Genome-wide association study identifies a gene responsible for temperature-dependent rice germination

Genome-wide association analysis and transgenic characterization for amylose content regulating gene in tuber of Dioscorea zingiberensis

BACKGROUND

Amylose, a prebiotic found in yams is known to be beneficial for the gut microflora and is particularly advantageous for diabetic patients' diet. However, the genetic machinery underlying amylose production remains elusive. A comprehensive characterization of the genetic basis of amylose content in yam tubers is a prerequisite for accelerating the genetic engineering of yams with respect to amylose content variation.

RESULTS

To uncover the genetic variants underlying variation in amylose content, we evaluated amylose content in freshly harvested tubers from 150 accessions of Dioscorea zingibensis. With 30,000 high-quality single nucleotide polymorphisms (SNP), we performed a genome-wide association analysis (GWAS). The population structure analysis classified the D. zingiberensis accessions into three groups. A total of 115 significant loci were detected on four chromosomes. Of these, 112 significant SNPs (log10(p) = 5, q-value < 0.004) were clustered in a narrow window on the chromosome 6 (chr6). The peak SNP at the position 75,609,202 on chr6 could explain 63.15% of amylose variation in the population and fell into the first exon of the ADP-glucose pyrophosphorylase (AGPase) small subunit gene, causing a non-synonymous modification of the resulting protein sequence. Allele segregation analysis showed that accessions with the rare G allele had a higher amylose content than those harboring the common A allele. However, AGPase, a key enzyme precursor of amylose biosynthesis, was not expressed differentially between accessions with A and G alleles. Overexpression of the two variants of AGPase in Arabidopsis thaliana resulted in a significantly higher amylose content in lines transformed with the AGPase-G allele.

CONCLUSIONS

Overall, this study showed that a major genetic variant in AGPase probably enhances the enzyme activity leading to high amylose content in D. zingiberensis tuber. The results provide valuable insights for the development of amylose-enriched genotypes.

Spatiotemporal transcriptomic landscape of rice embryonic cells during seed germination

Characterizing cellular features during seed germination is crucial for understanding the complex biological functions of different embryonic cells in regulating seed vigor and seedling establishment. We performed spatially enhanced resolution omics sequencing (Stereo-seq) and single-cell RNA sequencing (scRNA-seq) to capture spatially resolved single-cell transcriptomes of germinating rice embryos. An automated cell-segmentation model, employing deep learning, was developed to accommodate the analysis requirements. The spatial transcriptomes of 6, 24, 36, and 48 h after imbibition unveiled both known and previously unreported embryo cell types, including two unreported scutellum cell types, corroborated by in situ hybridization and functional exploration of marker genes. Temporal transcriptomic profiling delineated gene expression dynamics in distinct embryonic cell types during seed germination, highlighting key genes involved in nutrient metabolism, biosynthesis, and signaling of phytohormones, reprogrammed in a cell-type-specific manner. Our study provides a detailed spatiotemporal transcriptome of rice embryo and presents a previously undescribed methodology for exploring the roles of different embryonic cells in seed germination.

Genome resources for three modern cotton lines guide future breeding efforts

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

TCP transcription factor StAST1 represses potato tuberization by regulating tuberigen complex activity

Potato (Solanum tuberosum L.) is cultivated worldwide for its underground tubers, which provide an important part of human nutrition and serve as a model system for belowground storage organ formation. Similar to flowering, stolon-expressed FLOWERING LOCUS T-like (FT-like) protein SELF-PRUNING 6A (StSP6A) plays an instrumental role in tuberization by binding to the bZIP transcription factors StABI5-like 1 (StABL1) and StFD-like 1 (StFDL1), causing transcriptional reprogramming at the stolon subapical apices. However, the molecular mechanism regulating the widely conserved FT-bZIP interactions remains largely unexplored. Here, we identified a TCP transcription factor StAST1 (StABL1 and StSP6A-associated TCP protein 1) binding to both StSP6A and StABL1. StAST1 is specifically expressed in the vascular tissue of leaves and developing stolons. Silencing of StAST1 leads to accelerated tuberization and a shortened life cycle. Molecular dissection reveals that the interaction of StAST1 with StSP6A and StABL1 attenuates the formation of the alternative tuberigen activation complex (aTAC). We also observed StAST1 directly activates the expression of potato GA 20-oxidase gene (StGA20ox1) to regulate GA responses. These results demonstrate StAST1 functions as a tuberization repressor by regulating plant hormone levels; our findings also suggest a mechanism by which the widely conserved FT-FD genetic module is fine-tuned.

