Genome-scale, single-cell-type resolution of microRNA activities within a whole plant organ.

Loaded into ARGONAUTE(AGO) proteins, eukaryotic micro(mi)RNAs regulate gene expression via cleavage, translational repression, and/or accelerated decay of sequence-complementary target transcripts. Despite their importance in development, cell identity maintenance and stress responses, how individual miRNAs contribute to spatial gene regulation within the complex cell mosaics formed in tissues/organs has remained inaccessible in any organism to date. We have developed a non-invasive methodology to examine, at single-cell-type resolution, the AGO-loading and activity patterns of entire miRNA cohorts in intact organs, applied here to the Arabidopsis root tip. A dual miRNAome-targetome analytical interface allowing intuitive data integration/visualization was developed as the basis for in-depth investigations via single-cell-type experimentation. These uncovered an array of so far speculative or hitherto unknown types of spatial miRNA-mediated gene regulation schemes, including via widespread cell-to-cell movement between contiguous layers of distinct identities. This study provides the proof of principle that minimally invasive, genome-scale analysis of miRNA activities within and between single-cell types of whole organs is achievable.

A new full-length circular DNA sequencing method for viral-sized genomes reveals that RNAi transgenic plants provoke a shift in geminivirus populations in the field.

We present a new method, CIDER-Seq (Circular DNA Enrichment sequencing) for the unbiased enrichment and long-read sequencing of viral-sized circular DNA molecules. We used CIDER-Seq to produce single-read full-length virus genomes for the first time. CIDER-Seq combines PCR-free virus enrichment with Single Molecule Real Time sequencing and a new sequence de-concatenation algorithm. We apply our technique to produce >1200 full-length, highly accurate geminivirus genomes from RNAi-transgenic and control plants in a field trial in Kenya. Using CIDER-Seq we can demonstrate for the first time that the expression of antiviral double-stranded RNA (dsRNA) in transgenic plants causes a consistent shift in virus populations towards species sharing low homology to the transgene derived dsRNA. Our method and its application in an economically important crop plant opens new possibilities in periodic virus sequence surveillance and accurate profiling of diverse circular DNA elements.

Diurnal changes in concerted plant protein phosphorylation and acetylation in Arabidopsis organs and seedlings.

Protein phosphorylation and acetylation are the two most abundant post-translational modifications (PTMs) that regulate protein functions in eukaryotes. In plants, these PTMs have been investigated individually; however, their co-occurrence and dynamics on proteins is currently unknown. Using Arabidopsis thaliana, we quantified changes in protein phosphorylation, acetylation and protein abundance in leaf rosettes, roots, flowers, siliques and seedlings at the end of day (ED) and at the end of night (EN). This identified 2549 phosphorylated and 909 acetylated proteins, of which 1724 phosphorylated and 536 acetylated proteins were also quantified for changes in PTM abundance between ED and EN. Using a sequential dual-PTM workflow, we identified significant PTM changes and intersections in these organs and plant developmental stages. In particular, cellular process-, pathway- and protein-level analyses reveal that the phosphoproteome and acetylome predominantly intersect at the pathway- and cellular process-level at ED versus EN. We found 134 proteins involved in core plant cell processes, such as light harvesting and photosynthesis, translation, metabolism and cellular transport, that were both phosphorylated and acetylated. Our results establish connections between PTM motifs, PTM catalyzing enzymes and putative substrate networks. We also identified PTM motifs for further characterization of the regulatory mechanisms that control cellular processes during the diurnal cycle in different Arabidopsis organs and seedlings. The sequential dual-PTM analysis expands our understanding of diurnal plant cell regulation by PTMs and provides a useful resource for future analyses, while emphasizing the importance of analyzing multiple PTMs simultaneously to elucidate when, where and how they are involved in plant cell regulation.

Haplotype-resolved genomes of geminivirus-resistant and geminivirus-susceptible African cassava cultivars.

