Transcriptomic signature of cancer cachexia by integration of machine learning, literature mining and meta-analysis

BACKGROUND

Cancer cachexia is a severe metabolic syndrome marked by skeletal muscle atrophy. A successful clinical intervention for cancer cachexia is currently lacking. The study of cachexia mechanisms is largely based on preclinical animal models and the availability of high-throughput transcriptomic datasets of cachectic mouse muscles is increasing through the extensive use of next generation sequencing technologies.

METHODS

Cachectic mouse muscle transcriptomic datasets of ten different studies were combined and mined by seven attribute weighting models, which analysed both categorical variables and numerical variables. The transcriptomic signature of cancer cachexia was identified by attribute weighting algorithms and was used to evaluate the performance of eleven pattern discovery models. The signature was employed to find the best combination of drugs (drug repurposing) for developing cancer cachexia treatment strategies, as well as to evaluate currently used cachexia drugs by literature mining.

RESULTS

Attribute weighting algorithms ranked 26 genes as the transcriptomic signature of muscle from mice with cancer cachexia. Deep Learning and Random Forest models performed better in differentiating cancer cachexia cases based on muscle transcriptomic data. Literature mining revealed that a combination of melatonin and infliximab has negative interactions with 2 key genes (Rorc and Fbxo32) upregulated in the transcriptomic signature of cancer cachexia in muscle.

CONCLUSIONS

The integration of machine learning, meta-analysis and literature mining was found to be an efficient approach to identifying a robust transcriptomic signature for cancer cachexia, with implications for improving clinical diagnosis and management of this condition.

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.

Transcription Factor Dynamics in Cross-Regulation of Plant Hormone Signaling Pathways

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

A Practical Guide to Inferring Multi-Omics Networks in Plant Systems

The inference of gene regulatory networks can reveal molecular connections underlying biological processes and improve our understanding of complex biological phenomena in plants. Many previous network studies have inferred networks using only one type of omics data, such as transcriptomics. However, given more recent work applying multi-omics integration in plant biology, such as combining (phospho)proteomics with transcriptomics, it may be advantageous to integrate multiple omics data types into a comprehensive network prediction. Here, we describe a state-of-the-art approach for integrating multi-omics data with gene regulatory network inference to describe signaling pathways and uncover novel regulators. We detail how to download and process transcriptomics and (phospho)proteomics data for network inference, using an example dataset from the plant hormone signaling field. We provide a step-by-step protocol for inference, visualization, and analysis of an integrative multi-omics network using currently available methods. This chapter serves as an accessible guide for novice and intermediate bioinformaticians to analyze their own datasets and reanalyze published work.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

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.

scCloudMine: A cloud-based app for visualization, comparison, and exploration of single-cell transcriptomic data

scCloudMine is a cloud-based application for visualization, comparison, and exploration of single-cell transcriptome data. It does not require an on-site, high-power computing server, installation, or associated expertise and expense. Users upload their own or publicly available scRNA-seq datasets after pre-processing for visualization using a web browser. The data can be viewed in two color modes-Cluster, representing cell identity, and Values, showing levels of expression-and data can be queried using keywords or gene identification number(s). Using the app to compare studies, we determined that some genes frequently used as cell-type markers are in fact study specific. The apparent cell-specific expression of PHO1;H3 differed between GFP-tagging and scRNA-seq studies. Some phosphate transporter genes were induced by protoplasting, but they retained cell specificity, suggesting that cell-specific responses to stress (i.e., protoplasting) can occur. Examination of the cell specificity of hormone response genes revealed that 132 hormone-responsive genes display restricted expression and that the jasmonate response gene TIFY8 is expressed in endodermal cells, in contrast to previous reports. It also appears that JAZ repressors have cell-type-specific functions. These features identified using scCloudMine highlight the need for resources to enable biological researchers to compare their datasets of interest under a variety of parameters. scCloudMine enables researchers to form new hypotheses and perform comparative studies and allows for the easy re-use of data from this emerging technology by a wide variety of users who may not have access or funding for high-performance on-site computing and support.

Applications of hyperspectral imaging in plant phenotyping

