Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

Conserved functions of chromatin regulators in basal Archaeplastida

Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb-group genes which mediate H3K27me3-based silencing and play a role in balancing gene dosage and regulating haploid-to-diploid transitions by tissue-specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein-coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida.

Novel synthetic inducible promoters controlling gene expression during water-deficit stress with green tissue specificity in transgenic poplar

Synthetic promoters may be designed using short cis-regulatory elements (CREs) and core promoter sequences for specific purposes. We identified novel conserved DNA motifs from the promoter sequences of leaf palisade and vascular cell type-specific expressed genes in water-deficit stressed poplar (Populus tremula × Populus alba), collected through low-input RNA-seq analysis using laser capture microdissection. Hexamerized sequences of four conserved 20-base motifs were inserted into each synthetic promoter construct. Two of these synthetic promoters (Syn2 and Syn3) induced GFP in transformed poplar mesophyll protoplasts incubated in 0.5 M mannitol solution. To identify effect of length and sequence from a valuable 20 base motif, 5' and 3' regions from a basic sequence (GTTAACTTCAGGGCCTGTGG) of Syn3 were hexamerized to generate two shorter synthetic promoters, Syn3-10b-1 (5': GTTAACTTCA) and Syn3-10b-2 (3': GGGCCTGTGG). These promoters' activities were compared with Syn3 in plants. Syn3 and Syn3-10b-1 were specifically induced in transient agroinfiltrated Nicotiana benthamiana leaves in water cessation for 3 days. In stable transgenic poplar, Syn3 presented as a constitutive promoter but had the highest activity in leaves. Syn3-10b-1 had stronger induction in green tissues under water-deficit stress conditions than mock control. Therefore, a synthetic promoter containing the 5' sequence of Syn3 endowed both tissue-specificity and water-deficit inducibility in transgenic poplar, whereas the 3' sequence did not. Consequently, we have added two new synthetic promoters to the poplar engineering toolkit: Syn3-10b-1, a green tissue-specific and water-deficit stress-induced promoter, and Syn3, a green tissue-preferential constitutive promoter.

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Paramutation at the maize pl1 locus is associated with RdDM activity at distal tandem repeats

Exceptions to Mendelian inheritance often highlight novel chromosomal behaviors. The maize Pl1-Rhoades allele conferring plant pigmentation can display inheritance patterns deviating from Mendelian expectations in a behavior known as paramutation. However, the chromosome features mediating such exceptions remain unknown. Here we show that small RNA production reflecting RNA polymerase IV function within a distal downstream set of five tandem repeats is coincident with meiotically-heritable repression of the Pl1-Rhoades transcription unit. A related pl1 haplotype with three, but not one with two, repeat units also displays the trans-homolog silencing typifying paramutations. 4C interactions, CHD3a-dependent small RNA profiles, nuclease sensitivity, and polyadenylated RNA levels highlight a repeat subregion having regulatory potential. Our comparative and mutant analyses show that transcriptional repression of Pl1-Rhoades correlates with 24-nucleotide RNA production and cytosine methylation at this subregion indicating the action of a specific DNA-dependent RNA polymerase complex. These findings support a working model in which pl1 paramutation depends on trans-chromosomal RNA-directed DNA methylation operating at a discrete cis-linked and copy-number-dependent transcriptional regulatory element.

Genome-wide characterization and functional analysis of the melon TGA gene family in disease resistance through ectopic overexpression in Arabidopsis thaliana

TGA-binding (TGA) transcription factors, characterized by the basic region/leucine zipper motif (bZIP), have been recognized as pivotal regulators in plant growth, development, and stress responses through their binding to the as-1 element. In this study, the TGA gene families in melon, watermelon, cucumber, pumpkin, and zucchini were comprehensively characterized, encompassing analyses of gene/protein structures, phylogenetic relationships, gene duplication events, and cis-acting elements in gene promoters. Upon transient expression in Nicotiana benthamiana, the melon CmTGAs, with typical bZIP and DOG1 domains, were observed to localize within the nucleus. Biochemical investigation revealed specific interactions between CmTGA2/3/5/8/9 and CmNPR3 or CmNPR4. The CmTGA genes exhibited differential expression patterns in melon plants in response to different hormones like salicylic acid, methyl jasmonate, and ethylene, as well as a fungal pathogen, Stagonosporopsis cucurbitacearum that causes gummy stem blight in melon. The overexpression of CmTGA3, CmTGA8, and CmTGA9 in Arabidopsis plants resulted in the upregulation of AtPR1 and AtPR5 expression, thereby imparting enhanced resistance to Pseudomonas syringae pv. Tomato DC3000. In contrast, the overexpression of CmTGA7 or CmTGA9 resulted in a compromised resistance to Botrytis cinerea, coinciding with a concomitant reduction in the expression levels of AtPDF1.2 and AtMYC2 following infection with B. cinerea. These findings shed light on the important roles of specific CmTGA genes in plant immunity, suggesting that genetic manipulation of these genes could be a promising avenue for enhancing plant immune responses.

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.

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.

The Arabidopsis Mitochondrial Nucleoid-Associated Protein WHIRLY2 Is Required for a Proper Response to Salt Stress

In the last years, plant organelles have emerged as central coordinators of responses to internal and external stimuli, which can induce stress. Mitochondria play a fundamental role as stress sensors being part of a complex communication network between the organelles and the nucleus. Among the different environmental stresses, salt stress poses a significant challenge and requires efficient signaling and protective mechanisms. By using the why2 T-DNA insertion mutant and a novel knock-out mutant prepared by CRISPR/Cas9-mediated genome editing, this study revealed that WHIRLY2 is crucial for protecting mitochondrial DNA (mtDNA) integrity during salt stress. Loss-of-function mutants show an enhanced sensitivity to salt stress. The disruption of WHIRLY2 causes the impairment of mtDNA repair that results in the accumulation of aberrant recombination products, coinciding with severe alterations in nucleoid integrity and overall mitochondria morphology besides a compromised redox-dependent response and misregulation of antioxidant enzymes. The results of this study revealed that WHIRLY2-mediated structural features in mitochondria (nucleoid compactness and cristae) are important for an effective response to salt stress.