BACKGROUND: Cassava is an important food crop in tropical and sub-tropical regions worldwide. In Africa, cassava production is widely affected by cassava mosaic disease (CMD), which is caused by the African cassava mosaic geminivirus that is transmitted by whiteflies. Cassava breeders often use a single locus, CMD2, for introducing CMD resistance into susceptible cultivars. The CMD2 locus has been genetically mapped to a 10-Mbp region, but its organization and genes as well as their functions are unknown. RESULTS: We report haplotype-resolved de novo assemblies and annotations of the genomes for the African cassava cultivar TME (tropical Manihot esculenta), which is the origin of CMD2, and the CMD-susceptible cultivar 60444. The assemblies provide phased haplotype information for over 80% of the genomes. Haplotype comparison identified novel features previously hidden in collapsed and fragmented cassava genomes, including thousands of allelic variants, inter-haplotype diversity in coding regions, and patterns of diversification through allele-specific expression. Reconstruction of the CMD2 locus revealed a highly complex region with nearly identical gene sets but limited microsynteny between the two cultivars. CONCLUSIONS: The genome maps of the CMD2 locus in both 60444 and TME3, together with the newly annotated genes, will help the identification of the causal genetic basis of CMD2 resistance to geminiviruses. Our de novo cassava genome assemblies will also facilitate genetic mapping approaches to narrow the large CMD2 region to a few candidate genes for better informed strategies to develop robust geminivirus resistance in susceptible cassava cultivars.

Genome-scale analysis of regulatory protein acetylation enzymes from photosynthetic eukaryotes.

BACKGROUND: Reversible protein acetylation occurring on Lys-Ne has emerged as a key regulatory post-translational modification in eukaryotes. It is mediated by two groups of enzymes: lysine acetyltransferases (KATs) and lysine deacetylases (KDACs) that catalyze the addition and removal of acetyl groups from target proteins. Estimates indicate that protein acetylation is second to protein phosphorylation in abundance, with thousands of acetylated sites now identified in different subcellular compartments. Considering the important regulatory role of protein phosphorylation, elucidating the diversity of KATs and KDACs across photosynthetic eukaryotes is essential in furthering our understanding of the impact of reversible protein acetylation on plant cell processes. RESULTS: We report a genome-scale analysis of lysine acetyltransferase (KAT)- and lysine deacetylase (KDAC)-families from 53 photosynthetic eukaryotes. KAT and KDAC orthologs were identified in sequenced genomes ranging from glaucophytes and algae to land plants and then analyzed for evolutionary relationships. Based on consensus molecular phylogenetic and subcellular localization data we found new sub-classes of enzymes in established KAT- and KDAC-families. Specifically, we identified a non-photosynthetic origin of the HD-tuin family KDACs, a new monocot-specific Class I HDA-family sub-class, and a phylogenetically distinct Class II algal/heterokont sub-class which maintains an ankyrin domain not conserved in land plant Class II KDACs. Protein structure analysis showed that HDA- and SRT-KDACs exist as bare catalytic subunits with highly conserved median protein length, while all KATs maintained auxiliary domains, with CBP- and TAFII250-KATs displaying protein domain gain and loss over the course of photosynthetic eukaryote evolution in addition to variable protein length. Lastly, promoter element enrichment analyses across species revealed conserved cis-regulatory sequences that support KAT and KDAC involvement in the regulation of plant development, cold/drought stress response, as well as cellular processes such as the circadian clock. CONCLUSIONS: Our results reveal new evolutionary, structural, and biological insights into the KAT- and KDAC-families of photosynthetic eukaryotes, including evolutionary parallels to protein kinases and protein phosphatases. Further, we provide a comprehensive annotation framework through our extensive phylogenetic analysis, from which future research investigating aspects of protein acetylation in plants can use to position new findings in a broader context.

Large-Scale Proteomics of the Cassava Storage Root and Identification of a Target Gene to Reduce Postharvest Deterioration.

