UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice

Genome-wide association study uncovers major genetic loci associated with seed flooding tolerance in soybean

BACKGROUND

Seed flooding stress is one of the threatening environmental stressors that adversely limits soybean at the germination stage across the globe. The knowledge on the genetic basis underlying seed-flooding tolerance is limited. Therefore, we performed a genome-wide association study (GWAS) using 34,718 single nucleotide polymorphism (SNPs) in a panel of 243 worldwide soybean collections to identify genetic loci linked to soybean seed flooding tolerance at the germination stage.

RESULTS

In the present study, GWAS was performed with two contrasting models, Mixed Linear Model (MLM) and Multi-Locus Random-SNP-Effect Mixed Linear Model (mrMLM) to identify significant SNPs associated with electrical conductivity (EC), germination rate (GR), shoot length (ShL), and root length (RL) traits at germination stage in soybean. With MLM, a total of 20, 40, 4, and 9 SNPs associated with EC, GR, ShL and RL, respectively, whereas in the same order mrMLM detected 27, 17, 13, and 18 SNPs. Among these SNPs, two major SNPs, Gm_08_11971416, and Gm_08_46239716 were found to be consistently connected with seed-flooding tolerance related traits, namely EC and GR across two environments. We also detected two SNPs, Gm_05_1000479 and Gm_01_53535790 linked to ShL and RL, respectively. Based on Gene Ontology enrichment analysis, gene functional annotations, and protein-protein interaction network analysis, we predicted eight candidate genes and three hub genes within the regions of the four SNPs with Cis-elements in promoter regions which may be involved in seed-flooding tolerance in soybeans and these warrant further screening and functional validation.

CONCLUSIONS

Our findings demonstrate that GWAS based on high-density SNP markers is an efficient approach to dissect the genetic basis of complex traits and identify candidate genes in soybean. The trait associated SNPs could be used for genetic improvement in soybean breeding programs. The candidate genes could help researchers better understand the molecular mechanisms underlying seed-flooding stress tolerance in soybean.

A rice QTL GS3.1 regulates grain size through metabolic-flux distribution between flavonoid and lignin metabolons without affecting stress tolerance

Grain size is a key component trait of grain weight and yield. Numbers of quantitative trait loci (QTLs) have been identified in various bioprocesses, but there is still little known about how metabolism-related QTLs influence grain size and yield. The current study report GS3.1, a QTL that regulates rice grain size via metabolic flux allocation between two branches of phenylpropanoid metabolism. GS3.1 encodes a MATE (multidrug and toxic compounds extrusion) transporter that regulates grain size by directing the transport of p-coumaric acid from the p-coumaric acid biosynthetic metabolon to the flavonoid biosynthetic metabolon. A natural allele of GS3.1 was identified from an African rice with enlarged grains, reduced flavonoid content and increased lignin content in the panicles. Notably, the natural allele of GS3.1 caused no alterations in other tissues and did not affect stress tolerance, revealing an ideal candidate for breeding efforts. This study uncovers insights into the regulation of grain size though metabolic-flux distribution. In this way, it supports a strategy of enhancing crop yield without introducing deleterious side effects on stress tolerance mechanisms.

Meta-QTLs, ortho-MQTLs and candidate genes for nitrogen use efficiency and root system architecture in bread wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs), ortho-MQTLs, and candidate genes (CGs) were identified for nitrogen use efficiency and root system architecture. For this purpose, 1788 QTLs were available from 24 studies published during 2006-2020. Of these, 1098 QTLs were projected onto the consensus map resulting in 118 MQTLs. The average confidence interval (CI) of MQTLs was reduced up to 8.56 folds in comparison to the average CI of QTLs. Of the 118 MQTLs, 112 were anchored to the physical map of the wheat reference genome. The physical interval of MQTLs ranged from 0.02 to 666.18 Mb with a mean of 94.36 Mb. Eighty-eight of these 112 MQTLs were verified by marker-trait associations (MTAs) identified in published genome-wide association studies (GWAS); the MQTLs that were verified using GWAS also included 9 most robust MQTLs, which are particularly useful for breeders; we call them 'Breeder's QTLs'. Some selected wheat MQTLs were further utilized for the identification of ortho-MQTLs for wheat and maize; 9 such ortho-MQTLs were available. As many as 1991 candidate genes (CGs) were also detected, which included 930 CGs with an expression level of > 2 transcripts per million in relevant organs/tissues. Among the CGs, 97 CGs with functions previously reported as important for the traits under study were selected. Based on homology analysis and expression patterns, 49 orthologues of 35 rice genes were also identified in MQTL regions. The results of the present study may prove useful for the improvement of selection strategy for yield potential, stability, and performance under N-limiting conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01085-0.

