eVOC: a controlled vocabulary for unifying gene expression data.

Expression data contribute significantly to the biological value of the sequenced human genome, providing extensive information about gene structure and the pattern of gene expression. ESTs, together with SAGE libraries and microarray experiment information, provide a broad and rich view of the transcriptome. However, it is difficult to perform large-scale expression mining of the data generated by these diverse experimental approaches. Not only is the data stored in disparate locations, but there is frequent ambiguity in the meaning of terms used to describe the source of the material used in the experiment. Untangling semantic differences between the data provided by different resources is therefore largely reliant on the domain knowledge of a human expert. We present here eVOC, a system which associates labelled target cDNAs for microarray experiments, or cDNA libraries and their associated transcripts with controlled terms in a set of hierarchical vocabularies. eVOC consists of four orthogonal controlled vocabularies suitable for describing the domains of human gene expression data including Anatomical System, Cell Type, Pathology and Developmental Stage. We have curated and annotated 7016 cDNA libraries represented in dbEST, as well as 104 SAGE libraries,with expression information,and provide this as an integrated, public resource that allows the linking of transcripts and libraries with expression terms. Both the vocabularies and the vocabulary-annotated libraries can be retrieved from http://www.sanbi.ac.za/evoc/. Several groups are involved in developing this resource with the aim of unifying transcript expression information.

Comprehensive spatiotemporal transcriptomic analyses of the ganglionic eminences demonstrate the uniqueness of its caudal subdivision.

The elucidation of mechanisms underlying telencephalic neural development has been limited by the lack of knowledge regarding the molecular and cellular aspects of the ganglionic eminence (GE), an embryonic structure that supplies the brain with diverse sets of GABAergic neurons. Here, we report a comprehensive transcriptomic analysis of this structure including its medial (MGE), lateral (LGE) and caudal (CGE) subdivisions and its temporal dynamics in 12.5 to 16 day-old rat embryos. Surprisingly, comparison across subdivisions showed that CGE gene expression was the most unique providing unbiased genetic evidence for its differentiation from MGE and LGE. The molecular signature of the CGE comprised a large set of genes, including Rwdd3, Cyp26b1, Nr2f2, Egr3, Cpta1, Slit3, and Hod, of which several encode cell signaling and migration molecules such as WNT5A, DOCK9, VSNL1 and PRG1. Temporal analysis of the MGE revealed differential expression of unique sets of cell specification and migration genes, with early expression of Hes1, Lhx2, Ctgf and Mdk, and late enrichment of Olfm3, SerpinE2 and Wdr44. These GE profiles reveal new candidate regulators of spatiotemporally governed GABAergic neuronogenesis.

Genomic Study of RNA Polymerase II and III SNAP(c)-Bound Promoters Reveals a Gene Transcribed by Both Enzymes and a Broad Use of Common Activators.

SNAP(c) is one of a few basal transcription factors used by both RNA polymerase (pol) II and pol III. To define the set of active SNAP(c)-dependent promoters in human cells, we have localized genome-wide four SNAP(c) subunits, GTF2B (TFIIB), BRF2, pol II, and pol III. Among some seventy loci occupied by SNAP(c) and other factors, including pol II snRNA genes, pol III genes with type 3 promoters, and a few un-annotated loci, most are primarily occupied by either pol II and GTF2B, or pol III and BRF2. A notable exception is the RPPH1 gene, which is occupied by significant amounts of both polymerases. We show that the large majority of SNAP(c)-dependent promoters recruit POU2F1 and/or ZNF143 on their enhancer region, and a subset also recruits GABP, a factor newly implicated in SNAP(c)-dependent transcription. These activators associate with pol II and III promoters in G1 slightly before the polymerase, and ZNF143 is required for efficient transcription initiation complex assembly. The results characterize a set of genes with unique properties and establish that polymerase specificity is not absolute in vivo.

Radiation therapy and concurrent plus adjuvant temsirolimus (CCI-779) versus chemoirradiation with temozolomide in newly diagnosed glioblastoma without methylation of the MGMT gene promoter.

