Developing a systems-level understanding of gene expression

A report on the meeting 'Systems Biology: Global Regulation of Gene Expression' at the Cold Spring Harbor Laboratory, New York, USA, 28 March-1 April 2007

New gadolinium contrast agents for MRI

A collaboration between Bracco Imaging S.p.A. and Durham University allowed the work described in this thesis on the design and synthesis of new contrast agents for MRI. Significant enhancements in the relaxivity of contrast agents for MRI can be gained by increasing the complex rotational correlation time ((_TR)). Incorporating a Gd(^111) ion within a ligand structure possessing a suitably large dendritic framework, inspired the first part of this project. Thus, the periphery of a Gd-DOTA derivative was adorned with carbohydrate containing wedges. The symmetry of the mono- aqua tetra-substituted structure places the gadolinium-water vector at the centre of any tumbling motion, allowing a coherent tumbling of the macromolecule and an optimization of its rotational correlation time. The carbohydrates ensured high water solubility and favoured a large second sphere hydration contribution to the relaxivity. An increase in the hydration around the metal centre and a rapid exchange of the water molecules with the bulk solvent can also significantly increase the contrast agent efficacy, by efficiently transmitting the paramagnetic effect from the Gd(^111) to the solvent. In a second part of the work, the development was undertaken of diaqua systems based on the seven-membered heterocycle 6-amino-6-methyl- perhydro-1,4-diazepin (AMPED). The three N-positions were substituted with different phosphinate and carboxylate groups and lanthanide complexes (Eu(^111), Gd(^111), Yb(^111)) prepared and studied by multinuclear NMR methods The alkylation of the amino groups with chiral 1,5-dicarboxylate pendant arms led to complex diastereisomers, possessing different water exchange rates. The individual water exchange rates of each isomer were determined, and differed by a factor of six. Furthermore, the periphery of the isomer possessing a faster water exchange rate was adorned with carbohydrate containing wedges, and the relaxation properties studied

Ab initio genotype–phenotype association reveals intrinsic modularity in genetic networks

Microbial species express an astonishing diversity of phenotypic traits, behaviors, and metabolic capacities. However, our molecular understanding of these phenotypes is based almost entirely on studies in a handful of model organisms that together represent only a small fraction of this phenotypic diversity. Furthermore, many microbial species are not amenable to traditional laboratory analysis because of their exotic lifestyles and/or lack of suitable molecular genetic techniques. As an adjunct to experimental analysis, we have developed a computational information-theoretic framework that produces high-confidence gene–phenotype predictions using cross-species distributions of genes and phenotypes across 202 fully sequenced archaea and eubacteria. In addition to identifying the genetic basis of complex traits, our approach reveals the organization of these genes into generic preferentially co-inherited modules, many of which correspond directly to known enzymatic pathways, molecular complexes, signaling pathways, and molecular machines

Fast and systematic genome-wide discovery of conserved regulatory elements using a non-alignment based approach

We describe a powerful new approach for discovering globally conserved regulatory elements between two genomes. The method is fast, simple and comprehensive, without requiring alignments. Its application to pairs of yeasts, worms, flies and mammals yields a large number of known and novel putative regulatory elements. Many of these are validated by independent biological observations, have spatial and/or orientation biases, are co-conserved with other elements and show surprising conservation across large phylogenetic distances

Use of RNA sequencing to evaluate rheumatic disease patients

Studying the factors that control gene expression is of substantial importance for rheumatic diseases with poorly understood etiopathogenesis. In the past, gene expression microarrays have been used to measure transcript abundance on a genome-wide scale in a particular cell, tissue or organ. Microarray analysis has led to gene signatures that differentiate rheumatic diseases, and stages of a disease, as well as response to treatments. Nowadays, however, with the advent of next-generation sequencing methods, massive parallel sequencing of RNA tends to be the technology of choice for gene expression profiling, due to several advantages over microarrays, as well as for the detection of non-coding transcripts and alternative splicing events. In this review, we describe how RNA sequencing enables unbiased interrogation of the abundance and complexity of the transcriptome, and present a typical experimental workflow and bioinformatics tools that are often used for RNA sequencing analysis. We also discuss different uses of this next-generation sequencing technology to evaluate rheumatic disease patients and investigate the pathogenesis of rheumatic diseases such as rheumatoid arthritis, systemic lupus erythematosus, juvenile idiopathic arthritis and Sjögren\u27s syndrome

