A genome-wide association study using a Vietnamese landrace panel of rice (Oryza sativa) reveals new QTLs controlling panicle morphological traits

Optimized nitrogen management improves grain yield of rice by regulating panicle architecture in South China

Optimized nitrogen (N) management (OPT), with reduced total N input and more N applied during panicle development, has been proved to increase grain yield of rice through panicle enlargement. However, the changes in panicle architecture and source of variation are not well understood. A hybrid rice variety named Tianyou 3618 was subjected to OPT and farmer's fertilizer practice (FFP) in early cropping seasons of 2016 and 2017. With 16.7 % less N input, OPT increased panicle size by 8.6 % and 27.4 %, and grain yield by 13.8 % and 12.3 % for 2016 and 2017, respectively. OPT had greater dry matter accumulation and N uptake from panicle initiation to heading, which bolstered panicle enlargement. The number of surviving florets per branch was quite constant under different N treatments for all primary, secondary, and tertiary branches, implying that panicle size was mainly determined by the number of branches rather than the number of florets per branch. Little change was observed between OPT and FFP in differentiation, degeneration and survival of primary branches and their florets. Surviving secondary and tertiary branches and their florets were significantly more under OPT than those under FFP. The increase in surviving secondary branches under OPT resulted from both enhanced differentiation and reduced degeneration. While the increase in surviving tertiary branches under OPT was merely from enhanced differentiation though their degeneration was also dramatically increased. Among the increased differentiated florets under OPT, 32.4%-36.3 % and 61.6%-67.7 % came from secondary and tertiary branches, respectively. Among the increased surviving florets under OPT, 62.2%-65.2 % and 32.5%-37.8 % came from secondary and tertiary branches, respectively. Both secondary branches and tertiary branches were principal contributors to the increase in panicle size of OPT. To our knowledge, this is the first report on the detailed changes in panicle architecture and their involvement in panicle enlargement and yield gain under OPT.

New insights into QTNs and potential candidate genes governing rice yield via a multi-model genome-wide association study

BACKGROUND

Rice (Oryza sativa L.) is one of the globally important staple food crops, and yield-related traits are prerequisites for improved breeding efficiency in rice. Here, we used six different genome-wide association study (GWAS) models for 198 accessions, with 553,229 single nucleotide markers (SNPs) to identify the quantitative trait nucleotides (QTNs) and candidate genes (CGs) governing rice yield.

RESULTS

Amongst the 73 different QTNs in total, 24 were co-localized with already reported QTLs or loci in previous mapping studies. We obtained fifteen significant QTNs, pathway analysis revealed 10 potential candidates within 100kb of these QTNs that are predicted to govern plant height, days to flowering, and plot yield in rice. Based on their superior allelic information in 20 elite and 6 inferior genotypes, we found a higher percentage of superior alleles in the elite genotypes in comparison to inferior genotypes. Further, we implemented expression analysis and enrichment analysis enabling the identification of 73 candidate genes and 25 homologues of Arabidopsis, 19 of which might regulate rice yield traits. Of these candidate genes, 40 CGs were found to be enriched in 60 GO terms of the studied traits for instance, positive regulator metabolic process (GO:0010929), intracellular part (GO:0031090), and nucleic acid binding (GO:0090079). Haplotype and phenotypic variation analysis confirmed that LOC_OS09G15770, LOC_OS02G36710 and LOC_OS02G17520 are key candidates associated with rice yield.

CONCLUSIONS

Overall, we foresee that the QTNs, putative candidates elucidated in the study could summarize the polygenic regulatory networks controlling rice yield and be useful for breeding high-yielding varieties.

Genomic introgressions from African rice (Oryza glaberrima) in Asian rice (O. sativa) lead to the identification of key QTLs for panicle architecture

BACKGROUND

Developing high yielding varieties is a major challenge for breeders tackling the challenges of climate change in agriculture. The panicle (inflorescence) architecture of rice is one of the key components of yield potential and displays high inter- and intra-specific variability. The genus Oryza features two different crop species: Asian rice (Oryza sativa L.) and the African rice (O. glaberrima Steud.). One of the main morphological differences between the two independently domesticated species is the structure (or complexity) of the panicle, with O. sativa displaying a highly branched panicle, which in turn produces a larger number of grains than that of O. glaberrima. The gene regulatory network that governs intra- and interspecific panicle diversity is still under-studied.

RESULTS

