Assessment of Five Chilling Tolerance Traits and GWAS Mapping in Rice Using the USDA Mini-Core Collection

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.

Unlocking ABA's role in rice cold tolerance: insights from Zhonghua 11 and Kasalath

KEY MESSAGE

Unraveling key ABA pathways, including OsWRKY71-OsABA8ox1 and OsbZIP73-OsNCED5, provides valuable insights for improving cold tolerance in rice breeding for cold-prone regions. Cold stress limits rice (Oryza sativa L.) production in cooler climates. This study uncovers how abscisic acid (ABA) signaling enhances cold tolerance in the rice variety Zhonghua 11 (ZH11) compared to the cold-sensitive Kasalath. Under cold stress, ZH11 rapidly accumulates ABA through efficient regulation of key genes. The transcription factor OsWRKY71ZH11 represses the ABA catabolism gene OsABA8ox1 during early stress, enabling quick ABA accumulation. Additionally, OsbZIP73 regulates the ABA synthesis gene OsNCED5 to maintain ABA balance during prolonged stress. Transgenic ZH11 plants overexpressing OsWRKY71ZH11 exhibited enhanced cold tolerance, while overexpression of OsWRKY71Ka did not confer benefits. Haplotype analysis linked allelic variations in OsWRKY71 and OsNCED5 to differences in cold tolerance. Our findings highlight critical ABA signaling pathways that enhance cold tolerance in rice. Targeting these pathways offers promising strategies for breeding cold-resistant rice varieties, improving resilience in cold-prone regions.

Analysis of genome-wide association studies of low-temperature germination in Xian and Geng rice

Rice is the leading global staple crop. Low temperatures pose negative impacts on rice's optimal growth and development. Rice cultivars acclimating to low temperatures exhibited improved seedling emergence under direct-seeded sowing conditions, yet little is known about the genes that regulate germination at low temperatures (LTG). In this research investigation, we've performed whole genome sequencing for the 273 rice plant materials. Using the best linear unbiased prediction (BLUP) values for each rice material, we identified 7 LTG-related traits and performed the efficient genetic analysis and genome-wide association study (GWAS). As a result of this, 95 quantitative trait loci (QTLs) and 1001 candidate genes associated with LTG in rice were identified. Haplotype analysis and functional annotation of the candidate genes resulted in the identification of three promising candidate genes (LOC_Os08g30520 for regulating LTG4 and LTG5, LOC_Os10g02625 for regulating LTG6, LTg7 and LTG8, and LOC_Os12g31460 for regulating LTG7, LTg8 and LTG9) involving in the regulation of LTG in rice. This research provides a solid foundation for addressing the LTG issue in rice and will be valuable in future direct-seeded rice breeding programs.

OsEIN2-OsEIL1/2 pathway negatively regulates chilling tolerance by attenuating OsICE1 function in rice

Low temperature severely affects rice development and yield. Ethylene signal is essential for plant development and stress response. Here, we reported that the OsEIN2-OsEIL1/2 pathway reduced OsICE1-dependent chilling tolerance in rice. The overexpressing plants of OsEIN2, OsEIL1 and OsEIL2 exhibited severe stress symptoms with excessive reactive oxygen species (ROS) accumulation under chilling, while the mutants (osein2 and oseil1) and OsEIL2-RNA interference plants (OsEIL2-Ri) showed the enhanced chilling tolerance. We validated that OsEIL1 and OsEIL2 could form a heterxodimer and synergistically repressed OsICE1 expression by binding to its promoter. The expression of OsICE1 target genes, ROS scavenging- and photosynthesis-related genes were downregulated by OsEIN2 and OsEIL1/2, which were activated by OsICE1, suggesting that OsEIN2-OsEIL1/2 pathway might mediate ROS accumulation and photosynthetic capacity under chilling by attenuating OsICE1 function. Moreover, the association analysis of the seedling chilling tolerance with the haplotype showed that the lower expression of OsEIL1 and OsEIL2 caused by natural variation might confer chilling tolerance on rice seedlings. Finally, we generated OsEIL2-edited rice with an enhanced chilling tolerance. Taken together, our findings reveal a possible mechanism integrating OsEIN2-OsEIL1/2 pathway with OsICE1-dependent cascade in regulating chilling tolerance, providing a practical strategy for breeding chilling-tolerant rice.

Mapping of Candidate Genes for Nitrogen Uptake and Utilization in Japonica Rice at Seedling Stage

Nitrogen (N) is one of the essential nutrients for the growth and development of crops. The adequate application of N not only increases the yield of crops but also improves the quality of agricultural products, but the excessive application of N can cause many adverse effects on ecology and the environment. In this study, genome-wide association analysis (GWAS) was performed under low- and high-N conditions based on 788,396 SNPs and phenotypic traits relevant to N uptake and utilization (N content and N accumulation). A total of 75 QTLs were obtained using GWAS, which contained 811 genes. Of 811 genes, 281 genes showed different haplotypes, and 40 genes had significant phenotypic differences among different haplotypes. Of these 40 genes, 5 differentially expressed genes (Os01g0159250, Os02g0618200, Os02g0618400, Os02g0630300, and Os06g0619000) were finally identified as the more valuable candidate genes based on the transcriptome data sequenced from Longjing31 (low-N-tolerant variety) and Songjing 10 (low-N-sensitive variety) under low- and high-N treatments. These new findings enrich the genetic resources for N uptake and utilization in rice, as well as lay a theoretical foundation for improving the efficiency of N uptake and utilization in rice.

Pseudomonas mosselii improves cold tolerance of Asian rice (Oryza sativa L.) in a genotype-dependent manner by increasing proline in japonica and reduced glutathione in indica varieties

Cold stress is an important factor limiting rice production and distribution. Identifying factors that contribute to cold tolerance in rice is of primary importance. While some plant specific genetic factors involved in cold tolerance have been identified, the role of the rice microbiome remains unexplored. In this study, we evaluated the influence of plant growth promoting bacteria (PGPB) with the ability of phosphate solubilization on rice cold tolerance and survival. To reach this goal, inoculated and uninoculated 2-week-old seedlings were cold stressed and evaluated for survival and other phenotypes such as electrolyte leakage (EL) and necessary elements for cold tolerance. The results of this study showed that of the five bacteria, Pseudomonas mosselii, improved both indica and japonica varietal plants' survival and decreased EL, indicating increased membrane integrity. We observed different possible cold tolerance mechanisms in japonica and indica plants such as increases in proline and reduced glutathione levels, respectively. This bacterium also improved the shoot growth of cold exposed indica plants during the recovery period. This study confirmed the host genotype dependent activity of P. mosselii and indicated that there is an interaction between specific plant genes and bacterial genes that causes different plant responses to cold stress.

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.

AP2/EREBP Pathway Plays an Important Role in Chaling Wild Rice Tolerance to Cold Stress

