Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape.

The cistrome is the complete set of transcription factor (TF) binding sites (cis-elements) in an organism, while an epicistrome incorporates tissue-specific DNA chemical modifications and TF-specific chemical sensitivities into these binding profiles. Robust methods to construct comprehensive cistrome and epicistrome maps are critical for elucidating complex transcriptional networks that underlie growth, behavior, and disease. Here, we describe DNA affinity purification sequencing (DAP-seq), a high-throughput TF binding site discovery method that interrogates genomic DNA with in-vitro-expressed TFs. Using DAP-seq, we defined the Arabidopsis cistrome by resolving motifs and peaks for 529 TFs. Because genomic DNA used in DAP-seq retains 5-methylcytosines, we determined that >75% (248/327) of Arabidopsis TFs surveyed were methylation sensitive, a property that strongly impacts the epicistrome landscape. DAP-seq datasets also yielded insight into the biology and binding site architecture of numerous TFs, demonstrating the value of DAP-seq for cost-effective cistromic and epicistromic annotation in any organism.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

The retrograde signaling regulator ANAC017 recruits the MKK9-MPK3/6, ethylene, and auxin signaling pathways to balance mitochondrial dysfunction with growth.

In plant cells, mitochondria are ideally positioned to sense and balance changes in energy metabolism in response to changing environmental conditions. Retrograde signaling from mitochondria to the nucleus is crucial for adjusting the required transcriptional responses. We show that ANAC017, the master regulator of mitochondrial stress, directly recruits a signaling cascade involving the plant hormones ethylene and auxin as well as the MAP KINASE KINASE (MKK) 9-MAP KINASE (MPK) 3/6 pathway in Arabidopsis thaliana. Chromatin immunoprecipitation followed by sequencing and overexpression demonstrated that ANAC017 directly regulates several genes of the ethylene and auxin pathways, including MKK9, 1-AMINO-CYCLOPROPANE-1-CARBOXYLATE SYNTHASE 2, and YUCCA 5, in addition to genes encoding transcription factors regulating plant growth and stress responses such as BASIC REGION/LEUCINE ZIPPER MOTIF (bZIP) 60, bZIP53, ANAC081/ATAF2, and RADICAL-INDUCED CELL DEATH1. A time-resolved RNA-seq experiment established that ethylene signaling precedes the stimulation of auxin signaling in the mitochondrial stress response, with a large part of the transcriptional regulation dependent on ETHYLENE-INSENSITIVE 3. These results were confirmed by mutant analyses. Our findings identify the molecular components controlled by ANAC017, which integrates the primary stress responses to mitochondrial dysfunction with whole plant growth via the activation of regulatory and partly antagonistic feedback loops.

Spatially resolved transcriptomic analysis of the germinating barley grain.

Seeds are a vital source of calories for humans and a unique stage in the life cycle of flowering plants. During seed germination, the embryo undergoes major developmental transitions to become a seedling. Studying gene expression in individual seed cell types has been challenging due to the lack of spatial information or low throughput of existing methods. To overcome these limitations, a spatial transcriptomics workflow was developed for germinating barley grain. This approach enabled high-throughput analysis of spatial gene expression, revealing specific spatial expression patterns of various functional gene categories at a sub-tissue level. This study revealed over 14 000 genes differentially regulated during the first 24 h after imbibition. Individual genes, such as the aquaporin gene family, starch degradation, cell wall modification, transport processes, ribosomal proteins and transcription factors, were found to have specific spatial expression patterns over time. Using spatial autocorrelation algorithms, we identified auxin transport genes that had increasingly focused expression within subdomains of the embryo over time, suggesting their role in establishing the embryo axis. Overall, our study provides an unprecedented spatially resolved cellular map for barley germination and identifies specific functional genomics targets to better understand cellular restricted processes during germination. The data can be viewed at https://spatial.latrobe.edu.au/.

Spatiotemporal deposition of cell wall polysaccharides in oat endosperm during grain development.

