Predicting the functional repertoire of an organism from unassembled RNA-seq data.

BACKGROUND: The annotation of biomolecular functions is an essential step in the analysis of newly sequenced organisms. Usually, the functions are inferred from predicted genes on the genome using homology search techniques. A high quality genomic sequence is an important prerequisite which, however, is difficult to achieve for certain organisms, such as hybrids or organisms with a large genome. For functional analysis it is also possible to use a de novo transcriptome assembly but the computational requirements can be demanding. Up to now, it is unclear how much of the functional repertoire of an organism can be reliably predicted from unassembled RNA-seq short reads alone. RESULTS: We have conducted a study to investigate to what degree it is possible to reconstruct the functional profile of an organism from unassembled transcriptome data. We simulated the de novo prediction of biomolecular functions for Arabidopsis thaliana using a comprehensive RNA-seq data set. We evaluated the prediction performance using several homology search methods in combination with different evidence measures. For the decision on the presence or absence of a particular function under noisy conditions we propose a statistical mixture model enabling unsupervised estimation of a detection threshold. Our results indicate that the prediction of the biomolecular functions from the KEGG database is possible with a high sensitivity up to 94 percent. In this setting, the application of the mixture model for automatic threshold calibration allowed the reduction of the falsely predicted functions down to 4 percent. Furthermore, we found that our statistical approach even outperforms the prediction from a de novo transcriptome assembly. CONCLUSION: The analysis of an organism's transcriptome can provide a solid basis for the prediction of biomolecular functions. Using RNA-seq short reads directly, the functional profile of an organism can be reconstructed in a computationally efficient way to provide a draft annotation in cases where the classical genome-based approaches cannot be applied.

Soluble phenylpropanoids are involved in the defense response of Arabidopsis against Verticillium longisporum.

Verticillium longisporum is a soil-borne vascular pathogen causing economic loss in rape. Using the model plant Arabidopsis this study analyzed metabolic changes upon fungal infection in order to identify possible defense strategies of Brassicaceae against this fungus. Metabolite fingerprinting identified infection-induced metabolites derived from the phenylpropanoid pathway. Targeted analysis confirmed the accumulation of sinapoyl glucosides, coniferin, syringin and lignans in leaves from early stages of infection on. At later stages, the amounts of amino acids increased. To test the contribution of the phenylpropanoid pathway, mutants in the pathway were analyzed. The sinapate-deficient mutant fah1-2 showed stronger infection symptoms than wild-type plants, which is most likely due to the lack of sinapoyl esters. Moreover, the coniferin accumulating transgenic plant UGT72E2-OE was less susceptible. Consistently, sinapoyl glucose, coniferyl alcohol and coniferin inhibited fungal growth and melanization in vitro, whereas sinapyl alcohol and syringin did not. The amount of lignin was not significantly altered supporting the notion that soluble derivatives of the phenylpropanoid pathway contribute to defense. These data show that soluble phenylpropanoids are important for the defense response of Arabidopsis against V. longisporum and that metabolite fingerprinting is a valuable tool to identify infection-relevant metabolic markers.

Verticillium transcription activator of adhesion Vta2 suppresses microsclerotia formation and is required for systemic infection of plant roots.

Six transcription regulatory genes of the Verticillium plant pathogen, which reprogrammed nonadherent budding yeasts for adhesion, were isolated by a genetic screen to identify control elements for early plant infection. Verticillium transcription activator of adhesion Vta2 is highly conserved in filamentous fungi but not present in yeasts. The Magnaporthe grisea ortholog conidiation regulator Con7 controls the formation of appressoria which are absent in Verticillium species. Vta2 was analyzed by using genetics, cell biology, transcriptomics, secretome proteomics and plant pathogenicity assays. Nuclear Vta2 activates the expression of the adhesin-encoding yeast flocculin genes FLO1 and FLO11. Vta2 is required for fungal growth of Verticillium where it is a positive regulator of conidiation. Vta2 is mandatory for accurate timing and suppression of microsclerotia as resting structures. Vta2 controls expression of 270 transcripts, including 10 putative genes for adhesins and 57 for secreted proteins. Vta2 controls the level of 125 secreted proteins, including putative adhesins or effector molecules and a secreted catalase-peroxidase. Vta2 is a major regulator of fungal pathogenesis, and controls host-plant root infection and H2 O2 detoxification. Verticillium impaired in Vta2 is unable to colonize plants and induce disease symptoms. Vta2 represents an interesting target for controlling the growth and development of these vascular pathogens.