To identify genetic factors linked to panicle architecture diversity in the two species, we used a set of 60 Chromosome Segment Substitution Lines (CSSLs) issued from third generation backcross (BC3DH) and carrying genomic segments from O. glaberrima cv. MG12 in the genetic background of O. sativa Tropical Japonica cv. Caiapó. Phenotypic data were collected for rachis and primary branch length, primary, secondary and tertiary branch number and spikelet number. A total of 15 QTLs were localized on chromosomes 1, 2, 3, 7, 11 and 12, QTLs associated with enhanced secondary and tertiary branch numbers were detected in two CSSLs. Furthermore, BC4F3:5 lines carrying different combinations of substituted segments were produced to decipher the effects of the identified QTL regions on variations in panicle architecture. A detailed analysis of phenotypes versus genotypes was carried out between the two parental genomes within these regions in order to understand how O. glaberrima introgression events may lead to alterations in panicle traits.

CONCLUSION

Our analysis led to the detection of genomic variations between O. sativa cv. Caiapó and O. glaberrima cv. MG12 in regions associated with enhanced panicle traits in specific CSSLs. These regions contain a number of key genes that regulate panicle development in O. sativa and their interspecific genomic variations may explain the phenotypic effects observed.

Genome-wide association analysis identifies natural allelic variants associated with panicle architecture variation in African rice, Oryza glaberrima Steud

African rice (Oryza glaberrima Steud), a short-day cereal crop closely related to Asian rice (Oryza sativa L.), has been cultivated in Sub-Saharan Africa for ∼ 3,000 years. Although less cultivated globally, it is a valuable genetic resource in creating high-yielding cultivars that are better adapted to diverse biotic and abiotic stresses. While inflorescence architecture, a key trait for rice grain yield improvement, has been extensively studied in Asian rice, the morphological and genetic determinants of this complex trait are less understood in African rice. In this study, using a previously developed association panel of 162 O. glaberrima accessions and new SNP variants characterized through mapping to a new version of the O. glaberrima reference genome, we conducted a genome-wide association study of four major morphological panicle traits. We have found a total of 41 stable genomic regions that are significantly associated with these traits, of which 13 co-localized with previously identified QTLs in O. sativa populations and 28 were unique for this association panel. Additionally, we found a genomic region of interest on chromosome 3 that was associated with the number of spikelets and primary and secondary branches. Within this region was localized the O. sativa ortholog of the PHYTOCHROME B gene (Oglab_006903/OgPHYB). Haplotype analysis revealed the occurrence of natural sequence variants at the OgPHYB locus associated with panicle architecture variation through modulation of the flowering time phenotype, whereas no equivalent alleles were found in O. sativa. The identification in this study of genomic regions specific to O. glaberrima indicates panicle-related intra-specific genetic variation in this species, increasing our understanding of the underlying molecular processes governing panicle architecture. Identified candidate genes and major haplotypes may facilitate the breeding of new African rice cultivars with preferred panicle traits.

Crop germplasm: Current challenges, physiological-molecular perspective, and advance strategies towards development of climate-resilient crops

Germplasm is a long-term resource management mission and investment for civilization. An estimated ∼7.4 million accessions are held in 1750 plant germplasm centres around the world; yet, only 2% of these assets have been utilized as plant genetic resources (PGRs). According to recent studies, the current food yield trajectory will be insufficient to feed the world's population in 2050. Additionally, possible negative effects in terms of crop failure because of climate change are already being experienced across the world. Therefore, it is necessary to reconciliation of research advancement and innovation of practices for further exploration of the potential of crop germplasm especially for the complex traits associated with yield such as water- and nitrogen use efficiency. In this review, we tried to address current challenges, research gaps, physiological and molecular aspects of two broad spectrum complex traits such as water- and nitrogen-use efficiency, and advanced integrated strategies that could provide a platform for combined stress management for climate-smart crop development. Additionally, recent development in technologies that are directly related to germplasm characterization was highlighted for further molecular utilization towards the development of elite varieties.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum)

We investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment. Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

Genes determining panicle morphology and grain quality in rice (Oryza sativa)

The world's increase in rice (Oryza sativa L.) production is not keeping up with the increase in its population. To boost the introduction of new high-yielding cultivars, knowledge is being gained on the genes and quantitative trait loci (QTLs) determining the panicle phenotype. The important are those determining yield of the crop, such as grain numbers per panicle and size and weight of the grains. Biochemical and molecular functions of many of them are understood in some details. Among these, OsCKX2 and OsSPL14 have been shown to increase panicle branching and grain numbers when overexpressed. Furthermore, miRNAs appear to play an important role in determining the panicle morphology by regulating the expressions of the genes like OsSPL14 and GRF4 involved in panicle branching and grain numbers and length. Mutations also greatly influence the grain shape and size. However, the information gained so far on the genetic regulation of grain filling and panicle morphology has not been successfully put into commercial application. Furthermore, the identification of the gene(s)/QTLs regulating panicle compactness is still lacking, which may enable the researchers to convert a compact-panicle cultivar into a lax/open one, and thereby increasing the chances of enhancing the yield of a desired compact-panicle cultivar obtained by the breeding effort.

