Genome-Wide Association Mapping Reveals Multiple QTLs Governing Tolerance Response for Seedling Stage Chilling Stress in Indica Rice

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.

The G protein-coupled receptor COLD1 promotes chilling tolerance in maize during germination

The geographic range and yield of the staple crop maize (Zea mays L.) are both strongly limited by low-temperature conditions. One of the most economical and effective measures for improvement of maize production is chilling tolerance enhancement. In this study, a chilling-tolerance gene in maize, ZmCOLD1, was cloned and characterized. This gene encodes a G protein-coupled receptor that is localized to the plasma membrane and the endoplasmic reticulum. A single nucleotide polymorphism (SNP) in ZmCOLD1, SNP2738, was found to confer chilling tolerance and to have promoted maize adaptations during speciation from teosinte. Overexpression of the excellent haplotype ZmCOLD1Hap11 significantly enhanced chilling tolerance, whereas knocking down ZmCOLD1 increased sensitivity to low temperatures during the germination and seedling stages. ZmCOLD1 was associated with an influx of extracellular Ca2+, increases in abscisic acid content, and decreases in gibberellic acid and indole-3-acetic acid content under low temperatures during the germination stage. ZmCOLD1 interacted with the G protein α subunit ZmCT2 at the plasma membrane, and ZmCT2 interacted with ZmLanCL in the nucleus. These proteins are components of the chilling tolerance signaling pathway in maize that are triggered by abscisic acid and photosynthesis. These results offer novel strategies for improvement of chilling tolerance in key crop species.

Genome-wide association studies in rice germplasm reveal significant genomic regions for root and yield-related traits under aerobic and irrigated conditions

The development of nutrient-use efficient rice lines is a priority amidst the changing climate and depleting resources viz., water, land, and labor for achieving sustainability in rice cultivation. Along with the traditional transplanted irrigated system of cultivation, the dry direct-seeded aerobic system is gaining ground nationwide. The root-related traits play a crucial role in nutrient acquisition, adaptation and need to be concentrated along with the yield-attributing traits. We phenotyped an association panel of 118 rice lines for seedling vigour index (SVI) traits at 14 and 21 days after sowing (DAS), root-related traits at panicle initiation (PI) stage in polythene bags under controlled aerobic condition, yield and yield-related traits under the irrigated condition at ICAR-IIRR, Hyderabad, Telangana; irrigated and aerobic conditions at ARS, Dhadesugur, Raichur, Karnataka. The panel was genotyped using simple sequence repeats (SSR) markers and genome-wide association studies were conducted for identifying marker-trait associations (MTAs). Significant correlations were recorded for root length, root dry weight with SVI, root volume at the PI stage, number of productive tillers per plant, spikelet fertility, the total number of grains per panicle with grain yield per plant under irrigated conditions, and the total number of grains per panicle with grain yield per plant under aerobic condition. The panel was divided into three sub-groups (K = 3) and correlated with the principal component analysis. The maximum number of MTAs were found on chromosomes 2, 3, and 12 with considerable phenotypic variability. Consistent MTAs were recorded for SVI traits at 14 and 21 DAS (RM25310, RM80, RM22961, RM1385), yield traits under irrigated conditions (RM2584, RM5179, RM410, RM20698, RM14753) across years at ICAR-IIRR, grain yield per plant (RM22961, RM1146) under the aerobic condition, grain yield per plant at irrigated ICAR-IIRR and SVI (RM5501), root traits at PI stage (RM2584, RM80, RM410, RM1146, RM18472). Functionally relevant genes near the MTAs through in-silico expression analysis in root and panicle tissues viz., HBF2 bZIP transcription factor, WD40 repeat-like domain, OsPILS6a auxin efflux carrier, WRKY108, OsSCP42, OsMADS80, nodulin-like domain-containing protein, amino acid transporter using various rice expression databases were identified. The identified MTAs and rice lines having high SVI traits (Langphou, TI-128, Mouli, TI-124, JBB-631-1), high yield under aerobic (Phouren, NPK-43, JBB-684, Ratnamudi, TI-112), irrigated conditions (KR-209, KR-262, Phouren, Keibi-Phou, TI-17), robust root traits like root length (MoirangPhou-Angouba, Wangoo-Phou, JBB-661, Dissi, NPK-45), root volume (Ratnachudi, KJ-221, Mow, Heimang-Phou, PUP-229) can be further employed in breeding programs for the targeted environments aimed at improving seedling vigour, yield-related traits under irrigated condition, aerobic condition as adaptability to water-saving technology.

Transcriptomic Profiling of Cold Stress-Induced Differentially Expressed Genes in Seedling Stage of Indica Rice

