High-resolution QTL mapping for grain appearance traits and co-localization of chalkiness-associated differentially expressed candidate genes in rice

OsCRLK2, a Receptor-Like Kinase Identified by QTL Analysis, is Involved in the Regulation of Rice Quality

The quality of rice (Oryza sativa L) is determined by a combination of appearance, flavor, aroma, texture, storage characteristics, and nutritional composition. Rice quality directly influences acceptance by consumers and commercial value. The genetic mechanism underlying rice quality is highly complex, and is influenced by genotype, environment, and chemical factors such as starch type, protein content, and amino acid composition. Minor variations in these chemical components may lead to substantial differences in rice quality. Among these components, starch is the most crucial and influential factor in determining rice quality. In this study, quantitative trait loci (QTLs) associated with eight physicochemical properties related to the rapid viscosity analysis (RVA) profile were identified using a high-density sequence map constructed using recombinant inbred lines (RILs). Fifty-nine QTLs were identified across three environments, among which qGT6.4 was a novel locus co-located across all three environments. By integrating RNA-seq data, we identified the differentially expressed candidate gene OsCRLK2 within the qGT6.4 interval. osclrk2 mutants exhibited decreased gelatinization temperature (GT), apparent amylose content (AAC) and viscosity, and increased chalkiness. Furthermore, osclrk2 mutants exhibited downregulated expression of the majority of starch biosynthesis-related genes compared to wild type (WT) plants. In summary, OsCRLK2, which encodes a receptor-like protein kinase, appears to consistently influence rice quality across different environments. This discovery provides a new genetic resource for use in the molecular breeding of rice cultivars with improved quality.

Genetic dissection of cold tolerance at the budding stage of rice in an indica-japonica recombination inbred line population

Rice is highly cold-sensitive, and thus, the promotion of cold resistance in buds is essential. In this study, we conducted a mapping analysis to identify quantitative trait loci (QTLs) associated with cold tolerance in buds. The analysis was performed using a recombinant inbred line (RIL) population consisting of 192 lines derived from the cold-tolerant strain 02428 and the cold-sensitive strain YZX. Seven additive loci on chromosomes 1, 4, 5, and 6 were identified, of which loci 3 and 7 were found in two crop seasons, indicating stability. Three epistatic interactions, one present over two seasons, were found. Loci 3 and 7 pyramided with two main-effect QTLs observed to control the rate of low-temperature germination in our previous study. Two materials with good cold resistance at the germination and bud stages were obtained, namely, G93 and G146. Transcriptome sequencing analysis of the two parent buds after cold treatment found that genes expressed differentially between the two parents were related to photosynthesis, energy metabolism, and reactive oxygen scavenging. Five candidate genes, namely, Os01g0385400, Os01g0388000, Os06g0287700, Os06g0289200, and Os06g0291100, were selected in the two stable intervals based on gene expression profiles and annotations. These genetic loci exhibit strong potential as targets for breeding cold tolerance in buds and require additional investigation. In conclusion, this work provides valuable genetic resources that can be utilized to improve the cold tolerance of rice.

Construction of High-Density Genetic Map and QTL Mapping for Grain Shape in the Rice RIL Population

Grain shape is an important agronomic trait directly associated with yield in rice. In order to explore new genes related to rice grain shape, a high-density genetic map containing 2193 Bin markers (526957 SNP) was constructed by whole-genome resequencing of 208 recombinant inbred (RILs) derived from a cross between ZP37 and R8605, with a total genetic distance of 1542.27 cM. The average genetic distance between markers was 0.76 cM, and the physical distance was 201.29 kb. Quantitative trait locus (QTL) mapping was performed for six agronomic traits related to rice grain length, grain width, length-to-width ratio, thousand-grain weight, grain cross-sectional area, and grain perimeter under three different environments. A total of 39 QTLs were identified, with mapping intervals ranging from 8.1 kb to 1781.6 kb and an average physical distance of 517.5 kb. Among them, 15 QTLs were repeatedly detected in multiple environments. Analysis of the genetic effects of the identified QTLs revealed 14 stable genetic loci, including three loci that overlapped with previously reported gene positions, and the remaining 11 loci were newly identified loci associated with two or more environments or traits. Locus 1, Locus 3, Locus 10, and Locus 14 were novel loci exhibiting pleiotropic effects on at least three traits and were detected in multiple environments. Locus 14, with a contribution rate greater than 10%, influenced grain width, length-to-width ratio, and grain cross-sectional area. Furthermore, pyramiding effects analysis of three stable genetic loci showed that increasing the number of QTL could effectively improve the phenotypic value of grain shape. Collectively, our findings provided a theoretical basis and genetic resources for the cloning, functional analysis, and molecular breeding of genes related to rice grain shape.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Mapping of qChalk1 controlling grain chalkiness in japonica rice

BACKGROUND

Rice grain chalkiness is an undesirable characteristic that affects grain quality. The aim of this study was to map QTLs controlling grain chalkiness in japonica rice.

METHODS AND RESULTS

In this study, two japonica rice cultivars with similar grain shapes but different grain chalkiness rates were crossed and the F₂ and BC₁F₂ populations were subjected to QTL-seq analysis to map the QTLs controlling the grain chalkiness rate. QTL-seq analysis revealed SNP index differences on chromosome 1 in both of the segregating populations. Using polymorphic markers between the two parents, QTL mapping was conducted on 213 individual plants in the BC₁F₂ population. QTL mapping confined a QTL controlling grain chalkiness, qChalk1, to a 1.1 Mb genomic region on chromosome 1. qChalk1 explained 19.7% of the phenotypic variation.

CONCLUSION

A QTL controlling grain chalkiness qChalk1 was detected in both F₂ and BC₁F₂ segregating populations by QTL-Seq and QTL mapping methods. This result would be helpful for further cloning of the genes controlling grain chalkiness in japonica rice.