Effect of Panicle Morphology on Grain Filling and Rice Yield: Genetic Control and Molecular Regulation

The demand for rice is likely to increase approximately 1.5 times by the year 2050. In contrast, the rice production is stagnant since the past decade as the ongoing rice breeding program is unable to increase the production further, primarily because of the problem in grain filling. Investigations have revealed several reasons for poor filling of the grains in the inferior spikelets of the compact panicle, which are otherwise genetically competent to develop into well-filled grains. Among these, the important reasons are 1) poor activities of the starch biosynthesizing enzymes, 2) high ethylene production leading to inhibition in expressions of the starch biosynthesizing enzymes, 3) insufficient division of the endosperm cells and endoreduplication of their nuclei, 4) low accumulation of cytokinins and indole-3-acetic acid (IAA) that promote grain filling, and 5) altered expressions of the miRNAs unfavorable for grain filling. At the genetic level, several genes/QTLs linked to the yield traits have been identified, but the information so far has not been put into perspective toward increasing the rice production. Keeping in view the genetic competency of the inferior spikelets to develop into well-filled grains and based on the findings from the recent research studies, improving grain filling in these spikelets seems plausible through the following biotechnological interventions: 1) spikelet-specific knockdown of the genes involved in ethylene synthesis and overexpression of β-CAS (β-cyanoalanine) for enhanced scavenging of CN⁻ formed as a byproduct of ethylene biosynthesis; 2) designing molecular means for increased accumulation of cytokinins, abscisic acid (ABA), and IAA in the caryopses; 3) manipulation of expression of the transcription factors like MYC and OsbZIP58 to drive the expression of the starch biosynthesizing enzymes; 4) spikelet-specific overexpression of the cyclins like CycB;1 and CycH;1 for promoting endosperm cell division; and 5) the targeted increase in accumulation of ABA in the straw during the grain filling stage for increased carbon resource remobilization to the grains. Identification of genes determining panicle compactness could also lead to an increase in rice yield through conversion of a compact-panicle into a lax/open one. These efforts have the ability to increase rice production by as much as 30%, which could be more than the set production target by the year 2050.

Genetic Diversity and Population Structure of Myanmar Rice (Oryza sativa L.) Varieties Using DArTseq-Based SNP and SilicoDArT Markers

Myanmar is well known as a primary center of plant genetic resources for rice. A considerable number of genetic diversity studies have been conducted in Myanmar using various DNA markers. However, this is the first report using DArTseq technology for exploring the genetic diversity of Myanmar rice. In our study, two ultra-high-throughput diversity array technology markers were employed to investigate the genetic diversity and population structure of local rice varieties in the Ayeyarwady delta, the major region of rice cultivation. The study was performed using 117 rice genotypes with 7643 SNP and 4064 silicoDArT markers derived from the DArT platform. Genetic variance among the genotypes ranged from 0 to 0.753 in SNPs, and from 0.001 to 0.954 in silicoDArT. Two distinct population groups were identified from SNP data analysis. Cluster analysis with both markers clearly separated traditional Pawsan varieties and modern high-yielding varieties. A significant divergence was found between populations according to the Fst values (0.737) obtained from the analysis of molecular variance, which revealed 74% genetic variation at the population level. These findings support rice researchers in identifying useful DNA polymorphisms in genes and pinpointing specific genes conferring desirable phenotypic traits for further genome-wide association studies and parental selection for recombination breeding to enhance rice varietal development and release.

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.

OsTPR boosts the superior grains through increase in upper secondary rachis branches without incurring a grain quality penalty

To address the future food security in Asia, we need to improve the genetic gain of grain yield while ensuring the consumer acceptance. This study aimed to identify novel genes influencing the number of upper secondary rachis branches (USRB) to elevate superior grains without compromising grain quality by studying the genetic variance of 310 diverse O. sativa var. indica panel using single- and multi-locus genome-wide association studies (GWAS), gene set analyses and gene regulatory network analysis. GWAS of USRB identified 230 significant (q-value < 0.05) SNPs from chromosomes 1 and 2. GWAS targets narrowed down using gene set analyses identified large effect association on an important locus LOC_Os02g50790/LOC_Os02g50799 encoding a nuclear-pore anchor protein (OsTPR). The superior haplotype derived from non-synonymous SNPs identified in OsTPR was specifically associated with increase in USRB with superior grains being low chalk. Through haplotype mining, we further demonstrated the synergy of offering added yield advantage due to superior allele of OsTPR in elite materials with low glycaemic index (GI) property. We further validated the importance of OsTPR using recombinant inbred lines (RILs) population by introgressing a superior allele of OsTPR into elite materials resulted in raise in productivity in high amylose background. This confirmed a critical role for OsTPR in influencing yield while maintaining grain and nutritional quality.

