Mining metabolic pathways through gene expression

Multi-omics network model reveals key genes associated with p-coumaric acid stress response in an industrial yeast strain

The production of ethanol from lignocellulosic sources presents increasingly difficult issues for the global biofuel scenario, leading to increased production costs of current second-generation (2G) ethanol when compared to first-generation (1G) plants. Among the setbacks encountered in industrial processes, the presence of chemical inhibitors from pre-treatment processes severely hinders the potential of yeasts in producing ethanol at peak efficiency. However, some industrial yeast strains have, either naturally or artificially, higher tolerance levels to these compounds. Such is the case of S. cerevisiae SA-1, a Brazilian fuel ethanol industrial strain that has shown high resistance to inhibitors produced by the pre-treatment of cellulosic complexes. Our study focuses on the characterization of the transcriptomic and physiological impact of an inhibitor of this type, p-coumaric acid (pCA), on this strain under chemostat cultivation via RNAseq and quantitative physiological data. It was found that strain SA-1 tend to increase ethanol yield and production rate while decreasing biomass yield when exposed to pCA, in contrast to pCA-susceptible strains, which tend to decrease their ethanol yield and fermentation efficiency when exposed to this substance. This suggests increased metabolic activity linked to mitochondrial and peroxisomal processes. The transcriptomic analysis also revealed a plethora of differentially expressed genes located in co-expressed clusters that are associated with changes in biological pathways linked to biosynthetic and energetical processes. Furthermore, it was also identified 20 genes that act as interaction hubs for these clusters, while also having association with altered pathways and changes in metabolic outputs, potentially leading to the discovery of novel targets for metabolic engineering toward a more robust industrial yeast strain.

Extended mining of the oil biosynthesis pathway in biofuel plant Jatropha curcas by combined analysis of transcriptome and gene interactome data

BACKGROUND

Jatropha curcas L. is an important non-edible oilseed crop with a promising future in biodiesel production. However, little is known about the molecular biology of oil biosynthesis in this plant when compared with other established oilseed crops, resulting in the absence of agronomically improved varieties of Jatropha. To extensively discover the potentially novel genes and pathways associated with the oil biosynthesis in J. curcas, new strategy other than homology alignment is on the demand.

RESULTS

In this study, we proposed a multi-step computational framework that integrates transcriptome and gene interactome data to predict functional pathways in non-model organisms in an extended process, and applied it to study oil biosynthesis pathway in J. curcas. Using homologous mapping against Arabidopsis and transcriptome profile analysis, we first constructed protein-protein interaction (PPI) and co-expression networks in J. curcas. Then, using the homologs of Arabidopsis oil-biosynthesis-related genes as seeds, we respectively applied two algorithm models, random walk with restart (RWR) in PPI network and negative binomial distribution (NBD) in co-expression network, to further extend oil-biosynthesis-related pathways and genes in J. curcas. At last, using k-nearest neighbors (KNN) algorithm, the predicted genes were further classified into different sub-pathways according to their possible functional roles.

CONCLUSIONS

Our method exhibited a highly efficient way of mining the extended oil biosynthesis pathway of J. curcas. Overall, 27 novel oil-biosynthesis-related gene candidates were predicted and further assigned to 5 sub-pathways. These findings can help better understanding of the oil biosynthesis pathway of J. curcas, as well as paving the way for the following J. curcas breeding application.

Coupled analysis of transcriptome and BCR mutations reveals role of OXPHOS in affinity maturation

Antigen-activated B cells diversify variable regions of B cell antigen receptors by somatic hypermutation in germinal centers (GCs). The positive selection of GC B cells that acquire high-affinity mutations enables antibody affinity maturation. In spite of considerable progress, the genomic states underlying this process remain to be elucidated. Single-cell RNA sequencing and topic modeling revealed increased expression of the oxidative phosphorylation (OXPHOS) module in GC B cells undergoing mitoses. Coupled analysis of somatic hypermutation in immunoglobulin heavy chain (Igh) variable gene regions showed that GC B cells acquiring higher-affinity mutations had further elevated expression of OXPHOS genes. Deletion of mitochondrial Cox10 in GC B cells resulted in reduced cell division and impaired positive selection. Correspondingly, augmentation of OXPHOS activity with oltipraz promoted affinity maturation. We propose that elevated OXPHOS activity promotes B cell clonal expansion and also positive selection by tuning cell division times.

Biological signaling pathways and potential mathematical network representations: biological discovery through optimization

Establishing the role that different genes play in the development of cancer is a daunting task. A step toward this end is the detection of genes that are important in the illness from high-throughput biological experiments. Furthermore, it is safe to say that it is highly unlikely that these show expression changes independently, even with a list of potentially important genes. A biological signaling pathway is a more plausible underlying mechanism as favored in the literature. This work attempts to build a mathematical network problem through the analysis of microarray experiments. A preselection of genes is carried out with a multiple criteria optimization framework previously published by our research group . Afterward, application of the Traveling Salesperson Problem and Minimum Spanning Tree network optimization models are proposed to identify potential signaling pathways via the most correlated path among the genes of interest. Biological evidencing is provided to assess the effectiveness of the proposed methods. The capability of our analysis strategy is also demonstrated through the undertaking of meta-analysis studies. Three important aspects are novel in this work: (1) our joint analyses of different groups of lung cancer states reveal new correlations, biologically evidenced, and previously undocumented; (2) computation of the correlation coefficients from expression differences leads to an effective use of network optimization methods; and (3) the methods yield mathematically optimal correlation structures: no other configuration is better correlated using the available information.