Identification of a Unique Genomic Region in Sweet Chestnut (Castanea sativa Mill.) That Controls Resistance to Asian Chestnut Gall Wasp Dryocosmus kuriphilus Yasumatsu

The Asian chestnut gall wasp (ACGW) (Hymenoptera Dryocosmus kuriphilus Yasumatsu) is a severe pest of sweet chestnut (Castanea sativa Mill.) with a strong impact on growth and nut production. A comparative field trial in Central Italy, including provenances from Spain, Italy, and Greece, was screened for ACGW infestation over consecutive years. The Greek provenance Hortiatis expressed a high proportion of immune plants and was used to perform a genome-wide association study based on DNA pool sequencing (Pool-GWAS) by comparing two DNA pools from 25 susceptible and 25 resistant plants. DNA pools were sequenced with 50X coverage depth. Sequence reads were aligned to a C. mollissima reference genome and the pools were compared to identify SNPs associated with resistance. Twenty-one significant SNPs were identified and highlighted a small genomic region on pseudochromosome 3 (Chr 3), containing 12 candidate genes of three gene families: Cytochrome P450, UDP-glycosyltransferase, and Rac-like GTP-binding protein. Functional analyses revealed a putative metabolic gene cluster related to saccharide biosynthesis in the genomic regions associated with resistance that could be involved in the production of a toxic metabolite against parasites. The comparison with previous genetic studies confirmed the involvement of Chr 3 in the control of resistance to ACGW.

Genome-wide profiling of angiogenic cis-regulatory elements unravels cis-regulatory SNPs for vascular abnormality

Angiogenesis is extensively involved in embryonic development and requires complex regulation networks, whose defects can cause a variety of vascular abnormalities. Cis-regulatory elements control gene expression at all developmental stages, but they have not been studied or profiled in angiogenesis yet. In this study, we exploited public DNase-seq and RNA-seq datasets from a VEGFA-stimulated in vitro angiogenic model, and carried out an integrated analysis of the transcriptome and chromatin accessibility across the entire process. Totally, we generated a bank of 47,125 angiogenic cis-regulatory elements with promoter (marker by H3K4me3) and/or enhancer (marker by H3K27ac) activities. Motif enrichment analysis revealed that these angiogenic cis-regulatory elements interacted preferentially with ETS family TFs. With this tool, we performed an association study using our WES data of TAPVC and identified rs199530718 as a cis-regulatory SNP associated with disease risk. Altogether, this study generated a genome-wide bank of angiogenic cis-regulatory elements and illustrated its utility in identifying novel cis-regulatory SNPs for TAPVC, expanding new horizons of angiogenesis as well as vascular abnormality genetics.

Characterization of activating cis-regulatory elements from the histone genes of Chlamydomonas reinhardtii

Regulation of gene expression in eukaryotes is controlled by cis-regulatory modules (CRMs). A major class of CRMs are enhancers which are composed of activating cis-regulatory elements (CREs) responsible for upregulating transcription. To date, most enhancers and activating CREs have been studied in angiosperms; in contrast, our knowledge about these key regulators of gene expression in green algae is limited. In this study, we aimed at characterizing putative activating CREs/CRMs from the histone genes of the unicellular model alga Chlamydomonas reinhardtii. To test the activity of four candidates, reporter constructs consisting of a tetramerized CRE, an established promoter, and a gene for the mCerulean3 fluorescent protein were incorporated into the nuclear genome of C. reinhardtii, and their activity was quantified by flow cytometry. Two tested candidates, Eupstr and Ehist cons, significantly upregulated gene expression and were characterized in detail. Eupstr, which originates from highly expressed genes of C. reinhardtii, is an orientation-independent CRE capable of activating both the RBCS2 and β2-tubulin promoters. Ehist cons, which is a CRM from histone genes of angiosperms, upregulates the β2-tubulin promoter in C. reinhardtii over a distance of at least 1.5 kb. The octamer motif present in Ehist cons was identified in C. reinhardtii and the related green algae Chlamydomonas incerta, Chlamydomonas schloesseri, and Edaphochlamys debaryana, demonstrating its high evolutionary conservation. The results of this investigation expand our knowledge about the regulation of gene expression in green algae. Furthermore, the characterized activating CREs/CRMs can be applied as valuable genetic tools.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Predicting gene expression responses to environment in Arabidopsis thaliana using natural variation in DNA sequence

The evolution of gene expression responses are a critical component of adaptation to variable environments. Predicting how DNA sequence influences expression is challenging because the genotype to phenotype map is not well resolved for cis regulatory elements, transcription factor binding, regulatory interactions, and epigenetic features, not to mention how these factors respond to environment. We tested if flexible machine learning models could learn some of the underlying cis- regulatory genotype to phenotype map. We tested this approach using cold-responsive transcriptome profiles in 5 diverse Arabidopsis thaliana accessions. We first tested for evidence that cis regulation plays a role in environmental response, finding 14 and 15 motifs that were significantly enriched within the up- and down-stream regions of cold-responsive differentially regulated genes (DEGs). We next applied convolutional neural networks (CNNs), which learn de novo cis- regulatory motifs in DNA sequences to predict expression response to environment. We found that CNNs predicted differential expression with moderate accuracy, with evidence that predictions were hindered by biological complexity of regulation and the large potential regulatory code. Overall, DEGs between specific environments can be predicted based on variation in cis- regulatory sequences, although more information needs to be incorporated and better models may be required.

Unraveling the role of PlARF2 in regulating deed formancy in Paeonia lactiflora

MAIN CONCLUSION