Cassava (Manihot esculenta) is the most important root crop in the tropics, but rapid postharvest physiological deterioration (PPD) of the root is a major constraint to commercial cassava production. We established a reliable method for image-based PPD symptom quantification and used label-free quantitative proteomics to generate an extensive cassava root and PPD proteome. Over 2600 unique proteins were identified in the cassava root, and nearly 300 proteins showed significant abundance regulation during PPD. We identified protein abundance modulation in pathways associated with oxidative stress, phenylpropanoid biosynthesis (including scopoletin), the glutathione cycle, fatty acid α-oxidation, folate transformation, and the sulfate reduction II pathway. Increasing protein abundances and enzymatic activities of glutathione-associated enzymes, including glutathione reductases, glutaredoxins, and glutathione S-transferases, indicated a key role for ascorbate/glutathione cycles. Based on combined proteomics data, enzymatic activities, and lipid peroxidation assays, we identified glutathione peroxidase as a candidate for reducing PPD. Transgenic cassava overexpressing a cytosolic glutathione peroxidase in storage roots showed delayed PPD and reduced lipid peroxidation as well as decreased H2O2 accumulation. Quantitative proteomics data from ethene and phenylpropanoid pathways indicate additional gene candidates to further delay PPD. Cassava root proteomics data are available at www.pep2pro.ethz.ch for easy access and comparison with other proteomics data.

Diurnal changes in the histone H3 signature H3K9ac|H3K27ac|H3S28p are associated with diurnal gene expression in Arabidopsis.

Post-translational chromatin modifications are an important regulatory mechanism in light signalling and circadian clock function. The regulation of diurnal transcript level changes requires fine-tuning of the expression of generally active genes depending on the prevailing environmental conditions. We investigated the association of histone modifications H3K4me3, H3K9ac, H3K9me2, H3S10p, H3K27ac, H3K27me3 and H3S28p with diurnal changes in transcript expression using chromatin immunoprecipitations followed by sequencing (ChIP-Seq) in fully expanded leaves 6 of Arabidopsis thaliana grown in short-day optimal and water-deficit conditions. We identified a differential H3K9ac, H3K27ac and H3S28p signature between end-of-day and end-of-night that is correlated with changes in diurnal transcript levels. Genes with this signature have particular over-represented promoter elements and encode proteins that are significantly enriched for transcription factors, circadian clock and starch catabolic process. Additional activating modifications were prevalent in optimally watered (H3S10p) and in water-deficit (H3K4me3) plants. The data suggest a mechanism for diurnal transcript level regulation in which reduced binding of repressive transcription factors facilitates activating H3K9ac, H3K27ac and H3S28p chromatin modifications. The presence of activating chromatin modification patterns on genes only at times of the day when their expression is required can explain why some genes are differentially inducible during the diurnal cycle.

Protein abundance changes and ubiquitylation targets identified after inhibition of the proteasome with syringolin A.

As proteins are the main effectors inside cells, their levels need to be tightly regulated. This is partly achieved by specific protein degradation via the Ubiquitin-26S proteasome system (UPS). In plants, an exceptionally high number of proteins are involved in Ubiquitin-26S proteasome system-mediated protein degradation and it is known to regulate most, if not all, important cellular processes. Here, we investigated the response to the inhibition of the proteasome at the protein level treating leaves with the specific inhibitor Syringolin A (SylA) in a daytime specific manner and found 109 accumulated and 140 decreased proteins. The patterns of protein level changes indicate that the accumulating proteins cause proteotoxic stress that triggers various responses. Comparing protein level changes in SylA treated with those in a transgenic line over-expressing a mutated ubiquitin unable to form polyubiquitylated proteins produced little overlap pointing to different response pathways. To distinguish between direct and indirect targets of the UPS we also enriched and identified ubiquitylated proteins after inhibition of the proteasome, revealing a total of 1791 ubiquitylated proteins in leaves and roots and 1209 that were uniquely identified in our study. The comparison of the ubiquitylated proteins with those changing in abundance after SylA-mediated inhibition of the proteasome confirmed the complexity of the response and revealed that some proteins are regulated both at transcriptional and post-transcriptional level. For the ubiquitylated proteins that accumulate in the cytoplasm but are targeted to the plastid or the mitochondrion, we often found peptides in their target sequences, demonstrating that the UPS is involved in controlling organellar protein levels. Attempts to identify the sites of ubiquitylation revealed that the specific properties of this post-translational modification can lead to incorrect peptide spectrum assignments in complex peptide mixtures in which only a small fraction of peptides is expected to carry the ubiquitin footprint. This was confirmed with measurements of synthetically produced peptides and calculating the similarities between the different spectra.

