Setaria viridis: a model for C4 photosynthesis

Utilizing machine learning and bioinformatics analysis to identify drought-responsive genes affecting yield in foxtail millet

Drought stress is a major constraint on crop development, potentially causing huge yield losses and threatening global food security. Improving Crop's stress tolerance is usually associated with a yield penalty. One way to balance yield and stress tolerance is modification specific gene by emerging precision genome editing technology. However, our knowledge of yield-related drought-tolerant genes is still limited. Foxtail millet (Setaria italica) has a remarkable tolerance to drought and is considered to be a model C4 crop that is easy to engineer. Here, we have identified 46 drought-responsive candidate genes by performing a machine learning-based transcriptome study on two drought-tolerant and two drought-sensitive foxtail millet cultivars. A total of 12 important drought-responsive genes were screened out by principal component analysis and confirmed experimentally by qPCR. Significantly, by investigating the haplotype of these genes based on 1844 germplasm resources, we found two genes (Seita.5G251300 and Seita.8G036300) exhibiting drought-tolerant haplotypes that possess an apparent advantage in 1000 grain weight and main panicle grain weight without penalty in grain weight per plant. These results demonstrate the potential of Seita.5G251300 and Seita.8G036300 for breeding drought-tolerant high-yielding foxtail millet. It provides important insights for the breeding of drought-tolerant high-yielding crop cultivars through genetic manipulation technology.

Plasma membrane aquaporins regulate root hydraulic conductivity in the model plant Setaria viridis

The high rate of productivity observed in panicoid crops is in part due to their extensive root system. Recently, green foxtail (Setaria viridis) has emerged as a genetic model system for panicoid grasses. Natural accessions of S. viridis originating from different parts of the world, with differential leaf physiological behavior, have been identified. This work focused on understanding the physiological and molecular mechanisms controlling root hydraulic conductivity and root-to-shoot gas exchange signaling in S. viridis. We identified 2 accessions, SHA and ZHA, with contrasting behavior at the leaf, root, and whole-plant levels. Our results indicated a role for root aquaporin (AQP) plasma membrane (PM) intrinsic proteins in the differential behavior of SHA and ZHA. Moreover, a different root hydraulic response to low levels of abscisic acid between SHA and ZHA was observed, which was associated with root AQPs. Using cell imaging, biochemical, and reverse genetic approaches, we identified PM intrinsic protein 1;6 (PIP1;6) as a possible PIP1 candidate that regulates radial root hydraulics and root-to-shoot signaling of gas exchange in S. viridis. In heterologous systems, PIP1;6 localized in the endoplasmic reticulum, and upon interaction with PIP2s, relocalization to the PM was observed. PIP1;6 was predominantly expressed at the root endodermis. Generation of knockout PIP1;6 plants (KO-PIP1;6) in S. viridis showed altered root hydraulic conductivity, altered gas exchange, and alteration of root transcriptional patterns. Our results indicate that PIPs are essential in regulating whole-plant water homeostasis in S. viridis. We conclude that root hydraulic conductivity and gas exchange are positively associated and are regulated by AQPs.

Genome-Wide Characterization and Haplotypic Variation Analysis of the YUC Gene Family in Foxtail Millet (Setaria italica)

Panicle development and grain production in crop species are essential breeding characteristics affected by the synthesis of auxin, which is influenced by flavin monooxygenase-encoding genes such as YUC (YUCCA) family members. In this trial, fourteen YUCs were identified and named uniformly in foxtail millet, an ancient crop species cultivated across the world. The phylogenetic analysis revealed that the SiYUCs were clustered into four subgroups; protein motif and gene structure analyses suggested that the closely clustered SiYUC genes were relatively conserved within each subgroup; while genome mapping analysis indicated that the SiYUC genes were unevenly distributed on foxtail millet chromosomes and colinear with other grass species. Transcription analysis revealed that the SiYUC genes differed greatly in expression pattern in different tissues and contained hormonal/light/stress-responding cis-elements. The haplotype characterization of SiYUC genes indicated many superior haplotypes of SiYUCs correlated with higher panicle and grain weight could be favorably selected by breeding. These results will be useful for the further study of the functional characteristics of SiYUC genes, particularly with regard to the marker-assisted pyramiding of beneficial haplotypes in foxtail millet breeding programs.

The osmotic stress-activated receptor-like kinase DPY1 mediates SnRK2 kinase activation and drought tolerance in Setaria

Plant genomes encode many receptor-like kinases (RLKs) that localize to the cell surface and perceive a wide variety of environmental cues to initiate downstream signaling cascades. Whether these RLKs participate in dehydration stress signaling in plants is largely unknown. DROOPY LEAF1 (DPY1), a leucine-rich repeat (LRR)-RLK, was recently shown to regulate plant architecture by orchestrating early brassinosteroid signaling in foxtail millet (Setaria italica). Here, we show that DPY1 is essential for the acclimation of foxtail millet to drought stress. DPY1 can be phosphorylated and activated in response to osmotic stress and is required for more than half of osmotic stress-induced global phosphorylation events, including the phosphorylation of sucrose nonfermenting kinase 2s (SnRK2s), the central kinases involved in osmotic stress. DPY1 acts upstream of STRESS-ACTIVATED PROTEIN KINASE 6 (SAPK6, a subclass I SnRK2) and is required for full SAPK6 activation, thereby allowing regulation of downstream genes to mount a response against drought stress. These signaling events are largely independent of DPY1-mediated brassinosteroid signaling. The DPY1-SAPK6 module is specific to seed plants and is absent in ancestral nonseed plants. Our findings reveal a dehydration stress-activated RLK that plays an indispensable role in osmotic stress signaling and mediates SnRK2 activation at the cell surface.

Genome-wide identification and expression profiling analysis of DIR gene family in Setaria italica

Dirigent (DIR) proteins play essential roles in regulating plant growth and development, as well as enhancing resistance to abiotic and biotic stresses. However, the whole-genome identification and expression profiling analysis of DIR gene family in millet (Setaria italica (Si)) have not been systematically understood. In this study, we conducted genome-wide identification and expression analysis of the S. italica DIR gene family, including gene structures, conserved domains, evolutionary relationship, chromosomal locations, cis-elements, duplication events, gene collinearity and expression patterns. A total of 38 SiDIR members distributed on nine chromosomes were screened and identified. SiDIR family members in the same group showed higher sequence similarity. The phylogenetic tree divided the SiDIR proteins into six subfamilies: DIR-a, DIR-b/d, DIR-c, DIR-e, DIR-f, and DIR-g. According to the tertiary structure prediction, DIR proteins (like SiDIR7/8/9) themselves may form a trimer to exert function. The result of the syntenic analysis showed that tandem duplication may play the major driving force during the evolution of SiDIRs. RNA-seq data displayed higher expression of 16 SiDIR genes in root tissues, and this implied their potential functions during root development. The results of quantitative real-time PCR (RT-qPCR) assays revealed that SiDIR genes could respond to the stress of CaCl2, CdCl, NaCl, and PEG6000. This research shed light on the functions of SiDIRs in responding to abiotic stress and demonstrated their modulational potential during root development. In addition, the membrane localization of SiDIR7/19/22 was confirmed to be consistent with the forecast. The results above will provide a foundation for further and deeper investigation of DIRs.

Identification of Two Flip-Over Genes in Grass Family as Potential Signature of C4 Photosynthesis Evolution

In flowering plants, C4 photosynthesis is superior to C3 type in carbon fixation efficiency and adaptation to extreme environmental conditions, but the mechanisms behind the assembly of C4 machinery remain elusive. This study attempts to dissect the evolutionary divergence from C3 to C4 photosynthesis in five photosynthetic model plants from the grass family, using a combined comparative transcriptomics and deep learning technology. By examining and comparing gene expression levels in bundle sheath and mesophyll cells of five model plants, we identified 16 differentially expressed signature genes showing cell-specific expression patterns in C3 and C4 plants. Among them, two showed distinctively opposite cell-specific expression patterns in C3 vs. C4 plants (named as FOGs). The in silico physicochemical analysis of the two FOGs illustrated that C3 homologous proteins of LHCA6 had low and stable pI values of ~6, while the pI values of LHCA6 homologs increased drastically in C4 plants Setaria viridis (7), Zea mays (8), and Sorghum bicolor (over 9), suggesting this protein may have different functions in C3 and C4 plants. Interestingly, based on pairwise protein sequence/structure similarities between each homologous FOG protein, one FOG PGRL1A showed local inconsistency between sequence similarity and structure similarity. To find more examples of the evolutionary characteristics of FOG proteins, we investigated the protein sequence/structure similarities of other FOGs (transcription factors) and found that FOG proteins have diversified incompatibility between sequence and structure similarities during grass family evolution. This raised an interesting question as to whether the sequence similarity is related to structure similarity during C4 photosynthesis evolution.

