Conserved pleiotropy of an ancient plant homeobox gene uncovered by cis-regulatory dissection

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

In-locus gene silencing in plants using genome editing

Gene silencing is crucial in crop breeding for desired trait development. RNA interference (RNAi) has been used widely but is limited by ectopic expression of transgenes and genetic instability. Introducing an upstream start codon (uATG) into the 5'untranslated region (5'UTR) of a target gene may 'silence' the target gene by inhibiting protein translation from the primary start codon (pATG). Here, we report an efficient gene silencing method by introducing a tailor-designed uATG-containing element (ATGE) into the 5'UTR of genes in plants, occupying the original start site to act as a new pATG. Using base editing to introduce new uATGs failed to silence two of the tested three rice genes, indicating complex regulatory mechanisms. Precisely inserting an ATGE adjacent to pATG achieved significant target protein downregulation. Through extensive optimization, we demonstrated this strategy substantially and consistently downregulated target protein expression. By designing a bidirectional multifunctional ATGE4, we enabled tunable knockdown from 19% to 89% and observed expected phenotypes. Introducing ATGE into Waxy, which regulates starch synthesis, generated grains with lower amylose, revealing the value for crop breeding. Together, we have developed a programmable and robust method to knock down gene expression in plants, with potential for biological mechanism exploration and crop enhancement.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BACKGROUND

Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.

PURPOSE

Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.

RESULTS

Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.

CONCLUSION

Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

A genome-wide association study uncovers a ZmRap2.7-ZCN9/ZCN10 module to regulate ABA signalling and seed vigour in maize

Seed vigour, including rapid, uniform germination and robust seedling establishment under various field conditions, is becoming an increasingly essential agronomic trait for achieving high yield in crops. However, little is known about this important seed quality trait. In this study, we performed a genome-wide association study to identify a key transcription factor ZmRap2.7, which regulates seed vigour through transcriptionally repressing expressions of three ABA signalling genes ZmPYL3, ZmPP2C and ZmABI5 and two phosphatidylethanolamine-binding genes ZCN9 and ZCN10. In addition, ZCN9 and ZCN10 proteins could interact with ZmPYL3, ZmPP2C and ZmABI5 proteins, and loss-of-function of ZmRap2.7 and overexpression of ZCN9 and ZCN10 reduced ABA sensitivity and seed vigour, suggesting a complex regulatory network for regulation of ABA signalling mediated seed vigour. Finally, we showed that four SNPs in ZmRap2.7 coding region influenced its transcriptionally binding activity to the downstream gene promoters. Together with previously identified functional variants within and surrounding ZmRap2.7, we concluded that the distinct allelic variations of ZmRap2.7 were obtained independently during maize domestication and improvement, and responded separately for the diversities of seed vigour, flowering time and brace root development. These results provide novel genes, a new regulatory network and an evolutional mechanism for understanding the molecular mechanism of seed vigour.

The Evolution of the WUSCHEL-Related Homeobox Gene Family in Dendrobium Species and Its Role in Sex Organ Development in D. chrysotoxum

The WUSCHEL-related homeobox (WOX) transcription factor plays a vital role in stem cell maintenance and organ morphogenesis, which are essential processes for plant growth and development. Dendrobium chrysotoxum, D. huoshanense, and D. nobile are valued for their ornamental and medicinal properties. However, the specific functions of the WOX gene family in Dendrobium species are not well understood. In our study, a total of 30 WOX genes were present in the genomes of the three Dendrobium species (nine DchWOXs, 11 DhuWOXs, and ten DnoWOXs). These 30 WOXs were clustered into ancient clades, intermediate clades, and WUS/modern clades. All 30 WOXs contained a conserved homeodomain, and the conserved motifs and gene structures were similar among WOXs belonging to the same branch. D. chrysotoxum and D. huoshanense had one pair of fragment duplication genes and one pair of tandem duplication genes, respectively; D. nobile had two pairs of fragment duplication genes. The cis-acting regulatory elements (CREs) in the WOX promoter region were mainly enriched in the light response, stress response, and plant growth and development regulation. The expression pattern and RT-qPCR analysis revealed that the WOXs were involved in regulating the floral organ development of D. chrysotoxum. Among them, the high expression of DchWOX3 suggests that it might be involved in controlling lip development, whereas DchWOX5 might be involved in controlling ovary development. In conclusion, this work lays the groundwork for an in-depth investigation into the functions of WOX genes and their regulatory role in Dendrobium species' floral organ development.

Contribution of the transcription factor SfGATAe to Bt Cry toxin resistance in Spodoptera frugiperda through reduction of ABCC2 expression

Insect resistance evolution poses a significant threat to the advantages of biopesticides and transgenic crops utilizing insecticidal Cry-toxins from Bacillus thuringiensis (Bt). However, there is limited research on the relationship between transcriptional regulation of specific toxin receptors in lepidopteran insects and their resistance to Bt toxins. Here, we report the positive regulatory role of the SfGATAe transcription factor on the expression of the ABCC2 gene in Spodoptera frugiperda. DNA regions in the SfABCC2 promoter that are vital for regulation by SfGATAe, utilizing DAP-seq technology and promoter deletion mapping. Through yeast one-hybrid assays, DNA pull-down experiments, and site-directed mutagenesis, we confirmed that the transcription factor SfGATAe regulates the core control site PBS2 in the ABCC2 target gene. Tissue-specific expression analysis has revealed that SfGATAe is involved in the regulation and expression of midgut cells in the fall armyworm. Silencing SfGATAe in fall armyworm larvae resulted in reduced expression of SfABCC2 and decreased sensitivity to Cry1Ac toxin. Overall, this study elucidated the regulatory mechanism of the transcription factor SfGATAe on the expression of the toxin receptor gene SfABCC2 and this transcriptional control mechanism impacts the resistance of the fall armyworm to Bt toxins.