Our ability to interrogate and manipulate the genome far exceeds our capacity to measure the effects of genetic changes on plant traits. Much effort has been made recently by the plant science research community to address this imbalance. The responses of plants to environmental conditions can now be defined using a variety of imaging approaches. Hyperspectral imaging (HSI) has emerged as a promising approach to measure traits using a wide range of wavebands simultaneously in 3D to capture information in lab, glasshouse, or field settings. HSI has been applied to define abiotic, biotic, and quality traits for optimisation of crop management.

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.

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.

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.

Developmental normalization of phenomics data generated by high throughput plant phenotyping systems

BACKGROUND

Sowing time is commonly used as the temporal reference for Arabidopsis thaliana (Arabidopsis) experiments in high throughput plant phenotyping (HTPP) systems. This relies on the assumption that germination and seedling establishment are uniform across the population. However, individual seeds have different development trajectories even under uniform environmental conditions. This leads to increased variance in quantitative phenotyping approaches. We developed the Digital Adjustment of Plant Development (DAPD) normalization method. It normalizes time-series HTPP measurements by reference to an early developmental stage and in an automated manner. The timeline of each measurement series is shifted to a reference time. The normalization is determined by cross-correlation at multiple time points of the time-series measurements, which may include rosette area, leaf size, and number.

RESULTS

The DAPD method improved the accuracy of phenotyping measurements by decreasing the statistical dispersion of quantitative traits across a time-series. We applied DAPD to evaluate the relative growth rate in Arabidopsis plants and demonstrated that it improves uniformity in measurements, permitting a more informative comparison between individuals. Application of DAPD decreased variance of phenotyping measurements by up to 2.5 times compared to sowing-time normalization. The DAPD method also identified more outliers than any other central tendency technique applied to the non-normalized dataset.

CONCLUSIONS

DAPD is an effective method to control for temporal differences in development within plant phenotyping datasets. In principle, it can be applied to HTPP data from any species/trait combination for which a relevant developmental scale can be defined.

Laser-Capture Microdissection RNA-Sequencing for Spatial and Temporal Tissue-Specific Gene Expression Analysis in Plants

The development of a complex multicellular organism is governed by distinct cell types that have different transcriptional profiles. To identify transcriptional regulatory networks that govern developmental processes it is necessary to measure the spatial and temporal gene expression profiles of these individual cell types. Therefore, insight into the spatio-temporal control of gene expression is essential to gain understanding of how biological and developmental processes are regulated. Here, we describe a laser-capture microdissection (LCM) method to isolate small number of cells from three barley embryo organs over a time-course during germination followed by transcript profiling. The method consists of tissue fixation, tissue processing, paraffin embedding, sectioning, LCM and RNA extraction followed by real-time PCR or RNA-seq. This method has enabled us to obtain spatial and temporal profiles of seed organ transcriptomes from varying numbers of cells (tens to hundreds), providing much greater tissue-specificity than typical bulk-tissue analyses. From these data we were able to define and compare transcriptional regulatory networks as well as predict candidate regulatory transcription factors for individual tissues. The method should be applicable to other plant tissues with minimal optimization.

Publisher Correction: Integrated multi-omics framework of the plant response to jasmonic acid

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Integrated multi-omics framework of the plant response to jasmonic acid

Understanding the systems-level actions of transcriptional responses to hormones provides insight into how the genome is reprogrammed in response to environmental stimuli. Here, we investigated the signalling pathway of the hormone jasmonic acid (JA), which controls a plethora of critically important processes in plants and is orchestrated by the transcription factor MYC2 and its closest relatives in Arabidopsis thaliana. We generated an integrated framework of the response to JA, which spans from the activity of master and secondary regulatory transcription factors, through gene expression outputs and alternative splicing, to protein abundance changes, protein phosphorylation and chromatin remodelling. We integrated time-series transcriptome analysis with (phospho)proteomic data to reconstruct gene regulatory network models. These enabled us to predict previously unknown points of crosstalk of JA to other signalling pathways and to identify new components of the JA regulatory mechanism, which we validated through targeted mutant analysis. These results provide a comprehensive understanding of how a plant hormone remodels cellular functions and plant behaviour, the general principles of which provide a framework for analyses of cross-regulation between other hormone and stress signalling pathways.

Integrated multi-omics framework of the plant response to jasmonic acid