SHORT ROOT and INDETERMINATE DOMAIN family members govern PIN-FORMED expression to regulate minor vein differentiation in rice

C3 and C4 grasses directly and indirectly provide the vast majority of calories to the human diet, yet our understanding of the molecular mechanisms driving photosynthetic productivity in grasses is largely unexplored. Ground meristem cells divide to form mesophyll or vascular initial cells early in leaf development in C3 and C4 grasses. Here we define a genetic circuit composed of SHORT ROOT (SHR), INDETERMINATE DOMAIN (IDD), and PIN-FORMED (PIN) family members that specifies vascular identify and ground cell proliferation in leaves of both C3 and C4 grasses. Ectopic expression and loss-of-function mutant studies of SHR paralogs in the C3 plant Oryza sativa (rice) and the C4 plant Setaria viridis (green millet) revealed the roles of these genes in both minor vein formation and ground cell differentiation. Genetic and in vitro studies further suggested that SHR regulates this process through its interactions with IDD12 and 13. We also revealed direct interactions of these IDD proteins with a putative regulatory element within the auxin transporter gene PIN5c. Collectively, these findings indicate that a SHR-IDD regulatory circuit mediates auxin transport by negatively regulating PIN expression to modulate minor vein patterning in the grasses.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

An E2-E3 pair contributes to seed size control in grain crops

Understanding the molecular mechanisms that regulate grain yield is important for improving agricultural productivity. Protein ubiquitination controls various aspects of plant growth but lacks understanding on how E2-E3 enzyme pairs impact grain yield in major crops. Here, we identified a RING-type E3 ligase SGD1 and its E2 partner SiUBC32 responsible for grain yield control in Setaria italica. The conserved role of SGD1 was observed in wheat, maize, and rice. Furthermore, SGD1 ubiquitinates the brassinosteroid receptor BRI1, stabilizing it and promoting plant growth. Overexpression of an elite SGD1 haplotype improved grain yield by about 12.8% per plant, and promote complex biological processes such as protein processing in endoplasmic reticulum, stress responses, photosystem stabilization, and nitrogen metabolism. Our research not only identifies the SiUBC32-SGD1-BRI1 genetic module that contributes to grain yield improvement but also provides a strategy for exploring key genes controlling important traits in Poaceae crops using the Setaria model system.

The role of SWEET4 proteins in the post-phloem sugar transport pathway of Setaria viridis sink tissues

In the developing seeds of all higher plants, filial cells are symplastically isolated from the maternal tissue supplying photosynthate to the reproductive structure. Photoassimilates must be transported apoplastically, crossing several membrane barriers, a process facilitated by sugar transporters. Sugars Will Eventually be Exported Transporters (SWEETs) have been proposed to play a crucial role in apoplastic sugar transport during phloem unloading and the post-phloem pathway in sink tissues. Evidence for this is presented here for developing seeds of the C4 model grass Setaria viridis. Using immunolocalization, SvSWEET4 was detected in various maternal and filial tissues within the seed along the sugar transport pathway, in the vascular parenchyma of the pedicel, and in the xylem parenchyma of the stem. Expression of SvSWEET4a in Xenopus laevis oocytes indicated that it functions as a high-capacity glucose and sucrose transporter. Carbohydrate and transcriptional profiling of Setaria seed heads showed that there were some developmental shifts in hexose and sucrose content and consistent expression of SvSWEET4 homologues. Collectively, these results provide evidence for the involvement of SWEETs in the apoplastic transport pathway of sink tissues and allow a pathway for post-phloem sugar transport into the seed to be proposed.

MDSi: Multi-omics Database for Setaria italica

BACKGROUND

Foxtail millet (Setaria italica) harbors the small diploid genome (~ 450 Mb) and shows the high inbreeding rate and close relationship to several major foods, feed, fuel and bioenergy grasses. Previously, we created a mini foxtail millet, xiaomi, with an Arabidopsis-like life cycle. The de novo assembled genome data with high-quality and an efficient Agrobacterium-mediated genetic transformation system made xiaomi an ideal C₄ model system. The mini foxtail millet has been widely shared in the research community and as a result there is a growing need for a user-friendly portal and intuitive interface to perform exploratory analysis of the data.

RESULTS