PlARF2 can positively regulate the seed dormancy in Paeonia lactiflora Pall. and bind the RY cis-element. Auxin, a significant phytohormone influencing seed dormancy, has been demonstrated to be regulated by auxin response factors (ARFs), key transcriptional modulators in the auxin signaling pathway. However, the role of this class of transcription factors (TFs) in perennials with complex seed dormancy mechanisms remains largely unexplored. Here, we cloned and characterized an ARF gene from Paeonia lactiflora, named PlARF2, which exhibited differential expression levels in the seeds during the process of seed dormancy release. The deduced amino acid sequence of PlARF2 had high homology with those of other plants and contained typical conserved Auxin_resp domain of the ARF family. Phylogenetic analysis revealed that PlARF2 was closely related to VvARF3 in Vitis vinifera. The subcellular localization and transcriptional activation assay showed that PlARF2 is a nuclear protein possessing transcriptional activation activity. The expression levels of dormancy-related genes in transgenic callus indicated that PlARF2 was positively correlated with the contents of PlABI3 and PlDOG1. The germination assay showed that PlARF2 promoted seed dormancy. Moreover, TF Centered Yeast one-hybrid assay (TF-Centered Y1H), electrophoretic mobility shift assay (EMSA) and dual-luciferase reporter assay analysis (Dual-Luciferase) provided evidence that PlARF2 can bind to the 'CATGCATG' motif. Collectively, our findings suggest that PlARF2, as TF, could be involved in the regulation of seed dormancy and may act as a repressor of germination.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

Multiome in the Same Cell Reveals the Impact of Osmotic Stress on Arabidopsis Root Tip Development at Single-Cell Level

Cell-specific transcriptional regulatory networks (TRNs) play vital roles in plant development and response to environmental stresses. However, traditional single-cell mono-omics techniques are unable to directly capture the relationships and dynamics between different layers of molecular information within the same cells. While advanced algorithm facilitates merging scRNA-seq and scATAC-seq datasets, accurate data integration remains a challenge, particularly when investigating cell-type-specific TRNs. By examining gene expression and chromatin accessibility simultaneously in 16,670 Arabidopsis root tip nuclei, the TRNs are reconstructed that govern root tip development under osmotic stress. In contrast to commonly used computational integration at cell-type level, 12,968 peak-to-gene linkage is captured at the bona fide single-cell level and construct TRNs at an unprecedented resolution. Furthermore, the unprecedented datasets allow to more accurately reconstruct the coordinated changes of gene expression and chromatin states during cellular state transition. During root tip development, chromatin accessibility of initial cells precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for subsequent differentiation steps. Pseudo-time trajectory analysis reveal that osmotic stress can shift the functional differentiation of trichoblast. Candidate stress-related gene-linked cis-regulatory elements (gl-cCREs) as well as potential target genes are also identified, and uncovered large cellular heterogeneity under osmotic stress.

Increasing deletion sizes and the efficiency of CRISPR/Cas9-mediated mutagenesis by SunTag-mediated TREX1 recruitment

Previously, it has been shown that mutagenesis frequencies can be improved by directly fusing the human exonuclease TREX2 to Cas9, resulting in a strong increase in the frequency of smaller deletions at the cut site. Here, we demonstrate that, by using the SunTag system for recruitment of TREX2, the mutagenesis efficiency can be doubled in comparison to the direct fusion in Arabidopsis thaliana. Therefore, we also tested the efficiency of the system for targeted deletion formation by recruiting two other 3'-5' exonucleases, namely the human TREX1 and E. coli ExoI. It turns out that SunTag-mediated recruitment of TREX1 not only improved the general mutation induction efficiency slightly in comparison to TREX2, but that, more importantly, the mean size of the induced deletions was also enhanced, mainly via an increase of deletions of 25 bp or more. EcExoI also yielded a higher amount of larger deletions. However, only in the case of TREX1 and TREX2, the effect was predominately SunTag-dependent, indicating efficient target-specific recruitment. Using SunTag-mediated TREX1 recruitment at other genomic sites, we were able to obtain similar deletion patterns. Thus, we were able to develop an attractive novel editing tool that is especially useful for obtaining deletions in the range from 20 to 40 bp around the cut site. Such sizes are often required for the manipulation of cis-regulatory elements. This feature is closing an existing gap as previous approaches, based on single nucleases or paired nickases or nucleases, resulted in either shorter or longer deletions, respectively.

Heritable, multinucleotide deletions in plants using viral delivery of a repair exonuclease and guide RNAs

CRISPR/Cas9-mediated mutagenesis typically results in short insertion/deletion mutations, which are often too small to disrupt the function of cis-acting regulatory elements. Here, we describe a highly efficient in planta gene editing approach called VirTREX2-HLDel that achieves heritable multinucleotide deletions in both protein-coding genes and noncoding DNA regulatory elements. VirTREX2-HLDel uses RNA viruses to deliver both the 3 prime repair exonuclease 2 (TREX2) and single-guide RNAs. Our method enables recovery of multiplexed heritable deletions and increases the heritable gene editing frequency at poorly edited sites. We identified functional conservation and divergence of MICRORNA164 (miR164) in Nicotiana benthamiana and tomato (Solanum lycopersicum) using VirTREX2-HLDel and observed previously uncharacterized phenotypes in plants with large deletions at this locus. Our viral delivery method reduces the need for tissue culture and will accelerate the understanding of protein-coding and regulatory regions in plants.

Genetic dissection of cis-regulatory control of ZmWUSCHEL1 expression by type B RESPONSE REGULATORS

Mutations in cis-regulatory regions play an important role in the domestication and improvement of crops by altering gene expression. However, assessing the in vivo impact of cis-regulatory elements (CREs) on transcriptional regulation and phenotypic outcomes remains challenging. Previously, we showed that the dominant Barren inflorescence3 (Bif3) mutant of maize (Zea mays) contains a duplicated copy of the homeobox transcription factor gene ZmWUSCHEL1 (ZmWUS1), named ZmWUS1-B. ZmWUS1-B is controlled by a spontaneously generated novel promoter region that dramatically increases its expression and alters patterning and development of young ears. Overexpression of ZmWUS1-B is caused by a unique enhancer region containing multimerized binding sites for type B RESPONSE REGULATORs (RRs), key transcription factors in cytokinin signaling. To better understand how the enhancer increases the expression of ZmWUS1 in vivo, we specifically targeted the ZmWUS1-B enhancer region by CRISPR-Cas9-mediated editing. A series of deletion events with different numbers of type B RR DNA binding motifs (AGATAT) enabled us to determine how the number of AGATAT motifs impacts in vivo expression of ZmWUS1-B and consequently ear development. In combination with dual-luciferase assays in maize protoplasts, our analysis reveals that AGATAT motifs have an additive effect on ZmWUS1-B expression, while the distance separating AGATAT motifs does not appear to have a meaningful impact, indicating that the enhancer activity derives from the sum of individual CREs. These results also suggest that in maize inflorescence development, there is a threshold of buffering capacity for ZmWUS1 overexpression.