Understanding the systems-level actions of transcriptional responses to hormones provides insight into how the genome is reprogrammed in response to environmental stimuli. Here, we investigated the signalling pathway of the hormone jasmonic acid (JA), which controls a plethora of critically important processes in plants and is orchestrated by the transcription factor MYC2 and its closest relatives in Arabidopsis thaliana. We generated an integrated framework of the response to JA, which spans from the activity of master and secondary regulatory transcription factors, through gene expression outputs and alternative splicing, to protein abundance changes, protein phosphorylation and chromatin remodelling. We integrated time-series transcriptome analysis with (phospho)proteomic data to reconstruct gene regulatory network models. These enabled us to predict previously unknown points of crosstalk of JA to other signalling pathways and to identify new components of the JA regulatory mechanism, which we validated through targeted mutant analysis. These results provide a comprehensive understanding of how a plant hormone remodels cellular functions and plant behaviour, the general principles of which provide a framework for analyses of cross-regulation between other hormone and stress signalling pathways.

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.

Epigenetic silencing of a multifunctional plant stress regulator

The central regulator of the ethylene (ET) signaling pathway, which controls a plethora of developmental programs and responses to environmental cues in plants, is ETHYLENE-INSENSITIVE2 (EIN2). Here we identify a chromatin-dependent regulatory mechanism at EIN2 requiring two genes: ETHYLENE-INSENSITIVE6 (EIN6), which is a H3K27me3 demethylase also known as RELATIVE OF EARLY FLOWERING6 (REF6), and EIN6 ENHANCER (EEN), the Arabidopsis homolog of the yeast INO80 chromatin remodeling complex subunit IES6 (INO EIGHTY SUBUNIT). Strikingly, EIN6 (REF6) and the INO80 complex redundantly control the level and the localization of the repressive histone modification H3K27me3 and the histone variant H2A.Z at the 5’ untranslated region (5’UTR) intron of EIN2. Concomitant loss of EIN6 (REF6) and the INO80 complex shifts the chromatin landscape at EIN2 to a repressive state causing a dramatic reduction of EIN2 expression. These results uncover a unique type of chromatin regulation which safeguards the expression of an essential multifunctional plant stress regulator.

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.

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.

Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination

BACKGROUND

Seed germination involves progression from complete metabolic dormancy to a highly active, growing seedling. Many factors regulate germination and these interact extensively, forming a complex network of inputs that control the seed-to-seedling transition. Our understanding of the direct regulation of gene expression and the dynamic changes in the epigenome and small RNAs during germination is limited. The interactions between genome, transcriptome and epigenome must be revealed in order to identify the regulatory mechanisms that control seed germination.

RESULTS

We present an integrated analysis of high-resolution RNA sequencing, small RNA sequencing and MethylC sequencing over ten developmental time points in Arabidopsis thaliana seeds, finding extensive transcriptomic and epigenomic transformations associated with seed germination. We identify previously unannotated loci from which messenger RNAs are expressed transiently during germination and find widespread alternative splicing and divergent isoform abundance of genes involved in RNA processing and splicing. We generate the first dynamic transcription factor network model of germination, identifying known and novel regulatory factors. Expression of both microRNA and short interfering RNA loci changes significantly during germination, particularly between the seed and the post-germinative seedling. These are associated with changes in gene expression and large-scale demethylation observed towards the end of germination, as the epigenome transitions from an embryo-like to a vegetative seedling state.

CONCLUSIONS

This study reveals the complex dynamics and interactions of the transcriptome and epigenome during seed germination, including the extensive remodelling of the seed DNA methylome from an embryo-like to vegetative-like state during the seed-to-seedling transition. Data are available for exploration in a user-friendly browser at https://jbrowse.latrobe.edu.au/germination_epigenome .

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.

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.

Mobile small RNAs and their role in regulating cytosine methylation of DNA

Small (s)RNAs of 21 to 24 nucleotides are associated with RNA silencing and methylation of DNA cytosine residues. All sizes can move from cell-to-cell and long distance in plants, directing RNA silencing in destination cells. Twenty-four nucleotide sRNAs are the predominant long-distance mobile species. Thousands move from shoot to root, where they target methylation of transposable elements both directly and indirectly. We derive several classes of interaction between small RNAs and methylation and use these to explore the mechanisms of methylation and gene expression that associate with mobile sRNA signaling.