Here, we built a Multi-omics Database for Setaria italica (MDSi, http://sky.sxau.edu.cn/MDSi.htm ), that contains xiaomi genome of 161,844 annotations, 34,436 protein-coding genes and their expression information in 29 different tissues of xiaomi (6) and JG21 (23) samples that can be showed as an Electronic Fluorescent Pictograph (xEFP) in-situ. Moreover, the whole-genome resequencing (WGS) data of 398 germplasms, including 360 foxtail millets and 38 green foxtails and the corresponding metabolic data were available in MDSi. The SNPs and Indels of these germplasms were called in advance and can be searched and compared in an interactive manner. Common tools including BLAST, GBrowse, JBrowse, map viewer, and data downloads were implemented in MDSi.

CONCLUSION

The MDSi constructed in this study integrated and visualized data from three levels of genomics, transcriptomics and metabolomics, and also provides information on the variation of hundreds of germplasm resources that can satisfies the mainstream requirements and supports the corresponding research community.

Increased sedoheptulose-1,7-bisphosphatase content in Setaria viridis does not affect C4 photosynthesis

Sedoheptulose-1,7-bisphosphatase (SBPase) is one of the rate-limiting enzymes of the Calvin cycle, and increasing the abundance of SBPase in C3 plants provides higher photosynthetic rates and stimulates biomass and yield. C4 plants usually have higher photosynthetic rates because they operate a biochemical CO2-concentrating mechanism between mesophyll and bundle sheath cells. In the C4 system, SBPase and other enzymes of the Calvin cycle are localized to the bundle sheath cells. Here we tested what effect increasing abundance of SBPase would have on C4 photosynthesis. Using green foxtail millet (Setaria viridis), a model C4 plant of NADP-ME subtype, we created transgenic plants with 1.5 to 3.2 times higher SBPase content compared to wild-type plants. Transcripts of the transgene were found predominantly in the bundle sheaths suggesting the correct cellular localization of the protein. The abundance of ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit was not affected in transgenic plants overexpressing SBPase, and neither was leaf chlorophyll content or photosynthetic electron transport parameters. We found no association between SBPase content in S. viridis and saturating rates of CO2 assimilation. Moreover, a detailed analysis of CO2 assimilation rates at different CO2 partial pressures, irradiances, and leaf temperatures showed no improvement of photosynthesis in plants overexpressing SBPase. We discuss the potential implications of these results for understanding the role of SBPase in regulation of C4 photosynthesis.

Identification and haplotype analysis of SiCHLI: a gene for yellow-green seedling as morphological marker to accelerate foxtail millet (Setaria italica) hybrid breeding

We cloned and developed functional markers for the SiCHLI gene, which is responsible for the yellow-green color of leaves in foxtail millet, a frequently used marker trait in the hybrid breeding of foxtail millet by using bulked segregant analysis sequencing and haplotype analysis on the F₂ and core-collected nature populations. The color of leaves has been widely used as a marker for the hybrid breeding of foxtail millet; however, few related gene have been cloned to date. Here, we used two F₂ populations generated from crosses between the highly male-sterile material 125A with yellow-green leaves, and CG58 and S410, which have green leaves, to identify the genes underlying the yellow-green color of the leaves of foxtail millet. The leaves of 125A seedlings were yellow-green, but they became green at the heading stage. The content of chlorophyll a and chlorophyll b was lower, the number of thylakoid lamellae and grana was reduced, and the chloroplasts was more rounded in 125A than in S410 at the yellow-green leaf stage; however, no differences were observed between 125A and S410 in these traits and photosynthetic at the heading stage. Bulked segregant analysis and map-based cloning revealed that the SiCHLI gene is responsible for the leaf colors of 125A. A nonsynonymous mutation (C/T) in exon 3 causes yellow-green leaves in 125A at the seedling stage. Haplotype analysis of the SiCHLI gene in 596 core collected accessions revealed a new haplotype associated with high photosynthetic metabolic potential at the heading and mature stages, which could be used to enhance sterile lines with yellow-green leaves. We developed a functional marker that will facilitate the identification of foxtail millet accessions with the different types of yellow-green leaves. Generally, our study provides new genetic resources to guide the future marker-assisted or target-base editing in foxtail millet hybrid breeding.

Computational Docking Reveals Co-Evolution of C4 Carbon Delivery Enzymes in Diverse Plants

Proteins are modular functionalities regulating multiple cellular activities in prokaryotes and eukaryotes. As a consequence of higher plants adapting to arid and thermal conditions, C4 photosynthesis is the carbon fixation process involving multi-enzymes working in a coordinated fashion. However, how these enzymes interact with each other and whether they co-evolve in parallel to maintain interactions in different plants remain elusive to date. Here, we report our findings on the global protein co-evolution relationship and local dynamics of co-varying site shifts in key C4 photosynthetic enzymes. We found that in most of the selected key C4 photosynthetic enzymes, global pairwise co-evolution events exist to form functional couplings. Besides, protein-protein interactions between these enzymes may suggest their unknown functionalities in the carbon delivery process. For PEPC and PPCK regulation pairs, pocket formation at the interactive interface are not necessary for their function. This feature is distinct from another well-known regulation pair in C4 photosynthesis, namely, PPDK and PPDK-RP, where the pockets are necessary. Our findings facilitate the discovery of novel protein regulation types and contribute to expanding our knowledge about C4 photosynthesis.

Progress and prospects of C4 trait engineering in plants

Incorporating C₄ photosynthetic traits into C₃ crops is a rational approach for sustaining future demands for crop productivity. Using classical plant breeding, engineering this complex trait is unlikely to achieve its target. Therefore, it is critical and timely to implement novel biotechnological crop improvement strategies to accomplish this goal. However, a fundamental understanding of C₃ , C₄ , and C₃ -C₄ intermediate metabolism is crucial for the targeted use of biotechnological tools. This review assesses recent progress towards engineering C₄ photosynthetic traits in C₃ crops. We also discuss lessons learned from successes and failures of recent genetic engineering attempts in C₃ crops, highlighting the pros and cons of using rice as a model plant for short-, medium- and long-term goals of genetic engineering. This review provides an integrated approach towards engineering improved photosynthetic efficiency in C₃ crops for sustaining food, fibre and fuel production around the globe.

Pyruvate, phosphate dikinase regulatory protein impacts light response of C4 photosynthesis in Setaria viridis

In C4 plants, the pyruvate (Pyr), phosphate dikinase regulatory protein (PDRP) regulates the activity of the C4 pathway enzyme Pyr, phosphate dikinase (PPDK) in a light-/dark-dependent manner. The importance of this regulatory action to C4 pathway function and overall C4 photosynthesis is unknown. To resolve this question, we assessed in vivo PPDK phospho-regulation and whole leaf photophysiology in a CRISPR-Cas9 PDRP knockout (KO) mutant of the NADP-ME C4 grass green millet (Setaria viridis). PDRP enzyme activity was undetectable in leaf extracts from PDRP KO lines. Likewise, PPDK phosphorylated at the PDRP-regulatory Thr residue was immunologically undetectable in leaf extracts. PPDK enzyme activity in rapid leaf extracts was constitutively high in the PDRP KO lines, irrespective of light or dark pretreatment of leaves. Gas exchange analysis of net CO2 assimilation revealed PDRP KO leaves had markedly slower light induction kinetics when leaves transition from dark to high-light or low-light to high-light. In the initial 30 min of the light induction phase, KO leaves had an ∼15% lower net CO2 assimilation rate versus the wild-type (WT). Despite the impaired slower induction kinetics, we found growth and vigor of the KO lines to be visibly indistinguishable from the WT when grown in normal air and under standard growth chamber conditions. However, the PDRP KO plants grown under a fluctuating light regime exhibited a gradual multi-day decline in Fv/Fm, indicative of progressive photosystem II damage due to the absence of PDRP. Collectively, our results demonstrate that one of PDRP's functions in C4 photosynthesis is to ensure optimal photosynthetic light induction kinetics during dynamic changes in incident light.

Reduction of bundle sheath size boosts cyclic electron flow in C4 Setaria viridis acclimated to low light

When C₄ leaves are exposed to low light, the CO₂ concentration in the bundle sheath (BS) cells decreases, causing an increase in photorespiration relative to assimilation, and a consequent reduction in biochemical efficiency. These effects can be mitigated by complex acclimation syndromes, which are of primary importance for crop productivity but are not well studied. We unveil an acclimation strategy involving the coordination of electron transport processes. First, we characterize the anatomy, gas exchange and electron transport of C₄ Setaria viridis grown under low light. Through a purposely developed biochemical model, we resolve the photon fluxes and reaction rates to explain how the concerted acclimation strategies sustain photosynthetic efficiency. Our results show that a smaller BS in low-light-grown plants limited leakiness (the ratio of CO₂ leak rate out of the BS over the rate of supply via C₄ acid decarboxylation) but sacrificed light harvesting and ATP production. To counter ATP shortage and maintain high assimilation rates, plants facilitated light penetration through the mesophyll and upregulated cyclic electron flow in the BS. This shade tolerance mechanism, based on the optimization of light reactions, is possibly more efficient than the known mechanisms involving the rearrangement of carbon metabolism, and could potentially lead to innovative strategies for crop improvement.

Genome‑wide identification, phylogenetic and expression pattern analysis of GATA family genes in foxtail millet (Setaria italica)

Background

Transcription factors (TFs) play important roles in plants. Among the major TFs, GATA plays a crucial role in plant development, growth, and stress responses. However, there have been few studies on the GATA gene family in foxtail millet (Setaria italica). The release of the foxtail millet reference genome presents an opportunity for the genome-wide characterization of these GATA genes.

Results

In this study, we identified 28 GATA genes in foxtail millet distributed on seven chromosomes. According to the classification method of GATA members in Arabidopsis, SiGATA was divided into four subfamilies, namely subfamilies I, II, III, and IV. Structural analysis of the SiGATA genes showed that subfamily III had more introns than other subfamilies, and a large number of cis-acting elements were abundant in the promoter region of the SiGATA genes. Three tandem duplications and five segmental duplications were found among SiGATA genes. Tissue-specific results showed that the SiGATA genes were mainly expressed in foxtail millet leaves, followed by peels and seeds. Many genes were significantly induced under the eight abiotic stresses, such as SiGATA10, SiGATA16, SiGATA18, and SiGATA25, which deserve further attention.

Conclusions

Collectively, these findings will be helpful for further in-depth studies of the biological function of SiGATA, and will provide a reference for the future molecular breeding of foxtail millet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08786-0.

MicroRNAs expression profiles in early responses to different levels of water deficit in Setaria viridis

Water deficit is a major constraint for crops of economic importance in almost all agricultural regions. However, plants have an active defense system to adapt to these adverse conditions, acting in the reprogramming of gene expression responsible for encoding microRNAs (miRNAs). These miRNAs promote the regulation to the target gene expression by the post-transcriptional (PTGS) and transcriptional gene silencing (TGS), modulating several pathways including defense response to water deficit. The broader knowledge of the miRNA expression profile and its regulatory networks in response to water deficit can provide evidence for the development of new biotechnological tools for genetic improvement of several important crops. In this study, we used Setaria viridis accession A10.1 as a C4 model plant to widely investigate the miRNA expression profile in early responses to different levels of water deficit. Ecophysiological studies in Setaria viridis under water deficit and after rewatering demonstrated a drought tolerant accession, capable of a rapid recovery from the stress. Deep small RNA sequencing and degradome studies were performed in plants submitted to drought to identify differentially expressed miRNA genes and their predicted targets, using in silico analysis. Our findings showed that several miRNAs were differentially modulated in response to distinctive levels of water deficit and after rewatering. The predicted mRNA targets mainly corresponded to genes related to cell wall remodeling, antioxidant system and drought-related transcription factors, indicating that these genes are rapidly regulated in early responses to drought stress. The implications of these modulations are extensively discussed, and higher-effect miRNAs are suggested as major players for potential use in genetic engineering to improve drought tolerance in economically important crops, such as sugarcane, maize, and sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01226-z.

SvSTL1 in the large subunit family of ribonucleotide reductases plays a major role in chloroplast development of Setaria viridis

Ribonucleotide reductases (RNRs) are essential enzymes in DNA synthesis. However, little is known about the RNRs in plants. Here, we identified a svstl1 mutant from the self-created ethyl methanesulfonate (EMS) mutant library of Setaria viridis. The mutant leaves exhibited a bleaching phenotype at the heading stage. Paraffin section analysis showed the destruction of the C₄ Kranz anatomy. Transmission electron microscopy results further demonstrated the severely disturbed development of some chloroplasts. MutMap analysis revealed that the SvSTL1 gene is the primary candidate, encoding a large subunit of RNRs. Complementation experiments confirmed that SvSTL1 is responsible for the phenotype of svstl1. There are two additional RNR large subunit homologs in S. viridis, SvSTL2 and SvSTL3. To further understand the functions of these three RNR large subunit genes, a series of mutants were generated via CRISPR/Cas9 technology. In striking contrast to the finding that all three SvSTLs interact with the RNR small subunit, the phenotype varied along with the copies of chloroplast genome among different svstl single mutants: the svstl1 mutant exhibited pronounced chloroplast development and significantly fewer copies of the chloroplast genome than the svstl2 or svstl3 single mutants. These results suggested that SvSTL1 plays a major role in the optimal function of RNRs and is essential for chloroplast development. Furthermore, through the analysis of double and triple mutants, the study provides new insights into the finely tuned coordination among SvSTLs to maintain normal chloroplast development in the emerging C₄ model plant S. viridis.

Isolation of gametes and zygotes from Setaria viridis

Setaria viridis, the wild ancestor of foxtail millet (Setaria italica), is an effective model plant for larger C₄ crops because S. viridis has several desirable traits, such as short generation time, prolific seed production and a small genome size. These advantages are well suited for investigating molecular mechanisms in angiosperms, especially C₄ crop species. Here, we report a procedure for isolating gametes and zygotes from S. viridis flowers. To isolate egg cells, ovaries were harvested from unpollinated mature flowers and cut transversely, which allowed direct access to the embryo sac. Thereafter, an egg cell was released from the cut end of the basal portion of the dissected ovary. To isolate sperm cells, pollen grains released from anthers were immersed in a mannitol solution, resulting in pollen-grain bursting, which released sperm cells. Additionally, S. viridis zygotes were successfully isolated from freshly pollinated flowers. Isolated zygotes cultured in a liquid medium developed into globular-like embryos and cell masses. Thus, isolated S. viridis gametes, zygotes and embryos are attainable for detailed observations and investigations of fertilization and developmental events in angiosperms.

Characterizing serotonin biosynthesis in Setaria viridis leaves and its effect on aphids

A combined transcriptomic and metabolic analysis of Setaria viridis leaves responding to aphid infestation was used to identify genes related to serotonin biosynthesis. Setaria viridis (green foxtail), a short life-cycle C4 plant in the Poaceae family, is the wild ancestor of Setaria italica (foxtail millet), a resilient crop that provides good yields in dry and marginal land. Although S. viridis has been studied extensively in the last decade, the molecular mechanisms of insect resistance in this species remain under-investigated. To address this issue, we performed a metabolic analysis of S. viridis and discovered that these plants accumulate the tryptophan-derived compounds tryptamine and serotonin. To elucidate the defensive functions of serotonin, Rhophalosiphum padi (bird cherry-oat aphids) were exposed to this compound, either by exogenous application to the plant medium or with artificial diet bioassays. In both cases, exposure to serotonin increased aphid mortality. To identify genes that are involved in serotonin biosynthesis, we conducted a transcriptome analysis and identified several predicted S. viridis tryptophan decarboxylase (TDC) and tryptamine 5-hydroxylase (T5H) genes. Two candidate genes were ectopically expressed in Nicotiana tabacum, where SvTDC1 (Sevir.6G066200) had tryptophan decarboxylase activity, and SvT5H1 (Sevir.8G219600) had tryptamine hydroxylase activity. Moreover, the function of the SvTDC1 gene was validated using virus-induced gene silencing in S. italica, which caused a reduction in serotonin levels. This study provides the first evidence of serotonin biosynthesis in Setaria leaves. The biosynthesis of serotonin may play an important role in defense responses and could prove to be useful for developing more pest-tolerant Setaria italica cultivars.

Genome-wide loss of CHH methylation with limited transcriptome changes in Setaria viridis DOMAINS REARRANGED METHYLTRANSFERASE (DRM) mutants

The DOMAINS REARRANGED METHYLTRANSFERASEs (DRMs) are crucial for RNA-directed DNA methylation (RdDM) in plant species. Setaria viridis is a model monocot species with a relatively compact genome that has limited transposable element (TE) content. CRISPR-based genome editing approaches were used to create loss-of-function alleles for the two putative functional DRM genes in S. viridis to probe the role of RdDM. Double mutant (drm1ab) plants exhibit some morphological abnormalities but are fully viable. Whole-genome methylation profiling provided evidence for the widespread loss of methylation in CHH sequence contexts, particularly in regions with high CHH methylation in wild-type plants. Evidence was also found for the locus-specific loss of CG and CHG methylation, even in some regions that lack CHH methylation. Transcriptome profiling identified genes with altered expression in the drm1ab mutants. However, the majority of genes with high levels of CHH methylation directly surrounding the transcription start site or in nearby promoter regions in wild-type plants do not have altered expression in the drm1ab mutant, even when this methylation is lost, suggesting limited regulation of gene expression by RdDM. Detailed analysis of the expression of TEs identified several transposons that are transcriptionally activated in drm1ab mutants. These transposons are likely to require active RdDM for the maintenance of transcriptional repression.

Genome-wide identification and characterization of the SPL gene family and its expression in the various developmental stages and stress conditions in foxtail millet (Setaria italica)

BACKGROUND

Among the major transcription factors, SPL plays a crucial role in plant growth, development, and stress response. Foxtail millet (Setaria italica), as a C4 crop, is rich in nutrients and is beneficial to human health. However, research on the foxtail millet SPL (SQUAMOSA PROMOTER BINDING-LIKE) gene family is limited.  RESULTS: In this study, a total of 18 SPL genes were identified for the comprehensive analysis of the whole genome of foxtail millet. These SiSPL genes were divided into seven subfamilies (I, II, III, V, VI, VII, and VIII) according to the classification of the Arabidopsis thaliana SPL gene family. Structural analysis of the SiSPL genes showed that the number of introns in subfamilies I and II were much larger than others, and the promoter regions of SiSPL genes were rich in different cis-acting elements. Among the 18 SiSPL genes, nine genes had putative binding sites with foxtail millet miR156. No tandem duplication events were found between the SiSPL genes, but four pairs of segmental duplications were detected. The SiSPL genes expression were detected in different tissues, which was generally highly expressed in seeds development process, especially SiSPL6 and SiSPL16, which deserve further study. The results of the expression levels of SiSPL genes under eight types of abiotic stresses showed that many stress responsive genes, especially SiSPL9, SiSPL10, and SiSPL16, were highly expressed under multiple stresses, which deserves further attention.

CONCLUSIONS

In this research, 18 SPL genes were identified in foxtail millet, and their phylogenetic relationships, gene structural features, duplication events, gene expression and potential roles in foxtail millet development were studied. The findings provide a new perspective for the mining of the excellent SiSPL gene and the molecular breeding of foxtail millet.

Time-series transcriptomics reveals a drought-responsive temporal network and crosstalk between drought stress and the circadian clock in foxtail millet

Drought stress is a serious factor affecting crop growth and production worldwide. The circadian clock has been identified as key to improving regional adaptability of plants. However, our understanding of the contribution of the circadian clock to drought response and the impacts of drought stress on the circadian clock in plants is still limited. To explore the interactions between the circadian clock and drought stress, foxtail millet seedlings were treated with simulated drought (20% polyethylene glycol-6000) treatment starting at the day (DD) onset zeitgeber time 0 (ZT0, lights on) and at the night (DN) onset zeitgeber time 16 (ZT16, lights off). A high temporal-resolution transcriptomic investigation was performed using DD and DN samples collected at intervals of 2 or 4 h within a 24-h drought-treatment period. Overall, we identified 13 294 drought-responsive genes (DRGs). Among these DRGs, 7931 were common between DD and DN samples, 2638 were specific to DD, and 2725 were specific to DN. Additionally, we identified 1257 circadian genes, of which 67% were DRGs. Interestingly, with drought treatment starting at the day for 8, 12 or 16 h, the circadian phase shifted to 12 h. We also found that the circadian clock led to different day and night drought-responsive pathways. The identification of DRG_Clock (DRG and circadian clock) and DRG_NonClock (DRG and not circadian clock) genes provides a reference for selecting candidate drought resistance genes. Our work reveals the temporal drought-response process and crosstalk between drought stress and the circadian clock in foxtail millet.

Genome-wide identification of quantitative trait loci for morpho-agronomic and yield-related traits in foxtail millet (Setaria italica) across multi-environments

Foxtail millet (Setaria italica) is an ideal model of genetic system for functional genomics of the Panicoideae crop. Identification of QTL responsible for morpho-agronomic and yield-related traits facilitates dissection of genetic control and breeding in cereal crops. Here, based on a Yugu1 × Longgu7 RIL population and genome-wide resequencing data, an updated linkage map harboring 2297 bin and 74 SSR markers was constructed, spanning 1315.1 cM with an average distance of 0.56 cM between adjacent markers. A total of 221 QTL for 17 morpho-agronomic and yield-related traits explaining 5.5 ~ 36% of phenotypic variation were identified across multi-environments. Of these, 109 QTL were detected in two to nine environments, including the most stable qLMS6.1 harboring a promising candidate gene Seita.6G250500, of which 70 were repeatedly identified in different trials in the same geographic location, suggesting that foxtail millet has more identical genetic modules under the similar ecological environment. One hundred-thirty QTL with overlapping intervals formed 22 QTL clusters. Furthermore, six superior recombinant inbred lines, RIL35, RIL48, RIL77, RIL80, RIL115 and RIL125 with transgressive inheritance and enrichment of favorable alleles in plant height, tiller, panicle morphology and yield related-traits were screened by hierarchical cluster. These identified QTL, QTL clusters and superior lines lay ground for further gene-trait association studies and breeding practice in foxtail millet.

Leaf physiological traits of plants from the Qinghai-Tibet Plateau and other arid sites in China: Identifying susceptible species and well-adapted extremophiles

Extreme environments, such as deserts and high-elevation ecosystems, are very important from biodiversity and ecological perspectives. However, plant physiology at those sites has been scarcely studied, likely due to logistic difficulties. In the present study, leaf physiological traits in native plants were analyzed from arid zones across an elevational transect in Western China, from Turpan Basin to the Qinghai-Tibet Plateau (QTP) at Delingha. The aim of this study was to use leaf physiological traits to help identifying potentially threatened species and true extremophiles. Physiological measurements in the field, and particularly in situ measurements of gas exchange and chlorophyll fluorescence, have been determined to be useful to determine the current state of plants at a given environment. Using this approach plus a combination of leaf traits, several species performing particularly well at the QTP were identified, e.g. Hedysarum multijugum, as well as at Manas drylands, e.g. Peganum harmala and Setaria viridis. On the other hand, several species showed marked signs of severe stress, in particular a very low photosynthetic rate over its potential maximum, as well as other negative traits, like low water and/or nitrogen-use-efficiency, which should be considered in conservation plans. Interestingly, all C₄ species studied except Setaria viridis were among the most stressed species. Despite their higher water use efficiency and drought-tolerance reputation, they presented a much larger photosynthesis depression than most C₃ species. This is an intriguing and interesting observation that deserves further studies.

Future Outlook of Transferring Biological Nitrogen Fixation (BNF) to Cereals and Challenges to Retard Achieving this Dream

BNF is a fascinating phenomenon which contributes to protect the nature from environmental pollution that can be happened as a result of heavy nitrogen applications. The importance of BNF is due to its supply of the agricultural lands with about 200 million tons of N annually. In this biological process, a specific group of bacteria collectively called rhizobia fix the atmospheric N in symbiosis with legumes called symbiotic nitrogen fixation and others (free living) fix nitrogen gas from the atmosphere termed asymbiotic. Several trials were done by scientists around the world to make cereals more benefited from nitrogen gas through different approaches. The first approach is to engineer cereals to form nodulated roots. Secondly is to transfer nif genes directly to cereals and fix N without Rhizobium partner. The other two approaches are maximizing the inoculation of cereals with both of diazotrophs or endophytes. Recently, scientists solved some challenges that entangle engineering cereals with nif genes directly and they confirmed the suitability of mitochondria and plastids as a suitable place for better biological function of nif genes expression in cereals. Fortunately, this article is confirming the success of scientists not only to transfer synthetic nitrogenase enzyme to Escherichia coli that gave 50% of its activity of expression, but also move it to plants as Nicotiana benthamiana. This mini review aims at explaining the future outlook of BNF and the challenges limiting its transfer to cereals and levels of success to make cereals self nitrogen fixing.

The role of AUX1 during lateral root development in the domestication of the model C4 grass Setaria italica

C4 photosynthesis increases the efficiency of carbon fixation by spatially separating high concentrations of molecular oxygen from Rubisco. The specialized leaf anatomy required for this separation evolved independently many times. The morphology of C4 root systems is also distinctive and adapted to support high rates of photosynthesis; however, little is known about the molecular mechanisms that have driven the evolution of C4 root system architecture. Using a mutant screen in the C4 model plant Setaria italica, we identify Siaux1-1 and Siaux1-2 as root system architecture mutants. Unlike in S. viridis, AUX1 promotes lateral root development in S. italica. A cell by cell analysis of the Siaux1-1 root apical meristem revealed changes in the distribution of cell volumes in all cell layers and a dependence of the frequency of protophloem and protoxylem strands on SiAUX1. We explore the molecular basis of the role of SiAUX1 in seedling development using an RNAseq analysis of wild-type and Siaux1-1 plants and present novel targets for SiAUX1-dependent gene regulation. Using a selection sweep and haplotype analysis of SiAUX1, we show that Hap-2412TT in the promoter region of SiAUX1 is an allele which is associated with lateral root number and has been strongly selected for during Setaria domestication.

Elbaiomy R, Chen L, Wang Z, Jiao J, Zhu J, Zhou W, Chen B, Soaud SA, Abbas M, Lin N, El-Sappah AH. Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae)