Molecular mapping and characterization of QTLs for grain quality traits in a RIL population of US rice under high nighttime temperature stress

Elevated nighttime temperatures resulting from climate change significantly impact the rice crop worldwide. The rice (Oryza sativa L.) plant is highly sensitive to high nighttime temperature (HNT) during grain-filling (reproductive stage). HNT stress negatively affects grain quality traits and has a major impact on the value of the harvested rice crop. In addition, along with grain dimensions determining rice grain market classes, the grain appearance and quality traits determine the rice grain market value. During the last few years, there has been a major concern for rice growers and the rice industry over the prevalence of rice grains opacity and the reduction of grain dimensions affected by HNT stress. Hence, the improvement of heat-stress tolerance to maintain grain quality of the rice crop under HNT stress will bolster future rice value in the market. In this study, 185 F₁₂-recombinant inbred lines (RILs) derived from two US rice cultivars, Cypress (HNT-tolerant) and LaGrue (HNT-sensitive) were screened for the grain quality traits grain length (GL), grain width (GW), and percent chalkiness (%chalk) under control and HNT stress conditions and evaluated to identify the genomic regions associated with the grain quality traits. In total, there were 15 QTLs identified; 6 QTLs represented under control condition explaining 3.33% to 8.27% of the phenotypic variation, with additive effects ranging from - 0.99 to 0.0267 on six chromosomes and 9 QTLs represented under HNT stress elucidating 6.39 to 51.53% of the phenotypic variation, with additive effects ranging from - 8.8 to 0.028 on nine chromosomes for GL, GW, and % chalk. These 15 QTLs were further characterized and scanned for natural genetic variation in a japonica diversity panel (JDP) to identify candidate genes for GL, GW, and %chalk. We found 6160 high impact single nucleotide polymorphisms (SNPs) characterized as such depending on their type, region, functional class, position, and proximity to the gene and/or gene features, and 149 differentially expressed genes (DEGs) in the 51 Mbp genomic region comprising of the 15 QTLs. Out of which, 11 potential candidate genes showed high impact SNP associations. Therefore, the analysis of the mapped QTLs and their genetic dissection in the US grown Japonica rice genotypes at genomic and transcriptomic levels provide deep insights into genetic variation beneficial to rice breeders and geneticists for understanding the mechanisms related to grain quality under heat stress in rice.

Fine mapping of the grain chalkiness quantitative trait locus qCGP6 reveals the involvement of Wx in grain chalkiness formation

Grain chalkiness is an important index of rice appearance quality and is negatively associated with rice processing and eating qualities. However, the genetic mechanism underlying chalkiness formation is largely unknown. To identify the genetic basis of chalkiness, 410 recombinant inbred lines (RILs) derived from two representative indica rice varieties, Shuhui498 (R498) and Yihui3551 (R3551), were used to discover quantitative trait loci (QTL). The two parental lines and RILs were grown in three locations in China under three controlled fertilizer application level. Analyses indicated that chalkiness was significantly affected by genotype, the environment, and the interaction between the two, and that heritability was high. Several QTLs were isolated, including the two stable QTLs, i.e., qCGP6 and qCGP8. Fine mapping and candidate gene verification of qCGP6 showed that Wx may play a key role in chalkiness formation. Chromosomal segment substitution lines (CSSLs) and near-isogenic lines (NILs) carrying the Wxa or Wxin allele produced more chalky grain than the R498 parent. A similar result was also observed in the 3611 background. Notably, the effect of the Wx genotype on rice chalkiness was shown to be dependent on environmental conditions and Wx alleles exhibited different sensitivities to shading treatment. Using CRISPR/Cas9, the Wxa promoter region was successfully edited, down-regulating Wx alleviates chalkiness formation in NILR498-Wxa. This study developed a new strategy for synergistic improvement of eating and appearance qualities in rice, and created a novel Wx allele with great potential in breeding applications.

Meta-QTL and haplo-pheno analysis reveal superior haplotype combinations associated with low grain chalkiness under high temperature in rice

Chalk, an undesirable grain quality trait in rice, is primarily formed due to high temperatures during the grain-filling process. Owing to the disordered starch granule structure, air spaces and low amylose content, chalky grains are easily breakable during milling thereby lowering head rice recovery and its market price. Availability of multiple QTLs associated with grain chalkiness and associated attributes, provided us an opportunity to perform a meta-analysis and identify candidate genes and their alleles contributing to enhanced grain quality. From the 403 previously reported QTLs, 64 Meta-QTLs encompassing 5262 non-redundant genes were identified. MQTL analysis reduced the genetic and physical intervals and nearly 73% meta-QTLs were narrower than 5cM and 2Mb, revealing the hotspot genomic regions. By investigating expression profiles of 5262 genes in previously published datasets, 49 candidate genes were shortlisted on the basis of their differential regulation in at least two of the datasets. We identified non-synonymous allelic variations and haplotypes in 39 candidate genes across the 3K rice genome panel. Further, we phenotyped a subset panel of 60 rice accessions by exposing them to high temperature stress under natural field conditions over two Rabi cropping seasons. Haplo-pheno analysis uncovered haplotype combinations of two starch synthesis genes, GBSSI and SSIIa, significantly contributing towards the formation of grain chalk in rice. We, therefore, report not only markers and pre-breeding material, but also propose superior haplotype combinations which can be introduced using either marker-assisted breeding or CRISPR-Cas based prime editing to generate elite rice varieties with low grain chalkiness and high HRY traits.

High-resolution quantitative trait locus mapping for rice grain quality traits using genotyping by sequencing