Salicylic acid treatment and expression of an RNA-dependent RNA polymerase 1 transgene inhibit lethal symptoms and meristem invasion during tobacco mosaic virus infection in Nicotiana benthamiana

BACKGROUND

Host RNA-dependent RNA polymerases (RDRs) 1 and 6 contribute to antiviral RNA silencing in plants. RDR6 is constitutively expressed and was previously shown to limit invasion of Nicotiana benthamiana meristem tissue by potato virus X and thereby inhibit disease development. RDR1 is inducible by salicylic acid (SA) and several other phytohormones. But although it contributes to basal resistance to tobacco mosaic virus (TMV) it is dispensable for SA-induced resistance in inoculated leaves. The laboratory accession of N. benthamiana is a natural rdr1 mutant and highly susceptible to TMV. However, TMV-induced symptoms are ameliorated in transgenic plants expressing Medicago truncatula RDR1.

RESULTS

In MtRDR1-transgenic N. benthamiana plants the spread of TMV expressing the green fluorescent protein (TMV.GFP) into upper, non-inoculated, leaves was not inhibited. However, in these plants exclusion of TMV.GFP from the apical meristem and adjacent stem tissue was greater than in control plants and this exclusion effect was enhanced by SA. TMV normally kills N. benthamiana plants but although MtRDR1-transgenic plants initially displayed virus-induced necrosis they subsequently recovered. Recovery from disease was markedly enhanced by SA treatment in MtRDR1-transgenic plants whereas in control plants SA delayed but did not prevent systemic necrosis and death. Following SA treatment of MtRDR1-transgenic plants, extractable RDR enzyme activity was increased and Western blot analysis of RDR extracts revealed a band cross-reacting with an antibody raised against MtRDR1. Expression of MtRDR1 in the transgenic N. benthamiana plants was driven by a constitutive 35S promoter derived from cauliflower mosaic virus, confirmed to be non-responsive to SA. This suggests that the effects of SA on MtRDR1 are exerted at a post-transcriptional level.

CONCLUSIONS

MtRDR1 inhibits severe symptom development by limiting spread of virus into the growing tips of infected plants. Thus, RDR1 may act in a similar fashion to RDR6. MtRDR1 and SA acted additively to further promote recovery from disease symptoms in MtRDR1-transgenic plants. Thus it is possible that SA promotes MtRDR1 activity and/or stability through post-transcriptional effects.

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.

Domains of the cucumber mosaic virus 2b silencing suppressor protein affecting inhibition of salicylic acid-induced resistance and priming of salicylic acid accumulation during infection

The cucumber mosaic virus (CMV) 2b silencing suppressor protein allows the virus to overcome resistance to replication and local movement in inoculated leaves of plants treated with salicylic acid (SA), a resistance-inducing plant hormone. In Arabidopsis thaliana plants systemically infected with CMV, the 2b protein also primes the induction of SA biosynthesis during this compatible interaction. We found that CMV infection of susceptible tobacco (Nicotiana tabacum) also induced SA accumulation. Utilization of mutant 2b proteins expressed during infection of tobacco showed that the N- and C-terminal domains, which had previously been implicated in regulation of symptom induction, were both required for subversion of SA-induced resistance, while all mutants tested except those affecting the putative phosphorylation domain had lost the ability to prime SA accumulation and expression of the SA-induced marker gene PR-1.

Interference with jasmonic acid-regulated gene expression is a general property of viral suppressors of RNA silencing but only partly explains virus-induced changes in plant-aphid interactions

The cucumber mosaic virus (CMV) 2b viral suppressor of RNA silencing (VSR) inhibits host responses to jasmonic acid (JA), a chemical signal regulating resistance to insects. Previous experiments with a CMV subgroup IA strain and its 2b gene deletion mutant suggested that VSRs might neutralize aphid (Myzus persicae) resistance by inhibiting JA-regulated gene expression. To further investigate this, we examined JA-regulated gene expression and aphid performance in Nicotiana benthamiana infected with Potato virus X, Potato virus Y, Tobacco mosaic virus and a subgroup II CMV strain, as well as in transgenic plants expressing corresponding VSRs (p25, HC-Pro, 126 kDa and 2b). All the viruses or their VSRs inhibited JA-induced gene expression. However, this did not always correlate with enhanced aphid performance. Thus, VSRs are not the sole viral determinants of virus-induced changes in host-aphid interactions and interference with JA-regulated gene expression cannot completely explain enhanced aphid performance on virus-infected plants.

