Gene expression databases and data mining.

The DNA microarray technology has arguably caught the attention of the worldwide life science community and is now systematically supporting major discoveries in many fields of study. The majority of the initial technical challenges of conducting experiments are being resolved, only to be replaced with new informatics hurdles, including statistical analysis, data visualization, interpretation, and storage. Two systems of databases, one containing expression data and one containing annotation data are quickly becoming essential knowledge repositories of the research community. This present paper surveys several databases, which are considered "pillars" of research and important nodes in the network. This paper focuses on a generalized workflow scheme typical for microarray experiments using two examples related to cancer research. The workflow is used to reference appropriate databases and tools for each step in the process of array experimentation. Additionally, benefits and drawbacks of current array databases are addressed, and suggestions are made for their improvement.

SwissDock, a protein-small molecule docking web service based on EADock DSS.

Most life science processes involve, at the atomic scale, recognition between two molecules. The prediction of such interactions at the molecular level, by so-called docking software, is a non-trivial task. Docking programs have a wide range of applications ranging from protein engineering to drug design. This article presents SwissDock, a web server dedicated to the docking of small molecules on target proteins. It is based on the EADock DSS engine, combined with setup scripts for curating common problems and for preparing both the target protein and the ligand input files. An efficient Ajax/HTML interface was designed and implemented so that scientists can easily submit dockings and retrieve the predicted complexes. For automated docking tasks, a programmatic SOAP interface has been set up and template programs can be downloaded in Perl, Python and PHP. The web site also provides an access to a database of manually curated complexes, based on the Ligand Protein Database. A wiki and a forum are available to the community to promote interactions between users. The SwissDock web site is available online at http://www.swissdock.ch. We believe it constitutes a step toward generalizing the use of docking tools beyond the traditional molecular modeling community.

The IC3D classification of the corneal dystrophies.

BACKGROUND: The recent availability of genetic analyses has demonstrated the shortcomings of the current phenotypic method of corneal dystrophy classification. Abnormalities in different genes can cause a single phenotype, whereas different defects in a single gene can cause different phenotypes. Some disorders termed corneal dystrophies do not appear to have a genetic basis. PURPOSE: The purpose of this study was to develop a new classification system for corneal dystrophies, integrating up-to-date information on phenotypic description, pathologic examination, and genetic analysis. METHODS: The International Committee for Classification of Corneal Dystrophies (IC3D) was created to devise a current and accurate nomenclature. RESULTS: This anatomic classification continues to organize dystrophies according to the level chiefly affected. Each dystrophy has a template summarizing genetic, clinical, and pathologic information. A category number from 1 through 4 is assigned, reflecting the level of evidence supporting the existence of a given dystrophy. The most defined dystrophies belong to category 1 (a well-defined corneal dystrophy in which a gene has been mapped and identified and specific mutations are known) and the least defined belong to category 4 (a suspected dystrophy where the clinical and genetic evidence is not yet convincing). The nomenclature may be updated over time as new information regarding the dystrophies becomes available. CONCLUSIONS: The IC3D Classification of Corneal Dystrophies is a new classification system that incorporates many aspects of the traditional definitions of corneal dystrophies with new genetic, clinical, and pathologic information. Standardized templates provide key information that includes a level of evidence for there being a corneal dystrophy. The system is user-friendly and upgradeable and can be retrieved on the website www.corneasociety.org/ic3d.

Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.

A Novel Pulse-Chase SILAC Strategy Measures Changes in Protein Decay and Synthesis Rates Induced by Perturbation of Proteostasis with an Hsp90 Inhibitor.