Cold stress is the main factor limiting rice production and distribution. Chaling wild rice can survive in cold winters. AP2/EREBP is a known transcription factor family associated with abiotic stress. We identified the members of the AP2/EREBP transcription factor family in rice, maize, and Arabidopsis, and conducted collinearity analysis and gene family analysis. We used Affymetrix array technology to analyze the expression of AP2/EREBP family genes in Chaling wild rice and cultivated rice cultivar Pei'ai64S, which is sensitive to cold. According to the GeneChip results, the expression levels of AP2/EREBP genes in Chaling wild rice were different from those in Pei'ai64S; and the increase rate of 36 AP2/EREBP genes in Chaling wild rice was higher than that in Pei'ai64S. Meanwhile, the MYC elements in cultivated rice and Chaling wild rice for the Os01g49830, Os03g08470, and Os03g64260 genes had different promoter sequences, resulting in the high expression of these genes in Chaling wild rice under low-temperature conditions. Furthermore, we analyzed the upstream and downstream genes of the AP2/EREBP transcription factor family and studied the conservation of these genes. We found that the upstream transcription factors were more conserved, indicating that these upstream transcription factors may be more important in regulating cold stress. Meanwhile, we found the expression of AP2/EREBP pathway genes was significantly increased in recombinant inbred lines from Nipponbare crossing with Chaling wild rice, These results suggest that the AP2/EREBP signaling pathway plays an important role in Chaling wild rice tolerance to cold stress.

The basic helix-loop-helix transcription factor gene, OsbHLH38, plays a key role in controlling rice salt tolerance

The plant hormone abscisic acid (ABA) is crucial for plant seed germination and abiotic stress tolerance. However, the association between ABA sensitivity and plant abiotic stress tolerance remains largely unknown. In this study, 436 rice accessions were assessed for their sensitivity to ABA during seed germination. The considerable diversity in ABA sensitivity among rice germplasm accessions was primarily reflected by the differentiation between the Xian (indica) and Geng (japonica) subspecies and between the upland-Geng and lowland-Geng ecotypes. The upland-Geng accessions were most sensitive to ABA. Genome-wide association analyses identified four major quantitative trait loci containing 21 candidate genes associated with ABA sensitivity of which a basic helix-loop-helix transcription factor gene, OsbHLH38, was the most important for ABA sensitivity. Comprehensive functional analyses using knockout and overexpression transgenic lines revealed that OsbHLH38 expression was responsive to multiple abiotic stresses. Overexpression of OsbHLH38 increased seedling salt tolerance, while knockout of OsbHLH38 increased sensitivity to salt stress. A salt-responsive transcription factor, OsDREB2A, interacted with OsbHLH38 and was directly regulated by OsbHLH38. Moreover, OsbHLH38 affected rice abiotic stress tolerance by mediating the expression of a large set of transporter genes of phytohormones, transcription factor genes, and many downstream genes with diverse functions, including photosynthesis, redox homeostasis, and abiotic stress responsiveness. These results demonstrated that OsbHLH38 is a key regulator in plant abiotic stress tolerance.

Genetic analysis of protein content and oil content in soybean by genome-wide association study

Soybean seed protein content (PC) and oil content (OC) have important economic value. Detecting the loci/gene related to PC and OC is important for the marker-assisted selection (MAS) breeding of soybean. To detect the stable and new loci for PC and OC, a total of 320 soybean accessions collected from the major soybean-growing countries were used to conduct a genome-wide association study (GWAS) by resequencing. The PC ranged from 37.8% to 46.5% with an average of 41.1% and the OC ranged from 16.7% to 22.6% with an average of 21.0%. In total, 23 and 29 loci were identified, explaining 3.4%-15.4% and 5.1%-16.3% of the phenotypic variations for PC and OC, respectively. Of these, eight and five loci for PC and OC, respectively, overlapped previously reported loci and the other 15 and 24 loci were newly identified. In addition, nine candidate genes were identified, which are known to be involved in protein and oil biosynthesis/metabolism, including lipid transport and metabolism, signal transduction, and plant development pathway. These results uncover the genetic basis of soybean protein and oil biosynthesis and could be used to accelerate the progress in enhancing soybean PC and OC.

GWAS reveals novel loci and identifies a pentatricopeptide repeat-containing protein (CsPPR) that improves low temperature germination in cucumber

Low temperatures (LTs) negatively affect the percentage and rate of cucumber (Cucumis sativus L.) seed germination, which has deleterious effects on yield. Here, a genome-wide association study (GWAS) was used to identify the genetic loci underlying low temperature germination (LTG) in 151 cucumber accessions that represented seven diverse ecotypes. Over two years, phenotypic data for LTG i.e., relative germination rate (RGR), relative germination energy (RGE), relative germination index (RGI) and relative radical length (RRL), were collected in two environments, and 17 of the 151 accessions were found to be highly cold tolerant using cluster analysis. A total of 1,522,847 significantly associated single-nucleotide polymorphism (SNP) were identified, and seven loci associated with LTG, on four chromosomes, were detected: gLTG1.1, gLTG1.2, gLTG1.3, gLTG4.1, gLTG5.1, gLTG5.2, and gLTG6.1 after resequencing of the accessions. Of the seven loci, three, i.e., gLTG1.2, gLTG4.1, and gLTG5.2, showed strong signals that were consistent over two years using the four germination indices, and are thus strong and stable for LTG. Eight candidate genes associated with abiotic stress were identified, and three of them were potentially causal to LTG: CsaV3_1G044080 (a pentatricopeptide repeat-containing protein) for gLTG1.2, CsaV3_4G013480 (a RING-type E3 ubiquitin transferase) for gLTG4.1, and CsaV3_5G029350 (a serine/threonine-protein kinase) for gLTG5.2. The function for CsPPR (CsaV3_1G044080) in regulating LTG was confirmed, as Arabidopsis lines ectopically expressing CsPPR showed higher germination and survival rates at 4°C compared to the wild-type, which preliminarily illustrates that CsPPR positively regulates cucumber cold tolerance at the germination stage. This study will provide insights into cucumber LT-tolerance mechanisms and further promote cucumber breeding development.

Genome-wide association study identifies GhSAL1 affects cold tolerance at the seedling emergence stage in upland cotton (Gossypium hirsutum L.)