FRIZZY PANICLE (FZP), an essential gene that controls spikelet differentiation and development in the grass family (Poaceae), prevents the formation of axillary bud meristems and is closely associated with crop yields. It is unclear whether the FZP gene or its orthologs were selected during the evolutionary process of grass species, which possess diverse spike morphologies. In the present study, we adopted bioinformatics methods for the evolutionary analysis of FZP orthologs in species of the grass family. Thirty-five orthologs with protein sequences identical to that of the FZP gene were identified from 29 grass species. Analysis of conserved domains revealed that the AP2/ERF domains were highly conserved with almost no amino acid mutations. However, species of the tribe Triticeae, genus Oryza, and C4 plants exhibited more significant amino acid mutations in the acidic C-terminus region. Results of the phylogenetic analysis showed that the 29 grass species could be classified into three groups, namely, Triticeae, Oryza, and C4 plants. Within the Triticeae group, the FZP genes originating from the same genome were classified into the same sub-group. When selection pressure analysis was performed, significant positive selection sites were detected in species of the Triticeae and Oryza groups. Our results show that the FZP gene was selected during the grass family's evolutionary process, and functional divergence may have already occurred among the various species. Therefore, researchers investigating the FZP gene's functions should take note of the possible presence of various roles in other grass species.

Role of Cytosolic Malic Enzyme in Oleaginicity of High-Lipid-Producing Fungal Strain Mucor circinelloides WJ11