AMPK activator 991 specifically activates SnRK1 and thereby affects seed germination in rice

Sucrose non-fermenting-1-related protein kinase 1 (SnRK1) and AMP-activated protein kinase (AMPK) are highly conserved. Compound 991 is an AMPK activator in mammals. However, whether 991 also activates SnRK1 remains unknown. The addition of 991 significantly increased SnRK1 activity in desalted extracts from germinating rice seeds in vitro. To determine whether 991 has biological activity, rice seeds were treated with different concentrations of 991. Germination was promoted at low concentrations but inhibited at high concentrations. The effects of 991 on germination were similar to those of OsSnRK1a overexpression. To explore whether 991 affects germination by specifically affecting SnRK1, germination of an snrk1a mutant and the wild type under 1 μM 991 treatment was compared. The snrk1a mutant was insensitive to 991. Phosphoproteomic analysis showed that the differential phosphopeptides induced by 991 and OsSnRK1a overexpression largely overlapped. Furthermore, SnRK1 might regulate rice germination in a dosage-dependent manner by regulating the phosphorylation of three phosphosites, namely S285-PIP2;4, S1013-SOS1, and S110-ABI5. These results indicate that 991 is a specific SnRK1 activator in rice. The promotion and inhibition of germination by 991 also occurred in wheat seeds. Thus, 991 is useful for exploring SnRK1 function and the chemical regulation of growth and development in crops.

Redefining awn development in rice through the breeding history of Japanese awn reduction

The study challenges the conventional understanding of awn loss as a domestication syndrome, showing instead that many awned varieties continued to be widely grown in Japan until the early twentieth century and that selection for awn reduction was active at that time, demonstrating that awn loss is not a domestication syndrome but "a trait that emerged during crop improvement". Although selection for awnless mutants was carried out independently using different types of awned cultivars in the early twentieth century in Japan, awn loss was caused by the mutation in OsEPFL1. This suggests that a single mutant haplotype of OsEPFL1 was conserved in the genomes of different cultivars and subsequently selected within each line to meet the demand for awnless varieties. The study also conducts phylogenetic analyses of EPFL1 in 48 grass plants, revealing its unique involvement in awn formation in rice while potentially playing a different role in the domestication of other grass plants. Finally, an attempt is made to isolate an awn-forming gene that has not been identified from the awned rice cultivar "Omachi", which is still cultivated in Japan. The results presented in this paper provide a new perspective on domestication against the conventional understanding of awn development, shedding light on its potential as a useful organ for breeding to mitigate environmental stress.

The Function of Florigen in the Vegetative-to-Reproductive Phase Transition in and around the Shoot Apical Meristem

Plants undergo a series of developmental phases throughout their life-cycle, each characterized by specific processes. Three critical features distinguish these phases: the arrangement of primordia (phyllotaxis), the timing of their differentiation (plastochron) and the characteristics of the lateral organs and axillary meristems. Identifying the unique molecular features of each phase, determining the molecular triggers that cause transitions and understanding the molecular mechanisms underlying these transitions are keys to gleaning a complete understanding of plant development. During the vegetative phase, the shoot apical meristem (SAM) facilitates continuous leaf and stem formation, with leaf development as the hallmark. The transition to the reproductive phase induces significant changes in these processes, driven mainly by the protein FT (FLOWERING LOCUS T) in Arabidopsis and proteins encoded by FT orthologs, which are specified as 'florigen'. These proteins are synthesized in leaves and transported to the SAM, and act as the primary flowering signal, although its impact varies among species. Within the SAM, florigen integrates with other signals, culminating in developmental changes. This review explores the central question of how florigen induces developmental phase transition in the SAM. Future research may combine phase transition studies, potentially revealing the florigen-induced developmental phase transition in the SAM.

Low mutation rate of spontaneous mutants enables detection of causative genes by comparing whole genome sequences