MOTIFATOR: detection and characterization of regulatory motifs using prokaryote transcriptome data

Summary: Unraveling regulatory mechanisms (e.g. identification of motifs in cis-regulatory regions) remains a major challenge in the analysis of transcriptome experiments. Existing applications identify putative motifs from gene lists obtained at rather arbitrary cutoff and require additional manual processing steps. Our standalone application MOTIFATOR identifies the most optimal parameters for motif discovery and creates an interactive visualization of the results. Discovered putative motifs are functionally characterized, thereby providing valuable insight in the biological processes that could be controlled by the motif.

Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and near Stop Codons

SummaryMethylation of the N6 position of adenosine (m6A) is a posttranscriptional modification of RNA with poorly understood prevalence and physiological relevance. The recent discovery that FTO, an obesity risk gene, encodes an m6A demethylase implicates m6A as an important regulator of physiological processes. Here, we present a method for transcriptome-wide m6A localization, which combines m6A-specific methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-Seq). We use this method to identify mRNAs of 7,676 mammalian genes that contain m6A, indicating that m6A is a common base modification of mRNA. The m6A modification exhibits tissue-specific regulation and is markedly increased throughout brain development. We find that m6A sites are enriched near stop codons and in 3′ UTRs, and we uncover an association between m6A residues and microRNA-binding sites within 3′ UTRs. These findings provide a resource for identifying transcripts that are substrates for adenosine methylation and reveal insights into the epigenetic regulation of the mammalian transcriptome

Rapid identification of BCR/ABL1-like acute lymphoblastic leukaemia patients using a predictive statistical model based on quantitative real time-polymerase chain reaction: clinical, prognostic and therapeutic implications.

BCR/ABL1-like acute lymphoblastic leukaemia (ALL) is a subgroup of B-lineage acute lymphoblastic leukaemia that occurs within cases without recurrent molecular rearrangements. Gene expression profiling (GEP) can identify these cases but it is expensive and not widely available. Using GEP, we identified 10 genes specifically overexpressed by BCR/ABL1-like ALL cases and used their expression values - assessed by quantitative real time-polymerase chain reaction (Q-RT-PCR) in 26 BCR/ABL1-like and 26 non-BCR/ABL1-like cases to build a statistical "BCR/ABL1-like predictor", for the identification of BCR/ABL1-like cases. By screening 142 B-lineage ALL patients with the "BCR/ABL1-like predictor", we identified 28/142 BCR/ABL1-like patients (19·7%). Overall, BCR/ABL1-like cases were enriched in JAK/STAT mutations (P < 0·001), IKZF1 deletions (P < 0·001) and rearrangements involving cytokine receptors and tyrosine kinases (P = 0·001), thus corroborating the validity of the prediction. Clinically, the BCR/ABL1-like cases identified by the BCR/ABL1-like predictor achieved a lower rate of complete remission (P = 0·014) and a worse event-free survival (P = 0·0009) compared to non-BCR/ABL1-like ALL. Consistently, primary cells from BCR/ABL1-like cases responded in vitro to ponatinib. We propose a simple tool based on Q-RT-PCR and a statistical model that is capable of easily, quickly and reliably identifying BCR/ABL1-like ALL cases at diagnosis

Epigenetic profiles signify cell fate plasticity in unipotent spermatogonial stem and progenitor cells

Spermatogonial stem and progenitor cells (SSCs) generate adult male gametes. During in vitro expansion, these unipotent murine cells spontaneously convert to multipotent adult spermatogonial-derived stem cells (MASCs). Here we investigate this conversion process through integrative transcriptomic and epigenomic analyses. We find in SSCs that promoters essential to maintenance and differentiation of embryonic stem cells (ESCs) are enriched with histone H3-lysine4 and -lysine 27 trimethylations. These bivalent modifications are maintained at most somatic promoters after conversion, bestowing MASCs an ESC-like promoter chromatin. At enhancers, the core pluripotency circuitry is activated partially in SSCs and completely in MASCs, concomitant with loss of germ cell-specific gene expression and initiation of embryonic-like programs. Furthermore, SSCs in vitro maintain the epigenomic characteristics of germ cells in vivo. Our observations suggest that SSCs encode innate plasticity through the epigenome and that both conversion of promoter chromatin states and activation of cell type-specific enhancers are prominent features of reprogramming