Integrative analyses of targeted metabolome and transcriptome of Isatidis Radix autotetraploids highlighted key polyploidization-responsive regulators

Background

Isatidis Radix, the root of Isatis indigotica Fort. (Chinese woad) can produce a variety of efficacious compound with medicinal properties. The tetraploid I. indigotica plants exhibit superior phenotypic traits, such as greater yield, higher bioactive compounds accumulation and enhanced stress tolerance. In this study, a comparative transcriptomic and metabolomic study on Isatidis Radix autotetraploid and its progenitor was performed.

Results

Through the targeted metabolic profiling, 283 metabolites were identified in Isatidis Radix, and 70 polyploidization-altered metabolites were obtained. Moreover, the production of lignans was significantly increased post polyploidization, which implied that polyploidization-modulated changes in lignan biosynthesis. Regarding the transcriptomic shift, 2065 differentially expressed genes (DEGs) were identified as being polyploidy-responsive genes, and the polyploidization-altered DEGs were enriched in phenylpropanoid biosynthesis and plant hormone signal transduction. The further integrative analysis of polyploidy-responsive metabolome and transcriptome showed that 1584 DEGs were highly correlated with the 70 polyploidization-altered metabolites, and the transcriptional factors TFs-lignans network highlighted 10 polyploidy-altered TFs and 17 fluctuated phenylpropanoid pathway compounds.

Conclusions

These results collectively indicated that polyploidization contributed to the high content of active compounds in autotetraploid roots, and the gene–lignan pathway network analysis highlighted polyploidy–responsive key functional genes and regulators.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07980-w.

Verification and dissection of one quantitative trait locus for grain size and weight on chromosome 1 in rice

Grain size and weight are the key traits determining rice quality and yield and are mainly controlled by quantitative trait loci (QTL). In this study, one minor QTL that was previously mapped in the marker interval of JD1009-JD1019 using the Huanghuazhan/Jizi1560 (HHZ/JZ1560) recombinant inbred line (RIL) population, qTGW1-2, was validated to regulate grain size and weight across four rice-growing seasons using twenty-one near isogenic line (NIL)-F₂ populations. The twenty-one populations were in two types of genetic background that were derived from the same parents HHZ and JZ1560. Twelve F₉, F₁₀ or F₁₁ NIL-F₂ populations with the sequential residual heterozygous regions covering JD1009-RM6840 were developed from one residual heterozygote (RH) in the HHZ/JZ1560 RIL population, and the remaining nine BC₃F₃, BC₃F₄ or BC₃F₅ NIL-F₂ populations with the sequential residual heterozygous regions covering JD1009-RM6840 were constructed through consecutive backcrosses to the recurrent parent HHZ followed with marker assistant selection in each generation. Based on the QTL analysis of these genetic populations, qTGW1-2 was successfully confirmed to control grain length, width and weight and further dissected into two QTLs, qTGW1-2a and qTGW1-2b, which were respectively narrowed down to the marker intervals of JD1139-JD1127 (~ 978.2-kb) and JD1121-JD1102 (~ 54.8-kb). Furthermore, the two types of NIL-F₂ populations were proved to be able to decrease the genetic background noise and increase the detection power of minor QTL. These results provided an important basis for further map-based cloning and molecular design breeding with the two QTLs in rice.

