Genetic Architecture of Natural Variation in Rice Chlorophyll Content Revealed by a Genome-Wide Association Study

Genome-wide association study and genotypic variation for the major tocopherol content in rice grain

Rice tocopherols, vitamin E compounds with antioxidant activity, play essential roles in human health. Even though the key genes involved in vitamin E biosynthetic pathways have been identified in plants, the genetic architecture of vitamin E content in rice grain remains unclear. A genome-wide association study (GWAS) on 179 genotypically diverse rice accessions with 34,323 SNP markers was conducted to detect QTLs that define total and α- tocopherol contents in rice grains. Total and α-tocopherol contents had a strong positive correlation and varied greatly across the accessions, ranging from 0.230-31.76 and 0.011-30.83 (μg/g), respectively. A total of 13 QTLs were identified, which were spread across five of the rice chromosomes. Among the 13 QTLs, 11 were considered major with phenotypic variation explained (PVE) greater than 10%. Twelve transcription factor (TF) genes, one microprotein (miP), and a transposon were found to be associated with the QTLs with putative roles in controlling tocopherol contents. Moreover, intracellular transport proteins, ABC transporters, nonaspanins, and SNARE, were identified as associated genes on chromosomes 1 and 8. In the vicinity of seven QTLs, protein kinases were identified as key signaling factors. Haplotype analysis revealed the QTLs qAlph1.1, qTot1.1, qAlph2.1, qAlph6.1, qTot6.1, and qTot8.3 to have significant haplogroups. Quantitative RT-PCR validated the expression direction and magnitude of WRKY39 (Os02g0265200), PIP5Ks (Os08g0450800), and MADS59 (Os06g0347700) in defining the major tocopherol contents. This study provides insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in rice and other cereals.

Genome-wide association study reveals genetic basis and candidate genes for chlorophyll content of leaves in maize (Zea mays L.)

The chlorophyll content (CC) directly affects photosynthesis, growth, and yield. However, the genetic basis of CC is still unclear in maize (Zea mays L.). Here, we conducted a genome-wide association study using mixed linear model for CC of the fifth leaves at seedling stage (CCFSS) and the ear leaves at filling stage (CCEFS) for 334 maize inbred lines. The heritability estimates for CCFSS and CCEFS, obtained via variance components analysis using the lme4 package in R, were 70.84% and 78.99%, respectively, indicating that the CC of leaves is primarily controlled by genetic factors. A total of 15 CC-related SNPs and 177 candidate genes were identified with a p-value < 4.49 × 10-5, which explained 4.98-7.59% of the phenotypic variation. Lines with more favorable gene variants showed higher CC. Meanwhile, Gene Ontology (GO) analysis implied that these candidate genes were probably related to chlorophyll biosynthesis. In addition, gene-based association analyses revealed that six variants in GRMZM2G037152, GRMZM5G816561, GRMZM2G324462, and GRMZM2G064657 genes were significantly (p-value < 0.01) correlated with CC, of which GRMZM2G064657 (encodes a phosphate transporter protein) and GRMZM5G816561 (encodes a cytochrome P450 protein) were specifically highly expressed in leaves tissues. Interestingly, these candidate genes were previously reported to involve in the regulation of the contents of chlorophyll in plants or Chlamydomonas. These results may contribute to the understanding of genetic basis and molecular mechanisms of maize CC and the selection of maize varieties with improved CC.

Genome-wide analysis of the CCT gene family and functional characterization of SlCCT6 in response to drought stress in tomato

CCT transcription factors are important for photoperiod and abiotic stress regulation in Arabidopsis and rice. However, the CCT gene family has not been reported in tomato. Here, we systematically analyzed this. Thirty-one SlCCT genes were identified and divided into five groups (CMF, TIFY, PRR, S8, and COL), with members unevenly distributed across 12 chromosomes and the third chromosome exhibiting the most distribution. SlCCT was found to interact with an interacting protein (SlGI), transcription factor (MYB), and non-coding RNA (sly-miR156-5p) to jointly regulate the tomato stress response. cis-Acting element analysis of the SlCCT promoter region indicated large stress- and hormone-response elements in this family. Real-time PCR results indicated that SlPRR subfamily genes respond to various abiotic stresses and hormones. Tissue expression analysis revealed that several PRR subfamily genes are highly expressed in flowers, and subcellular localization analysis indicated an SlCCT6 nuclear location. Notably, SlCCT6 expression was significantly induced by drought, and its silencing reduced drought stress tolerance. Moreover, SlCCT6 overexpression enhanced tomato drought resistance by increasing antioxidant enzyme activity and activating stress-related genes, whereas SlCCT6 knockout decreased drought resistance. In conclusion, this provides valuable insights for future research on SlCCT functions.

The atypical Dof transcriptional factor OsDes1 contributes to stay-green, grain yield, and disease resistance in rice

The regulation of leaf senescence and disease resistance plays a crucial role in determining rice grain yield and quality, which are important to meet the ever-increasing food demands of the world. Here, we identified an atypical Dof transcriptional factor OsDes1 that contributes to the stay-green phenotype, grain yield, and disease resistance in rice. The expression level of OsDes1 is positively associated with stay-green in natural variations of japonica rice, suggesting that OsDes1 would be alternatively used in breeding programs. Mechanistically, OsDes1 targets the promoter of the Rieske FeS protein gene OsPetC to activate its expression and interacts with OsPetC to protect against its degradation, thus promoting stay-green and ultimately improving the grain yield. OsDes1 also binds to the promoter region of defense-related genes, such as OsPR1b, and activates their expression, leading to enhanced disease resistance. These findings offer a potential strategy for breeding rice to enhance grain yield and disease resistance.

Genome-wide identification of the sorghum OVATE gene family and revelation of its expression characteristics in sorghum seeds and leaves

The OVATE gene family plays an important role in regulating the development of plant organs and resisting stress, but its expression characteristics and functions in sorghum have not been revealed. In this study, we identified 26 OVATE genes in the sorghum BTx623 genome, which were divided into four groups and distributed unevenly across 9 chromosomes. Evolutionary analysis showed that after differentiation between sorghum and Arabidopsis, the OVATE gene family may have experienced unique expansion events, and all OVATE family members were negatively selected. Transcriptome sequencing and RT-qPCR results showed that OVATE genes in sorghum showed diverse expression characteristics, such as gene SORBl_3001G468900 and SORBl_3009G173400 were significantly expressed in seeds, while SORBI_3005G042700 and SORBI_3002G417700 were only highly expressed in L1. Meantime, in the promoter region, a large number of hormone-associated cis-acting elements were identified, and these results suggest that members of the OVATE gene family may be involved in regulating specific development of sorghum leaves and seeds. This study improves the understanding of the OVATE gene family of sorghum and provides important clues for further exploration of the function of the OVATE gene family.

Identification of Single Nucleotide Polymorphic Loci and Candidate Genes for Seed Germination Percentage in Okra under Salt and No-Salt Stresses by Genome-Wide Association Study

Excessive soil salinity is a major stressor inhibiting crops' growth, development, and yield. Seed germination is a critical stage of crop growth and development, as well as one of the most salt-sensitive stages. Salt stress has a significant inhibitory effect on seed germination. Okra is a nutritious vegetable, but its seed germination percentage (GP) is low, whether under salt stress conditions or suitable conditions. In this study, we used 180 okra accessions and conducted a genome-wide association study (GWAS) on the germination percentage using 20,133,859 single nucleotide polymorphic (SNP) markers under 0 (CK, diluted water), 70 (treatment 1, T1), and 140 mmol/L (treatment 2, T2) NaCl conditions. Using the mixed linear model (MLM) in Efficient Mixed-model Association eXpedated (EMMAX) and Genome-wide Efficient Mixed Model Association (GEMMA) software, 511 SNP loci were significantly associated during germination, of which 167 SNP loci were detected simultaneously by both programs. Among the 167 SNPs, SNP2619493 on chromosome 59 and SNP2692266 on chromosome 44 were detected simultaneously under the CK, T1, and T2 conditions, and were key SNP loci regulating the GP of okra seeds. Linkage disequilibrium block analysis revealed that nsSNP2626294 (C/T) in Ae59G004900 was near SNP2619493, and the amino acid changes caused by nsSNP2626294 led to an increase in the phenotypic values in some okra accessions. There was an nsSNP2688406 (A/G) in Ae44G005470 near SNP2692266, and the amino acid change caused by nsSNP2688406 led to a decrease in phenotypic values in some okra accessions. These results indicate that Ae59G004900 and Ae44G005470 regulate the GP of okra seeds under salt and no-salt stresses. The gene expression analysis further demonstrated these results. The SNP markers and genes that were identified in this study will provide reference for further research on the GP of okra, as well as new genetic markers and candidate genes for cultivating new okra varieties with high GPs under salt and no-salt stress conditions.