Importance of Landraces in Cereal Breeding for Stress Tolerance

The renewed focus on cereal landraces is a response to some negative consequences of modern agriculture and conventional breeding which led to a reduction of genetic diversity. Cereal landraces are still cultivated on marginal lands due to their adaptability to unfavourable conditions, constituting an important source of genetic diversity usable in modern plant breeding to improve the adaptation to abiotic or biotic stresses, yield performance and quality traits in limiting environments. Traditional agricultural production systems have played an important role in the evolution and conservation of wide variability in gene pools within species. Today, on-farm and ex situ conservation in gene bank collections, together with data sharing among researchers and breeders, will greatly benefit cereal improvement. Many efforts are usually made to collect, organize and phenotypically and genotypically analyse cereal landrace collections, which also utilize genomic approaches. Their use in breeding programs based on genomic selection, and the discovery of beneficial untapped QTL/genes/alleles which could be introgressed into modern varieties by MAS, pyramiding or biotechnological tools, increase the potential for their better deployment and exploitation in breeding for a more sustainable agricultural production, particularly enhancing adaptation and productivity in stress-prone environments to cope with current climate changes.

Resequencing of 672 Native Rice Accessions to Explore Genetic Diversity and Trait Associations in Vietnam

BACKGROUND

Vietnam possesses a vast diversity of rice landraces due to its geographical situation, latitudinal range, and a variety of ecosystems. This genetic diversity constitutes a highly valuable resource at a time when the highest rice production areas in the low-lying Mekong and Red River Deltas are enduring increasing threats from climate changes, particularly in rainfall and temperature patterns.

RESULTS

We analysed 672 Vietnamese rice genomes, 616 newly sequenced, that encompass the range of rice varieties grown in the diverse ecosystems found throughout Vietnam. We described four Japonica and five Indica subpopulations within Vietnam likely adapted to the region of origin. We compared the population structure and genetic diversity of these Vietnamese rice genomes to the 3000 genomes of Asian cultivated rice. The named Indica-5 (I5) subpopulation was expanded in Vietnam and contained lowland Indica accessions, which had very low shared ancestry with accessions from any other subpopulation and were previously overlooked as admixtures. We scored phenotypic measurements for nineteen traits and identified 453 unique genotype-phenotype significant associations comprising twenty-one QTLs (quantitative trait loci). The strongest associations were observed for grain size traits, while weaker associations were observed for a range of characteristics, including panicle length, heading date and leaf width.

CONCLUSIONS

We showed how the rice diversity within Vietnam relates to the wider Asian rice diversity by using a number of approaches to provide a clear picture of the novel diversity present within Vietnam, mainly around the Indica-5 subpopulation. Our results highlight differences in genome composition and trait associations among traditional Vietnamese rice accessions, which are likely the product of adaption to multiple environmental conditions and regional preferences in a very diverse country. Our results highlighted traits and their associated genomic regions that are a potential source of novel loci and alleles to breed a new generation of low input sustainable and climate resilient rice.

A cluster of Ankyrin and Ankyrin-TPR repeat genes is associated with panicle branching diversity in rice

The number of grains per panicle is an important yield-related trait in cereals which depends in part on panicle branching complexity. One component of this complexity is the number of secondary branches per panicle. Previously, a GWAS site associated with secondary branch and spikelet numbers per panicle in rice was identified. Here we combined gene capture, bi-parental genetic population analysis, expression profiling and transgenic approaches in order to investigate the functional significance of a cluster of 6 ANK and ANK-TPR genes within the QTL. Four of the ANK and ANK-TPR genes present a differential expression associated with panicle secondary branch number in contrasted accessions. These differential expression patterns correlate in the different alleles of these genes with specific deletions of potential cis-regulatory sequences in their promoters. Two of these genes were confirmed through functional analysis as playing a role in the control of panicle architecture. Our findings indicate that secondary branching diversity in the rice panicle is governed in part by differentially expressed genes within this cluster encoding ANK and ANK-TPR domain proteins that may act as positive or negative regulators of panicle meristem’s identity transition from indeterminate to determinate state.

Grain yield is one of the most important indexes in rice breeding, which is controlled in part by panicle branching complexity. A new QTL with co-location of spikelet number (SpN) and secondary branch number (SBN) traits was identified by genome-wide association study in a Vietnamese rice landrace panel. A set of four Ankyrin and Tetratricopeptide repeat domain-encoding genes was identified from this QTL based on their difference of expression levels between two contrasted haplotypes for the SpN and SBN traits. The differential expression is correlated with deletions in the promoter regions of these genes. Two of the genes act as negative regulators of the panicle meristem’s identity transition from indeterminate to determinate state while the other two act as positive regulators of this meristem fate transition. Based on the different phenotypes between overexpressed and mutant plants, two of these genes were confirmed as playing a role in the control of panicle architecture. These findings can be directly used to assist selection for grain yield improvement.