Breaking the silence: protein stabilization uncovers silenced biosynthetic gene clusters in the fungus Aspergillus nidulans.

The genomes of filamentous fungi comprise numerous putative gene clusters coding for the biosynthesis of chemically and structurally diverse secondary metabolites (SMs), which are rarely expressed under laboratory conditions. Previous approaches to activate these genes were based primarily on artificially targeting the cellular protein synthesis apparatus. Here, we applied an alternative approach of genetically impairing the protein degradation apparatus of the model fungus Aspergillus nidulans by deleting the conserved eukaryotic csnE/CSN5 deneddylase subunit of the COP9 signalosome. This defect in protein degradation results in the activation of a previously silenced gene cluster comprising a polyketide synthase gene producing the antibiotic 2,4-dihydroxy-3-methyl-6-(2-oxopropyl)benzaldehyde (DHMBA). The csnE/CSN5 gene is highly conserved in fungi, and therefore, the deletion is a feasible approach for the identification of new SMs.

Arabidopsis mutants of sphingolipid fatty acid α-hydroxylases accumulate ceramides and salicylates.

In Arabidopsis, the fatty acid moiety of sphingolipids is mainly α-hydroxylated. The consequences of a reduction in this modification were analysed. Mutants of both Fatty Acid Hydroxylase genes (AtFAH1 and AtFAH2) were analysed for sphingolipid profiles. To elucidate further consequences of the mutations, metabolic analyses were performed and the influence on pathogen defence was determined. Ceramide and glucosylceramide profiles of double-mutant plants showed a reduction in sphingolipids with α-hydroxylated fatty acid moieties, and an accumulation of sphingolipids without these moieties. In addition, the free trihydroxylated long-chain bases and ceramides were increased by five- and ten-fold, respectively, whereas the amount of glucosylceramides was decreased by 25%. Metabolite analysis of the double mutant revealed salicylates as enriched metabolites. Infection experiments supported the metabolic changes, as the double mutant showed an enhanced disease-resistant phenotype for infection with the obligate biotrophic pathogen Golovinomyces cichoracearum. In summary, these results suggest that fatty acid hydroxylation of ceramides is important for the biosynthesis of complex sphingolipids. Its absence leads to the accumulation of long-chain bases and ceramides as their precursors. This increases salicylate levels and resistance towards obligate biotrophic fungal pathogens, confirming a role of sphingolipids in salicylic acid-dependent defence reactions.

MarVis-Filter: ranking, filtering, adduct and isotope correction of mass spectrometry data.

Statistical ranking, filtering, adduct detection, isotope correction, and molecular formula calculation are essential tasks in processing mass spectrometry data in metabolomics studies. In order to obtain high-quality data sets, a framework which incorporates all these methods is required. We present the MarVis-Filter software, which provides well-established and specialized methods for processing mass spectrometry data. For the task of ranking and filtering multivariate intensity profiles, MarVis-Filter provides the ANOVA and Kruskal-Wallis tests with adjustment for multiple hypothesis testing. Adduct and isotope correction are based on a novel algorithm which takes the similarity of intensity profiles into account and allows user-defined ionization rules. The molecular formula calculation utilizes the results of the adduct and isotope correction. For a comprehensive analysis, MarVis-Filter provides an interactive interface to combine data sets deriving from positive and negative ionization mode. The software is exemplarily applied in a metabolic case study, where octadecanoids could be identified as markers for wounding in plants.

A 20-kb lineage-specific genomic region tames virulence in pathogenic amphidiploid Verticillium longisporum.

Amphidiploid fungal Verticillium longisporum strains Vl43 and Vl32 colonize the plant host Brassica napus but differ in their ability to cause disease symptoms. These strains represent two V. longisporum lineages derived from different hybridization events of haploid parental Verticillium strains. Vl32 and Vl43 carry same-sex mating-type genes derived from both parental lineages. Vl32 and Vl43 similarly colonize and penetrate plant roots, but asymptomatic Vl32 proliferation in planta is lower than virulent Vl43. The highly conserved Vl43 and Vl32 genomes include less than 1% unique genes, and the karyotypes of 15 or 16 chromosomes display changed genetic synteny due to substantial genomic reshuffling. A 20 kb Vl43 lineage-specific (LS) region apparently originating from the Verticillium dahliae-related ancestor is specific for symptomatic Vl43 and encodes seven genes, including two putative transcription factors. Either partial or complete deletion of this LS region in Vl43 did not reduce virulence but led to induction of even more severe disease symptoms in rapeseed. This suggests that the LS insertion in the genome of symptomatic V. longisporum Vl43 mediates virulence-reducing functions, limits damage on the host plant, and therefore tames Vl43 from being even more virulent.