Exploration of chlorophyll fluorescence characteristics gene regulatory in rice (Oryza sativa L.): a genome-wide association study

Chlorophyll content and fluorescence parameters are crucial indicators to evaluate the light use efficiency in rice; however, the correlations among these parameters and the underlying genetic mechanisms remain poorly understood. Here, to clarify these issues, we conducted a genome-wide association study (GWAS) on 225 rice accessions. In the phenotypic and Mendelian randomization (MR) analysis, a weak negative correlation was observed between the chlorophyll content and actual quantum yield of photosystem II (Φ I I ). The phenotypic diversity observed in SPAD, N P Q t , Φ N P Q , and F v / F m among accessions was affected by genetic background. Furthermore, the GWAS identified 78 SNPs and 17 candidate genes significantly associated with SPAD, N P Q t , Φ I I , Φ N P Q , q L and q P . Combining GWAS on 225 rice accessions with transcriptome analysis of two varieties exhibiting distinct fluorescence characteristics revealed two potential candidate genes (Os03g0583000 from Φ I I & q P traits and Os06g0587200 from N P Q t trait), which are respectively associated with peroxisomes, and protein kinase catalytic domains might involve in regulating the chlorophyll content and chlorophyll fluorescence. This study provides novel insights into the correlation among chlorophyll content and fluorescence parameters and the genetic mechanisms in rice, and offers valuable information for the breeding of rice with enhanced photosynthetic efficiency.

Genetic dissection of maize (Zea mays L.) chlorophyll content using multi-locus genome-wide association studies

BACKGROUND

The chlorophyll content (CC) is a key factor affecting maize photosynthetic efficiency and the final yield. However, its genetic basis remains unclear. The development of statistical methods has enabled researchers to design and apply various GWAS models, including MLM, MLMM, SUPER, FarmCPU, BLINK and 3VmrMLM. Comparative analysis of their results can lead to more effective mining of key genes.

RESULTS

The heritability of CC was 0.86. Six statistical models (MLM, BLINK, MLMM, FarmCPU, SUPER, and 3VmrMLM) and 1.25 million SNPs were used for the GWAS. A total of 140 quantitative trait nucleotides (QTNs) were detected, with 3VmrMLM and MLM detecting the most (118) and fewest (3) QTNs, respectively. The QTNs were associated with 481 genes and explained 0.29-10.28% of the phenotypic variation. Additionally, 10 co-located QTNs were detected by at least two different models or methods, three co-located QTNs were identified in at least two different environments, and six co-located QTNs were detected by different models or methods in different environments. Moreover, 69 candidate genes within or near these stable QTNs were screened based on the B73 (RefGen_v2) genome. GRMZM2G110408 (ZmCCS3) was identified by multiple models and in multiple environments. The functional characterization of this gene indicated the encoded protein likely contributes to chlorophyll biosynthesis. In addition, the CC differed significantly between the haplotypes of the significant QTN in this gene, and CC was higher for haplotype 1.

CONCLUSION

This study's results broaden our understanding of the genetic basis of CC, mining key genes related to CC and may be relevant for the ideotype-based breeding of new maize varieties with high photosynthetic efficiency.

Serine protease NAL1 exerts pleiotropic functions through degradation of TOPLESS-related corepressor in rice

NARROW LEAF 1 (NAL1) is a breeding-valuable pleiotropic gene that affects multiple agronomic traits in rice, although the molecular mechanism is largely unclear. Here, we report that NAL1 is a serine protease and displays a novel hexameric structure consisting of two ATP-mediated doughnut-shaped trimeric complexes. Moreover, we identified TOPLESS-related corepressor OsTPR2 involved in multiple growth and development processes as the substrate of NAL1. We found that NAL1 degraded OsTPR2, thus modulating the expression of downstream genes related to hormone signalling pathways, eventually achieving its pleiotropic physiological function. An elite allele, NAL1A, which may have originated from wild rice, could increase grain yield. Furthermore, the NAL1 homologues in different crops have a similar pleiotropic function to NAL1. Our study uncovers a NAL1-OsTPR2 regulatory module and provides gene resources for the design of high-yield crops.

Photosynthesis, Water Status and K+/Na+ Homeostasis of Buchoe dactyloides Responding to Salinity

Soil salinization is one of the most serious abiotic stresses restricting plant growth. Buffalograss is a C₄ perennial turfgrass and forage with an excellent resistance to harsh environments. To clarify the adaptative mechanisms of buffalograss in response to salinity, we investigated the effects of NaCl treatments on photosynthesis, water status and K⁺/Na⁺ homeostasis of this species, then analyzed the expression of key genes involved in these processes using the qRT-PCR method. The results showed that NaCl treatments up to 200 mM had no obvious effects on plant growth, photosynthesis and leaf hydrate status, and even substantially stimulated root activity. Furthermore, buffalograss could retain a large amount of Na⁺ in roots to restrict Na⁺ overaccumulation in shoots, and increase leaf K⁺ concentration to maintain a high K⁺/Na⁺ ratio under NaCl stresses. After 50 and 200 mM NaCl treatments, the expressions of several genes related to chlorophyll synthesis, photosynthetic electron transport and CO₂ assimilation, as well as aquaporin genes (BdPIPs and BdTIPs) were upregulated. Notably, under NaCl treatments, the increased expression of BdSOS1, BdHKT1 and BdNHX1 in roots might have helped Na⁺ exclusion by root tips, retrieval from xylem sap and accumulation in root cells, respectively; the upregulation of BdHAK5 and BdSKOR in roots likely enhanced K⁺ uptake and long-distance transport from roots to shoots, respectively. This work finds that buffalograss possesses a strong ability to sustain high photosynthetic capacity, water balance and leaf K⁺/Na⁺ homeostasis under salt stress, and lays a foundation for elucidating the molecular mechanism underlying the salt tolerance of buffalograss.

Physiological and transcriptomic analysis of OsLHCB3 knockdown lines in rice

The photosystem II (PSII) outer antenna LHCB3 protein plays critical roles in distributing the excitation energy and modulating the rate of state transition for photosynthesis. Here, OsLHCB3 knockdown mutants were produced using the RNAi system. Phenotypic analyses showed that OsLHCB3 knockdown led to pale green leaves and lower chlorophyll contents at both tillering and heading stages. In addition, mutant lines exhibited decreased non-photochemical quenching (NPQ) capacity and net photosynthetic rate (Pn) by downregulating the expression of PSII-related genes. Moreover, RNA-seq experiments were performed at both tillering and heading stages. The differentially expressed genes (DEGs) mainly involved in chlorophyll binding response to abscisic acid, photosystem II, response to chitin, and DNA-binding transcription factor. Besides, our transcriptomic and physiological data indicated that OsLHCB3 was essential for binding chlorophyll, but not for the metabolism of chlorophyll in rice. OsLHCB3 RNAi knockdown plants affected the expression of PS II-related genes, but not PS I-related genes. Overall, these results suggest that OsLHCB3 also plays vital roles in regulating photosynthesis and antenna proteins in rice as well as responses to environment stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01387-z.