Standard proteomics methods allow the relative quantitation of levels of thousands of proteins in two or more samples. While such methods are invaluable for defining the variations in protein concentrations which follow the perturbation of a biological system, they do not offer information on the mechanisms underlying such changes. Expanding on previous work [1], we developed a pulse-chase (pc) variant of SILAC (stable isotope labeling by amino acids in cell culture). pcSILAC can quantitate in one experiment and for two conditions the relative levels of proteins newly synthesized in a given time as well as the relative levels of remaining preexisting proteins. We validated the method studying the drug-mediated inhibition of the Hsp90 molecular chaperone, which is known to lead to increased synthesis of stress response proteins as well as the increased decay of Hsp90 "clients". We showed that pcSILAC can give information on changes in global cellular proteostasis induced by treatment with the inhibitor, which are normally not captured by standard relative quantitation techniques. Furthermore, we have developed a mathematical model and computational framework that uses pcSILAC data to determine degradation constants kd and synthesis rates Vs for proteins in both control and drug-treated cells. The results show that Hsp90 inhibition induced a generalized slowdown of protein synthesis and an increase in protein decay. Treatment with the inhibitor also resulted in widespread protein-specific changes in relative synthesis rates, together with variations in protein decay rates. The latter were more restricted to individual proteins or protein families than the variations in synthesis. Our results establish pcSILAC as a viable workflow for the mechanistic dissection of changes in the proteome which follow perturbations. Data are available via ProteomeXchange with identifier PXD000538.

High-density electric source imaging (esi) in focal epilepsy: a prospective study of 150 operated patients

Purpose: EEG is mandatory in the diagnosis of the epilepsy syndrome. However, its potential as imaging tool is still under estimated. In the present study, we aim to determine the prerequisites of maximal benefit of electric source imaging (ESI) to localize the irritative zone in patients with focal epilepsy. Methods: One hundred fifty patients suffering from focal epilepsy and with minimum 1 year postoperative follow-up were studied prospectively and blinded to the underlying diagnosis. We evaluated the influence of two important factors on sensitivity and specificity of ESI: the number of electrodes (low resolution, LR-ESI: <30 versus high resolution, HR-ESI: 128-256 electrodes), and the use of individual MRI (i-MRI) versus template MRI (t-MRI) as the head model. Findings: ESI had a sensitivity of 85% and a specificity of 87% when HR-ESI with i-MRI was used. Using LR-ESI, sensitivity decreased to 68%, or even 57% when only t-MRI was available. The sensitivity of HR-ESI/i-MRI compared favorably with those of MRI (76%), PET (69%) and ictal/interictal SPECT (64%). Interpretation: This study on a large patient group shows excellent sensitivity and specificity of ESI if 128 EEG channels or more are used for ESI and if the results are coregistered to the patient's individual MRI. Localization precision is as high as or even higher than established brain imagery techniques. HR-ESI appears to be a valuable additional imaging tool, given that larger electrode arrays are easily and rapidly applied with modern EEG equipment and that structural MRI is nearly always available for these patients.

Remodeling of DNA replication structures by the branch point translocase FANCM

Fanconi anemia (FA) is a genetically heterogeneous chromosome instability syndrome associated with congenital abnormalities, bone marrow failure, and cancer predisposition. Eight FA proteins form a nuclear core complex, which promotes tolerance of DNA lesions in S phase, but the underlying mechanisms are still elusive. We reported recently that the FA core complex protein FANCM can translocate Holliday junctions. Here we show that FANCM promotes reversal of model replication forks via concerted displacement and annealing of the nascent and parental DNA strands. Fork reversal by FANCM also occurs when the lagging strand template is partially single-stranded and bound by RPA. The combined fork reversal and branch migration activities of FANCM lead to extensive regression of model replication forks. These observations provide evidence that FANCM can remodel replication fork structures and suggest a mechanism by which FANCM could promote DNA damage tolerance in S phase

Histology-driven data mining of lipid signatures from multiple imaging mass spectrometry analyses: application to human colorectal cancer liver metastasis biopsies.

