RiceVarMap: a comprehensive database of rice genomic variations

A combination of upstream alleles involved in rice heading hastens natural long-day responses

BACKGROUND

The female parental line Jinbuol (JBO, early heading) and two recombinant isogenic lines, JSRIL1 and JSRIL2, have been shown to flower 44, 34 and 16 days earlier, respectively, than the male parental line Samgwang (SG, late heading) in paddy fields.

OBJECTIVE

To explore how photoperiodicity-related genes are involved in differential heading among these lines.

METHODS

Deep sequencing was conducted for these lines, photoperiodicity-related genes (71) were categorized, and qRT-PCR was performed for some key genes.

RESULTS

Deep sequencing revealed a nearly even contribution of parental groups, with 48.5% and 45% of the chromosomes in JSRIL1 and JSRIL2, respectively, inherited from the female parent JBO; however, Chr6 contained the most biased parental contribution, with 99.4% inherited from the female parent. The variation in single-nucleotide polymorphisms (SNPs) among many known flower-inducing genes, including rice GIGANTEA (OsGI); grain number, plant height and heading date 7 (Ghd7); and EARLY HEADING DATE 1 (Ehd1), was minimal. In the JSRILs, HEADING DATE 1 (Hd1) and VERNALIZATION INSENSITIVE 3-LIKE 1 (OsVIL2) originated from JBO, whereas FLAVIN-BINDING, KELCH REPEAT, F BOX 1 (OsFKF1) originated from SG. Interestingly, HEN1 suppressor 1 (OsHESO1) originated from SG in JSRIL1 and JBO in JSRIL2. RNA sequencing and qRT‒PCR analyses of plants at the floral meristem stage revealed that transcriptional regulation through chromosomal restructuring and posttranscriptional regulation might control minute gene regulation, resulting in delayed heading in JSRILs.

CONCLUSION

Our gene expression and SNP analyses of elite recombinant isogenic lines could be helpful in understanding how photoperiodicity-related genes in rice are modulated.

The cation/H+ exchanger OsCAX2 is involved in cadmium uptake and contributes to differential grain cadmium accumulation between Indica and Japonica rice

Rice is a major source of dietary cadmium (Cd), a toxic heavy metal that poses serious threat to human health. How rice takes up and accumulates Cd is not fully understood. Here, we characterize the function of a cation/H+ exchanger, OsCAX2, in Cd uptake in roots and Cd accumulation in shoots and grains. OsCAX2 exhibited Cd and calcium (Ca) transport activities when was heterologously expressed in yeast. OsCAX2 was mainly expressed in roots, particularly in lateral roots, and in the exodermis and endodermis of primary roots. OsCAX2 is localized at the plasma membrane. Knockout of OsCAX2 significantly decreased Cd uptake in roots and Cd accumulation in shoots and grains. Knockout of OsCAX2 also decreased the Ca concentration in roots, but not in shoots or grains. Surprisingly, overexpression of OsCAX2 also resulted in a significant decrease in the Cd concentrations in roots and shoots. We further reveal that the variation in the coding sequence of OsCAX2 contributes to differential grain Cd accumulation between two major rice subspecies, Indica and Japonica. Our results demonstrate that OsCAX2 functions in Cd/Ca uptake in roots and could be a useful target for breeding or genetic engineering low Cd rice varieties.

Natural variation in OsMADS1 transcript splicing affects rice grain thickness and quality by influencing monosaccharide loading to the endosperm

Grain size, which encompasses grain length, width, and thickness, is a critical determinant of both grain weight and quality in rice. Despite the extensive regulatory networks known to determine grain length and width, the pathway(s) that regulate grain thickness remain to be clarified. Here, we present the map-based cloning and characterization of qGT3, a major quantitative trait locus for grain thickness in rice that encodes the MADS-domain transcription factor OsMADS1. Our findings demonstrate that OsMADS1 regulates grain thickness by affecting sugar delivery during grain filling, and we show that OsMADS1 modulates expression of the downstream monosaccharide transporter gene MST4. A natural variant leads to alternative splicing and thus to a truncated OsMADS1 protein with attenuated transcriptional repressor activity. The truncated OsMADS1 protein results in increased expression of MST4, leading to enhanced loading of monosaccharides into the developing endosperm and thereby increasing grain thickness and improving grain quality. In addition, our results reveal that NF-YB1 and NF-YC12 interact directly with OsMADS1, acting as cofactors to enhance its transcriptional activity toward MST4. Collectively, these findings reveal a novel molecular mechanism underlying grain thickness regulation that is controlled by the OsMADS1-NF-YB1-YC12 complex and has great potential for synergistic improvement of grain yield and quality in rice.

Thousands of trait-specific KASP markers designed for diverse breeding applications in rice (Oryza sativa)

This study aimed to broaden applicability of KASP for Oryza sativa across diverse genotypes through incorporation of ambiguous (degenerate) bases into their primer designs and to validate 4,000 of them for genotyping applications. A bioinformatics pipeline was used to compare 129 rice genomes from 89 countries with the indica reference genome R498 and generate ∼1.6 million KASP designs for the more common variants between R498 and the other genomes. Of the designs, 98,238 were for predicted functional markers. Up to 5 KASP each for 1,024 breeder-selected loci were assayed in a panel of 178 diverse rice varieties, generating 3,366 validated KASP. The 84% success rate was within the normal range for KASP demonstrating that the ambiguous bases do not compromise efficacy. The 3,366-trait-specific marker panel was applied for population structure analysis in the diversity panel and resolved them into 4 expected groups. Target variations in 13 genomes used for designs were compared with the corresponding KASP genotypes in different accessions of the same 13 varieties in the diversity panel. There was agreement for 79% or more markers in 12 varieties; 10 having agreement >88%. One variety, a selection from a landrace, had only 46.5% marker agreement. Breeders can search for the validated KASP and more than a million so-far untested designs in three reference genomes (including Niponbare MSU7) with a search tool, that includes designs in proximity to previously published microsatellite markers, and retrieve target variations for 129 rice genomes plus their genomic locations with ±25 bp flanking sequences.

The RAD6-like Ubiquitin Conjugase Gene OsUBC7 Has a Positive Role in the Early Cold Stress Tolerance Response of Rice

BACKGROUND/OBJECTIVES

Cold stress poses a significant threat to Asian rice cultivation, disrupting important physiological processes crucial for seedling establishment and overall plant growth. It is, thus, crucial to elucidate genetic pathways involved in cold stress tolerance response mechanisms.

METHODS

We mapped OsUBC7, a Radiation-sensitive 6 (RAD6)-type homolog of rice, to a low-temperature seedling survivability (LTSS) QTL and used genomics, molecular genetics, and physiological assays to assess its role in plant resilience against low-temperature stress.

RESULTS

OsUBC7 is cold responsive and has higher expression levels in cold-tolerant japonica than cold-sensitive indica. Overexpression of OsUBC7 enhances LTSS of indica and freezing tolerance of Arabidopsis, increases levels of soluble sugars and chlorophyll A, boosts leaf development after cold exposure, and increases leaf cell numbers and plants size, but it does not affect membrane stability after cold stress exposure. Additionally, OsUBC7 has a positive role for germinability in the presence of salt and for flowering and yield-related traits. The OsUBC7 protein physically interacts with the developmental stage-specific and histone-modifying E3 ligases OsRFPH2-12 and OsHUB1/2, respectively, and potential target genes such as cell cycle dependent kinases were identified.

CONCLUSIONS

OsUBC7 might contribute to cold resilience by regulating sugar metabolism to provide energy for promoting cellular homeostasis restoration after cold stress exposure via new cell growth, particularly in leaf cells crucial for photosynthesis and metabolic activity, possibly by interacting with cell cycle regulating proteins. Overall, the present study suggests that OsUBC7 may be involved in plant development, reproduction, and stress adaptation, and contributes to a deeper understanding of rice plant cold stress tolerance response mechanisms. OsUBC7 may be a promising candidate for improving crop productivity and resilience to stressful environments.