Rice is a major food crop that sustains approximately half of the world population. Recent worldwide improvements in the standard of living have increased the demand for high-quality rice. Accurate identification of quantitative trait loci (QTLs) for rice grain quality traits will facilitate rice quality breeding and improvement. In the present study, we performed high-resolution QTL mapping for rice grain quality traits using a genotyping-by-sequencing approach. An F₂ population derived from a cross between an elite japonica variety, Koshihikari, and an indica variety, Nona Bokra, was used to construct a high-density genetic map. A total of 3,830 single nucleotide polymorphism markers were mapped to 12 linkage groups spanning a total length of 2,456.4 cM, with an average genetic distance of 0.82 cM. Seven grain quality traits-the percentage of whole grain, percentage of head rice, percentage of area of head rice, transparency, percentage of chalky rice, percentage of chalkiness area, and degree of chalkiness-of the F₂ population were investigated. In total, 15 QTLs with logarithm of the odds (LOD) scores >4 were identified, which mapped to chromosomes 6, 7, and 9. These loci include four QTLs for transparency, four for percentage of chalky rice, four for percentage of chalkiness area, and three for degree of chalkiness, accounting for 0.01%-61.64% of the total phenotypic variation. Of these QTLs, only one overlapped with previously reported QTLs, and the others were novel. By comparing the major QTL regions in the rice genome, several key candidate genes reported to play crucial roles in grain quality traits were identified. These findings will expedite the fine mapping of these QTLs and QTL pyramiding, which will facilitate the genetic improvement of rice grain quality.

Utilization of Genotyping-by-Sequencing (GBS) for Rice Pre-Breeding and Improvement: A Review

Molecular markers play a crucial role in the improvement of rice. To benefit from these markers, genotyping is carried out to identify the differences at a specific position in the genome of individuals. The advances in sequencing technologies have led to the development of different genotyping techniques such as genotyping-by-sequencing. Unlike PCR-fragment-based genotyping, genotyping-by-sequencing has enabled the parallel sequencing and genotyping of hundreds of samples in a single run, making it more cost-effective. Currently, GBS is being used in several pre-breeding programs of rice to identify beneficial genes and QTL from different rice genetic resources. In this review, we present the current advances in the utilization of genotyping-by-sequencing for the development of rice pre-breeding materials and the improvement of existing rice cultivars. The challenges and perspectives of using this approach are also highlighted.

Fine Mapping of Two Major Quantitative Trait Loci for Rice Chalkiness With High Temperature-Enhanced Additive Effects

Chalkiness is a crucial determinant of rice quality. During seed filling period, high temperature usually increases grain chalkiness, resulting in poor grain quality. Rice chalkiness was controlled by quantitative trait loci (QTLs) and influenced by environmental conditions. In this study, we identified two single-segment substitution lines (SSSLs) 22-05 and 15-06 with significantly lower percentage of grain chalkiness (PGC) than recipient Huajingxian 74 (HJX74) over 6 cropping seasons. Two major QTLs for chalkiness, qPGC5 and qPGC6, were located by substitution mapping of SSSLs 22-05 and 15-06, respectively. qPGC5 was located in the 876.5 kb interval of chromosome 5 and qPGC6 was located in the 269.1 kb interval of chromosome 6. Interestingly, the PGC of HJX74 was significantly different between the two cropping seasons per year, with 25.8% in the first cropping season (FCS) and 16.6% in the second cropping season (SCS), while the PGC of SSSLs 22-05 and 15-06 did not significantly differ between FCS and SCS. The additive effects of qPGC5 and qPGC6 on chalkiness in the SSSLs were significantly greater in FCS than in SCS. These results showed that qPGC5 and qPGC6 had major effects on chalkiness and the SSSL alleles were more effective in reducing chalkiness under high temperature condition in FCS. The fine-mapping of the two QTLs will facilitate the cloning of genes for chalkiness and provide new genetic resources to develop new cultivars with low chalkiness even under high temperature condition.

Substitution Mapping of Two Closely Linked QTLs on Chromosome 8 Controlling Grain Chalkiness in Rice

Rice varieties are required to have high yield and good grain quality. Grain chalkiness and grain shape are two important traits of rice grain quality. Low chalkiness slender grains are preferred by most rice consumers. Here, we dissected two closely linked quantitative trait loci (QTLs) controlling grain chalkiness and grain shape on rice chromosome 8 by substitution mapping. Two closely linked QTLs controlling grain chalkiness and grain shape were identified using single-segment substitution lines (SSSLs). The two QTLs were then dissected on rice chromosome 8 by secondary substitution mapping. qPGC8.1 was located in an interval of 1382.6 kb and qPGC8.2 was mapped in a 2057.1 kb region. The maximum distance of the two QTLs was 4.37 Mb and the space distance of two QTL intervals was 0.72 Mb. qPGC8.1 controlled grain chalkiness and grain width. qPGC8.2 was responsible for grain chalkiness, grain length and width. The additive effects of qPGC8.1 and qPGC8.2 on grain chalkiness were not affected by higher temperature. Two closely linked QTLs qPGC8.1 and qPGC8.2 were dissected on rice chromosome 8. They controlled the phenotypes of grain chalkiness and grain shape. The two QTLs were insensitive to higher temperature.

Identification and candidate gene screening of qCIR9.1, a novel QTL associated with anther culturability in rice (Oryza sativa L.)