Imaging mass spectrometry (IMS) represents an innovative tool in the cancer research pipeline, which is increasingly being used in clinical and pharmaceutical applications. The unique properties of the technique, especially the amount of data generated, make the handling of data from multiple IMS acquisitions challenging. This work presents a histology-driven IMS approach aiming to identify discriminant lipid signatures from the simultaneous mining of IMS data sets from multiple samples. The feasibility of the developed workflow is evaluated on a set of three human colorectal cancer liver metastasis (CRCLM) tissue sections. Lipid IMS on tissue sections was performed using MALDI-TOF/TOF MS in both negative and positive ionization modes after 1,5-diaminonaphthalene matrix deposition by sublimation. The combination of both positive and negative acquisition results was performed during data mining to simplify the process and interrogate a larger lipidome into a single analysis. To reduce the complexity of the IMS data sets, a sub data set was generated by randomly selecting a fixed number of spectra from a histologically defined region of interest, resulting in a 10-fold data reduction. Principal component analysis confirmed that the molecular selectivity of the regions of interest is maintained after data reduction. Partial least-squares and heat map analyses demonstrated a selective signature of the CRCLM, revealing lipids that are significantly up- and down-regulated in the tumor region. This comprehensive approach is thus of interest for defining disease signatures directly from IMS data sets by the use of combinatory data mining, opening novel routes of investigation for addressing the demands of the clinical setting.

Biochemical characterization of recombinant human telomerase