DNA Methylation and RNA-Sequencing Analysis Show Epigenetic Function During Grain Filling in Foxtail Millet (Setaria italica L.)

Grain filling is a crucial process for crop yield and quality. Certain studies already gained insight into the molecular mechanism of grain filling. However, it is unclear whether epigenetic modifications are associated with grain filling in foxtail millet. Global DNA methylation and transcriptome analysis were conducted in foxtail millet spikelets during different stages to interpret the epigenetic effects of the grain filling process. The study employed the whole-genome bisulfite deep sequencing and advanced bioinformatics to sequence and identify all DNA methylation during foxtail millet grain filling; the DNA methylation-mediated gene expression profiles and their involved gene network and biological pathway were systematically studied. One context of DNA methylation, namely, CHH methylation, was accounted for the largest percentage, and it was gradually increased during grain filling. Among all developmental stages, the methylation levels were lowest at T2, followed by T4, which mainly occurred in CHG. The distribution of differentially methylated regions (DMR) was varied in the different genetic regions for three contexts. In addition, gene expression was negatively associated with DNA methylation. Evaluation of the interconnection of the DNA methylome and transcriptome identified some stage-specific differentially expressed genes associated with the DMR at different stages compared with the T1 developmental stage, indicating the potential function of epigenetics on the expression regulation of genes related to the specific pathway at different stages of grain development. The results demonstrated that the dynamic change of DNA methylation plays a crucial function in gene regulation, revealing the potential function of epigenetics in grain development in foxtail millet.

Tissue-specific expression of GhnsLTPs identified via GWAS sophisticatedly coordinates disease and insect resistance by regulating metabolic flux redirection in cotton

Cotton (Gossypium hirsutum) is constantly attacked by pathogens and insects. The most efficient control strategy is to develop resistant varieties using broad-spectrum gene resources. Several resistance loci harboured by superior varieties have been identified through genome-wide association studies. However, the key genes and/or loci have not been functionally identified. In this study, we identified a locus significantly associated with Verticillium wilt (VW) resistance, and within a 145.5-kb linkage disequilibrium, two non-specific lipid transfer protein genes (named GhnsLTPsA10) were highly expressed under Verticillium pathogen stress. The expression of GhnsLTPsA10 significantly increased in roots upon Verticillium dahliae stress but significantly decreased in leaves under insect attack. Furthermore, GhnsLTPsA10 played antagonistic roles in positively regulating VW and Fusarium wilt resistance and negatively mediating aphid and bollworm resistance in transgenic Arabidopsis and silenced cotton. By combining transcriptomic, histological and physiological analyses, we determined that GhnsLTPsA10-mediated phenylpropanoid metabolism further affected the balance of the downstream metabolic flux of flavonoid and lignin biosynthesis. The divergent expression of GhnsLTPsA10 in roots and leaves coordinated resistance of cotton against fungal pathogens and insects via the redirection of metabolic flux. In addition, GhnsLTPsA10 contributed to reactive oxygen species accumulation. Therefore, in this study, we elucidated the novel function of GhnsLTP and the molecular association between disease resistance and insect resistance, balanced by GhnsLTPsA10. This broadens our knowledge of the biological function of GhnsLTPsA10 in crops and provides a useful locus for genetic improvement of cotton.

Genetic and Molecular Factors Determining Grain Weight in Rice

Grain weight is one of the major factors determining single plant yield production of rice and other cereal crops. Research has begun to reveal the regulatory mechanisms underlying grain weight as well as grain size, highlighting the importance of this research for plant molecular biology. The developmental trait of grain weight is affected by multiple molecular and genetic aspects that lead to dynamic changes in cell division, expansion and differentiation. Additionally, several important biological pathways contribute to grain weight, such as ubiquitination, phytohormones, G-proteins, photosynthesis, epigenetic modifications and microRNAs. Our review integrates early and more recent findings, and provides future perspectives for how a more complete understanding of grain weight can optimize strategies for improving yield production. It is surprising that the acquired wealth of knowledge has not revealed more insights into the underlying molecular mechanisms. To accelerating molecular breeding of rice and other cereals is becoming an emergent and critical task for agronomists. Lastly, we highlighted the importance of leveraging gene editing technologies as well as structural studies for future rice breeding applications.