Genetic Analysis and Fine Mapping of QTLs for Stigma Exsertion Rate in Rice

Stigma exsertion rate (SER) is a crucial trait that influences the seed production of hybrid rice by determining the outcrossing ability of male sterile lines (MSLs). However, the molecular genetic mechanisms underlying SER are still poorly understood. In this study, we identified 14 quantitative trait loci (QTLs) using a recombinant inbred line (RIL) population derived from B805D-MR-16-8-3 (B805D) and Hua6S. Two major QTLs, qSE1 and qSE9, were validated for their effects in the residual heterozygous line (RHL) background. The RHL carrying homozygous qSE1 region from Hua6S increased dual stigma exsertion rate (DSE) by 14.67% and 15.04%, and increased total stigma exsertion rate (TSE) by 11.73% and 13.04%, in F10 and F11 progeny, respectively. Conversely, the RHL carrying homozygous qSE9 region from B805D showed a substantial increase of 22.72% and 14.45% in single stigma exsertion rate (SSE), an increase of 13.46% and 8.30% in TSE, and an increase in percentage of spikelets with exserted stigma (PSE) by 24.82% and 15.57%, respectively, in F10 and F11 progeny. Furthermore, examination of floral organ traits revealed that both the Hua6S allele of qSE1 and the B805D allele of qSE9 increased pistil size to improve SER, but they had contrasting effects on spikelet shape. Subsequently, qSE1 and qSE9 were fine-mapped to intervals of 246.5 kb and 341.4 kb, respectively. A combination of sequencing, expression and haplotype analysis revealed that a single nucleotide variation (T to C) in the 5'UTR region of LOC_Os01g72020 (OsBOP1) was likely to be the functional variation for qSE1. Collectively, our work has laid a foundation for cloning the genes responsible for SER, and demonstrated that the Hua6S allele of qSE1 and the B805D allele of qSE9 can effectively increase SER, which could make important contributions to the genetic improvement of MSLs aimed at improving hybrid seed production.

Volatilome-based GWAS identifies OsWRKY19 and OsNAC021 as key regulators of rice aroma

Aromatic rice is globally favored for its distinctive scent, which not only increases its nutritional value but also enhances its economic importance. However, apart from 2-acetyl-1-pyrroline (2-AP), the metabolic basis of aroma remains to be clarified, and the genetic basis of the accumulation of fragrance metabolites is largely unknown. In this study, we revealed 2-AP and fatty acid-derived volatiles (FAVs) as key contributors to rice aroma by combining aroma rating with molecular docking. Using a volatilome-based genome-wide association study, we identified two regulatory genes that determine the natural variation of these fragrance metabolites. Genetic and molecular analyses showed that OsWRKY19 not only enhances fragrance by negatively regulating OsBADH2 but also improves agricultural traits in rice. Furthermore, we revealed that OsNAC021 negatively regulates FAV contents via the lipoxygenase pathway, and its knockout resulted in over-accumulation of grain FAVs without a yield penalty. Collectively, our study not only identifies two key regulators of rice aroma but also provides a compelling example about how to deciphering the genetic regulatory mechanisms that underlie rice fragrance, thereby paving the way for the creation of aromatic rice varieties.

OsIAA23 Promotes Heading by Directly Downregulating Ghd7 in rice

Ghd7 is a central regulator to multiple growth and development processes in rice. While it is not clear how Ghd7 is regulated by upstream factors. To identify its upstream regulator, the truncated Ghd7 promoter fragments were used to screen cis elements binding to rice total nuclear proteins. Electrophoretic mobility shift assays screened one truncated fragment f3 binding to the proteins. Subsequently, the fragment f3 was employed to screen a yeast one-hybrid library, and a transcription factor OsIAA23 was screened as a direct upstream regulator of Ghd7. Dual-luciferase transient assay demonstrated the transcriptional repression effect of OsIAA23 on the activity of Ghd7, and the location of the cis elements binding to OsIAA23 in the region 1264 to 1255 bp upstream of ATG. Genetic analysis between the wild type Ghd7-OsIAA23 and single/double mutants further verified that OsIAA23 downregulated Ghd7 expression and led to a delayed heading under long day conditions. Moreover, natural variations in fragment f3 were associated with heading and geographic distribution in rice. This study sheds light on the direct regulatory mechanism of OsIAA23 on Ghd7, which enriches the understanding of the Ghd7 involved flowering regulatory network in rice.

EARLY FLOWERING3-1 represses Grain number, plant height, and heading date7 to promote ABC1 REPRESSOR1 and regulate nitrogen uptake in rice

The extensive use of nitrogen fertilizer boosts rice (Oryza sativa) production but also harms ecosystems. Therefore, enhancing crop nitrogen use efficiency is crucial. Here, we performed map-based cloning and identified the EARLY FLOWERING3 (ELF3) like protein-encoding gene OsELF3-1, which confers enhanced nitrogen uptake in rice. OsELF3-1 forms a ternary complex (OsEC) with OsELF4s and OsLUX, the putative orthologs of ELF4 and LUX ARRHYTHMO (LUX) in Arabidopsis (Arabidopsis thaliana), respectively. OsEC directly binds to the promoter of Grain number, plant height, and heading date7 (Ghd7) and represses its expression. Ghd7 encodes a transcription factor that has major effects on multiple agronomic traits. Ghd7 is also a transcriptional repressor and directly suppresses the expression of ABC1 REPRESSOR1 (ARE1), a negative regulator of nitrogen use efficiency. Therefore, targeting the OsEC-Ghd7-ARE1 module offers an approach to enhance nitrogen uptake, presenting promising avenues for sustainable agriculture.

Genome-Wide Association-Based Identification of Alleles, Genes and Haplotypes Influencing Yield in Rice (Oryza sativa L.) Under Low-Phosphorus Acidic Lowland Soils

Rice provides poor yields in acidic soils due to several nutrient deficiencies and metal toxicities. The low availability of phosphorus (P) in acidic soils offers a natural condition for screening genotypes for grain yield and phosphorus utilization efficiency (PUE). The objective of this study was to phenotype a subset of indica rice accessions from 3000 Rice Genome Project (3K-RGP) under acidic soils and find associated genes and alleles. A panel of 234 genotypes, along with checks, were grown under low-input acidic soils for two consecutive seasons, followed by a low-P-based hydroponic screening experiment. The heritability of the agro-morphological traits was high across seasons, and Ward's clustering method identified 46 genotypes that can be used as low-P-tolerant donors in acidic soil conditions. Genotypes ARC10145, RPA5929, and K1559-4, with a higher grain yield than checks, were identified. Over 29 million SNPs were retrieved from the Rice SNP-Seek database, and after quality control, they were utilized for a genome-wide association study (GWAS) with seventeen traits. Ten quantitative trait nucleotides (QTNs) for three yield traits and five QTNs for PUE were identified. A set of 34 candidate genes for yield-related traits was also identified. An association study using this indica panel for an already reported 1.84 Mbp region on chromosome 2 identified genes Os02g09840 and Os02g08420 for yield and PUE, respectively. A haplotype analysis for the candidate genes identified favorable allelic combinations. Donors carrying the superior haplotypic combinations for the identified genes could be exploited in future breeding programs.

A chloroplast localized heavy metal-associated domain containing protein regulates grain calcium accumulation in rice

Calcium (Ca) is an essential mineral nutrient and plays a crucial signaling role in all living organisms. Increasing Ca content in staple foods such as rice is vital for improving Ca nutrition of humans. Here we map a quantitative trait locus that controls Ca concentration in rice grains and identify the causal gene as GCSC1 (Grain Ca and Sr Concentrations 1), which encodes a chloroplast vesicle localized homo-oligomeric protein. GCSC1 exhibits Ca2+ transport activity in heterologous assays in yeast and Xenopus laevis oocytes and is involved in the efflux of Ca2+ from the chloroplast to the cytosol. Knockout of GCSC1 results in increased chloroplast Ca concentration, lower stomatal conductance in leaves and enhanced Ca allocation to grains. Natural variation in grain Ca concentration is attributed to the variable expression of GCSC1 resulting from its promoter sequence variation. Our study identifies a chloroplast localized heavy metal-associated domain containing protein that regulates chloroplast Ca2+ efflux and provides a way to biofortify Ca in rice to benefit human nutrition.