Genome-Wide Association Mapping of Salinity Tolerance at the Seedling Stage in a Panel of Vietnamese Landraces Reveals New Valuable QTLs for Salinity Stress Tolerance Breeding in Rice

Rice tolerance to salinity stress involves diverse and complementary mechanisms, such as the regulation of genome expression, activation of specific ion-transport systems to manage excess sodium at the cell or plant level, and anatomical changes that avoid sodium penetration into the inner tissues of the plant. These complementary mechanisms can act synergistically to improve salinity tolerance in the plant, which is then interesting in breeding programs to pyramidize complementary QTLs (quantitative trait loci), to improve salinity stress tolerance of the plant at different developmental stages and in different environments. This approach presupposes the identification of salinity tolerance QTLs associated with different mechanisms involved in salinity tolerance, which requires the greatest possible genetic diversity to be explored. To contribute to this goal, we screened an original panel of 179 Vietnamese rice landraces genotyped with 21,623 SNP markers for salinity stress tolerance under 100 mM NaCl treatment, at the seedling stage, with the aim of identifying new QTLs involved in the salinity stress tolerance via a genome-wide association study (GWAS). Nine salinity tolerance-related traits, including the salt injury score, chlorophyll and water content, and K⁺ and Na⁺ contents were measured in leaves. GWAS analysis allowed the identification of 26 QTLs. Interestingly, ten of them were associated with several different traits, which indicates that these QTLs act pleiotropically to control the different levels of plant responses to salinity stress. Twenty-one identified QTLs colocalized with known QTLs. Several genes within these QTLs have functions related to salinity stress tolerance and are mainly involved in gene regulation, signal transduction or hormone signaling. Our study provides promising QTLs for breeding programs to enhance salinity tolerance and identifies candidate genes that should be further functionally studied to better understand salinity tolerance mechanisms in rice.

Reference gene identification for gene expression analysis in rice under different metal stress

Real-time quantitative polymerase chain reaction (RT-qPCR) is the most common approach to quantify changes in gene expression. Appropriate internal reference genes are essential for normalization of data of RT-qPCR. In the present study, we identified suitable reference genes for analysis of gene expression in rice seedlings subjected to different heavy metal stresses such as deficiencies of iron and zinc and toxicities of cobalt, cadmium and nickel. First, from publically available RNA-Seq data we identified 10 candidate genes having stable expression. We also included commonly used house-keeping gene OsUBQ5 (Ubiquitin 5) in our analysis. Expression stability of all the 11 genes was determined by two independent tools, NormFinder and geNorm. Our results show that selected candidate reference genes have higher stability in their expression compared to that of OsUBQ5. Genes with locus ID LOC_Os03g16690, encoding an oxysterol-binding protein (OsOBP) and LOC_Os01g56580, encoding Casein Kinase_1a.3 (OsCK1a.3) were identified to be the most stably expressed reference genes under most of the conditions tested. Finally, the study reveals that it is better to use a specific reference gene for a specific heavy metal stress condition rather than using a common reference gene for multiple heavy metal stress conditions. The reference genes identified here would be very useful for gene expression studies under heavy metal stresses in rice.

Genome-wide association study and Mendelian randomization analysis provide insights for improving rice yield potential

Rice yield per plant has a complex genetic architecture, which is mainly determined by its three component traits: the number of grains per panicle (GPP), kilo-grain weight (KGW), and tillers per plant (TP). Exploring ideotype breeding based on selection for genetically less complex component traits is an alternative route for further improving rice production. To understand the genetic basis of the relationship between rice yield and component traits, we investigated the four traits of two rice hybrid populations (575 + 1495 F₁) in different environments and conducted meta-analyses of genome-wide association study (meta-GWAS). In total, 3589 significant loci for three components traits were detected, while only 3 loci for yield were detected. It indicated that rice yield is mainly controlled by minor-effect loci and hardly to be identified. Selecting quantitative trait locus/gene affected component traits to further enhance yield is recommended. Mendelian randomization design is adopted to investigate the genetic effects of loci on yield through component traits and estimate the genetic relationship between rice yield and its component traits by these loci. The loci for GPP or TP mainly had a positive genetic effect on yield, but the loci for KGW with different direction effects (positive effect or negative effect). Additionally, TP (Beta = 1.865) has a greater effect on yield than KGW (Beta = 1.016) and GPP (Beta = 0.086). Five significant loci for component traits that had an indirect effect on yield were identified. Pyramiding superior alleles of the five loci revealed improved yield. A combination of direct and indirect effects may better contribute to the yield potential of rice. Our findings provided a rationale for using component traits as indirect indices to enhanced rice yield, which will be helpful for further understanding the genetic basis of yield and provide valuable information for improving rice yield potential.