Mucor circinelloides, an oleaginous filamentous fungus, is gaining popularity due to its ability to synthesize significant amounts of lipids containing γ-linolenic acid (GLA) that have important health benefits. Malic enzyme (ME), which serves as the main source of NADPH in some fungi, has been found to regulate lipid accumulation in oleaginous fungi. In the present study, the role of two cytosolic ME genes, cmalA and cmalB, in the lipid accumulation of the M. circinelloides high-lipid-producing strain WJ11, was evaluated. Strains overexpressing cmalA and cmalB showed a 9.8- and 6.4-fold rise in specific ME activity, respectively, and an elevation of the lipid content by 23.2% and 5.8%, respectively, suggesting that these genes are involved in lipid biosynthesis. Due to increased lipid accumulation, overall GLA content in biomass was observed to be elevated by 11.42% and 16.85% in cmalA and cmalB overexpressing strains, respectively. Our study gives an important insight into different studies exploring the role of the cmalA gene, while we have for the first time investigated the role of the cmalB gene in the M. circinelloides WJ11 strain.

Improvement of millets in the post-genomic era

Millets are food and nutrient security crops in the semi-arid tropics of developing countries. Crop improvement using modern tools is one of the priority areas of research in millets. The whole-genome sequence (WGS) of millets provides new insight into understanding and studying the genes, genome organization and genomic-assisted improvement of millets. The WGS of millets helps to carry out genome-wide comparison and co-linearity studies among millets and other cereal crops. This approach might lead to the identification of genes underlying biotic and abiotic stress tolerance in millets. The available genome sequence of millets can be used for SNP identification, allele discovery, association and linkage mapping, identification of valuable candidate genes, and marker-assisted breeding (MAB) programs. Next generation sequencing (NGS) technology provides opportunities for genome-assisted breeding (GAB) through genomic selection (GS) and genome-wide association studies (GAWS) for crop improvement. Clustered, regularly interspaced, short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) genome editing (GE) system provides new opportunities for millet improvement. In this review, we discuss the details on the WGS available for millets and highlight the importance of utilizing such resources in the post-genomic era for millet improvement. We also draw inroads on the utilization of various approaches such as GS, GWAS, functional genomics, gene validation and GE for millet improvement. This review might be helpful for understanding the developments in the post-genomic era of millet improvement.