Cold stress significantly constrains the growth, development, productivity, and distribution of rice, particularly the indica cultivar, known for its susceptibility to cold, limiting its cultivation to specific regions. This study investigated the genes associated with cold responsiveness in the roots of two indica cultivars, SQSL (cold-tolerant) and XZX45 (cold-susceptible), through transcriptome dynamics analysis during the seedling stage. The analysis identified 8144 and 6427 differentially expressed genes (DEGs) in XZX45 and SQSL, respectively. Among these DEGs, 4672 (G2) were shared by both cultivars, while 3472 DEGs (G1) were specific to XZX45, and 1755 DEGs (G3) were specific to SQSL. Additionally, 572 differentially expressed transcription factors (TFs) from 48 TF families, including WRKY, NAC, bHLH, ERF, bZIP, MYB, C2H2, and GRAS, were identified. Gene Ontology (GO) enrichment analysis revealed significant enrichment of DEGs in the G3 group, particularly in the "response to cold" category, highlighting the crucial role of these specific genes in response to cold stress in SQSL. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated pronounced enrichment of DEGs in the G3 group in metabolic pathways such as "Pyruvate metabolism", "Glycolysis/Gluconeogenesis", and "Starch and sucrose metabolism", contributing to cold tolerance mechanisms in SQSL. Overall, this study provides comprehensive insights into the molecular mechanisms underlying cold responses in the indica cultivar, informing future genetic improvement strategies to enhance cold tolerance in susceptible indica rice cultivars.

Transfer of Stress Resilient QTLs and Panicle Traits into the Rice Variety, Reeta through Classical and Marker-Assisted Breeding Approaches

Reeta is a popular late-maturing high-yielding rice variety recommended for cultivation in the eastern Indian states. The cultivar is highly sensitive to submergence stress. Phosphorus deficiency is an additional constraint for realizing high yield. The quantitative trait loci (QTLs), Sub1, for submergence and Pup1 for low phosphorus stress tolerance along with narrow-grained trait, GW5 were introgressed into the variety from the donor parent, Swarna-Sub1 through marker-assisted breeding. In addition, phenotypic selections for higher panicle weight, grain number, and spikelet fertility were performed in each segregating generation. Foreground selection detected the 3 target QTLs in 9, 8 and 7 progenies in the BC₁F₁, BC₂F₁, and BC₃F₁ generation, respectively. Recurrent parent's genome recovery was analyzed using 168 SSR polymorphic markers. The foreground analysis in 452 BC₃F₂ progenies showed five pyramided lines in homozygous condition for the target QTLs. No donor fragment drag was noticed in the Sub1 and GW5 QTLs carrier while a segmentwas observed in the Pup1 carrier chromosome. The developed lines were higher yielding, had submergence, and had low phosphorus stress-tolerance alongwith similar to the recipient parent in the studied morpho-quality traits. A promising pyramided line is released in the name of Reeta-Panidhan (CR Dhan 413) for the flood-prone areas of Odisha state.

Mapping the Genomic Regions Controlling Germination Rate and Early Seedling Growth Parameters in Rice

Seed vigor is the key performance parameter of good quality seed. A panel was prepared by shortlisting genotypes from all the phenotypic groups representing seedling growth parameters from a total of 278 germplasm lines. A wide variation was observed for the traits in the population. The panel was classified into four genetic structure groups. Fixation indices indicated the existence of linkage disequilibrium in the population. A moderate to high level of diversity parameters was assessed using 143 SSR markers. Principal component, coordinate, neighbor-joining tree and cluster analyses showed subpopulations with a fair degree of correspondence with the growth parameters. Marker-trait association analysis detected eight novel QTLs, namely qAGR4.1, qAGR6.1, qAGR6.2 and qAGR8.1 for absolute growth rate (AGR); qRSG6.1, qRSG7.1 and qRSG8.1 for relative shoot growth (RSG); and qRGR11.1 for relative growth rate (RGR), as analyzed by GLM and MLM. The reported QTL for germination rate (GR), qGR4-1, was validated in this population. Additionally, QTLs present on chromosome 6 controlling RSG and AGR at 221 cM and RSG and AGR on chromosome 8 at 27 cM were detected as genetic hotspots for the parameters. The QTLs identified in the study will be useful for improvement of the seed vigor trait in rice.

Integrated transcriptomics and metabolomics analysis reveals key regulatory network that response to cold stress in common Bean (Phaseolus vulgaris L.)

Cold temperatures can be detrimental to crop survival and productivity. Breeding progress can be improved by understanding the molecular basis of low temperature tolerance. We investigated the key routes and critical metabolites related to low temperature resistance in cold-tolerant and -sensitive common bean cultivars 120 and 093, respectively. Many potential genes and metabolites implicated in major metabolic pathways during the chilling stress response were identified through transcriptomics and metabolomics research. Under chilling stress, the expression of many genes involved in lipid, amino acid, and flavonoid metabolism, as well as metabolite accumulation increased in the two bean types. Malondialdehyde (MDA) content was lower in 120 than in 093. Regarding amino acid metabolism, 120 had a higher concentration of acidic amino acids than 093, whereas 093 had a higher concentration of basic amino acids. Methionine accumulation was clearly higher in 120 than in 093. In addition, 120 had a higher concentration of many types of flavonoids than 093. Flavonoids, methionine and malondialdehyde could be used as biomarkers of plant chilling injury. Transcriptome analysis of hormone metabolism revealed considerably greater, expression of abscisic acid (ABA), gibberellin (GA), and jasmonic acid (JA) in 093 than in 120 during chilling stress, indicating that hormone regulation modes in 093 and 120 were different. Thus, chilling stress tolerance is different between 093 and 120 possibly due to transcriptional and metabolic regulation.

