Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis

CRISPR-mediated iron and folate biofortification in crops: advances and perspectives

Micronutrient deficiency conditions, such as anemia, are the most prevalent global health problem due to inadequate iron and folate in dietary sources. Biofortification advancements can propel the rapid amelioration of nutritionally beneficial components in crops that are required to combat the adverse effects of micronutrient deficiencies on human health. To date, several strategies have been proposed to increase micronutrients in plants to improve food quality, but very few approaches have intrigued `clustered regularly interspaced short palindromic repeats' (CRISPR) modules for the enhancement of iron and folate concentration in the edible parts of plants. In this review, we discuss two important approaches to simultaneously enhance the bioavailability of iron and folate concentrations in rice endosperms by utilizing advanced CRISPR-Cas9-based technology. This includes the 'tuning of cis-elements' and 'enhancer re-shuffling' in the regulatory components of genes that play a vital role in iron and folate biosynthesis/transportation pathways. In particular, base-editing and enhancer re-installation in native promoters of selected genes can lead to enhanced accumulation of iron and folate levels in the rice endosperm. The re-distribution of micronutrients in specific plant organs can be made possible using the above-mentioned contemporary approaches. Overall, the present review discusses the possible approaches for synchronized iron and folate biofortification through modification in regulatory gene circuits employing CRISPR-Cas9 technology.

Comprehensive mapping and modelling of the rice regulome landscape unveils the regulatory architecture underlying complex traits

Unraveling the regulatory mechanisms that govern complex traits is pivotal for advancing crop improvement. Here we present a comprehensive regulome atlas for rice (Oryza sativa), charting the chromatin accessibility across 23 distinct tissues from three representative varieties. Our study uncovers 117,176 unique open chromatin regions (OCRs), accounting for ~15% of the rice genome, a notably higher proportion compared to previous reports in plants. Integrating RNA-seq data from matched tissues, we confidently predict 59,075 OCR-to-gene links, with enhancers constituting 69.54% of these associations, including many known enhancer-to-gene links. Leveraging this resource, we re-evaluate genome-wide association study results and discover a previously unknown function of OsbZIP06 in seed germination, which we subsequently confirm through experimental validation. We optimize deep learning models to decode regulatory grammar, achieving robust modeling of tissue-specific chromatin accessibility. This approach allows to predict cross-variety regulatory dynamics from genomic sequences, shedding light on the genetic underpinnings of cis-regulatory divergence and morphological disparities between varieties. Overall, our study establishes a foundational resource for rice functional genomics and precision molecular breeding, providing valuable insights into regulatory mechanisms governing complex traits.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

Intrapuparial stage aging and PMI estimation based on the developmental transcriptomes of forensically important Aldrichina grahami (Diptera: Calliphoridae) gene expression

Background

The expression profiles of differentially expressed genes (DEGs) during pupal development have been demonstrated to be vital in age estimation of forensic entomological study. Here, using forensically important Aldrichina grahami (Diptera: Calliphoridae), we aimed to explore the potential of intrapuparial stage aging and postmortem interval (PMI) estimation based on characterization of successive developmental transcriptomes and gene expression patterns.

Methods

We collected A. grahami pupae at 11 successive intrapuparial stages at 20 °C and used the RNA-seq technique to build the transcriptome profiles of their intrapuparial stages. The DEGs were identified during the different intrapuparial stages using comparative transcriptome analysis. The selected marker DEGs were classified and clustered for intrapuparial stage aging and PMI estimation and then further verified for transcriptome data validation. Ultimately, we categorized the overall gene expression levels as the dependent variable and the age of intrapuparial A. grahami as the independent variable to conduct nonlinear regression analysis.

Results

We redefined the intrapuparial stages of A. grahami into five key successive substages (I, II, III, IV, and V), based on the overall gene expression patterns during pupal development. We screened 99 specific time-dependent expressed genes (stage-specific DEGs) to determine the different intrapuparial stages based on comparison of the gene expression levels during the 11 developmental intrapuparial stages of A. grahami. We observed that 55 DEGs showed persistent upregulation during the development of intrapuparial A. grahami. We then selected four DEGs (act79b, act88f, up and ninac) which presented consistent upregulation using RT-qPCR (quantitative real-time PCR) analysis, along with consideration of the maximum fold changes during the pupal development. We conducted nonlinear regression analysis to simulate the calculations of the relationships between the expression levels of the four selected DEGs and the developmental time of intrapuparial A. grahami and constructed fitting curves. The curves demonstrated that act79b and ninac showed continuous relatively increasing levels.

Conclusions

This study redefined the intrapuparial stages of A. grahami based on expression profiles of developmental transcriptomes for the first time. The stage-specific DEGs and those with consistent tendencies of expression were found to have potential in age estimation of intrapuparial A. grahami and could be supplementary to a more accurate prediction of PMI.

Plant enhancers exhibit both cooperative and additive interactions among their functional elements

Enhancers are cis-regulatory elements that shape gene expression in response to numerous developmental and environmental cues. In animals, several models have been proposed to explain how enhancers integrate the activity of multiple transcription factors. However, it remains largely unclear how plant enhancers integrate transcription factor activity. Here, we use Plant STARR-seq to characterize 3 light-responsive plant enhancers-AB80, Cab-1, and rbcS-E9-derived from genes associated with photosynthesis. Saturation mutagenesis revealed mutations, many of which clustered in short regions, that strongly reduced enhancer activity in the light, in the dark, or in both conditions. When tested in the light, these mutation-sensitive regions did not function on their own; rather, cooperative interactions with other such regions were required for full activity. Epistatic interactions occurred between mutations in adjacent mutation-sensitive regions, and the spacing and order of mutation-sensitive regions in synthetic enhancers affected enhancer activity. In contrast, when tested in the dark, mutation-sensitive regions acted independently and additively in conferring enhancer activity. Taken together, this work demonstrates that plant enhancers show evidence for both cooperative and additive interactions among their functional elements. This knowledge can be harnessed to design strong, condition-specific synthetic enhancers.

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

A distal enhancer guides the negative selection of toxic glycoalkaloids during tomato domestication

Steroidal glycoalkaloids (SGAs) are major plant defense metabolites against pests, while they are considered poisonous in food. The genetic basis that guides negative selection of SGAs production during tomato domestication remains poorly understood. Here, we identify a distal enhancer, GAME Enhancer 1 (GE1), as the key regulator of SGAs metabolism in tomato. GE1 recruits MYC2-GAME9 transcriptional complex to regulate the expression of GAME cluster genes via the formation of chromatin loops located in the neighboring DNA region. A naturally occurring GE176 allelic variant is found to be more active in stimulating GAME expression. We show that the weaker GE1 allele has been the main driver for selecting reduced SGAs levels during tomato domestication. Unravelling the "TFs-Enhancer-Promoter" regulatory mechanism operating in SGAs metabolism opens unprecedented prospects for SGAs manipulation in Solanaceae via precision breeding strategies.

Beyond heat waves: Unlocking epigenetic heat stress memory in Arabidopsis

Plants remember their exposure to environmental changes and respond more effectively the next time they encounter a similar change by flexibly altering gene expression. Epigenetic mechanisms play a crucial role in establishing such memory of environmental changes and fine-tuning gene expression. With the recent advancements in biochemistry and sequencing technologies, it has become possible to characterize the dynamics of epigenetic changes on scales ranging from short term (minutes) to long term (generations). Here, our main focus is on describing the current understanding of the temporal regulation of histone modifications and chromatin changes during exposure to short-term recurring high temperatures and reevaluating them in the context of natural environments. Investigations of the dynamics of histone modifications and chromatin structural changes in Arabidopsis after repeated exposure to heat at short intervals have revealed the detailed molecular mechanisms of short-term heat stress memory, which include histone modification enzymes, chromatin remodelers, and key transcription factors. In addition, we summarize the spatial regulation of heat responses. Based on the natural temperature patterns during summer, we discuss how plants cope with recurring heat stress occurring at various time intervals by utilizing 2 distinct types of heat stress memory mechanisms. We also explore future research directions to provide a more precise understanding of the epigenetic regulation of heat stress memory.

Dynamic transcriptome analysis provides molecular insights into underground floral differentiation in Adonis Amurensis Regel & Radde

Understanding flower developmental processes is a prerequisite for improving flowering 'plants' production. Adonis amurensis is a fascinating spring ephemeral plant that develops its flower organs underground. Nevertheless, knowledge of the molecular mechanisms driving this particular process is scarce. Herein, we examined transcriptional changes during underground flower differentiation in A. amurensis and unveiled key differently regulated genes and pathways. High-throughput RNA sequencing of meristems at different flower developmental stages, including flower primordium (FP), sepal stage (SE), perianth primordium (PE), stamen stage (ST), and pistil stage (PI), identified 303,234 unigenes that showed 44.79% similarity with sequences in Aquilegia coerulea. Correlations, principal component, and differentially expressed genes (DEGs) analyses revealed that few molecular changes occurred during the transition from PE to ST. Many DEGs exhibited stage-specific regulations. Transcription factor (TF) and phytohormone family genes are critical regulators of the floral differentiation process in A. amurensis. The most differentially regulated TFs were MADS, FAR1, MYBs, AP2/ERF, B3, C2H2, and LOBs. We filtered out 186 candidate genes for future functional studies, including 18 flowering/circadian-related, 32 phytohormone-related, and TF family genes. Our findings deepen our understanding of the underground flower differentiation process and offer critical resources to dissect its regulatory network in A. amurensis. These findings establish a foundational platform for researchers dedicated to exploring the unique phenotypic characteristics of this specific flowering modality and delving into the intricate molecular mechanisms underpinning its regulation and expression.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

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.

Genome-wide enhancer identification by massively parallel reporter assay in Arabidopsis

Enhancers are critical cis-regulatory elements controlling gene expression during cell development and differentiation. However, genome-wide enhancer characterization has been challenging due to the lack of a well-defined relationship between enhancers and genes. Function-based methods are the gold standard for determining the biological function of cis-regulatory elements; however, these methods have not been widely applied to plants. Here, we applied a massively parallel reporter assay on Arabidopsis to measure enhancer activities across the genome. We identified 4327 enhancers with various combinations of epigenetic modifications distinctively different from animal enhancers. Furthermore, we showed that enhancers differ from promoters in their preference for transcription factors. Although some enhancers are not conserved and overlap with transposable elements forming clusters, enhancers are generally conserved across thousand Arabidopsis accessions, suggesting they are selected under evolution pressure and could play critical roles in the regulation of important genes. Moreover, comparison analysis reveals that enhancers identified by different strategies do not overlap, suggesting these methods are complementary in nature. In sum, we systematically investigated the features of enhancers identified by functional assay in A. thaliana, which lays the foundation for further investigation into enhancers' functional mechanisms in plants.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Low temperature-induced regulatory network rewiring via WRKY regulators during banana peel browning

Banana (Musa spp.) fruits, as typical tropical fruits, are cold sensitive, and lower temperatures can disrupt cellular compartmentalization and lead to severe browning. How tropical fruits respond to low temperature compared to the cold response mechanisms of model plants remains unknown. Here, we systematically characterized the changes in chromatin accessibility, histone modifications, distal cis-regulatory elements, transcription factor binding, and gene expression levels in banana peels in response to low temperature. Dynamic patterns of cold-induced transcripts were generally accompanied by concordant chromatin accessibility and histone modification changes. These upregulated genes were enriched for WRKY binding sites in their promoters and/or active enhancers. Compared to banana peel at room temperature, large amounts of banana WRKYs were specifically induced by cold and mediated enhancer-promoter interactions regulating critical browning pathways, including phospholipid degradation, oxidation, and cold tolerance. This hypothesis was supported by DNA affinity purification sequencing, luciferase reporter assays, and transient expression assay. Together, our findings highlight widespread transcriptional reprogramming via WRKYs during banana peel browning at low temperature and provide an extensive resource for studying gene regulation in tropical plants in response to cold stress, as well as potential targets for improving cold tolerance and shelf life of tropical fruits.

Promoter editing for the genetic improvement of crops

Gene expression plays a fundamental role in the regulation of agronomically important traits in crop plants. The genetic manipulation of plant promoters through genome editing has emerged as an effective strategy to create favorable traits in crops by altering the expression pattern of the pertinent genes. Promoter editing can be applied in a directed manner, where nucleotide sequences associated with favorable traits are precisely generated. Alternatively, promoter editing can also be exploited as a random mutagenic approach to generate novel genetic variations within a designated promoter, from which elite alleles are selected based on their phenotypic effects. Pioneering studies have demonstrated the potential of promoter editing in engineering agronomically important traits as well as in mining novel promoter alleles valuable for plant breeding. In this review, we provide an update on the application of promoter editing in crops for increased yield, enhanced tolerance to biotic and abiotic stresses, and improved quality. We also discuss several remaining technical bottlenecks and how this strategy may be better employed for the genetic improvement of crops in the future.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

AGAMOUS regulates various target genes via cell cycle-coupled H3K27me3 dilution in floral meristems and stamens

The MADS domain transcription factor AGAMOUS (AG) regulates floral meristem termination by preventing maintenance of the histone modification lysine 27 of histone H3 (H3K27me3) along the KNUCKLES (KNU) coding sequence. At 2 d after AG binding, cell division has diluted the repressive mark H3K27me3, allowing activation of KNU transcription prior to floral meristem termination. However, how many other downstream genes are temporally regulated by this intrinsic epigenetic timer and what their functions are remain unknown. Here, we identify direct AG targets regulated through cell cycle-coupled H3K27me3 dilution in Arabidopsis thaliana. Expression of the targets KNU, AT HOOK MOTIF NUCLEAR LOCALIZED PROTEIN18 (AHL18), and PLATZ10 occurred later in plants with longer H3K27me3-marked regions. We established a mathematical model to predict timing of gene expression and manipulated temporal gene expression using the H3K27me3-marked del region from the KNU coding sequence. Increasing the number of del copies delayed and reduced KNU expression in a polycomb repressive complex 2- and cell cycle-dependent manner. Furthermore, AHL18 was specifically expressed in stamens and caused developmental defects when misexpressed. Finally, AHL18 bound to genes important for stamen growth. Our results suggest that AG controls the timing of expression of various target genes via cell cycle-coupled dilution of H3K27me3 for proper floral meristem termination and stamen development.

Structural diversity and stress regulation of the plant immunity-associated CALMODULIN-BINDING PROTEIN 60 (CBP60) family of transcription factors in Solanum lycopersicum (tomato)

Cellular signaling generates calcium (Ca2+) ions, which are ubiquitous secondary messengers decoded by calcium-dependent protein kinases, calcineurins, calreticulin, calmodulins (CAMs), and CAM-binding proteins. Previous studies in the model plant Arabidopsis thaliana have shown the critical roles of the CAM-BINDING PROTEIN 60 (CBP60) protein family in plant growth, stress responses, and immunity. Certain CBP60 factors can regulate plant immune responses, like pattern-triggered immunity, effector-triggered immunity, and synthesis of major plant immune-activating metabolites salicylic acid (SA) and N-hydroxypipecolic acid (NHP). Although homologous CBP60 sequences have been identified in the plant kingdom, their function and regulation in most species remain unclear. In this paper, we specifically characterized 11 members of the CBP60 family in the agriculturally important crop tomato (Solanum lycopersicum). Protein sequence analyses revealed that three CBP60 homologs have the closest amino acid identity to Arabidopsis CBP60g and SARD1, master transcription factors involved in plant immunity. Strikingly, AlphaFold deep learning-assisted prediction of protein structures highlighted close structural similarity between these tomato and Arabidopsis CBP60 homologs. Conserved domain analyses revealed that they possess CAM-binding domains and DNA-binding domains, reflecting their potential involvement in linking Ca2+ signaling and transcriptional regulation in tomato plants. In terms of their gene expression profiles under biotic (Pseudomonas syringae pv. tomato DC3000 pathogen infection) and/or abiotic stress (warming temperatures), five tomato CBP60 genes were pathogen-responsive and temperature-sensitive, reminiscent of Arabidopsis CBP60g and SARD1. Overall, we present a genome-wide identification of the CBP60 gene/protein family in tomato plants, and we provide evidence on their regulation and potential function as Ca2+-sensing transcriptional regulators.

Isolation of 5' regulatory region of COLD1 gene and its functional characterization through transient expression analysis in tobacco and sugarcane

Chilling Tolerant Divergence 1 (COLD1) gene consists of Golgi pH Receptor (GPHR) as well as Abscisic Acid-linked G Protein-Coupled Receptor (ABA_GPCR), which are the major transmembrane proteins in plants. This gene expression has been found to be differentially regulated, under various stress conditions, in wild Saccharum-related genera, Erianthus arundinaceus, compared to commercial sugarcane variety. In this study, Rapid Amplification of Genomic Ends (RAGE) technique was employed to isolate the 5' upstream region of COLD1 gene to gain knowledge about the underlying stress regulatory mechanism. The current study established the cis-acting elements, main promoter regions, and Transcriptional Start Site (TSS) present within the isolated 5' upstream region (Cold1P) of COLD1, with the help of specific bioinformatics techniques. Phylogenetic analysis results revealed that the isolated Cold1P promoter is closely related to the species, Sorghum bicolor. Cold1P promoter-GUS gene construct was generated in pCAMBIA 1305.1 vector that displayed a constitutive expression of the GUS reporter gene in both monocot as well as dicot plants. The histochemical GUS assay outcomes confirmed that Cold1P can drive expression in both monocot as well as dicot plants. Cold1P's activities under several abiotic stresses such as cold, heat, salt, and drought, revealed its differential expression profile in commercial sugarcane variety. The highest activity of the GUS gene was found after 24 h of cold stress, driven by the isolated Cold1P promoter. The outcomes from GUS fluorimetric assay correlated with that of the GUS expression findings. This is the first report on Cold1P isolated from the species, E. arundinaceus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03650-8.

Enhancer activation via TCP and HD-ZIP and repression by Dof transcription factors mediate giant cell-specific expression

Proper cell-type identity relies on highly coordinated regulation of gene expression. Regulatory elements such as enhancers can produce cell type-specific expression patterns, but the mechanisms underlying specificity are not well understood. We previously identified an enhancer region capable of driving specific expression in giant cells, which are large, highly endoreduplicated cells in the Arabidopsis thaliana sepal epidermis. In this study, we use the giant cell enhancer as a model to understand the regulatory logic that promotes cell type-specific expression. Our dissection of the enhancer revealed that giant cell specificity is mediated primarily through the combination of two activators and one repressor. HD-ZIP and TCP transcription factors are involved in the activation of expression throughout the epidermis. High expression of HD-ZIP transcription factor genes in giant cells promoted higher expression driven by the enhancer in giant cells. Dof transcription factors repressed the activity of the enhancer such that only giant cells maintained enhancer activity. Thus, our data are consistent with a conceptual model whereby cell type-specific expression emerges from the combined activities of three transcription factor families activating and repressing expression in epidermal cells.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Predicting enhancer-promoter interaction based on epigenomic signals

Introduction: The physical interactions between enhancers and promoters are often involved in gene transcriptional regulation. High tissue-specific enhancer-promoter interactions (EPIs) are responsible for the differential expression of genes. Experimental methods are time-consuming and labor-intensive in measuring EPIs. An alternative approach, machine learning, has been widely used to predict EPIs. However, most existing machine learning methods require a large number of functional genomic and epigenomic features as input, which limits the application to different cell lines. Methods: In this paper, we developed a random forest model, HARD (H3K27ac, ATAC-seq, RAD21, and Distance), to predict EPI using only four types of features. Results: Independent tests on a benchmark dataset showed that HARD outperforms other models with the fewest features. Discussion: Our results revealed that chromatin accessibility and the binding of cohesin are important for cell-line-specific EPIs. Furthermore, we trained the HARD model in the GM12878 cell line and performed testing in the HeLa cell line. The cross-cell-lines prediction also performs well, suggesting it has the potential to be applied to other cell lines.

PlantCADB: A Comprehensive Plant Chromatin Accessibility Database

Chromatin accessibility landscapes are essential for detecting regulatory elements, illustrating the corresponding regulatory networks, and, ultimately, understanding the molecular basis underlying key biological processes. With the advancement of sequencing technologies, a large volume of chromatin accessibility data has been accumulated and integrated for humans and other mammals. These data have greatly advanced the study of disease pathogenesis, cancer survival prognosis, and tissue development. To advance the understanding of molecular mechanisms regulating plant key traits and biological processes, we developed a comprehensive plant chromatin accessibility database (PlantCADB) from 649 samples of 37 species. These samples are abiotic stress-related (such as heat, cold, drought, and salt; 159 samples), development-related (232 samples), and/or tissue-specific (376 samples). Overall, 18,339,426 accessible chromatin regions (ACRs) were compiled. These ACRs were annotated with genomic information, associated genes, transcription factor footprint, motif, and single-nucleotide polymorphisms (SNPs). Additionally, PlantCADB provides various tools to visualize ACRs and corresponding annotations. It thus forms an integrated, annotated, and analyzed plant-related chromatin accessibility resource, which can aid in better understanding genetic regulatory networks underlying development, important traits, stress adaptations, and evolution.PlantCADB is freely available at https://bioinfor.nefu.edu.cn/PlantCADB/.

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

Development of homozygous tetraploid potato and whole genome doubling-induced the enrichment of H3K27ac and potentially enhanced resistance to cold-induced sweetening in tubers

Polyploid plants typically display advantages on some agronomically important traits over their diploid counterparts. Extensive studies have shown genetic, transcriptomic, and epigenetic dynamics upon polyploidization in multiple plant species. However, few studies have unveiled those alternations imposed only by ploidy level, without any interference from heterozygosity. Cultivated potato is highly heterozygous. Thus, in this study, we developed two homozygous autotetraploid lines and one homozygous diploid line in parallel from a homozygous diploid potato. We confirmed their ploidy levels using chloroplast counting and karyotyping. Oligo-FISH and genome re-sequencing validated that these potato lines are nearly homozygous. We investigated variations in phenotypes, transcription, and histone modifications between two ploidies. Both autotetraploid lines produced larger but fewer tubers than the diploid line. Interestingly, each autotetraploid line displayed ploidy-related differential expression for various genes. We also discovered a genome-wide enrichment of H3K27ac in genic regions upon whole-genome doubling (WGD). However, such enrichment was not associated with the differential gene expression between two ploidies. The tetraploid lines may exhibit better resistance to cold-induced sweetening (CIS) than the diploid line in tubers, potentially regulated through the expression of CIS-related key genes, which seems to be associated with the levels of H3K4me3 in cold-stored tubers. These findings will help to understand the impacts of autotetraploidization on dynamics of phenotypes, transcription, and histone modifications, as well as on CIS-related genes in response to cold storage.

eRNA profiling uncovers the enhancer landscape of oesophageal adenocarcinoma and reveals new deregulated pathways

Cancer is driven by both genetic and epigenetic changes that impact on gene expression profiles and the resulting tumourigenic phenotype. Enhancers are transcriptional regulatory elements that are key to our understanding of how this rewiring of gene expression is achieved in cancer cells. Here, we have harnessed the power of RNA-seq data from hundreds of patients with oesophageal adenocarcinoma (OAC) or its precursor state Barrett's oesophagus coupled with open chromatin maps to identify potential enhancer RNAs and their associated enhancer regions in this cancer. We identify ~1000 OAC-specific enhancers and use these data to uncover new cellular pathways that are operational in OAC. Among these are enhancers for JUP, MYBL2, and CCNE1, and we show that their activity is required for cancer cell viability. We also demonstrate the clinical utility of our dataset for identifying disease stage and patient prognosis. Our data therefore identify an important set of regulatory elements that enhance our molecular understanding of OAC and point to potential new therapeutic directions.

Comparative transcriptome analyses reveal insights into catkin bloom patterns in pecan protogynous and protandrous cultivars

In perennial plants such as pecan, once reproductive maturity is attained, there are genetic switches that are regulated and required for flower development year after year. Pecan trees are heterodichogamous with both pistillate and staminate flowers produced on the same tree. Therefore, defining genes exclusively responsible for pistillate inflorescence and staminate inflorescence (catkin) initiation is challenging at best. To understand these genetic switches and their timing, this study analyzed catkin bloom and gene expression of lateral buds collected from a protogynous (Wichita) and a protandrous (Western) pecan cultivar in summer, autumn and spring. Our data showed that pistillate flowers in the current season on the same shoot negatively impacted catkin production on the protogynous 'Wichita' cultivar. Whereas fruit production the previous year on 'Wichita' had a positive effect on catkin production on the same shoot the following year. However, fruiting the previous year nor current year pistillate flower production had no significant effect on catkin production on 'Western' (protandrous cultivar) cultivar. The RNA-Seq results present more significant differences between the fruiting and non-fruiting shoots of the 'Wichita' cultivar compared to the 'Western' cultivar, revealing the genetic signals likely responsible for catkin production. Our data presented here, indicates the genes showing expression for the initiation of both types of flowers the season before bloom.

The ChIP-Hub Resource: Toward plantEncode

Recent advances in sequencing technologies lead to the generation of an enormous amount of regulome and epigenome data in a variety of plant species. However, a comprehensive standardized resource is so far not available. In this chapter, we present ChIP-Hub, an integrative platform that has been developed based on the ENCODE standards by collecting and reanalyzing regulatory genomic datasets from 41 plant species. The ChIP-hub website is introduced in this chapter, including information on detailed steps of searching, data download, and online analyses, which facilitates users to explore ChIP-Hub. We also provide a cross-species comparison of chromatin accessibility information that gives a thorough view of evolutionary regulatory networks in plants.

Floral Induction Systems for the Study of Arabidopsis Flower Development

Assessing the molecular changes that occur over the course of flower development is hampered by difficulties in isolating sufficient amounts of floral tissue at specific developmental stages. This is especially problematic when investigating molecular events at early stages of Arabidopsis flower development, as floral buds are minute and are initiated sequentially so that a single flower on an inflorescence is at a given developmental stage. Moreover, young floral buds are hidden by older flowers, which presents an additional challenge for dissection. To circumvent these issues, floral induction systems that allow the simultaneous induction of a large number of flowers on the inflorescence of a single plant were developed. To allow the plant community to avail of the full benefits of these systems, we address some common problems that can be encountered when growing these plants and collecting floral buds for analysis.

Systematic annotation of conservation states provides insights into regulatory regions in rice

Plant genomes contain a large fraction of noncoding sequences. The discovery and annotation of conserved noncoding sequences (CNSs) in plants is an ongoing challenge. Here we report the application of comparative genomics to systematically identify CNSs in 50 well-annotated Gramineae genomes using rice (Oryza sativa) as the reference. We conduct multiple-way whole-genome alignments to the rice genome. The rice genome is annotated as 20 conservation states (CSs) at single-nucleotide resolution using a multivariate hidden Markov model (ConsHMM) based on the multiple-genome alignments. Different states show distinct enrichments for various genomic features, and the conservation scores of CSs are highly correlated with the level of associated chromatin accessibility. We find that at least 33.5% of the rice genome is highly under selection, with more than 70% of the sequence lying outside of coding regions. A catalog of 855,366 regulatory CNSs is generated, and they significantly overlapped with putative active regulatory elements such as promoters, enhancers, and transcription factor binding sites. Collectively, our study provides a resource for elucidating functional noncoding regions of the rice genome and an evolutionary aspect of regulatory sequences in higher plants.

Comprehensive characterization of three classes of Arabidopsis SWI/SNF chromatin remodelling complexes

Although SWI/SNF chromatin remodelling complexes are known to regulate diverse biological functions in plants, the classification, compositions and functional mechanisms of the complexes remain to be determined. Here we comprehensively characterized SWI/SNF complexes by affinity purification and mass spectrometry in Arabidopsis thaliana, and found three classes of SWI/SNF complexes, which we termed BAS, SAS and MAS (BRM-, SYD- and MINU1/2-associated SWI/SNF complexes). By investigating multiple developmental phenotypes of SWI/SNF mutants, we found that three classes of SWI/SNF complexes have both overlapping and specific functions in regulating development. To investigate how the three classes of SWI/SNF complexes differentially regulate development, we mapped different SWI/SNF components on chromatin at the whole-genome level and determined their effects on chromatin accessibility. While all three classes of SWI/SNF complexes regulate chromatin accessibility at proximal promoter regions, SAS is a major SWI/SNF complex that is responsible for mediating chromatin accessibility at distal promoter regions and intergenic regions. Histone modifications are related to both the association of SWI/SNF complexes with chromatin and the SWI/SNF-dependent chromatin accessibility. Three classes of SWI/SNF-dependent accessibility may enable different sets of transcription factors to access chromatin. These findings lay a foundation for further investigation of the function of three classes of SWI/SNF complexes in plants.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

The chromatin accessibility landscape of pistils and anthers in rice

Transcription activation is tightly associated with the openness of chromatin and allows direct contact between transcriptional regulators and their targeted DNA for gene expression. However, there are limited studies on the annotation of open chromatin regions (OCRs) in rice (Oryza sativa), especially those in reproductive organs. Here, we characterized OCRs in rice pistils and anthers with an assay for transposase-accessible chromatin using sequencing. Despite a large overlap, we found more OCRs in pistils than in anthers. These OCRs were enriched in gene transcription start sites (TSSs) and showed tight associations with gene expression. Transcription factor (TF) binding motifs were enriched at these OCRs as validated by TF chromatin immunoprecipitation followed by sequencing. Pistil-specific OCRs provided potential regulatory networks by binding directly to the targets, indicating that pistil-specific OCRs may be indicators of cis-regulatory elements in regulating pistil development, which are absent in anthers. We also found that open chromatin of pistils and anthers responded differently to low temperature (LT). These data offer a comprehensive overview of OCRs regulating reproductive organ development and LT responses in rice.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Dynamic enhancer transcription associates with reprogramming of immune genes during pattern triggered immunity in Arabidopsis

BACKGROUND

Enhancers are cis-regulatory elements present in eukaryote genomes, which constitute indispensable determinants of gene regulation by governing the spatiotemporal and quantitative expression dynamics of target genes, and are involved in multiple life processes, for instance during development and disease states. The importance of enhancer activity has additionally been highlighted for immune responses in animals and plants; however, the dynamics of enhancer activities and molecular functions in plant innate immunity are largely unknown. Here, we investigated the involvement of distal enhancers in early innate immunity in Arabidopsis thaliana.

RESULTS

A group of putative distal enhancers producing low-abundance transcripts either unidirectionally or bidirectionally are identified. We show that enhancer transcripts are dynamically modulated in plant immunity triggered by microbe-associated molecular patterns and are strongly correlated with open chromatin, low levels of methylated DNA, and increases in RNA polymerase II targeting and acetylated histone marks. Dynamic enhancer transcription is correlated with target early immune gene expression patterns. Cis motifs that are bound by immune-related transcription factors, such as WRKYs and SARD1, are highly enriched within upregulated enhancers. Moreover, a subset of core pattern-induced enhancers are upregulated by multiple patterns from diverse pathogens. The expression dynamics of putative immunity-related enhancers and the importance of WRKY binding motifs for enhancer function were also validated.

CONCLUSIONS

Our study demonstrates the general occurrence of enhancer transcription in plants and provides novel information on the distal regulatory landscape during early plant innate immunity, providing new insights into immune gene regulation and ultimately improving the mechanistic understanding of the plant immune system.

ChIP-Hub provides an integrative platform for exploring plant regulome

Plant genomes encode a complex and evolutionary diverse regulatory grammar that forms the basis for most life on earth. A wealth of regulome and epigenome data have been generated in various plant species, but no common, standardized resource is available so far for biologists. Here, we present ChIP-Hub, an integrative web-based platform in the ENCODE standards that bundles >10,000 publicly available datasets reanalyzed from >40 plant species, allowing visualization and meta-analysis. We manually curate the datasets through assessing ~540 original publications and comprehensively evaluate their data quality. As a proof of concept, we extensively survey the co-association of different regulators and construct a hierarchical regulatory network under a broad developmental context. Furthermore, we show how our annotation allows to investigate the dynamic activity of tissue-specific regulatory elements (promoters and enhancers) and their underlying sequence grammar. Finally, we analyze the function and conservation of tissue-specific promoters, enhancers and chromatin states using comparative genomics approaches. Taken together, the ChIP-Hub platform and the analysis results provide rich resources for deep exploration of plant ENCODE. ChIP-Hub is available at https://biobigdata.nju.edu.cn/ChIPHub/ .

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

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Mapping open chromatin by ATAC-seq in bread wheat

Gene transcription is largely regulated by cis-regulatory elements. Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is an emerging technology that can accurately map cis-regulatory elements in animals and plants. However, the presence of cell walls and chloroplasts in plants hinders the extraction of high-quality nuclei, thereby affects the quality of ATAC-seq data. Meanwhile, it is tricky to perform ATAC-seq with different tissue types, especially for those with limited size and amount. Moreover, with rapid growth of ATAC-seq datasets from plants, powerful and easy-to-use data analysis pipelines for ATAC-seq, especially for wheat is lacking. Here, we provided an all-in-one solution for mapping open chromatin in wheat including both experimental and data analysis procedure. We efficiently obtained nuclei with less cell debris from various wheat tissues. High-quality ATAC-seq data from young spike and ovary, which are hard to harvest were generated. We determined that the saturation sequencing depth of wheat ATAC-seq is about 16 Gb. Particularly, we developed a powerful and easy-to-use online pipeline to analyze the wheat ATAC-seq data and this pipeline can be easily extended to other plant species. The method developed here will facilitate plant regulatory genome study not only for wheat but also for other plant species.

Building a Mammalian Retina: An Eye on Chromatin Structure

Regulation of gene expression by chromatin structure has been under intensive investigation, establishing nuclear organization and genome architecture as a potent and effective means of regulating developmental processes. The substantial growth in our knowledge of the molecular mechanisms underlying retinogenesis has been powered by several genome-wide based tools that mapped chromatin organization at multiple cellular and biochemical levels. Studies profiling the retinal epigenome and transcriptome have allowed the systematic annotation of putative cis-regulatory elements associated with transcriptional programs that drive retinal neural differentiation, laying the groundwork to understand spatiotemporal retinal gene regulation at a mechanistic level. In this review, we outline recent advances in our understanding of the chromatin architecture in the mammalian retina during development and disease. We focus on the emerging roles of non-coding regulatory elements in controlling retinal cell-type specific transcriptional programs, and discuss potential implications in untangling the etiology of eye-related disorders.

Transcriptional regulation in plants: Using omics data to crack the cis-regulatory code

Innovative omics technologies, advanced bioinformatics, and machine learning methods are rapidly becoming integral tools for plant functional genomics, with tremendous recent advances made in this field. In transcriptional regulation, an initial lag in the accumulation of plant omics data relative to that of animals stimulated the development of computational methods capable of extracting maximum information from the available data sets. Recent comprehensive studies of transcription factor-binding profiles in Arabidopsis and maize and the accumulation of uniformly processed omics data in public databases have brought plant biologists into the big leagues, with many cutting-edge methods available. Here, we summarize the state-of-the-art bioinformatics approaches used to predict or infer the cis-regulatory code behind transcriptional gene regulation, focusing on their plant research applications.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana

Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

Chicken C/EBPζ gene: Expression profiles, association analysis, and identification of functional variants for abdominal fat

CCAAT enhancer binding protein ζ (C/EBPζ) plays an important role in adipose proliferation and differentiation in humans. However, very little is known about the effect of C/EBPζ on the growth and development of adipose tissues in domesticated animals. The present study attempted to investigate the mRNA expression profiles of chicken C/EBPζ in a variety of tissues; analyze the association of its variants with abdominal fat; and identify the functional variants for abdominal fat. The tissue expression profiles revealed that C/EBPζ was highly expressed in 19 tissues obtained from broilers. The expression level of C/EBPζ in fat broilers was significantly lower than that in lean broilers in the duodenum, ileum, cecum, kidney, pectoral muscle, and liver (P < 0.05). Among 170 polymorphic loci of C/EBPζ, 9 single nucleotide polymorphisms (SNPs) demonstrated a significant association with chicken abdominal fat traits (P < 0.05) as well as significant discrepancies in their allelic frequencies between fat and lean birds. Particularly, only C/EBPζ g.7085A>C exhibited significant correlation with abdominal fat traits (P < 0.00015) using the Bonferroni method. The results revealed that, in preadipocyte immortalized cells (ICPI), the luciferase activity of the A allele of g.7085A>C locus was remarkably stronger than that of the C allele (P < 0.05). In silico analysis showed that g.7085A>C locus was located in the binding region of the transcription factor SOX5, which possesses the ability to transform C/EBPζ transcription efficiency through binding with SOX5. In summary, the data obtained from this study suggested that C/EBPζ is a potential candidate gene responsible for abdominal fat deposition in chicken and that g.7085A>C is a functional SNP that can be promisingly leveraged for marker assisted selection (MAS) in future chicken breeding programs.

H3.1K27me1 maintains transcriptional silencing and genome stability by preventing GCN5-mediated histone acetylation

Epigenetic mechanisms play diverse roles in the regulation of genome stability in eukaryotes. In Arabidopsis thaliana, genome stability is maintained during DNA replication by the H3.1K27 methyltransferases ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) and ATXR6, which catalyze the deposition of K27me1 on replication-dependent H3.1 variants. The loss of H3.1K27me1 in atxr5 atxr6 double mutants leads to heterochromatin defects, including transcriptional de-repression and genomic instability, but the molecular mechanisms involved remain largely unknown. In this study, we identified the transcriptional co-activator and conserved histone acetyltransferase GCN5 as a mediator of transcriptional de-repression and genomic instability in the absence of H3.1K27me1. GCN5 is part of a SAGA-like complex in plants that requires the GCN5-interacting protein ADA2b and the chromatin remodeler CHR6 to mediate the heterochromatic defects in atxr5 atxr6 mutants. Our results also indicate that Arabidopsis GCN5 acetylates multiple lysine residues on H3.1 variants, but H3.1K27 and H3.1K36 play essential functions in inducing genomic instability in the absence of H3.1K27me1. Finally, we show that H3.1K36 acetylation by GCN5 is negatively regulated by H3.1K27me1 in vitro. Overall, this work reveals a key molecular role for H3.1K27me1 in maintaining transcriptional silencing and genome stability in heterochromatin by restricting GCN5-mediated histone acetylation in plants.