Oat (Avena sativa) is a cereal crop whose grains are rich in (1,3;1,4)-ß-D-glucan (mixed-linkage glucan or MLG), a soluble dietary fiber. In our study, we analyzed oat endosperm development in 2 Canadian varieties with differing MLG content and nutritional value. We confirmed that oat undergoes a nuclear type of endosperm development but with a shorter cellularization phase than barley (Hordeum vulgare). Callose and cellulose were the first polysaccharides to be detected in the early anticlinal cell walls at 11 days postemergence (DPE) of the panicle. Other polysaccharides such as heteromannan and homogalacturonan were deposited early in cellularization around 12 DPE after the first periclinal walls are laid down. In contrast to barley, heteroxylan deposition coincided with completion of cellularization and was detected from 14 DPE but was only detectable after demasking. Notably, MLG was the last polysaccharide to be laid down at 18 DPE within the differentiation phase, rather than during cellularization. In addition, differences in the spatiotemporal patterning of MLG were also observed between the 2 varieties. The lower MLG-containing cultivar AC Morgan (3.5% w/w groats) was marked by the presence of a discontinuous pattern of MLG labeling, while labeling in the same walls in CDC Morrison (5.6% w/w groats) was mostly even and continuous. RNA-sequencing analysis revealed higher transcript levels of multiple MLG biosynthetic cellulose synthase-like F (CSLF) and CSLH genes during grain development in CDC Morrison compared with AC Morgan that likely contributes to the increased abundance of MLG at maturity in CDC Morrison. CDC Morrison was also observed to have smaller endosperm cells with thicker walls than AC Morgan from cellularization onwards, suggesting the processes controlling cell size and shape are established early in development. This study has highlighted that the molecular processes influencing MLG content and deposition are more complex than previously imagined.

Genomic and Cell-Specific Regulation of Benzylisoquinoline Alkaloid Biosynthesis in Opium Poppy.

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine and papaverine. Most enzymes involved in benzylisoquinoline alkaloids (BIAs) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.

Sink strength, nutrient allocation, cannabinoid yield and associated transcript profiles vary in two drug-type Cannabis chemovars.

Cannabis sativa L. is one of the oldest domesticated crops. Hemp-type cultivars, which predominantly produce non-intoxicating cannabidiol (CBD), have been selected for their fast growth, seed, and fibre production, while drug-type chemovars were bred for high accumulation of tetrahydrocannabinol (THC). We investigated how the generation of CBD-dominant chemovars by introgression of hemp- into drug-type Cannabis impacted plant performance. The THC-dominant chemovar showed superior sink strength, higher flower biomass and demand-driven control of nutrient uptake. By contrast, the CBD-dominant chemovar hyperaccumulated phosphate in sink organs leading to reduced carbon and nitrogen assimilation in leaves, which limited flower biomass and cannabinoid yield. RNA-seq analyses determined organ- and chemovar-specific differences in expression of genes associated with nitrate and phosphate homeostasis as well as growth-regulating transcription factors that were correlated with measured traits. Among these were genes positively selected for during Cannabis domestication encoding an inhibitor of the phosphate starvation response SPX DOMAIN GENE3, nitrate reductase and two nitrate transporters. Altered nutrient sensing, acquisition or distribution are likely a consequence of adaption to growth on marginal, low-nutrient input lands in hemp. Our data provide evidence that such ancestral traits may become detrimental for female flower development and consequently overall CBD yield in protected cropping environments.

RNA-seq analysis of laser microdissected Arabidopsis thaliana leaf epidermis, mesophyll and vasculature defines tissue-specific transcriptional responses to multiple stress treatments.

Acclimation of plants to adverse conditions requires the coordination of gene expression and signalling pathways between tissues and cell types. As the energy and carbon capturing organs, leaves are significantly affected by abiotic and biotic stresses. However, tissue- or cell type-specific analyses of stress responses have focussed on the Arabidopsis root. Here, we comparatively explore the transcriptomes of three leaf tissues (epidermis, mesophyll, vasculature) after induction of diverse stress pathways by chemical stimuli (antimycin A, 3-amino-1,2,4-triazole, methyl viologen, salicylic acid) and ultraviolet light in Arabidopsis using laser capture microdissection followed by RNA sequencing. Stimulation of stress pathways caused an overall reduction in the number of genes expressed in a tissue-specific manner, though a small subset gained or changed their tissue specificity. We find no evidence of a common stress response, with only a few genes consistently responsive to two or more treatments in the analysed tissues. However, differentially expressed genes overlap between tissues for individual treatments. A focussed analysis provided evidence for an interaction of auxin and ethylene that mediates retrograde signalling during mitochondrial dysfunction specifically in the epidermis, and a gene regulatory network defined the hierarchy of interactions. Taken together, we have generated an extensive reference dataset that will be valuable for future experiments analysing transcriptional responses on a tissue or single-cell level. Our results will enable the tailoring of the tissue-specific engineering of stress-tolerant plants.