A trio of viral proteins tunes aphid-plant interactions in Arabidopsis thaliana

Background

Virus-induced deterrence to aphid feeding is believed to promote plant virus transmission by encouraging migration of virus-bearing insects away from infected plants. We investigated the effects of infection by an aphid-transmitted virus, cucumber mosaic virus (CMV), on the interaction of Arabidopsis thaliana, one of the natural hosts for CMV, with Myzus persicae (common names: ‘peach-potato aphid’, ‘green peach aphid’).

Methodology/Principal Findings

Infection of Arabidopsis (ecotype Col-0) with CMV strain Fny (Fny-CMV) induced biosynthesis of the aphid feeding-deterrent 4-methoxy-indol-3-yl-methylglucosinolate (4MI3M). 4MI3M inhibited phloem ingestion by aphids and consequently discouraged aphid settling. The CMV 2b protein is a suppressor of antiviral RNA silencing, which has previously been implicated in altering plant-aphid interactions. Its presence in infected hosts enhances the accumulation of CMV and the other four viral proteins. Another viral gene product, the 2a protein (an RNA-dependent RNA polymerase), triggers defensive signaling, leading to increased 4MI3M accumulation. The 2b protein can inhibit ARGONAUTE1 (AGO1), a host factor that both positively-regulates 4MI3M biosynthesis and negatively-regulates accumulation of substance(s) toxic to aphids. However, the 1a replicase protein moderated 2b-mediated inhibition of AGO1, ensuring that aphids were deterred from feeding but not poisoned. The LS strain of CMV did not induce feeding deterrence in Arabidopsis ecotype Col-0.

Conclusions/Significance

Inhibition of AGO1 by the 2b protein could act as a booby trap since this will trigger antibiosis against aphids. However, for Fny-CMV the interplay of three viral proteins (1a, 2a and 2b) appears to balance the need of the virus to inhibit antiviral silencing, while inducing a mild resistance (antixenosis) that is thought to promote transmission. The strain-specific effects of CMV on Arabidopsis-aphid interactions, and differences between the effects of Fny-CMV on this plant and those seen previously in tobacco (inhibition of resistance to aphids) may have important epidemiological consequences.

Arabidopsis basic helix-loop-helix transcription factors MYC2, MYC3, and MYC4 regulate glucosinolate biosynthesis, insect performance, and feeding behavior

Arabidopsis thaliana plants fend off insect attack by constitutive and inducible production of toxic metabolites, such as glucosinolates (GSs). A triple mutant lacking MYC2, MYC3, and MYC4, three basic helix-loop-helix transcription factors that are known to additively control jasmonate-related defense responses, was shown to have a highly reduced expression of GS biosynthesis genes. The myc2 myc3 myc4 (myc234) triple mutant was almost completely devoid of GS and was extremely susceptible to the generalist herbivore Spodoptera littoralis. On the contrary, the specialist Pieris brassicae was unaffected by the presence of GS and preferred to feed on wild-type plants. In addition, lack of GS in myc234 drastically modified S. littoralis feeding behavior. Surprisingly, the expression of MYB factors known to regulate GS biosynthesis genes was not altered in myc234, suggesting that MYC2/MYC3/MYC4 are necessary for direct transcriptional activation of GS biosynthesis genes. To support this, chromatin immunoprecipitation analysis showed that MYC2 binds directly to the promoter of several GS biosynthesis genes in vivo. Furthermore, yeast two-hybrid and pull-down experiments indicated that MYC2/MYC3/MYC4 interact directly with GS-related MYBs. This specific MYC-MYB interaction plays a crucial role in the regulation of defense secondary metabolite production and underlines the importance of GS in shaping plant interactions with adapted and nonadapted herbivores.