Multistage Defense System Activated by Tetrachlorobiphenyl and its Hydroxylated and Methoxylated Derivatives in Oryza sativa

Crops can initiate various defense responses to environmental stresses. The process is often accompanied by extensive transcriptional and metabolic changes to reallocate metabolites. However, it remains unclear how organic pollutants activate the defense systems to reallocate metabolites in crops. The current study demonstrates that three defense systems, including the cytochrome P450s (CYP450s), glutathione S-transferases (GSTs), and phenylpropanoid biosynthesis, were sequentially activated after Oryza sativa was exposed to 2,3,4,5-tetrachlorobipheny l (PCB 61) and its derivatives 4'-hydroxy-2,3,4,5-tetrachlorobiphenyl (OH-PCB 61) and 4'-methoxy-2,3,4,5-tetrachlorobiphenyl (MeO-PCB 61), respectively. Genes encoding CYP76Ms and CYP72As were significantly upregulated after 0.5 h of exposure, followed by the GST-coding gene GSTU48, suggesting that the biotransformation and detoxification of PCB 61, OH-PCB 61, and MeO-PCB 61 occurred. Subsequently, CCR1 and CCR10 involved in phenylpropanoid biosynthesis were activated after 12 h, potentially reducing the oxidative stress induced by PCB 61 and its derivatives. Furthermore, β-d-glucan exohydrolase involved in both phenylpropanoid biosynthesis and starch and sucrose metabolism was significantly downregulated by 7.04-fold in the OH-PCB 61-treated group, potentially contributing to the inhibition of sugar hydrolysis. These findings provide insights into increasing rice adaptability to organic pollutants by reinforcing the enzyme-mediated defense systems, characterizing a novel and critical strategy that enables augmented crop outputs and quality in environments stressed by organic contaminants.

Drought Resistance in Qingke Involves a Reprogramming of the Phenylpropanoid Pathway and UDP-Glucosyltransferase Regulation of Abiotic Stress Tolerance Targeting Flavonoid Biosynthesis

Tibetan hulless barley (qingke) is an important food crop in the Tibetan plateau. However, it often suffers from drought stress resulting in reduction of food production because of the extreme plateau environment. To elucidate the molecular mechanisms underlying the drought resistance of qingke, the transcriptomic and metabolomic responses of drought-sensitive (D) and drought-resistant (XL) accessions were characterized in experiments with a time course design. The phenylpropanoid pathway was reprogrammed by downregulating the lignin pathway and increasing the biosynthesis of flavonoids and anthocyanins, and this regulation improved plant tolerance for drought stress. Besides, flavonoid glycosides have induced accumulation of metabolites that participated in drought stress resistance. HVUL7H11410 exhibited the activity of wide-spectrum glucosyltransferase and mediated flavonoid glycosylation to enhance drought stress resistance. Overall, the findings provide insights into the regulatory mechanism underlying drought stress tolerance associated with metabolic reprogramming. Furthermore, the flavonoid-enriched qingke is more tolerant to drought stress and can be used as a functional food to benefit human health.

Transcriptional Regulatory Networks Associate with Early Stages of Potato Virus X Infection of Solanum tuberosum