COG2 negatively regulates chilling tolerance through cell wall components altered in rice

Chilling-tolerant QTL gene COG2 encoded an extensin and repressed chilling tolerance by affecting the compositions of cell wall. Rice as a major crop is susceptible to chilling stress. Chilling tolerance is a complex trait controlled by multiple quantitative trait loci (QTLs). Here, we identify a QTL gene, COG2, that negatively regulates cold tolerance at seedling stage in rice. COG2 overexpression transgenic plants are sensitive to cold, whereas knockout transgenic lines enhance chilling tolerance. Natural variation analysis shows that Hap1 is a specific haplotype in japonica/Geng rice and correlates with chilling tolerance. The SNP1 in COG2 promoter is a specific divergency and leads to the difference in the expression level of COG2 between japonica/Geng and indica/Xian cultivars. COG2 encodes a cell wall-localized extensin and affects the compositions of cell wall, including pectin and cellulose, to defense the chilling stress. The results extend the understanding of the adaptation to the environment and provide an editing target for molecular design breeding of cold tolerance in rice.

Association mapping for protein, total soluble sugars, starch, amylose and chlorophyll content in rice

BACKGROUND

Protein, starch, amylose and total soluble sugars are basic metabolites of seed that influence the eating, cooking and nutritional qualities of rice. Chlorophyll is responsible for the absorption and utilization of the light energy influencing photosynthetic efficiency in rice plant. Mapping of these traits are very important for detection of more number of robust markers for improvement of these traits through molecular breeding approaches.

RESULTS

A representative panel population was developed by including 120 germplasm lines from the initial shortlisted 274 lines for mapping of the six biochemical traits using 136 microsatellite markers through association mapping. A wide genetic variation was detected for the traits, total protein, starch, amylose, total soluble sugars, chlorophyll a, and chlorophyll b content in the population. Specific allele frequency, gene diversity, informative markers and other diversity parameters obtained from the population indicated the effectiveness of utilization of the population and markers for mapping of these traits. The fixation indices values estimated from the population indicated the existence of linkage disequilibrium for the six traits. The population genetic structure at K = 3 showed correspondence with majority of the members in each group for the six traits. The reported QTL, qProt1, qPC6.2, and qPC8.2 for protein content; qTSS8.1 for total soluble sugar; qAC1.2 for amylose content; qCH2 and qSLCHH for chlorophyll a (Chl. a) while qChl5D for chlorophyll b (Chl. b) were validated in this population. The QTL controlling total protein content qPC1.2; qTSS7.1, qTSS8.2 and qTSS12.1 for total soluble sugars; qSC2.1, qSC2.2, qSC6.1 and qSC11.1 for starch content; qAC11.1, qAC11.2 and qAC11.3 for amylose content; qChla8.1 for Chl. a content and qChlb7.1 and qChlb8.1 for Chl. b identified by both Generalized Linear Model and Mixed Linear Model were detected as novel QTL. The chromosomal regions on chromosome 8 at 234 cM for grain protein content and total soluble sugars and at 363 cM for Chl. a and Chl. b along with the position at 48 cM on chromosome 11 for starch and amylose content are genetic hot spots for these traits.

CONCLUSION

The validated, co-localized and the novel QTL detected in this study will be useful for improvement of protein, starch, amylose, total soluble sugars and chlorophyll content in rice.

Comparison of anther transcriptomes in response to cold stress at the reproductive stage between susceptible and resistant Japonica rice varieties

BACKGROUND

Rice is one of the most important cereal crops in the world but is susceptible to cold stress (CS). In this study, we carried out parallel transcriptomic analysis at the reproductive stage on the anthers of two Japonica rice varieties with contrasting CS resistance: cold susceptible Longjing11 (LJ11) and cold resistant Longjing25 (LJ25).

RESULTS

According to the obtained results, a total of 16,762 differentially expressed genes (DEGs) were identified under CS, including 7,050 and 14,531 DEGs in LJ25 and LJ11, respectively. Examining gene ontology (GO) enrichment identified 35 up- and 39 down-regulated biological process BP GO terms were significantly enriched in the two varieties, with 'response to heat' and 'response to cold' being the most enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 33 significantly enriched pathways. Only the carbon metabolism and amino acid biosynthesis pathways with down-regulated DEGs were enriched considerably in LJ11, while the plant hormone signal transduction pathway (containing 153 DEGs) was dramatically improved. Eight kinds of plant hormones were detected in the pathway, while auxin, abscisic acid (ABA), salicylic acid (SA), and ethylene (ETH) signaling pathways were found to be the top four pathways with the most DEGs. Furthermore, the protein-protein interaction (PPI) network analysis identified ten hub genes (co-expressed gene number ≥ 30), including six ABA-related genes. Various DEGs (such as OsDREB1A, OsICE1, OsMYB2, OsABF1, OsbZIP23, OsCATC, and so on) revealed distinct expression patterns among rice types when the DEGs between LJ11 and LJ25 were compared, indicating that they are likely responsible for CS resistance of rice in cold region.

CONCLUSION

Collectively, our findings provide comprehensive insights into complex molecular mechanisms of CS response and can aid in CS resistant molecular breeding of rice in cold regions.