Discovery of new genetic determinants of morphological plasticity in rice roots and shoots under phosphate starvation using GWAS

Phosphorus is an essential nutrient for plants that is often in short supply. In rice (Oryza sativa L.), inorganic phosphate (P_i) deficiency leads to various physiological disorders that consequently affect plant productivity. In this study, a large-scale phenotyping experiment using 160 Vietnamese rice landraces was performed under greenhouse conditions, by employing an alpha lattice design with three replicates, to identify quantitative trait loci (QTLs) associated with plant growth inhibition caused by P_i deficiency. Rice plantlets were grown for six weeks in the PVC sand column (16 cm diameter × 80 cm height) supplied with P_i-deficient medium (10 μM P) or full-P_i Yoshida medium (320 μM P). The effects of P_i deficiency on the number of crown roots, root length, shoot length, root weight, shoot weight and total weight were studied. From 36 significant markers identified using a genome-wide association study, 21 QTLs associated with plant growth inhibition under P_i starvation were defined. In total, 158 candidate genes co-located with the defined QTLs were identified. Interestingly, one QTL (qRST9.14) was associated with all three weight-traits. The co-located gene GLYCEROPHOSPHODIESTER PHOSPHODIESTERASE 13 was found to be potentially involved in P_i transport. Understanding the molecular mechanisms of Pi-starvation responses, and identifying the potential QTLs responsible for low-P_i stress tolerance, will provide valuable information for developing new varieties tolerant of low-P_i conditions.

Cell cycle events and expression of cell cycle regulators are determining factors in differential grain filling in rice spikelets based on their spatial location on compact panicles

Rice being a staple crop for human, its production is required to be increased significantly, particularly keeping in view the expected world's population of 9.6 billion by the year 2050. In this context, although the rice breeding programs have been successful in increasing the number of spikelets per panicle, the basal spikelets remain poorly filled, undermining the yield potential. The present study also found the grain filling to bear negative correlation with the panicle grain density. The poorly filled basal spikelets of the compact-panicle cultivars showed a lower endosperm cell division rate and ploidy status of the endosperm nuclei coupled with no significant greater expression of CYCB;1 and CYCH;1 compared with the apical spikelets, unlike that observed in the lax-panicle cultivars, which might have prevented them from overcoming apical dominance. Significantly greater expression of CYCB2;2 in the basal spikelets than in the apical spikelets might also have prevented the former to enter into endoreduplication. Furthermore, expression studies of KRPs in the caryopses revealed that a higher expression of KRP;1 and KRP;4 in the basal spikelets than in the apical spikelets of the compact-panicle cultivars could also be detrimental to grain filling in the former, as KRPs form complex primarily with CDKA-CYCD that promotes S-phase activity and G1/S transition, and thus inhibits endosperm cell division. The study indicates that targeted manipulation of expression of CYCB1;1, CYCB2;2, CYCH1;1, KRP;1 and KRP4 in the basal spikelets of the compact-panicle cultivars may significantly improve their yield performance.

Genome Wide Association Study Pinpoints Key Agronomic QTLs in African Rice Oryza glaberrima

BACKGROUND

African rice, Oryza glaberrima, is an invaluable resource for rice cultivation and for the improvement of biotic and abiotic resistance properties. Since its domestication in the inner Niger delta ca. 2500 years BP, African rice has colonized a variety of ecologically and climatically diverse regions. However, little is known about the genetic basis of quantitative traits and adaptive variation of agricultural interest for this species.

RESULTS

Using a reference set of 163 fully re-sequenced accessions, we report the results of a Genome Wide Association Study carried out for African rice. We investigated a diverse panel of traits, including flowering date, panicle architecture and resistance to Rice yellow mottle virus. For this, we devised a pipeline using complementary statistical association methods. First, using flowering time as a target trait, we found several association peaks, one of which co-localised with a well described gene in the Asian rice flowering pathway, OsGi, and identified new genomic regions that would deserve more study. Then we applied our pipeline to panicle- and resistance-related traits, highlighting some interesting genomic regions and candidate genes. Lastly, using a high-resolution climate database, we performed an association analysis based on climatic variables, searching for genomic regions that might be involved in adaptation to climatic variations.

CONCLUSION

Our results collectively provide insights into the extent to which adaptive variation is governed by sequence diversity within the O. glaberrima genome, paving the way for in-depth studies of the genetic basis of traits of interest that might be useful to the rice breeding community.

