Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity.

Spatial patterns of coherent activity across different brain areas have been identified during the resting-state fluctuations of the brain. However, recent studies indicate that resting-state activity is not stationary, but shows complex temporal dynamics. We were interested in the spatiotemporal dynamics of the phase interactions among resting-state fMRI BOLD signals from human subjects. We found that the global phase synchrony of the BOLD signals evolves on a characteristic ultra-slow (<0.01Hz) time scale, and that its temporal variations reflect the transient formation and dissolution of multiple communities of synchronized brain regions. Synchronized communities reoccurred intermittently in time and across scanning sessions. We found that the synchronization communities relate to previously defined functional networks known to be engaged in sensory-motor or cognitive function, called resting-state networks (RSNs), including the default mode network, the somato-motor network, the visual network, the auditory network, the cognitive control networks, the self-referential network, and combinations of these and other RSNs. We studied the mechanism originating the observed spatiotemporal synchronization dynamics by using a network model of phase oscillators connected through the brain's anatomical connectivity estimated using diffusion imaging human data. The model consistently approximates the temporal and spatial synchronization patterns of the empirical data, and reveals that multiple clusters that transiently synchronize and desynchronize emerge from the complex topology of anatomical connections, provided that oscillators are heterogeneous.

EGASP: the human ENCODE Genome Annotation Assessment Project.

BACKGROUND: We present the results of EGASP, a community experiment to assess the state-of-the-art in genome annotation within the ENCODE regions, which span 1% of the human genome sequence. The experiment had two major goals: the assessment of the accuracy of computational methods to predict protein coding genes; and the overall assessment of the completeness of the current human genome annotations as represented in the ENCODE regions. For the computational prediction assessment, eighteen groups contributed gene predictions. We evaluated these submissions against each other based on a 'reference set' of annotations generated as part of the GENCODE project. These annotations were not available to the prediction groups prior to the submission deadline, so that their predictions were blind and an external advisory committee could perform a fair assessment. RESULTS: The best methods had at least one gene transcript correctly predicted for close to 70% of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into account alternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotide level, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programs relying on mRNA and protein sequences were the most accurate in reproducing the manually curated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could be verified. CONCLUSION: This is the first such experiment in human DNA, and we have followed the standards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe the results presented here contribute to the value of ongoing large-scale annotation projects and should guide further experimental methods when being scaled up to the entire human genome sequence.

p21 as a transcriptional co-repressor of S-phase and mitotic control genes.

It has been previously described that p21 functions not only as a CDK inhibitor but also as a transcriptional co-repressor in some systems. To investigate the roles of p21 in transcriptional control, we studied the gene expression changes in two human cell systems. Using a human leukemia cell line (K562) with inducible p21 expression and human primary keratinocytes with adenoviral-mediated p21 expression, we carried out microarray-based gene expression profiling. We found that p21 rapidly and strongly repressed the mRNA levels of a number of genes involved in cell cycle and mitosis. One of the most strongly down-regulated genes was CCNE2 (cyclin E2 gene). Mutational analysis in K562 cells showed that the N-terminal region of p21 is required for repression of gene expression of CCNE2 and other genes. Chromatin immunoprecipitation assays indicated that p21 was bound to human CCNE2 and other p21-repressed genes gene in the vicinity of the transcription start site. Moreover, p21 repressed human CCNE2 promoter-luciferase constructs in K562 cells. Bioinformatic analysis revealed that the CDE motif is present in most of the promoters of the p21-regulated genes. Altogether, the results suggest that p21 exerts a repressive effect on a relevant number of genes controlling S phase and mitosis. Thus, p21 activity as inhibitor of cell cycle progression would be mediated not only by the inhibition of CDKs but also by the transcriptional down-regulation of key genes.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

ANIBAL, stable isotope-based quantitative proteomics by aniline and benzoic acid labeling of amino and carboxylic groups