A novel QTL, qCIR9.1, that controls callus induction rate in anther culture was identified on chromosome 9 in rice, and based on RNA-seq data, Os09g0551600 was the most promising candidate gene. Anther culture, a doubled haploid (DH) technique, has become an important technology in many plant-breeding programmes. Although anther culturability is the key factor in this technique, its genetic mechanisms in rice remain poorly understood. In this study, we mapped quantitative trait loci (QTLs) responsible for anther culturability by using 192 recombinant inbred lines (RILs) derived from YZX (Oryza sativa ssp. indica) × 02428 (Oryza sativa ssp. japonica) and a high-density bin map. A total of eight QTLs for anther culturability were detected in three environments. Among these QTLs, a novel major QTL for callus induction rate (CIR) named qCIR9.1 was repeatedly mapped to a ~ 100 kb genomic interval on chromosome 9 and explained 8.39-14.14% of the phenotypic variation. Additionally, RNA sequencing (RNA-seq) was performed for the parents (YZX and 02428), low- (L-Pool) and high-CIR RILs (H-Pool) after 16 and 26 days of culture. By using the RNA of the bulked RILs for background normalization, the number of differentially expressed genes (DEGs) both between the parents and between the bulked RILs after 26 days of culture was drastically reduced to only 78. Among these DEGs, only one gene, Os09g0551600, encoding a high-mobility group (HMG) protein, was located in the candidate region of qCIR9.1. qRT-PCR analysis of Os09g0551600 showed the same results as RNA-seq, and the expression of this gene was decreased in the low-callus-induction parent (YZX) and L-Pool. Our results provide a foundational step for further cloning of qCIR9.1 and will be very useful for improving anther culturability in rice.

Genome-wide association coupled gene to gene interaction studies unveil novel epistatic targets among major effect loci impacting rice grain chalkiness

Rice varieties whose quality is graded as excellent have a lower percent grain chalkiness (PGC) of two per cent and below with higher whole grain yields upon milling, leading to higher economic returns for farmers. We have conducted a genome-wide association study (GWAS) using a combined population panel of indica and japonica rice varieties, and identified a total of 746 single nucleotide polymorphisms (SNPs) that were strongly associated with the chalk phenotype, covered 78 Quantitative Trait Loci (QTL) regions. Among them, 21 were high-value QTLs, as they explained at least 10 % of the phenotypic variance for PGC. A combined epistasis and GWAS was applied to dissect the genetics of the complex chalkiness trait, and its regulatory cascades were validated using gene regulatory networks. Promising novel epistatic interactions were found between the loci of chromosomes 6 (PGC6.1) and 7 (PGC7.8) that contributed to lower PGC. Based on haplotype mining only a few modern rice varieties confounded with a lower chalkiness, and they possess several PGC QTLs. The importance of PGC6.1 was validated through multi-parent advanced generation intercrosses and several low-chalk lines possessing superior haplotypes were identified. The results of this investigation have deciphered the underlying genetic networks that can reduce PGC to 2%, and will thus support future breeding programs to improve the grain quality of elite genetic material with high-yielding potentials.

Fine mapping of two grain chalkiness QTLs sensitive to high temperature in rice

BACKGROUND

Grain chalkiness is one of important factors affected rice grain quality. It is known that chalkiness is affected by the high temperature during the seed filling period. Although a larger of QTLs for chalkiness were reported across all 12 chromosomes, only a few of the QTLs were fine mapped or cloned up to now. Here, we fine map two QTLs for chalkiness in two single-segment substitution lines (SSSLs), 11-09 with substitution segment from O. sativa and HP67-11 with substitution segment from O. glaberrima.

RESULTS

The grain chalkiness of SSSLs 11-09 and HP67-11 was significantly lower than that in the recipient Huajingxian 74 (HJX74) in consecutive 8 cropping seasons. The regression correlation analysis showed that percentage of chalky grain (PCG) and percentage of chalky area (PCA) were significantly and positively correlated with percentage of grain chalkiness (PGC). Two QTLs for grain chalkiness were located on two chromosomes by substitution mapping. qPGC9 was mapped on chromosome 9 with an estimated interval of 345.6 kb. qPGC11 was located on chromosome 11 and delimited to a 432.1 kb interval in the O. sativa genome and a 332.9 kb interval in the O. glaberrima genome. qPGC11 is a QTL for grain chalkiness from O. glaberrima and was mapped in a new region of chromosome 11. The effect of two QTLs was incomplete dominance. The additive effects of two QTLs on chalkiness in second cropping season (SCS) were significantly greater than that in first cropping season (FCS).

CONCLUSIONS

qPGC11 is a new QTL for grain chalkiness. The two QTLs were fine mapped. The donor alleles of qPGC9 and qPGC11 were sensitive to the high temperature of FCS.

Bin-based genome-wide association studies reveal superior alleles for improvement of appearance quality using a 4-way MAGIC population in rice

•

4-way Multiparental population covered the limitations of the biparental structure.

•

The combination of SNP and bin-GWAS showed a powerful tool for QTL mapping.

•

qPGWC8.2 harbored a novel predicted gene for rice chalkiness quality.

Introduction

The multiparental population provides us the chance to identify superior alleles controlling a trait for genetic improvement. Genome wide association studies at bin level (bin-GWAS) are expected to be more power in QTL mapping than GWAS at SNP level (SNP-GWAS).

Objectives

This study is to estimate genetic effects of QTL conferring grain appearance quality in rice by SNP-GWAS and bin-GWAS, compare their power in QTL mapping and identify the superior alleles of all detected QTL from 4 parents for genetic improvement.

Methods

A 4-way MAGIC population and its four founders were cultivated in two environments to dissect the genetic basis of rice grain appearance quality. Both SNP-GWAS and bin-GWAS were conducted for QTL mapping. Multiple comparison among 4 parental bin/alleles was used to identify the superior alleles.

Results

A total of 16 and 20 QTL associated with grain appearance quality were identified by SNP- and bin-GWAS, respectively. A minor chalkiness QTL qPGWC8.2/qDEC8 was assigned to a 30-kb genomic region, in which OsMH_08T0121900 is the potential candidate gene because its encoded protein, glucan endo-1,3-beta-glucosidase precursor is involved in the starch and sucrose metabolism pathway. The superior parental alleles for GS3, GL3.1, GW5, GW7, and Chalk5 and two QTLs were almost carried by the high-quality parents Cypress and Yuejingsimiao (YJSM), while the poor-quality parent Guichao-2 (GC2) always carried the inferior alleles. The top five recombinant inbred lines with the highest quality of grain shape and chalkiness traits all carried gene combinations of superior alleles.