Genomic analysis of upland cotton revealed that cold tolerance was associated with ecological distribution. GhSAL1 on chromosome D09 negatively regulated cold tolerance of upland cotton. Cotton can undergo low-temperature stress at the seedling emergence stage, which adversely affects growth and yield; however, the regulatory mechanism underlying cold tolerance remains nebulous. Here, we analyze the phenotypic and physiological parameters in 200 accessions from 5 ecological distributions under constant chilling (CC) and diurnal variation of chilling (DVC) stresses at the seedling emergence stage. All accessions were clustered into four groups, of which Group IV, with most germplasms from the northwest inland region (NIR), had better phenotypes than Groups I-III under the two kinds of chilling stresses. A total of 575 significantly associated single-nucleotide polymorphism (SNP) were identified, and 35 stable genetic quantitative trait loci (QTL) were obtained, of which 5 were associated with traits under CC and DVC stress, respectively, while the remaining 25 were co-associated. The accumulation of dry weight (DW) of seedling was associated with the flavonoid biosynthesis process regulated by Gh_A10G0500. The emergence rate (ER), DW, and total length of seedling (TL) under CC stress were associated with the SNPs variation of Gh_D09G0189 (GhSAL1). GhSAL1^HapB was the elite haplotype, which increased ER, DW, and TL by 19.04%, 11.26%, and 7.69%, respectively, compared with that of GhSAL1^HapA. The results of virus-induced gene silencing (VIGS) experiment and determination of metabolic substrate content preliminarily illustrated that GhSAL1 negatively regulated cotton cold tolerance through IP₃-Ca²⁺ signaling pathway. The elite haplotypes and candidate genes identified in this study could be used to improve cold tolerance at the seedling emergence stage in future upland cotton breeding.

Oligomeric Proanthocyanidins Confer Cold Tolerance in Rice through Maintaining Energy Homeostasis

Oligomeric proanthocyanidins (OPCs) are abundant polyphenols found in foods and botanicals that benefit human health, but our understanding of the functions of OPCs in rice plants is limited, particularly under cold stress. Two rice genotypes, named Zhongzao39 (ZZ39) and its recombinant inbred line RIL82, were subjected to cold stress. More damage was caused to RIL82 by cold stress than to ZZ39 plants. Transcriptome analysis suggested that OPCs were involved in regulating cold tolerance in the two genotypes. A greater increase in OPCs content was detected in ZZ39 than in RIL82 plants under cold stress compared to their respective controls. Exogenous OPCs alleviated cold damage of rice plants by increasing antioxidant capacity. ATPase activity was higher and poly (ADP-ribose) polymerase (PARP) activity was lower under cold stress in ZZ39 than in RIL82 plants. Importantly, improvements in cold tolerance were observed in plants treated with the OPCs and 3-aminobenzamide (PARP inhibitor, 3ab) combination compared to the seedling plants treated with H₂O, OPCs, or 3ab alone. Therefore, OPCs increased ATPase activity and inhibited PARP activity to provide sufficient energy for rice seedling plants to develop antioxidant capacity against cold stress.

Multi-omics revolution to promote plant breeding efficiency

Crop production is the primary goal of agricultural activities, which is always taken into consideration. However, global agricultural systems are coming under increasing pressure from the rising food demand of the rapidly growing world population and changing climate. To address these issues, improving high-yield and climate-resilient related-traits in crop breeding is an effective strategy. In recent years, advances in omics techniques, including genomics, transcriptomics, proteomics, and metabolomics, paved the way for accelerating plant/crop breeding to cope with the changing climate and enhance food production. Optimized omics and phenotypic plasticity platform integration, exploited by evolving machine learning algorithms will aid in the development of biological interpretations for complex crop traits. The precise and progressive assembly of desire alleles using precise genome editing approaches and enhanced breeding strategies would enable future crops to excel in combating the changing climates. Furthermore, plant breeding and genetic engineering ensures an exclusive approach to developing nutrient sufficient and climate-resilient crops, the productivity of which can sustainably and adequately meet the world's food, nutrition, and energy needs. This review provides an overview of how the integration of omics approaches could be exploited to select crop varieties with desired traits.

Global identification of quantitative trait loci and candidate genes for cold stress and chilling acclimation in rice through GWAS and RNA-seq

Associated analysis of GWAS with RNA-seq had detected candidate genes responsible for cold stress and chilling acclimation in rice. Haplotypes of two candidate genes and geographic distribution were analyzed. To explore new candidate genes and genetic resources for cold tolerance improvement in rice, genome-wide association study (GWAS) mapping experiments with 351 rice core germplasms was performed for three traits (survival rate, shoot length and chlorophyll content) under three temperature conditions (normal temperature, cold stress and chilling acclimation), yielding a total of 134 QTLs, of which 54, 59 and 21 QTLs were responsible for normal temperature, cold stress and chilling acclimation conditions, respectively. Integrated analysis of significant SNPs in 134 QTLs further identified 116 QTLs for three temperature treatments, 53, 43 and 18 QTLs responsible for normal temperature, cold stress and chilling acclimation, respectively, and 2 QTLs were responsible for both cold stress and chilling acclimation. Matching differentially expressed genes from RNA-seq to 43 and 18 QTLs for cold stress and chilling acclimation, we identified 69 and 44 trait-associated candidate genes, respectively, to be classified into six and five groups, particularly involved in metabolisms, reactive oxygen species scavenging and hormone signaling. Interestingly, two candidate genes LOC_Os01g04814, encoding a vacuolar protein sorting-associating protein 4B, and LOC_Os01g48440, encoding glycosyltransferase family 43 protein, showed the highest expression levels under chilling acclimation. Haplotype analysis revealed that both genes had a distinctive differentiation with subpopulation. Haplotypes of both genes with more japonica accessions have higher latitude distribution and higher chilling tolerance than the chilling sensitive indica accessions. These findings reveal the new insight into the molecular mechanism and candidate genes for cold stress and chilling acclimation in rice.

Assessment of Rice Sheath Blight Resistance Including Associations with Plant Architecture, as Revealed by Genome-Wide Association Studies

BACKGROUND

Sheath blight (ShB) disease caused by Rhizoctonia solani Kühn, is one of the most economically damaging rice (Oryza sativa L.) diseases worldwide. There are no known major resistance genes, leaving only partial resistance from small-effect QTL to deploy for cultivar improvement. Many ShB-QTL are associated with plant architectural traits detrimental to yield, including tall plants, late maturity, or open canopy from few or procumbent tillers, which confound detection of physiological resistance.

RESULTS

To identify QTL for ShB resistance, 417 accessions from the Rice Diversity Panel 1 (RDP1), developed for association mapping studies, were evaluated for ShB resistance, plant height and days to heading in inoculated field plots in Arkansas, USA (AR) and Nanning, China (NC). Inoculated greenhouse-grown plants were used to evaluate ShB using a seedling-stage method to eliminate effects from height or maturity, and tiller (TN) and panicle number (PN) per plant. Potted plants were used to evaluate the RDP1 for TN and PN. Genome-wide association (GWA) mapping with over 3.4 million SNPs identified 21 targeted SNP markers associated with ShB which tagged 18 ShB-QTL not associated with undesirable plant architecture traits. Ten SNPs were associated with ShB among accessions of the Indica subspecies, ten among Japonica subspecies accessions, and one among all RDP1 accessions. Across the 18 ShB QTL, only qShB4-1 was not previously reported in biparental mapping studies and qShB9 was not reported in the GWA ShB studies. All 14 PN QTL overlapped with TN QTL, with 15 total TN QTL identified. Allele effects at the five TN QTL co-located with ShB QTL indicated that increased TN does not inevitably increase disease development; in fact, for four ShB QTL that overlapped TN QTL, the alleles increasing resistance were associated with increased TN and PN, suggesting a desirable coupling of alleles at linked genes.

