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

Meta-QTL and ortho analysis unravels the genetic architecture and key candidate genes for cold tolerance at seedling stage in rice

Rice, a critical cereal crop, grapples with productivity challenges due to its inherent sensitivity to low temperatures, primarily during the seedling and booting stages. Recognizing the polygenic complexity of cold stress signaling in rice, a meta-analysis was undertaken, focusing on 20 physiological traits integral to cold tolerance. This initiative allowed the consolidation of genetic data from 242 QTLs into 58 meta-QTLs, thereby significantly constricting the genetic and physical intervals, with 84% of meta-QTLs (MQTLs) being reduced to less than 2 Mb. The list of 10,505 genes within these MQTLs, was further refined utilizing expression datasets to pinpoint 46 pivotal genes exhibiting noteworthy differential regulation during cold stress. The study underscored the presence of several TFs such as WRKY, NAC, CBF/DREB, MYB, and bHLH, known for their roles in cold stress response. Further, ortho-analysis involving maize, barley, and Arabidopsis identified OsWRKY71, among others, as a prospective candidate for enhancing cold tolerance in diverse crop plants. In conclusion, our study delineates the intricate genetic architecture underpinning cold tolerance in rice and propounds significant candidate genes, offering crucial insights for further research and breeding strategies focused on fortifying crops against cold stress, thereby bolstering global food resilience.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01412-1.

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.

Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice

The resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.

Biotechnological Advances to Improve Abiotic Stress Tolerance in Crops

The major challenges that agriculture is facing in the twenty-first century are increasing droughts, water scarcity, flooding, poorer soils, and extreme temperatures due to climate change. However, most crops are not tolerant to extreme climatic environments. The aim in the near future, in a world with hunger and an increasing population, is to breed and/or engineer crops to tolerate abiotic stress with a higher yield. Some crop varieties display a certain degree of tolerance, which has been exploited by plant breeders to develop varieties that thrive under stress conditions. Moreover, a long list of genes involved in abiotic stress tolerance have been identified and characterized by molecular techniques and overexpressed individually in plant transformation experiments. Nevertheless, stress tolerance phenotypes are polygenetic traits, which current genomic tools are dissecting to exploit their use by accelerating genetic introgression using molecular markers or site-directed mutagenesis such as CRISPR-Cas9. In this review, we describe plant mechanisms to sense and tolerate adverse climate conditions and examine and discuss classic and new molecular tools to select and improve abiotic stress tolerance in major crops.

Construction of a high-density genetic map and localization of grazing-tolerant QTLs in Medicago falcata L

BACKGROUND

Using genomic DNA from 79 F1 plants resulted from a crossing between parents with strong and weak grazing tolerance in Medicago falcata L., we generated an EcoRI restriction site-associated DNA (RAD) sequencing library. After sequencing and assembly, a high-density genetic map with high-quality SNP markers was constructed, with a total length of 1312.238 cM and an average density of 0.844 SNP/cM.

METHODS

The phenotypic traits of 79 F1 families were observed and the QTLS of 6 traits were analyzed by interval mapping.

RESULTS

Sixty three QTLs were identified for seven traits with LOD values from 3 to 6 and the contribution rates from 15% to 30%. Among the 63 QTLs, 17 were for natural shoot height, 12 for rhizome Length, 10 for Shoot canopy diameter, 9 for Basal plant diameter, 6 for stem number, 5 for absolute shoot height, and 4 for rhizome width. These QTLs were concentrated on LG2, LG4, LG5, LG7, and LG8. LG6 had only 6 QTLs. According to the results of QTL mapping, comparison of reference genomes, and functional annotation, 10 candidate genes that may be related to grazing tolerance were screened. qRT-PCR analysis showed that two candidate genes (LOC11412291 and LOC11440209) may be the key genes related to grazing tolerance of M. falcata.

CONCLUSION

The identified trait-associated QTLs and candidate genes in this study will provide a solid foundation for future molecular breeding for enhanced grazing-tolerance in M. falcata.

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.

Plant Growth-Promoting Activity of Bacteria Isolated from Asian Rice (Oryza sativa L.) Depends on Rice Genotype

Asian rice is one of the most important crops because it is a staple food for almost half of the world's population. To have production of rice keep pace with a growing world population, it is anticipated that the use of fertilizers will also need to increase, which may cause environmental damage through runoff impacts. An alternative strategy to increase crop yield is the use of plant growth-promoting bacteria. Thousands of microbial species can exist in association with plant roots and shoots, and some are critical to the plant's survival. We isolated 140 bacteria from two distantly related rice accessions and investigated whether their impact on the growth of four different rice accessions. The bacterial isolates were screened for their ability to solubilize phosphate, a known plant growth-promoting characteristic, and 25 isolates were selected for further analysis. These 25 phosphate-solubilizing isolates were also able to produce other potentially growth-promoting factors. Five of the most promising bacterial isolates were chosen for whole-genome sequencing. Four of these bacteria, isolates related to Pseudomonas mosselii, a Microvirga sp., Paenibacillus rigui, and Paenibacillus graminis, improved root and shoot growth in a rice genotype-dependent manner. This indicates that while bacteria have several known plant growth-promoting functions, their effects on growth parameters are rice genotype dependent and suggest a close relationship between plants and their microbial partners. IMPORTANCE In this study, endophytic bacterial isolates from roots and shoots of two distantly related rice accessions were characterized phenotypically and genotypically. From the isolated bacterial species, five of the most promising plant growth-promoting bacteria were selected to test their abilities to enhance growth of the four rice accessions. Interestingly, plant growth enhancement was both bacterial isolate specific and plant genotype specific. However, the positive interactions between plant and bacteria could not easily be predicted because rice growth-promoting bacteria isolated from their host plants did not necessarily stimulate growth of their own host.