In the early 1900s, mutation breeding to select varieties with desirable traits using spontaneous mutation was actively conducted around the world, including Japan. In rice, the number of fixed mutations per generation was estimated to be 1.38-2.25. Although this low mutation rate was a major problem for breeding in those days, in the modern era with the development of next-generation sequencing (NGS) technology, it was conversely considered to be an advantage for efficient gene identification. In this paper, we proposed an in silico approach using NGS to compare the whole genome sequence of a spontaneous mutant with that of a closely related strain with a nearly identical genome, to find polymorphisms that differ between them, and to identify the causal gene by predicting the functional variation of the gene caused by the polymorphism. Using this approach, we found four causal genes for the dwarf mutation, the round shape grain mutation and the awnless mutation. Three of these genes were the same as those previously reported, but one was a novel gene involved in awn formation. The novel gene was isolated from Bozu-Aikoku, a mutant of Aikoku with the awnless trait, in which nine polymorphisms were predicted to alter gene function by their whole-genome comparison. Based on the information on gene function and tissue-specific expression patterns of these candidate genes, Os03g0115700/LOC_Os03g02460, annotated as a short-chain dehydrogenase/reductase SDR family protein, is most likely to be involved in the awnless mutation. Indeed, complementation tests by transformation showed that it is involved in awn formation. Thus, this method is an effective way to accelerate genome breeding of various crop species by enabling the identification of useful genes that can be used for crop breeding with minimal effort for NGS analysis.

GWAS unravels acid phosphatase ACP2 as a photosynthesis regulator under phosphate starvation conditions through modulating serine metabolism in rice

Inorganic phosphorus (Pi) deficiency significantly impacts plant growth, development, and photosynthetic efficiency. This study evaluated 206 rice accessions from a MiniCore population under both Pi-sufficient (Pi+) and Pi-starvation (Pi-) conditions in the field to assess photosynthetic phosphorus use efficiency (PPUE), defined as the ratio of AsatPi- to AsatPi+. A genome-wide association study and differential gene expression analyses identified an acid phosphatase gene (ACP2) that responds strongly to phosphate availability. Overexpression and knockout of ACP2 led to a 67% increase and 32% decrease in PPUE, respectively, compared with wild type. Introduction of an elite allele A, by substituting the v5 SNP G with A, resulted in an 18% increase in PPUE in gene-edited ACP2 rice lines. The phosphate-responsive gene PHR2 was found to transcriptionally activate ACP2 in parallel with PHR2 overexpression, resulting in an 11% increase in PPUE. Biochemical assays indicated that ACP2 primarily catalyzes the hydrolysis of phosphoethanolamine and phospho-L-serine. In addition, serine levels increased significantly in the ACP2v8G-overexpression line, along with a concomitant decrease in the expression of all nine genes involved in the photorespiratory pathway. Application of serine enhanced PPUE and reduced photorespiration rates in ACP2 mutants under Pi-starvation conditions. We deduce that ACP2 plays a crucial role in promoting photosynthesis adaptation to Pi starvation by regulating serine metabolism in rice.

A transposon insertion in the promoter of OsUBC12 enhances cold tolerance during japonica rice germination

Low-temperature germination (LTG) is an important agronomic trait for rice (Oryza sativa). Japonica rice generally has greater capacity for germination at low temperatures than the indica subpopulation. However, the genetic basis and molecular mechanisms underlying this complex trait are poorly understood. Here, we report that OsUBC12, encoding an E2 ubiquitin-conjugating enzyme, increases low-temperature germinability in japonica, owing to a transposon insertion in its promoter enhancing its expression. Natural variation analysis reveals that transposon insertion in the OsUBC12 promoter mainly occurs in the japonica lineage. The variation detected in eight representative two-line male sterile lines suggests the existence of this allele introgression by indica-japonica hybridization breeding, and varieties carrying the japonica OsUBC12 locus (transposon insertion) have higher low-temperature germinability than varieties without the locus. Further molecular analysis shows that OsUBC12 negatively regulate ABA signaling. OsUBC12-regulated seed germination and ABA signaling mainly depend on a conserved active site required for ubiquitin-conjugating enzyme activity. Furthermore, OsUBC12 directly associates with rice SUCROSE NON-FERMENTING 1-RELATED PROTEIN KINASE 1.1 (OsSnRK1.1), promoting its degradation. OsSnRK1.1 inhibits LTG by enhancing ABA signaling and acts downstream of OsUBC12. These findings shed light on the underlying mechanisms of UBC12 regulating LTG and provide genetic reference points for improving LTG in indica rice.