Genome wide association mapping of yield and various desirable agronomic traits in Rice

BACKGROUND

Rice (Oryza sativa L.) is one of the staple foods worldwide. To feed the growing population, the improvement of rice cultivars is important. To make the improvement in the rice breeding program, it is imperative to understand the similarities and differences of the existing rice accessions to find out the genetic diversity. Previous studies demonstrated the existence of abundant elite genes in rice landraces. A genome-wide association study (GWAS) was performed for yield and yield related traits to find the genetic diversity.

DESIGN

Experimental study.

METHODS AND RESULTS

A total of 204 SSRs markers were used among 17 SSRs found to be located on each chromosome in the rice genome. The diversity was analyzed using different genetic characters i.e., the total number of alleles (TNA), polymorphic information content (PIC), and gene diversity by Power markers, and the values for each genetic character per marker ranged from 2 to 9, 0.332 to 0.887 and 0.423 to 0.900 respectively across the whole genome. The results of population structure identified four main groups. MTA identified several markers associated with many agronomically important traits. These results will be very useful for the selection of potential parents, recombinants, and MTAs that govern the improvements and developments of new high yielding rice varieties.

CONCLUSIONS

Analysis of diversity in germplasm is important for the improvement of cultivars in the breeding program. In the present study, the diversity was analyzed with different methods and found that enormous diversity was present in the studied rice germplasm. The structure analysis found the presence of 4 genetic groups in the existing germplasm. A total of 129 marker-trait associations (MTAs) have been found in this study.

Unraveling the genomic regions controlling the seed vigour index, root growth parameters and germination per cent in rice

High seed vigour ensures good quality seed and higher productivity. Early seedling growth parameters indicate seed vigour in rice. Seed vigour via physiological growth parameters is a complex trait controlled by many quantitative trait loci. A panel was prepared representing a population of 274 rice landraces by including genotypes from all the phenotypic groups of sixseedling stage physiological parameters including germination % for association mapping. Wide variations for the six studiedtraits were observed in the population. The population was classified into 3 genetic groups. Fixation indices indicated the presence of linkage disequilibrium in the population. The population was classified into subpopulations and each subpopulation showed correspondence with the 6 physiological traits. A total of 5 reported QTLs viz., qGP8.1 for germination % (GP); qSVII2.1, qSVII6.1 and qSVII6.2 for seed vigour index II (SVII), and qRSR11.1 for root-shoot ratio (RSR) were validated in this mapping population. In addition, 13 QTLs regulating the physiological parameters such as qSVI 11.1 for seed vigour index I; qSVI11.1 and qSVI12.1 for seed vigour index II; qRRG10.1, qRRG8.1, qRRG8.2, qRRG6.1 and qRRG4.1 for rate of root growth (RRG); qRSR2.1, qRSR3.1 and qRSR5.1 for root-shoot ratio (RSR) while qGP6.2 and qGP6.3 for germination %were identified. Additionally, co-localization or co-inheritance of QTLs, qGP8.1 and qSVI8.1 for GP and SVI-1; qGP6.2 and qRRG6.1 for GP and RRG, and qSVI11.1 and qRSR11.1 for SVI and RSR were detected. The QTLs identified in this study will be useful for improvement of seed vigour trait in rice.

Combined QTL mapping and RNA-Seq pro-filing reveal candidate genes related to low-temperature tolerance in maize

Maize (Zea mays L.) is the most important food crop in the world, with significant acreage and production across the globe. However, it is affected by low temperatures throughout its growth process, especially during germination. Therefore, it is important to identify more QTLs or genes associated with germination under low-temperature conditions. For the QTL analysis of traits related to low-temperature germination, we used a high-res genetic map of 213 lines of the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, which had 6618 bin markers. We detected 28 QTLs of eight phenotypic characteristics associated with low-temperature germination, while they explained the phenotypic contribution rate of 5.4 ~ 13.34%. Additionally, 14 overlapping QTLs produced six QTL clusters on every chromosome, except for 8 and 10. RNA-Seq found six genes related to low-temperature tolerance in these QTLs, while qRT-PCR analysis demonstrated that the expression trends of the Zm00001d045568 gene in the LT_BvsLT_M group and the CK_BvsCK_M group were highly significantly different at all four-time points (P < 0.01), and encoded the RING zinc finger protein. It was located on qRTL9-2 and qRSVI9-1 and is related to the total length and simple vitality index. These results provided potential candidate genes for further gene cloning and improving the low-temperature tolerance of maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01297-6.

Detection of Genomic Regions Controlling the Antioxidant Enzymes, Phenolic Content, and Antioxidant Activities in Rice Grain through Association Mapping