Conclusions

Both SNP- and bin-GWAS methods are encouraged for joint QTL mapping with MAGIC population. qPGWC8.2/qDEC8 is a novel candidate gene strongly associated with chalkiness. The superior alleles of GS3, GW5, GL3.1, GW7, Chalk5 and qPGWC8.2 were identified, and the pyramiding of these superior alleles is helpful to improve rice appearance quality.

Genome-wide association study reveals novel genetic loci contributing to cold tolerance at the germination stage in indica rice

Low temperature at the germination stage is one of the major abiotic stresses limiting rice (Oryza sativa L.) production, especially in regions where rice seeds are sown directly. However, few relevant genetic loci and genes have been identified. In this study, we report the phenotypic analysis of low temperature germination (LTG) in 200 indica rice varieties and a genome-wide association study (GWAS) of LTG in this collection using 161,657 high-quality SNPs, which were identified via genotyping-by-sequencing (GBS) of all the rice varieties. A total of 159 genetic loci were detected, and they were evenly distributed on all 12 chromosomes. Among them, 51 loci were detected more than twice; in particular, 23 loci were detected repeatedly in both the wet and dry seasons, and 569 genes were predicted in the 200-kb genomic region harbouring these 23 loci. Furthermore, 14,742 differentially expressed genes (DEGs) were identified using RNA sequencing. By integrating GWAS and RNA sequencing, 179 candidate DEGs were obtained. Sequence variation in the region of loci 95 was analyzed using 20 varieties with extreme phenotype. The polymorphisms of three DEGs (Os07g0585500, Os07g0585700, Os07g0585900) were associated with their phenotypes. Haplotype analysis of the three genes demonstrated that almost all the varieties with the same haplotype as japonica Nipponbare on the three DEGs showed high LTG ability. These findings provide valuable information for understanding the genetic control of LTG and performing molecular breeding with marker-assisted selection in indica rice.

Genome-Wide SNP Discovery and Mapping QTLs for Seed Iron and Zinc Concentrations in Chickpea (Cicer arietinum L.)

Biofortification through plant breeding is a cost-effective and sustainable approach towards addressing micronutrient malnutrition prevailing across the globe. Screening cultivars for micronutrient content and identification of quantitative trait loci (QTLs)/genes and markers help in the development of biofortified varieties in chickpea (Cicer arietinum L.). With the aim of identifying the genomic regions controlling seed Fe and Zn concentrations, the F2:3 population derived from a cross between MNK-1 and Annigeri 1 was genotyped using genotyping by sequencing approach and evaluated for Fe and Zn concentration. An intraspecific genetic linkage map comprising 839 single nucleotide polymorphisms (SNPs) spanning a total distance of 1,088.04 cM with an average marker density of 1.30 cM was constructed. By integrating the linkage map data with the phenotypic data of the F2:3 population, a total of 11 QTLs were detected for seed Fe concentration on CaLG03, CaLG04, and CaLG05, with phenotypic variation explained ranging from 7.2% (CaqFe3.4) to 13.4% (CaqFe4.2). For seed Zn concentration, eight QTLs were identified on CaLG04, CaLG05, and CaLG08. The QTLs individually explained phenotypic variations ranging between 5.7% (CaqZn8.1) and 13.7% (CaqZn4.3). Three QTLs for seed Fe and Zn concentrations (CaqFe4.4, CaqFe4.5, and CaqZn4.1) were colocated in the "QTL-hotspot" region on CaLG04 that harbors several drought tolerance-related QTLs. We identified genes in the QTL regions that encode iron-sulfur metabolism and zinc-dependent alcohol dehydrogenase activity on CaLG03, iron ion binding oxidoreductase on CaLG04, and zinc-induced facilitator-like protein and ZIP zinc/iron transport family protein on CaLG05. These genomic regions and the associated markers can be used in marker-assisted selection to increase seed Fe and Zn concentrations in agronomically superior chickpea varieties.

Dynamic transcriptome and metabolome analyses of two types of rice during the seed germination and young seedling growth stages

Background

Seed germination and young seedling growth are important agricultural traits for developing populations of both irrigated and directly seeded rice. Previous studies have focused on the identification of QTLs. However, there are few studies on the metabolome or transcriptome of germination and young seedling growth in rice.

Results

Here, an indica rice and a japonica rice were used as materials, and the transcripts and metabolites were detected during the germination and young seedling growth periods on a large scale by using RNA sequencing and a widely targeted metabolomics method, respectively. Fourteen shared transcripts and 15 shared metabolites that were continuously differentially expressed in the two materials were identified and may be essential for seed germination and young seedling growth. Enrichment analysis of differentially expressed genes in transcriptome expression profiles at different stages indicated that cell wall metabolism, lipid metabolism, nucleotide degradation, amino acid, etc., were enriched at 0–2 days, and most of the results are consistent with those of previous reports. Specifically, phenylpropanoid biosynthesis and glutathione metabolism were continuously enriched during the seed germination and young seedling growth stages. Next, KO enrichment analysis was conducted by using the differentially expressed genes of the two materials at 2, 3 and 4 days. Fourteen pathways were enriched. Additionally, 44 differentially expressed metabolites at 2, 3 and 4 days were identified. These metabolites may be responsible for the differences in germination and young seedling growth between the two materials. Further attention was focused on the ascorbate–glutathione pathway, and it was found that differences in ROS-scavenging abilities mediated by some APX, GPX and GST genes may be directly involved in mediating differences in the germination and young seedling growth speed of the two materials.

Conclusions