Studying Salt-Induced Shifts in Gene Expression Patterns of Glucosinolate Transporters and Glucosinolate Accumulation in Two Contrasting Brassica Species

Brassica crops are well known for the accumulation of glucosinolates-secondary metabolites crucial for plants' adaptation to various stresses. Glucosinolates also functioning as defence compounds pose challenges to food quality due to their goitrogenic properties. Their disruption leaves plants susceptible to insect pests and diseases. Hence, a targeted reduction in seed glucosinolate content is of paramount importance to increase food acceptance. GLUCOSINOLATE TRANSPORTERS (GTRs) present a promising avenue for selectively reducing glucosinolate concentrations in seeds while preserving biosynthesis elsewhere. In this study, 54 putative GTR protein sequences found in Brassica were retrieved, employing Arabidopsis GTR1 and GTR2 templates. Comprehensive bioinformatics analyses, encompassing gene structure organization, domain analysis, motif assessments, promoter analysis, and cis-regulatory elements, affirmed the existence of transporter domains and stress-related regulatory elements. Phylogenetic analysis revealed patterns of conservation and divergence across species. Glucosinolates have been shown to increase under stress conditions, indicating a potential role in stress response. To elucidate the role of GTRs in glucosinolate transportation under NaCl stress in two distinct Brassica species, B. juncea and B. napus, plants were subjected to 0, 100, or 200 mM NaCl. Based on the literature, key GTR genes were chosen and their expression across various plant parts was assessed. Both species displayed divergent trends in their biochemical profiles as well as glucosinolate contents under elevated salt stress conditions. Statistical modelling identified significant contributors to glucosinolate variations, guiding the development of targeted breeding strategies for low-glucosinolate varieties. Notably, GTR2A2 exhibited pronounced expressions in stems, contributing approximately 52% to glucosinolate content variance, while GTR2B1/C2 displayed significant expression in flowers. Additionally, GTR2A1 and GTR1A2/B1 demonstrated noteworthy expression in roots. This study enhances our understanding of glucosinolate regulation under stress conditions, offering avenues to improve Brassica crop quality and resilience.

Genome-wide identification, phylogenetic classification of histone acetyltransferase genes, and their expression analysis in sugar beet (Beta vulgaris L.) under salt stress

MAIN CONCLUSION

This study identified seven histone acetyltransferase-encoding genes (HATs) from Beta vulgaris L. (sugar beet) genome through bioinformatics tools and analyzed their expression profiles under salt stress. Sugar beet HATs are phylogenetically divided into four families: GNAT, MYST, CBP, and TAFII250. The BvHAT genes were differentially transcribed in leaves, stems, and roots of B. vulgaris salt-resistant (Casino) and -sensitive (Bravo) cultivars under salt stress. Histone acetylation is regulated by histone acetyltransferases (HATs), which catalyze ɛ-amino bond formation between lysine residues and acetyl groups with a cofactor, acetyl-CoA. Even though the HATs are known to participate in stress response and development in model plants, little is known about the functions of HATs in crops. In sugar beet (Beta vulgaris L.), they have not yet been identified and characterized. Here, an in silico analysis of the HAT gene family in sugar beet was performed, and their expression patterns in leaves, stems, and roots of B. vulgaris were analyzed under salt stress. Salt-resistant (Casino) and -sensitive (Bravo) beet cultivars were used for gene expression assays. Seven HATs were identified from sugar beet genome, and named BvHAG1, BvHAG2, BvHAG3, BvHAG4, BvHAC1, BvHAC2, and BvHAF1. The HAT proteins were divided into 4 groups including MYST, GNAT (GCN5, HAT1, ELP3), CBP and TAFII250. Analysis of cis-acting elements indicated that the BvHAT genes might be involved in hormonal regulation, light response, plant development, and abiotic stress response. The BvHAT genes were differentially expressed in leaves, stems, and roots under control and 300 mM NaCl. In roots of B. vulgaris cv. Bravo, the BvHAG1, BvHAG2, BvHAG4, BvHAF1, and BvHAC1 genes were dramatically expressed after 7 and 14 days of salt stress. Interestingly, the BvHAC2 gene was not expressed under both control and stress conditions. However, the expression of BvHAG2, BvHAG3, BvHAG4, BvHAC1, BvHAC2 genes showed a significant increase in response to salt stress in the roots of cv. Casino. This study provides new insights into the potential roles of histone acetyltransferases in sugar beet.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

Optimized barley phytase gene expression by focused FIND-IT screening for mutations in cis-acting regulatory elements

Introduction

Induced modification of plant gene expression is of both fundamental and applied importance. Cis-acting regulatory elements (CREs) are major determinants of the spatiotemporal strength of gene expression. Yet, there are few examples where induced genetic variation in predetermined CREs has been exploited to improve or investigate crop plants.

Methods

The digital PCR based FIND-IT technology was applied to discover barley mutants with CRE variants in the promoter of the nutritional important barley grain phytase (PAPhy_a) gene.

Results and discussion

Mutants with higher or lower gene expression and ultimately higher or lower mature grain phytase activity (MGPA), respectively, were discovered. Field trials and inositol phosphate profiling during germination showed that PAPhy_a does not influence agronomic performance under the trial conditions but it does shorten the lag time of phosphate mobilization during germination. Higher endogenous MGPA is an improvement of grain quality for feed use as it improves the phosphate bioavailability for monogastric animals. Moreover, as the targeted CRE motifs of the PAPhy_a promoter are shared with a range of seed expressed genes like key cereal and legume storage genes, the current results demonstrates a concept for modulating individual gene expression levels of a range of seed genes.