Because it is rich in antioxidant compounds, the staple food of rice provides many health benefits. Four antioxidant traits in rice grain, viz., catalase, CUPRAC, DPPH, FRAP and peroxidase, were mapped in a representative panel population containing 117 germplasm lines using 131 SSR markers through association mapping. Donor lines rich in multiple antioxidant properties were identified from the mapping population. The population was classified into three genetic groups and each group showed reasonable correspondence with the antioxidant traits. The presence of linkage disequilibrium in the population was confirmed from the estimated Fst values. A strong positive correlation of DPPH was established with TPC, FRAP and CUPRAC. A moderate to high mean gene diversity was observed in the panel population. Eleven significant marker-trait associations for antioxidant traits were mapped, namely, qACD2.1, qACD11.1 and qACD12.2 for DPPH; qCAT8.1 and qCAT11.1 for catalase; qFRAP11.1, qFRAP12.1 and qFRAP12.2 for FRAP; and qCUPRAC3.1, qCUPRAC11.1 and qCUPRA12.1 regulating CUPRAC. Co-localization of the QTLs for qACD11.1, qFRAP11.1 and qCUPRAC11.1 were detected, which may act as antioxidant hotspots regulating DPPH, FRAP and CUPRAC activities, respectively, while qACD12.2 and qFRAP12.1 remained close on the chromosome 12. These detected QTLs will be useful in antioxidant improvement programs in rice.

Phenotypic and Genotypic Characterization and Correlation Analysis of Pea (Pisum sativum L.) Diversity Panel

Phenotypic and genotypic characterization were performed to assess heritability, variability, and seed yield stability of pea genotypes used in breeding to increase the pea production area. A European pea diversity panel, including genotypes from North America, Asia, and Australia consisting of varieties, breeding lines, pea, and landraces was examined in 2019 and 2020 in Serbia and Belgium using augmented block design. The highest heritability was for thousand seed weight; the highest coefficient of variation was for seed yield. The highest positive correlation was between number of seeds per plant and number of pods per plant; the highest negative correlation was between seed yield and protein content. Hierarchical clustering separated pea germplasm based on use and type. Different Principal component analysis grouping of landraces, breeding lines, and varieties, as well as forage types and garden and dry peas, confirms that there was an apparent decrease in similarity between the genotypes, which can be explained by their different purposes. Pea breeding should be focused on traits with consistent heritability and a positive effect on seed yield when selecting high-yielding genotypes, and on allowing for more widespread use of pea in various agricultural production systems.

An Integration of MicroRNA and Transcriptome Sequencing Analysis Reveal Regulatory Roles of miRNAs in Response to Chilling Stress in Wild Rice

A chromosome single segment substitution line (CSSL) DC90, which was generated by introgressing CTS-12, a locus derived from common wild rice (Oryza rufipogon Griff.), into the 9311 (Oryza sativa L. ssp. indica) background, exhibits a chilling tolerance phenotype under chilling stress. Here, an integration of microRNA (miRNA) deep sequencing and transcriptomic sequencing analysis was performed to explore the expression profiles of miRNAs and their target genes mediated by CTS-12 under chilling stress, and to reveal the possible regulatory mechanisms of miRNAs that are involved in chilling tolerance. Integration analysis revealed that a number of differentially expressed miRNAs (DEMs) and putative target genes with different expression patterns and levels were identified in 9311 and DC90 under chilling stress. KEGG enrichment analysis revealed that the target genes that are regulated by chilling-induced miRNAs are involved in the regulation of various biological processes/pathways, including protein biosynthesis, redox process, photosynthetic process, and chloroplast development in two genotypes. CRISPR/Cas9 editing of the target genes of the key DEMs in a chilling tolerant rice variety Zhonghua 11 (ZH11) found that LOC_Os11g48020 (OsGL1-11), one of the putative target genes of osa-miR1846a/b-5p and encoding a wax synthesis protein, is correlated with a chilling stress tolerance phenotype, implying osa-miR1846a/b-5p/OsGL1-11 plays an important role in CTS-12-mediated chilling stress tolerance regulatory pathway(s). Therefore, we speculate that the CTS-12 may regulate the key miRNA target genes in response to chilling stress by differential regulation of miRNAs in wild rice, thereby resulting in the variation of chilling tolerance phenotype between 9311 and DC90.

Population genetic structure and association mapping for iron toxicity tolerance in rice

Iron (Fe) toxicity is a major abiotic stress which severely reduces rice yield in many countries of the world. Genetic variation for this stress tolerance exists in rice germplasms. Mapping of gene(s)/QTL controlling the stress tolerance and transfer of the traits into high yielding rice varieties are essential for improvement against the stress. A panel population of 119 genotypes from 352 germplasm lines was constituted for detecting the candidate gene(s)/QTL through association mapping. STRUCTURE, GenAlEx and Darwin softwares were used to classify the population. The marker-trait association was detected by considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) analyses. Wide genetic variation was observed among the genotypes present in the panel population for the stress tolerance. Linkage disequilibrium was detected in the population for iron toxicity tolerance. The population was categorized into three genetic structure groups. Marker-trait association study considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) showed significant association of leaf browning index (LBI) with markers RM471, RM3, RM590 and RM243. Three novel QTL controlling Fe-toxicity tolerance were detected and designated as qFeTox4.3, qFeTox6.1 and qFeTox10.1. A QTL reported earlier in the marker interval of C955-C885 on chromosome 1 is validated using this panel population. The present study showed that QTL controlling Fe-toxicity tolerance to be co-localized with the QTL for Fe-biofortification of rice grain indicating involvement of common pathway for Fe toxicity tolerance and Fe content in rice grain. Fe-toxicity tolerance QTL qFeTox6.1 was co-localized with grain Fe-biofortification QTLs qFe6.1 and qFe6.2 on chromosome 6, whereas qFeTox10.1 was co-localized with qFe10.1 on chromosome 10. The Fe-toxicity tolerance QTL detected from this mapping study will be useful in marker-assisted breeding programs.

