Predmoter-cross-species prediction of plant promoter and enhancer regions.

Motivation: Identifying cis-regulatory elements (CREs) is crucial for analyzing gene regulatory networks. Next generation sequencing methods were developed to identify CREs but represent a considerable expenditure for targeted analysis of few genomic loci. Thus, predicting the outputs of these methods would significantly cut costs and time investment. Results: We present Predmoter, a deep neural network that predicts base-wise Assay for Transposase Accessible Chromatin using sequencing (ATAC-seq) and histone Chromatin immunoprecipitation DNA-sequencing (ChIP-seq) read coverage for plant genomes. Predmoter uses only the DNA sequence as input. We trained our final model on 21 species for 13 of which ATAC-seq data and for 17 of which ChIP-seq data was publicly available. We evaluated our models on Arabidopsis thaliana and Oryza sativa. Our best models showed accurate predictions in peak position and pattern for ATAC- and histone ChIP-seq. Annotating putatively accessible chromatin regions provides valuable input for the identification of CREs. In conjunction with other in silico data, this can significantly reduce the search space for experimentally verifiable DNA-protein interaction pairs. Availability and implementation: The source code for Predmoter is available at: https://github.com/weberlab-hhu/Predmoter. Predmoter takes a fasta file as input and outputs h5, and optionally bigWig and bedGraph files.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity.

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene ß-cyclase whose function we demonstrate.

Helixer: cross-species gene annotation of large eukaryotic genomes using deep learning.

MOTIVATION: Current state-of-the-art tools for the de novo annotation of genes in eukaryotic genomes have to be specifically fitted for each species and still often produce annotations that can be improved much further. The fundamental algorithmic architecture for these tools has remained largely unchanged for about two decades, limiting learning capabilities. Here, we set out to improve the cross-species annotation of genes from DNA sequence alone with the help of deep learning. The goal is to eliminate the dependency on a closely related gene model while also improving the predictive quality in general with a fundamentally new architecture. RESULTS: We present Helixer, a framework for the development and usage of a cross-species deep learning model that improves significantly on performance and generalizability when compared to more traditional methods. We evaluate our approach by building a single vertebrate model for the base-wise annotation of 186 animal genomes and a separate land plant model for 51 plant genomes. Our predictions are shown to be much less sensitive to the length of the genome than those of a current state-of-the-art tool. We also present two novel post-processing techniques that each worked to further strengthen our annotations and show in-depth results of an RNA-Seq based comparison of our predictions. Our method does not yet produce comprehensive gene models but rather outputs base pair wise probabilities. AVAILABILITY AND IMPLEMENTATION: The source code of this work is available at https://github.com/weberlab-hhu/Helixer under the GNU General Public License v3.0. The trained models are available at https://doi.org/10.5281/zenodo.3974409. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phosphate starvation causes different stress responses in the lipid metabolism of tomato leaves and roots.

Plants have evolved various acclimation responses to cope with phosphate depletion, including several changes in lipid metabolism. Thereby membrane phospholipids are dephosphorylated and can be used as an internal phosphate source, while galactolipids are incorporated into the membrane to maintain membrane functionality. Still little is known about the lipidomic and transcriptomic response of plants other than Arabidopsis thaliana upon phosphate starvation. Therefore, we employed lipidomics and transcriptomics to characterize the phosphate starvation response of lipid metabolism in tomato leaves and roots. Overall, phospholipid levels decreased and galactolipids increased during the acclimation response. In addition, an early increase of triacylglycerol was observed. Interestingly, there were major differences in the acclimation response of tomato leaves and roots: leaves mainly accumulated polyunsaturated triacylglycerol, while roots showed a massive increase in galactolipid content. In line with these results, we observed transcriptional induction of phospholipid degradation and galactolipid synthesis pathways in both analyzed tissues. In contrast, other aspects of the transcriptional response, in particular, the induction of phospholipid degradation, ER-localized fatty acid desaturation and triacylglycerol assembly differed between tomato leaves and roots. These results suggest a different modulation of degraded phospholipids toward triacylglycerols and galactolipids in phosphate-starved tomato leaves and roots. Possibly the availability and composition of acyl-CoA pools and ER-derived precursors trigger the synthesis of triacylglycerols or galactolipids. As the mechanism of triacylglycerol accumulation is poorly characterized outside of seed oil formation, these findings enhance our understanding of the phosphate starvation response and of how storage lipids accumulate under stress in vegetative tissue.

Modeling Protein Destiny in Developing Fruit.