CONCLUSIONS

Nineteen accessions identified as containing the most SNP alleles associated with ShB resistance for each subpopulation were resistant in both AR and NC field trials. Rice breeders can utilize these accessions and SNPs to develop cultivars with enhanced ShB resistance along with increased TN and PN for improved yield potential.

Phenotypic and Molecular Characterization of Rice Genotypes’ Tolerance to Cold Stress at the Seedling Stage

Rice plants are affected by low-temperature stress during germination, vegetative growth, and reproductive stages. Thirty-nine rice genotypes including 36 near-isogenic lines (NILs) of BRRI dhan29 were evaluated to investigate the level of cold tolerance under artificially induced low temperature at the seedling stage. Three cold-related traits, leaf discolouration (LD), survivability, and recovery rate, were measured to determine the level of cold tolerance. Highly significant variation among the genotypes was observed for LD, survivability, and recovery rate. Three NILs, IR90688-74-1-1-1-1-1, IR90688-81-1-1-1-1-1, and IR90688-103-1-1-1-1-1, showed tolerance in all three traits, while IR90688-118-1-1-1-1-1 showed cold tolerance with LD and recovery rate. IR90688-92-1-1-1-1-1, IR90688-125-1-1-1-1-1, IR90688-104-1-1-1-1-1, IR90688-124-1-1-1-1-P2, IR90688-15-1-1-1-1-1, and IR90688-27-1-1-1-1-1 showed significantly higher yield coupled with short growth duration and good grain quality. Genetic analysis with SSRs markers revealed that the high-yielding NILs were genetically 67% similar to BRRI dhan28 and possessed cold tolerance at the seedling stage. These cold-tolerant NILs could be used as potential resources to broaden the genetic base of the breeding germplasm to develop high-yielding cold-tolerant rice varieties.

Natural variation in the fast phase of chlorophyll a fluorescence induction curve (OJIP) in a global rice minicore panel

Photosynthesis can be probed through Chlorophyll a fluorescence induction (FI), which provides detailed insight into the electron transfer process in Photosystem II, and beyond. Here, we have systematically studied the natural variation of the fast phase of the FI, i.e. the OJIP phase, in rice. The OJIP phase of the Chl a fluorescence induction curve is referred to as "fast transient" lasting for less than a second; it is obtained after a dark-adapted sample is exposed to saturating light. In the OJIP curve, "O" stands for "origin" (minimal fluorescence), "P" for "peak" (maximum fluorescence), and J and I for inflection points between the O and P levels. Further, F_o is the fluorescence intensity at the "O" level, whereas F_m is the intensity at the P level, and F_v (= F_m - F_o) is the variable fluorescence. We surveyed a set of quantitative parameters derived from the FI curves of 199 rice accessions, grown under both field condition (FC) and growth room condition (GC). Our results show a significant variation between Japonica (JAP) and Indica (IND) subgroups, under both the growth conditions, in almost all the parameters derived from the OJIP curves. The ratio of the variable to the maximum (F_v/F_m) and of the variable to the minimum (F_v/F_o) fluorescence, the performance index (PI_abs), as well as the amplitude of the I-P phase (A_I-P) show higher values in JAP compared to that in the IND subpopulation. In contrast, the amplitude of the O-J phase (A_O-J) and the normalized area above the OJIP curve (S_m) show an opposite trend. The performed genetic analysis shows that plants grown under GC appear much more affected by environmental factors than those grown in the field. We further conducted a genome-wide association study (GWAS) using 11 parameters derived from plants grown in the field. In total, 596 non-unique significant loci based on these parameters were identified by GWAS. Several photosynthesis-related proteins were identified to be associated with different OJIP parameters. We found that traits with high correlation are usually associated with similar genomic regions. Specifically, the thermal phase of FI, which includes the amplitudes of the J-I and I-P subphases (A_J-I and A_I-P) of the OJIP curve, is, in turn, associated with certain common genomic regions. Our study is the first one dealing with the natural variations in rice, with the aim to characterize potential candidate genes controlling the magnitude and half-time of each of the phases in the OJIP FI curve.

Novel QTLs for salinity tolerance revealed by genome-wide association studies of biomass, chlorophyll and tissue ion content in 176 rice landraces from Bangladesh

Farmland is on the decline and worldwide food security is at risk. Rice is the staple of choice for over half the Earth's people. To sustain current demands and ascertain a food secure future, substandard farmland affected by abiotic stresses must be utilized. For rapid crop improvement, a broader understanding of polygenic traits like stress tolerance and crop yield is indispensable. To this end, the hidden diversity of resilient and neglected wild varieties must be traced back to their genetic roots. In this study, we separately assayed 11 phenotypes in a panel of 176 diverse accessions predominantly comprised of local landraces from Bangladesh. We compiled high resolution sequence data for these accessions. We collectively studied the ties between the observed phenotypic differences and the examined additive genetic effects underlying these variations. We applied a fixed effect model to associate phenotypes with genotypes on a genomic scale. Discovered QTLs were mapped to known genes. Our explorations yielded 13 QTLs related to various traits in multiple trait classes. 10 identified QTLs were equivalent to findings from previous studies. Integrative analysis assumes potential novel functionality for a number of candidate genes. These findings will usher novel avenues for the bioengineering of high yielding crops of the future fortified with genetic defenses against abiotic stressors.

Low Temperature Antioxidant Activity QTL Associate with Genomic Regions Involved in Physiological Cold Stress Tolerance Responses in Rice (Oryza sativa L.)

Boosting cold stress tolerance in crop plants can minimize stress-mediated yield losses. Asian rice (Oryza sativa L.), one of the most consumed cereal crops, originated from subtropical regions and is generally sensitive to low temperature environments. Within the two subspecies of rice, JAPONICA, and INDICA, the cold tolerance potential of its accessions is highly variable and depends on their genetic background. Yet, cold stress tolerance response mechanisms are complex and not well understood. This study utilized 370 accessions from the Rice Diversity Panel 1 (RDP1) to investigate and correlate four cold stress tolerance response phenotypes: membrane damage, seedling survivability, and catalase and anthocyanin antioxidative activity. Most JAPONICA accessions, and admixed accessions within JAPONICA, had lower membrane damage, higher antioxidative activity, and overall, higher seedling survivability compared to INDICA accessions. Genome-wide association study (GWAS) mapping was done using the four traits to find novel quantitative trait loci (QTL), and to validate and fine-map previously identified QTL. A total of 20 QTL associated to two or more traits were uncovered by our study. Gene Ontology (GO) term enrichment analyses satisfying four layers of filtering retrieved three potential pathways: signal transduction, maintenance of plasma membrane and cell wall integrity, and nucleic acids metabolism as general mechanisms of cold stress tolerance responses involving antioxidant activity.

