Evaluation and integration of functional annotation pipelines for newly sequenced organisms: the potato genome as a test case

Temporally distinct regulatory pathways coordinate thermo-responsive storage organ formation in potato

Tuberization is an important developmental process in potatoes, but it is highly affected by environmental conditions. Temperature is a major environmental factor affecting tuberization, with high temperatures suppressing tuber development. However, the temporal aspects of thermo-responsive tuberization remain elusive. In this study, we show that FT homolog StSP6A is suppressed by temporally distinct regulatory pathways. Experiments using StSP6A-overexpressing plants show that post-transcriptional regulation plays a major role at the early stage, while transcriptional regulation is an important late-stage factor, suppressing StSP6A at high temperatures in leaves. Overexpression of StSP6A in leaves restores tuber formation but does not recover tuber yield at the late stage, possibly because of suppressed sugar transport at high temperatures. Transcriptome analyses lead to the identification of potential regulators that may be involved in thermo-responsive tuberization at different stages. Our work shows that potato has temporally distinct molecular mechanisms that finely control tuber development at high temperatures.

Heat-stress-induced sprouting and differential gene expression in growing potato tubers: Comparative transcriptomics with that induced by postharvest sprouting

Crops face increased risk from heat stress due to climate change. Potato (Solanum tuberosum L.) tubers grown in hot summers often have defects including pre-harvest sprouting ("heat sprouts"). We have used 18 potato cultivars to investigate whether heat stress (HS) conditions alone could cause heat sprouting and dormancy changes in tubers. We also examined transcriptomic responses of potato to HS and whether these responses are like those induced by postharvest sprouting. We demonstrated that HS alone caused heat sprouts and shortened postharvest dormancy period, heat-sprouted tubers became dormant after harvest, and cultivars varied substantially for producing heat spouts but there was no clear association with cultivar maturity earliness. Cultivar Innovator did not show any heat sprouts and still had long dormancy. Dormancy-associated genes (DOG1 and SLP) were downregulated in HS tubers like in postharvest sprouting tubers. We have identified 1201 differentially expressed genes, 14 enriched GO terms and 12 enriched KEGG pathways in response to HS in growing tubers of 'Russet Burbank'. Transcriptomic response of 'Russet Burbank' to HS showed significant similarities to that of postharvest non-HS sprouted tubers. Gibberellin biosynthesis pathway was enriched in heat-stressed tubers and was likely involved in heat sprouting and dormancy release. Heat sprouting and postharvest sprouting shared common candidate genes and had significant similarity in gene expression. Our study has significance for selecting potato cultivars for farming, planning storage and utilization of heat-stressed tubers, identifying sprouting-related genes, understanding heat-stress biology, and breeding heat-tolerant potato cultivars, especially for sustainable potato production under climate change.

Gene Ontology Meta Annotator for Plants (GOMAP)

Annotating gene structures and functions to genome assemblies is necessary to make assembly resources useful for biological inference. Gene Ontology (GO) term assignment is the most used functional annotation system, and new methods for GO assignment have improved the quality of GO-based function predictions. The Gene Ontology Meta Annotator for Plants (GOMAP) is an optimized, high-throughput, and reproducible pipeline for genome-scale GO annotation of plants. We containerized GOMAP to increase portability and reproducibility and also optimized its performance for HPC environments. Here we report on the pipeline's availability and performance for annotating large, repetitive plant genomes and describe how GOMAP was used to annotate multiple maize genomes as a test case. Assessment shows that GOMAP expands and improves the number of genes annotated and annotations assigned per gene as well as the quality (based on [Formula: see text]) of GO assignments in maize. GOMAP has been deployed to annotate other species including wheat, rice, barley, cotton, and soy. Instructions and access to the GOMAP Singularity container are freely available online at https://bioinformapping.com/gomap/ . A list of annotated genomes and links to data is maintained at https://dill-picl.org/projects/gomap/ .

TOA: A software package for automated functional annotation in non-model plant species