Genome-wide association study reveals the genetic basis of rice resistance to three herbicides

Crop resistance to herbicides is crucial for agricultural productivity and sustainability amidst escalating challenges of weed resistance. Uncovering herbicide resistant genes is particularly important for rice production. In this study, we tested the resistance to three commonly used herbicides: glufosinate, glyphosate and mesotrione of 421 diverse rice cultivars and employed genome-wide association studies (GWAS) to unravel the genetic underpinnings of resistance to these three herbicides in rice. We discovered that cultivated rice exhibited rich variation in resistance to the three herbicides, and the differences among subpopulations were significant. Six identified associations harboring candidate genes for resistance to these herbicides were significant. Among them, RGlu6 and RGly8 were the major QTL for resistance to glufosinate and glyphosate, respectively. The favorable alleles of RGlu6 and RGly8 were primarily present in japonica cultivars that originated from Europe, highlighting the geographic and genetic diversity of herbicide resistance and emphasizing the localized selection pressures in European rice varieties. Moreover, our findings might suggest that traditional target genes may not contain tolerant alleles in nature, and alternative mechanisms with novel loci associated with resistance may work. By mapping the genes for herbicide resistance, our results may help develop new strategies to combat the dual challenges on effective weed management and herbicide sustainability.

The OsZIP2 transporter is involved in root-to-shoot translocation and intervascular transfer of cadmium in rice

Cadmium (Cd) is a toxic metal that poses serious threats to human health. Rice is a major source of dietary Cd but how rice plants transport Cd to the grain is not fully understood. Here, we characterize the function of the ZIP (ZRT, IRT-like protein) family protein, OsZIP2, in the root-to-shoot translocation of Cd and intervascular transfer of Cd in nodes. OsZIP2 is localized at the plasma membrane and exhibited Cd2+ transport activity when heterologously expressed in yeast. OsZIP2 is strongly expressed in xylem parenchyma cells in roots and in enlarged vascular bundles in nodes. Knockout of OsZIP2 significantly enhanced root-to-shoot translocation of Cd and alleviated the inhibition of root elongation by excess Cd stress; whereas overexpression of OsZIP2 decreased Cd translocation to shoots and resulted in Cd sensitivity. Knockout of OsZIP2 increased Cd allocation to the flag leaf but decreased Cd allocation to the panicle and grain. We further reveal that the variation of OsZIP2 expression level contributes to grain Cd concentration among rice germplasms. Our results demonstrate that OsZIP2 functions in root-to-shoot translocation of Cd in roots and intervascular transfer of Cd in nodes, which can be used for breeding low Cd rice varieties.

Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many quantitative trait loci (QTLs) and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that was associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five single nucleotide polymophysim (SNP) variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity toward LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice and suggest the potential of WLGhap.B in improving rice yield.

Haplotype-based breeding: A new insight in crop improvement

Haplotype-based breeding (HBB) is one of the cutting-edge technologies in the realm of crop improvement due to the increasing availability of Single Nucleotide Polymorphisms identified by Next Generation Sequencing technologies. The complexity of the data can be decreased with fewer statistical tests and a lower probability of spurious associations by combining thousands of SNPs into a few hundred haplotype blocks. The presence of strong genomic regions in breeding lines of most crop species facilitates the use of haplotypes to improve the efficiency of genomic and marker-assisted selection. Haplotype-based breeding as a Genomic Assisted Breeding (GAB) approach harnesses the genome sequence data to pinpoint the allelic variation used to hasten the breeding cycle and circumvent the challenges associated with linkage drag. This review article demonstrates ways to identify candidate genes, superior haplotype identification, haplo-pheno analysis, and haplotype-based marker-assisted selection. The crop improvement strategies that utilize superior haplotypes will hasten the breeding progress to safeguard global food security.

Dissecting the genetic basis of UV-B responsive metabolites in rice

BACKGROUND

UV-B, an important environmental factor, has been shown to affect the yield and quality of rice (Oryza sativa) worldwide. However, the molecular mechanisms underlying the response to UV-B stress remain elusive in rice.

RESULTS

We perform comprehensive metabolic profiling of leaves from 160 diverse rice accessions under UV-B and normal light conditions using a widely targeted metabolomics approach. Our results reveal substantial differences in metabolite accumulation between the two major rice subspecies indica and japonica, especially after UV-B treatment, implying the possible role and mechanism of metabolome changes in subspecies differentiation and the stress response. We next conduct a transcriptome analysis from four representative rice varieties under UV-B stress, revealing genes from amino acid and flavonoid pathways involved in the UV-B response. We further perform a metabolite-based genome-wide association study (mGWAS), which reveals 3307 distinct loci under UV-B stress. Identification and functional validation of candidate genes show that OsMYB44 regulates tryptamine accumulation to mediate UV-B tolerance, while OsUVR8 interacts with OsMYB110 to promote flavonoid accumulation and UV-B tolerance in a coordinated manner. Additionally, haplotype analysis suggests that natural variation of OsUVR8groupA contributes to UV-B resistance in rice.

CONCLUSIONS

Our study reveals the complex biochemical and genetic foundations that govern the metabolite dynamics underlying the response, tolerance, and adaptive strategies of rice to UV-B stress. These findings provide new insights into the biochemical and genetic basis of the metabolome underlying the crop response, tolerance, and adaptation to UV-B stress.

Variation in a Poaceae-conserved fatty acid metabolic gene cluster controls rice yield by regulating male fertility

A wide variety of metabolic gene clusters exist in eukaryotic genomes, but fatty acid metabolic gene clusters have not been discovered. Here, combining with metabolic and phenotypic genome-wide association studies, we identify a major locus containing a six-gene fatty acid metabolic gene cluster on chromosome 3 (FGC3) that controls the cutin monomer hydroxymonoacylglycerols (HMGs) contents and rice yield, possibly through variation in the transcription of FGC3 members. We show that HMGs are sequentially synthesized in the endoplasmic reticulum by OsFAR2, OsKCS11, OsGPAT6, OsCYP704B2 and subsequently transported to the apoplast by OsABCG22 and OsLTPL82. Mutation of FGC3 members reduces HMGs, leading to defective male reproductive development and a significant decrease in yield. OsMADS6 and OsMADS17 directly regulate FGC3 and thus influence male reproduction and yield. FGC3 is conserved in Poaceae and likely formed prior to the divergence of Pharus latifolius. The eukaryotic fatty acid and plant primary metabolic gene cluster we identified show a significant impact on the origin and evolution of Poaceae and has potential for application in hybrid crop breeding.

Stress-induced nuclear translocation of ONAC023 improves drought and heat tolerance through multiple processes in rice

Drought and heat are major abiotic stresses frequently coinciding to threaten rice production. Despite hundreds of stress-related genes being identified, only a few have been confirmed to confer resistance to multiple stresses in crops. Here we report ONAC023, a hub stress regulator that integrates the regulations of both drought and heat tolerance in rice. ONAC023 positively regulates drought and heat tolerance at both seedling and reproductive stages. Notably, the functioning of ONAC023 is obliterated without stress treatment and can be triggered by drought and heat stresses at two layers. The expression of ONAC023 is induced in response to stress stimuli. We show that overexpressed ONAC23 is translocated to the nucleus under stress and evidence from protoplasts suggests that the dephosphorylation of the remorin protein OSREM1.5 can promote this translocation. Under drought or heat stress, the nuclear ONAC023 can target and promote the expression of diverse genes, such as OsPIP2;7, PGL3, OsFKBP20-1b, and OsSF3B1, which are involved in various processes including water transport, reactive oxygen species homeostasis, and alternative splicing. These results manifest that ONAC023 is fine-tuned to positively regulate drought and heat tolerance through the integration of multiple stress-responsive processes. Our findings provide not only an underlying connection between drought and heat responses, but also a promising candidate for engineering multi-stress-resilient rice.