Association of molecular markers with physio-biochemical traits related to seed vigour in rice

Eighteen physio-biochemical traits influencing seed vigour were studied for their association with molecular markers using a mini core set constituted from 120 germplasm lines. High genetic variation was detected in the parameters namely chlrophyll a, Chlrophyll b, total chlorophyll, carotenoids, total anthocyanin content, gamma-oryzanols, total phenolics content, superoxide dismutase, catalase, guaicol peroxidase, total soluble sugar, total protein, seed vigour index -I and seed vigour index -II. Strong positive correlation of seed vigour index II was observed with amylose content, total anthocyanin content, catalase, total phenolic content and total flavonoid content while a negative association was observed for gamma-oryzanol content. High gene diversity (0.7169) and informative markers value (0.6789) were estimated from the investigation. Three genetic structure groups were observed in the panel population and genotypes were grouped in the subpopulations based on the seed vigour trait. Differences in the fixation indices of the three sub populations indicated existence of linkage disequilibrium in the studied panel population. Association of the traits namely total flavonoids, superoxide dismutase, catalase, chlorophyll a, Chlorophyll b, total chlorophyll, carotenoids, starch, amylose, total anthocyanin, gamma-oryzanol, total phenolics with the molecular markers were detected by Generalized Linear Model and Mixed Linear Model showing > 0.10 R² value. Association of the trait, total flavonoids with marker RM7364 located on chromosome 8 reported in earlier study was validated in this investigation. The validated markers and the novel markers detected showing higher R² value will be useful for improvement of seed vigour in rice.

Identification of QTLs for high grain yield and component traits in new plant types of rice

A panel of 60 genotypes comprising New Plant Types (NPTs) along with indica, tropical and temperate japonica genotypes was phenotypically evaluated for four seasons in irrigated situation for grain yield per se and component traits. Twenty NPT genotypes were found promising with an average grain yield varying from 5.45 to 8.8 t/ha. A total of 85 SSR markers were used in the study to identify QTLs associated with grain yield per se and related traits. Sixty-six (77.65%) markers were found to be polymorphic. The PIC values varied from 0.516 to 0.92 with an average of 0.704. A moderate level of genetic diversity (0.39) was detected among genotypes. Variation to the tune of 8% within genotypes, 68% among the genotypes within the population and 24% among the populations were observed (AMOVA). This information may help in identification of potential parents for development of transgressive segregants with very high yield. The association analysis using GLM and MLM models led to the identification of 30 and 10 SSR markers associated with 70 and 16 QTLs, respectively. Thirty novel QTLs linked with 16 SSRs were identified to be associated with eleven traits, namely tiller number (qTL-6.1, qTL-11.1, qTL-4.1), panicle length (qPL-1.1, qPL-5.1, qPL-7.1, qPL-8.1), flag leaf length (qFLL-8.1, qFLL-9.1), flag leaf width (qFLW-6.2, qFLW-5.1, qFLW-8.1, qFLW-7.1), total no. of grains (qTG-2.2, qTG-a7.1), thousand-grain weight (qTGW-a1.1, qTGW-a9.2, qTGW-5.1, qTGW-8.1), fertile grains (qFG-7.1), seed length-breadth ratio (qSlb-3.1), plant height (qPHT-6.1, qPHT-9.1), days to 50% flowering (qFD-1.1) and grain yield per se (qYLD-5.1, qYLD-6.1a, qYLD-11.1).Some of the SSRs were co-localized with more than two traits. The highest co-localization was identified with RM5709 linked to nine traits, followed by RM297 with five traits. Similarly, RM5575, RM204, RM168, RM112, RM26499 and RM22899 were also recorded to be co-localized with more than one trait and could be rated as important for marker-assisted backcross breeding programs, for pyramiding of these QTLs for important yield traits, to produce new-generation rice for prospective increment in yield potentiality and breaking yield ceiling.

Integrated RNA-Seq Analysis and Meta-QTLs Mapping Provide Insights into Cold Stress Response in Rice Seedling Roots

Rice (Oryza sativa L.) is a widely cultivated food crop around the world, especially in Asia. However, rice seedlings often suffer from cold stress, which affects their growth and yield. Here, RNA-seq analysis and Meta-QTLs mapping were performed to understand the molecular mechanisms underlying cold tolerance in the roots of 14-day-old seedlings of rice (RPY geng, cold-tolerant genotype). A total of 4779 of the differentially expressed genes (DEGs) were identified, including 2457 up-regulated and 2322 down-regulated DEGs. The GO, COG, KEEG, and Mapman enrichment results of DEGs revealed that DEGs are mainly involved in carbohydrate transport and metabolism, signal transduction mechanisms (plant hormone signal transduction), biosynthesis, transport and catabolism of secondary metabolites (phenylpropanoid biosynthesis), defense mechanisms, and large enzyme families mechanisms. Notably, the AP2/ERF-ERF, NAC, WRKY, MYB, C2H2, and bHLH transcription factors participated in rice’s cold–stress response and tolerance. On the other hand, we mapped the identified DEGs to 44 published cold–stress-related genes and 41 cold-tolerant Meta-QTLs regions. Of them, 12 DEGs were the published cold–stress-related genes and 418 DEGs fell into the cold-tolerant Meta-QTLs regions. In this study, the identified DEGs and the putative molecular regulatory network can provide insights for understanding the mechanism of cold stress tolerance in rice. In addition, DEGs in KEGG term-enriched terms or cold-tolerant Meta-QTLs will help to secure key candidate genes for further functional studies on the molecular mechanism of cold stress response in rice.