Molecular and genetic dissection of the USDA rice mini-core collection using high-density SNP markers

Molecular tools and knowledge of crop germplasm are vital for their effective utilization. In this study, we developed 40,866 high-quality and well distributed SNPs for a rice mini-core collection (RMC) developed by the United States Department of Agriculture (USDA). The high-quality SNPs clustered the USDA-RMC into five subpopulations (Ind, indica; Aus, aus; Afr, African rice; TeJ, temperate japonica; TrJ, tropical japonica) and one admixture (Adm). This classification was further confirmed by phylogenetic and principal component analyses. The rice ARO (aromatic) subpopulation of previous studies was re-assigned with Adm and the WD (wild-type) subpopulation was re-defined to the Afr subpopulation because most of its accessions are African cultivated rice. The Aus and Ind subpopulations had a substantially wider genetic variation than the TrJ and TeJ subpopulations. The genetic diversities were much larger between the Ind or Aus subpopulation and the TrJ or TeJ subpopulation than between the Afr subpopulation and the Ind, Aus, TrJ or TeJ subpopulation. Comparative agronomic trait analysis between the subpopulations also supported the genetic structure and variation of the RMC, and suggested the existence of extensive variation in the genes controlling agronomic traits among them. Furthermore, analysis of ancestral membership of the RMC accessions revealed that reproductive barrier or wide incompatibility existed between the Indica and Japonica groups, while gene flow occurred between them. These results provide high-quality SNPs and knowledge of genetic structure and diversity of the USDA-RMC necessary for enhanced rice research and breeding.

Identification of Candidate Genes Conferring Cold Tolerance to Rice (Oryza sativa L.) at the Bud-Bursting Stage Using Bulk Segregant Analysis Sequencing and Linkage Mapping

Low-temperature tolerance during the bud-bursting stage is an important characteristic of direct-seeded rice. The identification of cold-tolerance quantitative trait loci (QTL) in species that can stably tolerate cold environments is crucial for the molecular breeding of rice with such traits. In our study, high-throughput QTL-sequencing analyses were performed in a 460-individual F_{2 : 3} mapping population to identify the major QTL genomic regions governing cold tolerance at the bud-bursting (CTBB) stage in rice. A novel major QTL, qCTBB9, which controls seed survival rate (SR) under low-temperature conditions of 5°C/9 days, was mapped on the 5.40-Mb interval on chromosome 9. Twenty-six non-synonymous single-nucleotide polymorphism (nSNP) markers were designed for the qCTBB9 region based on re-sequencing data and local QTL mapping conducted using traditional linkage analysis. We mapped qCTBB9 to a 483.87-kb region containing 58 annotated genes, among which six predicted genes contained nine nSNP loci. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis revealed that only Os09g0444200 was strongly induced by cold stress. Haplotype analysis further confirmed that the SNP 1,654,225 bp in the Os09g0444200 coding region plays a key role in regulating the cold tolerance of rice. These results suggest that Os09g0444200 is a potential candidate for qCTBB9. Our results are of great significance to explore the genetic mechanism of rice CTBB and to improve the cold tolerance of rice varieties by marker-assisted selection.

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.

The Candidate Genes Underlying a Stably Expressed QTL for Low Temperature Germinability in Rice (Oryza sativa L.)

BACKGROUND

Direct seeding is an efficient cultivation technique in rice. However, poor low temperature germinability (LTG) of modern rice cultivars limits its application. Identifying the genes associated with LTG and performing molecular breeding is the fundamental way to address this issue. However, few LTG QTLs have been fine mapped and cloned so far.

RESULTS

In the present study, the LTG evaluation of 375 rice accessions selected from the Rice Diversity Panel 2 showed that there were large LTG variations within the population, and the LTG of Indica group was significantly higher than that of Japonica and Aus groups (p < 0.01). In total, eleven QTLs for LTG were identified through genome-wide association study (GWAS). Among them, qLTG_sRDP2-3/qLTG_JAP-3, qLTG_AUS-3 and qLTG_sRDP2-12 are first reported in the present study. The QTL on chromosome 10, qLTG_sRDP2-10a had the largest contribution to LTG variations in 375 rice accessions, and was further validated using single segment substitution line (SSSL). The presence of qLTG_sRDP2-10a could result in 59.8% increase in LTG under 15 °C low temperature. The expression analysis of the genes within qLTG_sRDP2-10a region indicated that LOC_Os10g22520 and LOC_Os10g22484 exhibited differential expression between the high and low LTG lines. Further sequence comparisons revealed that there were insertion and deletion sequence differences in the promoter and intron region of LOC_Os10g22520, and an about 6 kb variation at the 3' end of LOC_Os10g22484 between the high and low LTG lines, suggesting that the sequence variations of the two genes could be the cause for their differential expression in high and low LTG lines.

CONCLUSION

Among the 11 QTLs identified in this study, qLTG_sRDP2-10a could also be detected in other three studies using different germplasm under different cold environments. Its large effect and stable expression make qLTG_sRDP2-10a particularly valuable in rice breeding. The two genes, LOC_Os10g22484 and LOC_Os10g22520, were considered as the candidate genes underlying qLTG_sRDP2-10a. Our results suggest that integrating GWAS and SSSL can facilitate identification of QTL for complex traits in rice. The identification of qLTG_sRDP2-10a and its candidate genes provide a promising source for gene cloning of LTG and molecular breeding for LTG in rice.

Multiple Cold Tolerance Trait Phenotyping Reveals Shared Quantitative Trait Loci in Oryza sativa

BACKGROUND

Developing chilling tolerant accessions of domesticated Asian rice is a potential source of significant crop improvement. The uniquely chilling sensitive nature of the tropically originating Oryza sativa make it the most important cereal crop that can gain significantly from improved tolerance to low temperatures. However, mechanisms underlying this complex trait are not fully understood. Oryza sativa has two subspecies with different levels of chilling tolerance, JAPONICA and INDICA, providing an ideal tool to investigate mechanistic differences in the chilling stress tolerance responses within this important crop species.

RESULTS