Genome-Wide Association Study Identified Novel SNPs Associated with Chlorophyll Content in Maize

Chlorophyll is an essential component that captures light energy to drive photosynthesis. Chlorophyll content can affect photosynthetic activity and thus yield. Therefore, mining candidate genes of chlorophyll content will help increase maize production. Here, we performed a genome-wide association study (GWAS) on chlorophyll content and its dynamic changes in 378 maize inbred lines with extensive natural variation. Our phenotypic assessment showed that chlorophyll content and its dynamic changes were natural variations with a moderate genetic level of 0.66/0.67. A total of 19 single-nucleotide polymorphisms (SNPs) were found associated with 76 candidate genes, of which one SNP, 2376873-7-G, co-localized in chlorophyll content and area under the chlorophyll content curve (AUCCC). Zm00001d026568 and Zm00001d026569 were highly associated with SNP 2376873-7-G and encoded pentatricopeptide repeat-containing protein and chloroplastic palmitoyl-acyl carrier protein thioesterase, respectively. As expected, higher expression levels of these two genes are associated with higher chlorophyll contents. These results provide a certain experimental basis for discovering the candidate genes of chlorophyll content and finally provide new insights for cultivating high-yield and excellent maize suitable for planting environment.

Genome-wide association study identifies a gene conferring high physiological phosphorus use efficiency in rice

Phosphate (Pi) is indispensable for the growth and development of plant, and low-Pi stress is a major limitation for crop growth and yield worldwide. The tolerance to low-Pi stress varied among rice germplasm resources. However, the mechanisms underlying the tolerance of rice to low-Pi stress, as a complex quantitative trait, are not clear. We performed a genome-wide association study (GWAS) through a diverse worldwide collection of 191 rice accessions in the field under normal-Pi and low-Pi supply in two years. Twenty and three significant association loci were identified for biomass and grain yield per plant under low-Pi supply respectively. The expression level of OsAAD as a candidate gene from a associated locus was significantly up-regulated after low-Pi stress treatment for five days and tended to return to normal levels after Pi re-supply in shoots. Suppression of OsAAD expression could improve the physiological phosphorus use efficiency (PPUE) and grain yields through affecting the expression of several genes associated with GA biosynthesis and metabolism. OsAAD would be a promising gene for increasing PPUE and grain yield in rice under normal- and low-Pi supply via genome editing.

Unraveling the genetics underlying micronutrient signatures of diversity panel present in brown rice through genome-ionome linkages

Rice (Oryza sativa) is an important staple crop to address the Hidden Hunger problem not only in Asia but also in Africa where rice is fast becoming an important source of calories. The brown rice (whole grain with bran) is known to be more nutritious due to elevated mineral composition. The genetics underlying brown rice ionome (sum total of such mineral composition) remains largely unexplored. Hence, we conducted a comprehensive study to dissect the genetic architecture of the brown rice ionome. We used genome-wide association studies, gene set analysis, and targeted association analysis for 12 micronutrients in the brown rice grains. A diverse panel of 300 resequenced indica accessions, with more than 1.02 million single nucleotide polymorphisms, was used. We identified 109 candidate genes with 5-20% phenotypic variation explained for the 12 micronutrients and identified epistatic interactions with multiple micronutrients. Pooling all candidate genes per micronutrient exhibited phenotypic variation explained values ranging from 11% to almost 40%. The key donor lines with larger concentrations for most of the micronutrients possessed superior alleles, which were absent in the breeding lines. Through gene regulatory networks we identified enriched functional pathways for central regulators that were detected as key candidate genes through genome-wide association studies. This study provided important insights on the ionome variations in rice, on the genetic basis of the genome-ionome relationships and on the molecular mechanisms underlying micronutrient signatures.

Comparative analysis of physiological variations and genetic architecture for cold stress response in soybean germplasm

Soybean (Glycine max L.) is susceptible to low temperatures. Increasing lines of evidence indicate that abiotic stress-responsive genes are involved in plant low-temperature stress response. However, the involvement of photosynthesis, antioxidants and metabolites genes in low temperature response is largely unexplored in Soybean. In the current study, a genetic panel of diverse soybean varieties was analyzed for photosynthesis, chlorophyll fluorescence and leaf injury parameters under cold stress and control conditions. This helps us to identify cold tolerant (V100) and cold sensitive (V45) varieties. The V100 variety outperformed for antioxidant enzymes activities and relative expression of photosynthesis (Glyma.08G204800.1, Glyma.12G232000.1), GmSOD (GmSOD01, GmSOD08), GmPOD (GmPOD29, GmPOD47), trehalose (GmTPS01, GmTPS13) and cold marker genes (DREB1E, DREB1D, SCOF1) than V45 under cold stress. Upon cold stress, the V100 variety showed reduced accumulation of H₂O₂ and MDA levels and subsequently showed lower leaf injury compared to V45. Together, our results uncovered new avenues for identifying cold tolerant soybean varieties from a large panel. Additionally, we identified the role of antioxidants, osmo-protectants and their posttranscriptional regulators miRNAs such as miR319, miR394, miR397, and miR398 in Soybean cold stress tolerance.

Multi-model genome-wide association studies of leaf anatomical traits and vein architecture in rice

Introduction

The anatomy of rice leaves is closely related to photosynthesis and grain yield. Therefore, exploring insight into the quantitative trait loci (QTLs) and alleles related to rice flag leaf anatomical and vein traits is vital for rice improvement.

Methods

Here, we aimed to explore the genetic architecture of eight flag leaf traits using one single-locus model; mixed-linear model (MLM), and two multi-locus models; fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage disequilibrium iteratively nested keyway (BLINK). We performed multi-model GWAS using 329 rice accessions of RDP1 with 700K single-nucleotide polymorphisms (SNPs) markers.

Results

The phenotypic correlation results indicated that rice flag leaf thickness was strongly correlated with leaf mesophyll cells layer (ML) and thickness of both major and minor veins. All three models were able to identify several significant loci associated with the traits. MLM identified three non-synonymous SNPs near NARROW LEAF 1 (NAL1) in association with ML and the distance between minor veins (IVD) traits.

Discussion

Several numbers of significant SNPs associated with known gene function in leaf development and yield traits were detected by multi-model GWAS performed in this study. Our findings indicate that flag leaf traits could be improved via molecular breeding and can be one of the targets in high-yield rice development.

Improving crop yield potential: Underlying biological processes and future prospects

The growing world population and global increases in the standard of living both result in an increasing demand for food, feed and other plant-derived products. In the coming years, plant-based research will be among the major drivers ensuring food security and the expansion of the bio-based economy. Crop productivity is determined by several factors, including the available physical and agricultural resources, crop management, and the resource use efficiency, quality and intrinsic yield potential of the chosen crop. This review focuses on intrinsic yield potential, since understanding its determinants and their biological basis will allow to maximize the plant's potential in food and energy production. Yield potential is determined by a variety of complex traits that integrate strictly regulated processes and their underlying gene regulatory networks. Due to this inherent complexity, numerous potential targets have been identified that could be exploited to increase crop yield. These encompass diverse metabolic and physical processes at the cellular, organ and canopy level. We present an overview of some of the distinct biological processes considered to be crucial for yield determination that could further be exploited to improve future crop productivity.

eQTLs play critical roles in regulating gene expression and identifying key regulators in rice