A route-based pathway analysis framework integrating mutation information and gene expression data

We propose a new way of analyzing biological pathways in which the analysis combines both transcriptome data and mutation information and uses the outcome to identify "routes" of aberrant pathways potentially responsible for the etiology of disease. Each pathway route is encoded as a Bayesian Network which is initialized with a sequence of conditional probabilities which are designed to encode directionality of regulatory relationships encoded in the pathways, i.e. activation and inhibition relationships. First, we demonstrate the effectiveness of our model through simulation in which the model was able to easily separate Test samples from Control samples using fictitiously perturbed pathway routes. Second, we apply our model to analyze the Breast Cancer data set, available from TCGA, against many cancer pathways available from KEGG and rank the significance of identified pathways. The outcome is consistent with what have already been reported in the literature. Third, survival analysis has been carried out on the same data set by using pathway routes as features. Overall, we envision that our model of using pathway routes for analysis can further refine the conventional ways of subtyping cancer patients as it can discover additional characteristics specific to individual's tumor.

NetPathMiner: R/Bioconductor package for network path mining through gene expression

NetPathMiner is a general framework for mining, from genome-scale networks, paths that are related to specific experimental conditions. NetPathMiner interfaces with various input formats including KGML, SBML and BioPAX files and allows for manipulation of networks in three different forms: metabolic, reaction and gene representations. NetPathMiner ranks the obtained paths and applies Markov model-based clustering and classification methods to the ranked paths for easy interpretation. NetPathMiner also provides static and interactive visualizations of networks and paths to aid manual investigation.

AVAILABILITY

The package is available through Bioconductor and from Github at http://github.com/ahmohamed/NetPathMiner.

What mRNA Abundances Can Tell us about Metabolism

Inferring decreased or increased metabolic functions from transcript profiles is at first sight a bold and speculative attempt because of the functional layers in between: proteins, enzymatic activities, and reaction fluxes. However, the growing interest in this field can easily be explained by two facts: the high quality of genome-scale metabolic network reconstructions and the highly developed technology to obtain genome-covering RNA profiles. Here, an overview of important algorithmic approaches is given by means of criteria by which published procedures can be classified. The frontiers of the methods are sketched and critical voices are being heard. Finally, an outlook for the prospects of the field is given.

Genomic analysis of the hydrocarbon-producing, cellulolytic, endophytic fungus Ascocoryne sarcoides

The microbial conversion of solid cellulosic biomass to liquid biofuels may provide a renewable energy source for transportation fuels. Endophytes represent a promising group of organisms, as they are a mostly untapped reservoir of metabolic diversity. They are often able to degrade cellulose, and they can produce an extraordinary diversity of metabolites. The filamentous fungal endophyte Ascocoryne sarcoides was shown to produce potential-biofuel metabolites when grown on a cellulose-based medium; however, the genetic pathways needed for this production are unknown and the lack of genetic tools makes traditional reverse genetics difficult. We present the genomic characterization of A. sarcoides and use transcriptomic and metabolomic data to describe the genes involved in cellulose degradation and to provide hypotheses for the biofuel production pathways. In total, almost 80 biosynthetic clusters were identified, including several previously found only in plants. Additionally, many transcriptionally active regions outside of genes showed condition-specific expression, offering more evidence for the role of long non-coding RNA in gene regulation. This is one of the highest quality fungal genomes and, to our knowledge, the only thoroughly annotated and transcriptionally profiled fungal endophyte genome currently available. The analyses and datasets contribute to the study of cellulose degradation and biofuel production and provide the genomic foundation for the study of a model endophyte system.

A renewable source of energy is a pressing global need. The biological conversion of lignocellulose to biofuels by microorganisms presents a promising avenue, but few organisms have been studied thoroughly enough to develop the genetic tools necessary for rigorous experimentation. The filamentous-fungal endophyte A. sarcoides produces metabolites when grown on a cellulose-based medium that include eight-carbon volatile organic compounds, which are potential biofuel targets. Here we use broadly applicable methods including genomics, transcriptomics, and metabolomics to explore the biofuel production of A. sarcoides. These data were used to assemble the genome into 16 scaffolds, to thoroughly annotate the cellulose-degradation machinery, and to make predictions for the production pathway for the eight-carbon volatiles. Extremely high expression of the gene swollenin when grown on cellulose highlights the importance of accessory proteins in addition to the enzymes that catalyze the breakdown of the polymers. Correlation of the production of the eight-carbon biofuel-like metabolites with the expression of lipoxygenase pathway genes suggests the catabolism of linoleic acid as the mechanism of eight-carbon compound production. This is the first fungal genome to be sequenced in the family Helotiaceae, and A. sarcoides was isolated as an endophyte, making this work also potentially useful in fungal systematics and the study of plant–fungus relationships.

Identifying neighborhoods of coordinated gene expression and metabolite profiles

In this paper we investigate how metabolic network structure affects any coordination between transcript and metabolite profiles. To achieve this goal we conduct two complementary analyses focused on the metabolic response to stress. First, we investigate the general size of any relationship between metabolic network gene expression and metabolite profiles. We find that strongly correlated transcript-metabolite profiles are sustained over surprisingly long network distances away from any target metabolite. Secondly, we employ a novel pathway mining method to investigate the structure of this transcript-metabolite relationship. The objective of this method is to identify a minimum set of metabolites which are the target of significantly correlated gene expression pathways. The results reveal that in general, a global regulation signature targeting a small number of metabolites is responsible for a large scale metabolic response. However, our method also reveals pathway specific effects that can degrade this global regulation signature and complicates the observed coordination between transcript-metabolite profiles.