Résumé Les télomères sont les structures ADN-protéines des extrémités des chromosomes des eucaryotes. L'ADN télomérique est constitué de courtes séquences répétitives. L'intégrité des télomères est essentielle pour protéger les extrémités des chromosomes contre les systèmes de dégradations et pour les distinguer des cassures de l'ADN double brin. Parce que la machinerie de la réplication de l'ADN n'est pas capable de répliquer l'extrémité des chromosomes, les télomères raccourcissent au fur et à mesure des cycles de réplication. Dès que les télomères atteignent une longueur critique, leur structure protectrice est perdue. Cela induit un signal de dommage de l'ADN et l'arrêt du cycle cellulaire. Pour contrebalancer le raccourcissement des télomères, les cellules qui s'auto régénèrent, dont les cellules de la moelle osseuse, les lymphocytes activés et 80-90% des cellules cancéreuses, expriment la télomérase. C'est une ribonucléoprotéine qui a la capacité de synthétiser des séquences télomériques par transcription inverse d'une courte séquence contenue dans sa propre sous-unité ARN avec laquelle elle est associée. La télomérase humaine est une enzyme processive au niveau de l'addition des nucléotides et aussi des répétitions télomériques. La télomérase de levure et la télomérase humaine sont toutes deux dimériques et il a été montré que la télomérase humaine recombinante contient deux ARN qui coopèrent pour fonctionner ainsi que deux sous-unités catalytiques. Cependant, il n'a pas encore été montré quel est le rôle de la dimérisation dans l'activité de la télomérase. Afin d'élucider ce rôle, nous avons exprimé, reconstitué et purifié la télomérase humaine dimérique recombinante. Et pour étudier l'effet d'ARN mutants sur l'activité de la télomérase, nous avons développé une méthode pour reconstituer et enrichir en hétérodimères de télomérase. Les hétérodimères contiennent une sous-unité ARN sauvage et une sous-unité ARN mutée au niveau de la séquence de la matrice. Sur l'ARN muté nous avons introduit une étiquette aptamer ARN-S1 puis nous avons purifié la télomérase via l'etiquette Si. Nous avons montré que la dimérisation est essentielle pour l'activité de la télomérase. Nos données indiquent que chaque télomérase du dimère allonge leur substrat, l'ADN télomérique, indépendamment l'une de l'autre à chaque cycle d'élongation mais que l'addition itérative de répétitions télomériques nécessite une coopération entre les deux télomérases du dimère. Nous proposons donc un modèle dans lequel les deux télomérases du dimères se lient et allongent deux substrats télomères et que pendant l'élongation processive les deux enzymes subissent un changement de conformation de manière coordonnée, ce changement va permettre le repositionnement des substrats pour d'autres cycles d'additions de répétitions télomériques. Dyskeratosis congenita est une maladie mortelle due majoritairement au disfonctionnement de la moelle osseuse. Dans la forme autosomale de la maladie, l'ARN de la télomérase contient des mutations. En utilisant notre système de reconstitution, nous avons montré que ces ARN mutés, qui ont perdu leur activité enzymatique dans le cas d'un homodimère de mutants, sont dominant négatifs quand ils sont présents dans les hétérodimères sauvage/mutant. Cet effet trans-dominant négatif pourrait contribuer à la progression de la maladie. Abstract Telomeres are protein-DNA structures at the ends of linear eukaryotic chromosomes. The telomeric DNA consists of tandemly repeated sequences. Telomeric integrity is essential to protect chromosomal ends from nucleolytic degradation and to prevent their recognition as DNA double strand breaks. Due to the inability of the conventional DNA replication machinery to replicate terminal DNA stretches, telomeres shorten with continuous rounds of DNA replication. As soon as telomeres reach a critical length, their protective structure is lost and the deprotected telomeres will induce a DNA damage response leading to cell cycle arrest. To counteract telomere shortening, self-renewing cells, including bone marrow cells, activated lymphocytes and 80-90% of cancer cells express the cellular reverse transcriptase telomerase, which has the capacity to synthesize telomeric repeats by reverse transcription of a short template sequence encoded by its stably associated RNA subunit. Human telomerase is a processive enzyme for nucleotide as well as repeat addition. Both yeast and human telomerase are dimeric enzymes and recombinant human telomerase has been shown to contain two functionally cooperating RNAs and most probably also two protein subunits. However, it has remained unclear how dimerization may contribute to telomerase activity. To study the role of dimerization, we expressed, reconstituted and purified recombinant human telomerase. We also developed a new method to reconstitute and enrich for telomerase heterodimers containing wild-type (wt) and mutant telomerase RNA subunits. To this end we introduced an S1-RNA-aptamer tag into telomerase RNA and purified telomerase reconstituted with a mixture of untagged and tagged RNA via the S1-tag. Using this experimental system, we introduced template mutations in the tagged RNA subunit and examined the effect of mutant RNAs on wt telomerase activity in wt/mutant heterodimers. We obtained evidence that dimerization is essential for telomerase activity. Our data indicate that the two subunits elongate telomere substrates independently of each other during single rounds of elongation, but that iterative addition of telomeric repeats requires cooperation between the two subunits. We suggest a model, in which dimeric telomerases bind and elongate two telomere substrates and that the two subunits undergo coordinated conformational changes during processive elongation that enable repositioning the substrates for subsequent rounds of repeat addition. Dyskeratosis congenita is a multisystemic disease with bone marrow failure as the major cause of death. The autosomal form of this disease was found to harbor mutations in the telomerase RNA. Using our reconstitution system, we tested whether mutant dyskeratosis telomerase RNAs behaved in a dominant negative manner. We observed that dyskeratosis telomerase RNA mutants, which lacked enzymatic activity were dominant negative, when present in wt/ mutant heterodimers. The transdominant negative effect of these mutants may contribute to disease progression.

An immuno-electron microscopical analysis of transcribing multinucleosomal templates: what happens to the histones?

Immuno-electron microscopy was used to visualize the structure of reconstituted chromatin after in vitro transcription by purified T7 RNA polymerase. T7 RNA polymerase disrupts the nucleosomal structure in the transcribed region. This disruption is not influenced by the template, linear or supercoiled, and the presence or absence of nucleosomal positioning sequences in the transcribed region. In this study, we used monoclonal autoantibodies reacting with the nucleosome core particles and epitopes within several regions of the four different core histones. Some of the residues recognized by the autoantibodies are accessible on the surface of the nucleosomes and some are more internal and therefore less exposed at the surface. We show that the loss of the nucleosomal configuration during transcription is due to the loss of histone/DNA binding and that at least part of the histones are transferred to the nascent RNA chains. Consequently, after in vitro transcription by T7 RNA polymerase, the nucleosomal template does not conserve its original configuration, and no interaction of antigen/antibodies is observed anymore in the region that has been transcribed. Therefore, we conclude that in our in vitro transcription assay, nucleosomes are detached from the template, and not simply unfolded with histones remaining attached to the DNA.