The Rice Diversity Panel 1 (RDP1) was used to investigate a core set of Oryza sativa accessions. The tools available for this panel allowed for a comprehensive analysis of two chilling tolerance traits at multiple temperatures across a 354-cultivar subset of the RDP1. Chilling tolerance trait values were distributed as mostly subpopulation specific clusters of Tolerant, Intermediate, and Sensitive accessions. Genome-wide association study (GWAS) mapping approaches using all 354 accessions yielded a total of 245 quantitative trait loci (QTL), containing 178 unique QTL covering 25% of the rice genome, while 40 QTL were identified by multiple traits. QTL mappings using subsets of rice accession clusters yielded another 255 QTL, for a total of 500 QTL. The genes within these multiple trait QTL were analyzed for Gene Ontology (GO) term and potential pathway enrichments. Terms related to "carbohydrate biosynthesis", "carbohydrate transmembrane transport", "small molecule protein modification", and "plasma membrane" were enriched from this list. Filtering was done to identify more likely candidate pathways involved in conferring chilling tolerance, resulting in enrichment of terms related to "Golgi apparatus", "stress response", "transmembrane transport", and "signal transduction".

CONCLUSIONS

Taken together, these GO term clusters revealed a likely involvement of Golgi-mediated subcellular and extracellular vesicle and intracellular carbohydrate transport as a general cold stress tolerance response mechanism to achieve cell and metabolic homeostasis under chilling stress.

Glycosyltransferase OsUGT90A1 helps protect the plasma membrane during chilling stress in rice

Due to its subtropical origins, rice (Oryza sativa) is sensitive to low-temperature stress. In this study, we identify LOC_Os04g24110, annotated to encode the UDP-glycosyltransferase enzyme UGT90A1, as a gene associated with the low-temperature seedling survivability (LTSS) quantitative trait locus qLTSS4-1. Differences between haplotypes in the control region of OsUGT90A1 correlate with chilling tolerance phenotypes, and reflect differential expression between tolerant and sensitive accessions rather than differences in protein sequences. Expression of OsUGT90A1 is initially enhanced by low temperature, and its overexpression helps to maintain membrane integrity during cold stress and promotes leaf growth during stress recovery, which are correlated with reduced levels of reactive oxygen species due to increased activities of antioxidant enzymes. In addition, overexpression of OsUGT90A1 in Arabidopsis improves freezing survival and tolerance to salt stress, again correlated with enhanced activities of antioxidant enzymes. Overexpression of OsUGT90A1 in rice decreases root lengths in 3-week-old seedlings while gene-knockout increases the length, indicating that its differential expression may affect phytohormone activities. We conclude that higher OsUGT90A1 expression in chilling-tolerant accessions helps to maintain cell membrane integrity as an abiotic stress-tolerance mechanism that prepares plants for the resumption of growth and development during subsequent stress recovery.

Genome-Wide Association Mapping to Identify Genetic Loci for Cold Tolerance and Cold Recovery During Germination in Rice

Low temperature significantly affects rice growth and yield. Temperatures lower than 15°C are generally detrimental for germination and uniform seedling stand. To investigate the genetic architecture underlying cold tolerance during germination in rice, we conducted a genome-wide association study using a novel diversity panel of 257 rice accessions from around the world and the 7K SNP marker array. Phenotyping was conducted in controlled growth chambers under dark conditions at 13°C. The rice accessions were measured for low-temperature germinability, germination index, coleoptile length under cold stress, plumule length at 4-day recovery, and plumule length recovery rate. A total of 51 QTLs were identified at p < 0.001 and 17 QTLs were identified using an FDR < 0.05 across the different chilling indices with the whole panel of accessions. At the threshold of p < 0.001, a total of 20 QTLs were identified in the subset of japonica accessions, while 9 QTLs were identified in the subset of indica accessions. Considering the recurring SNPs and linked SNPs across different chilling indices, we identified 31 distinct QTL regions in the whole panel, 13 QTL regions in the japonica subset, and 7 distinct QTL regions in the indica subset. Among these QTL regions, three regions were common between the whole panel and japonica, three regions were common between the whole panel and indica, and one region was common between indica and japonica. A subset of QTL regions was potentially colocalized with previously identified genes and QTLs, including 10 from the japonica subset, 4 from the indica subset, and 6 from the whole panel. On the other hand, a total of 21 potentially novel QTL regions from the whole panel, 10 from the japonica subset, and 1 from the indica subset were identified. The results of our study provide useful information on the genetic architecture underlying cold tolerance during germination in rice, which in turn can be used for further molecular study and crop improvement for low-temperature stressed environments.

Identification of main-effect quantitative trait loci (QTLs) for low-temperature stress tolerance germination- and early seedling vigor-related traits in rice (Oryza sativa L.)

An attempt was made in the current study to identify the main-effect and co-localized quantitative trait loci (QTLs) for germination and early seedling growth traits under low-temperature stress (LTS) conditions in rice. The plant material used in this study was an early backcross population of 230 introgression lines (ILs) in BCIF7 generation derived from the Weed Tolerant Rice-1 (WTR-1) (as the recipient) and Haoannong (HNG) (as the donor). Genetic analyses of LTS tolerance revealed a total of 27 main-effect quantitative trait loci (M-QTLs) mapped on 12 chromosomes. These QTLs explained more than 10% of phenotypic variance (PV), and average PV of 12.71% while employing 704 high-quality SNP markers. Of these 27 QTLs distributed on 12 chromosomes, 11 were associated with low-temperature germination (LTG), nine with low-temperature germination stress index (LTGS), five with root length stress index (RLSI), and two with biomass stress index (BMSI) QTLs, shoot length stress index (SLSI) and root length stress index (RLSI), seven with seed vigor index (SVI), and single QTL with root length (RL). Among them, five significant major QTLs (qLTG(I) 1 , qLTGS(I) 1-2 , qLTG(I) 5 , qLTGS(I) 5 , and qLTG(I) 7 ) mapped on chromosomes 1, 5, and 7 were associated with LTG and LTGS traits and the PV explained ranged from 16 to 23.3%. The genomic regions of these QTLs were co-localized with two to six QTLs. Most of the QTLs were growth stage-specific and found to harbor QTLs governing multiple traits. Eight chromosomes had more than four QTLs and were clustered together and designated as promising LTS tolerance QTLs (qLTTs), as qLTT 1 , qLTT 2 , qLTT 3 , qLTT 5 , qLTT 6 , qLTT 8 , qLTT 9 , and qLTT 11 . A total of 16 putative candidate genes were identified in the major M-QTLs and co-localized QTL regions distributed on different chromosomes. Overall, these significant genomic regions of M-QTLs are responsible for multiple traits and this suggested that these could serve as the best predictors of LTS tolerance at germination and early seedling growth stages. Furthermore, it is necessary to fine-map these regions and to find functional markers for marker-assisted selection in rice breeding programs for cold tolerance.

Metabolic profile analysis and identification of key metabolites during rice seed germination under low-temperature stress