Potato virus X (PVX) belongs to genus Potexvirus. This study characterizes the cellular transcriptome responses to PVX infection in Russet potato at 2 and 3 days post infection (dpi). Among the 1242 differentially expressed genes (DEGs), 268 genes were upregulated, and 37 genes were downregulated at 2 dpi while 677 genes were upregulated, and 265 genes were downregulated at 3 dpi. DEGs related to signal transduction, stress response, and redox processes. Key stress related transcription factors were identified. Twenty-five pathogen resistance gene analogs linked to effector triggered immunity or pathogen-associated molecular pattern (PAMP)-triggered immunity were identified. Comparative analysis with Arabidopsis unfolded protein response (UPR) induced DEGs revealed genes associated with UPR and plasmodesmata transport that are likely needed to establish infection. In conclusion, this study provides an insight on major transcriptional regulatory networked involved in early response to PVX infection and establishment.

Oryza sativa cv. Nipponbare and Oryza barthii as Unexpected Tolerance and Susceptibility Sources Against Schizotetranychus oryzae (Acari: Tetranychidae) Mite Infestation

Cultivated rice (Oryza sativa L.) is frequently exposed to multiple stresses, including Schizotetranychus oryzae mite infestation. Rice domestication has narrowed the genetic diversity of the species, leading to a wide susceptibility. This work aimed to analyze the response of two African rice species (Oryza barthii and Oryza glaberrima), weedy rice (O. sativa f. spontanea), and O. sativa cv. Nipponbare to S. oryzae infestation. Surprisingly, leaf damage, histochemistry, and chlorophyll concentration/fluorescence indicated that the African species present a higher level of leaf damage, increased accumulation of H₂O₂, and lower photosynthetic capacity when compared to O. sativa plants under infested conditions. Infestation decreased tiller number, except in Nipponbare, and caused the death of O. barthii and O. glaberrima plants during the reproductive stage. While infestation did not affect the weight of 1,000 grains in both O. sativa, the number of panicles per plant was affected only in O. sativa f. spontanea, and the percentage of full seeds per panicle and seed length were increased only in Nipponbare. Using proteomic analysis, we identified 195 differentially abundant proteins when comparing susceptible (O. barthii) and tolerant (Nipponbare) plants under control and infested conditions. O. barthii presents a less abundant antioxidant arsenal and is unable to modulate proteins involved in general metabolism and energy production under infested condition. Nipponbare presents high abundance of detoxification-related proteins, general metabolic processes, and energy production, suggesting that the primary metabolism is maintained more active compared to O. barthii under infested condition. Also, under infested conditions, Nipponbare presents higher levels of proline and a greater abundance of defense-related proteins, such as osmotin, ricin B-like lectin, and protease inhibitors (PIs). These differentially abundant proteins can be used as biotechnological tools in breeding programs aiming at increased tolerance to mite infestation.

Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions

Phenylpropanoid metabolism is one of the most important metabolisms in plants, yielding more than 8,000 metabolites contributing to plant development and plant-environment interplay. Phenylpropanoid metabolism materialized during the evolution of early freshwater algae that were initiating terrestrialization and land plants have evolved multiple branches of this pathway, which give rise to metabolites including lignin, flavonoids, lignans, phenylpropanoid esters, hydroxycinnamic acid amides, and sporopollenin. Recent studies have revealed that many factors participate in the regulation of phenylpropanoid metabolism, and modulate phenylpropanoid homeostasis when plants undergo successive developmental processes and are subjected to stressful environments. In this review, we summarize recent progress on elucidating the contribution of phenylpropanoid metabolism to the coordination of plant development and plant-environment interaction, and metabolic flux redirection among diverse metabolic routes. In addition, our review focuses on the regulation of phenylpropanoid metabolism at the transcriptional, post-transcriptional, post-translational, and epigenetic levels, and in response to phytohormones and biotic and abiotic stresses.

A moso bamboo transcription factor, Phehdz1, positively regulates the drought stress response of transgenic rice