In summary, these results may enhance the understanding of the overall mechanism of seed germination and young seedling growth, and the outcome of this study is expected to facilitate rice breeding for direct seeding.

Genetic Analysis of Agronomic Traits and Grain Iron and Zinc Concentrations in a Doubled Haploid Population of Rice (Oryza sativa L.)

The development of micronutrient dense rice varieties with good agronomic traits is one of the sustainable and cost-effective approaches for reducing malnutrition. Identification of QTLs for high grain Fe and Zn, yield and yield components helps in precise and faster development of high Fe and Zn rice. We carried out a three-season evaluation using IR05F102 x IR69428 derived doubled-haploid population at IRRI. Inclusive composite interval mapping was carried out using SNP markers and Best Linear Unbiased Estimates of the phenotypic traits. A total of 23 QTLs were identified for eight agronomic traits and grain Fe and Zn concentration that explained 7.2 to 22.0% PV. A QTL by environment interaction analysis confirmed the stability of nine QTLs, including two QTLs for Zn on chromosomes 5 and 12. One epistatic interaction for plant height was significant with 28.4% PVE. Moreover, five QTLs were identified for Fe and Zn that harbor several candidate genes, e.g. OsZIP6 on QTL qZn_5.1. A number of QTLs were associated with a combination of greater yield and increased grain Zn levels. These results are useful for development of new rice varieties with good agronomic traits and high grain Zn using MAS, and identification of genetic resources with the novel QTLs for grain Zn.

QTLs identification for nitrogen and phosphorus uptake-related traits using ultra-high density SNP linkage

To understand the genetic basis of nitrogen and phosphorus uptake in the cultivated rice, quantitative trait loci (QTL) analysis for 7 nitrogen and phosphorus uptake-related traits including above-ground biomass (AGB), leaf colour value (SPAD) in heading stage, grain nitrogen concentration (GNC), grain nitrogen content of the plant, total nitrogen content (TNC), grain phosphorus concentration, total phosphorus content (TPC) were conducted using SNP markers in a F2 population derived from a cross between GH128 and W6827. A total of 21 QTLs for nitrogen and phosphorus uptake-related traits distributed in 16 regions along 6 chromosomes were detected using a high density genetic map consisting of 1582 bin markers, with QTLs maximum explaining 8.19% of the phenotypic variation. Nine QTLs (42.9% of total QTLs) were detected on chromosome 2. Among them, two QTL clusters including AGB, TNC, TPC and GNC were also detected in the region bin 140 and bin 146 on the chromosome 2. The distance between the two clusters was only 4.1 cM. The presence of QTL clusters has important significance and could be useful in molecular marker assisted breeding. These genomic regions might be deployed for the simultaneous improving the use efficiency of nitrogen and phosphorus in rice breeding.

Dissecting the genome-wide genetic variants of milling and appearance quality traits in rice

Higher head rice yield (HRY), which represents the proportion of intact grains that survive milling, and lower grain chalkiness (opacity) are key quality traits. We investigated the genetic basis of HRY and chalkiness in 320 diverse resequenced accessions of indica rice with integrated single- and multi-locus genome-wide association studies using 2.26 million single-nucleotide polymorphisms. We identified novel haplotypes that underly higher HRY on chromosomes 3, 6, 8, and 11, and that lower grain chalkiness in a fine-mapped region on chromosome 5. Whole-genome sequencing of 92 IRRI breeding lines was performed to identify the genetic variants of HRY and chalkiness. Rare and novel haplotypes were found for lowering chalkiness, but missing alleles hindered progress towards enhancing HRY in breeding material. The novel haplotypes that we identified have potential use in breeding programs aimed at improving these important traits in the rice crop.

Molecular insights into the regulation of rice kernel elongation

A large number of rice agronomic traits are complex, multi factorial and polygenic. As the mechanisms and genes determining grain size and yield are largely unknown, the identification of regulatory genes related to grain development remains a preeminent approach in rice genetic studies and breeding programs. Genes regulating cell proliferation and expansion in spikelet hulls and participating in endosperm development are the main controllers of rice kernel elongation and grain size. We review here and discuss recent findings on genes controlling rice grain size and the mechanisms, epialleles, epigenomic variation, and assessment of controlling genes using genome-editing tools relating to kernel elongation.

Identification of stable QTLs and candidate genes involved in anaerobic germination tolerance in rice via high-density genetic mapping and RNA-Seq

Background

Anaerobic germination tolerance is an important trait for direct-seeded rice varieties. Understanding the genetic basis of anaerobic germination is a key for breeding direct-seeded rice varieties.

Results

In this study, a recombinant inbred line (RIL) population derived from a cross between YZX and 02428 exhibited obvious coleoptile phenotypic differences. Mapping analysis using a high-density bin map indicated that a total of 25 loci were detected across two cropping seasons, including 10 previously detected loci and a total of 13 stable loci. Analysis of the 13 stable loci demonstrated that the more elite alleles that were pyramided in an individual, the higher the values of these traits were in the two cropping seasons. Furthermore, some anaerobic germination-tolerant recombinant inbred lines, namely G9, G10, G16, and G151, were identified. A total of 84 differentially expressed genes were obtained from the 13 stable loci via genome-wide expression analysis of the two parents at three key periods. Among them, Os06g0110200, Os07g0638300, Os07g0638400, Os09g0532900, Os09g0531701 and Os12g0539751 constitute the best candidates associated with anaerobic germination.

Conclusions

Both the anaerobic germination-tolerant recombinant inbred lines and the loci identified in this study will provide new genetic resources for improving the anaerobic germination tolerance of rice using molecular breeding strategies, as well as will broaden our understanding of the genetic control of germination tolerance under anaerobic conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5741-y) contains supplementary material, which is available to authorized users.