The regulation of gene expression plays an essential role in both the phenotype and adaptation of plants. Transcriptome sequencing enables simultaneous identification of exonic variants and quantification of gene expression. Here, we sequenced the leaf transcriptomes of 287 rice accessions from around the world and obtained a total of 177 853 high-quality single nucleotide polymorphisms after filtering. Genome-wide association study identified 44 354 expression quantitative trait loci (eQTLs), which regulate the expression of 13 201 genes, as well as 17 local eQTL hotspots and 96 distant eQTL hotspots. Furthermore, a transcriptome-wide association study screened 21 candidate genes for starch content in the flag leaves at the heading stage. HS002 was identified as a significant distant eQTL hotspot with five downstream genes enriched for diterpene antitoxin synthesis. Co-expression analysis, eQTL analysis, and linkage mapping together demonstrated that bHLH026 acts as a key regulator to activate the expression of downstream genes. The transgenic assay revealed that bHLH026 is an important regulator of diterpenoid antitoxin synthesis and enhances the disease resistance of rice. These findings improve our knowledge of the regulatory mechanisms of gene expression variation and complex regulatory networks of the rice genome and will facilitate genetic improvement of cultivated rice varieties.

Effect of Silica-Based Nanomaterials on Seed Germination and Seedling Growth of Rice (Oryza sativa L.)

The application of nanomaterials (NMs) in agriculture has become a global concern in recent years. However, studies on their effects on plants are still limited. Here, we conducted a seed germination experiment for 5 days and a hydroponics experiment for 14 days to study the effects of silicon dioxide NMs(nSiO₂) and silicon carbide NMs(nSiC) (0,10, 50, 200 mg/L) on rice (Oryza sativa L.). Bulk SiO₂ (bSiO₂) and sodium silicate (Na₂SiO₃) were used as controls. The results showed that nSiO₂ and nSiC increased the shoot length (11-37%, 6-25%) and root length (17-87%, 59-207%) of germinating seeds, respectively, compared with the control. Similarly, inter-root exposure to nSiO₂, bSiO_2, and nSiC improved the activity of aboveground catalase (10-55%, 31-34%, and 13-51%) and increased the content of trace elements magnesium, copper, and zinc, thus promoting the photosynthesis of rice. However, Na₂SiO₃ at a concentration of 200 mg/L reduced the aboveground and root biomass of rice by 27-51% and 4-17%, respectively. This may be because excess silicon not only inhibited the activity of root antioxidant enzymes but also disrupted the balance of mineral elements. This finding provides a new basis for the effect of silica-based NMs promotion on seed germination and rice growth.

A Ghd7-centered regulatory network provides a mechanistic approximation to optimal heterosis in an elite rice hybrid

Heterosis refers to the superior performance of hybrids over their parents, which is a general phenomenon occurring in diverse organisms. Many commercial hybrids produce high yield without delayed flowering, which we refer to as optimal heterosis and is desired in hybrid breeding. Here, we attempted to illustrate the genomic basis of optimal heterosis by reinvestigating the single-locus quantitative trait loci and digenic interactions of two traits, the number of spikelets per panicle (SP) and heading date (HD), using recombinant inbred lines and 'immortalized F₂ s' derived from the elite rice (Oryza sativa) hybrid Shanyou 63. Our analysis revealed a regulatory network that may provide an approximation to the genetic constitution of the optimal heterosis observed in this hybrid. In this network, Ghd7 works as the core element, and three other genes, Ghd7.1, Hd1, and Hd3a/RFT1, also have major roles. The effects of positive dominance by Ghd7 and Ghd7.1 and negative dominance by Hd1 and Hd3a/RFT1 in the hybrid background contribute the major part to the high SP without delaying HD; numerous epistatic interactions, most of which involve Ghd7, also play important roles collectively. The results expand our understanding of the genic interaction networks underlying hybrid rice breeding programs, which may be very useful in future crop genetic improvement.

OsbHLHq11, the Basic Helix-Loop-Helix Transcription Factor, Involved in Regulation of Chlorophyll Content in Rice

Photosynthesis is an important factor in determining the yield of rice. In particular, the size and efficiency of the photosynthetic system after the heading has a great impact on the yield. Research related to high-efficiency photosynthesis is essential to meet the growing demands of crops for the growing population. Chlorophyll is a key molecule in photosynthesis, a pigment that acts as an antenna to absorb light energy. Improvement of chlorophyll content characteristics has been emphasized in rice breeding for several decades. It is expected that an increase in chlorophyll content may increase photosynthetic efficiency, and understanding the genetic basis involved is important. In this study, we measured leaf color (CIELAB), chlorophyll content (SPAD), and chlorophyll fluorescence, and quantitative trait loci (QTL) mapping was performed using 120 Cheongcheong/Nagdong double haploid (CNDH) line after the heading date. A major QTL related to chlorophyll content was detected in the RM26981-RM287 region of chromosome 11. OsbHLHq11 was finally selected through screening of genes related to chlorophyll content in the RM26981-RM287 region. The relative expression level of the gene of OsbHLHq11 was highly expressed in cultivars with low chlorophyll content, and is expected to have a similar function to BHLH62 of the Gramineae genus. OsbHLHq11 is expected to increase photosynthetic efficiency by being involved in the chlorophyll content, and is expected to be utilized as a new genetic resource for breeding high-yield rice.

Partially functional NARROW LEAF1 balances leaf photosynthesis and plant architecture for greater rice yield

NARROW LEAF1 (NAL1) is an elite gene in rice (Oryza sativa), given its close connection to leaf photosynthesis, hybrid vigor, and yield-related agronomic traits; however, the underlying mechanism by which this gene affects these traits remains elusive. In this study, we systematically measured leaf photosynthetic parameters, leaf anatomical parameters, architectural parameters, and agronomic traits in indica cultivar 9311, in 9311 with the native NAL1 replaced by the Nipponbare NAL1 (9311-NIL), and in 9311 with the NAL1 fully mutated (9311-nal1). Leaf length, width, and spikelet number gradually increased from lowest to highest in 9311-nal1, 9311, and 9311-NIL. In contrast, the leaf photosynthetic rate on a leaf area basis, leaf thickness, and panicle number gradually decreased from highest to lowest in 9311-nal1, 9311, and 9311-NIL. RNA-seq analysis showed that NAL1 negatively regulates the expression of photosynthesis-related genes; NAL1 also influenced expression of many genes related to phytohormone signaling, as also shown by different leaf contents of 3-Indoleacetic acid, jasmonic acid, Gibberellin A3, and isopentenyladenine among these genotypes. Furthermore, field experiments with different planting densities showed that 9311 had a larger biomass and yield advantage under low planting density compared to either 9311-NIL or 9311-nall. This study shows both direct and indirect effects of NAL1 on leaf photosynthesis; furthermore, we show that a partially functional NAL1 allele helps maintain a balanced leaf photosynthesis and plant architecture for increased biomass and grain yield in the field.

A model-guided holistic review of exploiting natural variation of photosynthesis traits in crop improvement

Breeding for improved leaf photosynthesis is considered as a viable approach to increase crop yield. Whether it should be improved in combination with other traits has not been assessed critically. Based on the quantitative crop model GECROS that interconnects various traits to crop productivity, we review natural variation in relevant traits, from biochemical aspects of leaf photosynthesis to morpho-physiological crop characteristics. While large phenotypic variations (sometimes >2-fold) for leaf photosynthesis and its underlying biochemical parameters were reported, few quantitative trait loci (QTL) were identified, accounting for a small percentage of phenotypic variation. More QTL were reported for sink size (that feeds back on photosynthesis) or morpho-physiological traits (that affect canopy productivity and duration), together explaining a much greater percentage of their phenotypic variation. Traits for both photosynthetic rate and sustaining it during grain filling were strongly related to nitrogen-related traits. Much of the molecular basis of known photosynthesis QTL thus resides in genes controlling photosynthesis indirectly. Simulation using GECROS demonstrated the overwhelming importance of electron transport parameters, compared with the maximum Rubisco activity that largely determines the commonly studied light-saturated photosynthetic rate. Exploiting photosynthetic natural variation might significantly improve crop yield if nitrogen uptake, sink capacity, and other morpho-physiological traits are co-selected synergistically.

Genome-wide analysis of the CCT gene family in Chinese white pear (Pyrus bretschneideri Rehd.) and characterization of PbPRR2 in response to varying light signals

BACKGROUND