The metabolic profile of rice (Oryza sativa) during germination under low temperature (LT) has not been reported. In this study, the rice varieties 02428 (japonica) and YZX (indica) were subjected to experiments consisting of treatments including LT, normal temperature (NT) and a transition from LT to NT, and tissues were sampled at different time points during germination. A total of 730 metabolites were detected by a liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based widely targeted metabolomics method. On the basis of the screening criteria of increased contents under LT and decreased contents under NT, we identified 35 different metabolites that responded to LT stress among the 730 metabolites. Furthermore, the content differences of the 35 metabolites were compared when the samples were transferred from LT to NT. According to a fold change <0.5 or a variable importance in projection (VIP) score>1 at the transition point, 7 out of the 35 metabolites responded significantly to LT stress and were defined as key metabolites. A partial least squares (PLS) regression model of seven key metabolites with seedling length (SL), seedling area (SSA), and seedling volume (SV) was constructed, and the fitting effect was good. These seven key metabolites participate in the biosynthesis of amino acids and phenylpropanoids and in the metabolism of glutathione and inositol phosphate. This study laid a foundation for an improved understanding of the LT-germination mechanism of rice seeds.

Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection

Salt stress is a major constraint to rice acreage and production worldwide. The purpose of this study was to evaluate the natural genetic variation available in the United States Department of Agriculture (USDA) rice mini-core collection (URMC) for early vigor traits under salt stress and identify quantitative trait loci (QTLs) for seedling-stage salt tolerance via a genome-wide association study (GWAS). Using a hydroponic system, the seedlings of 162 accessions were subjected to electrical conductivity (EC) 6.0 dS m⁻¹ salt stress at the three-to-four leaf stage. After completion of the study, 59.4% of the accessions were identified as sensitive, 23.9% were identified as moderately tolerant, and 16.7% were identified as highly tolerant. Pokkali was the most tolerant variety, while Nerica-6 was the most sensitive. Adapting standard International Rice Research Institute (IRRI) protocols, eight variables associated with salt tolerance were determined. The GWAS of the URMC, using over three million single-nucleotide polymorphisms (SNPs), identified nine genomic regions associated with salt tolerance that were mapped to five different chromosomes. Of these, none were in the known Saltol QTL region, suggesting different probable genes and mechanisms responsible for salt tolerance in the URMC. The study uncovered genetic loci that explained a large portion of the variation in salt tolerance at the seedling stage. Fourteen highly salt-tolerant accessions, six novel loci, and 16 candidate genes in their vicinity were identified that may be useful in breeding for salt stress tolerance. Identified QTLs can be targeted for fine mapping, candidate gene verification, and marker-assisted breeding in future studies.

Association Analysis of Three Diverse Rice (Oryza sativa L.) Germplasm Collections for Loci Regulating Grain Quality Traits

Rice ( L.) end-use cooking quality is vital for producers and billions of consumers worldwide. Grain quality is a complex trait with interacting genetic and environmental factors. Deciphering the complex genetic architecture associated with grain quality provides essential information for improved breeding strategies to enhance desirable traits that are stable across variable climatic and environmental conditions. In this study, genome-wide association (GWA) analysis of three rice diversity panels, the USDA rice core subset (1364 accessions), the minicore (MC) (173 accessions after removing non-), and the high density rice array-MC (HDMC) (383 accessions), with simple sequence repeats, single nucleotide polymorphic markers, or both, revealed large- and small-effect loci associated with known genes and previously uncharacterized genomic regions. Clustering of the significant regions in the GWA results suggests that multiple grain quality traits are inherited together. The 11 novel candidate loci for grain quality traits and the seven candidates for grain chalk identified are involved in the starch biosynthesis pathway. This study highlights the intricate pleiotropic relationships that exist in complex genotype-phenotypic associations and gives a greater insight into effective breeding strategies for grain quality improvement.

Translational genomics and multi-omics integrated approaches as a useful strategy for crop breeding

Recent next generation sequencing-driven mass production of genomic data and multi-omics-integrated approaches have significantly contributed to broadening and deepening our knowledge on the molecular system of living organisms. Accordingly, translational genomics (TG) approach can play a pivotal role in creating an informational bridge between model systems and relatively less studied plants. This review focuses mainly on addressing recent advancement in omics-related technologies, a diverse array of bioinformatic resources and potential applications of TG for the crop breeding. To accomplish above objectives, information on omics data production, various DBs and high throughput technologies was collected, integrated, and used to analyze current status and future perspectives towards omics-assisted crop breeding. Various omics data and resources have been organized and integrated into the databases and/or bioinformatic infrastructures, and thereby serve as the ome's information center for cross-genome translation of biological data. Although the size of accumulated omics data and availability of reference genomes are different among plant families, translational approaches have been actively progressing to access particular biological characteristics. When multi-layered omics data are integrated in a synthetic manner, it will allow providing a stereoscopic view of dynamic molecular behavior and interacting networks of genes occurring in plants. Consequently, TG approach will lead us to broader and deeper insights into target traits for the plant breeding. Furthermore, such systems approach will renovate conventional breeding programs and accelerate precision crop breeding in the future.

Dissecting the Genetic Architecture of Melon Chilling Tolerance at the Seedling Stage by Association Mapping and Identification of the Elite Alleles

Low temperature is an important abiotic stress that negatively affects morphological growth and fruit development in melon (Cucumis melo L.). Chilling stress at the seedling stage causes seedling injury and poor stand establishment, prolonging vegetative growth and delaying fruit harvest. In this study, association mapping was performed for chilling tolerance at the seedling stage on an expanded melon core collection containing 212 diverse accessions by 272 SSRs and 27 CAPSs. Chilling tolerance of the melon seedlings was evaluated by calculating the chilling injury index (CII) in 2016 and 2017. Genetic diversity analysis of the whole accession panel presented two main groups, which corresponded to the two subspecies of C. melo, melo, and agrestis. Both the subspecies were sensitive to chilling but with agrestis being more tolerant. Genome-wide association study (GWAS) was conducted, respectively, on the whole panel and the two subspecies, totally detecting 51 loci that contributed to 74 marker-trait associations. Of these associations, 35 were detected in the whole panel, 21 in melo, and 18 in agrestis. About half of the associations identified in the two subspecies were also observed in the whole panel, and seven associations were shared by both the subspecies. CMCT505_Chr.1 was repeatedly detected in different populations with high phenotypic contribution and could be a key locus controlling chilling tolerance in C. melo. Nine loci were selected for evaluation of the phenotypic effects related to their alleles, which identified 11 elite alleles contributing to seedling chilling tolerance. Four such alleles existed in both the subspecies and six in either of the two subspecies. Analysis of 20 parental combinations for their allelic status and phenotypic values showed that the elite alleles collectively contributed to enhancement of the chilling tolerance. Tagging the loci responsible for chilling tolerance may simultaneously favor dissecting the complex adaptability traits and elevate the efficiency to improve chilling tolerance using marker-assisted selection in melon.