78 HD-Zip family genes in Phyllostachys edulis were analyzed. Overexpression of Phehdz1 can improve the drought tolerance of transgenic rice and affect its secondary metabolism. Many studies suggested homeodomain-leucine zipper (HD-Zip) transcription factors are important regulators of plant growth and development, signal transduction, and responses to environmental stresses. In this study, 78 moso bamboo (Phyllostachys edulis) HD-Zip genes were investigated and classified into four subfamilies (HD-Zip I-IV). Additionally, Phehdz1 (HD-Zip I gene) was isolated and confirmed to be highly expressed in the roots. A quantitative real-time PCR analysis indicated Phehdz1 expression was significantly induced by drought, high salinity, and abscisic acid (ABA). A transient expression assay proved that Phehdz1 was localized in the nucleus of tobacco cells. Moreover, it could bind to the core region encoded by the H-box sequence (CAATAATTG) in yeast. In response to mannitol treatments, the Phehdz1-overexpressing transgenic rice had a higher germination rate and longer shoots than the wild-type controls. Moreover, Phehdz1-overexpressing rice plants had a higher survival rate as well as higher relative water and proline contents, but a lower malondialdehyde content, than the WT plants after a 30% polyethylene glycol 6000 treatment. Accordingly, the overexpression of Phehdz1 enhances the drought tolerance of transgenic rice. Many of the differentially expressed genes identified by a transcriptome analysis are involved in MAPK signal transduction and the biosynthesis of secondary metabolites. Thus, the overexpression of Phehdz1 enhances the drought stress tolerance of transgenic rice, while also potentially modulating the expression of metabolism-related genes.

Rice Glycosyltransferase Gene UGT85E1 Is Involved in Drought Stress Tolerance Through Enhancing Abscisic Acid Response

Drought is one of the most important environmental constraints affecting plant growth and development and ultimately leads to yield loss. Uridine diphosphate (UDP)-dependent glycosyltransferases (UGTs) are believed to play key roles in coping with environmental stresses. In rice, it is estimated that there are more than 200 UGT genes. However, most of them have not been identified as their physiological significance. In this study, we reported the characterization of a putative glycosyltransferase gene UGT85E1 in rice. UGT85E1 gene is significantly upregulated by drought stress and abscisic acid (ABA) treatment. The overexpression of UGT85E1 led to an enhanced tolerance in transgenic rice plants to drought stress, while the ugt85e1 mutants of rice showed a more sensitive phenotype to drought stress. Further studies indicated that UGT85E1 overexpression induced ABA accumulation, stomatal closure, enhanced reactive oxygen species (ROS) scavenging capacity, increased proline and sugar contents, and upregulated expression of stress-related genes under drought stress conditions. Moreover, when UGT85E1 was ectopically overexpressed in Arabidopsis, the transgenic plants showed increased tolerance to drought as well as in rice. Our findings suggest that UGT85E1 plays an important role in mediating plant response to drought and oxidative stresses. This work may provide a promising candidate gene for cultivating drought-tolerant crops both in dicots and monocots.

Differential Contributions of MYCs to Insect Defense Reveals Flavonoids Alleviating Growth Inhibition Caused by Wounding in Arabidopsis

In Arabidopsis, basic helix-loop-helix transcription factors (TFs) MYC2, MYC3, and MYC4 are involved in many biological processes, such as defense against insects. We found that despite functional redundancy, MYC-related mutants displayed different resistance to cotton bollworm (Helicoverpa armigera). To screen out the most likely genes involved in defense against insects, we analyzed the correlation of gene expression with cotton bollworm resistance in wild-type (WT) and MYC-related mutants. In total, the expression of 94 genes in untreated plants and 545 genes in wounded plants were strongly correlated with insect resistance, and these genes were defined as MGAIs (MYC-related genes against insects). MYC3 had the greatest impact on the total expression of MGAIs. Gene ontology (GO) analysis revealed that besides the biosynthesis pathway of glucosinolates (GLSs), MGAIs, which are well-known defense compounds, were also enriched in flavonoid biosynthesis. Moreover, MYC3 dominantly affected the gene expression of flavonoid biosynthesis. Weighted gene co-expression network analysis (WGCNA) revealed that AAE18, which is involved in activating auxin precursor 2,4-dichlorophenoxybutyric acid (2,4-DB) and two other auxin response genes, was highly co-expressed with flavonoid biosynthesis genes. With wounding treatment, the WT plants exhibited better growth performance than chalcone synthase (CHS), which was defective in flavonoid biosynthesis. The data demonstrated dominant contributions of MYC3 to cotton bollworm resistance and imply that flavonoids might alleviate the growth inhibition caused by wounding in Arabidopsis.