CHD8 associates with human Staf and contributes to efficient U6 RNA polymerase III transcription.

Chromatin remodeling and histone modification are essential for eukaryotic transcription regulation, but little is known about chromatin-modifying activities acting on RNA polymerase III (Pol III)-transcribed genes. The human U6 small nuclear RNA promoter, located 5' of the transcription start site, consists of a core region directing basal transcription and an activating region that recruits the transcription factors Oct-1 and Staf (ZNF143). Oct-1 activates transcription in part by helping recruit core binding factors, but nothing is known about the mechanisms of transcription activation by Staf. We show that Staf activates U6 transcription from a preassembled chromatin template in vitro and associates with several proteins linked to chromatin modification, among them chromodomain-helicase-DNA binding protein 8 (CHD8). CHD8 binds to histone H3 di- and trimethylated on lysine 4. It resides on the human U6 promoter as well as the mRNA IRF3 promoter in vivo and contributes to efficient transcription from both these promoters. Thus, Pol III transcription from type 3 promoters uses some of the same factors used for chromatin remodeling at Pol II promoters.

Eeny meeny miny moe, catch a transcript by the toe, or how to enumerate eukaryotic transcripts.

In this issue of Genes & Development, Revyakin and colleagues (pp. 1691-1702) measure the relation between individual RNA polymerase II transcription events and transcription factor assembly by counting RNA transcripts retained on the template DNA using single-molecule fluorescence.

Event-related potential analyses via single-subject topographic classification.

Introduction: Responses to external stimuli are typically investigated by averaging peri-stimulus electroencephalography (EEG) epochs in order to derive event-related potentials (ERPs) across the electrode montage, under the assumption that signals that are related to the external stimulus are fixed in time across trials. We demonstrate the applicability of a single-trial model based on patterns of scalp topographies (De Lucia et al, 2007) that can be used for ERP analysis at the single-subject level. The model is able to classify new trials (or groups of trials) with minimal a priori hypotheses, using information derived from a training dataset. The features used for the classification (the topography of responses and their latency) can be neurophysiologically interpreted, because a difference in scalp topography indicates a different configuration of brain generators. An above chance classification accuracy on test datasets implicitly demonstrates the suitability of this model for EEG data. Methods: The data analyzed in this study were acquired from two separate visual evoked potential (VEP) experiments. The first entailed passive presentation of checkerboard stimuli to each of the four visual quadrants (hereafter, "Checkerboard Experiment") (Plomp et al, submitted). The second entailed active discrimination of novel versus repeated line drawings of common objects (hereafter, "Priming Experiment") (Murray et al, 2004). Four subjects per experiment were analyzed, using approx. 200 trials per experimental condition. These trials were randomly separated in training (90%) and testing (10%) datasets in 10 independent shuffles. In order to perform the ERP analysis we estimated the statistical distribution of voltage topographies by a Mixture of Gaussians (MofGs), which reduces our original dataset to a small number of representative voltage topographies. We then evaluated statistically the degree of presence of these template maps across trials and whether and when this was different across experimental conditions. Based on these differences, single-trials or sets of a few single-trials were classified as belonging to one or the other experimental condition. Classification performance was assessed using the Receiver Operating Characteristic (ROC) curve. Results: For the Checkerboard Experiment contrasts entailed left vs. right visual field presentations for upper and lower quadrants, separately. The average posterior probabilities, indicating the presence of the computed template maps in time and across trials revealed significant differences starting at ~60-70 ms post-stimulus. The average ROC curve area across all four subjects was 0.80 and 0.85 for upper and lower quadrants, respectively and was in all cases significantly higher than chance (unpaired t-test, p<0.0001). In the Priming Experiment, we contrasted initial versus repeated presentations of visual object stimuli. Their posterior probabilities revealed significant differences, which started at 250ms post-stimulus onset. The classification accuracy rates with single-trial test data were at chance level. We therefore considered sub-averages based on five single trials. We found that for three out of four subjects' classification rates were significantly above chance level (unpaired t-test, p<0.0001). Conclusions: The main advantage of the present approach is that it is based on topographic features that are readily interpretable along neurophysiologic lines. As these maps were previously normalized by the overall strength of the field potential on the scalp, a change in their presence across trials and between conditions forcibly reflects a change in the underlying generator configurations. The temporal periods of statistical difference between conditions were estimated for each training dataset for ten shuffles of the data. Across the ten shuffles and in both experiments, we observed a high level of consistency in the temporal periods over which the two conditions differed. With this method we are able to analyze ERPs at the single-subject level providing a novel tool to compare normal electrophysiological responses versus single cases that cannot be considered part of any cohort of subjects. This aspect promises to have a strong impact on both basic and clinical research.