Identification of hub genes involved in gibberellin-regulated elongation of coleoptiles of rice seeds germinating under submerged conditions

Rapid elongation of coleoptiles from rice seeds to reach the water surface enables plants to survive submergence stress and therefore plays a crucial role in allowing direct seeding in rice cultivation. Gibberellin (GA) positively influences growth in rice, but the molecular mechanisms underlying its regulation of coleoptile elongation under submerged conditions remain unclear. In this study, we performed a weighted gene co-expression network analysis to conduct a preliminarily examination of the mechanisms. Four key modules were identified with high correlations to the GA regulation of submergence tolerance. The genes within these modules were mainly involved in the Golgi apparatus and carbohydrate metabolic pathways, suggesting their involvement in enhancing submergence tolerance. Further analysis of natural variation revealed that the specific hub genes Os03g0337900, Os03g0355600, and Os07g0638400 exhibited strong correlations with subspecies divergence of the coleoptile elongation phenotype. Consistent with this analysis, mutation of Os07g0638400 resulted in a lower germination potential and a stronger inhibition of coleoptile elongation under submerged conditions. The hub genes identified in this study provide new insights into the molecular mechanisms underlying GA-dependent tolerance to submergence stress in rice, and a potential basis for future modification of rice germplasm to allow for direct seeding.

A natural variation in OsDSK2a modulates plant growth and salt tolerance through phosphorylation by SnRK1A in rice

Plants grow rapidly for maximal production under optimal conditions; however, they adopt a slower growth strategy to maintain survival when facing environmental stresses. As salt stress restricts crop architecture and grain yield, identifying genetic variations associated with growth and yield responses to salinity is critical for breeding optimal crop varieties. OsDSK2a is a pivotal modulator of plant growth and salt tolerance via the modulation of gibberellic acid (GA) metabolism; however, its regulation remains unclear. Here, we showed that OsDSK2a can be phosphorylated at the second amino acid (S2) to maintain its stability. The gene-edited mutant osdsk2aS2G showed decreased plant height and enhanced salt tolerance. SnRK1A modulated OsDSK2a-S2 phosphorylation and played a substantial role in GA metabolism. Genetic analysis indicated that SnRK1A functions upstream of OsDSK2a and affects plant growth and salt tolerance. Moreover, SnRK1A activity was suppressed under salt stress, resulting in decreased phosphorylation and abundance of OsDSK2a. Thus, SnRK1A preserves the stability of OsDSK2a to maintain plant growth under normal conditions, and reduces the abundance of OsDSK2a to limit growth under salt stress. Haplotype analysis using 3 K-RG data identified a natural variation in OsDSK2a-S2. The allele of OsDSK2a-G downregulates plant height and improves salt-inhibited grain yield. Thus, our findings revealed a new mechanism for OsDSK2a stability and provided a valuable target for crop breeding to overcome yield limitations under salinity stress.

The OsCLV2s-OsCRN1 co-receptor regulates grain shape in rice

The highly conserved CLV-WUS negative feedback pathway plays a decisive role in regulating stem cell maintenance in shoot and floral meristems in higher plants, including Arabidopsis, rice, maize, and tomato. Here, we find significant natural variations in the OsCLV2c, OsCLV2d, and OsCRN1 loci in a genome-wide association study of grain shape in rice. OsCLV2a, OsCLV2c, OsCLV2d, and OsCRN1 negatively regulate grain length-width ratio and show distinctive geographical distribution, indica-japonica differentiation, and artificial selection signatures. Notably, OsCLV2a and OsCRN1 interact biochemically and genetically, suggesting that the two components function in a complex to regulate grain shape of rice. Furthermore, the genetic contributions of the haplotypes combining OsCLV2a, OsCLV2c, and OsCRN1 are significantly higher than those of each single gene alone in controlling key yield traits. These findings identify two groups of receptor-like kinases that may function as distinct co-receptors to control grain size in rice, thereby revealing a previously unrecognized role of the CLV class genes in regulating seed development and proposing a framework to understand the molecular mechanisms of the CLV-WUS pathway in rice and other crops.

A genome-wide association study of panicle blast resistance to Magnaporthe oryzae in rice

Rice blast, caused by Magnaporthe oryzae (M. oryzae), is one of the most serious diseases worldwide. Developing blast-resistant rice varieties is an effective strategy to control the spread of rice blast and reduce the reliance on chemical pesticides. In this study, 477 sequenced rice germplasms from 48 countries were inoculated and assessed at the booting stage. We found that 23 germplasms exhibited high panicle blast resistance against M. oryzae. Genome-wide association analysis (GWAS) identified 43 quantitative trait loci (QTLs) significantly associated (P < 1.0 × 10-4) with resistance to rice panicle blast. These QTL intervals encompass four genes (OsAKT1, OsRACK1A, Bsr-k1 and Pi25/Pid3) previously reported to contribute to rice blast resistance. We selected QTLs with -Log10 (P-value) greater than 6.0 or those detected in two-year replicates, amounting to 12 QTLs, for further candidate gene analysis. Three blast resistance candidate genes (Os06g0316800, Os06g0320000, Pi25/Pid3) were identified based on significant single nucleotide polymorphisms (SNP) distributions within annotated gene sequences across these 12 QTLs and the differential expression levels among blast-resistant varieties after 72 h of inoculation. Os06g0316800 encodes a glycine-rich protein, OsGrp6, an important component of plant cell walls involved in cellular stress responses and signaling. Os06g0320000 encodes a protein with unknown function (DUF953), part of the thioredoxin-like family, which is crucial for maintaining reactive oxygen species (ROS) homeostasis in vivo, named as OsTrxl1. Lastly, Pi25/Pid3 encodes a disease resistance protein, underscoring its potential importance in plant biology. By analyzing the haplotypes of these three genes, we identified favorable haplotypes for blast resistance, providing valuable genetic resources for future rice blast resistance breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01486-5.

The GW5-WRKY53-SGW5 module regulates grain size variation in rice

Grain size is a crucial agronomic trait that affects stable yield, appearance, milling quality, and domestication in rice. However, the molecular and genetic relationships among QTL genes (QTGs) underlying natural variation for grain size remain elusive. Here, we identified a novel QTG SGW5 (suppressor of gw5) by map-based cloning using an F2 segregation population by fixing same genotype of the master QTG GW5. SGW5 positively regulates grain width by influencing cell division and cell size in spikelet hulls. Two nearly isogenic lines exhibited a significant differential expression of SGW5 and a 12.2% increase in grain yield. Introducing the higher expression allele into the genetic background containing the lower expression allele resulted in increased grain width, while its knockout resulted in shorter grain hulls and dwarf plants. Moreover, a cis-element variation in the SGW5 promoter influenced its differential binding affinity for the WRKY53 transcription factor, causing the differential SGW5 expression, which ultimately leads to grain size variation. GW5 physically and genetically interacts with WRKY53 to suppress the expression of SGW5. These findings elucidated a new pathway for grain size regulation by the GW5-WRKY53-SGW5 module and provided a novel case for generally uncovering QTG interactions underlying the genetic diversity of an important trait in crops.

Pan-transcriptomic analysis reveals alternative splicing control of cold tolerance in rice

Plants have evolved complex mechanisms to adapt to harsh environmental conditions. Rice (Oryza sativa) is a staple food crop that is sensitive to low temperatures. However, its cold stress responses remain poorly understood, thus limiting possibilities for crop engineering to achieve greater cold tolerance. In this study, we constructed a rice pan-transcriptome and characterized its transcriptional regulatory landscape in response to cold stress. We performed Iso-Seq and RNA-Seq of 11 rice cultivars subjected to a time-course cold treatment. Our analyses revealed that alternative splicing-regulated gene expression plays a significant role in the cold stress response. Moreover, we identified CATALASE C (OsCATC) and Os03g0701200 as candidate genes for engineering enhanced cold tolerance. Importantly, we uncovered central roles for the 2 serine-arginine-rich proteins OsRS33 and OsRS2Z38 in cold tolerance. Our analysis of cold tolerance and resequencing data from a diverse collection of 165 rice cultivars suggested that OsRS2Z38 may be a key selection gene in japonica domestication for cold adaptation, associated with the adaptive evolution of rice. This study systematically investigated the distribution, dynamic changes, and regulatory mechanisms of alternative splicing in rice under cold stress. Overall, our work generates a rich resource with broad implications for understanding the genetic basis of cold response mechanisms in plants.