Developmental pleiotropy of SDP1 from seedling to mature stages in B. napus

Rapeseed, an important oil crop, relies on robust seedling emergence for optimal yields. Seedling emergence in the field is vulnerable to various factors, among which inadequate self-supply of energy is crucial to limiting seedling growth in early stage. SUGAR-DEPENDENT1 (SDP1) initiates triacylglycerol (TAG) degradation, yet its detailed function has not been determined in B. napus. Here, we focused on the effects of plant growth during whole growth stages and energy mobilization during seedling establishment by mutation in BnSDP1. Protein sequence alignment and haplotypic analysis revealed the conservation of SDP1 among species, with a favorable haplotype enhancing oil content. Investigation of agronomic traits indicated bnsdp1 had a minor impact on vegetative growth and no obvious developmental defects when compared with wild type (WT) across growth stages. The seed oil content was improved by 2.0-2.37% in bnsdp1 lines, with slight reductions in silique length and seed number per silique. Furthermore, bnsdp1 resulted in lower seedling emergence, characterized by a shrunken hypocotyl and poor photosynthetic capacity in the early stages. Additionally, impaired seedling growth, especially in yellow seedlings, was not fully rescued in medium supplemented with exogenous sucrose. The limited lipid turnover in bnsdp1 was accompanied by induced amino acid degradation and PPDK-dependent gluconeogenesis pathway. Analysis of the metabolites in cotyledons revealed active amino acid metabolism and suppressed lipid degradation, consistent with the RNA-seq results. Finally, we proposed strategies for applying BnSDP1 in molecular breeding. Our study provides theoretical guidance for understanding trade-off between oil accumulation and seedling energy mobilization in B. napus.

Engineering Plant Cell Fates and Functions for Agriculture and Industry

Many plant species are grown to enable access to specific organs or tissues, such as seeds, fruits, or stems. In some cases, a value is associated with a molecule that accumulates in a single type of cell. Domestication and subsequent breeding have often increased the yields of these target products by increasing the size, number, and quality of harvested organs and tissues but also via changes to overall plant growth architecture to suit large-scale cultivation. Many of the mutations that underlie these changes have been identified in key regulators of cellular identity and function. As key determinants of yield, these regulators are key targets for synthetic biology approaches to engineer new forms and functions. However, our understanding of many plant developmental programs and cell-type specific functions is still incomplete. In this Perspective, we discuss how advances in cellular genomics together with synthetic biology tools such as biosensors and DNA-recording devices are advancing our understanding of cell-specific programs and cell fates. We then discuss advances and emerging opportunities for cell-type-specific engineering to optimize plant morphology, responses to the environment, and the production of valuable compounds.

Disruption of Zfh3 abolishes mulberry-specific monophagy in silkworm larvae

Feeding behavior is critical for insect survival and fitness. Most researchers have explored the molecular basis of feeding behaviors by identifying and elucidating the function of olfactory receptors (ORs) and gustatory receptors (GRs). Other types of genes, such as transcription factors, have rarely been investigated, and little is known about their potential roles. The silkworm (Bombyx mori) is a well-studied monophagic insect which primarily feeds on mulberry leaves, but the genetic basis of its monophagy is still not understood. In this report, we focused on a transcription factor encoded by the Zfh3 gene, which is highly expressed in the silkworm central and peripheral nervous systems, including brain, antenna, and maxilla. To investigate its function, Zfh3 was abrogated using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) mutagenesis. Since Zfh3 knockout homozygotes are not viable, we studied feeding behavior in heterozygotes, and found that disruption of Zfh3 affects both gustation and olfaction. Mutant larvae lose preference for mulberry leaves, acquire the ability to consume an expanded range of diets, and exhibit improved adaptation to the M0 artificial diet, which contains no mulberry leaves. These results provide the first demonstration that a transcription factor modulates feeding behaviors in an insect.

WUSCHEL-related homeobox transcription factor SlWOX13 regulates tomato fruit ripening

Fruit ripening is a complex, genetically programmed process involving the action of critical transcription factors (TFs). Despite the established importance of WUSCHEL-related homeobox (WOX) TFs in plant development, the involvement of WOX and its underlying mechanism in the regulation of fruit ripening remain unclear. Here, we demonstrate that SlWOX13 regulates fruit ripening in tomato (Solanum lycopersicum). Overexpression of SlWOX13 accelerates fruit ripening, whereas loss-of-function mutation in SlWOX13 delays this process. Moreover, ethylene synthesis and carotenoid accumulation are significantly inhibited in slwox13 mutant fruit but accelerated in SlWOX13 transgenic fruit. Integrated analyses of RNA-seq and chromatin immunoprecipitation (ChIP)-seq identified 422 direct targets of SlWOX13, of which 243 genes are negatively regulated and 179 are positively regulated by SlWOX13. Electrophoretic mobility shift assay, RT-qPCR, dual-luciferase reporter assay, and ChIP-qPCR analyses demonstrated that SlWOX13 directly activates the expression of several genes involved in ethylene synthesis and signaling and carotenoid biosynthesis. Furthermore, SlWOX13 modulates tomato fruit ripening through key ripening-related TFs, such as RIPENING INHIBITOR (RIN), NON-RIPENING (NOR), and NAM, ATAF1, 2, and CUC2 4 (NAC4). Consequently, these effects promote fruit ripening. Taken together, these results demonstrate that SlWOX13 positively regulates tomato fruit ripening via both ethylene synthesis and signaling and by transcriptional regulation of key ripening-related TFs.

EARLY FLOWERING is a dominant gain-of-function allele of FANTASTIC FOUR 1/2c that promotes early flowering in tomato