A large sequenced mutant library - valuable reverse genetic resource that covers 98% of sorghum genes

Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum, a drought and heat-tolerant C4 species, requires a vast, large-scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large-scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)-induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis-regulatory elements. A comparison of the distribution of mutations in the large-scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS-induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5' UTR, gene body, and 3'-UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the Ka /Ks value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB (https://www.depts.ttu.edu/igcast/sorbmutdb.php) and SorghumBase (https://sorghumbase.org/). This current large-scale sequence-indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

Prioritizing Metabolic Gene Regulators through Multi-Omic Network Integration in Maize

Elucidating gene regulatory networks (GRNs) is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. In this study, we integrated publicly available co-expression networks derived from more than 6,000 RNA-seq samples, 283 protein-DNA interaction assays, and 16 million of SNPs used to identify expression quantitative loci (eQTL), to construct TF-target networks. In total, we analyzed ~4.6M interactions to generate four distinct types of TF-target networks: co-expression, protein-DNA interaction (PDI), trans-expression quantitative loci (trans-eQTL), and cis-eQTL combined with PDIs. To improve the functional annotation of TFs based on its target genes, we implemented three different strategies to integrate these four types of networks. We subsequently evaluated the effectiveness of our method through loss-of function mutant and random networks. The multi-network integration allowed us to identify transcriptional regulators of hormone-, metabolic- and development-related processes. Finally, using the topological properties of the fully integrated network, we identified potentially functional redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.)

BACKGROUND

N6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy.

RESULTS

In this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance.

CONCLUSIONS

This study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.

A rapid and scalable approach to build synthetic repetitive hormone-responsive promoters

Advancement of DNA-synthesis technologies has greatly facilitated the development of synthetic biology tools. However, high-complexity DNA sequences containing tandems of short repeats are still notoriously difficult to produce synthetically, with commercial DNA synthesis companies usually rejecting orders that exceed specific sequence complexity thresholds. To overcome this limitation, we developed a simple, single-tube reaction method that enables the generation of DNA sequences containing multiple repetitive elements. Our strategy involves commercial synthesis and PCR amplification of padded sequences that contain the repeats of interest, along with random intervening sequence stuffers that include type IIS restriction enzyme sites. GoldenBraid molecular cloning technology is then employed to remove the stuffers, rejoin the repeats together in a predefined order, and subclone the tandem(s) in a vector using a single-tube digestion-ligation reaction. In our hands, this new approach is much simpler, more versatile and efficient than previously developed solutions to this problem. As a proof of concept, two different phytohormone-responsive, synthetic, repetitive proximal promoters were generated and tested in planta in the context of transcriptional reporters. Analysis of transgenic lines carrying the synthetic ethylene-responsive promoter 10x2EBS-S10 fused to the GUS reporter gene uncovered several developmentally regulated ethylene response maxima, indicating the utility of this reporter for monitoring the involvement of ethylene in a variety of physiologically relevant processes. These encouraging results suggest that this reporter system can be leveraged to investigate the ethylene response to biotic and abiotic factors with high spatial and temporal resolution.

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.

GhBRX.1, GhBRX.2, and GhBRX4.3 improve resistance to salt and cold stress in upland cotton

Introduction

Abiotic stress during growth readily reduces cotton crop yield. The different survival tactics of plants include the activation of numerous stress response genes, such as BREVIS RADIX (BRX).

Methods

In this study, the BRX gene family of upland cotton was identified and analyzed by bioinformatics method, three salt-tolerant and cold-resistant GhBRX genes were screened. The expression of GhBRX.1, GhBRX.2 and GhBRXL4.3 in upland cotton was silenced by virus-induced gene silencing (VIGS) technique. The physiological and biochemical indexes of plants and the expression of related stress-response genes were detected before and after gene silencing. The effects of GhBRX.1, GhBRX.2 and GhBRXL4.3 on salt and cold resistance of upland cotton were further verified.

Results and discussion

We discovered 12, 6, and 6 BRX genes in Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively. Chromosomal localization indicated that the retention and loss of GhBRX genes on homologous chromosomes did not have a clear preference for the subgenomes. Collinearity analysis suggested that segmental duplications were the main force for BRX gene amplification. The upland cotton genes GhBRX.1, GhBRX.2 and GhBRXL4.3 are highly expressed in roots, and GhBRXL4.3 is also strongly expressed in the pistil. Transcriptome data and qRT‒PCR validation showed that abiotic stress strongly induced GhBRX.1, GhBRX.2 and GhBRXL4.3. Under salt stress and low-temperature stress conditions, the activities of superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) and the content of soluble sugar and chlorophyll decreased in GhBRX.1-, GhBRX.2- and GhBRXL4.3-silenced cotton plants compared with those in the control (TRV: 00). Moreover, GhBRX.1-, GhBRX.2- and GhBRXL4.3-silenced cotton plants exhibited greater malondialdehyde (MDA) levels than did the control plants. Moreover, the expression of stress marker genes (GhSOS1, GhSOS2, GhNHX1, GhCIPK6, GhBIN2, GhSnRK2.6, GhHDT4D, GhCBF1 and GhPP2C) decreased significantly in the three target genes of silenced plants following exposure to stress. These results imply that the GhBRX.1, GhBRX.2 and GhBRXL4.3 genes may be regulators of salt stress and low-temperature stress responses in upland cotton.

Predicting the unpredictable: the regulatory nature and promiscuity of herbicide cross resistance