Canopy architecture is critical in determining the light environment and subsequently the photosynthetic productivity of fruit crops. Numerous CCT domain-containing genes are crucial for plant adaptive responses to diverse environmental cues. Two CCT genes, the orthologues of AtPRR5 in pear, have been reported to be strongly correlated with photosynthetic performance under distinct canopy microclimates. However, knowledge concerning the specific expression patterns and roles of pear CCT family genes (PbCCTs) remains very limited. The key roles played by PbCCTs in the light response led us to examine this large gene family in more detail.

RESULTS

Genome-wide sequence analysis identified 42 putative PbCCTs in the genome of pear (Pyrus bretschneideri Rehd.). Phylogenetic analysis indicated that these genes were divided into five subfamilies, namely, COL (14 members), PRR (8 members), ZIM (6 members), TCR1 (6 members) and ASML2 (8 members). Analysis of exon-intron structures and conserved domains provided support for the classification. Genome duplication analysis indicated that whole-genome duplication/segmental duplication events played a crucial role in the expansion of the CCT family in pear and that the CCT family evolved under the effect of purifying selection. Expression profiles exhibited diverse expression patterns of PbCCTs in various tissues and in response to varying light signals. Additionally, transient overexpression of PbPRR2 in tobacco leaves resulted in inhibition of photosynthetic performance, suggesting its possible involvement in the repression of photosynthesis.

CONCLUSIONS

This study provides a comprehensive analysis of the CCT gene family in pear and will facilitate further functional investigations of PbCCTs to uncover their biological roles in the light response.

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study’s novel candidate genes and allele combinations will provide a valuable reference for future breeding work.

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

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

Genomic prediction and QTL mapping of root system architecture and above-ground agronomic traits in rice (Oryza sativa L.) with a multitrait index and Bayesian networks

Root system architecture (RSA) is a crucial factor in resource acquisition and plant productivity. Roots are difficult to phenotype in the field, thus new tools for predicting phenotype from genotype are particularly valuable for plant breeders aiming to improve RSA. This study identifies quantitative trait loci (QTLs) for RSA and agronomic traits in a rice (Oryza sativa) recombinant inbred line (RIL) population derived from parents with contrasting RSA traits (PI312777 × Katy). The lines were phenotyped for agronomic traits in the field, and separately grown as seedlings on agar plates which were imaged to extract RSA trait measurements. QTLs were discovered from conventional linkage analysis and from a machine learning approach using a Bayesian network (BN) consisting of genome-wide SNP data and phenotypic data. The genomic prediction abilities (GPAs) of multi-QTL models and the BN analysis were compared with the several standard genomic prediction (GP) methods. We found GPAs were improved using multitrait (BN) compared to single trait GP in traits with low to moderate heritability. Two groups of individuals were selected based on GPs and a modified rank sum index (GSRI) indicating their divergence across multiple RSA traits. Selections made on GPs did result in differences between the group means for numerous RSA. The ranking accuracy across RSA traits among the individual selected RILs ranged from 0.14 for root volume to 0.59 for lateral root tips. We conclude that the multitrait GP model using BN can in some cases improve the GPA of RSA and agronomic traits, and the GSRI approach is useful to simultaneously select for a desired set of RSA traits in a segregating population.

Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: cross-species 2D GWAS reveals the underlying genetics

A 1D/2D genome-wide association study strategy was adopted to investigate the genetic systems underlying the reciprocal adaptation of rice (Oryza sativa) and its bacterial pathogen, Xanthomonas oryzae pv. oryzae (Xoo) using the whole-genome sequencing and large-scale phenotyping data of 701 rice accessions and 23 diverse Xoo strains. Forty-seven Xoo virulence-related genes and 318 rice quantitative resistance genes (QR-genes) mainly located in 41 genomic regions, and genome-wide interactions between the detected virulence-related genes and QR genes were identified, including well-known resistance genes/virulence genes plus many previously uncharacterized ones. The relationship between rice and Xoo was characterized by strong differentiation among Xoo races corresponding to the subspecific differentiation of rice, by strong shifts toward increased resistance/virulence of rice/Xoo populations and by rich genetic diversity at the detected rice QR-genes and Xoo virulence genes, and by genome-wide interactions between many rice QR-genes and Xoo virulence genes in a multiple-to-multiple manner, presumably resulting either from direct protein-protein interactions or from genetic epistasis. The observed complex genetic interaction system between rice and Xoo likely exists in other crop-pathogen systems that would maintain high levels of diversity at their QR-loci/virulence-loci, resulting in dynamic coevolutionary consequences during their reciprocal adaptation.

CO2 -responsive CCT protein interacts with 14-3-3 proteins and controls the expression of starch synthesis-related genes

CO₂ -responsive CCT protein (CRCT) is a positive regulator of starch synthesis-related genes such as ADP-glucose pyrophosphorylase large subunit 1 and starch branching enzyme I particularly in the leaf sheath of rice (Oryza sativa L.). The promoter GUS analysis revealed that CRCT expressed exclusively in the vascular bundle, whereas starch synthesis-related genes were expressed in different sites such as mesophyll cell and starch storage parenchyma cell. However, the chromatin immunoprecipitation (ChIP) using a FLAG-CRCT overexpression line and subsequent qPCR analyses showed that the 5'-flanking regions of these starch synthesis-related genes tended to be enriched by ChIP, suggesting that CRCT can bind to the promoter regions of these genes. The monomer of CRCT is 34.2 kDa; however, CRCT was detected at 270 kDa via gel filtration chromatography, suggesting that CRCT forms a complex in vivo. Immunoprecipitation and subsequent MS analysis pulled down several 14-3-3-like proteins. A yeast two-hybrid analysis and bimolecular fluorescence complementation assays confirmed the interaction between CRCT and 14-3-3-like proteins. Although there is an inconsistency in the place of expression, this study provides important findings regarding the molecular function of CRCT to control the expression of key starch synthesis-related genes.

Genome-wide association studies: assessing trait characteristics in model and crop plants

GWAS involves testing genetic variants across the genomes of many individuals of a population to identify genotype–phenotype association. It was initially developed and has proven highly successful in human disease genetics. In plants genome-wide association studies (GWAS) initially focused on single feature polymorphism and recombination and linkage disequilibrium but has now been embraced by a plethora of different disciplines with several thousand studies being published in model and crop species within the last decade or so. Here we will provide a comprehensive review of these studies providing cases studies on biotic resistance, abiotic tolerance, yield associated traits, and metabolic composition. We also detail current strategies of candidate gene validation as well as the functional study of haplotypes. Furthermore, we provide a critical evaluation of the GWAS strategy and its alternatives as well as future perspectives that are emerging with the emergence of pan-genomic datasets.

The Ghd7 transcription factor represses ARE1 expression to enhance nitrogen utilization and grain yield in rice

The genetic improvement of nitrogen use efficiency (NUE) of crops is vital for grain productivity and sustainable agriculture. However, the regulatory mechanism of NUE remains largely elusive. Here, we report that the rice Grain number, plant height, and heading date7 (Ghd7) gene genetically acts upstream of ABC1 REPRESSOR1 (ARE1), a negative regulator of NUE, to positively regulate nitrogen utilization. As a transcriptional repressor, Ghd7 directly binds to two Evening Element-like motifs in the promoter and intron 1 of ARE1, likely in a cooperative manner, to repress its expression. Ghd7 and ARE1 display diurnal expression patterns in an inverse oscillation manner, mirroring a regulatory scheme based on these two loci. Analysis of a panel of 2656 rice varieties suggests that the elite alleles of Ghd7 and ARE1 have undergone diversifying selection during breeding. Moreover, the allelic distribution of Ghd7 and ARE1 is associated with the soil nitrogen deposition rate in East Asia and South Asia. Remarkably, the combination of the Ghd7 and ARE1 elite alleles substantially improves NUE and yield performance under nitrogen-limiting conditions. Collectively, these results define a Ghd7-ARE1-based regulatory mechanism of nitrogen utilization, providing useful targets for genetic improvement of rice NUE.