Flowering time, an important factor in plant adaptability and genetic improvement, is regulated by various genes in tomato (Solanum lycopersicum). In this study, we characterized a tomato mutant, EARLY FLOWERING (EF), that developed flowers much earlier than its parental control. EF is a dominant gain-of-function allele with a T-DNA inserted 139 bp downstream of the stop codon of FANTASTIC FOUR 1/2c (FAF1/2c). The transcript of SlFAF1/2c was at elevated levels in the EF mutant. Overexpressing SlFAF1/2c in tomato plants phenocopied the early flowering trait of the EF mutant. Knocking out SlFAF1/2c in the EF mutant reverted the early flowering phenotype of the mutant to the normal flowering time of the wild-type tomato plants. SlFAF1/2c promoted the floral transition by shortening the vegetative phase rather than by reducing the number of leaves produced before the emergence of the first inflorescence. The COP9 signalosome subunit 5B (CSN5B) was shown to interact with FAF1/2c, and knocking out CSN5B led to an early flowering phenotype in tomato. Interestingly, FAF1/2c was found to reduce the accumulation of the CSN5B protein by reducing its protein stability. These findings imply that FAF1/2c regulates flowering time in tomato by reducing the accumulation and stability of CSN5B, which influences the expression of SINGLE FLOWER TRUSS (SFT), JOINTLESS (J) and UNIFLORA (UF). Thus, a new allele of SlFAF1/2c was discovered and found to regulate flowering time in tomato.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

Precise fine-turning of GhTFL1 by base editing tools defines ideal cotton plant architecture

BACKGROUND

CRISPR/Cas-derived base editor enables precise editing of target sites and has been widely used for basic research and crop genetic improvement. However, the editing efficiency of base editors at different targets varies greatly.

RESULTS

Here, we develop a set of highly efficient base editors in cotton plants. GhABE8e, which is fused to conventional nCas9, exhibits 99.9% editing efficiency, compared to GhABE7.10 with 64.9%, and no off-target editing is detected. We further replace nCas9 with dCpf1, which recognizes TTTV PAM sequences, to broaden the range of the target site. To explore the functional divergence of TERMINAL FLOWER 1 (TFL1), we edit the non-coding and coding regions of GhTFL1 with 26 targets to generate a comprehensive allelic population including 300 independent lines in cotton. This allows hidden pleiotropic roles for GhTFL1 to be revealed and allows us to rapidly achieve directed domestication of cotton and create ideotype germplasm with moderate height, shortened fruiting branches, compact plant, and early-flowering. Further, by exploring the molecular mechanism of the GhTFL1L86P and GhTFL1K53G+S78G mutations, we find that the GhTFL1L86P mutation weakens the binding strength of the GhTFL1 to other proteins but does not lead to a complete loss of GhTFL1 function.

CONCLUSIONS

This strategy provides an important technical platform and genetic information for the study and creation of ideal plant architecture.

Evolution of plant cell-type-specific cis -regulatory elements

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Application of genome editing techniques to regulate gene expression in crops

BACKGROUND

Enhanced agricultural production is urgently required to meet the food demands of the increasing global population. Abundant genetic diversity is expected to accelerate crop development. In particular, the development of the CRISPR/Cas genome editing technology has greatly enhanced our ability to improve crop's genetic diversity through direct artificial gene modification. However, recent studies have shown that most crop improvement efforts using CRISPR/Cas techniques have mainly focused on the coding regions, and there is a relatively lack of studies on the regulatory regions of gene expression.

RESULTS

This review briefly summarizes the development of CRISPR/Cas system in the beginning. Subsequently, the importance of gene regulatory regions in plants is discussed. The review focuses on recent developments and applications of mutations in regulatory regions via CRISPR/Cas techniques in crop breeding.

CONCLUSION

Finally, an outline of perspectives for future crop breeding using genome editing technologies is provided. This review provides new research insights for crop improvement using genome editing techniques.

PHB3 interacts with BRI1 and BAK1 to mediate brassinosteroid signal transduction in Arabidopsis and tomato

Brassinosteroids (BRs) are plant hormones that are essential in plant growth and development. BRASSINOSTEROID-INSENSITIVE 1 (BRI1) and BRI1 ASSOCIATED RECEPTOR KINASE 1 (BAK1), which are located on the plasma membrane, function as co-receptors that accept and transmit BR signals. PROHIBITIN 3 (PHB3) was identified in both BRI1 and BAK1 complexes by affinity purification and LC-MS/MS analysis. Biochemical data showed that BRI1/BAK1 interacted with PHB3 in vitro and in vivo. BRI1/BAK1 phosphorylated PHB3 in vitro. When the Thr-80 amino acid in PHB3 was mutated to Ala, the mutant protein was not phosphorylated by BRI1 and the mutant protein interaction with BRI1 was abolished in the yeast two-hybrid assay. BAK1 did not phosphorylate the mutant protein PHB3T54A . The loss-of-function phb3 mutant showed a weaker BR signal than the wild-type. Genetic analyses revealed that PHB3 is a BRI1/BAK1 downstream substrate that participates in BR signalling. PHB3 has five homozygous in tomato, and we named the closest to AtPHB3 as SlPHB3.1. Biochemical data showed that SlBRI1/SlSERK3A/SlSERK3B interacted with SlPHB3.1 and SlPHB3.3. The CRISPR-Cas9 method generated slphb3.1 mutant led to a BR signal stunted relatively in tomatoes. PHB3 is a new component of the BR signal pathway in both Arabidopsis and tomato.

New Advances in the Study of Regulation of Tomato Flowering-Related Genes Using Biotechnological Approaches