GWAS analysis reveals the genetic basis of blast resistance associated with heading date in rice

Rice blast is a destructive fungal disease affecting rice plants at various growth stages, significantly threatening global yield stability. Development of resistant rice cultivars stands as a practical means of disease control. Generally, association mapping with a diversity panel powerfully identifies new alleles controlling trait of interest. On the other hand, utilization of a breeding panel has its advantage that can be directly applied in a breeding program. In this study, we conducted a genome-wide association study (GWAS) for blast resistance using 296 commercial rice cultivars with low population structure but large phenotypic diversity. We attempt to answer the genetic basis behind rice blast resistance among early maturing cultivars by subdividing the population based on its Heading date 1 (Hd1) functionality. Subpopulation-specific GWAS using the mixed linear model (MLM) based on blast nursery screening conducted in three years revealed a total of 26 significant signals, including three nucleotide-binding site leucine-rich repeat (NBS-LRR) genes (Os06g0286500, Os06g0286700, and Os06g0287500) located at Piz locus on chromosome 6, and one at the Pi-ta locus (Os12g0281300) on chromosome 12. Haplotype analysis revealed blast resistance associated with Piz locus was exclusively specific to Type 14 hd1 among japonica rice. Our findings provide valuable insights for breeding blast resistant rice and highlight the applicability of our elite cultivar panel to detect superior alleles associated with important agronomic traits.

New Insights into the Genetic Basis of Lysine Accumulation in Rice Revealed by Multi-Model GWAS

Lysine is an essential amino acid that cannot be synthesized in humans. Rice is a global staple food for humans but has a rather low lysine content. Identification of the quantitative trait nucleotides (QTNs) and genes underlying lysine content is crucial to increase lysine accumulation. In this study, five grain and three leaf lysine content datasets and 4,630,367 single nucleotide polymorphisms (SNPs) of 387 rice accessions were used to perform a genome-wide association study (GWAS) by ten statistical models. A total of 248 and 71 common QTNs associated with grain/leaf lysine content were identified. The accuracy of genomic selection/prediction RR-BLUP models was up to 0.85, and the significant correlation between the number of favorable alleles per accession and lysine content was up to 0.71, which validated the reliability and additive effects of these QTNs. Several key genes were uncovered for fine-tuning lysine accumulation. Additionally, 20 and 30 QTN-by-environment interactions (QEIs) were detected in grains/leaves. The QEI-sf0111954416 candidate gene LOC_Os01g21380 putatively accounted for gene-by-environment interaction was identified in grains. These findings suggested the application of multi-model GWAS facilitates a better understanding of lysine accumulation in rice. The identified QTNs and genes hold the potential for lysine-rich rice with a normal phenotype.

Oschib1 gene encoding a GH18 chitinase confers resistance against sheath blight disease of rice caused by Rhizoctonia solani AG1-IA

Sheath blight disease of rice caused by Rhizoctonia solani AG1-IA, is a major fungal disease responsible for huge loss to grain yield and quality. The major limitation of achieving persistent and reliable resistance against R. solani is the governance of disease resistance trait by many genes. Therefore, functional characterization of new genes involved in sheath blight resistance is necessary to understand the mechanism of resistance as well as evolving effective strategies to manage the disease through host-plant resistance. In this study, we performed RNA sequencing of six diverse rice genotypes (TN1, BPT5204, Vandana, N22, Tetep, and Pankaj) from sheath and leaf tissue of control and fungal infected samples. The approach for identification of candidate resistant genes led to identification of 352 differentially expressed genes commonly present in all the six genotypes. 23 genes were analyzed for RT-qPCR expression which helped identification of Oschib1 showing differences in expression level in a time-course manner between susceptible and resistant genotypes. The Oschib1 encoding classIII chitinase was cloned from resistant variety Tetep and over-expressed in susceptible variety Taipei 309. The over-expression lines showed resistance against R. solani, as analyzed by detached leaf and whole plant assays. Interestingly, the resistance response was correlated with the level of transgene expression suggesting that the enzyme functions in a dose dependent manner. We report here the classIIIb chitinase from chromosome10 of rice showing anti-R. solani activity to combat the dreaded sheath blight disease.

OsUGT88C3 Encodes a UDP-Glycosyltransferase Responsible for Biosynthesis of Malvidin 3-O-Galactoside in Rice

The diversity of anthocyanins is largely due to the action of glycosyltransferases, which add sugar moieties to anthocyanidins. Although a number of glycosyltransferases have been identified to glycosylate anthocyanidin in plants, the enzyme that catalyzes malvidin galactosylation remains unclear. In this study, we identified three rice varieties with different leaf color patterns, different anthocyanin accumulation patterns, and different expression patterns of anthocyanin biosynthesis genes (ABGs) to explore uridine diphosphate (UDP)-glycosyltransferases (UGTs) responsible for biosynthesis of galactosylated malvidin. Based on correlation analysis of transcriptome data, nine candidate UGT genes coexpressed with 12 ABGs were identified (r values range from 0.27 to 1.00). Further analysis showed that the expression levels of one candidate gene, OsUGT88C3, were highly correlated with the contents of malvidin 3-O-galactoside, and recombinant OsUGT88C3 catalyzed production of malvidin 3-O-galactoside using UDP-galactose and malvidin as substrates. OsUGT88C3 was closely related to UGTs with flavone and flavonol glycosylation activities in phylogeny. Its plant secondary product glycosyltransferase (PSPG) motif ended with glutamine. Haplotype analysis suggested that the malvidin galactosylation function of OsUGT88C3 was conserved among most of the rice germplasms. OsUGT88C3 was highly expressed in the leaf, pistil, and embryo, and its protein was located in the endoplasmic reticulum and nucleus. Our findings indicate that OsUGT88C3 is responsible for the biosynthesis of malvidin 3-O-galactoside in rice and provide insight into the biosynthesis of anthocyanin in plants.

Insertion of a miniature inverted-repeat transposable element into the promoter of OsTCP4 results in more tillers and a lower grain size in rice

A class I PCF type protein, TCP4, was identified as a transcription factor associated with both grain size and tillering through a DNA pull-down-MS assay combined with a genome-wide association study. This transcription factor was found to have a significant role in the variations among the 533 rice accessions, dividing them into two main subspecies. A Tourist-like miniature inverted-repeat transposable element (MITE) was discovered in the promoter of TCP4 in japonica/geng accessions (TCP4M+), which was found to suppress the expression of TCP4 at the transcriptional level. The MITE-deleted haplotype (TCP4M-) was mainly found in indica/xian accessions. ChIP-qPCR and EMSA demonstrated the binding of TCP4 to promoters of grain reservoir genes such as SSIIa and Amy3D in vivo and in vitro, respectively. The introduction of the genomic sequence of TCP4M+ into different TCP4M- cultivars was found to affect the expression of TCP4 in the transgenic rice, resulting in decreased expression of its downstream target gene SSIIa, increased tiller number, and decreased seed length. This study revealed that a Tourist-like MITE contributes to subspecies divergence by regulating the expression of TCP4 in response to environmental pressure, thus influencing source-sink balance by regulating starch biosynthesis in rice.

Transcription factor encoding gene OsC1 regulates leaf sheath color through anthocyanidin metabolism in Oryza rufipogon and Oryza sativa