Genome-Wide Association Study Identified Novel Candidate Loci/Genes Affecting Lodging Resistance in Rice

Lodging reduces rice yield, but increasing lodging resistance (LR) usually limits yield potential. Stem strength and leaf type are major traits related to LR and yield, respectively. Hence, understanding the genetic basis of stem strength and leaf type is of help to reduce lodging and increase yield in LR breeding. Here, we carried out an association analysis to identify quantitative trait locus (QTLs) affecting stem strength-related traits (internode length/IL, stem wall thickness/SWT, stem outer diameter/SOD, and stem inner diameter/SID) and leaf type-associated traits (Flag leaf length/FLL, Flag leaf angle/FLA, Flag leaf width/FLW, leaf-rolling/LFR and SPAD/Soil, and plant analyzer development) using a diverse panel of 550 accessions and evaluated over two years. Genome-wide association study (GWAS) using 4,076,837 high-quality single-nucleotide polymorphisms (SNPs) identified 89 QTLs for the nine traits. Next, through "gene-based association analysis, haplotype analysis, and functional annotation", the scope was narrowed down step by step. Finally, we identified 21 candidate genes in 9 important QTLs that included four reported genes (TUT1, OsCCC1, CFL1, and ACL-D), and seventeen novel candidate genes. Introgression of alleles, which are beneficial for both stem strength and leaf type, or pyramiding stem strength alleles and leaf type alleles, can be employed for LR breeding. All in all, the experimental data and the identified candidate genes in this study provide a useful reference for the genetic improvement of rice LR.

Global analysis of CCT family knockout mutants identifies four genes involved in regulating heading date in rice

Many genes encoding CCT domain-containing proteins regulate flowering time. In rice (Oryza sativa), 41 such genes have been identified, but only a few have been shown to regulate heading date. Here, to test whether and how additional CCT family genes regulate heading date in rice, we classified these genes into five groups based on their diurnal expression patterns. The expression patterns of genes in the same subfamily or in close phylogenetic clades tended to be similar. We generated knockout mutants of the entire gene family via CRISPR/Cas9. The heading dates of knockout mutants of only 4 of 14 genes previously shown to regulate heading date were altered, pointing to functional redundancy of CCT family genes in regulating this trait. Analysis of mutants of four other genes showed that OsCCT22, OsCCT38, and OsCCT41 suppress heading under long-day conditions and promote heading under short-day conditions. OsCCT03 promotes heading under both conditions and upregulates the expression of Hd1 and Ehd1, a phenomenon not previously reported for other such genes. To date, at least 18 CCT domain-containing genes involved in regulating heading have been identified, providing diverse, flexible gene combinations for generating rice varieties with a given heading date.

Current Understanding of Leaf Senescence in Rice

Leaf senescence, which is the last developmental phase of plant growth, is controlled by multiple genetic and environmental factors. Leaf yellowing is a visual indicator of senescence due to the loss of the green pigment chlorophyll. During senescence, the methodical disassembly of macromolecules occurs, facilitating nutrient recycling and translocation from the sink to the source organs, which is critical for plant fitness and productivity. Leaf senescence is a complex and tightly regulated process, with coordinated actions of multiple pathways, responding to a sophisticated integration of leaf age and various environmental signals. Many studies have been carried out to understand the leaf senescence-associated molecular mechanisms including the chlorophyll breakdown, phytohormonal and transcriptional regulation, interaction with environmental signals, and associated metabolic changes. The metabolic reprogramming and nutrient recycling occurring during leaf senescence highlight the fundamental role of this developmental stage for the nutrient economy at the whole plant level. The strong impact of the senescence-associated nutrient remobilization on cereal productivity and grain quality is of interest in many breeding programs. This review summarizes our current knowledge in rice on (i) the actors of chlorophyll degradation, (ii) the identification of stay-green genotypes, (iii) the identification of transcription factors involved in the regulation of leaf senescence, (iv) the roles of leaf-senescence-associated nitrogen enzymes on plant performance, and (v) stress-induced senescence. Compiling the different advances obtained on rice leaf senescence will provide a framework for future rice breeding strategies to improve grain yield.

Genetic analyses of lodging resistance and yield provide insights into post-Green-Revolution breeding in rice

Lodging reduces grain yield in cereal crops. Understanding the genetic basis of lodging resistance (LR) benefits LR breeding. In the study, 524 accessions from a rice germplasm collection and 193 recombinant inbred lines were phenotyped for 17 LR-related traits. Height and culm strength (the magnitude of applied force necessary to break the culm) were two major factors affecting LR. We conducted genome-wide association study (GWAS) and identified 127 LR-associated loci. Significant phenotypic correlations between culm-strength traits and yield-related traits were observed. To reveal the genetic relationship between them, we conducted GWAS of culm-strength traits with adding yield-related trait as a covariate and detected 63 loci linking culm strength and yield. As a proof, a near-isogenic line for an association locus on chromosome 7 showed enhanced LR and yield. Strikingly, 58 additional loci were identified in the covariate-added GWAS. Several LR-associated loci had undergone divergent selection. Linkage analysis supported the GWAS results. We propose that introgression of alleles beneficial for both culm strength and panicle weight without negative effects on panicle number or pyramiding high-yielding alleles and lodging-resistant alleles without effects on yield can be employed for the post-Green-Revolution breeding.

Spatiotemporal Heterogeneity of Chlorophyll Content and Fluorescence Response Within Rice (Oryza sativa L.) Canopies Under Different Nitrogen Treatments

Accurate acquisition of plant phenotypic information has raised long-standing concerns in support of crop breeding programs. Different methods have been developed for high throughput plant phenotyping, while they mainly focused on the canopy level without considering the spatiotemporal heterogeneity at different canopy layers and growth stages. This study aims to phenotype spatiotemporal heterogeneity of chlorophyll (Chl) content and fluorescence response within rice leaves and canopies. Multipoint Chl content and high time-resolved Chl a fluorescence (ChlF) transient (OJIP transient) of rice plants were measured at different nitrogen levels and growth stages. Results showed that the Chl content within the upper leaves exhibited an increasing trend from the basal to the top portions but a decreasing pattern within the lower leaves at the most growth stages. Leaf Chl content within the rice canopy was higher in the lower leaves in the vegetative phase, while from the initial heading stage the pattern gradually reversed with the highest Chl content appearing in the upper leaves. Nitrogen supply mainly affects the occurrence time of the reverse vertical pattern. This could be the result of different nutritional demands of leaves transforming from sinks to sources, and it was further confirmed by the fall of the JI phase of OJIP transient in the vegetative phase and the rise in the reproductive phase. We further deduced that the vertical distribution of Chl content could have a defined pattern at a specific growth stage. Furthermore, the reduction of end acceptors at photosystem I (PSI) electron acceptor side per cross section (RE0/CS) was found to be a potential sensitive predictor for identifying the vertical heterogeneity of leaf Chl content. These findings provide prior knowledge on the vertical profiles of crop physiological traits, which explore the opportunity to develop more efficient plant phenotyping tools for crop breeding.

Genome-wide characterization and analysis of the CCT motif family genes in soybean (Glycine max)

MAIN CONCLUSION