The tomato is a convenient object for studying reproductive processes, which has become a classic. Such complex processes as flowering and fruit setting require an understanding of the fundamental principles of molecular interaction, the structures of genes and proteins, the construction of signaling pathways for transcription regulation, including the synchronous actions of cis-regulatory elements (promoter and enhancer), trans-regulatory elements (transcription factors and regulatory RNAs), and transposable elements and epigenetic regulators (DNA methylation and acetylation, chromatin structure). Here, we discuss the current state of research on tomatoes (2017-2023) devoted to studying the function of genes that regulate flowering and signal regulation systems using genome-editing technologies, RNA interference gene silencing, and gene overexpression, including heterologous expression. Although the central candidate genes for these regulatory components have been identified, a complete picture of their relationship has yet to be formed. Therefore, this review summarizes the latest achievements related to studying the processes of flowering and fruit set. This work attempts to display the gene interaction scheme to better understand the events under consideration.

The WOX Genes from the Intermediate Clade: Influence on the Somatic Embryogenesis in Medicago truncatula

Transcription factors from the WOX family are well-known regulators of cell proliferation and differentiation in plants. Herein, we focused on several WOX genes from the intermediate clade and checked their impact on somatic embryogenesis using the model legume object Medicago truncatula. As a result, we show that MtWOX9-1 overexpression not only stimulates somatic embryogenesis in the embryogenic M. truncatula line, as it was shown previously, but can also induce somatic embryogenesis in the non-embryogenic line. Other intermediate clade WOX, including the close paralog of MtWOX9-1, as well as WOX11 homologs, did not have any significant impact on somatic embryogenesis in our in vitro cultivation system. Together, our results give new information about the diversity of the WOX family proteins and their specific functions. These data can be used for the search of new regeneration stimulators.

PlantPAN 4.0: updated database for identifying conserved non-coding sequences and exploring dynamic transcriptional regulation in plant promoters

PlantPAN 4.0 (http://PlantPAN.itps.ncku.edu.tw/) is an integrative resource for constructing transcriptional regulatory networks for diverse plant species. In this release, the gene annotation and promoter sequences were expanded to cover 115 species. PlantPAN 4.0 can help users characterize the evolutionary differences and similarities among cis-regulatory elements; furthermore, this system can now help in identification of conserved non-coding sequences among homologous genes. The updated transcription factor binding site repository contains 3428 nonredundant matrices for 18305 transcription factors; this expansion helps in exploration of combinational and nucleotide variants of cis-regulatory elements in conserved non-coding sequences. Additionally, the genomic landscapes of regulatory factors were manually updated, and ChIP-seq data sets derived from a single-cell green alga (Chlamydomonas reinhardtii) were added. Furthermore, the statistical review and graphical analysis components were improved to offer intelligible information through ChIP-seq data analysis. These improvements included easy-to-read experimental condition clusters, searchable gene-centered interfaces for the identification of promoter regions' binding preferences by considering experimental condition clusters and peak visualization for all regulatory factors, and the 20 most significantly enriched gene ontology functions for regulatory factors. Thus, PlantPAN 4.0 can effectively reconstruct gene regulatory networks and help compare genomic cis-regulatory elements across plant species and experiments.

Boosting Triticeae crop grain yield by manipulating molecular modules to regulate inflorescence architecture: insights and knowledge from other cereal crops

One of the challenges for global food security is to reliably and sustainably improve the grain yield of cereal crops. One solution is to modify the architecture of the grain-bearing inflorescence to optimize for grain number and size. Cereal inflorescences are complex structures, with determinacy, branching patterns, and spikelet/floret growth patterns that vary by species. Recent decades have witnessed rapid advancements in our understanding of the genetic regulation of inflorescence architecture in rice, maize, wheat, and barley. Here, we summarize current knowledge on key genetic factors underlying the different inflorescence morphologies of these crops and model plants (Arabidopsis and tomato), focusing particularly on the regulation of inflorescence meristem determinacy and spikelet meristem identity and determinacy. We also discuss strategies to identify and utilize these superior alleles to optimize inflorescence architecture and, ultimately, improve crop grain yield.

Cultivating potential: Harnessing plant stem cells for agricultural crop improvement

Meristems are stem cell-containing structures that produce all plant organs and are therefore important targets for crop improvement. Developmental regulators control the balance and rate of cell divisions within the meristem. Altering these regulators impacts meristem architecture and, as a consequence, plant form. In this review, we discuss genes involved in regulating the shoot apical meristem, inflorescence meristem, axillary meristem, root apical meristem, and vascular cambium in plants. We highlight several examples showing how crop breeders have manipulated developmental regulators to modify meristem growth and alter crop traits such as inflorescence size and branching patterns. Plant transformation techniques are another innovation related to plant meristem research because they make crop genome engineering possible. We discuss recent advances on plant transformation made possible by studying genes controlling meristem development. Finally, we conclude with discussions about how meristem research can contribute to crop improvement in the coming decades.

The STF/WOX1 MD is required for physical interaction with MtWOX9 and leaf blade outgrowth in Medicago truncatula

Plant-specific WUSCHEL-related homeobox (WOX) family transcription factors play critical roles in maintaining meristems and lateral organ development. The WUS clade member STF/LAM1 physically interacts with the intermediate clade member WOX9. This interaction contributes to their antagonistical functions on leaf blade outgrowth by competing for the same cis-elements in the promoter of their common target in M. truncatula and N. sylvestris. Here, we identified the main interaction domains of STF and MtWOX9 in Medicago, shedding light on the mechanism of WOX gene function. The middle domain of STF and MtWOX9 are both critical for the interaction, while the conserved motif of STF in the C-terminal domain is also required. Deletion of the middle domain of STF partially rescued the leaf blade phenotypes of the stf null mutant, indicating that the middle domain plays an essential role during leaf blade expansion. This finding provides a new insight that the versatility of WOX function is not only caused by the conserved DNA binding and repression domains but also by the middle domain that recruits different partners.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Extreme restructuring of cis -regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis -regulatory sequences. It remains unclear how such drastic cis -regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis- regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ∼125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 ( CLV3 ) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis -regulatory sequences. Our results demonstrate remarkable malleability in cis -regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis -regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis -regulation to effectively engineer trait variation from conserved productivity genes in crops.

Author summary

We investigated the evolution of cis -regulatory elements (CREs) and their interactions in the regulation of a plant stem cell regulator gene, CLAVATA3 (CLV3) , in Arabidopsis and tomato. Despite diverging ∼125 million years ago, the function and expression of CLV3 is conserved in these species; however, cis -regulatory sequences upstream and downstream have drastically diverged, preventing identification of conserved non-coding sequences between them. We used CRISPR-Cas9 to engineer dozens of mutations within the cis -regulatory regions of Arabidopsis and tomato CLV3. In tomato, our results show that tomato CLV3 function primarily relies on interactions among CREs in the 5' non-coding region, unlike Arabidopsis CLV3 , which depends on a more balanced distribution of functional CREs between the 5' and 3' regions. Therefore, despite a high degree of functional conservation, our study demonstrates divergent regulatory strategies between two distantly related CLV3 orthologs, with substantial alterations in regulatory sequences, their spatial arrangement, and their relative effects on CLV3 regulation. These results suggest that regulatory regions are not only extremely robust to mutagenesis, but also that the sequences underlying this robustness can be lineage-specific for conserved genes, due to the complex and often redundant interactions among CREs that ensure proper gene function amidst large-scale sequence turnover.

GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) homologs share conserved roles in growth repression