Protein synthesis and degradation are essential processes that regulate cell status. Because labeling in bulky organs, such as fruits, is difficult, we developed a modeling approach to study protein turnover at the global scale in developing tomato (Solanum lycopersicum) fruit. Quantitative data were collected for transcripts and proteins during fruit development. Clustering analysis showed smaller changes in protein abundance compared to mRNA abundance. Furthermore, protein and transcript abundance were poorly correlated, and the coefficient of correlation decreased during fruit development and ripening, with transcript levels decreasing more than protein levels. A mathematical model with one ordinary differential equation was used to estimate translation (kt ) and degradation (kd ) rate constants for almost 2,400 detected transcript-protein pairs and was satisfactorily fitted for >1,000 pairs. The model predicted median values of ∼2 min for the translation of a protein, and a protein lifetime of ∼11 d. The constants were validated and inspected for biological relevance. Proteins involved in protein synthesis had higher kt and kd values, indicating that the protein machinery is particularly flexible. Our model also predicts that protein concentration is more strongly affected by the rate of translation than that of degradation.

Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris.

A parasitic lifestyle, where plants procure some or all of their nutrients from other living plants, has evolved independently in many dicotyledonous plant families and is a major threat for agriculture globally. Nevertheless, no genome sequence of a parasitic plant has been reported to date. Here we describe the genome sequence of the parasitic field dodder, Cuscuta campestris. The genome contains signatures of a fairly recent whole-genome duplication and lacks genes for pathways superfluous to a parasitic lifestyle. Specifically, genes needed for high photosynthetic activity are lost, explaining the low photosynthesis rates displayed by the parasite. Moreover, several genes involved in nutrient uptake processes from the soil are lost. On the other hand, evidence for horizontal gene transfer by way of genomic DNA integration from the parasite's hosts is found. We conclude that the parasitic lifestyle has left characteristic footprints in the C. campestris genome.

De Novo Assembly of a New Solanum pennellii Accession Using Nanopore Sequencing.

Updates in nanopore technology have made it possible to obtain gigabases of sequence data. Prior to this, nanopore sequencing technology was mainly used to analyze microbial samples. Here, we describe the generation of a comprehensive nanopore sequencing data set with a median read length of 11,979 bp for a self-compatible accession of the wild tomato species Solanum pennellii We describe the assembly of its genome to a contig N50 of 2.5 MB. The assembly pipeline comprised initial read correction with Canu and assembly with SMARTdenovo. The resulting raw nanopore-based de novo genome is structurally highly similar to that of the reference S. pennellii LA716 accession but has a high error rate and was rich in homopolymer deletions. After polishing the assembly with Illumina reads, we obtained an error rate of <0.02% when assessed versus the same Illumina data. We obtained a gene completeness of 96.53%, slightly surpassing that of the reference S. pennellii Taken together, our data indicate that such long read sequencing data can be used to affordably sequence and assemble gigabase-sized plant genomes.

Freeze-quenched maize mesophyll and bundle sheath separation uncovers bias in previous tissue-specific RNA-Seq data.