Identification and relative quantification of hundreds to thousands of proteins within complex biological samples have become realistic with the emergence of stable isotope labeling in combination with high throughput mass spectrometry. However, all current chemical approaches target a single amino acid functionality (most often lysine or cysteine) despite the fact that addressing two or more amino acid side chains would drastically increase quantifiable information as shown by in silico analysis in this study. Although the combination of existing approaches, e.g. ICAT with isotope-coded protein labeling, is analytically feasible, it implies high costs, and the combined application of two different chemistries (kits) may not be straightforward. Therefore, we describe here the development and validation of a new stable isotope-based quantitative proteomics approach, termed aniline benzoic acid labeling (ANIBAL), using a twin chemistry approach targeting two frequent amino acid functionalities, the carboxylic and amino groups. Two simple and inexpensive reagents, aniline and benzoic acid, in their (12)C and (13)C form with convenient mass peak spacing (6 Da) and without chromatographic discrimination or modification in fragmentation behavior, are used to modify carboxylic and amino groups at the protein level, resulting in an identical peptide bond-linked benzoyl modification for both reactions. The ANIBAL chemistry is simple and straightforward and is the first method that uses a (13)C-reagent for a general stable isotope labeling approach of carboxylic groups. In silico as well as in vitro analyses clearly revealed the increase in available quantifiable information using such a twin approach. ANIBAL was validated by means of model peptides and proteins with regard to the quality of the chemistry as well as the ionization behavior of the derivatized peptides. A milk fraction was used for dynamic range assessment of protein quantification, and a bacterial lysate was used for the evaluation of relative protein quantification in a complex sample in two different biological states

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Estimating dormant and active hematopoietic stem cell kinetics through extensive modeling of bromodeoxyuridine label-retaining cell dynamics.

Bone marrow hematopoietic stem cells (HSCs) are responsible for both lifelong daily maintenance of all blood cells and for repair after cell loss. Until recently the cellular mechanisms by which HSCs accomplish these two very different tasks remained an open question. Biological evidence has now been found for the existence of two related mouse HSC populations. First, a dormant HSC (d-HSC) population which harbors the highest self-renewal potential of all blood cells but is only induced into active self-renewal in response to hematopoietic stress. And second, an active HSC (a-HSC) subset that by and large produces the progenitors and mature cells required for maintenance of day-to-day hematopoiesis. Here we present computational analyses further supporting the d-HSC concept through extensive modeling of experimental DNA label-retaining cell (LRC) data. Our conclusion that the presence of a slowly dividing subpopulation of HSCs is the most likely explanation (amongst the various possible causes including stochastic cellular variation) of the observed long term Bromodeoxyuridine (BrdU) retention, is confirmed by the deterministic and stochastic models presented here. Moreover, modeling both HSC BrdU uptake and dilution in three stages and careful treatment of the BrdU detection sensitivity permitted improved estimates of HSC turnover rates. This analysis predicts that d-HSCs cycle about once every 149-193 days and a-HSCs about once every 28-36 days. We further predict that, using LRC assays, a 75%-92.5% purification of d-HSCs can be achieved after 59-130 days of chase. Interestingly, the d-HSC proportion is now estimated to be around 30-45% of total HSCs - more than twice that of our previous estimate.

Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease.

ObjectiveCandidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.Research Design and MethodsBy integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).Results273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.ConclusionsUsing a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.

Fourmidable: a database for ant genomics.

BACKGROUND: Fourmidable is an infrastructure to curate and share the emerging genetic, molecular, and functional genomic data and protocols for ants. DESCRIPTION: The Fourmidable assembly pipeline groups nucleotide sequences into clusters before independently assembling each cluster. Subsequently, assembled sequences are annotated via Interproscan and BLAST against general and insect-specific databases. Gene-specific information can be retrieved using gene identifiers, searching for similar sequences or browsing through inferred Gene Ontology annotations. The database will readily scale as ultra-high throughput sequence data and sequences from additional species become available. CONCLUSION: Fourmidable currently houses EST data from two ant species and microarray gene expression data for one of these. Fourmidable is publicly available at http://fourmidable.unil.ch.

Loss of mitochondrial functions associated with azole resistance in Candida glabrata results in enhanced virulence in mice.

Mitochondrial dysfunction is one of the possible mechanisms by which azole resistance can occur in Candida glabrata. Cells with mitochondrial DNA deficiency (so-called "petite mutants") upregulate ATP binding cassette (ABC) transporter genes and thus display increased resistance to azoles. Isolation of such C. glabrata mutants from patients receiving antifungal therapy or prophylaxis has been rarely reported. In this study, we characterized two sequential and related C. glabrata isolates recovered from the same patient undergoing azole therapy. The first isolate (BPY40) was azole susceptible (fluconazole MIC, 4 μg/ml), and the second (BPY41) was azole resistant (fluconazole MIC, >256 μg/ml). BPY41 exhibited mitochondrial dysfunction and upregulation of the ABC transporter genes C. glabrata CDR1 (CgCDR1), CgCDR2, and CgSNQ2. We next assessed whether mitochondrial dysfunction conferred a selective advantage during host infection by testing the virulence of BPY40 and BPY41 in mice. Surprisingly, even with in vitro growth deficiency compared to BPY40, BPY41 was more virulent (as judged by mortality and fungal tissue burden) than BPY40 in both systemic and vaginal murine infection models. The increased virulence of the petite mutant correlated with a drastic gain of fitness in mice compared to that of its parental isolate. To understand this unexpected feature, genome-wide changes in gene expression driven by the petite mutation were analyzed by use of microarrays during in vitro growth. Enrichment of specific biological processes (oxido-reductive metabolism and the stress response) was observed in BPY41, all of which was consistent with mitochondrial dysfunction. Finally, some genes involved in cell wall remodelling were upregulated in BPY41 compared to BPY40, which may partially explain the enhanced virulence of BPY41. In conclusion, this study shows for the first time that mitochondrial dysfunction selected in vivo under azole therapy, even if strongly affecting in vitro growth characteristics, can confer a selective advantage under host conditions, allowing the C. glabrata mutant to be more virulent than wild-type isolates.