Crop engineering and de novo domestication using gene editing are new frontiers in agriculture. However, outside of well-studied crops and model systems, prioritizing engineering targets remains challenging. Evolution can guide us, revealing genes with deeply conserved roles that have repeatedly been selected in the evolution of plant form. Homologs of the transcription factor genes GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) have repeatedly been targets of selection in domestication and evolution, where they repress growth in many developmental contexts. This suggests a conserved role for these genes in regulating growth repression. To test this, we determined the roles of GT1 and VRS1 homologs in maize (Zea mays) and the distantly related grass brachypodium (Brachypodium distachyon) using gene editing and mutant analysis. In maize, gt1; vrs1-like1 (vrl1) mutants have derepressed growth of floral organs. In addition, gt1; vrl1 mutants bore more ears and more branches, indicating broad roles in growth repression. In brachypodium, Bdgt1; Bdvrl1 mutants have more branches, spikelets, and flowers than wild-type plants, indicating conserved roles for GT1 and VRS1 homologs in growth suppression over ca. 59 My of grass evolution. Importantly, many of these traits influence crop productivity. Notably, maize GT1 can suppress growth in arabidopsis (Arabidopsis thaliana) floral organs, despite ca. 160 My of evolution separating the grasses and arabidopsis. Thus, GT1 and VRS1 maintain their potency as growth regulators across vast timescales and in distinct developmental contexts. This work highlights the power of evolution to inform gene editing in crop improvement.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Genome editing for plant synthetic metabolic engineering and developmental regulation

Plant metabolism and development are a reflection of the orderly expression of genetic information intertwined with the environment interactions. Genome editing is the cornerstone for scientists to modify endogenous genes or introduce exogenous functional genes and metabolic pathways, holding immense potential applications in molecular breeding and biosynthesis. Over the course of nearly a decade of development, genome editing has advanced significantly beyond the simple cutting of double-stranded DNA, now enabling precise base and fragment replacements, regulation of gene expression and translation, as well as epigenetic modifications. However, the utilization of genome editing in plant synthetic metabolic engineering and developmental regulation remains exploratory. Here, we provide an introduction and a comprehensive overview of the editing attributes associated with various CRISPR/Cas tools, along with diverse strategies for the meticulous control of plant metabolic pathways and developments. Furthermore, we discuss the limitations of current approaches and future prospects for genome editing-driven plant breeding.

Tuning plant phenotypes by precise, graded downregulation of gene expression

The ability to control gene expression and generate quantitative phenotypic changes is essential for breeding new and desired traits into crops. Here we report an efficient, facile method for downregulating gene expression to predictable, desired levels by engineering upstream open reading frames (uORFs). We used base editing or prime editing to generate de novo uORFs or to extend existing uORFs by mutating their stop codons. By combining these approaches, we generated a suite of uORFs that incrementally downregulate the translation of primary open reading frames (pORFs) to 2.5-84.9% of the wild-type level. By editing the 5' untranslated region of OsDLT, which encodes a member of the GRAS family and is involved in the brassinosteroid transduction pathway, we obtained, as predicted, a series of rice plants with varied plant heights and tiller numbers. These methods offer an efficient way to obtain genome-edited plants with graded expression of traits.

Engineering homoeologs provide a fine scale for quantitative traits in polyploid

Numerous staple crops exhibit polyploidy and are difficult to genetically modify. However, recent advances in genome sequencing and editing have enabled polyploid genome engineering. The hexaploid black nightshade species Solanum nigrum has immense potential as a beneficial food supplement. We assembled its genome at the scaffold level. After functional annotations, we identified homoeologous gene sets, with similar sequence and expression profiles, based on comparative analyses of orthologous genes with close diploid relatives Solanum americanum and S. lycopersicum. Using CRISPR-Cas9-mediated mutagenesis, we generated various mutation combinations in homoeologous genes. Multiple mutants showed quantitative phenotypic changes based on the genotype, resulting in a broad-spectrum effect on the quantitative traits of hexaploid S. nigrum. Furthermore, we successfully improved the fruit productivity of Boranong, an orphan cultivar of S. nigrum suggesting that engineering homoeologous genes could be useful for agricultural improvement of polyploid crops.

Sea-ATI unravels novel vocabularies of plant active cistrome

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants.

Targeting editing of tomato SPEECHLESS cis-regulatory regions generates plants with altered stomatal density in response to changing climate conditions