Population Evolution, Genetic Diversity and Structure of the Medicinal Legume, Glycyrrhiza uralensis and the Effects of Geographical Distribution on Leaves Nutrient Elements and Photosynthesis

Glycyrrhiza uralensis is a valuable medicinal legume, which occurs widely in arid and semi-arid regions. G. uralensis demand has risen steeply due to its high medical and commercial value. Interpret genome-wide information can stimulate the G. uralensis development as far as its increased bioactive compound levels, and plant yield are concerned. In this study, leaf nutrient concentration and photosynthetic chlorophyll index of G. uralensis were evaluated to determine the G. uralensis growth physiology in three habitats. We observed that G. uralensis nutrient levels and photosynthesis differed significantly in three regions (p < 0.05). Whole-genome re-sequencing of the sixty G. uralensis populations samples from different habitats was performed using an Illumina HiSeq sequencing platform to elucidate the distribution patterns, population evolution, and genetic diversity of G. uralensis. 150.06 Gb high-quality clean data was obtained after strict filtering. The 895237686 reads were mapped against the reference genome, with an average 89.7% mapping rate and 87.02% average sample coverage rate. A total of 6985987 SNPs were identified, and 117970 high-quality SNPs were obtained after filtering, which were subjected to subsequent analysis. Principal component analysis (PCA) based on interindividual SNPs and phylogenetic analysis based on interindividual SNPs showed that the G. uralensis samples could be categorized into central, southern, and eastern populations, which reflected strong genetic differentiation due to long periods of geographic isolation. In this study, a total of 131 candidate regions were screened, and 145 candidate genes (such as Glyur001802s00036258, Glyur003702s00044485, Glyur001802s00036257, Glyur007364s00047495, Glyur000028s00003476, and Glyur000398s00034457) were identified by selective clearance analysis based on Fst and θπ values. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed significant enrichment of 110 GO terms including carbohydrate metabolic process, carbohydrate biosynthetic process, carbohydrate derivative biosynthetic process, and glucose catabolic process (p < 0.05). Alpha-linolenic acid metabolism, biosynthesis of unsaturated fatty acids, and fatty acid degradation pathways were significantly enriched (p < 0.05). This study provides information on the genetic diversity, genetic structure, and population adaptability of the medicinal legumes, G. uralensis. The data obtained in this study provide valuable information for plant development and future optimization of breeding programs for functional genes.

Chilling tolerance in rice: Past and present

Rice is generally sensitive to chilling stress, which seriously affects growth and yield. Since early in the last century, considerable efforts have been made to understand the physiological and molecular mechanisms underlying the response to chilling stress and improve rice chilling tolerance. Here, we review the research trends and advances in this field. The phenotypic and biochemical changes caused by cold stress and the physiological explanations are briefly summarized. Using published data from the past 20 years, we reviewed the past progress and important techniques in the identification of quantitative trait loci (QTL), novel genes, and cellular pathways involved in rice chilling tolerance. The advent of novel technologies has significantly advanced studies of cold tolerance, and the characterization of QTLs, key genes, and molecular modules have sped up molecular design breeding for cold tolerance in rice varieties. In addition to gene function studies based on overexpression or artificially generated mutants, elucidating natural allelic variation in specific backgrounds is emerging as a novel approach for the study of cold tolerance in rice, and the superior alleles identified using this approach can directly facilitate breeding.

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.

Can omics deliver temperature resilient ready-to-grow crops?

Plants are extensively well-thought-out as the main source for nourishing natural life on earth. In the natural environment, plants have to face several stresses, mainly heat stress (HS), chilling stress (CS) and freezing stress (FS) due to adverse climate fluctuations. These stresses are considered as a major threat for sustainable agriculture by hindering plant growth and development, causing damage, ultimately leading to yield losses worldwide and counteracting to achieve the goal of "zero hunger" proposed by the Food and Agricultural Organization (FAO) of the United Nations. Notably, this is primarily because of the numerous inequities happening at the cellular, molecular and/or physiological levels, especially during plant developmental stages under temperature stress. Plants counter to temperature stress via a complex phenomenon including variations at different developmental stages that comprise modifications in physiological and biochemical processes, gene expression and differences in the levels of metabolites and proteins. During the last decade, omics approaches have revolutionized how plant biologists explore stress-responsive mechanisms and pathways, driven by current scientific developments. However, investigations are still required to explore numerous features of temperature stress responses in plants to create a complete idea in the arena of stress signaling. Therefore, this review highlights the recent advances in the utilization of omics approaches to understand stress adaptation and tolerance mechanisms. Additionally, how to overcome persisting knowledge gaps. Shortly, the combination of integrated omics, genome editing, and speed breeding can revolutionize modern agricultural production to feed millions worldwide in order to accomplish the goal of "zero hunger."