The emergence of herbicide-resistant weeds is a significant threat to modern agriculture. Cross resistance, a phenomenon where resistance to one herbicide confers resistance to another, is a particular concern owing to its unpredictability. Nontarget-site (NTS) cross resistance is especially challenging to predict, as it arises from genes that encode enzymes that do not directly involve the herbicide target site and can affect multiple herbicides. Recent advancements in genomic and structural biology techniques could provide new venues for predicting NTS resistance in weed species. In this review, we present an overview of the latest approaches that could be used. We discuss the use of genomic and epigenomics techniques such as ATAC-seq and DAP-seq to identify transcription factors and cis-regulatory elements associated with resistance traits. Enzyme/protein structure prediction and docking analysis are discussed as an initial step for predicting herbicide binding affinities with key enzymes to identify candidates for subsequent in vitro validation. We also provide example analyses that can be deployed toward elucidating cross resistance and its regulatory patterns. Ultimately, our review provides important insights into the latest scientific advancements and potential directions for predicting and managing herbicide cross resistance in weeds. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Origin and early divergence of tandem duplicated sorbitol transporter genes in Rosaceae: insights from evolutionary analysis of the SOT gene family in angiosperms

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol transporters (SOTs) play a vital role in facilitating sorbitol allocation from source to sink organs and sugar accumulation in sink organs. While prior research has addressed gene duplications within the SOT gene family in Rosaceae, the precise origin and evolutionary dynamics of these duplications remain unclear, largely due to the complicated interplay of whole genome duplications and tandem duplications. Here, we investigated the synteny relationships among all identified Polyol/Monosaccharide Transporter (PLT) genes in 61 angiosperm genomes and SOT genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of PLTs and SOTs across diverse plant lineages. We found that Rosaceae SOTs, as PLT family members, originated from a pair of tandemly duplicated PLT genes within Class III-A. Furthermore, our investigation highlights the role of lineage-specific and synergistic duplications in Amygdaloideae in contributing to the expansion of SOTs in Rosaceae plants. Collectively, our findings provide insights into the genomic origins, duplication events, and subsequent divergence of SOT gene family members. Such insights lay a crucial foundation for comprehensive functional characterizations in future studies.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Comparative phylogenomics and phylotranscriptomics provide insights into the genetic complexity of nitrogen-fixing root-nodule symbiosis

Plant root-nodule symbiosis (RNS) with mutualistic nitrogen-fixing bacteria is restricted to a single clade of angiosperms, the Nitrogen-Fixing Nodulation Clade (NFNC), and is best understood in the legume family. Nodulating species share many commonalities, explained either by divergence from a common ancestor over 100 million years ago or by convergence following independent origins over that same time period. Regardless, comparative analyses of diverse nodulation syndromes can provide insights into constraints on nodulation-what must be acquired or cannot be lost for a functional symbiosis-and the latitude for variation in the symbiosis. However, much remains to be learned about nodulation, especially outside of legumes. Here, we employed a large-scale phylogenomic analysis across 88 species, complemented by 151 RNA-seq libraries, to elucidate the evolution of RNS. Our phylogenomic analyses further emphasize the uniqueness of the transcription factor NIN as a master regulator of nodulation and identify key mutations that affect its function across the NFNC. Comparative transcriptomic assessment revealed nodule-specific upregulated genes across diverse nodulating plants, while also identifying nodule-specific and nitrogen-response genes. Approximately 70% of symbiosis-related genes are highly conserved in the four representative species, whereas defense-related and host-range restriction genes tend to be lineage specific. Our study also identified over 900 000 conserved non-coding elements (CNEs), over 300 000 of which are unique to sampled NFNC species. NFNC-specific CNEs are enriched with the active H3K9ac mark and are correlated with accessible chromatin regions, thus representing a pool of candidate regulatory elements for genes involved in RNS. Collectively, our results provide novel insights into the evolution of nodulation and lay a foundation for engineering of RNS traits in agriculturally important crops.

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.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Genome-wide DNase I-hypersensitive site assay reveals distinct genomic distributions and functional features of open chromatin in autopolyploid sugarcane

The characterization of cis-regulatory DNA elements (CREs) is essential for deciphering the regulation of gene expression in eukaryotes. Although there have been endeavors to identify CREs in plants, the properties of CREs in polyploid genomes are still largely unknown. Here, we conducted the genome-wide identification of DNase I-hypersensitive sites (DHSs) in leaf and stem tissues of the auto-octoploid species Saccharum officinarum. We revealed that DHSs showed highly similar distributions in the genomes of these two S. officinarum tissues. Notably, we observed that approximately 74% of DHSs were located in distal intergenic regions, suggesting considerable differences in the abundance of distal CREs between S. officinarum and other plants. Leaf- and stem-dependent transcriptional regulatory networks were also developed by mining the binding motifs of transcription factors (TFs) from tissue-specific DHSs. Four TEOSINTE BRANCHED 1, CYCLOIDEA, and PCF1 (TCP) TFs (TCP2, TCP4, TCP7, and TCP14) and two ethylene-responsive factors (ERFs) (ERF109 and ERF03) showed strong causal connections with short binding distances from each other, pointing to their possible roles in the regulatory networks of leaf and stem development. Through functional validation in transiently transgenic protoplasts, we isolate a set of tissue-specific promoters. Overall, the DHS maps presented here offer a global view of the potential transcriptional regulatory elements in polyploid sugarcane and can be expected to serve as a valuable resource for both transcriptional network elucidation and genome editing in sugarcane breeding.

Distinct guard cell-specific remodeling of chromatin accessibility during abscisic acid- and CO2-dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures, thereby regulating gas exchange. Chromatin structure controls transcription factor (TF) access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remains unknown. Here, we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2 (carbon dioxide), regulate guard cell chromatin during stomatal movements. Our cell type-specific analyses uncover patterns of chromatin accessibility specific to guard cells and define cis-regulatory sequences supporting guard cell-specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell type specificity. DNA motif analyses uncover binding sites for distinct TFs enriched in ABA-induced and ABA-repressed chromatin. We identify the Abscisic Acid Response Element (ABRE) Binding Factor (ABF) bZIP-type TFs that are required for ABA-triggered chromatin opening in guard cells and roots and implicate the inhibition of a clade of bHLH-type TFs in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling, whereby elevated atmospheric CO2 had only minimal impact on chromatin dynamics. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

The long non-coding RNA LINDA restrains cellular collapse following DNA damage in Arabidopsis thaliana

The genomic integrity of every organism is endangered by various intrinsic and extrinsic stresses. To maintain genomic integrity, a sophisticated DNA damage response (DDR) network is activated rapidly after DNA damage. Notably, the fundamental DDR mechanisms are conserved in eukaryotes. However, knowledge about many regulatory aspects of the plant DDR is still limited. Important, yet little understood, regulatory factors of the DDR are the long non-coding RNAs (lncRNAs). In humans, 13 lncRNAs functioning in DDR have been characterized to date, whereas no such lncRNAs have been characterized in plants yet. By meta-analysis, we identified the putative long intergenic non-coding RNA induced by DNA damage (LINDA) that responds strongly to various DNA double-strand break-inducing treatments, but not to replication stress induced by mitomycin C. After DNA damage, LINDA is rapidly induced in an ATM- and SOG1-dependent manner. Intriguingly, the transcriptional response of LINDA to DNA damage is similar to that of its flanking hypothetical protein-encoding gene. Phylogenetic analysis of putative Brassicales and Malvales LINDA homologs indicates that LINDA lncRNAs originate from duplication of a flanking small protein-encoding gene followed by pseudogenization. We demonstrate that LINDA is not only needed for the regulation of this flanking gene but also fine-tuning of the DDR after the occurrence of DNA double-strand breaks. Moreover, Δlinda mutant root stem cells are unable to recover from DNA damage, most likely due to hyper-induced cell death.