Flexible developmental programs enable plants to customize their organ size and cellular composition. In leaves of eudicots, the stomatal lineage produces two essential cell types, stomata and pavement cells, but the total numbers and ratio of these cell types can vary. Central to this flexibility is the stomatal lineage initiating transcription factor, SPEECHLESS (SPCH). Here we show, by multiplex CRISPR/Cas9 editing of SlSPCH cis-regulatory sequences in tomato, that we can identify variants with altered stomatal development responses to light and temperature cues. Analysis of tomato leaf development across different conditions, aided by newly-created tools for live-cell imaging and translational reporters of SlSPCH and its paralogues SlMUTE and SlFAMA, revealed the series of cellular events that lead to the environmental change-driven responses in leaf form. Plants bearing the novel SlSPCH variants generated in this study are powerful resources for fundamental and applied studies of tomato resilience in response to climate change.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

Gene editing of non-coding regulatory DNA and its application in crop improvement

The development of the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein (CRISPR/Cas) system has provided precise and efficient strategies to edit target genes and generate transgene-free crops. Significant progress has been made in the editing of protein-coding genes; however, studies on the editing of non-coding DNA with regulatory roles lags far behind. Non-coding regulatory DNAs, including those which can be transcribed into long non-coding RNAs (lncRNAs), and miRNAs, together with cis-regulatory elements (CREs), play crucial roles in regulating plant growth and development. Therefore, the combination of CRISPR/Cas technology and non-coding regulatory DNA has great potential to generate novel alleles that affect various agronomic traits of crops, thus providing valuable genetic resources for crop breeding. Herein, we review recent advances in the roles of non-coding regulatory DNA, attempts to edit non-coding regulatory DNA for crop improvement, and potential application of novel editing tools in modulating non-coding regulatory DNA. Finally, the existing problems, possible solutions, and future applications of gene editing of non-coding regulatory DNA in modern crop breeding practice are also discussed.

Characterization and expression profiles of WUSCHEL-related homeobox (WOX) gene family in cultivated alfalfa (Medicago sativa L.)

The WUSCHEL-related homeobox (WOX) family members are plant-specific transcriptional factors, which function in meristem maintenance, embryogenesis, lateral organ development, as well as abiotic stress tolerance. In this study, 14 MsWOX transcription factors were identified and comprehensively analyzed in the cultivated alfalfa cv. Zhongmu No.1. Overall, 14 putative MsWOX members containing conserved structural regions were clustered into three clades according to phylogenetic analysis. Specific expression patterns of MsWOXs in different tissues at different levels indicated that the MsWOX genes play various roles in alfalfa. MsWUS, MsWOX3, MsWOX9, and MsWOX13-1 from the three subclades were localized in the nucleus, among which, MsWUS and MsWOX13-1 exhibited strong self-activations in yeast. In addition, various cis-acting elements related to hormone responses, plant growth, and stress responses were identified in the 3.0 kb promoter regions of MsWOXs. Expression detection of separated shoots and roots under hormones including auxin, cytokinin, GA, and ABA, as well as drought and cold stresses, showed that MsWOX genes respond to different hormones and abiotic stress treatments. Furthermore, transcript abundance of MsWOX3, and MsWOX13-2 were significantly increased after rhizobia inoculation. This study presented comprehensive data on MsWOX transcription factors and provided valuable insights into further studies of their roles in developmental processes and abiotic stress responses in alfalfa.

Frontiers and techniques in plant gene regulation

Understanding plant gene regulation has been a priority for generations of plant scientists. However, due to its complex nature, the regulatory code governing plant gene expression has yet to be deciphered comprehensively. Recently developed methods-often relying on next-generation sequencing technology and state-of-the-art computational approaches-have started to further our understanding of the gene regulatory logic used by plants. In this review, we discuss these methods and the insights into the regulatory code of plants that they can yield.

CRISPR/Cas-mediated plant genome editing: outstanding challenges a decade after implementation

The discovery of the CRISPR/Cas genome-editing system has revolutionized our understanding of the plant genome. CRISPR/Cas has been used for over a decade to modify plant genomes for the study of specific genes and biosynthetic pathways as well as to speed up breeding in many plant species, including both model and non-model crops. Although the CRISPR/Cas system is very efficient for genome editing, many bottlenecks and challenges slow down further improvement and applications. In this review we discuss the challenges that can occur during tissue culture, transformation, regeneration, and mutant detection. We also review the opportunities provided by new CRISPR platforms and specific applications related to gene regulation, abiotic and biotic stress response improvement, and de novo domestication of plants.

The CCT transcriptional activator Ghd2 constantly delays the heading date by upregulating CO3 in rice

CONSTANS, CO-like, and TOC1 (CCT) family genes play important roles in regulating heading date, which exerts a large impact on the regional and seasonal adaptation of rice. Previous studies have shown that Grain number, plant height, and heading date2 (Ghd2) exhibits a negative response to drought stress by directly upregulating Rubisco activase and exerting a negative effect on heading date. However, the target gene of Ghd2 regulating heading date is still unknown. In this study, CO3 is identified by analyzing Ghd2 ChIP-seq data. Ghd2 activates CO3 expression by binding to the CO3 promoter through its CCT domain. EMSA experiments show that the motif CCACTA in the CO3 promoter was recognized by Ghd2. A comparison of the heading dates among plants with CO3 knocked out or overexpressed and double-mutants with Ghd2 overexpressed and CO3 knocked out shows that CO3 negatively and constantly regulates flowering by repressing the transcription of Ehd1, Hd3a, and RFT1. In addition, the target genes of CO3 are explored via a comprehensive analysis of DAP-seq and RNA-seq data. Taken together, these results suggest that Ghd2 directly binds to the downstream gene CO3, and the Ghd2-CO3 module constantly delays heading date via the Ehd1-mediated pathway.

Advances in cis-element- and natural variation-mediated transcriptional regulation and applications in gene editing of major crops