Carbohydrates, proteins, lipids, minerals and vitamins are nutrient substances commonly seen in rice grains, but anthocyanidin, with benefit for plant growth and animal health, exists mainly in the common wild rice but hardly in the cultivated rice. To screen the rice germplasm with high intensity of anthocyanidins and identify the variations, we used metabolomics technique and detected significant different accumulation of anthocyanidins in common wild rice (Oryza rufipogon, with purple leaf sheath) and cultivated rice (Oryza sativa, with green leaf sheath). In this study, we identified and characterized a well-known MYB transcription factor, OsC1, through phenotypic (leaf sheath color) and metabolic (metabolite profiling) genome-wide association studies (pGWAS and mGWAS) in 160 common wild rice (O. rufipogon) and 151 cultivated (O. sativa) rice varieties. Transgenic experiments demonstrated that biosynthesis and accumulation of cyanidin-3-Galc, cyanidin 3-O-rutinoside and cyanidin O-syringic acid, as well as purple pigmentation in leaf sheath were regulated by OsC1. A total of 25 sequence variations of OsC1 constructed 16 functional haplotypes (higher accumulation of the three anthocyanidin types within purple leaf sheath) and 9 non-functional haplotypes (less accumulation of anthocyanidins within green leaf sheath). Three haplotypes of OsC1 were newly identified in our germplasm, which have potential values in functional genomics and molecular breeding of rice. Gene-to-metabolite analysis by mGWAS and pGWAS provides a useful and efficient tool for functional gene identification and omics-based crop genetic improvement.

A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset

BACKGROUND

Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable.

RESULTS

Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq).

CONCLUSIONS

This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

Uncovering Novel QTLs and Candidate Genes for Salt Tolerance at the Bud Burst Stage in Rice through Genome-Wide Association Study

Salt stress is one of the most important factors limiting rice growth and yield increase. Salt tolerance of rice at the bud burst (STB) stage determines whether germinated seeds can grow normally under salt stress, which is very important for direct seeding. However, reports on quantitative trait loci (QTLs) and candidate genes for STB in rice are very limited. In this study, a natural population of 130 indica and 81 japonica rice accessions was used to identify STB-related QTLs and candidate genes using a genome-wide association study (GWAS). Nine QTLs, including five for relative shoot length (RSL), two for relative root length (RRL), and two for relative root number (RRN), were identified. Five of these STB-related QTLs are located at the same site as the characterized salt tolerance genes, such as OsMDH1, OsSRFP1, and OsCDPK7. However, an important QTL related to RSL, qRSL1-2, has not been previously identified and was detected on chromosome 1. The candidate region for qRSL1-2 was identified by linkage disequilibrium analysis, 18 genes were found to have altered expression levels under salt stress through the RNA-seq database, and 10 of them were found to be highly expressed in the shoot. It was also found that, eight candidate genes (LOC_Os01g62980, LOC_Os01g63190, LOC_Os01g63230, LOC_Os01g63280, LOC_Os01g63400, LOC_Os01g63460, and LOC_Os01g63580) for qRSL1-2 carry different haplotypes between indica and japonica rice, which exactly corresponds to the significant difference in RSL values between indica and japonica rice in this study. Most of the accessions with elite haplotypes were indica rice, which had higher RSL values. These genes with indica-japonica specific haplotypes were identified as candidate genes. Rice accessions with elite haplotypes could be used as important resources for direct seeding. This study also provides new insights into the genetic mechanism of STB.

Identification of hub genes that variate the qCSS12-mediated cold tolerance between indica and japonica rice using WGCNA

KEY MESSAGE

Cold-tolerant QTL qCSS12-regulated 14 hub genes are involved in the chloroplastic biological processes and in the protein synthesis and degradation processes in japonica rice. Low temperature is a main constraint factor for rice growth and production. To better understand the regulatory mechanisms underlying the cold tolerance phenotype in rice, here, we selected a cold-sensitive nearly isogenic line (NIL) NIL(qcss12) as materials to identify hub genes that are mediated by the cold-tolerant locus qCSS12 through weighted gene co-expression network analysis (WGCNA). Fourteen cold-responsive genes were identified, of which, 6 are involved in regulating biological processes in chloroplasts, including the reported EF-Tu, Prk, and ChlD, and 8 are involved in the protein synthesis and degradation processes. Differential expression of these genes between NIL(qcss12) and its controls under cold stress may be responsible for qCSS12-mediated cold tolerance in japonica rice. Moreover, natural variations in 12 of these hub genes are highly correlated with the cold tolerance divergence in two rice subspecies. The results provide deep insights into a better understanding of the molecular basis of cold adaptation in rice and provide a theoretical basis for molecular breeding.

Navigating rice seedling cold resilience: QTL mapping in two inbred line populations and the search for genes

Due to global climate change resulting in extreme temperature fluctuations, it becomes increasingly necessary to explore the natural genetic variation in model crops such as rice to facilitate the breeding of climate-resilient cultivars. To uncover genomic regions in rice involved in managing cold stress tolerance responses and to identify associated cold tolerance genes, two inbred line populations developed from crosses between cold-tolerant and cold-sensitive parents were used for quantitative trait locus (QTL) mapping of two traits: degree of membrane damage after 1 week of cold exposure quantified as percent electrolyte leakage (EL) and percent low-temperature seedling survivability (LTSS) after 1 week of recovery growth. This revealed four EL QTL and 12 LTSS QTL, all overlapping with larger QTL regions previously uncovered by genome-wide association study (GWAS) mapping approaches. Within the QTL regions, 25 cold-tolerant candidate genes were identified based on genomic differences between the cold-tolerant and cold-sensitive parents. Of those genes, 20% coded for receptor-like kinases potentially involved in signal transduction of cold tolerance responses; 16% coded for transcription factors or factors potentially involved in regulating cold tolerance response effector genes; and 64% coded for protein chaperons or enzymes potentially serving as cold tolerance effector proteins. Most of the 25 genes were cold temperature regulated and had deleterious nucleotide variants in the cold-sensitive parent, which might contribute to its cold-sensitive phenotype.

Control of rice ratooning ability by a nucleoredoxin that inhibits histidine kinase dimerization to attenuate cytokinin signaling in axillary buds

Rice ratooning, the fast outgrowth of dormant buds on stubble, is an important cropping practice in rice production. However, the low ratooning ability (RA) of most rice varieties restricts the application of this cost-efficient system, and the genetic basis of RA remains unknown. In this study, we dissected the genetic architecture of RA by a genome-wide association study in a natural rice population. Rice ratooning ability 3 (RRA3), encoding a hitherto not characterized nucleoredoxin involved in reduction of disulfide bonds, was identified as the causal gene of a major locus controlling RA. Overexpression of RRA3 in rice significantly accelerated leaf senescence and reduced RA, whereas knockout of RRA3 significantly delayed leaf senescence and increased RA and ratoon yield. We demonstrated that RRA3 interacts with Oryza sativa histidine kinase 4 (OHK4), a cytokinin receptor, and inhibits the dimerization of OHK4 through disulfide bond reduction. This inhibition ultimately led to decreased cytokinin signaling and reduced RA. In addition, variations in the RRA3 promoter were identified to be associated with RA. Introgression of a superior haplotype with weak expression of RRA3 into the elite rice variety Guichao 2 significantly increased RA and ratoon yield by 23.8%. Collectively, this study not only uncovers an undocumented regulatory mechanism of cytokinin signaling through de-dimerization of a histidine kinase receptor-but also provides an eximious gene with promising value for ratoon rice breeding.

Yellow-Green Leaf 19 Encoding a Specific and Conservative Protein for Photosynthetic Organisms Affects Tetrapyrrole Biosynthesis, Photosynthesis, and Reactive Oxygen Species Metabolism in Rice