Deciphering the microbial and molecular responses of geographically diverse Setaria accessions grown in a nutrient-poor soil

The microbial and molecular characterization of the ectorhizosphere is an important step towards developing a more complete understanding of how the cultivation of biofuel crops can be undertaken in nutrient poor environments. The ectorhizosphere of Setaria is of particular interest because the plant component of this plant-microbe system is an important agricultural grain crop and a model for biofuel grasses. Importantly, Setaria lends itself to high throughput molecular studies. As such, we have identified important intra- and interspecific microbial and molecular differences in the ectorhizospheres of three geographically distant Setaria italica accessions and their wild ancestor S. viridis. All were grown in a nutrient-poor soil with and without nutrient addition. To assess the contrasting impact of nutrient deficiency observed for two S. italica accessions, we quantitatively evaluated differences in soil organic matter, microbial community, and metabolite profiles. Together, these measurements suggest that rhizosphere priming differs with Setaria accession, which comes from alterations in microbial community abundances, specifically Actinobacteria and Proteobacteria populations. When globally comparing the metabolomic response of Setaria to nutrient addition, plants produced distinctly different metabolic profiles in the leaves and roots. With nutrient addition, increases of nitrogen containing metabolites were significantly higher in plant leaves and roots along with significant increases in tyrosine derived alkaloids, serotonin, and synephrine. Glycerol was also found to be significantly increased in the leaves as well as the ectorhizosphere. These differences provide insight into how C₄ grasses adapt to changing nutrient availability in soils or with contrasting fertilization schemas. Gained knowledge could then be utilized in plant enhancement and bioengineering efforts to produce plants with superior traits when grown in nutrient poor soils.

Ground tissue circuitry regulates organ complexity in maize and Setaria

[Figure: see text].

Evaluation of the molecular and physiological response to dehydration of two accessions of the model plant Setaria viridis

Water deficits are responsible for countless agricultural losses. Among the affected crops, C4 plants are of special interest due to their high water and nitrogen use efficiency. Two accessions of Setaria viridis (Ast-1 and A10.1) with contrasting responses to water deficit were used in the current work to better understand the mechanisms behind drought tolerance in C4 species. Our results showed that although the A10.1 accession exhibited a reduced size and lower Rfd values in comparison to Ast-1, it had overall higher Fv/Fm and lower NPQ values in well-watered conditions. The water deficit induction was performed with PEG-8000 at the grain-filling stage using dehydration cycles. Analysis of physiological measurements showed the A10.1 accession as being more tolerant to multiple water deficit exposures. In addition, PCA identified a clear difference in the pattern of drought response of the accessions. Four drought marker genes previously described in the literature were chosen to evaluate the response at the molecular level: SvP5CS2, SvDHN1, SvNAC6, and SvWRKY1. Besides confirming that Ast-1 is a more sensitive accession, the expression analysis revealed that SvNAC1 might better monitor drought stress, while SvWRKY1 was able to differentiate the two accessions. Distinct evolutionary histories of each accession may be behind their differences in response to water deficits.

Expression of a CO2-permeable aquaporin enhances mesophyll conductance in the C4 species Setaria viridis

A fundamental limitation of photosynthetic carbon fixation is the availability of CO₂. In C₄ plants, primary carboxylation occurs in mesophyll cytosol, and little is known about the role of CO₂ diffusion in facilitating C₄ photosynthesis. We have examined the expression, localization, and functional role of selected plasma membrane intrinsic aquaporins (PIPs) from Setaria italica (foxtail millet) and discovered that SiPIP2;7 is CO₂-permeable. When ectopically expressed in mesophyll cells of Setaria viridis (green foxtail), SiPIP2;7 was localized to the plasma membrane and caused no marked changes in leaf biochemistry. Gas exchange and C¹⁸O¹⁶O discrimination measurements revealed that targeted expression of SiPIP2;7 enhanced the conductance to CO₂ diffusion from the intercellular airspace to the mesophyll cytosol. Our results demonstrate that mesophyll conductance limits C₄ photosynthesis at low pCO₂ and that SiPIP2;7 is a functional CO₂ permeable aquaporin that can improve CO₂ diffusion at the airspace/mesophyll interface and enhance C₄ photosynthesis.

Challenges and Approaches to Crop Improvement Through C3-to-C4 Engineering

With a rapidly growing world population and dwindling natural resources, we are now facing the enormous challenge of increasing crop yields while simultaneously improving the efficiency of resource utilization. Introduction of C4 photosynthesis into C3 crops is widely accepted as a key strategy to meet this challenge because C4 plants are more efficient than C3 plants in photosynthesis and resource usage, particularly in hot climates, where the potential for productivity is high. Lending support to the feasibility of this C3-to-C4 engineering, evidence indicates that C4 photosynthesis has evolved from C3 photosynthesis in multiple lineages. Nevertheless, C3-to-C4 engineering is not an easy task, as several features essential to C4 photosynthesis must be introduced into C3 plants. One such feature is the spatial separation of the two phases of photosynthesis (CO₂ fixation and carbohydrate synthesis) into the mesophyll and bundle sheath cells, respectively. Another feature is the Kranz anatomy, characterized by a close association between the mesophyll and bundle sheath (BS) cells (1:1 ratio). These anatomical features, along with a C4-specific carbon fixation enzyme (PEPC), form a CO₂-concentration mechanism that ensures a high photosynthetic efficiency. Much effort has been taken in the past to introduce the C4 mechanism into C3 plants, but none of these attempts has met with success, which is in my opinion due to a lack of system-level understanding and manipulation of the C3 and C4 pathways. As a prerequisite for the C3-to-C4 engineering, I propose that not only the mechanisms that control the Kranz anatomy and cell-type-specific expression in C3 and C4 plants must be elucidated, but also a good understanding of the gene regulatory network underlying C3 and C4 photosynthesis must be achieved. In this review, I first describe the past and current efforts to increase photosynthetic efficiency in C3 plants and their limitations; I then discuss a systems approach to tackling down this challenge, some practical issues, and recent technical innovations that would help us to solve these problems.

Updating the steady-state model of C4 photosynthesis

C4 plants play a key role in world agriculture. For example, C4 crops such as maize and sorghum are major contributors to food production in both developed and developing countries, and the C4 grasses sugarcane, miscanthus, and switchgrass are major plant sources of bioenergy. In the challenge to manipulate and enhance C4 photosynthesis, steady-state models of leaf photosynthesis provide an important tool for gas exchange analysis and thought experiments that can explore photosynthetic pathway changes. Here a previous C4 photosynthetic model developed by von Caemmerer and Furbank has been updated with new kinetic parameterization and temperature dependencies added. The parameterization was derived from experiments on the C4 monocot, Setaria viridis, which for the first time provides a cohesive parameterization. Mesophyll conductance and its temperature dependence have also been included, as this is an important step in the quantitative correlation between the initial slope of the CO2 response curve of CO2 assimilation and in vitro phosphoenolpyruvate carboxylase activity. Furthermore, the equations for chloroplast electron transport have been updated to include cyclic electron transport flow, and equations have been added to calculate the electron transport rate from measured CO2 assimilation rates.

Updating the steady-state model of C4 photosynthesis

C4 plants play a key role in world agriculture. For example, C4 crops such as maize and sorghum are major contributors to food production in both developed and developing countries, and the C4 grasses sugarcane, miscanthus, and switchgrass are major plant sources of bioenergy. In the challenge to manipulate and enhance C4 photosynthesis, steady-state models of leaf photosynthesis provide an important tool for gas exchange analysis and thought experiments that can explore photosynthetic pathway changes. Here a previous C4 photosynthetic model developed by von Caemmerer and Furbank has been updated with new kinetic parameterization and temperature dependencies added. The parameterization was derived from experiments on the C4 monocot, Setaria viridis, which for the first time provides a cohesive parameterization. Mesophyll conductance and its temperature dependence have also been included, as this is an important step in the quantitative correlation between the initial slope of the CO2 response curve of CO2 assimilation and in vitro phosphoenolpyruvate carboxylase activity. Furthermore, the equations for chloroplast electron transport have been updated to include cyclic electron transport flow, and equations have been added to calculate the electron transport rate from measured CO2 assimilation rates.