Verticillium longisporum Elicits Media-Dependent Secretome Responses With Capacity to Distinguish Between Plant-Related Environments.

Verticillia cause a vascular wilt disease affecting a broad range of economically valuable crops. The fungus enters its host plants through the roots and colonizes the vascular system. It requires extracellular proteins for a successful plant colonization. The exoproteomes of the allodiploid Verticillium longisporum upon cultivation in different media or xylem sap extracted from its host plant Brassica napus were compared. Secreted fungal proteins were identified by label free liquid chromatography-tandem mass spectrometry screening. V. longisporum induced two main secretion patterns. One response pattern was elicited in various non-plant related environments. The second pattern includes the exoprotein responses to the plant-related media, pectin-rich simulated xylem medium and pure xylem sap, which exhibited similar but additional distinct features. These exoproteomes include a shared core set of 221 secreted and similarly enriched fungal proteins. The pectin-rich medium significantly induced the secretion of 143 proteins including a number of pectin degrading enzymes, whereas xylem sap triggered a smaller but unique fungal exoproteome pattern with 32 enriched proteins. The latter pattern included proteins with domains of known pathogenicity factors, metallopeptidases and carbohydrate-active enzymes. The most abundant proteins of these different groups are the necrosis and ethylene inducing-like proteins Nlp2 and Nlp3, the cerato-platanin proteins Cp1 and Cp2, the metallopeptidases Mep1 and Mep2 and the carbohydrate-active enzymes Gla1, Amy1 and Cbd1. Their pathogenicity contribution was analyzed in the haploid parental strain V. dahliae. Deletion of the majority of the corresponding genes caused no phenotypic changes during ex planta growth or invasion and colonization of tomato plants. However, we discovered that the MEP1, NLP2, and NLP3 deletion strains were compromised in plant infections. Overall, our exoproteome approach revealed that the fungus induces specific secretion responses in different environments. The fungus has a general response to non-plant related media whereas it is able to fine-tune its exoproteome in the presence of plant material. Importantly, the xylem sap-specific exoproteome pinpointed Nlp2 and Nlp3 as single effectors required for successful V. dahliae colonization.

Evaluation of molecular subtypes and clonal selection during establishment of patient-derived tumor xenografts from gastric adenocarcinoma.

Patient-derived xenografts (PDX) have emerged as an important translational research tool for understanding tumor biology and enabling drug efficacy testing. They are established by transfer of patient tumor into immune compromised mice with the intent of using them as Avatars; operating under the assumption that they closely resemble patient tumors. In this study, we established 27 PDX from 100 resected gastric cancers and studied their fidelity in histological and molecular subtypes. We show that the established PDX preserved histology and molecular subtypes of parental tumors. However, in depth investigation of the entire cohort revealed that not all histological and molecular subtypes are established. Also, for the established PDX models, genetic changes are selected at early passages and rare subclones can emerge in PDX. This study highlights the importance of considering the molecular and evolutionary characteristics of PDX for a proper use of such models, particularly for Avatar trials.

Author Correction: Evaluation of molecular subtypes and clonal selection during establishment of patient-derived tumor xenografts from gastric adenocarcinoma.

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

Patient derived renal cell carcinoma xenografts exhibit distinct sensitivity patterns in response to antiangiogenic therapy and constitute a suitable tool for biomarker development.

Systemic treatment is necessary for one third of patients with renal cell carcinoma. No valid biomarker is currently available to tailor personalized therapy. In this study we established a representative panel of patient derived xenograft (PDX) mouse models from patients with renal cell carcinomas and determined serum levels of high mobility group B1 (HMGB1) protein under treatment with sunitinib, pazopanib, sorafenib, axitinib, temsirolimus and bevacizumab. Serum HMGB1 levels were significantly higher in a subset of the PDX collection, which exhibited slower tumor growth during subsequent passages than tumors with low HMGB1 serum levels. Pre-treatment PDX serum HMGB1 levels also correlated with response to systemic treatment: PDX models with high HMGB1 levels predicted response to bevacizumab. Taken together, we provide for the first time evidence that the damage associated molecular pattern biomarker HMGB1 can predict response to systemic treatment with bevacizumab. Our data support the future evaluation of HMGB1 as a predictive biomarker for bevacizumab sensitivity in patients with renal cell carcinoma.