Transcriptional regulation is crucial to control of gene expression. Both spatio-temporal expression patterns and expression levels of genes are determined by the interaction between cis-acting elements and trans-acting factors. Numerous studies have focused on the trans-acting factors that mediate transcriptional regulatory networks. However, cis-acting elements, such as enhancers, silencers, transposons, and natural variations in the genome, are also vital for gene expression regulation and could be utilized by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated gene editing to improve crop quality and yield. In this review, we discuss current understanding of cis-element-mediated transcriptional regulation in major crops, including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays), as well as the latest advancements in gene editing techniques and their applications in crops to highlight prospective strategies for crop breeding.

CRISPR/Cas9 and Agrobacterium tumefaciens virulence proteins synergistically increase efficiency of precise genome editing via homology directed repair in plants

CRISPR/Cas9 genome editing and Agrobacterium tumefaciens-mediated genetic transformation are widely-used plant biotechnology tools derived from bacterial immunity-related systems, each involving DNA modification. The Cas9 endonuclease introduces DNA double-strand breaks (DSBs), and the A. tumefaciens T-DNA is released by the VirD2 endonuclease assisted by VirDl and attached by VirE2, transferred to the plant nucleus and integrated into the genome. Here, we explored the potential for synergy between the two systems and found that Cas9 and three virulence (Vir) proteins achieve precise genome editing via the homology directed repair (HDR) pathway in tobacco and rice plants. Compared with Cas9T (Cas9, VirD1, VirE2) and CvD (Cas9-VirD2) systems, the HDR frequencies of a foreign GFPm gene in the CvDT system (Cas9-VirD2, VirD1, VirE2) increased 52-fold and 22-fold, respectively. Further optimization of the CvDT process with a donor linker (CvDTL) achieved a remarkable increase in the efficiency of HDR-mediated genome editing. Additionally, the HDR efficiency of the three rice endogenous genes ACETOLACTATE SYNTHASE (ALS), PHYTOENE DESATURASE (PDS), and NITROGEN TRANSPORTER 1.1 B (NRT1.1B) increased 24-, 32- and 16-fold, respectively, in the CvDTL system, compared with corresponding Cas9TL (Cas9T process with a donor linker). Our results suggest that collaboration between CRISPR/Cas9 and Agrobacterium-mediated genetic transformation can make great progress towards highly efficient and precise genome editing via the HDR pathway.

Fertilization controls tiller numbers via transcriptional regulation of a MAX1-like gene in rice cultivation

Fertilization controls various aspects of cereal growth such as tiller number, leaf size, and panicle size. However, despite such benefits, global chemical fertilizer use must be reduced to achieve sustainable agriculture. Here, based on field transcriptome data from leaf samples collected during rice cultivation, we identify fertilizer responsive genes and focus on Os1900, a gene orthologous to Arabidopsis thaliana MAX1, which is involved in strigolactone biosynthesis. Elaborate genetic and biochemical analyses using CRISPR/Cas9 mutants reveal that Os1900 together with another MAX1-like gene, Os5100, play a critical role in controlling the conversion of carlactone into carlactonoic acid during strigolactone biosynthesis and tillering in rice. Detailed analyses of a series of Os1900 promoter deletion mutations suggest that fertilization controls tiller number in rice through transcriptional regulation of Os1900, and that a few promoter mutations alone can increase tiller numbers and grain yields even under minor-fertilizer conditions, whereas a single defective os1900 mutation does not increase tillers under normal fertilizer condition. Such Os1900 promoter mutations have potential uses in breeding programs for sustainable rice production.

Antagonistic regulation of target genes by the SISTER OF TM3-JOINTLESS2 complex in tomato inflorescence branching

Inflorescence branch number is a yield-related trait controlled by cell fate determination in meristems. Two MADS-box transcription factors (TFs)-SISTER OF TM3 (STM3) and JOINTLESS 2 (J2)-have opposing regulatory roles in inflorescence branching. However, the mechanisms underlying their regulatory functions in inflorescence determinacy remain unclear. Here, we characterized the functions of these TFs in tomato (Solanum lycopersicum) floral meristem and inflorescence meristem (IM) through chromatin immunoprecipitation and sequencing analysis of their genome-wide occupancy. STM3 and J2 activate or repress the transcription of a set of common putative target genes, respectively, through recognition and binding to CArG box motifs. FRUITFULL1 (FUL1) is a shared putative target of STM3 and J2 and these TFs antagonistically regulate FUL1 in inflorescence branching. Moreover, STM3 physically interacts with J2 to mediate its cytosolic redistribution and restricts J2 repressor activity by reducing its binding to target genes. Conversely, J2 limits STM3 regulation of target genes by transcriptional repression of the STM3 promoter and reducing STM3-binding activity. Our study thus reveals an antagonistic regulatory relationship in which STM3 and J2 control tomato IM determinacy and branch number.

Long-distance control of potato storage organ formation by SELF PRUNING 3D and FLOWERING LOCUS T-like 1

Plants program their meristem-associated developmental switches for timely adaptation to a changing environment. Potato (Solanum tuberosum L.) tubers differentiate from specialized belowground branches or stolons through radial expansion of their terminal ends. During this process, the stolon apex and closest axillary buds enter a dormancy state that leads to tuber eyes, which are reactivated the following spring and generate a clonally identical plant. The potato FLOWERING LOCUS T homolog SELF-PRUNING 6A (StSP6A) was previously identified as the major tuber-inducing signal that integrates day-length cues to control the storage switch. However, whether some other long-range signals also act as tuber organogenesis stimuli remains unknown. Here, we show that the florigen SELF PRUNING 3D (StSP3D) and FLOWERING LOCUS T-like 1 (StFTL1) genes are activated by short days, analogously to StSP6A. Overexpression of StSP3D or StFTL1 promotes tuber formation under non-inductive long days, and the tuber-inducing activity of these proteins is graft transmissible. Using the non-tuber-bearing wild species Solanum etuberosum, a natural SP6A null mutant, we show that leaf-expressed SP6A is dispensable for StSP3D long-range activity. StSP3D and StFTL1 mediate secondary activation of StSP6A in stolon tips, leading to amplification of this tuberigen signal. StSP3D and StFTL1 were observed to bind the same protein partners as StSP6A, suggesting that they can also form transcriptionally active complexes. Together, our findings show that additional mobile tuber-inducing signals are regulated by the photoperiodic pathway.