Plant Inflorescence Architecture: The Formation, Activity, and Fate of Axillary Meristems

The above-ground plant body in different plant species can have very distinct forms or architectures that arise by recurrent redeployment of a finite set of building blocks-leaves with axillary meristems, stems or branches, and flowers. The unique architectures of plant inflorescences in different plant families and species, on which this review focuses, determine the reproductive success and yield of wild and cultivated plants. Major contributors to the inflorescence architecture are the activity and developmental trajectories adopted by axillary meristems, which determine the degree of branching and the number of flowers formed. Recent advances in genetic and molecular analyses in diverse flowering plants have uncovered both common regulatory principles and unique players and/or regulatory interactions that underlie inflorescence architecture. Modulating activity of these regulators has already led to yield increases in the field. Additional insight into the underlying regulatory interactions and principles will not only uncover how their rewiring resulted in altered plant form, but will also enhance efforts at optimizing plant architecture in desirable ways in crop species.

PI-Plat: a high-resolution image-based 3D reconstruction method to estimate growth dynamics of rice inflorescence traits

BACKGROUND

Recent advances in image-based plant phenotyping have improved our capability to study vegetative stage growth dynamics. However, more complex agronomic traits such as inflorescence architecture (IA), which predominantly contributes to grain crop yield are more challenging to quantify and hence are relatively less explored. Previous efforts to estimate inflorescence-related traits using image-based phenotyping have been limited to destructive end-point measurements. Development of non-destructive inflorescence phenotyping platforms could accelerate the discovery of the phenotypic variation with respect to inflorescence dynamics and mapping of the underlying genes regulating critical yield components.

RESULTS

The major objective of this study is to evaluate post-fertilization development and growth dynamics of inflorescence at high spatial and temporal resolution in rice. For this, we developed the Panicle Imaging Platform (PI-Plat) to comprehend multi-dimensional features of IA in a non-destructive manner. We used 11 rice genotypes to capture multi-view images of primary panicle on weekly basis after the fertilization. These images were used to reconstruct a 3D point cloud of the panicle, which enabled us to extract digital traits such as voxel count and color intensity. We found that the voxel count of developing panicles is positively correlated with seed number and weight at maturity. The voxel count from developing panicles projected overall volumes that increased during the grain filling phase, wherein quantification of color intensity estimated the rate of panicle maturation. Our 3D based phenotyping solution showed superior performance compared to conventional 2D based approaches.

CONCLUSIONS

For harnessing the potential of the existing genetic resources, we need a comprehensive understanding of the genotype-to-phenotype relationship. Relatively low-cost sequencing platforms have facilitated high-throughput genotyping, while phenotyping, especially for complex traits, has posed major challenges for crop improvement. PI-Plat offers a low cost and high-resolution platform to phenotype inflorescence-related traits using 3D reconstruction-based approach. Further, the non-destructive nature of the platform facilitates analyses of the same panicle at multiple developmental time points, which can be utilized to explore the genetic variation for dynamic inflorescence traits in cereals.

Characterizing 3D inflorescence architecture in grapevine using X-ray imaging and advanced morphometrics: implications for understanding cluster density

Inflorescence architecture provides the scaffold on which flowers and fruits develop, and consequently is a primary trait under investigation in many crop systems. Yet the challenge remains to analyse these complex 3D branching structures with appropriate tools. High information content datasets are required to represent the actual structure and facilitate full analysis of both the geometric and the topological features relevant to phenotypic variation in order to clarify evolutionary and developmental inflorescence patterns. We combined advanced imaging (X-ray tomography) and computational approaches (topological and geometric data analysis and structural simulations) to comprehensively characterize grapevine inflorescence architecture (the rachis and all branches without berries) among 10 wild Vitis species. Clustering and correlation analyses revealed unexpected relationships, for example pedicel branch angles were largely independent of other traits. We identified multivariate traits that typified species, which allowed us to classify species with 78.3% accuracy, versus 10% by chance. Twelve traits had strong signals across phylogenetic clades, providing insight into the evolution of inflorescence architecture. We provide an advanced framework to quantify 3D inflorescence and other branched plant structures that can be used to tease apart subtle, heritable features for a better understanding of genetic and environmental effects on plant phenotypes.

Unraveling the Genetic Elements Involved in Shoot and Root Growth Regulation by Jasmonate in Rice Using a Genome-Wide Association Study

BACKGROUND

Due to their sessile life style, plant survival is dependent on the ability to build up fast and highly adapted responses to environmental stresses by modulating defense response and organ growth. The phytohormone jasmonate plays an essential role in regulating these plant responses to stress.

RESULTS