MarVis-Pathway: integrative and exploratory pathway analysis of non-targeted metabolomics data.

A central aim in the evaluation of non-targeted metabolomics data is the detection of intensity patterns that differ between experimental conditions as well as the identification of the underlying metabolites and their association with metabolic pathways. In this context, the identification of metabolites based on non-targeted mass spectrometry data is a major bottleneck. In many applications, this identification needs to be guided by expert knowledge and interactive tools for exploratory data analysis can significantly support this process. Additionally, the integration of data from other omics platforms, such as DNA microarray-based transcriptomics, can provide valuable hints and thereby facilitate the identification of metabolites via the reconstruction of related metabolic pathways. We here introduce the MarVis-Pathway tool, which allows the user to identify metabolites by annotation of pathways from cross-omics data. The analysis is supported by an extensive framework for pathway enrichment and meta-analysis. The tool allows the mapping of data set features by ID, name, and accurate mass, and can incorporate information from adduct and isotope correction of mass spectrometry data. MarVis-Pathway was integrated in the MarVis-Suite (http://marvis.gobics.de), which features the seamless highly interactive filtering, combination, clustering, and visualization of omics data sets. The functionality of the new software tool is illustrated using combined mass spectrometry and DNA microarray data. This application confirms jasmonate biosynthesis as important metabolic pathway that is upregulated during the wound response of Arabidopsis plants.

Meta-analysis of pathway enrichment: combining independent and dependent omics data sets.

A major challenge in current systems biology is the combination and integrative analysis of large data sets obtained from different high-throughput omics platforms, such as mass spectrometry based Metabolomics and Proteomics or DNA microarray or RNA-seq-based Transcriptomics. Especially in the case of non-targeted Metabolomics experiments, where it is often impossible to unambiguously map ion features from mass spectrometry analysis to metabolites, the integration of more reliable omics technologies is highly desirable. A popular method for the knowledge-based interpretation of single data sets is the (Gene) Set Enrichment Analysis. In order to combine the results from different analyses, we introduce a methodical framework for the meta-analysis of p-values obtained from Pathway Enrichment Analysis (Set Enrichment Analysis based on pathways) of multiple dependent or independent data sets from different omics platforms. For dependent data sets, e.g. obtained from the same biological samples, the framework utilizes a covariance estimation procedure based on the nonsignificant pathways in single data set enrichment analysis. The framework is evaluated and applied in the joint analysis of Metabolomics mass spectrometry and Transcriptomics DNA microarray data in the context of plant wounding. In extensive studies of simulated data set dependence, the introduced correlation could be fully reconstructed by means of the covariance estimation based on pathway enrichment. By restricting the range of p-values of pathways considered in the estimation, the overestimation of correlation, which is introduced by the significant pathways, could be reduced. When applying the proposed methods to the real data sets, the meta-analysis was shown not only to be a powerful tool to investigate the correlation between different data sets and summarize the results of multiple analyses but also to distinguish experiment-specific key pathways.

Integrative study of Arabidopsis thaliana metabolomic and transcriptomic data with the interactive MarVis-Graph software.

State of the art high-throughput technologies allow comprehensive experimental studies of organism metabolism and induce the need for a convenient presentation of large heterogeneous datasets. Especially, the combined analysis and visualization of data from different high-throughput technologies remains a key challenge in bioinformatics. We present here the MarVis-Graph software for integrative analysis of metabolic and transcriptomic data. All experimental data is investigated in terms of the full metabolic network obtained from a reference database. The reactions of the network are scored based on the associated data, and sub-networks, according to connected high-scoring reactions, are identified. Finally, MarVis-Graph scores the detected sub-networks, evaluates them by means of a random permutation test and presents them as a ranked list. Furthermore, MarVis-Graph features an interactive network visualization that provides researchers with a convenient view on the results. The key advantage of MarVis-Graph is the analysis of reactions detached from their pathways so that it is possible to identify new pathways or to connect known pathways by previously unrelated reactions. The MarVis-Graph software is freely available for academic use and can be downloaded at: http://marvis.gobics.de/marvis-graph.