Chlorophyll is the main photosynthetic pigment and is crucial for plant photosynthesis. Leaf color mutants are widely used to identify genes involved in the synthesis or metabolism of chlorophyll. In this study, a spontaneous mutant, yellow-green leaf 19 (ygl19), was isolated from rice (Oryza sativa). This ygl19 mutant showed yellow-green leaves and decreased chlorophyll level and net photosynthetic rate. Brown necrotic spots appeared on the surface of ygl19 leaves at the tillering stage. And the agronomic traits of the ygl19 mutant, including the plant height, tiller number per plant, and total number of grains per plant, were significantly reduced. Map-based cloning revealed that the candidate YGL19 gene was LOC_Os03g21370. Complementation of the ygl19 mutant with the wild-type CDS of LOC_Os03g21370 led to the restoration of the mutant to the normal phenotype. Evolutionary analysis revealed that YGL19 protein and its homologues were unique for photoautotrophs, containing a conserved Ycf54 functional domain. A conserved amino acid substitution from proline to serine on the Ycf54 domain led to the ygl19 mutation. Sequence analysis of the YGL19 gene in 4726 rice accessions found that the YGL19 gene was conserved in natural rice variants with no resulting amino acid variation. The YGL19 gene was mainly expressed in green tissues, especially in leaf organs. And the YGL19 protein was localized in the chloroplast for function. Gene expression analysis via qRT-PCR showed that the expression levels of tetrapyrrole synthesis-related genes and photosynthesis-related genes were regulated in the ygl19 mutant. Reactive oxygen species (ROS) such as superoxide anions and hydrogen peroxide accumulated in spotted leaves of the ygl19 mutant at the tillering stage, accompanied by the regulation of ROS scavenging enzyme-encoding genes and ROS-responsive defense signaling genes. This study demonstrates that a novel yellow-green leaf gene YGL19 affects tetrapyrrole biosynthesis, photosynthesis, and ROS metabolism in rice.

Transcriptome-wide association analyses reveal the impact of regulatory variants on rice panicle architecture and causal gene regulatory networks

Panicle architecture is a key determinant of rice grain yield and is mainly determined at the 1-2 mm young panicle stage. Here, we investigated the transcriptome of the 1-2 mm young panicles from 275 rice varieties and identified thousands of genes whose expression levels were associated with panicle traits. Multimodel association studies suggested that many small-effect genetic loci determine spikelet per panicle (SPP) by regulating the expression of genes associated with panicle traits. We found that alleles at cis-expression quantitative trait loci of SPP-associated genes underwent positive selection, with a strong preference for alleles increasing SPP. We further developed a method that integrates the associations of cis- and trans-expression components of genes with traits to identify causal genes at even small-effect loci and construct regulatory networks. We identified 36 putative causal genes of SPP, including SDT (MIR156j) and OsMADS17, and inferred that OsMADS17 regulates SDT expression, which was experimentally validated. Our study reveals the impact of regulatory variants on rice panicle architecture and provides new insights into the gene regulatory networks of panicle traits.

The ribosomal protein P0A is required for embryo development in rice

BACKGROUND

The P-stalk is a conserved and vital structural element of ribosome. The eukaryotic P-stalk exists as a P0-(P1-P2)2 pentameric complex, in which P0 function as a base structure for incorporating the stalk onto 60S pre-ribosome. Prior studies have suggested that P0 genes are indispensable for survival in yeast and animals. However, the functions of P0 genes in plants remain elusive.

RESULTS

In the present study, we show that rice has three P0 genes predicted to encode highly conserved proteins OsP0A, OsP0B and OsP0C. All of these P0 proteins were localized both in cytoplasm and nucleus, and all interacted with OsP1. Intriguingly, the transcripts of OsP0A presented more than 90% of the total P0 transcripts. Moreover, knockout of OsP0A led to embryo lethality, while single or double knockout of OsP0B and OsP0C did not show any visible defects in rice. The genomic DNA of OsP0A could well complement the lethal phenotypes of osp0a mutant. Finally, sequence and syntenic analyses revealed that OsP0C evolved from OsP0A, and that duplication of genomic fragment harboring OsP0C further gave birth to OsP0B, and both of these duplication events might happen prior to the differentiation of indica and japonica subspecies in rice ancestor.

CONCLUSION

These data suggested that OsP0A functions as the predominant P0 gene, playing an essential role in embryo development in rice. Our findings highlighted the importance of P0 genes in plant development.

A Structure Variation in qPH8.2 Detrimentally Affects Plant Architecture and Yield in Rice

Plant height is an important agronomic trait associated with plant architecture and grain yield in rice (Oryza sativa L.). In this study, we report the identification of quantitative trait loci (QTL) for plant height using a chromosomal segment substitution line (CSSL) population with substituted segments from japonica variety Nipponbare (NIP) in the background of the indica variety 9311. Eight stable QTLs for plant height were identified in three environments. Among them, six loci were co-localized with known genes such as semidwarf-1 (sd1) and Grain Number per Panicle1 (GNP1) involved in gibberellin biosynthesis. A minor QTL qPH8.2 on chromosome 8 was verified and fine-mapped to a 74 kb region. Sequence comparison of the genomic region revealed the presence/absence of a 42 kb insertion between NIP and 9311. This insertion occurred predominantly in temperate japonica rice. Comparisons on the near-isogenic lines showed that the qPH8.2 allele from NIP exhibits pleiotropic effects on plant growth, including reduced plant height, leaf length, photosynthetic capacity, delayed heading date, decreased yield, and increased tiller angle. These results indicate that qPH8.2 from temperate japonica triggers adverse effects on plant growth and yield when introduced into the indica rice, highlighting the importance of the inter-subspecies crossing breeding programs.

Haplotype analysis and marker development of five salt-tolerant-related genes in rice (Oryza sativa L.)

Salinity stress is a great threat to the growth and productivity of crops, and development of salt-tolerant crops is of great necessity to ensure food security. Although a few genes with natural variations that confer salt tolerance at germination and seedling stage in rice have been cloned, effective intragenic markers for these genes are awaited to be developed, which hinder the use of these genes in genetic improvement of salt tolerance in rice. In this study, we first performed haplotype analysis of five rice salt-tolerant-related genes using 38 rice accessions with reference genome and 4,726 rice germplasm accessions with imputed genotypes and classified main haplotype groups and haplotypes. Subsequently, we identified unique variations for elite haplotypes reported in previous studies and developed 11 effective intragenic makers. Finally, we conducted genotyping of 533 of the 4,726 rice accessions from worldwide and 70 approved temperate geng/japonica cultivars in China using the developed markers. These results could provide effective donors and markers of salt-tolerant-related genes and thus could be of great use in genetic improvement of salt tolerance in rice.

Fantastic genes: where and how to find them? Exploiting rice genetic resources for the improvement of yield, tolerance, and resistance to a wide array of stresses in rice

Rice production is a critical component of global food security. To date, rice is grown in over 100 countries and is the primary source of food for more than 3 billion people. Despite its importance, rice production is facing numerous challenges that threaten its future viability. One of the primary problems is the advent of climate change. The changing climatic conditions greatly affect the growth and productivity of rice crop and the quality of rice yield. Similarly, biotic stresses brought about by pathogen and pest infestations are greatly affecting the productivity of rice. To address these issues, the utilization of rice genetic resources is necessary to map, identify, and understand the genetics of important agronomic traits. This review paper highlights the role of rice genetic resources for developing high-yielding and stress-tolerant rice varieties. The integration of genetic, genomic, and phenomic tools in rice breeding programs has led to the development of high-yielding and stress-tolerant rice varieties. The collaboration of multidisciplinary teams of experts, sustainable farming practices, and extension services for farmers is essential for accelerating the development of high-yielding and stress-tolerant rice varieties.

A phytoexclusion strategy for reducing contamination risk of rice based on low-Cd natural variations pyramid of root transporters

Cadmium pollution in rice is a threat to human health. Phytoexclusion is an effective strategy to reduce the Cd accumulation. Soil-to-root is the first step of Cd entering rice and plays a crucial role in Cd accumulation, so targeting root transporters could be an effective approach for phytoexclusion. This study utilized single-gene & multi-gene joint haplotype analysis to reveal the law of natural variations. The result showed that natural variations of rice root transporters assembled regularly following a certain pattern, rather than randomly. A total of 3 dominant nature variation combinations with 2 high-Cd combinations and 1 low-Cd combination were identified. In addition, indica-japonica differentiation was observed, with indica germplasms harboring high-Cd combinations while japonica germplasms harboring. In Chinese rice landraces, most of the collected indica landraces contained high-Cd combinations, indicating a high Cd contamination risk in indica landraces in terms of both phenotype and genotype. To address this issue, multiple superior low-Cd natural variations were pyramided to create two new low-Cd germplasms. In both pond and farmland trials, the ameliorated rice grain Cd did not exceed safety standards. This research provided a framework for future phytoexclusion, thus to reduce Cd-contamination risk in soil-rice system.