Soybean possesses 19 CMF genes which mainly arose from duplication events. Their features and motifs are highly conserved but transcriptional data indicated functional diversity in metabolism and stress responses. CCT [for CONSTANS, CONSTANS-like (CO-like), and timing of CAB expression1 (TOC1)] domain-containing genes play important roles in regulating flowering, plant growth, and grain yield and are also involved in stress responses. The CMF (CCT motif family) genes, included in the CCT family, contain a single CCT domain as the only identifiable domain in their predicted protein sequence and are interesting targets for breeding programs. In this study, we identified 19 putative GmCMF genes, based on the latest soybean (Glycine max) genome annotation. The predicted GmCMF proteins were characterized based on conserved structural features, and a phylogenetic tree was constructed including all CMF proteins from rice and Arabidopsis as representative examples of the monocotyledonous (monocot) and dicotyledonous (dicot) plants, respectively. High similarities in the conserved motifs of the protein sequences and the gene structures were found. In addition, by analyzing the CMF gene family in soybean, we identified seven pairs of genes that originated from segmental chromosomal duplication events attributable to the most recent whole-genome duplication (WGD) event in the Glycine lineage. Expression analysis of GmCMF genes in various tissues and after specific treatments demonstrated tissue and stress-response specific differential expression. Gene expression analysis was complemented by the identification of putative cis-elements present in the promoter regions of the genes through a bioinformatics approach, using the existing soybean reference genome sequence and gene models. Co-functional networks inferred from distinct types of genomics data-including microarrays and RNA-seq samples from soybean-revealed that GmCMF genes might play crucial roles in metabolism and transport processes. The results of this study, the first systematic analysis of the soybean CCT gene family, can serve as a strong foundation for further elucidation of their physiological functions and biological roles.

Genetic basis of vascular bundle variations in rice revealed by genome-wide association study

The vascular bundles play important roles in transportation of photoassimilate, and the number, size, and capacity of vascular bundles influence the transportation efficiency. Dissecting the genetic basis may help to make better use of naturally occurring vascular bundle variations. Here, we conducted a genome-wide association study (GWAS) of the vascular bundle variations in a worldwide collection of 529 Oryza sativa accessions. A total of 42 and 93 significant association loci were identified in the neck panicle and flag leaf, respectively. The introgression lines showing extreme values of the target traits harbored at least one GWAS signal, indicating the reliability of the GWAS loci. Based on the data of near-isogenic lines and transgenic plants, Grain number, plant height, and heading date7 (Ghd7) was identified as a major locus for the natural variation of vascular bundles in the neck panicle at the heading stage. In addition, Narrow leaf1 (NAL1) was found to influence the vascular bundles in both the neck panicle and flag leaf, and the effects of the major haplotypes of NAL1 were characterized. The loci or candidate genes identified would help to improve vascular bundle system in rice breeding.

Association Mapping and Transcriptome Analysis Reveal the Genetic Architecture of Maize Kernel Size

Kernel length, kernel width, and kernel thickness are important traits affecting grain yield and product quality. Here, the genetic architecture of the three kernel size traits was dissected in an association panel of 309 maize inbred lines using four statistical methods. Forty-two significant single nucleotide polymorphisms (SNPs; p < 1.72E-05) and 70 genes for the three traits were identified under five environments. One and eight SNPs were co-detected in two environments and by at least two methods, respectively, and they explained 5.87-9.59% of the phenotypic variation. Comparing the transcriptomes of two inbred lines with contrasting seed size, three and eight genes identified in the association panel showed significantly differential expression between the two genotypes at 15 and 39 days after pollination, respectively. Ten and 17 genes identified by a genome-wide association study were significantly differentially expressed between the two development stages in the two genotypes. Combining environment-/method-stable SNPs and differential expression analysis, ribosomal protein L7, jasmonate-regulated gene 21, serine/threonine-protein kinase RUNKEL, AP2-EREBP-transcription factor 16, and Zm00001d035222 (cell wall protein IFF6-like) were important candidate genes for maize kernel size and development.

Genome-wide association study of flowering time reveals complex genetic heterogeneity and epistatic interactions in rice

Since domestication, rice has cultivated in a wide range of latitudes with different day lengths. Selection of diverse natural variations in heading date and photoperiod sensitivity is critical for adaptation of rice to different geographical environments. To unravel the genetic architecture underlying natural variation of rice flowering time, we conducted a genome wide association study (GWAS) using several association analysis strategies with a diverse worldwide collection of 529 O. sativa accessions. Heading date was investigated in three environments under long-day or short-day conditions, and photosensitivity was evaluated. By dividing the whole association panel into subpopulations and performing GWAS with both linear mixed models and multi-locus mixed-models, we revealed hundreds of significant loci harboring novel candidate genes as well as most of the known flowering time genes. In total, 127 hotspots were detected in at least two GWAS. Universal genetic heterogeneity was found across subpopulations. We further detected abundant interactions between GWAS loci, especially in indica. Functional gene families were revealed from enrichment analysis of the 127 hotspots. The results demonstrated a rich of genetic interactions in rice flowering time genes and such epistatic interactions contributed to the large portions of missing heritability in GWAS. It suggests the increased complexity of genetic heterogeneity might discount the power of increasing the sample sizes in GWAS.

Alterations in stomatal response to fluctuating light increase biomass and yield of rice under drought conditions

The acceleration of stomatal closure upon high to low light transition could improve plant water use efficiency and drought tolerance. Herein, using genome-wide association study, we showed that the genetic variation in OsNHX1 was strongly associated with the changes in τ_cl , the time constant of stomatal closure, in 206 rice accessions. OsNHX1 overexpression in rice resulted in a decrease in τ_cl , and an increase in biomass, grain yield under drought. Conversely, OsNHX1 knockout by CRISPR/CAS9 shows opposite trends for these traits. We further found three haplotypes spanning the OsNHX1 promoter and CDS regions. Two among them, HapII and HapIII, were found to be associated with a high and low τ_cl , respectively. A near-isogenic line (NIL, S464) was developed through replacing the genomic region harboring HapII (~10 kb) from MH63 (recipient) rice cultivar by the same sized genomic region containing Hap III from 02428 (donor). Compared with MH63, S464 shows a reduction by 35% in τ_cl and an increase by 40% in the grain yield under drought. However, under normal conditions, S464 maintains closely similar grain yield as MH63. The global distribution of the two OsNHX1 haplotypes is associated with the local precipitation. Taken together, the natural variation in OsNHX1 could be utilized to manipulate the stomatal dynamics for an improved rice drought tolerance.

Whole-genome de novo assemblies reveal extensive structural variations and dynamic organelle-to-nucleus DNA transfers in African and Asian rice

Asian cultivated rice (Oryza sativa) and African cultivated rice (Oryza glaberrima) originated from the wild rice species Oryza rufipogon and Oryza barthii, respectively. The genomes of both cultivated species have undergone profound changes during domestication. Whole-genome de novo assemblies of O. barthii, O. glaberrima, O. rufipogon and Oryza nivara, produced using PacBio single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies, showed that Gypsy-like retrotransposons are the major contributors to genome size variation in African and Asian rice. Through the detection of genome-wide structural variations (SVs), we observed that besides 28 shared SV hot spots, another 67 hot spots existed in either the Asian or African rice genomes. Based on gene annotation information of the SVs, we established that organelle-to-nucleus DNA transfers resulted in numerous SVs that participated in the nuclear genome divergence of rice species and subspecies. We detected 52 giant nuclear integrants of organelle DNA (NORGs, defined as >10 kb) in six Oryza AA genomes. In addition, we developed an effective method to genotype giant NORGs, based on genome assembly, and first showed the dynamic change in the distribution of giant NORGs in rice natural population. Interestingly, 16 highly differentiated giant NORGs tended to accumulate in natural populations of Asian rice from higher latitude regions, grown at lower temperatures and light intensities. Our study provides new insight into the genome divergence of African and Asian rice, and establishes that organelle-to-nucleus DNA transfers, as potentially powerful contributors to environmental adaptation during rice evolution, play a major role in producing SVs in rice genomes.

CCT domain-containing genes in cereal crops: flowering time and beyond

The review summarizes the functions of the plant special transcription factors CCT family genes in multiple traits and discusses the molecular breeding strategies with CCT family genes in the future. Plants integrate circadian clock and external signals such as temperature and photoperiod to synchronize flowering with seasonal environmental changes. This process makes cereal crops including short-day crops, such as rice and maize, and long-day crops, such as wheat and barley, better adapt to varied growth zones from temperate to tropical regions. CCT family genes involve circadian clock and photoperiodic flowering pathways and help plants set a suitable flowering time to produce offspring. Beyond the flowering time, CCT family genes in cereal crops are associated with biomass and grain yield. Moreover, recent studies showed that they also associate with photosynthesis, nutrition use efficiency and stress tolerance. Here, we systematically review the progress in functional characterization of CCT family genes in flowering, geographical adaptation and grain yield formation, raise the core questions related to their molecular mechanisms and discuss how to practice them in genetic improvement in cereal crops by combining gene diagnosis and top-level design.