Comparative quantitative trait loci analysis framework reveals relationships between salt stress responsive phenotypes and pathways

Soil salinity is a complex abiotic stress that involves several biological pathways. Hence, focusing on a specific or a few salt-tolerant phenotypes is unlikely to provide comprehensive insights into the intricate and interwinding mechanisms that regulate salt responsiveness. In this study, we develop a heuristic framework for systematically integrating and comprehensively evaluating quantitative trait loci (QTL) analyses from multiple stress-related traits obtained by different studies. Making use of a combined set of 46 salinity-related traits from three independent studies that were based on the same chromosome segment substitution line (CSSL) population of rice (Oryza sativa), we demonstrate how our approach can address technical biases and limitations from different QTL studies and calling methods. This allows us to compile a comprehensive list of trait-specific and multi-trait QTLs, as well as salinity-related candidate genes. In doing so, we discover several novel relationships between traits that demonstrate similar trends of phenotype scores across the CSSLs, as well as the similarities between genomic locations that the traits were mapped to. Finally, we experimentally validate our findings by expression analyses and functional validations of several selected candidate genes from multiple pathways in rice and Arabidopsis orthologous genes, including OsKS7 (ENT-KAURENE SYNTHASE 7), OsNUC1 (NUCLEOLIN 1) and OsFRO1 (FERRIC REDUCTASE OXIDASE 1) to name a few. This work not only introduces a novel approach for conducting comparative analyses of multiple QTLs, but also provides a list of candidate genes and testable hypotheses for salinity-related mechanisms across several biological pathways.

A mediator of OsbZIP46 deactivation and degradation negatively regulates seed dormancy in rice

Preharvest sprouting (PHS) is a deleterious phenotype that occurs frequently in rice-growing regions where the temperature and precipitation are high. It negatively affects yield, quality, and downstream grain processing. Seed dormancy is a trait related to PHS. Longer seed dormancy is preferred for rice production as it can prevent PHS. Here, we map QTLs associated with rice seed dormancy and clone Seed Dormancy 3.1 (SDR3.1) underlying one major QTL. SDR3.1 encodes a mediator of OsbZIP46 deactivation and degradation (MODD). We show that SDR3.1 negatively regulates seed dormancy by inhibiting the transcriptional activity of ABIs. In addition, we reveal two critical amino acids of SDR3.1 that are critical for the differences in seed dormancy between the Xian/indica and Geng/japonica cultivars. Further, SDR3.1 has been artificially selected during rice domestication. We propose a two-line model for the process of rice seed dormancy domestication from wild rice to modern cultivars. We believe the candidate gene and germplasm studied in this study would be beneficial for the genetic improvement of rice seed dormancy.

Molecular evolution and interaction of 14-3-3 proteins with H+-ATPases in plant abiotic stresses

Environmental stresses severely affect plant growth and crop productivity. Regulated by 14-3-3 proteins (14-3-3s), H+-ATPases (AHAs) are important proton pumps that can induce diverse secondary transport via channels and co-transporters for the abiotic stress response of plants. Many studies demonstrated the roles of 14-3-3s and AHAs in coordinating the processes of plant growth, phytohormone signaling, and stress responses. However, the molecular evolution of 14-3-3s and AHAs has not been summarized in parallel with evolutionary insights across multiple plant species. Here, we comprehensively review the roles of 14-3-3s and AHAs in cell signaling to enhance plant responses to diverse environmental stresses. We analyzed the molecular evolution of key proteins and functional domains that are associated with 14-3-3s and AHAs in plant growth and hormone signaling. The results revealed evolution, duplication, contraction, and expansion of 14-3-3s and AHAs in green plants. We also discussed the stress-specific expression of those 14-3-3and AHA genes in a eudicotyledon (Arabidopsis thaliana), a monocotyledon (Hordeum vulgare), and a moss (Physcomitrium patens) under abiotic stresses. We propose that 14-3-3s and AHAs respond to abiotic stresses through many important targets and signaling components of phytohormones, which could be promising to improve plant tolerance to single or multiple environmental stresses.

Integrated genome-wide differentiation and association analyses identify causal genes underlying breeding-selected grain quality traits in japonica rice