Plastid proteome assembly without Toc159: photosynthetic protein import and accumulation of N-acetylated plastid precursor proteins.

Import of nuclear-encoded precursor proteins from the cytosol is an essential step in chloroplast biogenesis that is mediated by protein translocon complexes at the inner and outer envelope membrane (TOC). Toc159 is thought to be the main receptor for photosynthetic proteins, but lacking a large-scale systems approach, this hypothesis has only been tested for a handful of photosynthetic and nonphotosynthetic proteins. To assess Toc159 precursor specificity, we quantitatively analyzed the accumulation of plastid proteins in two mutant lines deficient in this receptor. Parallel genome-wide transcript profiling allowed us to discern the consequences of impaired protein import from systemic transcriptional responses that contribute to the loss of photosynthetic capacity. On this basis, we defined putative Toc159-independent and Toc159-dependent precursor proteins. Many photosynthetic proteins accumulate in Toc159-deficient plastids, and, surprisingly, several distinct metabolic pathways are negatively affected by Toc159 depletion. Lack of Toc159 furthermore affects several proteins that accumulate as unprocessed N-acetylated precursor proteins outside of plastids. Together, our data show an unexpected client protein promiscuity of Toc159 that requires a far more differentiated view of Toc159 receptor function and regulation of plastid protein import, in which cytosolic Met removal followed by N-terminal acetylation of precursors emerges as an additional regulatory step.

Psilocybin enhances insightfulness in meditation: a perspective on the global topology of brain imaging during meditation.

In this study, for the first time, we explored a dataset of functional magnetic resonance images collected during focused attention and open monitoring meditation before and after a five-day psilocybin-assisted meditation retreat using a recently established approach, based on the Mapper algorithm from topological data analysis. After generating subject-specific maps for two groups (psilocybin vs. placebo, 18 subjects/group) of experienced meditators, organizational principles were uncovered using graph topological tools, including the optimal transport (OT) distance, a geometrically rich measure of similarity between brain activity patterns. This revealed characteristics of the topology (i.e. shape) in space (i.e. abstract space of voxels) and time dimension of whole-brain activity patterns during different styles of meditation and psilocybin-induced alterations. Most interestingly, we found that (psilocybin-induced) positive derealization, which fosters insightfulness specifically when accompanied by enhanced open-monitoring meditation, was linked to the OT distance between open-monitoring and resting state. Our findings suggest that enhanced meta-awareness through meditation practice in experienced meditators combined with potential psilocybin-induced positive alterations in perception mediate insightfulness. Together, these findings provide a novel perspective on meditation and psychedelics that may reveal potential novel brain markers for positive synergistic effects between mindfulness practices and psilocybin.

Systems-based analysis of Arabidopsis leaf growth reveals adaptation to water deficit.