The high efficiency of C4 photosynthesis relies on spatial division of labor, classically with initial carbon fixation in the mesophyll and carbon reduction in the bundle sheath. By employing grinding and serial filtration over liquid nitrogen, we enriched C4 tissues along a developing leaf gradient. This method treats both C4 tissues in an integrity-preserving and consistent manner, while allowing complementary measurements of metabolite abundance and enzyme activity, thus providing a comprehensive data set. Meta-analysis of this and the previous studies highlights the strengths and weaknesses of different C4 tissue separation techniques. While the method reported here achieves the least enrichment, it is the only one that shows neither strong 3' (degradation) bias, nor different severity of 3' bias between samples. The meta-analysis highlighted previously unappreciated observations, such as an accumulation of evidence that aspartate aminotransferase is more mesophyll specific than expected from the current NADP-ME C4 cycle model, and a shift in enrichment of protein synthesis genes from bundle sheath to mesophyll during development. The full comparative dataset is available for download, and a web visualization tool (available at http://www.plant-biochemistry.hhu.de/resources.html) facilitates comparison of the the Z. mays bundle sheath and mesophyll studies, their consistencies and their conflicts.

Understanding metabolite transport and metabolism in C4 plants through RNA-seq.

RNA-seq, the measurement of steady-state RNA levels by next generation sequencing, has enabled quantitative transcriptome analyses of complex traits in many species without requiring the parallel sequencing of their genomes. The complex trait of C4 photosynthesis, which increases photosynthetic efficiency via a biochemical pump that concentrates CO2 around RubisCO, has evolved convergently multiple times. Due to these interesting properties, C4 photosynthesis has been analyzed in a series of comparative RNA-seq projects. These projects compared both species with and without the C4 trait and different tissues or organs within a C4 plant. The RNA-seq studies were evaluated by comparing to earlier single gene studies. The studies confirmed the marked changes expected for C4 signature genes, but also revealed numerous new players in C4 metabolism showing that the C4 cycle is more complex than previously thought, and suggesting modes of integration into the underlying C3 metabolism.

Retained duplicate genes in green alga Chlamydomonas reinhardtii tend to be stress responsive and experience frequent response gains.

BACKGROUND: Green algae belong to a group of photosynthetic organisms that occupy diverse habitats, are closely related to land plants, and have been studied as sources of food and biofuel. Although multiple green algal genomes are available, a global comparative study of algal gene families has not been carried out. To investigate how gene families and gene expression have evolved, particularly in the context of stress response that have been shown to correlate with gene family expansion in multiple eukaryotes, we characterized the expansion patterns of gene families in nine green algal species, and examined evolution of stress response among gene duplicates in Chlamydomonas reinhardtii. RESULTS: Substantial variation in domain family sizes exists among green algal species. Lineage-specific expansion of families occurred throughout the green algal lineage but inferred gene losses occurred more often than gene gains, suggesting a continuous reduction of algal gene repertoire. Retained duplicates tend to be involved in stress response, similar to land plant species. However, stress responsive genes tend to be pseudogenized as well. When comparing ancestral and extant gene stress response state, we found that response gains occur in 13% of duplicate gene branches, much higher than 6% in Arabidopsis thaliana. CONCLUSION: The frequent gains of stress response among green algal duplicates potentially reflect a high rate of innovation, resulting in a species-specific gene repertoire that contributed to adaptive response to stress. This could be further explored towards deciphering the mechanism of stress response, and identifying suitable green algal species for oil production.

Comparative transcriptome atlases reveal altered gene expression modules between two Cleomaceae C3 and C4 plant species.

C(4) photosynthesis outperforms the ancestral C(3) state in a wide range of natural and agro-ecosystems by affording higher water-use and nitrogen-use efficiencies. It therefore represents a prime target for engineering novel, high-yielding crops by introducing the trait into C(3) backgrounds. However, the genetic architecture of C(4) photosynthesis remains largely unknown. To define the divergence in gene expression modules between C(3) and C(4) photosynthesis during leaf ontogeny, we generated comprehensive transcriptome atlases of two Cleomaceae species, Gynandropsis gynandra (C(4)) and Tarenaya hassleriana (C(3)), by RNA sequencing. Overall, the gene expression profiles appear remarkably similar between the C(3) and C(4) species. We found that known C(4) genes were recruited to photosynthesis from different expression domains in C(3), including typical housekeeping gene expression patterns in various tissues as well as individual heterotrophic tissues. Furthermore, we identified a structure-related module recruited from the C(3) root. Comparison of gene expression patterns with anatomy during leaf ontogeny provided insight into genetic features of Kranz anatomy. Altered expression of developmental factors and cell cycle genes is associated with a higher degree of endoreduplication in enlarged C(4) bundle sheath cells. A delay in mesophyll differentiation apparent both in the leaf anatomy and the transcriptome allows for extended vein formation in the C(4) leaf.

C4 photosynthesis: from evolutionary analyses to strategies for synthetic reconstruction of the trait.

C4 photosynthesis represents the most productive modes of photosynthesis in land plants and some of the most productive crops on the planet, such as maize and sugarcane, and many ecologically important native plants use this type of photosynthesis. Despite its ecological and economic importance, the genetic basis of C4 photosynthesis remains largely unknown. Even many fundamental aspects of C4 biochemistry, such as the molecular identity of solute transporters, and many aspects of C4 plant leaf development, such as the Kranz anatomy, are currently not understood. Here, we review recent progress in gaining a mechanistic understanding of the complex C4 trait through comparative evolutionary analyses of C3 and C4 species.

Systems analysis of a maize leaf developmental gradient redefines the current C4 model and provides candidates for regulation.

We systematically analyzed a developmental gradient of the third maize (Zea mays) leaf from the point of emergence into the light to the tip in 10 continuous leaf slices to study organ development and physiological and biochemical functions. Transcriptome analysis, oxygen sensitivity of photosynthesis, and photosynthetic rate measurements showed that the maize leaf undergoes a sink-to-source transition without an intermediate phase of C(3) photosynthesis or operation of a photorespiratory carbon pump. Metabolome and transcriptome analysis, chlorophyll and protein measurements, as well as dry weight determination, showed continuous gradients for all analyzed items. The absence of binary on-off switches and regulons pointed to a morphogradient along the leaf as the determining factor of developmental stage. Analysis of transcription factors for differential expression along the leaf gradient defined a list of putative regulators orchestrating the sink-to-source transition and establishment of C(4) photosynthesis. Finally, transcriptome and metabolome analysis, as well as enzyme activity measurements, and absolute quantification of selected metabolites revised the current model of maize C(4) photosynthesis. All data sets are included within the publication to serve as a resource for maize leaf systems biology.