A MYC2/MYC3/MYC4-dependent transcription factor network regulates water spray-responsive gene expression and jasmonate levels.

Mechanical stimuli, such as wind, rain, and touch affect plant development, growth, pest resistance, and ultimately reproductive success. Using water spray to simulate rain, we demonstrate that jasmonic acid (JA) signaling plays a key role in early gene-expression changes, well before it leads to developmental changes in flowering and plant architecture. The JA-activated transcription factors MYC2/MYC3/MYC4 modulate transiently induced expression of 266 genes, most of which peak within 30 min, and control 52% of genes induced >100-fold. Chromatin immunoprecipitation-sequencing analysis indicates that MYC2 dynamically binds >1,300 promoters and trans-activation assays show that MYC2 activates these promoters. By mining our multiomic datasets, we identified a core MYC2/MYC3/MYC4-dependent "regulon" of 82 genes containing many previously unknown MYC2 targets, including transcription factors bHLH19 and ERF109 bHLH19 can in turn directly activate the ORA47 promoter, indicating that MYC2/MYC3/MYC4 initiate a hierarchical network of downstream transcription factors. Finally, we also reveal that rapid water spray-induced accumulation of JA and JA-isoleucine is directly controlled by MYC2/MYC3/MYC4 through a positive amplification loop that regulates JA-biosynthesis genes.

Temporal tissue-specific regulation of transcriptomes during barley (Hordeum vulgare) seed germination.

The distinct functions of individual cell types require cells to express specific sets of genes. The germinating seed is an excellent model to study genome regulation between cell types since the majority of the transcriptome is differentially expressed in a short period, beginning from a uniform, metabolically inactive state. In this study, we applied laser-capture microdissection RNA-sequencing to small numbers of cells from the plumule, radicle tip and scutellum of germinating barley seeds every 8 h, over a 48 h time course. Tissue-specific gene expression was notably common; 25% (910) of differentially expressed transcripts in plumule, 34% (1876) in radicle tip and 41% (2562) in scutellum were exclusive to that organ. We also determined that tissue-specific storage of transcripts occurs during seed development and maturation. Co-expression of genes had strong spatiotemporal structure, with most co-expression occurring within one organ and at a subset of specific time points during germination. Overlapping and distinct enrichment of functional categories were observed in the tissue-specific profiles. We identified candidate transcription factors amongst these that may be regulators of spatiotemporal gene expression programs. Our findings contribute to the broader goal of generating an integrative model that describes the structure and function of individual cells within seeds during germination.

The JA-pathway MYC transcription factors regulate photomorphogenic responses by targeting HY5 gene expression.

Jasmonates are key regulators of the balance between defence and growth in plants. However, the molecular mechanisms by which activation of defence reduces growth are not yet fully understood. Here, we analyze the role of MYC transcription factors (TFs) and jasmonic acid (JA) in photomorphogenic growth. We found that multiple myc mutants share light-associated phenotypes with mutants of the phytochrome B photoreceptor, such as delayed seed germination in the dark and long hypocotyl growth. Overexpression of MYC2 in a phyB background partially suppressed its long hypocotyl phenotype. Transcriptomic analysis of multiple myc mutants confirmed that MYCs are required for full expression of red (R) light-regulated genes, including the master regulator HY5. ChIP-seq analyses revealed that MYC2 and MYC3 bind directly to the promoter of HY5 and that HY5 gene expression and protein levels are compromised in multiple myc mutants. Altogether, our results pinpoint MYCs as photomorphogenic TFs that control phytochrome responses by activating HY5 expression. This has important implications in understanding the trade-off between growth and defence as the same TFs that activate defence responses are photomorphogenic growth regulators.

Recent advances in Cannabis sativa genomics research.

