Bayesian Modeling Approaches for Temporal Dynamics in RNA-seq Data

Analysis of differential expression has been a central role to address the variety of biological questions in the manner to characterize abnormal patterns of cellular and molecular functions for last decades. To date, identification of differentially expressed genes and isoforms has been more intensively focused on temporal dynamics over a series of time points. Bayesian strategies have been successfully employed to uncover the complexity of biological interest with the methodological and analytical perspectives for the various platforms of high-throughput data, for instance, methods in differential expression analysis and network modules in transcriptome data, peak-callers in ChipSeq data, target prediction in microRNA data and meta-methods between different platforms. In this chapter, we will discuss how our methodological works based on Bayesian models address important questions to arise in the architecture of temporal dynamics in RNA-seq data

Functional exploration of co-expression networks identifies a nexus for modulating protein and citric acid titres in Aspergillus niger submerged culture

Background: Filamentous fungal cell factories are used to produce numerous proteins, enzymes, and organic acids. Protein secretion and filamentous growth are tightly coupled at the hyphal tip. Additionally, both these processes require ATP and amino acid precursors derived from the citric acid cycle. Despite this interconnection of organic acid production and protein secretion/filamentous growth, few studies in fungi have identified genes which may concomitantly impact all three processes. Results: We applied a novel screen of a global co-expression network in the cell factory Aspergillus niger to identify candidate genes which may concomitantly impact macromorphology, and protein/organic acid fermentation. This identified genes predicted to encode the Golgi localized ArfA GTPase activating protein (GAP, AgeB), and ArfA guanine nucleotide exchange factors (GEFs SecG and GeaB) to be co-expressed with citric acid cycle genes. Consequently, we used CRISPR-based genome editing to place the titratable Tet-on expression system upstream of ageB, secG, and geaB in A. niger. Functional analysis revealed that ageB and geaB are essential whereas secG was dispensable for early filamentous growth. Next, gene expression was titrated during submerged cultivations under conditions for either protein or organic acid production. ArfA regulators played varied and culture-dependent roles on pellet formation. Notably, ageB or geaB expression levels had major impacts on protein secretion, whereas secG was dispensable. In contrast, reduced expression of each predicted ArfA regulator resulted in an absence of citric acid in growth media. Finally, titrated expression of either GEFs resulted in an increase in oxaloacetic acid concentrations in supernatants. Conclusion: Our data suggest that the Golgi may play an underappreciated role in modulating organic acid titres during industrial applications, and that this is SecG, GeaB and AgeB dependent in A. niger. These data may lead to novel avenues for strain optimization in filamentous fungi for improved protein and organic acid titres.TU Berlin, Open-Access-Mittel - 201

MicroRNA co-expression networks exhibit increased complexity in pancreatic ductal compared to Vater’s papilla adenocarcinoma

iRNA expression abnormalities in adenocarcinoma arising from pancreatic ductal system (PDAC) and Vater’s papilla (PVAC) could be associated with distinctive pathologic features and clinical cancer behaviours. Our previous miRNA expression profiling data on PDAC (n=9) and PVAC (n=4) were revaluated to define differences/ similarities in miRNA expression patterns. Afterwards, in order to uncover target genes and core signalling pathways regulated by specific miRNAs in these two tumour entities, miRNA interaction networks were wired for each tumour entity, and experimentally validated target genes underwent pathways enrichment analysis. One hundred and one miRNAs were altered, mainly over-expressed, in PDAC samples. Twenty-six miRNAs were deregulated in PVAC samples, where more miRNAs were down-expressed in tumours compared to normal tissues. Four miRNAs were significantly altered in both subgroups of patients, while 27 miRNAs were differentially expressed between PDAC and PVAC. Although miRNA interaction networks were more complex and dense in PDAC than in PVAC, pathways enrichment analysis uncovered a functional overlapping between PDAC and PVAC. However, shared signalling events were influenced by different miRNA and/or genes in the two tumour entities. Overall, specific miRNA expression patterns were involved in the regulation of a limited core signalling pathways in the biology landscape of PDAC and PVAC

Detecting Differentially Co-Expressed Gene Modules Via The Edge-Count Test

Background Gene expression profiling by microarray has been used to uncover molecular variations in many different diseases. Complementary to conventional differential expression analysis, differential co-expression analysis can identify gene markers from the systematic and granular level. There are three aspects for differential co-expression network analysis, including the network global topological comparison, differential co-expression cluster identification, and differential co-expressed genes and gene pair identification. To date, most of the methods available still rely on Pearson’s correlation coefficient despite its nonlinear insensitivity. Results Here we present an approach that is robust to nonlinearity by using the edge-count test for differential co-expression analysis. The performance of the new approach was tested with synthetic data and found to have significant results. For real data, we used a human cervical cancer data set prepared from 29 pairs of cervical tumor and matched normal tissue samples. Hierarchical cluster analysis resulted in the identification of clusters containing differentially co-expressed genes associated with the regulation of cervical cancer. Conclusion The proposed approach targets all different types of differential co-expression and it is sensitive to nonlinear relations. It is easy to implement and can be applied to any sequencing data to identify gene co-expression differences between multiple conditions

Cardiac expression of the microsomal triglyceride transport protein protects the heart function during ischemia

Aims: The microsomal triglyceride transport protein (MTTP) is critical for assembly and secretion of apolipoprotein B (apoB)-containing lipoproteins and is most abundant in the liver and intestine. Surprisingly, MTTP is also expressed in the heart. Here we tested the functional relevance of cardiac MTTP expression. Materials and methods: We combined clinical studies, advanced expression analysis of human heart biopsies and analyses in genetically modified mice lacking cardiac expression of the MTTP-A isoform of MTTP. Results: Our results indicate that lower cardiac MTTP expression in humans is associated with structural and perfusion abnormalities in patients with ischemic heart disease. MTTP-A deficiency in mice heart does not affect total MTTP expression, activity or lipid concentration in the heart. Despite this, MTTP-A deficient mice displayed impaired cardiac function after a myocardial infarction. Expression analysis of MTTP indicates that MTTP expression is linked to cardiac function and responses in the heart. Conclusions: Our results indicate that MTTP may play an important role for the heart function in conjunction to ischemic events