Genetic Diversity and Breeding Signatures for Regional Indica Rice Improvement in Guangdong of Southern China

As the pioneer of the Green Revolution in China, Guangdong province witnessed the improvement and spread of semi-dwarf Xian/Indica rice cultivars and possessed diverse rice germplasm of landrace and cultivars. A total of 517 accessions containing a core germplasm of 479 newly sequenced landraces and modern cultivars were used to reveal breeding signatures and key variations for regional genetic improvement of indica rice from Guangdong. Four subpopulations were identified in the collection, which including Ind IV as a novel subpopulation that not covered by previously released accessions. Modern cultivars of subpopulation Ind II were inferred to have less deleterious variations, especially in yield related genes. About 15 Mb genomic segments were identified as potential breeding signatures by cross-population likelihood method (XP-CLR) of modern cultivars and landraces. The selected regions spanning multiple yield related QTLs (quantitative trait locus) which identified by GWAS (genome-wide association studies) of the same population, and specific variations that fixed in modern cultivars of Ind II were characterized. This study highlights genetic differences between traditional landraces and modern cultivars, which revealed the potential molecular basis of regional genetic improvement for Guangdong indica rice from southern China.

geneHapR: an R package for gene haplotypic statistics and visualization

BACKGROUND

Together with application of next-generation sequencing technologies and increased accumulation of genomic variation data in different organism species, an opportunity for effectively identification of superior alleles of functional genes to facilitate marker-assisted selection is emerging, and the clarification of haplotypes of functional genes is becoming an essential target in recent study works.

RESULTS

In this paper, we describe an R package 'geneHapR' developed for haplotypes identification, statistics and visualization analysis of candidate genes. This package could integrate genotype data, genomic annotating information and phenotypic variation data to clarify genotype variations, evolutionary-ship, and morphological effects among haplotypes through variants visualization, network construction and phenotypic comparison. 'geneHapR' also provides functions for Linkage Disequilibrium block analysis and visualizing of haplotypes geo-distribution.

CONCLUSIONS

The R package 'geneHapR' provided an easy-to-use tool for haplotype identification, statistic and visualization for candidate gene and will provide useful clues for gene functional dissection and molecular-assistant pyramiding of beneficial alleles of functional locus in future breeding programs.

BOS1 is a basic helix-loop-helix transcription factor involved in regulating panicle development in rice

Panicle development is crucial to increase the grain yield of rice (Oryza sativa). The molecular mechanisms of the control of panicle development in rice remain unclear. In this study, we identified a mutant with abnormal panicles, termed branch one seed 1-1 (bos1-1). The bos1-1 mutant showed pleiotropic defects in panicle development, such as the abortion of lateral spikelets and the decreased number of primary panicle branches and secondary panicle branches. A combined map-based cloning and MutMap approach was used to clone BOS1 gene. The bos1-1 mutation was located in chromosome 1. A T-to-A mutation in BOS1 was identified, which changed the codon from TAC to AAC, resulting in the amino acid change from tyrosine to asparagine. BOS1 gene encoded a grass-specific basic helix-loop-helix transcription factor, which is a novel allele of the previously cloned LAX PANICLE 1 (LAX1) gene. Spatial and temporal expression profile analyses showed that BOS1 was expressed in young panicles and was induced by phytohormones. BOS1 protein was mainly localized in the nucleus. The expression of panicle development-related genes, such as OsPIN2, OsPIN3, APO1, and FZP, was changed by bos1-1 mutation, suggesting that the genes may be the direct or indirect targets of BOS1 to regulate panicle development. The analysis of BOS1 genomic variation, haplotype, and haplotype network showed that BOS1 gene had several genomic variations and haplotypes. These results laid the foundation for us to further dissect the functions of BOS1.

Double Mutation of Days to Heading 2 and CONSTANS 3 Improves Agronomic Performance of Japonica Rice under Short Daylight Conditions in Southern China

Some progress has been made in understanding the pathways related to rice heading, but their applications to breeding japonica rice varieties adapted to grow in low-latitude areas ("indica to japonica") are limited. We edited eight adaptation-related genes via a lab-established CRISPR/Cas9 system in a japonica variety, Shennong265 (SN265). All T₀ plants and their progeny bearing random mutation permutations were planted in southern China and screened for changes in heading date. We found that the double mutant of Days to heading 2 (DTH2) and CONSTANS 3 (OsCO3) (dth2-osco3), two CONSTANS-like (COL) genes, showed significantly delayed heading under both short-day (SD) and long-day (LD) conditions in Guangzhou and manifested great yield increase under SD conditions. We further demonstrated that the heading-related Hd3a-OsMADS14 pathway was down-regulated in the dth2-osco3 mutant lines. The editing of the COL genes DTH2 and OsCO3 greatly improves the agronomic performance of japonica rice in Southern China.

A minimal genome design to maximally guarantee fertile inter-subspecific hybrid rice

Hybrid rice has made considerable contributions to achieve the ambitious goal of food security for the world's population. Hybrid rice from indica/xian and japonica/geng subspecies shows much higher heterosis and is thereby an important innovation in promoting rice production in the next decade. However, such inter-subspecific hybrid rice has long suffered from serious hybrid sterility, which is a major challenge that needs to be addressed. In this study, we performed a genome design strategy to produce fertile inter-subspecific hybrid by creation of wide compatibility varieties that are able to overcome hybrid sterility. Based on combined genetic analyses in two indica-japonica crosses, we determined that four hybrid sterility loci, S5, f5, pf12 and Sc, are the major QTLs controlling inter-subspecific hybrid sterility and thus the minimal targets that can be manipulated for breeding sub-specific hybrid rice. We then cloned the pf12 locus, one of the most effective loci for hybrid male sterility, by map-based cloning, and showed that artificial disruption of pf12A gene at this locus could successfully rescue hybrid fertility. We further dissected the genetic basis of wide compatibility using three pairwise crosses from a wide-compatibility variety Dular and representative indica and japonica varieties. On this basis, we constructed and assembled different combinations of naturally compatible alleles of four loci, S5, Sc, pf12, and f5, and found that the improved lines could fully recover pollen and embryo sac fertility in test-crossed F₁s, thereby completely fulfilling the demands of inter-subspecific hybrid spikelet fertility in agricultural production. This breeding scheme would facilitate redesign of future inter-subspecific hybrid rice with a higher yield potential.

Natural variations of OsAUX5, a target gene of OsWRKY78, control the neutral essential amino acid content in rice grains

Grain essential amino acid (EAA) levels contribute to rice nutritional quality. However, the molecular mechanisms underlying EAA accumulation and natural variation in rice grains remain unclear. Here we report the identification of a previously unrecognized auxin influx carrier subfamily gene, OsAUX5, which encodes an amino acid transporter that functions in uptake of multiple amino acids. We identified an elite haplotype of Pro::OsAUX5^Hap2 that enhances grain EAA accumulation without an apparent negative effect on agronomic traits. Natural variations of OsAUX5 occur in the cis elements of its promoter, which are differentially activated because of the different binding affinity between OsWRKY78 and the W-box, contributing to grain EAA variation among rice varieties. The two distinct haplotypes were shown to have originated from different Oryza rufipogon progenitors, which contributed to the divergence between japonica and indica. Introduction of the indica-type Pro::OsAUX5^Hap2 genotype into japonica could significantly increase EAA levels, indicating that indica-type Pro::OsAUX5^Hap2 can be utilized to increase grain EAAs of japonica varieties. Collectively, our study uncovers an WRKY78-OsAUX5-based regulatory mechanism controlling grain EAA accumulation and provides a potential target for breeding EAA-rich rice.