Linkage disequilibrium mapping for grain Fe and Zn enhancing QTLs useful for nutrient dense rice breeding

BACKGROUND

High yielding rice varieties are usually low in grain iron (Fe) and zinc (Zn) content. These two micronutrients are involved in many enzymatic activities, lack of which cause many disorders in human body. Bio-fortification is a cheaper and easier way to improve the content of these nutrients in rice grain.

RESULTS

A population panel was prepared representing all the phenotypic classes for grain Fe-Zn content from 485 germplasm lines. The panel was studied for genetic diversity, population structure and association mapping of grain Fe-Zn content in the milled rice. The population showed linkage disequilibrium showing deviation of Hardy-Weinberg's expectation for Fe-Zn content in rice. Population structure at K = 3 categorized the panel population into distinct sub-populations corroborating with their grain Fe-Zn content. STRUCTURE analysis revealed a common primary ancestor for each sub-population. Novel quantitative trait loci (QTLs) namely qFe3.3 and qFe7.3 for grain Fe and qZn2.2, qZn8.3 and qZn12.3 for Zn content were detected using association mapping. Four QTLs, namely qFe3.3, qFe7.3, qFe8.1 and qFe12.2 for grain Fe content were detected to be co-localized with qZn3.1, qZn7, qZn8.3 and qZn12.3 QTLs controlling grain Zn content, respectively. Additionally, some Fe-Zn controlling QTLs were co-localized with the yield component QTLs, qTBGW, OsSPL14 and qPN. The QTLs qFe1.1, qFe3.1, qFe5.1, qFe7.1, qFe8.1, qZn6, qZn7 and gRMm9-1 for grain Fe-Zn content reported in earlier studies were validated in this study.

CONCLUSION

Novel QTLs, qFe3.3 and qFe7.3 for grain Fe and qZn2.2, qZn8.3 and qZn12.3 for Zn content were detected for these two traits. Four Fe-Zn controlling QTLs and few yield component QTLs were detected to be co-localized. The QTLs, qFe1.1, qFe3.1, qFe5.1, qFe7.1, qFe8.1, qFe3.3, qFe7.3, qZn6, qZn7, qZn2.2, qZn8.3 and qZn12.3 will be useful for biofortification of the micronutrients. Simultaneous enhancement of Fe-Zn content may be possible with yield component traits in rice.

Analysis of differential gene expression in cold-tolerant vs. cold-sensitive varieties of snap bean (Phaseolus vulgaris L.) in response to low temperature stress

BACKGROUND

Snap bean, Phaseolus vulgaris L., as a warm-season vegetable, low temperature stress seriously affect the yield and quality. At present, little is known about the genes and molecular regulation mechanism in cold response in snap bean exposed to low temperature.

OBJECTIVES

Our objectives were to identify the low temperature response genes in snap bean and to examine differences in the gene response between cold-tolerant and cold-sensitive genotypes.

METHODS

We used two highly inbred snap bean lines in this study, the cold-tolerant line '120', and the cold-sensitive line '093'. The plants were grown to the three leaf and one heart stage and exposed to 4 °C low temperature. We used RNA sequencing (RNA-seq) to analyze the differences of gene expression.

RESULTS

988 and 874 cold-responsive genes were identified in 'T120 vs CK120' and 'T093 vs CK093' ('T' stands for low temperature treatment, and 'CK' stands for control at room temperature), respectively. Of these, 555 and 442 genes were unique to cold-stressed lines '120' and '093', respectively compared to the control. Our analysis of these differentially expressed genes indicates that Ca²⁺, ROS, and hormones act as signaling molecules that play important roles in low temperature response in P. vulgaris. Altering the expression of genes in these signaling pathways activates expression of downstream response genes which can interact with other signaling regulatory networks. This may maintained the balance of ROS and hormones, making line '120' more cold-tolerant than line '093'.

CONCLUSION

Our results provide a preliminarily understanding of the molecular basis of low temperature response in snap bean, and also establish a foundation for the future genetic improvement of cold sensitivity in snap bean by incorporating genes for cold tolerance.

Development of flash-flood tolerant and durable bacterial blight resistant versions of mega rice variety 'Swarna' through marker-assisted backcross breeding