Genome-wide identification of PTI1 family in Setaria italica and salinity-responsive functional analysis of SiPTI1-5

BACKGROUND

PTI1 (Pto-interacting 1) protein kinase belongs to the receptor-like cytoplasmic kinase (RLCK) group of receptor-like protein kinases (RLK), but lack extracellular and transmembrane domains. PTI1 was first identified in tomato (Solanum lycopersicum) and named SlPTI1, which has been reported to interact with bacterial effector Pto, a serine/threonine protein kinase involved in plant resistance to bacterial disease. Briefly, the host PTI1 specifically recognizes and interacts with the bacterial effector AvrPto, which triggers hypersensitive cell death to inhibit the pathogen growth in the local infection site. Previous studies have demonstrated that PTI1 is associated with oxidative stress and hypersensitivity.

RESULTS

We identified 12 putative PTI1 genes from the genome of foxtail millet (Setaria italica) in this study. Gene replication analysis indicated that both segmental replication events played an important role in the expansion of PTI1 gene family in foxtail millet. The PTI1 family members of model plants, i.e. S. italica, Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), maize (Zea mays), S. lycopersicum, and soybean (Glycine max), were classified into six major categories according to the phylogenetic analysis, among which the PTI1 family members in foxtail millet showed higher degree of homology with those of rice and maize. The analysis of a complete set of SiPTI1 genes/proteins including classification, chromosomal location, orthologous relationships and duplication. The tissue expression characteristics revealed that SiPTI1 genes are mainly expressed in stems and leaves. Experimental qRT-PCR results demonstrated that 12 SiPTI1 genes were induced by multiple stresses. Subcellular localization visualized that all of foxtail millet SiPTI1s were localized to the plasma membrane. Additionally, heterologous expression of SiPTI1-5 in yeast and E. coli enhanced their tolerance to salt stress.

CONCLUSIONS

Our results contribute to a more comprehensive understanding of the roles of PTI1 protein kinases and will be useful in prioritizing particular PTI1 for future functional validation studies in foxtail millet.

Enhanced aluminum tolerance in sugarcane: evaluation of SbMATE overexpression and genome-wide identification of ALMTs in Saccharum spp

BACKGROUND

A major limiting factor for plant growth is the aluminum (Al) toxicity in acidic soils, especially in tropical regions. The exclusion of Al from the root apex through root exudation of organic acids such as malate and citrate is one of the most ubiquitous tolerance mechanisms in the plant kingdom. Two families of anion channels that confer Al tolerance are well described in the literature, ALMT and MATE family.

RESULTS

In this study, sugarcane plants constitutively overexpressing the Sorghum bicolor MATE gene (SbMATE) showed improved tolerance to Al when compared to non-transgenic (NT) plants, characterized by sustained root growth and exclusion of aluminum from the root apex based on the result obtained with hematoxylin staining. In addition, genome-wide analysis of the recently released sugarcane genome identified 11 ALMT genes and molecular studies showed potential new targets for aluminum tolerance.

CONCLUSIONS

Our results indicate that the transgenic plants overexpressing the Sorghum bicolor MATE has an improved tolerance to Al. The expression profile of ALMT genes revels potential candidate genes to be used has an alternative for agricultural expansion in Brazil and other areas with aluminum toxicity in poor and acid soils.

Engineering Posttranslational Regulation of Glutamine Synthetase for Controllable Ammonia Production in the Plant Symbiont Azospirillum brasilense

Nitrogen requirements for modern agriculture far exceed the levels of bioavailable nitrogen in most arable soils. As a result, the addition of nitrogen fertilizer is necessary to sustain productivity and yields, especially for cereal crops, the planet's major calorie suppliers. Given the unsustainability of industrial fertilizer production and application, engineering biological nitrogen fixation directly at the roots of plants has been a grand challenge for biotechnology. Here, we designed and tested a potentially broadly applicable metabolic engineering strategy for the overproduction of ammonia in the diazotrophic symbiont Azospirillum brasilense. Our approach is based on an engineered unidirectional adenylyltransferase (uAT) that posttranslationally modifies and deactivates glutamine synthetase (GS), a key regulator of nitrogen metabolism in the cell. We show that this circuit can be controlled inducibly, and we leveraged the inherent self-contained nature of our posttranslational approach to demonstrate that multicopy redundancy can improve strain evolutionary stability. uAT-engineered Azospirillum is capable of producing ammonia at rates of up to 500 μM h-1 unit of OD600 (optical density at 600 nm)-1. We demonstrated that when grown in coculture with the model monocot Setaria viridis, these strains increase the biomass and chlorophyll content of plants up to 54% and 71%, respectively, relative to the wild type (WT). Furthermore, we rigorously demonstrated direct transfer of atmospheric nitrogen to extracellular ammonia and then plant biomass using isotopic labeling: after 14 days of cocultivation with engineered uAT strains, 9% of chlorophyll nitrogen in Setaria seedlings was derived from diazotrophically fixed dinitrogen, whereas no nitrogen was incorporated in plants cocultivated with WT controls. This rational design for tunable ammonia overproduction is modular and flexible, and we envision that it could be deployable in a consortium of nitrogen-fixing symbiotic diazotrophs for plant fertilization. IMPORTANCE Nitrogen is the most limiting nutrient in modern agriculture. Free-living diazotrophs, such as Azospirillum, are common colonizers of cereal grasses and have the ability to fix nitrogen but natively do not release excess ammonia. Here, we used a rational engineering approach to generate ammonia-excreting strains of Azospirillum. Our design features posttranslational control of highly conserved central metabolism, enabling tunability and flexibility of circuit placement. We found that our strains promote the growth and health of the model grass S. viridis and rigorously demonstrated that in comparison to WT controls, our engineered strains can transfer nitrogen from 15N2 gas to plant biomass. Unlike previously reported ammonia-producing mutants, our rationally designed approach easily lends itself to further engineering opportunities and has the potential to be broadly deployable.

Efficient genome editing and gene knockout in Setaria viridis with CRISPR/Cas9 directed gene editing by the non-homologous end-joining pathway

The CRISPR/Cas9 system has been used for genome editing in several organisms, including higher plants. This system induces site-specific mutations in the genome based on the nucleotide sequence of engineered guide RNAs. The complex genomes of C4 grasses makes genome editing a challenge in key grass crops like maize (Zea mays), sorghum (Sorghum bicolor), Brachiaria spp., switchgrass (Panicum virgatum), and sugarcane (Saccharum spp.). Setaria viridis is a diploid C4 grass widely used as a model for these C4 crop plants. Here, an optimized CRISPR/Cas9 binary vector that exploits the non-homologous end joining (NHEJ) system was used to knockout a green fluorescent protein (gfp) transgene in S. viridis accession A10.1. Transformation of embryogenic callus by A. tumefaciens generated ten glufosinate-ammonium resistant transgenic events. In the T0 generation, 60% of the events were biallelic mutants in the gfp transgene with no detectable accumulation of GFP protein and without insertions or deletions in predicted off-target sites. The gfp mutations generated by CRISPR/Cas9 were stable and displayed Mendelian segregation in the T1 generation. Altogether, the system described here is a highly efficient genome editing system for S. viridis, an important model plant for functional genomics studies in C4 grasses. Also, this system is a potential tool for improvement of agronomic traits in C4 crop plants with complex genomes.

Correlation and co-localization of QTL for stomatal density, canopy temperature, and productivity with and without drought stress in Setaria

Mechanistic modeling indicates that stomatal conductance could be reduced to improve water use efficiency (WUE) in C4 crops. Genetic variation in stomatal density and canopy temperature was evaluated in the model C4 genus, Setaria. Recombinant inbred lines (RILs) derived from a Setaria italica×Setaria viridis cross were grown with ample or limiting water supply under field conditions in Illinois. An optical profilometer was used to rapidly assess stomatal patterning, and canopy temperature was measured using infrared imaging. Stomatal density and canopy temperature were positively correlated but both were negatively correlated with total above-ground biomass. These trait relationships suggest a likely interaction between stomatal density and the other drivers of water use such as stomatal size and aperture. Multiple quantitative trait loci (QTL) were identified for stomatal density and canopy temperature, including co-located QTL on chromosomes 5 and 9. The direction of the additive effect of these QTL on chromosome 5 and 9 was in accordance with the positive phenotypic relationship between these two traits. This, along with prior experiments, suggests a common genetic architecture between stomatal patterning and WUE in controlled environments with canopy transpiration and productivity in the field, while highlighting the potential of Setaria as a model to understand the physiology and genetics of WUE in C4 species.