To assess variation of plant growth responses and identify genetic determinants associated to JA treatment, we conducted a genome-wide association study (GWAS) using an original panel of Vietnamese rice accessions. The phenotyping results showed a high natural genetic variability of the 155 tested rice accessions in response to JA for shoot and root growth. The level of growth inhibition by JA is different according to the rice varieties tested. We conducted genome-wide association study and identified 28 significant associations for root length (RTL), shoot length (SHL), root weight (RTW), shoot weight (SHW) and total weight (TTW) in response to JA treatment. Three common QTLs were found for RTL, RTW and SHL. Among a list of 560 candidate genes found to co-locate with the QTLs, a transcriptome analysis from public database for the JA response allows us to identify 232 regulated genes including several JA-responsive transcription factors known to play a role in stress response.

CONCLUSION

Our genome-wide association study shows that common and specific genetic elements are associated with inhibition of shoot and root growth under JA treatment suggesting the involvement of a complex JA-dependent genetic control of rice growth inhibition at the whole plant level. Besides, numerous candidate genes associated to stress and JA response are co-located with the association loci, providing useful information for future studies on genetics and breeding to optimize the growth-defense trade-off in rice.

Genome-wide association mapping of leaf mass traits in a Vietnamese rice landrace panel

Leaf traits are often strongly correlated with yield, which poses a major challenge in rice breeding. In the present study, using a panel of Vietnamese rice landraces genotyped with 21,623 single-nucleotide polymorphism markers, a genome-wide association study (GWAS) was conducted for several leaf traits during the vegetative stage. Vietnamese landraces are often poorly represented in panels used for GWAS, even though they are adapted to contrasting agrosystems and can contain original, valuable genetic determinants. A panel of 180 rice varieties was grown in pots for four weeks with three replicates under nethouse conditions. Different leaf traits were measured on the second fully expanded leaf of the main tiller, which often plays a major role in determining the photosynthetic capacity of the plant. The leaf fresh weight, turgid weight and dry weight were measured; then, from these measurements, the relative tissue weight and leaf dry matter percentage were computed. The leaf dry matter percentage can be considered a proxy for the photosynthetic efficiency per unit leaf area, which contributes to yield. By a GWAS, thirteen QTLs associated with these leaf traits were identified. Eleven QTLs were identified for fresh weight, eleven for turgid weight, one for dry weight, one for relative tissue weight and one for leaf dry matter percentage. Eleven QTLs presented associations with several traits, suggesting that these traits share common genetic determinants, while one QTL was specific to leaf dry matter percentage and one QTL was specific to relative tissue weight. Interestingly, some of these QTLs colocalize with leaf- or yield-related QTLs previously identified using other material. Several genes within these QTLs with a known function in leaf development or physiology are reviewed.

PI-Plat: a high-resolution image-based 3D reconstruction method to estimate growth dynamics of rice inflorescence traits

BACKGROUND

Recent advances in image-based plant phenotyping have improved our capability to study vegetative stage growth dynamics. However, more complex agronomic traits such as inflorescence architecture (IA), which predominantly contributes to grain crop yield are more challenging to quantify and hence are relatively less explored. Previous efforts to estimate inflorescence-related traits using image-based phenotyping have been limited to destructive end-point measurements. Development of non-destructive inflorescence phenotyping platforms could accelerate the discovery of the phenotypic variation with respect to inflorescence dynamics and mapping of the underlying genes regulating critical yield components.

RESULTS

The major objective of this study is to evaluate post-fertilization development and growth dynamics of inflorescence at high spatial and temporal resolution in rice. For this, we developed the Panicle Imaging Platform (PI-Plat) to comprehend multi-dimensional features of IA in a non-destructive manner. We used 11 rice genotypes to capture multi-view images of primary panicle on weekly basis after the fertilization. These images were used to reconstruct a 3D point cloud of the panicle, which enabled us to extract digital traits such as voxel count and color intensity. We found that the voxel count of developing panicles is positively correlated with seed number and weight at maturity. The voxel count from developing panicles projected overall volumes that increased during the grain filling phase, wherein quantification of color intensity estimated the rate of panicle maturation. Our 3D based phenotyping solution showed superior performance compared to conventional 2D based approaches.

CONCLUSIONS

For harnessing the potential of the existing genetic resources, we need a comprehensive understanding of the genotype-to-phenotype relationship. Relatively low-cost sequencing platforms have facilitated high-throughput genotyping, while phenotyping, especially for complex traits, has posed major challenges for crop improvement. PI-Plat offers a low cost and high-resolution platform to phenotype inflorescence-related traits using 3D reconstruction-based approach. Further, the non-destructive nature of the platform facilitates analyses of the same panicle at multiple developmental time points, which can be utilized to explore the genetic variation for dynamic inflorescence traits in cereals.