Proteomic analysis of membrane rafts of melanoma cells identifies protein patterns characteristic of the tumor progression stage.

The molecular mechanisms controlling the progression of melanoma from a localized tumor to an invasive and metastatic disease are poorly understood. In the attempt to start defining a functional protein profile of melanoma progression, we have analyzed by LC-MS/MS the proteins associated with detergent resistant membranes (DRMs), which are enriched in cholesterol/sphingolipids-containing membrane rafts, of melanoma cell lines derived from tumors at different stages of progression. Since membrane rafts are involved in several biological processes, including signal transduction and protein trafficking, we hypothesized that the association of proteins with rafts can be regulated during melanoma development and affect protein function and disease progression. We have identified a total of 177 proteins in the DRMs of the cell lines examined. Among these, we have found groups of proteins preferentially associated with DRMs of either less malignant radial growth phase/vertical growth phase (VGP) cells, or aggressive VGP and metastatic cells suggesting that melanoma cells with different degrees of malignancy have different DRM profiles. Moreover, some proteins were found in DRMs of only some cell lines despite being expressed at similar levels in all the cell lines examined, suggesting the existence of mechanisms controlling their association with DRMs. We expect that understanding the mechanisms regulating DRM targeting and the activity of the proteins differentially associated with DRMs in relation to cell malignancy will help identify new molecular determinants of melanoma progression.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Automated scoring of AFLPs using RawGeno v 2.0, a free R CRAN library.

Amplified Fragment Length Polymorphisms (AFLPs) are a cheap and efficient protocol for generating large sets of genetic markers. This technique has become increasingly used during the last decade in various fields of biology, including population genomics, phylogeography, and genome mapping. Here, we present RawGeno, an R library dedicated to the automated scoring of AFLPs (i.e., the coding of electropherogram signals into ready-to-use datasets). Our program includes a complete suite of tools for binning, editing, visualizing, and exporting results obtained from AFLP experiments. RawGeno can either be used with command lines and program analysis routines or through a user-friendly graphical user interface. We describe the whole RawGeno pipeline along with recommendations for (a) setting the analysis of electropherograms in combination with PeakScanner, a program freely distributed by Applied Biosystems; (b) performing quality checks; (c) defining bins and proceeding to scoring; (d) filtering nonoptimal bins; and (e) exporting results in different formats.

Meta-analysis of heat- and chemically upregulated chaperone genes in plant and human cells.

Molecular chaperones are central to cellular protein homeostasis. In mammals, protein misfolding diseases and aging cause inflammation and progressive tissue loss, in correlation with the accumulation of toxic protein aggregates and the defective expression of chaperone genes. Bacteria and non-diseased, non-aged eukaryotic cells effectively respond to heat shock by inducing the accumulation of heat-shock proteins (HSPs), many of which molecular chaperones involved in protein homeostasis, in reducing stress damages and promoting cellular recovery and thermotolerance. We performed a meta-analysis of published microarray data and compared expression profiles of HSP genes from mammalian and plant cells in response to heat or isothermal treatments with drugs. The differences and overlaps between HSP and chaperone genes were analyzed, and expression patterns were clustered and organized in a network. HSPs and chaperones only partly overlapped. Heat-shock induced a subset of chaperones primarily targeted to the cytoplasm and organelles but not to the endoplasmic reticulum, which organized into a network with a central core of Hsp90s, Hsp70s, and sHSPs. Heat was best mimicked by isothermal treatments with Hsp90 inhibitors, whereas less toxic drugs, some of which non-steroidal anti-inflammatory drugs, weakly expressed different subsets of Hsp chaperones. This type of analysis may uncover new HSP-inducing drugs to improve protein homeostasis in misfolding and aging diseases.

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.