Beyond element-wise interactions: identifying complex interactions in biological processes

Background: Biological processes typically involve the interactions of a number of elements (genes, cells) acting on each others. Such processes are often modelled as networks whose nodes are the elements in question and edges pairwise relations between them (transcription, inhibition). But more often than not, elements actually work cooperatively or competitively to achieve a task. Or an element can act on the interaction between two others, as in the case of an enzyme controlling a reaction rate. We call “complex” these types of interaction and propose ways to identify them from time-series observations. Methodology: We use Granger Causality, a measure of the interaction between two signals, to characterize the influence of an enzyme on a reaction rate. We extend its traditional formulation to the case of multi-dimensional signals in order to capture group interactions, and not only element interactions. Our method is extensively tested on simulated data and applied to three biological datasets: microarray data of the Saccharomyces cerevisiae yeast, local field potential recordings of two brain areas and a metabolic reaction. Conclusions: Our results demonstrate that complex Granger causality can reveal new types of relation between signals and is particularly suited to biological data. Our approach raises some fundamental issues of the systems biology approach since finding all complex causalities (interactions) is an NP hard problem

Translational web robots for pathogen genome analysis

4 page(s

Identification and characterization of extensive intra-molecular associations between 3′-UTRs and their ORFs

During eukaryotic translation, mRNAs may form intra-molecular interactions between distant domains. The 5′-cap and the polyA tail were shown to interact through their associated proteins, and this can induce physical compaction of the mRNA in vitro. However, the stability of this intra-molecular association in translating mRNAs and whether additional contacts exist in vivo are largely unknown. To explore this, we applied a novel approach in which several endogenous polysomal mRNAs from Saccharomyces cerevisiae were cleaved near their stop codon and the resulting 3′-UTR fragments were tested either for co-sedimentation or co-immunoprecipitation (co-IP) with their ORFs. In all cases a significant fraction of the 3′-UTR fragments sedimented similarly to their ORF-containing fragments, yet the extent of co-sedimentation differed between mRNAs. Similar observations were obtained by a co-IP assay. Interestingly, various treatments that are expected to interfere with the cap to polyA interactions had no effect on the co-sedimentation pattern. Moreover, the 3′-UTR appeared to co-sediment with different regions from within the ORF. Taken together, these results indicate extensive physical associations between 3′-UTRs and their ORFs that vary between genes. This implies that polyribosomal mRNAs are in a compact configuration in vivo

Analyzing the genomic variation of microbial cell factories in the era of “New Biotechnology”

The application of genome-scale technologies, both experimental and in silico, to industrial biotechnology has allowed improving the conversion of biomass-derived feedstocks to chemicals, materials and fuels through microbial fermentation. In particular, due to rapidly decreasing costs and its suitability for identifying the genetic determinants of a phenotypic trait of interest, whole genome sequencing is expected to be one of the major driving forces in industrial biotechnology in the coming years. We present some of the recent studies that have successfully applied high-throughput sequencing technologies for finding the underlying molecular mechanisms for (a) improved carbon source utilization, (b) increased product formation, and (c) stress tolerance. We also discuss the strengths and weaknesses of different strategies for mapping industrially relevant genotype-to-phenotype links including exploiting natural diversity in natural isolates or crosses between isolates, classical mutagenesis and evolutionary engineering

Stimulation of translation by human Unr requires cold shock domains 2 and 4, and correlates with poly(A) binding protein interaction

The RNA binding protein Unr, which contains five cold shock domains, has several specific roles in post-transcriptional control of gene expression. It can act as an activator or inhibitor of translation initiation, promote mRNA turnover, or stabilise mRNA. Its role depends on the mRNA and other proteins to which it binds, which includes cytoplasmic poly(A) binding protein 1 (PABP1). Since PABP1 binds to all polyadenylated mRNAs, and is involved in translation initiation by interaction with eukaryotic translation initiation factor 4G (eIF4G), we investigated whether Unr has a general role in translational control. We found that Unr strongly stimulates translation in vitro, and mutation of cold shock domains 2 or 4 inhibited its translation activity. The ability of Unr and its mutants to stimulate translation correlated with its ability to bind RNA, and to interact with PABP1. We found that Unr stimulated the binding of PABP1 to mRNA, and that Unr was required for the stable interaction of PABP1 and eIF4G in cells. siRNA-mediated knockdown of Unr reduced the overall level of cellular translation in cells, as well as that of cap-dependent and IRES-dependent reporters. These data describe a novel role for Unr in regulating cellular gene expression

Time-dependent increase in ribosome processivity

We created a novel tripartite reporter RNA to separately and simultaneously examine ribosome translation rates at the 5′- and 3′-ends of a large open reading frame (ORF) in vitro in HeLa cell lysates. The construct contained Renilla luciferase (RLuc), β-galactosidase and firefly luciferase (FLuc) ORFs linked in frame and separated by a viral peptide sequence that causes cotranslational scission of emerging peptide chains. The length of the ORF contributed to low ribosome processivity, a low number of initiating ribosomes completing translation of the entire ORF. We observed a time-dependent increase in FLuc production rate that was dependent on a poly(A) tail and poly(A)-binding protein, but was independent of eIF4F function. Stimulation of FLuc production occurred earlier on shorter RNA templates. Cleavage of eIF4G at times after ribosome loading on templates occurred did not cause immediate cessation of 5′-RLuc translation; rather, a delay was observed that shortened when shorter templates were translated. Electron microscopic analysis of polysome structures in translation lysates revealed a time-dependent increase in ribosome packing and contact that correlated with increased processivity on the FLuc ORF. The results suggest that ORF transit combined with PABP function contribute to interactions between ribosomes that increase or sustain processivity on long ORFs

Into the unknown: expression profiling without genome sequence information in CHO by next generation sequencing

The arrival of next-generation sequencing (NGS) technologies has led to novel opportunities for expression profiling and genome analysis by utilizing vast amounts of short read sequence data. Here, we demonstrate that expression profiling in organisms lacking any genome or transcriptome sequence information is feasible by combining Illumina’s mRNA-seq technology with a novel bioinformatics pipeline that integrates assembled and annotated Chinese hamster ovary (CHO) sequences with information derived from related organisms. We applied this pipeline to the analysis of CHO cells which were chosen as a model system owing to its relevance in the production of therapeutic proteins. Specifically, we analysed CHO cells undergoing butyrate treatment which is known to affect cell cycle regulation and to increase the specific productivity of recombinant proteins. By this means, we identified sequences for >13 000 CHO genes which added sequence information of ∼5000 novel genes to the CHO model. More than 6000 transcript sequences are predicted to be complete, as they covered >95% of the corresponding mouse orthologs. Detailed analysis of selected biological functions such as DNA replication and cell cycle control, demonstrated the potential of NGS expression profiling in organisms without extended genome sequence to improve both data quantity and quality

Next-generation insights into regulatory T cells: expression profiling and FoxP3 occupancy in Human

Regulatory T-cells (Treg) play an essential role in the negative regulation of immune answers by developing an attenuated cytokine response that allows suppressing proliferation and effector function of T-cells (CD4+ Th). The transcription factor FoxP3 is responsible for the regulation of many genes involved in the Treg gene signature. Its ablation leads to severe immune deficiencies in human and mice. Recent developments in sequencing technologies have revolutionized the possibilities to gain insights into transcription factor binding by ChiP-seq and into transcriptome analysis by mRNA-seq. We combine FoxP3 ChiP-seq and mRNA-seq in order to understand the transcriptional differences between primary human CD4+ T helper and regulatory T-cells, as well as to study the role of FoxP3 in generating those differences. We show, that mRNA-seq allows analyzing the transcriptomal landscape of T-cells including the expression of specific splice variants at much greater depth than previous approaches, whereas 50% of transcriptional regulation events have not been described before by using diverse array technologies. We discovered splicing patterns like the expression of a kinase-dead isoform of IRAK1 upon T-cell activation. The immunoproteasome is up-regulated in both Treg and CD4+ Th cells upon activation, whereas the ‘standard’ proteasome is up-regulated in Tregs only upon activation

Nucleotide Discrimination with DNA Immobilized in the MspA Nanopore

Nanopore sequencing has the potential to become a fast and low-cost DNA sequencing platform. An ionic current passing through a small pore would directly map the sequence of single stranded DNA (ssDNA) driven through the constriction. The pore protein, MspA, derived from Mycobacterium smegmatis, has a short and narrow channel constriction ideally suited for nanopore sequencing. To study MspA's ability to resolve nucleotides, we held ssDNA within the pore using a biotin-NeutrAvidin complex. We show that homopolymers of adenine, cytosine, thymine, and guanine in MspA exhibit much larger current differences than in α-hemolysin. Additionally, methylated cytosine is distinguishable from unmethylated cytosine. We establish that single nucleotide substitutions within homopolymer ssDNA can be detected when held in MspA's constriction. Using genomic single nucleotide polymorphisms, we demonstrate that single nucleotides within random DNA can be identified. Our results indicate that MspA has high signal-to-noise ratio and the single nucleotide sensitivity desired for nanopore sequencing devices