Identification through fine mapping and verification using CRISPR/Cas9-targeted mutagenesis for a minor QTL controlling grain weight in rice

A minor QTL for grain weight in rice, qTGW1.2b, was fine-mapped. Its casual gene OsVQ4 was confirmed through CRISPR/Cas9-targeted mutagenesis, exhibiting an effect that was larger than the original QTL effect. The CRISPR/Cas system exhibits a great potential for rice improvement, but the application was severely hindered due to insufficient target genes, especial the lack of validated genes underlying quantitative trait loci having small effects. In this study, a minor QTL for grain weight, qTGW1.2b, was fine-mapped into a 44.0 kb region using seven sets of near isogenic lines (NILs) developed from the indica rice cross (Zhenshan 97)³/Milyang 46, followed by validation of the causal gene using CRISPR/Cas9-targeted mutagenesis. In the NIL populations, 1000-grain weight of the Zhenshan 97 homozygous lines decreased by 0.9-2.0% compared with the Milyang 46 homozygous lines. A gene encoding VQ-motif protein, OsVQ4, was identified as the candidate gene based on parental sequence differences. The effect of OsVQ4 was confirmed by creating CRISPR/Cas9 knockout lines, whose 1000-grain weight decreased by 2.8-9.8% compared with the wild-type transgenic line and the recipient. These results indicate that applying genome editing system could create novel alleles with large phenotypic variation at minor QTLs, which is an effective way to validate causal genes of minor QTLs. Our study establishes a strategy for cloning minor QTLs, which could also be used to identify a large number of potential target genes for the application of CRISPR/Cas system.

Programmed Editing of Rice (Oryza sativa L.) OsSPL16 Gene Using CRISPR/Cas9 Improves Grain Yield by Modulating the Expression of Pyruvate Enzymes and Cell Cycle Proteins

Rice (Oryza sativa L.) is one of the major crops in the world and significant increase in grain yield is constant demand for breeders to meet the needs of a rapidly growing population. The size of grains is one of major components determining rice yield and a vital trait for domestication and breeding. To increase the grain size in rice, OsSPL16/qGW8 was mutagenized through CRISPR/Cas9, and proteomic analysis was performed to reveal variations triggered by mutations. More specifically, mutants were generated with two separate guide RNAs targeting recognition sites on opposite strands and genomic insertions and deletions were characterized. Mutations followed Mendelian inheritance and homozygous and heterozygous mutants lacking any T-DNA and off-target effects were screened. The mutant lines showed a significant increase in grain yield without any change in other agronomic traits in T0, T1, and T2 generations. Proteomic screening found a total of 44 differentially expressed proteins (DEPs), out of which 33 and 11 were up and downregulated, respectively. Most of the DEPs related to pyruvate kinase, pyruvate dehydrogenase, and cell division and proliferation were upregulated in the mutant plants. Pathway analysis revealed that DEPs were enriched in the biosynthesis of secondary metabolites, pyruvate metabolism, glycolysis/gluconeogenesis, carbon metabolism, ubiquinone and other terpenoid-quinone biosynthesis, and citrate cycle. Gene Ontology (GO) analysis presented that most of the DEPs were involved in the pyruvate metabolic process and pyruvate dehydrogenase complex. Proteins related to pyruvate dehydrogenase E1 component subunit alpha-1 displayed higher interaction in the protein-protein interaction (PPI) network. Thus, the overall results revealed that CRISPR/Cas9-guided OsSPL16 mutations have the potential to boost the grain yield of rice. Additionally, global proteome analysis has broad applications for discovering molecular components and dynamic regulation underlying the targeted gene mutations.