Applications of the CRISPR/Cas9 System for Rice Grain Quality Improvement: Perspectives and Opportunities

Grain quality improvement is a key target for rice breeders, along with yield. It is a multigenic trait that is simultaneously influenced by many factors. Over the past few decades, breeding for semi-dwarf cultivars and hybrids has significantly contributed to the attainment of high yield demands but reduced grain quality, which thus needs the attention of researchers. The availability of rice genome sequences has facilitated gene discovery, targeted mutagenesis, and revealed functional aspects of rice grain quality attributes. Some success has been achieved through the application of molecular markers to understand the genetic mechanisms for better rice grain quality; however, researchers have opted for novel strategies. Genomic alteration employing genome editing technologies (GETs) like clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) for reverse genetics has opened new avenues of research in the life sciences, including for rice grain quality improvement. Currently, CRISPR/Cas9 technology is widely used by researchers for genome editing to achieve the desired biological objectives, because of its simple targeting. Over the past few years many genes that are related to various aspects of rice grain quality have been successfully edited via CRISPR/Cas9 technology. Interestingly, studies on functional genomics at larger scales have become possible because of the availability of GETs. In this review, we discuss the progress made in rice by employing the CRISPR/Cas9 editing system and its eminent applications. We also elaborate possible future avenues of research with this system, and our understanding regarding the biological mechanism of rice grain quality improvement.

A Cluster of MYB Transcription Factors Regulates Anthocyanin Biosynthesis in Carrot (Daucus carota L.) Root and Petiole

Purple carrots can accumulate large quantities of anthocyanins in their roots and -in some genetic backgrounds- petioles, and therefore they represent an excellent dietary source of antioxidant phytonutrients. In a previous study, using linkage analysis in a carrot F₂ mapping population segregating for root and petiole anthocyanin pigmentation, we identified a region in chromosome 3 with co-localized QTL for all anthocyanin pigments of the carrot root, whereas petiole pigmentation segregated as a single dominant gene and mapped to one of these "root pigmentation" regions conditioning anthocyanin biosynthesis. In the present study, we performed fine mapping combined with gene expression analyses (RNA-Seq and RT-qPCR) to identify candidate genes controlling anthocyanin pigmentation in the carrot root and petiole. Fine mapping was performed in four carrot populations with different genetic backgrounds and patterns of pigmentation. The regions controlling root and petiole pigmentation in chromosome 3 were delimited to 541 and 535 kb, respectively. Genome wide prediction of transcription factor families known to regulate the anthocyanin biosynthetic pathway coupled with orthologous and phylogenetic analyses enabled the identification of a cluster of six MYB transcription factors, denominated DcMYB6 to DcMYB11, associated with the regulation of anthocyanin biosynthesis. No anthocyanin biosynthetic genes were present in this region. Comparative transcriptome analysis indicated that upregulation of DcMYB7 was always associated with anthocyanin pigmentation in both root and petiole tissues, whereas DcMYB11 was only upregulated with pigmentation in petioles. In the petiole, the level of expression of DcMYB11 was higher than DcMYB7. DcMYB6, a gene previously suggested as a key regulator of carrot anthocyanin biosynthesis, was not consistently associated with pigmentation in either tissue. These results strongly suggest that DcMYB7 is a candidate gene for root anthocyanin pigmentation in all the genetic backgrounds included in this study. DcMYB11 is a candidate gene for petiole pigmentation in all the purple carrot sources in this study. Since DcMYB7 is co-expressed with DcMYB11 in purple petioles, the latter gene may act also as a co-regulator of anthocyanin pigmentation in the petioles. This study provides linkage-mapping and functional evidence for the candidacy of these genes for the regulation of carrot anthocyanin biosynthesis.

Down-Regulation of SSSII-2 Gene Expression Results in Novel Low-Amylose Rice with Soft, Transparent Grains

Although soft rice, with low amylose content (AC), has high eating and cooking quality (ECQ), its appearance is poor due to the opaque endosperm. Here, a novel soft rice with low AC but a transparent appearance was generated by knocking-down the expression of SSSII-2, a gene encoding one isoform of soluble starch synthase (SSS). The physicochemical properties of the SSSII-2 RNAi rice are quite different from the control but more like the popular soft rice "Nanjing 46". The taste value assay further demonstrated that the ECQ of SSSII-2 RNAi rice was as high as "Nanjing 46", but only SSSII-2 RNAi rice retained the transparent endosperm under low moisture conditions. Further examination showed that the different morphologies and fine structures of the starch granules may contribute to the specific properties of SSSII-2 RNAi rice. Therefore, SSSII-2 has potential application in future high quality rice breeding programs.

Colocalization of QTLs for hull-cracked rice and grain size in elite rice varieties in Japan

The control of insects that consume cereal grains is important for the production and storage of grains. Hull-cracked rice, which has splits in the hull, becomes more susceptible to insects both in the paddy field and during storage. The development of varieties with a low frequency of hull-cracked rice is the most economical and effective strategy to avoid insect damage and the environmental risks from agricultural chemical entering rice grains. In this study, we identified that QTLs for the frequency of hull-cracked rice and for grain width are located on the same chromosome using recombinant inbred lines derived from a cross between the elite rice varieties in Hokkaido, Japan, which are from the same pedigree and are genetically closely related. These QTLs were detected close to different molecular markers, which were separated by 1,101,675 bp, on chromosome 5 in the reference Nipponbare genome. In addition, low coefficient values of the phenotype were found between hull-cracked rice and grain size. These results suggested that the ratio of hull-cracked rice is independent of grain size. Using these QTLs, new varieties with low hull-cracked rice could be developed regardless of grain size.

Genetic Dissection of qPCG1 for a Quantitative Trait Locus for Percentage of Chalky Grain in Rice (Oryza sativa L.)