Bacterial blight (BB) disease and submergence due to flash flood are the two major constraints for achieving higher yield from rainfed lowland rice. Marker-assisted backcross breeding was followed to develop submergence tolerant and durable BB resistant variety in the background of popular cultivar ‘Swarna’. Four BB resistance genes viz., Xa4, xa5, xa13, Xa21 and Sub1 QTL for submergence tolerance were incorporated into the mega variety. Foreground selection for the five target genes was performed using closely linked markers and tracked in each backcross generations. Background selection in plants carrying the target genes was performed by using 100 simple sequence repeat markers. Amongst backcross derivatives, the plant carrying five target genes and maximum recurrent parent genome content was selected in each generation and hybridized with recipient parent. Eighteen BC₃F₂ plants were obtained by selfing the selected BC₃F₁ line. Amongst the pyramided lines, 3 lines were homozygous for all the target genes. Bioassay of the 18 pyramided lines containing BB resistance genes was conducted against different Xoo strains conferred very high levels of resistance to the predominant isolates. The pyramided lines also exhibited submergence tolerance for 14 days. The pyramided lines were similar to the recurrent parent in 14 morpho-quality traits.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

A Log-Ratio Biplot Approach for Exploring Genetic Relatedness Based on Identity by State

The detection of cryptic relatedness in large population-based cohorts is of great importance in genome research. The usual approach for detecting closely related individuals is to plot allele sharing statistics, based on identity-by-state or identity-by-descent, in a two-dimensional scatterplot. This approach ignores that allele sharing data across individuals has in reality a higher dimensionality, and neither regards the compositional nature of the underlying counts of shared genotypes. In this paper we develop biplot methodology based on log-ratio principal component analysis that overcomes these restrictions. This leads to entirely new graphics that are essentially useful for exploring relatedness in genetic databases from homogeneous populations. The proposed method can be applied in an iterative manner, acting as a looking glass for more remote relationships that are harder to classify. Datasets from the 1,000 Genomes Project and the Genomes For Life-GCAT Project are used to illustrate the proposed method. The discriminatory power of the log-ratio biplot approach is compared with the classical plots in a simulation study. In a non-inbred homogeneous population the classification rate of the log-ratio principal component approach outperforms the classical graphics across the whole allele frequency spectrum, using only identity by state. In these circumstances, simulations show that with 35,000 independent bi-allelic variants, log-ratio principal component analysis, combined with discriminant analysis, can correctly classify relationships up to and including the fourth degree.

Association mapping reveals multiple QTLs for grain protein content in rice useful for biofortification

Rice is the staple food for majority of the global population. But, rice grain has low protein content (PC). Mapping of QTLs controlling grain PC is essential for enhancement of the trait through breeding programs. A shortlisted panel population for grain protein content was studied for genetic diversity, population structure and association mapping for grain PC. Phenotyping results showed a wide variation for grain PC. The panel population showed a moderate level of genetic diversity estimated through 98 molecular markers. AMOVA and structure analysis indicated linkage disequilibrium for grain PC and deviation of Hardy-Weinberg's expectation. The analysis showed 15% of the variation among populations and 73% among individuals in the panel population. STRUCTURE analysis categorized the panel population into three subpopulations. The analysis also revealed a common primary ancestor for each subpopulation with few admix individuals. Marker-trait association using 98 molecular markers detected 7 strongly associated QTLs for grain PC by both MLM and GLM analysis. Three novel QTLs qPC3.1, qPC5.1 and qPC9.1 were detected for controlling the grain PC. Four reported QTLs viz., qPC3, QPC8, qPC6.1 and qPC12.1 were validated for use in breeding programs. Reported QTLs, qPC6, qPC6.1 and qPC6.2 may be same QTL controlling PC in rice. A very close marker RM407 near to protein controlling QTL, qProt8 and qPC8, was detected. The study provided clue for simultaneous improvement of PC with high grain yield in rice. The strongly associated markers with grain PC, namely qPC3, qPC3.1, qPC5.1, qPC6.1, qPC8, qPC9.1 and qPC12.1, will be useful for their pyramiding for developing protein rich high yielding rice.

QTL Mapping in Three Connected Populations Reveals a Set of Consensus Genomic Regions for Low Temperature Germination Ability in Zea mays L

Improving seed vigor in response to cold stress is an important breeding objective in maize that allows early sowing. Using two cold tolerant inbred lines 220 and P9-10 and two susceptible lines Y1518 and PH4CV, three connected F_2:3 populations were generated for detecting quantitative trait locus (QTL) related to seed low-temperature germination ability. At 10°C, two germination traits (emergence rate and germination index) were collected from a sand bed and three seedling traits (seedling root length, shoot length, and total length) were extracted from paper rolls. Significant correlations were found among all traits in all populations. Via single-population analysis, 43 QTL were detected with explained phenotypic variance of 0.62%∼39.44%. Seventeen QTL explained more than 10% phenotypic variance; of them sixteen (94.12%) inherited favorable alleles from the tolerant lines. After constructing a consensus map, three meta-QTL (mQTL) were identified to include at least two initial QTL from different populations. mQTL1-1 included seven initial QTL for both germination and seedling traits; with three explaining more than 30% phenotypic variance. mQTL2-1 and mQTL9-1 covered two to three initial QTL. The favorable alleles of the QTL within these three mQTL regions were all inherited from the tolerant line 220 and P9-10. These results provided a basis for cloning of genes underlying the mQTL regions to uncover the molecular mechanisms of maize cold tolerance during germination.