Cannabis (Cannabis sativa L.) is one of the oldest cultivated plants purported to have unique medicinal properties. However, scientific research of cannabis has been restricted by the Single Convention on Narcotic Drugs of 1961, an international treaty that prohibits the production and supply of narcotic drugs except under license. Legislation governing cannabis cultivation for research, medicinal and even recreational purposes has been relaxed recently in certain jurisdictions. As a result, there is now potential to accelerate cultivar development of this multi-use and potentially medically useful plant species by application of modern genomics technologies. Whilst genomics has been pivotal to our understanding of the basic biology and molecular mechanisms controlling key traits in several crop species, much work is needed for cannabis. In this review we provide a comprehensive summary of key cannabis genomics resources and their applications. We also discuss prospective applications of existing and emerging genomics technologies for accelerating the genetic improvement of cannabis.

ANAC017 Coordinates Organellar Functions and Stress Responses by Reprogramming Retrograde Signaling.

Mitochondria adjust their activities in response to external and internal stimuli to optimize growth via the mitochondrial retrograde response signaling pathway. The Arabidopsis (Arabidopsis thaliana) NAC domain transcription factor ANAC017 has previously been identified as a regulator of the mitochondrial retrograde response. We show here that overexpression of ANAC017 in Arabidopsis leads to growth retardation, altered leaf development with decreased cell size and viability, and early leaf senescence. RNA sequencing analyses revealed that increased ANAC017 expression leads to higher expression of genes related to mitochondrial stress, cell death/autophagy, and leaf senescence under nonlimiting growth conditions as well as extensive repression of chloroplast function. Gene regulatory network analysis indicated that a complex hierarchy of transcription factors exists downstream of ANAC017. These involve a set of up-regulated ANAC and WRKY transcription factors associated with organellar signaling and senescence. The network also includes a number of ethylene- and gibberellic acid-related transcription factors with established functions in stress responses and growth regulation, which down-regulate their target genes. A number of BASIC LEUCINE-ZIPPER MOTIF transcription factors involved in the endoplasmic reticulum unfolded protein response or balancing of energy homeostasis via the SNF1-RELATED PROTEIN KINASE1 were also down-regulated by ANAC017 overexpression. Our results show that the endoplasmic reticulum membrane tethering of the constitutively expressed ANAC017, and its controlled release, are crucial to fine-tune a fast reactive but potentially harmful signaling cascade. Thus, ANAC017 is a master regulator of cellular responses with mitochondria acting as central sensors.

Analysis of Spatio-Temporal Transcriptome Profiles of Soybean (Glycine max) Tissues during Early Seed Development.

Soybean (Glycine max) is an important crop providing oil and protein for both human and animal consumption. Knowing which biological processes take place in specific tissues in a temporal manner will enable directed breeding or synthetic approaches to improve seed quantity and quality. We analyzed a genome-wide transcriptome dataset from embryo, endosperm, endothelium, epidermis, hilum, outer and inner integument and suspensor at the global, heart and cotyledon stages of soybean seed development. The tissue specificity of gene expression was greater than stage specificity, and only three genes were differentially expressed in all seed tissues. Tissues had both unique and shared enriched functional categories of tissue-specifically expressed genes associated with them. Strong spatio-temporal correlation in gene expression was identified using weighted gene co-expression network analysis, with the most co-expression occurring in one seed tissue. Transcription factors with distinct spatiotemporal gene expression programs in each seed tissue were identified as candidate regulators of expression within those tissues. Gene ontology (GO) enrichment of orthogroup clusters revealed the conserved functions and unique roles of orthogroups with similar and contrasting expression patterns in transcript abundance between soybean and Arabidopsis during embryo proper and endosperm development. Key regulators in each seed tissue and hub genes connecting those networks were characterized by constructing gene regulatory networks. Our findings provide an important resource for describing the structure and function of individual soybean seed compartments during early seed development.

Mobile small RNAs regulate genome-wide DNA methylation.