Improving grain quality is a primary objective in contemporary rice breeding. Japanese modern rice breeding has developed two different types of rice, eating and sake-brewing rice, with different grain characteristics, indicating the selection of variant gene alleles during the breeding process. Given the critical importance of promptly and efficiently identifying genes selected in past breeding for future molecular breeding, we conducted genome scans for divergence, genome-wide association studies, and map-based cloning. Consequently, we successfully identified two genes, OsMnS and OsWOX9D, both contributing to rice grain traits. OsMnS encodes a mannan synthase that increases the white core frequency in the endosperm, a desirable trait for sake brewing but decreases the grain appearance quality. OsWOX9D encodes a grass-specific homeobox-containing transcription factor, which enhances grain width for better sake brewing. Furthermore, haplotype analysis revealed that their defective alleles were selected in East Asia, but not Europe, during modern improvement. In addition, our analyses indicate that a reduction in grain mannan content during African rice domestication may also be caused a defective OsMnS allele due to breeding selection. This study not only reveals the delicate balance between grain appearance quality and nutrition in rice but also provides a new strategy for isolating causal genes underlying complex traits, based on the concept of "breeding-assisted genomics" in plants.

OsMFT1 Inhibits Seed Germination by Modulating Abscisic Acid Signaling and Gibberellin Biosynthesis under Salt Stress in Rice

Seed dormancy and germination are regulated by endogenous gene expression as well as hormonal and environmental conditions, such as salinity, which greatly inhibits seed germination. MOTHER OF FT AND TFL1 (MFT), which encodes a phosphatidylethanolamine-binding protein, is a key regulator of seed germination in Arabidopsis thaliana. There are two orthologous genes of AtMFT in rice (Oryza sativa), namely, OsMFT1 and OsMFT2. However, the functions of these two genes in regulating rice seed germination under salt stress remain unknown. In this study, we found that seeds of loss-of-function osmft1 mutants germinated faster than wild-type (WT) seeds under salt stress, but this was not the case for loss-of-function osmft2 mutants. Overexpression of OsMFT1 (OsMFT1OE) or OsMFT2 increased the sensitivity to salt stress during seed germination. Transcriptome comparisons of osmft1 vs WT in the absence and presence of salt stress yielded several differentially expressed genes, which were associated with salt stress, plant hormone metabolism and signaling pathways, such as B-BOX ZINC FINGER 6, O. sativa bZIP PROTEIN 8 and GIBBERELLIN (GA) 20-oxidase 1. In addition, the sensitivity of OsMFT1OE seeds to GA and osmft1 seeds to abscisic acid (ABA) during seed germination increased under salt stress. Overall, our results indicate that ABA and GA metabolism and their signaling pathways are regulated by OsMFT1, modulating seed germination in rice under salt stress.

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.

Transcriptomic and QTL Analysis of Seed Germination Vigor under Low Temperature in Weedy Rice WR04-6

Low temperature is one of the major factors affecting rice germination, and low temperature germination (LTG) is an important agronomic trait. Although significant progress has been made in the study of rice LTG, the molecular mechanism of LTG remains poorly understood. To explore more rice LTG gene resources, we first demonstrated that weedy rice WR04-6 (Oryza sativa f. spontanea) had significantly higher LTG ability at 10 °C than the cultivated rice Qishanzhan (QSZ Oryza sativa L. ssp. indica). RNA-seq was used to investigate the gene expression of WR04-6 and QSZ at 10 °C for 10, 12 and 14 days after imbibition (DAI) of seed germination. The results of Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed that the differentially expressed genes (DEGs) between WR04-6 and QSZ were mainly concentrated on the response to starch catabolic processes and the response to abscisic acid (ABA). This is consistent with the results of α-amylase activity, ABA and gibberellins (GA) treatment. A recombinant inbred line (RIL) population derived from a cross between WR04-6 and QSZ and its high-density SNP genetic map were used to detect quantitative trait loci (QTL) for LTG rates. The results showed that two new QTLs were located on chromosome 3 and chromosome 12. Combined with the mapped QTLs and RNA-seq DEGs, sixteen candidate genes potentially associated with LTG were identified. Validation of the expression of the candidates by qRT-PCR were consistent with the RNA-seq data. These results will enable us to understand the genetic basis of LTG in weedy rice and provide new genetic resources for the generation of rice germplasm with improved LTG.