The increase of sequencing capacity provided by high-throughput platforms has made it possible to routinely obtain large sets of genomic and transcriptomic sequences from model and non-model organisms. Subsequent genomic analysis and gene discovery in next-generation sequencing experiments are, however, bottlenecked by functional annotation. One common way to perform functional annotation of sets of sequences obtained from next-generation sequencing experiments, is by searching for homologous sequences and accessing the related functional information deposited in genomic databases. Functional annotation is especially challenging for non-model organisms, like many plant species. In such cases, existing free and commercial general-purpose applications may not offer complete and accurate results. We present TOA (Taxonomy-oriented annotation), a Python-based user-friendly open source application designed to establish functional annotation pipelines geared towards non-model plant species that can run in Linux/Mac computers, HPCs and cloud servers. TOA performs homology searches against proteins stored in the PLAZA databases, NCBI RefSeq Plant, Nucleotide Database and Non-Redundant Protein Sequence Database, and outputs functional information from several ontology systems: Gene Ontology, InterPro, EC, KEGG, Mapman and MetaCyc. The software performance was validated by comparing the runtimes, total number of annotated sequences and accuracy of the functional information obtained for several plant benchmark data sets with TOA and other functional annotation solutions. TOA outperformed the other software in terms of number of annotated sequences and accuracy of the annotation and constitutes a good alternative to improve functional annotation in plants. TOA is especially recommended for gymnosperms or for low quality sequence data sets of non-model plants.

Linking crop traits to transcriptome differences in a progeny population of tetraploid potato

BACKGROUND

Potato is the third most consumed crop in the world. Breeding for traits such as yield, product quality and pathogen resistance are main priorities. Identifying molecular signatures of these and other important traits is important in future breeding efforts. In this study, a progeny population from a cross between a breeding line, SW93-1015, and a cultivar, Désirée, was studied by trait analysis and RNA-seq in order to develop understanding of segregating traits at the molecular level and identify transcripts with expressional correlation to these traits. Transcript markers with predictive value for field performance applicable under controlled environments would be of great value for plant breeding.

RESULTS

A total of 34 progeny lines from SW93-1015 and Désirée were phenotyped for 17 different traits in a field in Nordic climate conditions and controlled climate settings. A master transcriptome was constructed with all 34 progeny lines and the parents through a de novo assembly of RNA-seq reads. Gene expression data obtained in a controlled environment from the 34 lines was correlated to traits by different similarity indices, including Pearson and Spearman, as well as DUO, which calculates the co-occurrence between high and low values for gene expression and trait. Our study linked transcripts to traits such as yield, growth rate, high laying tubers, late and tuber blight, tuber greening and early flowering. We found several transcripts associated to late blight resistance and transcripts encoding receptors were associated to Dickeya solani susceptibility. Transcript levels of a UBX-domain protein was negatively associated to yield and a GLABRA2 expression modulator was negatively associated to growth rate.

CONCLUSION

In our study, we identify 100's of transcripts, putatively linked based on expression with 17 traits of potato, representing both well-known and novel associations. This approach can be used to link the transcriptome to traits. We explore the possibility of associating the level of transcript expression from controlled, optimal environments to traits in a progeny population with different methods introducing the application of DUO for the first time on transcriptome data. We verify the expression pattern for five of the putative transcript markers in another progeny population.

The effects of the neonicotinoid imidacloprid on gene expression and DNA methylation in the buff-tailed bumblebee Bombus terrestris

Neonicotinoids are effective insecticides used on many important arable and horticultural crops. They are nicotinic acetylcholine receptor agonists which disrupt the function of insect neurons and cause paralysis and death. In addition to direct mortality, there are numerous sublethal effects of low doses of neonicotinoids on bees. We hypothesize that some of these large array of effects could be a consequence of epigenetic changes in bees induced by neonicotinoids. We compared whole methylome (BS-seq) and RNA-seq libraries of the brains of buff-tailed bumblebee Bombus terrestris workers exposed to field-realistic doses of the neonicotinoid imidacloprid to libraries from control workers. We found numerous genes which show differential expression between neonicotinoid-treated bees and control bees, but no differentially methylated cytosines in any context. We found CpG methylation to be focused mainly in exons and associated with highly expressed genes. We discuss the implications of our results for future legislation.

Plant genome and transcriptome annotations: from misconceptions to simple solutions

Next-generation sequencing has triggered an explosion of available genomic and transcriptomic resources in the plant sciences. Although genome and transcriptome sequencing has become orders of magnitudes cheaper and more efficient, often the functional annotation process is lagging behind. This might be hampered by the lack of a comprehensive enumeration of simple-to-use tools available to the plant researcher. In this comprehensive review, we present (i) typical ontologies to be used in the plant sciences, (ii) useful databases and resources used for functional annotation, (iii) what to expect from an annotated plant genome, (iv) an automated annotation pipeline and (v) a recipe and reference chart outlining typical steps used to annotate plant genomes/transcriptomes using publicly available resources.