Early selection of bZIP73 facilitated adaptation of japonica rice to cold climates

Cold stress is a major factor limiting production and geographic distribution of rice (Oryza sativa). Although the growth range of japonica subspecies has expanded northward compared to modern wild rice (O. rufipogon), the molecular basis of the adaptation remains unclear. Here we report bZIP73, a bZIP transcription factor-coding gene with only one functional polymorphism (+511 G>A) between the two subspecies japonica and indica, may have facilitated japonica adaptation to cold climates. We show the japonica version of bZIP73 (bZIP73^Jap) interacts with bZIP71 and modulates ABA levels and ROS homeostasis. Evolutionary and population genetic analyses suggest bZIP73 has undergone balancing selection; the bZIP73^Jap allele has firstly selected from standing variations in wild rice and likely facilitated cold climate adaptation during initial japonica domestication, while the indica allele bZIP73^Ind was subsequently selected for reasons that remain unclear. Our findings reveal early selection of bZIP73^Jap may have facilitated climate adaptation of primitive rice germplasms.

Japonica rice can grow further north than wild or indica rice and is more tolerant of cold climates. Here, the authors show that bZIP73 likely underwent selection in the early phase of rice domestication to facilitate cold tolerance in japonica by modulating ABA and ROS homeostasis.

Genome-Wide Association Study of Cadmium Accumulation at the Seedling Stage in Rapeseed (Brassica napus L.)

Cadmium is a potentially toxic heavy metal to human health. Rapeseed (Brassica napus L.), a vegetable and oilseed crop, might also be a Cd hyperaccumulator, but there is little information on this trait in rapeseed. We evaluated Cd accumulation in different oilseed accessions and employed a genome-wide association study to identify quantitative trait loci (QTLs) related to Cd accumulation. A total of 419 B. napus accessions and inbred lines were genotyped with a 60K Illumina Infinium SNP array of Brassica. Wide genotypic variations in Cd concentration and translocation were found. Twenty-five QTLs integrated with 98 single-nucleotide polymorphisms (SNPs) located at 15 chromosomes were associated with Cd accumulation traits. These QTLs explained 3.49-7.57% of the phenotypic variation observed. Thirty-two candidate genes were identified in these genomic regions, and they were 0.33-497.97 kb away from the SNPs. We found orthologs of Arabidopsis thaliana located near the significant SNPs on the B. napus genome, including NRAMP6 (natural resistance-associated macrophage protein 6), IRT1 (iron-regulated transporter 1), CAD1 (cadmium-sensitive 1), and PCS2 (phytochelatin synthase 2). Of them, four candidate genes were verified by qRT-PCR, the expression levels of which were significantly higher after exposure to Cd than in the controls. Our results might facilitate the study of the genetic basis of Cd accumulation and the cloning of candidate Cd accumulation genes, which could be used to help reduce Cd levels in edible plant parts and/or create more efficient hyperaccumulators.

Genome-wide association study of cold tolerance of Chinese indica rice varieties at the bud burst stage

A region containing three genes on chromosome 1 of indica rice was associated with cold tolerance at the bud burst stage; these results may be useful for breeding cold-tolerant lines. Low temperature at the bud burst stage is one of the major abiotic stresses limiting rice growth, especially in regions where rice seeds are sown directly. In this study, we investigated cold tolerance of rice at the bud burst stage and conducted a genome-wide association study (GWAS) based on the 5K rice array of 249 indica rice varieties widely distributed in China. We improved the method to assess cold tolerance at the bud burst stage in indica rice, and used severity of damage (SD) and seed survival rate (SR) as the cold-tolerant indices. Population structure analysis demonstrated that the Chinese indica panel was divided into three subgroups. In total, 47 significant single-nucleotide polymorphism (SNP) loci associated with SD and SR, were detected by association mapping based on mixed linear model. Because some loci overlapped between SD and SR, the loci contained 13 genome intervals and most of them have been reported previously. A major QTL for cold tolerance on chromosome 1 at the position of 31.6 Mb, explaining 13.2% of phenotypic variation, was selected for further analysis. Through LD decay, GO enrichment, RNA-seq data, and gene expression pattern analyses, we identified three genes (LOC_Os01g55510, LOC_Os01g55350 and LOC_Os01g55560) that were differentially expressed between cold-tolerant and cold-sensitive varieties, suggesting they may be candidate genes for cold tolerance. Together, our results provide a new method to assess cold tolerance in indica rice, and establish the foundation for isolating genes related to cold tolerance that could be used in rice breeding.

Differential transcriptome profiling of chilling stress response between shoots and rhizomes of Oryza longistaminata using RNA sequencing

Rice (Oryza sativa) is very sensitive to chilling stress at seedling and reproductive stages, whereas wild rice, O. longistaminata, tolerates non-freezing cold temperatures and has overwintering ability. Elucidating the molecular mechanisms of chilling tolerance (CT) in O. longistaminata should thus provide a basis for rice CT improvement through molecular breeding. In this study, high-throughput RNA sequencing was performed to profile global transcriptome alterations and crucial genes involved in response to long-term low temperature in O. longistaminata shoots and rhizomes subjected to 7 days of chilling stress. A total of 605 and 403 genes were respectively identified as up- and down-regulated in O. longistaminata under 7 days of chilling stress, with 354 and 371 differentially expressed genes (DEGs) found exclusively in shoots and rhizomes, respectively. GO enrichment and KEGG pathway analyses revealed that multiple transcriptional regulatory pathways were enriched in commonly induced genes in both tissues; in contrast, only the photosynthesis pathway was prevalent in genes uniquely induced in shoots, whereas several key metabolic pathways and the programmed cell death process were enriched in genes induced only in rhizomes. Further analysis of these tissue-specific DEGs showed that the CBF/DREB1 regulon and other transcription factors (TFs), including AP2/EREBPs, MYBs, and WRKYs, were synergistically involved in transcriptional regulation of chilling stress response in shoots. Different sets of TFs, such as OsERF922, OsNAC9, OsWRKY25, and WRKY74, and eight genes encoding antioxidant enzymes were exclusively activated in rhizomes under long-term low-temperature treatment. Furthermore, several cis-regulatory elements, including the ICE1-binding site, the GATA element for phytochrome regulation, and the W-box for WRKY binding, were highly abundant in both tissues, confirming the involvement of multiple regulatory genes and complex networks in the transcriptional regulation of CT in O. longistaminata. Finally, most chilling-induced genes with alternative splicing exclusive to shoots were associated with photosynthesis and regulation of gene expression, while those enriched in rhizomes were primarily related to stress signal transduction; this indicates that tissue-specific transcriptional and post-transcriptional regulation mechanisms synergistically contribute to O. longistaminata long-term CT. Our findings provide an overview of the complex regulatory networks of CT in O. longistaminata.