De novo genome assembly of a foxtail millet cultivar Huagu11 uncovered the genetic difference to the cultivar Yugu1, and the genetic mechanism of imazethapyr tolerance

BACKGROUND

Setaria italica is the second-most widely planted species of millets in the world and an important model grain crop for the research of C4 photosynthesis and abiotic stress tolerance. Through three genomes assembly and annotation efforts, all genomes were based on next generation sequencing technology, which limited the genome continuity.

RESULTS

Here we report a high-quality whole-genome of new cultivar Huagu11, using single-molecule real-time sequencing and High-throughput chromosome conformation capture (Hi-C) mapping technologies. The total assembly size of the Huagu11 genome was 408.37 Mb with a scaffold N50 size of 45.89 Mb. Compared with the other three reported millet genomes based on the next generation sequencing technology, the Huagu11 genome had the highest genomic continuity. Intraspecies comparison showed about 94.97 and 94.66% of the Yugu1 and Huagu11 genomes, respectively, were able to be aligned as one-to-one blocks with four chromosome inversion. The Huagu11 genome contained approximately 19.43 Mb Presence/absence Variation (PAV) with 627 protein-coding transcripts, while Yugu1 genomes had 20.53 Mb PAV sequences encoding 737 proteins. Overall, 969,596 Single-nucleotide polymorphism (SNPs) and 156,282 insertion-deletion (InDels) were identified between these two genomes. The genome comparison between Huagu11 and Yugu1 should reflect the genetic identity and variation between the cultivars of foxtail millet to a certain extent. The Ser-626-Aln substitution in acetohydroxy acid synthase (AHAS) was found to be relative to the imazethapyr tolerance in Huagu11.

CONCLUSIONS

A new improved high-quality reference genome sequence of Setaria italica was assembled, and intraspecies genome comparison determined the genetic identity and variation between the cultivars of foxtail millet. Based on the genome sequence, it was inferred that the Ser-626-Aln substitution in AHAS was responsible for the imazethapyr tolerance in Huagu11. The new improved reference genome of Setaria italica will promote the genic and genomic studies of this species and be beneficial for cultivar improvement.

Developing Portulaca oleracea as a model system for functional genomics analysis of C4/CAM photosynthesis

Previously regarded as an intriguing photosynthetic curiosity, the occurrence of C4 and Crassulacean acid metabolism (CAM) photosynthesis within a single organism has recently emerged as a source of information for future biotechnological use. Among C4/CAM facultative species, Portulaca oleracea L. has been used as a model for biochemical and gene expression analysis of C4/CAM under field and laboratory conditions. In the present work, we focussed on developing molecular tools to facilitate functional genomics studies in this species, from the optimisation of RNA isolation protocols to a method for stable genetic transformation. Eleven variations of RNA extraction procedures were tested and compared for RNA quantity and quality. Also, 7 sample sets comprising total RNA from hormonal and abiotic stress treatments, distinct plant organs, leaf developmental stages, and subspecies were used to select, among 12 reference genes, the most stable reference genes for RT-qPCR analysis of each experimental condition. Furthermore, different explant sources, Agrobacterium tumefaciens strains, and regeneration and antibiotic selection media were tested in various combinations to optimise a protocol for stable genetic transformation of P. oleracea. Altogether, we provide essential tools for functional gene analysis in the context of C4/CAM photosynthesis, including an efficient RNA isolation method, preferred reference genes for RT-qPCR normalisation for a range of experimental conditions, and a protocol to produce P. oleracea stable transformants using A. tumefaciens.

Agrobacterium-Mediated Transformation of Setaria viridis, a Model System for Cereals and Bioenergy Crops

Setaria viridis is an emerging model system for the genetic and molecular characterization of cereals and bioenergy crops. Here, we describe a detailed procedure for genetic transformation of the S. viridis accession line ME034V-1. This method utilizes callus generated from mature seeds for infection with Agrobacterium tumefaciens strain AGL1 to regenerate hygromycin-resistant stable transgenic plants. It takes approximately 7 weeks to generate callus from mature seeds, 11-17 weeks from infection to the regeneration of transgenic lines, and an additional 3-4 weeks for plant growth in the greenhouse for seed collection. The protocol as presented consistently results in transformation frequency of approximately 25% for the generation of transgenic plants, with fewer escapes and higher survivability in soil for optimal seed collection. © 2021 Donald Danforth Plant Science Center. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Generation of S. viridis (Accession ME034V-1) callus from mature seeds Basic Protocol 2: Agrobacterium-mediated transformation of callus to generate transgenic plants Basic Protocol 3: Plantlet transplantation in soil, plant growth in greenhouse, and seed collection Support Protocol: Preparation of Agrobacterium culture for infection.

Diazotrophic Bacteria and Their Mechanisms to Interact and Benefit Cereals

Plant-growth-promoting bacteria (PGPB) stimulate plant growth through diverse mechanisms. In addition to biological nitrogen fixation, diazotrophic PGPB can improve nutrient uptake efficiency from the soil, produce and release phytohormones to the host, and confer resistance against pathogens. The genetic determinants that drive the success of biological nitrogen fixation in nonlegume plants are understudied. These determinants include recognition and signaling pathways, bacterial colonization, and genotype specificity between host and bacteria. This review presents recent discoveries of how nitrogen-fixing PGPB interact with cereals and promote plant growth. We suggest adopting an experimental model system, such as the Setaria-diazotrophic bacteria association, as a reliable way to better understand the associated mechanisms and, ultimately, increase the use of PGPB inoculants for sustainable agriculture.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Shedding Light on the Dynamic Role of the "Target of Rapamycin" Kinase in the Fast-Growing C4 Species Setaria viridis, a Suitable Model for Biomass Crops

The Target of Rapamycin (TOR) kinase pathway integrates energy and nutrient availability into metabolism promoting growth in eukaryotes. The overall higher efficiency on nutrient use translated into faster growth rates in C₄ grass plants led to the investigation of differential transcriptional and metabolic responses to short-term chemical TOR complex (TORC) suppression in the model Setaria viridis. In addition to previously described responses to TORC inhibition (i.e., general growth arrest, translational repression, and primary metabolism reprogramming) in Arabidopsis thaliana (C₃), the magnitude of changes was smaller in S. viridis, particularly regarding nutrient use efficiency and C allocation and partitioning that promote biosynthetic growth. Besides photosynthetic differences, S. viridis and A. thaliana present several specificities that classify them into distinct lineages, which also contribute to the observed alterations mediated by TOR. Indeed, cell wall metabolism seems to be distinctly regulated according to each cell wall type, as synthesis of non-pectic polysaccharides were affected in S. viridis, whilst assembly and structure in A. thaliana. Our results indicate that the metabolic network needed to achieve faster growth seems to be less stringently controlled by TORC in S. viridis.

Potassium supply promotes the mitigation of NaCl-induced effects on leaf photochemistry, metabolism and morphology of Setaria viridis

Soil salinity has the potential to severely affect crop performance. To maintain cell functioning and improve salt tolerance, the maintenance of K⁺ homeostasis is crucial in several plant metabolism processes. Besides, potassium fertilization can efficiently alleviate the perilous effects of salinity. We characterized impacts in Setaria viridis exposed to NaCl and KCl to underlying photochemistry mechanisms, K⁺ and Na⁺ shoot contents, enzymatic activity, electrolytic leakage, and morphological responses focusing on non-stomatal limitation of photosynthesis. Plants were exposed to sodium chloride (NaCl; 0, 150 and 250 mM) and potassium chloride (KCl; 0, 5, 9 mM). The exposure to NaCl affected S. viridis leaves morphological and physiologically. Plants submitted to 150 mM showed reductions in performance indexes (PI_abs and PI_total; JIP-test), and the presence of positive K- and L-bands. Plants exposed to 250 mM exhibited blockage in electron flow further than Q_A within 48h and permanent photoinhibition at 96 h. The presence of 9 and 5 mM of KCl counteracted the effects of NaCl on plants submitted to 150 mM, concomitant with increases in K⁺ accumulation and cell turgidity conservation, causing positive effects in plant growth and metabolism. Neither KCl concentrations were effective in reducing NaCl-induced effects on plants exposed to 250 mM of NaCl. Our results support the conclusion that greater availability of K⁺ alleviates the harmful effects of salinity in S. viridis under moderate stress and that application of KCl as means of lightning saline stress has a concentration and a salt level limit that must be experimentally determined.