An Updated Functional Annotation of Protein-Coding Genes in the Cucumber Genome

Background: Although the cucumber reference genome and its annotation were published several years ago, the functional annotation of predicted genes, particularly protein-coding genes, still requires further improvement. In general, accurately determining orthologous relationships between genes allows for better and more robust functional assignments of predicted genes. As one of the most reliable strategies, the determination of collinearity information may facilitate reliable orthology inferences among genes from multiple related genomes. Currently, the identification of collinear segments has mainly been based on conservation of gene order and orientation. Over the course of plant genome evolution, various evolutionary events have disrupted or distorted the order of genes along chromosomes, making it difficult to use those genes as genome-wide markers for plant genome comparisons. Results: Using the localized LASTZ/MULTIZ analysis pipeline, we aligned 15 genomes, including cucumber and other related angiosperm plants, and identified a set of genomic segments that are short in length, stable in structure, uniform in distribution and highly conserved across all 15 plants. Compared with protein-coding genes, these conserved segments were more suitable for use as genomic markers for detecting collinear segments among distantly divergent plants. Guided by this set of identified collinear genomic segments, we inferred 94,486 orthologous protein-coding gene pairs (OPPs) between cucumber and 14 other angiosperm species, which were used as proxies for transferring functional terms to cucumber genes from the annotations of the other 14 genomes. In total, 10,885 protein-coding genes were assigned Gene Ontology (GO) terms which was nearly 1,300 more than results collected in Uniprot-proteomic database. Our results showed that annotation accuracy would been improved compared with other existing approaches. Conclusions: In this study, we provided an alternative resource for the functional annotation of predicted cucumber protein-coding genes, which we expect will be beneficial for the cucumber's biological study, accessible from http://cmb.bnu.edu.cn/functional_annotation. Meanwhile, using the cucumber reference genome as a case study, we presented an efficient strategy for transferring gene functional information from previously well-characterized protein-coding genes in model species to newly sequenced or "non-model" plant species.

MorphDB: Prioritizing Genes for Specialized Metabolism Pathways and Gene Ontology Categories in Plants

Recent times have seen an enormous growth of "omics" data, of which high-throughput gene expression data are arguably the most important from a functional perspective. Despite huge improvements in computational techniques for the functional classification of gene sequences, common similarity-based methods often fall short of providing full and reliable functional information. Recently, the combination of comparative genomics with approaches in functional genomics has received considerable interest for gene function analysis, leveraging both gene expression based guilt-by-association methods and annotation efforts in closely related model organisms. Besides the identification of missing genes in pathways, these methods also typically enable the discovery of biological regulators (i.e., transcription factors or signaling genes). A previously built guilt-by-association method is MORPH, which was proven to be an efficient algorithm that performs particularly well in identifying and prioritizing missing genes in plant metabolic pathways. Here, we present MorphDB, a resource where MORPH-based candidate genes for large-scale functional annotations (Gene Ontology, MapMan bins) are integrated across multiple plant species. Besides a gene centric query utility, we present a comparative network approach that enables researchers to efficiently browse MORPH predictions across functional gene sets and species, facilitating efficient gene discovery and candidate gene prioritization. MorphDB is available at http://bioinformatics.psb.ugent.be/webtools/morphdb/morphDB/index/. We also provide a toolkit, named "MORPH bulk" (https://github.com/arzwa/morph-bulk), for running MORPH in bulk mode on novel data sets, enabling researchers to apply MORPH to their own species of interest.

Identification of Preferred DNA-Binding Sites for the Thermus thermophilus Transcriptional Regulator SbtR by the Combinatorial Approach REPSA