Advances in genome-wide association studies of complex traits in rice

Genome-wide association studies (GWAS), genetic surveys of the whole genome to detect variants associated with a trait in natural populations, are a powerful approach for dissecting complex traits. This genetic mapping approach has been applied in rice over the last 10 years. During the last decade, GWAS was used to identify the loci underlying tens of rice traits, and several important genes were detected in GWAS and further confirmed in follow-up functional experiments. In this review, we present an overview of the whole process in a typical GWAS, including population design, genotyping, phenotyping and analysis methods. Recent advances in rice GWAS are also provided, including several examples of the functional characterization of candidate genes. The possible breakthroughs of rice GWAS in the next decade are discussed with regard to their application in breeding, the consideration of epistatic interactions and in-depth functional annotations of DNA elements and genetic variants throughout the rice genome.

Analysis of genetic architecture and favorable allele usage of agronomic traits in a large collection of Chinese rice accessions

Genotyping and phenotyping large natural populations provide opportunities for population genomic analysis and genome-wide association studies (GWAS). Several rice populations have been re-sequenced in the past decade; however, many major Chinese rice cultivars were not included in these studies. Here, we report large-scale genomic and phenotypic datasets for a collection mainly comprised of 1,275 rice accessions of widely planted cultivars and parental hybrid rice lines from China. The population was divided into three indica/Xian and three japonica/Geng phylogenetic subgroups that correlate strongly with their geographic or breeding origins. We acquired a total of 146 phenotypic datasets for 29 agronomic traits under multi-environments for different subpopulations. With GWAS, we identified a total of 143 significant association loci, including three newly identified candidate genes or alleles that control heading date or amylose content. Our genotypic analysis of agronomically important genes in the population revealed that many favorable alleles are underused in elite accessions, suggesting they may be used to provide improvements in future breeding efforts. Our study provides useful resources for rice genetics research and breeding.

A xylan glucuronosyltransferase gene exhibits pleiotropic effects on cellular composition and leaf development in rice

Leaf chlorophyll content is an important physiological indicator of plant growth, metabolism and nutritional status, and it is highly correlated with leaf nitrogen content and photosynthesis. In this study, we report the cloning and identification of a xylan glucuronosyltransferase gene (OsGUX1) that affects relative chlorophyll content in rice leaf. Using a set of chromosomal segment substitution lines derived from a cross of wild rice accession ACC10 and indica variety Zhenshan 97 (ZS97), we identified numerous quantitative trait loci for relative chlorophyll content. One major locus of them for relative chlorophyll content was mapped to a 10.3-kb region that contains OsGUX1. The allele OsGUX1^AC from ACC10 significantly decreases nitrogen content and chlorophyll content of leaf compared with OsGUX1^ZS from ZS97. The overexpression of OsGUX1 reduced chlorophyll content, and the suppression of this gene increased chlorophyll content of rice leaf. OsGUX1 is located in Golgi apparatus, and highly expressed in seedling leaf and the tissues in which primary cell wall synthesis occurring. Our experimental data indicate that OsGUX1 is responsible for addition of glucuronic acid residues onto xylan and participates in accumulation of cellulose and hemicellulose in the cell wall deposition, thus thickening the primary cell wall of mesophyll cells, which might lead to reduced chlorophyll content in rice leaf. These findings provide insights into the association of cell wall components with leaf nitrogen content in rice.

Dissecting the Genetic Complexity of Fusarium Crown Rot Resistance in Wheat

Fusarium crown rot (FCR) is one of the most important diseases of wheat (Triticum aestivum L.). FCR is mainly caused by the fungal pathogens Fusarium culmorum and F. pseudograminearum. In order to identify new sources of resistance to FCR and to dissect the complexity of FCR resistance, a panel of 161 wheat accessions was phenotyped under growth room (GR) and greenhouse conditions (GH). Analysis of variance showed significant differences in crown rot development among wheat accessions and high heritability of genotype-environment interactions for GR (0.96) and GH (0.91). Mixed linear model analysis revealed seven novel quantitative trait loci (QTLs) linked to F. culmorum on chromosomes 2AL, 3AS, 4BS, 5BS, 5DS, 5DL and 6DS for GR and eight QTLs on chromosomes on 3AS, 3BS, 3DL, 4BS (2), 5BS, 6BS and 6BL for GH. Total phenotypic variances (R²) explained by the QTLs linked to GR and GH were 48% and 59%, respectively. In addition, five favorable epistasis interactions among the QTLs were detected for both GR and GH with and without main effects. Epistatic interaction contributed additional variation up to 21% under GR and 7% under GH indicating strong effects of environment on the expression of QTLs. Our results revealed FCR resistance responses in wheat to be complex and controlled by multiple QTLs.

Crop Phenomics and High-Throughput Phenotyping: Past Decades, Current Challenges, and Future Perspectives

Since whole-genome sequencing of many crops has been achieved, crop functional genomics studies have stepped into the big-data and high-throughput era. However, acquisition of large-scale phenotypic data has become one of the major bottlenecks hindering crop breeding and functional genomics studies. Nevertheless, recent technological advances provide us potential solutions to relieve this bottleneck and to explore advanced methods for large-scale phenotyping data acquisition and processing in the coming years. In this article, we review the major progress on high-throughput phenotyping in controlled environments and field conditions as well as its use for post-harvest yield and quality assessment in the past decades. We then discuss the latest multi-omics research combining high-throughput phenotyping with genetic studies. Finally, we propose some conceptual challenges and provide our perspectives on how to bridge the phenotype-genotype gap. It is no doubt that accurate high-throughput phenotyping will accelerate plant genetic improvements and promote the next green revolution in crop breeding.

Genetic dissection of yield-related traits via genome-wide association analysis across multiple environments in wild soybean (Glycine soja Sieb. and Zucc.)

MAIN CONCLUSION

A total of 41 SNPs were identified as significantly associated with five yield-related traits in wild soybean populations across multiple environments, and the candidate gene GsCID1 was found to be associated with seed weight. These results may facilitate improvements in cultivated soybean. Crop-related wild species contain new sources of genetic diversity for crop improvement. Wild soybean (Glycine soja Sieb. and Zucc.) is the progenitor of cultivated soybean [Glycine max (L.) Merr.] and can be used as an essential genetic resource for yield improvements. In this research, using genome-wide association study (GWAS) in 96 out of 113 wild soybean accessions with 114,090 single nucleotide polymorphisms (SNPs) (with minor allele frequencies ≤ 0.05), SNPs associated with five yield-related traits were identified across multiple environments. In total, 41 SNPs were significantly associated with the traits in two or more environments (significance threshold P ≤ 8.76 × 10-6), with 29, 7, 3, and 2 SNPs detected for 100-seed weight (SW), maturity time (MT), seed yield per plant (SY) and flowering time (FT), respectively. BLAST search against the Glycine soja W05 reference genome was performed, 20 candidate genes were identified based on these 41 significant SNPs. One candidate gene, GsCID1 (Glysoja.04g010563), harbored two significant SNPs-AX-93713187, with a non-synonymous mutation, and AX-93713188, with a synonymous mutation. GsCID1 was highly expressed during seed development based on public information resources. The polymorphisms in this gene were associated with SW. We developed a derived cleaved amplified polymorphic sequence (dCAPS) marker for GsCID1 that was highly associated with SW and was validated as a functional marker. In summary, the revealed SNPs/genes are useful for understanding the genetic architecture of yield-related traits in wild soybean, which could be used as a potential exotic resource to improve cultivated soybean yields.