Rice is a pivotal cereal crop that provides the staple food for more than half of the world's population. Along with improvements in the standard of living, people not only pay attention to the grain yield but also to the grain quality. Chalkiness is one of the most important index of grain quality. In this study, qPCG1, a QTL for percentage of chalky grain, was mapped in an interval with a physical distance about 139 kb on chromosome 1 by residual heterozygous line (RHL) method. qPCG1 was incomplete dominant and the additive effect plays a major role and explained 6.8-21.9% of phenotypic variance within the heterogeneous region on chromosome 1. The effect of allele from Zhonghui9308 was decreasing the percentage of chalky grains (PCG). Microscope observation results indicated that there are great differences in the shape, structure and arrangement of starch granule between the chalky part and transparent part. Analysis of starch physicochemical properties showed that the total starch content, amylose content and chain length distribution of amylopectin changed while the protein contents were not apparently affected with the changed chalkiness. qPCG1 had little influence on main agronomic traits and it might be useful in rice breeding for it did not bring negative effect on grain yield while reducing the chalkiness.

Selection of Transcripts Affecting Initial Growth Rate of Rice Backcrossed Inbred Lines Using RNA Sequencing Data

Seedling growth is an important factor for direct seeding of rice. However, the genetic and transcriptomic factors involved in this process are largely unknown. In this study, transcripts affecting shoot weight were identified in rice (Oryza sativa L.) using RNA sequencing (RNA-Seq) data from 20 backcrossed inbred lines (BILs) and their parental cultivars. The selection frequency of the genes for the regression model was determined using repeated analysis of random subsets of the transcriptome. The qLTG3-1gene, controlling low-temperature germinability, and short grain 1 gene (SG1), known to decrease organ elongation, showed high frequency. The quantitative trait loci (QTLs) analysis performed for BILs revealed that qLTG3-1 was included in the QTLs for shoot weight but SG1 was not. No nucleotide polymorphisms were found in the coding region of SG1 in either of the parental cultivars. Quantitative real-time PCR showed that SG1 expression was negatively correlated with shoot weight for all 104 BILs analyzed in this study. Expression QTL (eQTLs) analysis showed an eQTL for SG1 expression located in the same region as the QTL for shoot weight. However, no eQTLs were detected on the same chromosome as SG1, suggesting that nucleotide polymorphisms around the gene do not affect its expression in analyzed growth stage. Overall, these results indicate that RNA-Seq is a useful tool for identifying transcripts that can be related to seedling growth rate.

QTL mapping using an ultra-high-density SNP map reveals a major locus for grain yield in an elite rice restorer R998

To dissect the genetic basis of yield formation in restorer line of hybrid rice, we conducted QTL analysis for 6 yield traits including panicles per plant (PPP), grains per panicle (GPP), grain yield per plant (GY), thousand-grain weight (TGW), above-ground biomass (AGB), and harvest index (HI) using SNP markers in a recombinant inbred lines (RILs) population derived from a cross between a tropical japonica inbred Francis and an elite indica restorer Guanghui 998 (R998). A total of 26 QTLs were detected using a high density genetic map consisting of 3016 bin markers. Nineteen out of the 26 QTL alleles from R998 had a beneficial effect on yield traits. Most of the QTLs were co-located with previously reported rice QTLs. qAGB6 and qHI9, controlling AGB and HI respectively, were detected as novel QTLs. Four QTLs for GY were repeatedly detected across two years, with all the beneficial alleles from R998. Notably, qGY8 explained over 20% of the yield variance in both years. Moreover, qGY8 together with qTGW8 and qHI8 formed a QTL cluster. Markers tightly linked with qGY8 were developed. Cloning of qGY8 will facilitate its further exploitation in high-yield breeding.

Detection and validation of QTLs for milky-white grains caused by high temperature during the ripening period in Japonica rice

The occurrence of chalky rice (Oryza sativa L.) grains caused by high temperature is a serious problem in rice production. Of the several kinds of chalky grains, milky-white grains are not well analyzed. The milky-white rice grain phenomenon is caused by genetic factors as well as environmental and nutritional conditions. To analyze the genetic control system for rice grain quality, we raised recombinant inbred lines from progeny produced from 'Tsukushiroman' (high temperature sensitive) and 'Chikushi 52' (high temperature tolerant) cultivars. Quantitative trait locus (QTL) analysis revealed that the QTL on chromosome 4, linked to the simple sequence repeat marker RM16424, contributed substantially to the occurrence of milky-white grains, as it was detected over two experimental years. To validate the effect of the QTL, we developed near isogenic lines that have the 'Chikushi 52' segment on the short arm of chromosome 4 in the 'Tsukushiroman' genetic background, and that had a lower milky-white grain ratio than that of 'Tsukushiroman' when exposed to high temperatures during the ripening period. These results suggest that the 'Chikushi 52' allele on chromosome 4 suppresses the occurrence of milky-white grains and improves rice grain quality under heat stress during the grain ripening period.

A Natural Allele of a Transcription Factor in Rice Confers Broad-Spectrum Blast Resistance

Rice feeds half the world's population, and rice blast is often a destructive disease that results in significant crop loss. Non-race-specific resistance has been more effective in controlling crop diseases than race-specific resistance because of its broad spectrum and durability. Through a genome-wide association study, we report the identification of a natural allele of a C₂H₂-type transcription factor in rice that confers non-race-specific resistance to blast. A survey of 3,000 sequenced rice genomes reveals that this allele exists in 10% of rice, suggesting that this favorable trait has been selected through breeding. This allele causes a single nucleotide change in the promoter of the bsr-d1 gene, which results in reduced expression of the gene through the binding of the repressive MYB transcription factor and, consequently, an inhibition of H₂O₂ degradation and enhanced disease resistance. Our discovery highlights this novel allele as a strategy for breeding durable resistance in rice.