Cucumber mosaic virus and its 2b RNA silencing suppressor modify plant-aphid interactions in tobacco

The cucumber mosaic virus (CMV) 2b protein not only inhibits anti-viral RNA silencing but also quenches transcriptional responses of plant genes to jasmonic acid, a key signalling molecule in defence against insects. This suggested that it might affect interactions between infected plants and aphids, insects that transmit CMV. We found that infection of tobacco with a 2b gene deletion mutant (CMVΔ2b) induced strong resistance to aphids (Myzus persicae) while CMV infection fostered aphid survival. Using electrical penetration graph methodology we found that higher proportions of aphids showed sustained phloem ingestion on CMV-infected plants than on CMVΔ2b-infected or mock-inoculated plants although this did not increase the rate of growth of individual aphids. This indicates that while CMV infection or certain viral gene products might elicit aphid resistance, the 2b protein normally counteracts this during a wild-type CMV infection. Our findings suggest that the 2b protein could indirectly affect aphid-mediated virus transmission.

Transgenerational epigenetic instability is a source of novel methylation variants

Epigenetic information, which may affect an organism's phenotype, can be stored and stably inherited in the form of cytosine DNA methylation. Changes in DNA methylation can produce meiotically stable epialleles that affect transcription and morphology, but the rates of spontaneous gain or loss of DNA methylation are unknown. We examined spontaneously occurring variation in DNA methylation in Arabidopsis thaliana plants propagated by single-seed descent for 30 generations. We identified 114,287 CG single methylation polymorphisms and 2485 CG differentially methylated regions (DMRs), both of which show patterns of divergence compared with the ancestral state. Thus, transgenerational epigenetic variation in DNA methylation may generate new allelic states that alter transcription, providing a mechanism for phenotypic diversity in the absence of genetic mutation.

An antiviral defense role of AGO2 in plants

BACKGROUND

Argonaute (AGO) proteins bind to small-interfering (si)RNAs and micro (mi)RNAs to target RNA silencing against viruses, transgenes and in regulation of mRNAs. Plants encode multiple AGO proteins but, in Arabidopsis, only AGO1 is known to have an antiviral role.

METHODOLOGY/PRINCIPAL FINDINGS

To uncover the roles of specific AGOs in limiting virus accumulation we inoculated turnip crinkle virus (TCV) to Arabidopsis plants that were mutant for each of the ten AGO genes. The viral symptoms on most of the plants were the same as on wild type plants although the ago2 mutants were markedly hyper-susceptible to this virus. ago2 plants were also hyper-susceptible to cucumber mosaic virus (CMV), confirming that the antiviral role of AGO2 is not specific to a single virus. For both viruses, this phenotype was associated with transient increase in virus accumulation. In wild type plants the AGO2 protein was induced by TCV and CMV infection.

CONCLUSIONS/SIGNIFICANCE

Based on these results we propose that there are multiple layers to RNA-mediated defense and counter-defense in the interactions between plants and their viruses. AGO1 represents a first layer. With some viruses, including TCV and CMV, this layer is overcome by viral suppressors of silencing that can target AGO1 and a second layer involving AGO2 limits virus accumulation. The second layer is activated when the first layer is suppressed because AGO2 is repressed by AGO1 via miR403. The activation of the second layer is therefore a direct consequence of the loss of the first layer of defense.

Disruption of two defensive signaling pathways by a viral RNA silencing suppressor

The Cucumber mosaic virus (CMV) 2b counter-defense protein disrupts plant antiviral mechanisms mediated by RNA silencing and salicylic acid (SA). We used microarrays to investigate defensive gene expression in 2b-transgenic Arabidopsis thaliana plants. Surprisingly, 2b inhibited expression of few SA-regulated genes and, in some instances, enhanced the effect of SA on certain genes. Strikingly, the 2b protein inhibited changes in the expression of 90% of genes regulated by jasmonic acid (JA). Consistent with this, infection of plants with CMV, but not the 2b gene-deletion mutant CMVDelta2b, strongly inhibited JA-inducible gene expression. JA levels were unaffected by infection with either CMV or CMVDelta2b. Although the CMV-Arabidopsis interaction is a compatible one, SA accumulation, usually considered to be an indicator of plant resistance, was increased in CMV-infected plants but not in CMVDelta2b-infected plants. Thus, the 2b protein inhibits JA signaling at a step downstream of JA biosynthesis but it primes induction of SA biosynthesis by another CMV gene product or by the process of infection itself. Like many plant viruses, CMV is aphid transmitted. JA is important in plant defense against insects. This raises the possibility that disruption of JA-mediated gene expression by the 2b protein may influence CMV transmission by aphids.