The m6 A reader SiYTH1 enhances drought tolerance by affecting the messenger RNA stability of genes related to stomatal closure and reactive oxygen species scavenging in Setaria italica

Foxtail millet (Setaria italica), a vital drought-resistant crop, plays a significant role in ensuring food and nutritional security. However, its drought resistance mechanism is not fully understood. N6 -methyladenosine (m6 A) modification of RNA, a prevalent epi-transcriptomic modification in eukaryotes, provides a binding site for m6 A readers and affects plant growth and stress responses by regulating RNA metabolism. In this study, we unveiled that the YT521-B homology (YTH) family gene SiYTH1 positively regulated the drought tolerance of foxtail millet. Notably, the siyth1 mutant exhibited reduced stomatal closure and augmented accumulation of excessive H2 O2 under drought stress. Further investigations demonstrated that SiYTH1 positively regulated the transcripts harboring m6 A modification related to stomatal closure and reactive oxygen species (ROS) scavenging under drought stress. SiYTH1 was uniformly distributed in the cytoplasm of SiYTH1-GFP transgenic foxtail millet. It formed dynamic liquid-like SiYTH1 cytosol condensates in response to drought stress. Moreover, the cytoplasmic protein SiYTH1 was identified as a distinct m6 A reader, facilitating the stabilization of its directly bound SiARDP and ROS scavenging-related transcripts under drought stress. Furthermore, natural variation analysis revealed SiYTH1AGTG as the dominant allele responsible for drought tolerance in foxtail millet. Collectively, this study provides novel insights into the intricate mechanism of m6 A reader-mediated drought tolerance and presents a valuable genetic resource for improving drought tolerance in foxtail millet breeding.

Computational tools for inferring transcription factor activity

Transcription factors (TFs) are essential players in orchestrating the regulatory landscape in cells. Still, their exact modes of action and dependencies on other regulatory aspects remain elusive. Since TFs act cell type-specific and each TF has its own characteristics, untangling their regulatory interactions from an experimental point of view is laborious and convoluted. Thus, there is an ongoing development of computational tools that estimate transcription factor activity (TFA) from a variety of data modalities, either based on a mapping of TFs to their putative target genes or in a genome-wide, gene-unspecific fashion. These tools can help to gain insights into TF regulation and to prioritize candidates for experimental validation. We want to give an overview of available computational tools that estimate TFA, illustrate examples of their application, debate common result validation strategies, and discuss assumptions and concomitant limitations.

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.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters that lack any known promoter elements, thus designated as Coreless, tend to be expressed more uniformly. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression pattern. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Decoding the guardians of cotton resilience: A comprehensive exploration of the βCA genes and its role in Verticillium dahliae resistance

Plant Carbonic anhydrases (Cas) have been shown to be stress-responsive enzymes that may play a role in adapting to adverse conditions. Cotton is a significant economic crop in China, with upland cotton (Gossypium hirsutum) being the most widely cultivated species. We conducted genome-wide identification of the βCA gene in six cotton species and preliminary analysis of the βCA gene in upland cotton. In total, 73 βCA genes from six cotton species were identified, with phylogenetic analysis dividing them into five subgroups. GHβCA proteins were predominantly localized in the chloroplast and cytoplasm. The genes exhibited conserved motifs, with motifs 1, 2, and 3 being prominent. GHβCA genes were unevenly distributed across chromosomes and were associated with stress-responsive cis-regulatory elements, including those responding to light, MeJA, salicylic acid, abscisic acid, cell cycle regulation, and defence/stress. Expression analysis indicated that GHβCA6, GHβCA7, GHβCA10, GHβCA15, and GHβCA16 were highly expressed under various abiotic stress conditions, whereas GHβCA3, GHβCA9, GHβCA10, and GHβCA18 had higher expression patterns under Verticillium dahliae infection at different time intervals. In Gossypium thurberi, GthβCA1, GthβCA2, and GthβCA4 showed elevated expression across stress conditions and tissues. Silencing GHβCA10 through VIGS increased Verticillium wilt severity and reduced lignin deposition compared to non-silenced plants. GHβCA10 is crucial for cotton's defense against Verticillium dahliae. Further research is needed to understand the underlying mechanisms and develop strategies to enhance resistance against Verticillium wilt.

Differential gene expression during floral transition in pineapple

Pineapple (Ananas comosus var. comosus) and ornamental bromeliads are commercially induced to flower by treatment with ethylene or its analogs. The apex is transformed from a vegetative to a floral meristem and shows morphological changes in 8 to 10 days, with flowers developing 8 to 10 weeks later. During eight sampling stages ranging from 6 h to 8 days after treatment, 7961 genes were found to exhibit differential expression (DE) after the application of ethylene. In the first 3 days after treatment, there was little change in ethylene synthesis or in the early stages of the ethylene response. Subsequently, three ethylene response transcription factors (ERTF) were up-regulated and the potential gene targets were predicted to be the positive flowering regulator CONSTANS-like 3 (CO), a WUSCHEL gene, two APETALA1/FRUITFULL (AP1/FUL) genes, an epidermal patterning gene, and a jasmonic acid synthesis gene. We confirm that pineapple has lost the flowering repressor FLOWERING LOCUS C. At the initial stages, the SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was not significantly involved in this transition. Another WUSCHEL gene and a PHD homeobox transcription factor, though not apparent direct targets of ERTF, were up-regulated within a day of treatment, their predicted targets being the up-regulated CO, auxin response factors, SQUAMOSA, and histone H3 genes with suppression of abscisic acid response genes. The FLOWERING LOCUS T (FT), TERMINAL FLOWER (TFL), AGAMOUS-like APETELAR (AP2), and SEPETALA (SEP) increased rapidly within 2 to 3 days after ethylene treatment. Two FT genes were up-regulated at the apex and not at the leaf bases after treatment, suggesting that transport did not occur. These results indicated that the ethylene response in pineapple and possibly most bromeliads act directly to promote the vegetative to flower transition via APETALA1/FRUITFULL (AP1/FUL) and its interaction with SPL, FT, TFL, SEP, and AP2. A model based on AP2/ERTF DE and predicted DE target genes was developed to give focus to future research. The identified candidate genes are potential targets for genetic manipulation to determine their molecular role in flower transition.