Characterization of a selective antagonist of neuropeptide Y at the Y2 receptor. Synthesis and pharmacological evaluation of a Y2 antagonist.

Neuropeptide Y (NPY) is a potent inhibitor of neurotransmitter release through the Y2 receptor subtype. Specific antagonists for the Y2 receptors have not yet been described. Based on the concept of template-assembled synthetic proteins we have used a cyclic template molecule containing two beta-turn mimetics for covalent attachment of four COOH-terminal fragments RQRYNH2 (NPY 33-36), termed T4-[NPY(33-36)]4. This structurally defined template-assembled synthetic protein has been tested for binding using SK-N-MC and LN319 cell lines that express the Y1 and Y2 receptor, respectively. T4-[NPY(33-36)]4 binds to the Y2 receptor with high affinity (IC50 = 67.2 nM) and has poor binding to the Y1 receptor. This peptidomimetic tested on LN319 cells at concentrations up to 10 microM shows no inhibitory effect on forskolin-stimulated cAMP levels (IC50 for NPY = 2.5 nM). Furthermore, we used confocal microscopy to examine the NPY-induced increase in intracellular calcium in single LN319 cells. Preincubation of the cells with T4-[NPY(33-36)]4 shifted to the right the dose-response curves for intracellular mobilization of calcium induced by NPY at concentrations ranging from 0.1 nM to 10 microM. Finally, we assessed the competitive antagonistic properties of T4-[NPY(33-36)]4 at presynaptic peptidergic Y2 receptors modulating noradrenaline release. the compound T4-[NPY(33-36)]4 caused a marked shift to the right of the concentration-response curve of NPY 13-36, a Y2-selective fragment, yielding a pA2 value of 8.48. Thus, to our best knowledge, T4-[NPY(33-36)]4 represents the first potent and selective Y2 antagonist.

Not one odour but two: A new model for nestmate recognition.

Recognition systems play a key role in a range of biological processes, including mate choice, immune defence and altruistic behaviour. Social insects provide an excellent model for studying recognition systems because workers need to discriminate between nestmates and non-nestmates, enabling them to direct altruistic behaviour towards closer kin and to repel potential invaders. However, the level of aggression directed towards conspecific intruders can vary enormously, even among workers within the same colony. This is usually attributed to differences in the aggression thresholds of individuals or to workers having different roles within the colony. Recent evidence from the weaver ant Oecophylla smaragdina suggests that this does not tell the whole story. Here I propose a new model for nestmate recognition based on a vector template derived from both the individual's innate odour and the shared colony odour. This model accounts for the recent findings concerning weaver ants, and also provides an alternative explanation for why the level of aggression expressed by a colony decreases as the diversity within the colony increases, even when odour is well-mixed. The model makes additional predictions that are easily tested, and represents a significant advance in our conceptualisation of recognition systems.