RNA silencing in plants: Flash report!

Earlier this year plant scientists met in Santa Fe, New Mexico at the Keystone Symposium "RNA Silencing Mechanisms in Plants". Sessions included small RNA biogenesis and signalling, development and stress responses, small RNA-directed DNA methylation, and interaction with pathogens. This report highlights some of the prominent and recurring themes at the meeting and emerging arenas of future research.

Symptom induction and RNA silencing suppression by the cucumber mosaic virus 2b protein

The 2b protein encoded by Cucumber mosaic virus (CMV) and other cucumoviruses is multifunctional, having roles in local and systemic virus movement, symptom determination, evasion of defense mediated by salicylic acid, and in suppression of antiviral RNA silencing. It also perturbs silencing-mediated regulation of host transcripts, suggesting that another function of 2b protein is to manipulate host gene expression and physiology in a way that may aid the virus. The 2b proteins encoded by the various cucumoviruses (CMV strains, as well as Tomato aspermy virus and Peanut stunt virus) share conserved amino acid sequence domains, suggesting that these might determine specific functions of the protein. We analyzed the effect of mutations in these domains on functions of the 2b protein during viral infection. This revealed that binding of short RNAs, the key determinants of RNA silencing specificity, correlates with RNA silencing suppression activity. Two putative phosphorylation sites were found to be required for virus symptom induction, despite having no influence on RNA silencing suppression. This indicates that the ability to suppress silencing is not the only factor affecting symptom induction by the 2b protein. In accordance with this, our studies also revealed that the 2b protein acts synergistically with some other CMV product(s) to induce symptoms, and that the role of the 2b protein in symptom determination is host species specific.

Cucumber mosaic virus 2b protein subcellular targets and interactions: their significance to RNA silencing suppressor activity

The RNA silencing suppressor activity of the 2b protein of Cucumber mosaic virus (CMV) has been variously attributed to its nuclear targeting, its interaction with and inhibition of Argonaute 1 (AGO1), or its ability to bind small RNAs in vitro. In addition, the 2b ortholog of Tomato aspermy virus forms aggregates and binds RNAs in vitro. We have further studied the relationships between CMV 2b protein silencing suppressor activity and its subcellular distribution, protein-protein interactions in vivo, and interactions with small interfering RNAs in vitro. To do this, we tagged the protein with fluorescent markers and showed that it retained suppressor activity. We showed that the 2b protein is present in the nucleolus and that it self-interacts and interacts with AGO1 and AGO4 in vivo. Using a battery of mutants, we showed that the putative nuclear localization signals and phosphorylation motif of the 2b protein are not required for self-interaction or for interaction with AGO proteins. The occurrence of neither of these interactions or of nucleolar targeting was sufficient to provide local silencing-suppression activity. In contrast, the ability of the 2b protein to bind small RNAs appears to be indispensable for silencing suppressor function.

Effects of DICER-like proteins 2, 3 and 4 on cucumber mosaic virus and tobacco mosaic virus infections in salicylic acid-treated plants

Salicylic acid (SA)-mediated resistance and RNA silencing are both important plant antiviral defence mechanisms. To investigate overlap between these resistance phenomena, we examined the ability of mutant Arabidopsis thaliana plants lacking DICER-like (DCL) endoribonucleases 2, 3 and 4 to exhibit SA-induced defence. We found that in dcl2/3/4 triple mutant plants, treatment with exogenous SA stimulated resistance to two positive-sense RNA viruses: cucumber mosaic virus and tobacco mosaic virus. We conclude that DCLs 2, 3 and 4, which are the predominant DCL endoribonucleases involved in silencing of positive-sense RNA viruses, are not required for effective SA-induced resistance to these viruses. However, the findings do not exclude RNA silencing from making a contribution to SA-mediated resistance in wild-type plants.