Identification of clade-wide putative cis-regulatory elements from conserved non-coding sequences in Cucurbitaceae genomes

Cis-regulatory elements regulate gene expression and play an essential role in the development and physiology of organisms. Many conserved non-coding sequences (CNSs) function as cis-regulatory elements. They control the development of various lineages. However, predicting clade-wide cis-regulatory elements across several closely related species remains challenging. Based on the relationship between CNSs and cis-regulatory elements, we present a computational approach that predicts the clade-wide putative cis-regulatory elements in 12 Cucurbitaceae genomes. Using 12-way whole-genome alignment, we first obtained 632 112 CNSs in Cucurbitaceae. Next, we identified 16 552 Cucurbitaceae-wide cis-regulatory elements based on collinearity among all 12 Cucurbitaceae plants. Furthermore, we predicted 3 271 potential regulatory pairs in the cucumber genome, of which 98 were verified using integrative RNA sequencing and ChIP sequencing datasets from samples collected during various fruit development stages. The CNSs, Cucurbitaceae-wide cis-regulatory elements, and their target genes are accessible at http://cmb.bnu.edu.cn/cisRCNEs_cucurbit/. These elements are valuable resources for functionally annotating CNSs and their regulatory roles in Cucurbitaceae genomes.

Combined genome-wide association study and epistasis analysis reveal multifaceted genetic architectures of plant height in Asian cultivated rice

Plant height (PH) in rice (Oryza sativa) is an important trait for its adaptation and agricultural performance. Discovery of the semi-dwarf1 (SD1) mutation initiated the Green Revolution, boosting rice yield and fitness, but the underlying genetic regulation of PH in rice remains largely unknown. Here, we performed genome-wide association study (GWAS) and identified 12 non-repetitive QTL/genes regulating PH variation in 619 Asian cultivated rice accessions. One of these was an SD1 structural variant, not normally detected in standard GWAS analyses. Given the strong effect of SD1 on PH, we also divided 619 accessions into subgroups harbouring distinct SD1 haplotypes, and found a further 85 QTL/genes for PH, revealing genetic heterogeneity that may be missed by analysing a broad, diverse population. Moreover, we uncovered two epistatic interaction networks of PH-associated QTL/genes in the japonica (Geng)-dominant SD1^NIP subgroup. In one of them, the hub QTL/gene qphSN1.4/GAMYB interacted with qphSN3.1/OsINO80, qphSN3.4/HD16/EL1, qphSN6.2/LOC_Os06g11130, and qphSN10.2/MADS56. Sequence variations in GAMYB and MADS56 were associated with their expression levels and PH variations, and MADS56 was shown to physically interact with MADS57 to coregulate expression of gibberellin (GA) metabolic genes OsGA2ox3 and Elongated Uppermost Internode1 (EUI1). Our study uncovered the multifaceted genetic architectures of rice PH, and provided novel and abundant genetic resources for breeding semi-dwarf rice and new candidates for further mechanistic studies on regulation of PH in rice.

An efficient CRISPR-Cas12a promoter editing system for crop improvement

Promoter editing represents an innovative approach to introduce quantitative trait variation (QTV) in crops. However, an efficient promoter editing system for QTV needs to be established. Here we develop a CRISPR-Cas12a promoter editing (CAPE) system that combines a promoter key-region estimating model and an efficient CRISPR-Cas12a-based multiplexed or singular editing system. CAPE is benchmarked in rice to produce QTV continuums for grain starch content and size by targeting OsGBSS1 and OsGS3, respectively. We then apply CAPE for promoter editing of OsD18, a gene encoding GA3ox in the gibberellin biosynthesis pathway. The resulting lines carry a QTV continuum of semidwarfism without significantly compromising grain measures. Field trials demonstrated that the OsD18 promoter editing lines have the same yield performance and antilodging phenotype as the Green Revolution OsSD1 mutants in different genetic backgrounds. Hence, promoter editing of OsD18 generates a quantitative Green Revolution trait. Together, we demonstrate a CAPE-based promoter editing and tuning pipeline for efficient production of useful QTV continuum in crops.

Can the Wild Perennial, Rhizomatous Rice Species Oryza longistaminata be a Candidate for De Novo Domestication?

As climate change intensifies, the development of resilient rice that can tolerate abiotic stresses is urgently needed. In nature, many wild plants have evolved a variety of mechanisms to protect themselves from environmental stresses. Wild relatives of rice may have abundant and virtually untapped genetic diversity and are an essential source of germplasm for the improvement of abiotic stress tolerance in cultivated rice. Unfortunately, the barriers of traditional breeding approaches, such as backcrossing and transgenesis, make it challenging and complex to transfer the underlying resilience traits between plants. However, de novo domestication via genome editing is a quick approach to produce rice with high yields from orphans or wild relatives. African wild rice, Oryza longistaminata, which is part of the AA-genome Oryza species has two types of propagation strategies viz. vegetative propagation via rhizome and seed propagation. It also shows tolerance to multiple types of abiotic stress, and therefore O. longistaminata is considered a key candidate of wild rice for heat, drought, and salinity tolerance, and it is also resistant to lodging. Importantly, O. longistaminata is perennial and propagates also via rhizomes both of which are traits that are highly valuable for the sustainable production of rice. Therefore, O. longistaminata may be a good candidate for de novo domestication through genome editing to obtain rice that is more climate resilient than modern elite cultivars of O. sativa.