One of the first steps towards elucidating the biological function of a putative transcriptional regulator is to ascertain its preferred DNA-binding sequences. This may be rapidly and effectively achieved through the application of a combinatorial approach, one involving very large numbers of randomized oligonucleotides and reiterative selection and amplification steps to enrich for high-affinity nucleic acid-binding sequences. Previously, we had developed the novel combinatorial approach Restriction Endonuclease Protection, Selection and Amplification (REPSA), which relies not on the physical separation of ligand-nucleic acid complexes but instead selects on the basis of ligand-dependent inhibition of enzymatic template inactivation, specifically cleavage by type IIS restriction endonucleases. Thus, no prior knowledge of the ligand is required for REPSA, making it more amenable for discovery purposes. Here we describe using REPSA, massively parallel sequencing, and bioinformatics to identify the preferred DNA-binding sites for the transcriptional regulator SbtR, encoded by the TTHA0167 gene from the model extreme thermophile Thermus thermophilus HB8. From the resulting position weight matrix, we can identify multiple operons potentially regulated by SbtR and postulate a biological role for this protein in regulating extracellular transport processes. Our study provides a proof-of-concept for the application of REPSA for the identification of preferred DNA-binding sites for orphan transcriptional regulators and a first step towards determining their possible biological roles.

The nitrogen responsive transcriptome in potato (Solanum tuberosum L.) reveals significant gene regulatory motifs

Nitrogen (N) is the most important nutrient for the growth of potato (Solanum tuberosum L.). Foliar gene expression in potato plants with and without N supplementation at 180 kg N ha(-1) was compared at mid-season. Genes with consistent differences in foliar expression due to N supplementation over three cultivars and two developmental time points were examined. In total, thirty genes were found to be over-expressed and nine genes were found to be under-expressed with supplemented N. Functional relationships between over-expressed genes were found. The main metabolic pathway represented among differentially expressed genes was amino acid metabolism. The 1000 bp upstream flanking regions of the differentially expressed genes were analysed and nine overrepresented motifs were found using three motif discovery algorithms (Seeder, Weeder and MEME). These results point to coordinated gene regulation at the transcriptional level controlling steady state potato responses to N sufficiency.

Short interspersed nuclear elements (SINEs) are abundant in Solanaceae and have a family-specific impact on gene structure and genome organization

Short interspersed nuclear elements (SINEs) are highly abundant non-autonomous retrotransposons that are widespread in plants. They are short in size, non-coding, show high sequence diversity, and are therefore mostly not or not correctly annotated in plant genome sequences. Hence, comparative studies on genomic SINE populations are rare. To explore the structural organization and impact of SINEs, we comparatively investigated the genome sequences of the Solanaceae species potato (Solanum tuberosum), tomato (Solanum lycopersicum), wild tomato (Solanum pennellii), and two pepper cultivars (Capsicum annuum). Based on 8.5 Gbp sequence data, we annotated 82 983 SINE copies belonging to 10 families and subfamilies on a base pair level. Solanaceae SINEs are dispersed over all chromosomes with enrichments in distal regions. Depending on the genome assemblies and gene predictions, 30% of all SINE copies are associated with genes, particularly frequent in introns and untranslated regions (UTRs). The close association with genes is family specific. More than 10% of all genes annotated in the Solanaceae species investigated contain at least one SINE insertion, and we found genes harbouring up to 16 SINE copies. We demonstrate the involvement of SINEs in gene and genome evolution including the donation of splice sites, start and stop codons and exons to genes, enlargement of introns and UTRs, generation of tandem-like duplications and transduction of adjacent sequence regions.

CoExpNetViz: Comparative Co-Expression Networks Construction and Visualization Tool

MOTIVATION

Comparative transcriptomics is a common approach in functional gene discovery efforts. It allows for finding conserved co-expression patterns between orthologous genes in closely related plant species, suggesting that these genes potentially share similar function and regulation. Several efficient co-expression-based tools have been commonly used in plant research but most of these pipelines are limited to data from model systems, which greatly limit their utility. Moreover, in addition, none of the existing pipelines allow plant researchers to make use of their own unpublished gene expression data for performing a comparative co-expression analysis and generate multi-species co-expression networks.

RESULTS

We introduce CoExpNetViz, a computational tool that uses a set of query or "bait" genes as an input (chosen by the user) and a minimum of one pre-processed gene expression dataset. The CoExpNetViz algorithm proceeds in three main steps; (i) for every bait gene submitted, co-expression values are calculated using mutual information and Pearson correlation coefficients, (ii) non-bait (or target) genes are grouped based on cross-species orthology, and (iii) output files are generated and results can be visualized as network graphs in Cytoscape.

AVAILABILITY

The CoExpNetViz tool is freely available both as a PHP web server (link: http://bioinformatics.psb.ugent.be/webtools/coexpr/) (implemented in C++) and as a Cytoscape plugin (implemented in Java). Both versions of the CoExpNetViz tool support LINUX and Windows platforms.

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.