Leaves have a central role in plant energy capture and carbon conversion and therefore must continuously adapt their development to prevailing environmental conditions. To reveal the dynamic systems behaviour of leaf development, we profiled Arabidopsis leaf number six in depth at four different growth stages, at both the end-of-day and end-of-night, in plants growing in two controlled experimental conditions: short-day conditions with optimal soil water content and constant reduced soil water conditions. We found that the lower soil water potential led to reduced, but prolonged, growth and an adaptation at the molecular level without a drought stress response. Clustering of the protein and transcript data using a decision tree revealed different patterns in abundance changes across the growth stages and between end-of-day and end-of-night that are linked to specific biological functions. Correlations between protein and transcript levels depend on the time-of-day and also on protein localisation and function. Surprisingly, only very few of >1700 quantified proteins showed diurnal abundance fluctuations, despite strong fluctuations at the transcript level.

MASCP Gator: an aggregation portal for the visualization of Arabidopsis proteomics data.

Proteomics has become a critical tool in the functional understanding of plant processes at the molecular level. Proteomics-based studies have also contributed to the ever-expanding array of data in modern biology, with many generating Web portals and online resources that contain incrementally expanding and updated information. Many of these resources reflect specialist research areas with significant and novel information that is not currently captured by centralized repositories. The Arabidopsis (Arabidopsis thaliana) community is well served by a number of online proteomics resources that hold an abundance of functional information. These sites can be difficult to locate among a multitude of online resources. Furthermore, they can be difficult to navigate in order to identify specific features of interest without significant technical knowledge. Recently, members of the Arabidopsis proteomics community involved in developing many of these resources decided to develop a summary aggregation portal that is capable of retrieving proteomics data from a series of online resources on the fly. The Web portal is known as the MASCP Gator and can be accessed at the following address: http://gator.masc-proteomics.org/. Significantly, proteomics data displayed at this site retrieve information from the data repositories upon each request. This means that information is always up to date and displays the latest data sets. The site also provides hyperlinks back to the source information hosted at each of the curated databases to facilitate more in-depth analysis of the primary data.

Full-length sequencing of circular DNA viruses and extrachromosomal circular DNA using CIDER-Seq.

Circular DNA is ubiquitous in nature in the form of plasmids, circular DNA viruses, and extrachromosomal circular DNA (eccDNA) in eukaryotes. Sequencing of such molecules is essential to profiling virus distributions, discovering new viruses and understanding the roles of eccDNAs in eukaryotic cells. Circular DNA enrichment sequencing (CIDER-Seq) is a technique to enrich and accurately sequence circular DNA without the need for polymerase chain reaction amplification, cloning, and computational sequence assembly. The approach is based on randomly primed circular DNA amplification, which is followed by several enzymatic DNA repair steps and then by long-read sequencing. CIDER-Seq includes a custom data analysis package (CIDER-Seq Data Analysis Software 2) that implements the DeConcat algorithm to deconcatenate the long sequencing products of random circular DNA amplification into the intact sequences of the input circular DNA. The CIDER-Seq data analysis package can generate full-length annotated virus genomes, as well as circular DNA sequences of novel viruses. Applications of CIDER-Seq also include profiling of eccDNA molecules such as transposable elements (TEs) from biological samples. The method takes ~2 weeks to complete, depending on the computational resources available. Owing to the present constraints of long-read single-molecule sequencing, the accuracy of circular virus and eccDNA sequences generated by the CIDER-Seq method scales with sequence length, and the greatest accuracy is obtained for molecules <10 kb long.

Linking CRISPR-Cas9 interference in cassava to the evolution of editing-resistant geminiviruses.

BACKGROUND: Geminiviruses cause damaging diseases in several important crop species. However, limited progress has been made in developing crop varieties resistant to these highly diverse DNA viruses. Recently, the bacterial CRISPR/Cas9 system has been transferred to plants to target and confer immunity to geminiviruses. In this study, we use CRISPR-Cas9 interference in the staple food crop cassava with the aim of engineering resistance to African cassava mosaic virus, a member of a widespread and important family (Geminiviridae) of plant-pathogenic DNA viruses. RESULTS: Our results show that the CRISPR system fails to confer effective resistance to the virus during glasshouse inoculations. Further, we find that between 33 and 48% of edited virus genomes evolve a conserved single-nucleotide mutation that confers resistance to CRISPR-Cas9 cleavage. We also find that in the model plant Nicotiana benthamiana the replication of the novel, mutant virus is dependent on the presence of the wild-type virus. CONCLUSIONS: Our study highlights the risks associated with CRISPR-Cas9 virus immunity in eukaryotes given that the mutagenic nature of the system generates viral escapes in a short time period. Our in-depth analysis of virus populations also represents a template for future studies analyzing virus escape from anti-viral CRISPR transgenics. This is especially important for informing regulation of such actively mutagenic applications of CRISPR-Cas9 technology in agriculture.