RNA silencing at the transcriptional and posttranscriptional levels regulates endogenous gene expression, controls invading transposable elements (TEs), and protects the cell against viruses. Key components of the mechanism are small RNAs (sRNAs) of 21-24 nt that guide the silencing machinery to their nucleic acid targets in a nucleotide sequence-specific manner. Transcriptional gene silencing is associated with 24-nt sRNAs and RNA-directed DNA methylation (RdDM) at cytosine residues in three DNA sequence contexts (CG, CHG, and CHH). We previously demonstrated that 24-nt sRNAs are mobile from shoot to root in Arabidopsis thaliana and confirmed that they mediate DNA methylation at three sites in recipient cells. In this study, we extend this finding by demonstrating that RdDM of thousands of loci in root tissues is dependent upon mobile sRNAs from the shoot and that mobile sRNA-dependent DNA methylation occurs predominantly in non-CG contexts. Mobile sRNA-dependent non-CG methylation is largely dependent on the DOMAINS REARRANGED METHYLTRANSFERASES 1/2 (DRM1/DRM2) RdDM pathway but is independent of the CHROMOMETHYLASE (CMT)2/3 DNA methyltransferases. Specific superfamilies of TEs, including those typically found in gene-rich euchromatic regions, lose DNA methylation in a mutant lacking 22- to 24-nt sRNAs (dicer-like 2, 3, 4 triple mutant). Transcriptome analyses identified a small number of genes whose expression in roots is associated with mobile sRNAs and connected to DNA methylation directly or indirectly. Finally, we demonstrate that sRNAs from shoots of one accession move across a graft union and target DNA methylation de novo at normally unmethylated sites in the genomes of root cells from a different accession.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species.

BACKGROUND: The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops. DESCRIPTION: AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications. CONCLUSION: Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Regulation of genome-wide DNA methylation by mobile small RNAs.

Contents Summary 540 I. Introduction 540 II. There are different types of sRNA mobility 541 III. Mechanisms of sRNA movement 541 IV. Long-distance, shoot-root, mobile siRNAs influence DNA methylation in recipient tissues 541 V. Classes of interactions between shoot-root mobile siRNAs and DNA methylation 542 VI. Loci targeted directly and indirectly by shoot-root mobile siRNAs are associated with different histone modifications 543 VII. Is mobile siRNA-regulated DNA methylation important in specific tissues or under specific conditions? 543 VIII. Mobile sRNAs can be used to modify plant traits 544 IX. Conclusions 544 Acknowledgements 544 References 544 SUMMARY: RNA-directed DNA methylation (RdDM) at cytosine residues regulates gene expression, silences transposable elements and influences genome stability. The mechanisms responsible for RdDM are guided to target loci by small RNAs (sRNAs) that can move within plants cell to cell and long distance. Here we discuss recent advances in the understanding of interactions between mobile sRNAs and DNA methylation. We describe the mechanisms of sRNA movement, the differences between known classes of mobile sRNA-DNA methylation interactions and the limits of current knowledge. Finally, we discuss potential applications of mobile sRNAs in modifying plant traits.

JAZ2 controls stomata dynamics during bacterial invasion.

Coronatine (COR) facilitates entry of bacteria into the plant apoplast by stimulating stomata opening. COR-induced signaling events at stomata remain unclear. We found that the COR and jasmonate isoleucine (JA-Ile) co-receptor JAZ2 is constitutively expressed in guard cells and modulates stomatal dynamics during bacterial invasion We analyzed tissue expression patterns of AtJAZ genes and measured stomata opening and pathogen resistance in loss- and gain-of-function mutants. Arabidopsis jaz2 mutants are partially impaired in pathogen-induced stomatal closing and more susceptible to Pseudomonas. Gain-of-function mutations in JAZ2 prevent stomatal reopening by COR and are highly resistant to bacterial penetration. The JAZ2 targets MYC2, MYC3 and MYC4 directly regulate the expression of ANAC19, ANAC55 and ANAC72 to modulate stomata aperture. Due to the antagonistic interactions between the salicylic acid (SA) and JA defense pathways, efforts to increase resistance to biotrophs result in enhanced susceptibility to necrotrophs, and vice versa. Remarkably, dominant jaz2Δjas mutants are resistant to Pseudomonas syringae but retain unaltered resistance against necrotrophs. Our results demonstrate the existence of a COI1-JAZ2-MYC2,3,4-ANAC19,55,72 module responsible for the regulation of stomatal aperture that is hijacked by bacterial COR to promote infection. They also provide novel strategies for crop protection against biotrophs without compromising resistance to necrotrophs.