Genetic diversity, transcript heterogeneity and allele mining of TaSKP1-6B-4 gene variants across diverse genotypes under terminal heat stress and genome wide characterization of TaSKP1 gene family from bread wheat (Triticum aestivum L.)

SKP1 (S-phase kinase protein1) is an essential regulatory component of SCF (Skp1-cullin-F-box) E3 ubiquitin ligases involved in maintenance of cellular protein homeostasis through ubiquitin mediated proteasome system (UPS). UPS play a key role in stress response and grain yield. Earlier, we isolated TaSKP1-6B-4, highly induced in flag leaf tissues (Accession No. KJ830759.1) of developing wheat caryopses under heat stress. To further assess the functional role of SKP1, genetic variability analysis was carried out in a panel of 25 contrasting germplasm through extensive phenotyping and transcript profiling of TaSKP1-6B-4 during anthesis under ambient and terminal heat stress (THS) in field experiments for two consecutive years. The analysis of variance revealed significant variations for all the traits studied. Higher H2(%), GCV, PCV, GA and GA% mean observed in tiller number per plant (23.81, 17.65, 5.71, 28, 30.86%) and grain number per head (30.27, 82.79, 60.16, 105.00, 108.64%) under THS over ambient temperature. Higher fold induction of TaSKP1-6B-4 transcripts was recorded in 10 genotypes viz. HD2967 (9.9), IC145456 (6.18) in flag leaf; while C-306 (15.88), RAJ3765 (8.37) in ear head. Allele mining of SKP1-6B-4 showed genotypic sequence variations. Whole genome wide search of SKP1 gene family identified 95 SKP1 genes which were structurally characterized. Grain yield, leaf senescence and other agronomic-morpho-physiological parameters combined with transcript profiling, cvHD2967, was found to be the best positively responsive to THS which by pedigree was not heat tolerant. We report a novel 2 year comprehensive field based analysis on collective genetic variability and SKP1/UPS modulation under a natural environmental setting. The data reveals potential functional role of UPS under THS and tolerant cultivars can be further utilized for clarifying the role of UPS mechanistically at the molecular level and for developing terminal heat stress tolerant wheat.

Distinct guard cell specific remodeling of chromatin accessibility during abscisic acid and CO2 dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures thereby regulating gas exchange. Chromatin structure controls transcription factor access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remain unknown. Here we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2, regulate guard cell chromatin during stomatal movements. Our cell type specific analyses uncover patterns of chromatin accessibility specific to guard cells and define novel cis-regulatory sequences supporting guard cell specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell-type specificity. DNA motif analyses uncover binding sites for distinct transcription factors enriched in ABA-induced and ABA-repressed chromatin. We identify the ABF/AREB bZIP-type transcription factors that are required for ABA-triggered chromatin opening in guard cells and implicate the inhibition of a set of bHLH-type transcription factors in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Engineering custom morpho- and chemotypes of Populus for sustainable production of biofuels, bioproducts, and biomaterials

Humans have been modifying plant traits for thousands of years, first through selection (i.e., domestication) then modern breeding, and in the last 30 years, through biotechnology. These modifications have resulted in increased yield, more efficient agronomic practices, and enhanced quality traits. Precision knowledge of gene regulation and function through high-resolution single-cell omics technologies, coupled with the ability to engineer plant genomes at the DNA sequence, chromatin accessibility, and gene expression levels, can enable engineering of complex and complementary traits at the biosystem level. Populus spp., the primary genetic model system for woody perennials, are among the fastest growing trees in temperate zones and are important for both carbon sequestration and global carbon cycling. Ample genomic and transcriptomic resources for poplar are available including emerging single-cell omics datasets. To expand use of poplar outside of valorization of woody biomass, chassis with novel morphotypes in which stem branching and tree height are modified can be fabricated thereby leading to trees with altered leaf to wood ratios. These morphotypes can then be engineered into customized chemotypes that produce high value biofuels, bioproducts, and biomaterials not only in specific organs but also in a cell-type-specific manner. For example, the recent discovery of triterpene production in poplar leaf trichomes can be exploited using cell-type specific regulatory sequences to synthesize high value terpenes such as the jet fuel precursor bisabolene specifically in the trichomes. By spatially and temporally controlling expression, not only can pools of abundant precursors be exploited but engineered molecules can be sequestered in discrete cell structures in the leaf. The structural diversity of the hemicellulose xylan is a barrier to fully utilizing lignocellulose in biomaterial production and by leveraging cell-type-specific omics data, cell wall composition can be modified in a tailored and targeted specific manner to generate poplar wood with novel chemical features that are amenable for processing or advanced manufacturing. Precision engineering poplar as a multi-purpose sustainable feedstock highlights how genome engineering can be used to re-imagine a crop species.

Transcription factors KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4 regulate long intergenic noncoding RNAs expressed in Arabidopsis roots

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of 8 annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis (Arabidopsis thaliana). For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation (ChIP) data were used to construct a TF-lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation, and regulation differ between plants and animals. Experimental validation confirmed the role of 3 TFs, KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4, as key regulators controlling root-specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell type-specific expression or chromatin accessibility, which are linked with genome-wide association studies genetic associations related to root system development and growth. The multilevel genome-wide characterization covering chromatin state information, promoter conservation, and chromatin immunoprecipitation-based TF binding, for all detectable lincRNAs across 769 expression samples, permits rapidly defining the biological context and relevance of Arabidopsis lincRNAs through regulatory networks.

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.