Design and Use of a Digitally Controlled Device for Accurate, Multiplexed Gas Exchange Measurements of the Complete Foliar Parts of Plants.

Performing accurate measurements of photosynthetic and respiration rates is vital to a large proportion of plant-based studies. While several commercial systems exist to perform such measurements, few are ideal for whole-plant measurements of small herbaceous plants such as Arabidopsis and none offer the capacity for simultaneous analysis of multiple plants. We, therefore, designed a multi-chamber, computer-controlled, infrared gas analyzer-coupled system for the continuous measurement of gas exchange in whole-plant shoots or rosettes. This system was called ETH Gas Exchange System-1 (EGES-1). We have subsequently expanded the device to accommodate a wider variety of species while providing precise control over environmental parameters. Critically, we have (1) increased the flow rates through each of the eight chambers, (2) introduced a computer-controlled feedback loop for the precise introduction of CO2, and (3) added an additional feedback loop for the introduction and control of humidity. The advantages of this new system (EGES-2) are illustrated here in the context of a variety of physiological experiments.

Web-based analysis of the mouse transcriptome using Genevestigator.

BACKGROUND: Gene function analysis often requires a complex and laborious sequence of laboratory and computer-based experiments. Choosing an effective experimental design generally results from hypotheses derived from prior knowledge or experimentation. Knowledge obtained from meta-analyzing compendia of expression data with annotation libraries can provide significant clues in understanding gene and network function, resulting in better hypotheses that can be tested in the laboratory. DESCRIPTION: Genevestigator is a microarray database and analysis system allowing context-driven queries. Simple but powerful tools allow biologists with little computational background to retrieve information about when, where and how genes are expressed. We manually curated and quality-controlled 3110 mouse Affymetrix arrays from public repositories. Data queries can be run against an annotation library comprising 160 anatomy categories, 12 developmental stage groups, 80 stimuli, and 182 genetic backgrounds or modifications. The quality of results obtained through Genevestigator is illustrated by a number of biological scenarios that are substantiated by other types of experimentation in the literature. CONCLUSION: The Genevestigator-Mouse database effectively provides biologically meaningful results and can be accessed at https://www.genevestigator.ethz.ch.

pep2pro: the high-throughput proteomics data processing, analysis, and visualization tool.

The pep2pro database was built to support effective high-throughput proteome data analysis. Its database schema allows the coherent integration of search results from different database-dependent search algorithms and filtering of the data including control for unambiguous assignment of peptides to proteins. The capacity of the pep2pro database has been exploited in data analysis of various Arabidopsis proteome datasets. The diversity of the datasets and the associated scientific questions required thorough querying of the data. This was supported by the relational format structure of the data that links all information on the sample, spectrum, search database, and algorithm to peptide and protein identifications and their post-translational modifications. After publication of datasets they are made available on the pep2pro website at www.pep2pro.ethz.ch. Further, the pep2pro data analysis pipeline also handles data export do the PRIDE database (http://www.ebi.ac.uk/pride) and data retrieval by the MASCP Gator (http://gator.masc-proteomics.org/). The utility of pep2pro will continue to be used for analysis of additional datasets and as a data warehouse. The capacity of the pep2pro database for proteome data analysis has now also been made publicly available through the release of pep2pro4all, which consists of a database schema and a script that will populate the database with mass spectrometry data provided in mzIdentML format.