Background: Preclinical data indicate activity of mammalian target of rapamycin inhibitors and synergistic activity together with radiotherapy in glioblastoma. The aim of this trial is to assess the therapeutic activity of temsirolimus (CCI-779), an intravenous mTOR inhibitor, in patients with newly diagnosed glioblastoma with unmethylated O6 methlyguanine-DNA-methlytransferase (MGMT)promoter. Methods: Patients (n=257) with newly diagnosed glioblastoma after open surgical biopsy or resection fulfilling basic eligibility criteria underwent a central MGMT promoter analysis using quantitative methylation specific PCR. Patients with glioblastoma harboring an unmethylated MGMT promoter (n=111) were randomized 1:1 between radiotherapy (60 Gy; 5 times 2 Gy per week) plus concomitant and six cycles of maintenance temozolomide or radiotherapy plus weekly temsirolimus at 25 mg flat dose to be continued until progression or undue toxicity. Primary endpoint was overall survival at 12 months (OS12). Sample size of the investigational treatment arm required 54 patients to assess adequacy of temsirolimus activity set at 80%. More than 38 patients alive at 12 months in the per protocol population was considered a positive signal. A control arm of 54 patients treated with the standard of care was implemented to evaluate the assumptions on OS12. Results: Between December 2009 and October 2012, 111 pts in 14 centers were randomized and treated. Median age was 55 and 58 years in the temsirolimus and standard arm, respectively. Most patients (95.5%) had a WHO performance status of 0 or 1. Both therapies were properly administered with a median of 13 cycles of maintenance temsirolimus. In the per protocolpopulation, exactly 38 patients treated with temsirolimus (out of 54 eligible) reached OS12. In the intention to treat population OS12 was 72.2% [95% CI (58.2, 82.2)] in the temozolomide arm and 69.6% [95% CI (55.8, 79.9) in the temsirolimus arm [HR=1.16 95% CI (0.77, 1.76), p=0.47]. Conclusions: The therapeutic activity of temsirolimus in patients with newly diagnosed glioblastoma with an unmethylated MGMT promoter is too low.

Molecular phylogenetics of soricid shrews (Mammalia) based on mitochondrial cytochrome b gene sequences: with special reference to the Soricinae

Molecular phylogeny of soricid shrews (Soricidae, Eulipotyphla, Mammalia) based on 1140 bp mitochondrial cytochrome b gene (cytb) sequences was inferred by the maximum likelihood (ML) method. All 13 genera of extant Soricinae and two genera of Crocidurinae were included in the analyses. Anourosorex was phylogenetically distant from the main groupings within Soricinae and Crocidurinae in the ML tree. Thus, it could not be determined to which subfamily Anourosorex should be assigned: Soricinae, Crocidurinae or a new subfamily. Soricinae (excluding Anourosorex) should be divided into four tribes: Neomyini, Notiosoricini, Soricini and Blarinini. However, monophyly of Blarinini was not robust in the present data set. Also, branching orders among tribes of Soricinae and those among genera of Neomyini could not be determined because of insufficient phylogenetic information of the cytb sequences. For water shrews of Neomyini (Chimarrogale, Nectogale and Neomys), monophyly of Neomys and the Chimarrogale-Nectogale group could not be verified, which implies the possibility of multiple origins for the semi-aquatic mode of living among taxa within Neomyini. Episoriculus may contain several separate genera. Blarinella was included in Blarinini not Soricini, based on the cytb sequences, but the confidence level was rather low; hence more phylogenetic information is needed to determine its phylogenetic position. Furthermore, some specific problems of taxonomy of soricid shrews were clarified, for example phylogeny of local populations of Notiosorex crawfordi, Chimarrogale himalayica and Crocidura attenuata.

Plastid genome characterisation in Brassica and Brassicaceae using a new set of nine SSRs.

We report a new set of nine primer pairs specifically developed for amplification of Brassica plastid SSR markers. The wide utility of these markers is demonstrated for haplotype identification and detection of polymorphism in B. napus, B. nigra, B. oleracea, B. rapa and in related genera Arabidopsis, Camelina, Raphanus and Sinapis. Eleven gene regions (ndhB-rps7 spacer, rbcL-accD spacer, rpl16 intron, rps16 intron, atpB-rbcL spacer, trnE-trnT spacer, trnL intron, trnL-trnF spacer, trnM-atpE spacer, trnR-rpoC2 spacer, ycf3-psaA spacer) were sequenced from a range of Brassica and related genera for SSR detection and primer design. Other sequences were obtained from GenBank/EMBL. Eight out of nine selected SSR loci showed polymorphism when amplified using the new primers and a combined analysis detected variation within and between Brassica species, with the number of alleles detected per locus ranging from 5 (loci MF-6, MF-1) to 11 (locus MF-7). The combined SSR data were used in a neighbour-joining analysis (SMM, D (DM) distances) to group the samples based on the presence and absence of alleles. The analysis was generally able to separate plastid types into taxon-specific groups. Multi-allelic haplotypes were plotted onto the neighbour joining tree. A total number of 28 haplotypes were detected and these differentiated 22 of the 41 accessions screened from all other accessions. None of these haplotypes was shared by more than one species and some were not characteristic of their predicted type. We interpret our results with respect to taxon differentiation, hybridisation and introgression patterns relating to the 'Triangle of U'.

MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Very low prevalence of meticillin-resistant Staphylococcus aureus carrying the mecC gene in western Switzerland.

We report the first case of meticillin-resistant Staphylococcus aureus (MRSA) with the mecC gene in a patient in western Switzerland. After this first identification, a polymerase chain reaction protocol was established to investigate the occurrence of this new mecC gene in the population of this region. Enrichment broths were investigated from 1062 patients screened for MRSA, meticillin-susceptible Staphylococcus aureus isolates from clinical specimens from 475 patients, and 80 MRSA isolates (from 2005 to 2011) showing discrepancies between genotypic and phenotypic meticillin resistance. None was positive for mecC, suggesting that it is rare in the patient population of this region.

Detection of plant-modulated alterations in antifungal gene expression in Pseudomonas fluorescens CHA0 on roots by flow cytometry.

The biocontrol activity of the root-colonizing Pseudomonas fluorescens strain CHA0 is largely determined by the production of antifungal metabolites, especially 2,4-diacetylphloroglucinol. The expression of these metabolites depends on abiotic and biotic environmental factors, in particular, elements present in the rhizosphere. In this study, we have developed a new method for the in situ analysis of antifungal gene expression using flow cytometry combined with green fluorescent protein (GFP)-based reporter fusions to the phlA and prnA genes essential for the production of the antifungal compounds 2,4-diacetylphloroglucinol and pyrrolnitrin, respectively, in strain CHA0. Expression of phlA-gfp and prnA-gfp in CHA0 cells harvested from the rhizosphere of a set of plant species as well as from the roots of healthy, leaf pathogen-attacked, and physically stressed plants were analyzed using a FACSCalibur. After subtraction of background fluorescence emitted by plant-derived particles and CHA0 cells not carrying the gfp reporters, the average gene expression per bacterial cell could be calculated. Levels of phlA and prnA expression varied significantly in the rhizospheres of different plant species. Physical stress and leaf pathogen infection lowered phlA expression levels in the rhizosphere of cucumber. Our results demonstrate that the newly developed approach is suitable to monitor differences in levels of antifungal gene expression in response to various plant-derived factors. An advantage of the method is that it allows quantification of bacterial gene expression in rhizosphere populations at a single-cell level. To our best knowledge, this is the first study using flow cytometry for the in situ analysis of biocontrol gene expression in a plant-beneficial bacterium in the rhizosphere.

In vivo gene transfer into the ocular ciliary muscle mediated by ultrasound and microbubbles.

This study aimed to assess application of ultrasound (US) combined with microbubbles (MB) to transfect the ciliary muscle of rat eyes. Reporter DNA plasmids encoding for Gaussia luciferase, β-galactosidase or the green fluorescent protein (GFP), alone or mixed with 50% Artison MB, were injected into the ciliary muscle, with or without US exposure (US set at 1 MHz, 2 W/cm(2), 50% duty cycle for 2 min). Luciferase activity was measured in ocular fluids at 7 and 30 days after sonoporation. At 1 week, the US+MB treatment showed a significant increase in luminescence compared with control eyes, injected with plasmid only, with or without MB (×2.6), and, reporter proteins were localized in the ciliary muscle by histochemical analysis. At 1 month, a significant decrease in luciferase activity was observed in all groups. A rise in lens and ciliary muscle temperature was measured during the procedure but did not result in any observable or microscopic damages at 1 and 8 days. The feasibility to transfer gene into the ciliary muscle by US and MB suggests that sonoporation may allow intraocular production of proteins for the treatment of inflammatory, angiogenic and/or degenerative retinal diseases.

The Gene Ontology in 2010: extensions and refinements.

The Gene Ontology (GO) Consortium (http://www.geneontology.org) (GOC) continues to develop, maintain and use a set of structured, controlled vocabularies for the annotation of genes, gene products and sequences. The GO ontologies are expanding both in content and in structure. Several new relationship types have been introduced and used, along with existing relationships, to create links between and within the GO domains. These improve the representation of biology, facilitate querying, and allow GO developers to systematically check for and correct inconsistencies within the GO. Gene product annotation using GO continues to increase both in the number of total annotations and in species coverage. GO tools, such as OBO-Edit, an ontology-editing tool, and AmiGO, the GOC ontology browser, have seen major improvements in functionality, speed and ease of use.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Relationship between estrogen receptor alpha location and gene induction reveals the importance of downstream sites and cofactors.

BACKGROUND: To understand cancer-related modifications to transcriptional programs requires detailed knowledge about the activation of signal-transduction pathways and gene expression programs. To investigate the mechanisms of target gene regulation by human estrogen receptor alpha (hERalpha), we combine extensive location and expression datasets with genomic sequence analysis. In particular, we study the influence of patterns of DNA occupancy by hERalpha on expression phenotypes. RESULTS: We find that strong ChIP-chip sites co-localize with strong hERalpha consensus sites and detect nucleotide bias near hERalpha sites. The localization of ChIP-chip sites relative to annotated genes shows that weak sites are enriched near transcription start sites, while stronger sites show no positional bias. Assessing the relationship between binding configurations and expression phenotypes, we find binding sites downstream of the transcription start site (TSS) to be equally good or better predictors of hERalpha-mediated expression as upstream sites. The study of FOX and SP1 cofactor sites near hERalpha ChIP sites shows that induced genes frequently have FOX or SP1 sites. Finally we integrate these multiple datasets to define a high confidence set of primary hERalpha target genes. CONCLUSION: Our results support the model of long-range interactions of hERalpha with the promoter-bound cofactor SP1 residing at the promoter of hERalpha target genes. FOX motifs co-occur with hERalpha motifs along responsive genes. Importantly we show that the spatial arrangement of sites near the start sites and within the full transcript is important in determining response to estrogen signaling.

Gene-disease network analysis reveals functional modules in mendelian, complex and environmental diseases

AbstractBACKGROUND: Scientists have been trying to understand the molecular mechanisms of diseases to design preventive and therapeutic strategies for a long time. For some diseases, it has become evident that it is not enough to obtain a catalogue of the disease-related genes but to uncover how disruptions of molecular networks in the cell give rise to disease phenotypes. Moreover, with the unprecedented wealth of information available, even obtaining such catalogue is extremely difficult.PRINCIPAL FINDINGS: We developed a comprehensive gene-disease association database by integrating associations from several sources that cover different biomedical aspects of diseases. In particular, we focus on the current knowledge of human genetic diseases including mendelian, complex and environmental diseases. To assess the concept of modularity of human diseases, we performed a systematic study of the emergent properties of human gene-disease networks by means of network topology and functional annotation analysis. The results indicate a highly shared genetic origin of human diseases and show that for most diseases, including mendelian, complex and environmental diseases, functional modules exist. Moreover, a core set of biological pathways is found to be associated with most human diseases. We obtained similar results when studying clusters of diseases, suggesting that related diseases might arise due to dysfunction of common biological processes in the cell.CONCLUSIONS: For the first time, we include mendelian, complex and environmental diseases in an integrated gene-disease association database and show that the concept of modularity applies for all of them. We furthermore provide a functional analysis of disease-related modules providing important new biological insights, which might not be discovered when considering each of the gene-disease association repositories independently. Hence, we present a suitable framework for the study of how genetic and environmental factors, such as drugs, contribute to diseases.AVAILABILITY: The gene-disease networks used in this study and part of the analysis are available at http://ibi.imim.es/DisGeNET/DisGeNETweb.html#Download

Gene expression profiling in human pathogenic dermatophytes in vitro and in vivo

Objectives: Dermatophytes are highly specialized fungi which are the most common agents of superficial mycoses in humans and animals. The particular ability of these microorganisms to invade and multiply within keratinized host structures is presumably linked to their secreted keratinolytic activity, which is therefore a major putative virulence attribute of these fungi. The overall adaptation and transcriptional response of dermatophytes during protein degradation and/or infection is largely unknown. Methods: A Trichophyton rubrum cDNA microarray was developed and used for the transcriptional analysis of T. rubrum and Arthroderma benhamiae cells during growth on protein substrates. Moreover, the gene expression profile in A. benhamiae cells was monitored during infection of guinea pigs. Results: T. rubrum and A. benhamiae cells activate a large set of genes encoding secreted endo- and exoproteases during growth on soy and keratin. In addition, other specifically induced factors with potential implication in protein utilization were identified, e.g. multiple transporters, metabolic enzymes, transcription factors and hypothetical proteins with unknown function. Notably however, the protease gene expression profile in the fungal cells during infection was significantly different from the pattern elicited during in vitro growth on keratin. Conclusions: Our results suggest specific functions of individual proteases during infection, which may not be restricted to the degradation of keratin. This first, broad in vivo transcriptional profiling approach in dermatophytes gives new molecular insights into pathogenicity associated adaptation mechanisms that make these microorganisms the most successful causitive agents of superficial mycoses.