Integrative analysis of neuroblastoma by single-cell RNA sequencing identifies the NECTIN2-TIGIT axis as a target for immunotherapy

Pediatric patients with high-risk neuroblastoma have poor survival rates and urgently need more effective treatment options with less side effects. Since novel and improved immunotherapies may fill this need, we dissect the immunoregulatory interactions in neuroblastoma by single-cell RNA-sequencing of 24 tumors (10 pre- and 14 post-chemotherapy, including 5 pairs) to identify strategies for optimizing immunotherapy efficacy. Neuroblastomas are infiltrated by natural killer (NK), T and B cells, and immunosuppressive myeloid populations. NK cells show reduced cytotoxicity and T cells have a dysfunctional profile. Interaction analysis reveals a vast immunoregulatory network and identifies NECTIN2-TIGIT as a crucial immune checkpoint. Combined blockade of TIGIT and PD-L1 significantly reduces neuroblastoma growth, with complete responses (CR) in vivo. Moreover, addition of TIGIT+PD-L1 blockade to standard relapse treatment in a chemotherapy-resistant Th-ALKF1174L/MYCN 129/SvJ syngeneic model induces CR. In conclusion, our integrative analysis provides promising targets and a rationale for immunotherapeutic combination strategies.

An ERK-dependent molecular switch antagonizes fibrosis and promotes regeneration in spiny mice (Acomys)

Although most mammals heal injured tissues and organs with scarring, spiny mice (Acomys) naturally regenerate skin and complex musculoskeletal tissues. Now, the core signaling pathways driving mammalian tissue regeneration are poorly characterized. Here, we show that, while immediate extracellular signal-regulated kinase (ERK) activation is a shared feature of scarring (Mus) and regenerating (Acomys) injuries, ERK activity is only sustained at high levels during complex tissue regeneration. Following ERK inhibition, ear punch regeneration in Acomys shifted toward fibrotic repair. Using single-cell RNA sequencing, we identified ERK-responsive cell types. Loss- and gain-of-function experiments prompted us to uncover fibroblast growth factor and ErbB signaling as upstream ERK regulators of regeneration. The ectopic activation of ERK in scar-prone injuries induced a pro-regenerative response, including cell proliferation, extracellular matrix remodeling, and hair follicle neogenesis. Our data detail an important distinction in ERK activity between regenerating and poorly regenerating adult mammals and open avenues to redirect fibrotic repair toward regenerative healing.

Lysosomal glycogen accumulation in Pompe disease results in disturbed cytoplasmic glycogen metabolism

Pompe disease is an inherited metabolic myopathy caused by deficiency of acid alpha-glucosidase (GAA), resulting in lysosomal glycogen accumulation. Residual GAA enzyme activity affects disease onset and severity, although other factors, including dysregulation of cytoplasmic glycogen metabolism, are suspected to modulate the disease course. In this study, performed in mice and patient biopsies, we found elevated protein levels of enzymes involved in glucose uptake and cytoplasmic glycogen synthesis in skeletal muscle from mice with Pompe disease, including glycogenin (GYG1), glycogen synthase (GYS1), glucose transporter 4 (GLUT4), glycogen branching enzyme 1 (GBE1), and UDP-glucose pyrophosphorylase (UGP2). Expression levels were elevated before the loss of muscle mass and function. For first time, quantitative mass spectrometry in skeletal muscle biopsies from five adult patients with Pompe disease showed increased expression of GBE1 protein relative to healthy controls at the group level. Paired analysis of individual patients who responded well to treatment with enzyme replacement therapy (ERT) showed reduction of GYS1, GYG1, and GBE1 in all patients after start of ERT compared to baseline. These results indicate that metabolic changes precede muscle wasting in Pompe disease, and imply a positive feedforward loop in Pompe disease, in which lysosomal glycogen accumulation promotes cytoplasmic glycogen synthesis and glucose uptake, resulting in aggravation of the disease phenotype.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing

PURPOSE

Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events.

METHODS

We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing.

RESULTS

We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents.

CONCLUSION

We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

LGR6 marks nephron progenitor cells

Background

Nephron progenitor cells (NPCs) undergo a stepwise process to generate all mature nephron structures. Mesenchymal to epithelial transition (MET) is considered a multistep process of NPC differentiation to ensure progressive establishment of new nephrons. However, despite this important role, to date, no marker for NPCs undergoing MET in the nephron exists.

Results

Here, we identify LGR6 as a NPC marker, expressed in very early cap mesenchyme, pre‐tubular aggregates, renal vesicles, and in segments of S‐shaped bodies, following the trajectory of MET. By using a lineage tracing approach in embryonic explants in combination with confocal imaging and single‐cell RNA sequencing, we provide evidence for the multiple fates of LGR6+ cells during embryonic nephrogenesis. Moreover, by using long‐term in vivo lineage tracing, we show that postnatal LGR6+ cells are capable of generating the multiple lineages of the nephrons.

Conclusions

Given the profound early mesenchymal expression and MET signature of LGR6⁺ cells, together with the lineage tracing of mesenchymal LGR6⁺ cells, we conclude that LGR6+ cells contribute to all nephrogenic segments by undergoing MET. LGR6+ cells can therefore be considered an early committed NPC population during embryonic and postnatal nephrogenesis with potential regenerative capability.

Lgr6 is expressed in the earliest cap mesenchyme pool, a niche where nephrogenic progenitor cells (NPCs) are found.

Lgr6 marks NPCs undergoing mesenchymal to epithelial transition, following the main process of nephron development.

Using ex vivo and vivo lineage tracing, we show that mesenchymal Lgr6 expressing cells give rise to multiple types of mesenchymal derived nephron segments, including specialized glomerular epithelium, such as podocytes.

Tumor to normal single-cell mRNA comparisons reveal a pan-neuroblastoma cancer cell

Neuroblastoma is a childhood cancer that resembles developmental stages of the neural crest. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer (n = 19,723) with normal fetal adrenal single-cell transcriptomes (n = 57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients, and clinical phenotypes. We substantiated our findings in 650 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that a pan-neuroblastoma cancer cell state exists, which may be attractive for novel immunotherapeutic and targeted avenues.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

Abstract IA27: Patient-derived organoids in pediatric cancer research

Recent advances in in vitro culture technologies, such as adult stem cell-derived organoids, have opened up new avenues for the development of novel, more physiologic human cancer models. Such preclinical models are essential for efficient translation of basic cancer research into novel treatment regimens. We succeeded in growing organoids from a range of pediatric solid tumors, including Wilms’ tumors, renal cell carcinomas, and different types of rhabdoid tumors (i.e., AT/RT, MRT). Tumor organoids retain many characteristics of parental tumor tissue. For instance, Wilms’ tumor organoids retain the cellular heterogeneity of tumors, as they are composed of an intricate network of different cell types. Moreover, we demonstrate that tumor organoids are amenable to gene editing and high-throughput drug screens. In conclusion, our pediatric cancer organoids capture disease and tissue heterogeneity and provide a platform for basic cancer research, drug screening, and personalized medicine.

Citation Format: Camilla Calandrini, Frans Schutgens, Rurika Oka, Thanasis Margaritis, Tito Candelli, Luka Mathijsen, Carola Ammerlaan, Ravian van Ineveld, Sepideh Derakhshan, Lars Custers, Philip Lijnzaad, Harry Begthel, Hinri Kerstens, Maarten Rookmaker, Marianne Verhaar, Patrick Kemmeren, Ronald de Krijger, Kathy Pritchard-Jones, Anne Rios, Marry van den Heuvel-Eibrink, Frank Holstege, Ruben van Boxtel, Hans Clevers, Jarno Drost. Patient-derived organoids in pediatric cancer research [abstract]. In: Proceedings of the AACR Special Conference on the Advances in Pediatric Cancer Research; 2019 Sep 17-20; Montreal, QC, Canada. Philadelphia (PA): AACR; Cancer Res 2020;80(14 Suppl):Abstract nr IA27.

An organoid biobank for childhood kidney cancers that captures disease and tissue heterogeneity

Kidney tumours are among the most common solid tumours in children, comprising distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first paediatric cancer organoid biobank. It contains tumour and matching normal kidney organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumours, malignant rhabdoid tumours, renal cell carcinomas, and congenital mesoblastic nephromas. Paediatric kidney tumour organoids retain key properties of native tumours, useful for revealing patient-specific drug sensitivities. Using single cell RNA-sequencing and high resolution 3D imaging, we further demonstrate that organoid cultures derived from Wilms tumours consist of multiple different cell types, including epithelial, stromal and blastemal-like cells. Our organoid biobank captures the heterogeneity of paediatric kidney tumours, providing a representative collection of well-characterised models for basic cancer research, drug-screening and personalised medicine.

Pre-clinical cell culture models capturing the heterogeneity of childhood kidney tumours are limited. Here, the authors establish and characterise an organoid biobank of tumour and matched normal organoid cultures from over 50 children with different subtypes of kidney cancer.

CHETAH: a selective, hierarchical cell type identification method for single-cell RNA sequencing

Cell type identification is essential for single-cell RNA sequencing (scRNA-seq) studies, currently transforming the life sciences. CHETAH (CHaracterization of cEll Types Aided by Hierarchical classification) is an accurate cell type identification algorithm that is rapid and selective, including the possibility of intermediate or unassigned categories. Evidence for assignment is based on a classification tree of previously available scRNA-seq reference data and includes a confidence score based on the variance in gene expression per cell type. For cell types represented in the reference data, CHETAH’s accuracy is as good as existing methods. Its specificity is superior when cells of an unknown type are encountered, such as malignant cells in tumor samples which it pinpoints as intermediate or unassigned. Although designed for tumor samples in particular, the use of unassigned and intermediate types is also valuable in other exploratory studies. This is exemplified in pancreas datasets where CHETAH highlights cell populations not well represented in the reference dataset, including cells with profiles that lie on a continuum between that of acinar and ductal cell types. Having the possibility of unassigned and intermediate cell types is pivotal for preventing misclassification and can yield important biological information for previously unexplored tissues.

The ability of transcription factors to differentially regulate gene expression is a crucial component of the mechanism underlying inversion, a frequently observed genetic interaction pattern

Genetic interactions, a phenomenon whereby combinations of mutations lead to unexpected effects, reflect how cellular processes are wired and play an important role in complex genetic diseases. Understanding the molecular basis of genetic interactions is crucial for deciphering pathway organization as well as understanding the relationship between genetic variation and disease. Several hypothetical molecular mechanisms have been linked to different genetic interaction types. However, differences in genetic interaction patterns and their underlying mechanisms have not yet been compared systematically between different functional gene classes. Here, differences in the occurrence and types of genetic interactions are compared for two classes, gene-specific transcription factors (GSTFs) and signaling genes (kinases and phosphatases). Genome-wide gene expression data for 63 single and double deletion mutants in baker's yeast reveals that the two most common genetic interaction patterns are buffering and inversion. Buffering is typically associated with redundancy and is well understood. In inversion, genes show opposite behavior in the double mutant compared to the corresponding single mutants. The underlying mechanism is poorly understood. Although both classes show buffering and inversion patterns, the prevalence of inversion is much stronger in GSTFs. To decipher potential mechanisms, a Petri Net modeling approach was employed, where genes are represented as nodes and relationships between genes as edges. This allowed over 9 million possible three and four node models to be exhaustively enumerated. The models show that a quantitative difference in interaction strength is a strict requirement for obtaining inversion. In addition, this difference is frequently accompanied with a second gene that shows buffering. Taken together, these results provide a mechanistic explanation for inversion. Furthermore, the ability of transcription factors to differentially regulate expression of their targets provides a likely explanation why inversion is more prevalent for GSTFs compared to kinases and phosphatases.

Developing a strategy for computational lab skills training through Software and Data Carpentry: Experiences from the ELIXIR Pilot action

Quality training in computational skills for life scientists is essential to allow them to deliver robust, reproducible and cutting-edge research. A pan-European bioinformatics programme, ELIXIR, has adopted a well-established and progressive programme of computational lab and data skills training from Software and Data Carpentry, aimed at increasing the number of skilled life scientists and building a sustainable training community in this field. This article describes the Pilot action, which introduced the Carpentry training model to the ELIXIR community.

Molecular mechanisms that distinguish TFIID housekeeping from regulatable SAGA promoters

An important distinction is frequently made between constitutively expressed housekeeping genes versus regulated genes. Although generally characterized by different DNA elements, chromatin architecture and cofactors, it is not known to what degree promoter classes strictly follow regulatability rules and which molecular mechanisms dictate such differences. We show that SAGA‐dominated/TATA‐box promoters are more responsive to changes in the amount of activator, even compared to TFIID/TATA‐like promoters that depend on the same activator Hsf1. Regulatability is therefore an inherent property of promoter class. Further analyses show that SAGA/TATA‐box promoters are more dynamic because TATA‐binding protein recruitment through SAGA is susceptible to removal by Mot1. In addition, the nucleosome configuration upon activator depletion shifts on SAGA/TATA‐box promoters and seems less amenable to preinitiation complex formation. The results explain the fundamental difference between housekeeping and regulatable genes, revealing an additional facet of combinatorial control: an activator can elicit a different response dependent on core promoter class.

A high-resolution gene expression atlas of epistasis between gene-specific transcription factors exposes potential mechanisms for genetic interactions

BACKGROUND

Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease.

RESULTS

To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression".

CONCLUSIONS

These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

Cell cycle population effects in perturbation studies

Growth condition perturbation or gene function disruption are commonly used strategies to study cellular systems. Although it is widely appreciated that such experiments may involve indirect effects, these frequently remain uncharacterized. Here, analysis of functionally unrelated Saccharyomyces cerevisiae deletion strains reveals a common gene expression signature. One property shared by these strains is slower growth, with increased presence of the signature in more slowly growing strains. The slow growth signature is highly similar to the environmental stress response (ESR), an expression response common to diverse environmental perturbations. Both environmental and genetic perturbations result in growth rate changes. These are accompanied by a change in the distribution of cells over different cell cycle phases. Rather than representing a direct expression response in single cells, both the slow growth signature and ESR mainly reflect a redistribution of cells over different cell cycle phases, primarily characterized by an increase in the G1 population. The findings have implications for any study of perturbation that is accompanied by growth rate changes. Strategies to counter these effects are presented and discussed.

A novel microscopy-based high-throughput screening method to identify proteins that regulate global histone modification levels

Posttranslational modifications of histones play an important role in the regulation of gene expression and chromatin structure in eukaryotes. The balance between chromatin factors depositing (writers) and removing (erasers) histone marks regulates the steady-state levels of chromatin modifications. Here we describe a novel microscopy-based screening method to identify proteins that regulate histone modification levels in a high-throughput fashion. We named our method CROSS, for Chromatin Regulation Ontology SiRNA Screening. CROSS is based on an siRNA library targeting the expression of 529 proteins involved in chromatin regulation. As a proof of principle, we used CROSS to identify chromatin factors involved in histone H3 methylation on either lysine-4 or lysine-27. Furthermore, we show that CROSS can be used to identify chromatin factors that affect growth in cancer cell lines. Taken together, CROSS is a powerful method to identify the writers and erasers of novel and known chromatin marks and facilitates the identification of drugs targeting epigenetic modifications.

Gcn4 misregulation reveals a direct role for the evolutionary conserved EKC/KEOPS in the t6A modification of tRNAs

The EKC/KEOPS complex is universally conserved in Archaea and Eukarya and has been implicated in several cellular processes, including transcription, telomere homeostasis and genomic instability. However, the molecular function of the complex has remained elusive so far. We analyzed the transcriptome of EKC/KEOPS mutants and observed a specific profile that is highly enriched in targets of the Gcn4p transcriptional activator. GCN4 expression was found to be activated at the translational level in mutants via the defective recognition of the inhibitory upstream ORFs (uORFs) present in its leader. We show that EKC/KEOPS mutants are defective for the N6-threonylcarbamoyl adenosine modification at position 37 (t⁶A₃₇) of tRNAs decoding ANN codons, which affects initiation at the inhibitory uORFs and provokes Gcn4 de-repression. Structural modeling reveals similarities between Kae1 and bacterial enzymes involved in carbamoylation reactions analogous to t⁶A₃₇ formation, supporting a direct role for the EKC in tRNA modification. These findings are further supported by strong genetic interactions of EKC mutants with a translation initiation factor and with threonine biosynthesis genes. Overall, our data provide a novel twist to understanding the primary function of the EKC/KEOPS and its impact on several essential cellular functions like transcription and telomere homeostasis.

Functional overlap and regulatory links shape genetic interactions between signaling pathways

To understand relationships between phosphorylation-based signaling pathways, we analyzed 150 deletion mutants of protein kinases and phosphatases in S. cerevisiae using DNA microarrays. Downstream changes in gene expression were treated as a phenotypic readout. Double mutants with synthetic genetic interactions were included to investigate genetic buffering relationships such as redundancy. Three types of genetic buffering relationships are identified: mixed epistasis, complete redundancy, and quantitative redundancy. In mixed epistasis, the most common buffering relationship, different gene sets respond in different epistatic ways. Mixed epistasis arises from pairs of regulators that have only partial overlap in function and that are coupled by additional regulatory links such as repression of one by the other. Such regulatory modules confer the ability to control different combinations of processes depending on condition or context. These properties likely contribute to the evolutionary maintenance of paralogs and indicate a way in which signaling pathways connect for multiprocess control.

A consensus of core protein complex compositions for Saccharomyces cerevisiae

Analyses of biological processes would benefit from accurate definitions of protein complexes. High-throughput mass spectrometry data offer the possibility of systematically defining protein complexes; however, the predicted compositions vary substantially depending on the algorithm applied. We determine consensus compositions for 409 core protein complexes from Saccharomyces cerevisiae by merging previous predictions with a new approach. Various analyses indicate that the consensus is comprehensive and of high quality. For 85 out of 259 complexes not recorded in GO, literature search revealed strong support in the form of coprecipitation. New complexes were verified by an independent interaction assay and by gene expression profiling of strains with deleted subunits, often revealing which cellular processes are affected. The consensus complexes are available in various formats, including a merge with GO, resulting in 518 protein complex compositions. The utility is further demonstrated by comparison with binary interaction data to reveal interactions between core complexes.

Genome-wide location of the coactivator mediator: Binding without activation and transient Cdk8 interaction on DNA

Mediator is a general coactivator of RNA polymerase II (Pol II) transcription. Genomic location analyses of different Mediator subunits indicate a uniformly composed core complex upstream of active genes but unexpectedly also upstream of inactive genes and on the coding regions of some highly active genes. The repressive Cdk8 submodule is associated with core Mediator at all sites but with a lower degree of occupancy, indicating transient interaction, regardless of promoter activity. This suggests gene-specific regulation of Cdk8 activity, rather than regulated Cdk8 recruitment. Mediator presence is not necessarily linked to transcription. This goes beyond Cdk8-repressed genes, indicating that Mediator can mark some regulatory regions ahead of additional signals. Overlap with intergenic Pol II location in stationary phase points to a role as a binding platform for inactive Pol II during quiescence. These results shed light on Cdk8 repression, suggest additional roles for Mediator, and query models of recruitment-coupled regulation.

Genome-wide analyses reveal RNA polymerase II located upstream of genes poised for rapid response upon S. cerevisiae stationary phase exit

The resting state of eukaryotic cells (G0) is relatively uncharacterized. We have applied DNA microarray expression profiling of S. cerevisiae to reveal multiple transitions during a complete 9-day growth cycle between stationary phase (SP) exit and entry. The findings include distinct waves of transcription after the diauxic shift (DS), identification of genes active in SP, and upregulation of over 2500 genes during the first minutes of lag phase. This provides a framework for analyzing large-scale reprogramming of gene expression. Despite global repression, the general transcription machinery is found to be present in quiescent cells but is largely inactive. Genome-wide location analysis by chromatin immunoprecipitation (ChIP on chip) reveals that RNA polymerase II is more predominantly bound at intergenic regions in SP, upstream of hundreds of genes immediately induced upon exit. In contrast to current models of activation-coupled recruitment, the results show that RNA polymerase II is located and maintained upstream of many inactive genes in quiescence.

An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas

Metastasis is the process by which cancers spread to distinct sites in the body. It is the principal cause of death in individuals suffering from cancer. For some types of cancer, early detection of metastasis at lymph nodes close to the site of the primary tumor is pivotal for appropriate treatment. Because it can be difficult to detect lymph node metastases reliably, many individuals currently receive inappropriate treatment. We show here that DNA microarray gene-expression profiling can detect lymph node metastases for primary head and neck squamous cell carcinomas that arise in the oral cavity and oropharynx. The predictor, established with an 82-tumor training set, outperforms current clinical diagnosis when independently validated. The 102 predictor genes offer unique insights into the processes underlying metastasis. The results show that the metastatic state can be deciphered from the primary tumor gene-expression pattern and that treatment can be substantially improved.

The Ensembl genome database project

The Ensembl (http://www.ensembl.org/) database project provides a bioinformatics framework to organise biology around the sequences of large genomes. It is a comprehensive source of stable automatic annotation of the human genome sequence, with confirmed gene predictions that have been integrated with external data sources, and is available as either an interactive web site or as flat files. It is also an open source software engineering project to develop a portable system able to handle very large genomes and associated requirements from sequence analysis to data storage and visualisation. The Ensembl site is one of the leading sources of human genome sequence annotation and provided much of the analysis for publication by the international human genome project of the draft genome. The Ensembl system is being installed around the world in both companies and academic sites on machines ranging from supercomputers to laptops.

RHdb: the Radiation Hybrid database

Since July 1995, the European Bioinformatics Institute (EBI) has maintained RHdb (http://www.ebi.ac.uk/RHdb), a public database for radiation hybrid data. Radiation hybrid mapping is an important technique for determining high resolution maps. RHdb is also served by CORBA servers. The EBI is an Outstation of the European Molecular Biology Laboratory (EMBL).

Accessing and distributing EMBL data using CORBA (common object request broker architecture)

BACKGROUND

The EMBL Nucleotide Sequence Database is a comprehensive database of DNA and RNA sequences and related information traditionally made available in flat-file format. Queries through tools such as SRS (Sequence Retrieval System) also return data in flat-file format. Flat files have a number of shortcomings, however, and the resources therefore currently lack a flexible environment to meet individual researchers' needs. The Object Management Group's common object request broker architecture (CORBA) is an industry standard that provides platform-independent programming interfaces and models for portable distributed object-oriented computing applications. Its independence from programming languages, computing platforms and network protocols makes it attractive for developing new applications for querying and distributing biological data.

RESULTS

A CORBA infrastructure developed by EMBL-EBI provides an efficient means of accessing and distributing EMBL data. The EMBL object model is defined such that it provides a basis for specifying interfaces in interface definition language (IDL) and thus for developing the CORBA servers. The mapping from the object model to the relational schema in the underlying Oracle database uses the facilities provided by PersistenceTM, an object/relational tool. The techniques of developing loaders and 'live object caching' with persistent objects achieve a smart live object cache where objects are created on demand. The objects are managed by an evictor pattern mechanism.

CONCLUSIONS

The CORBA interfaces to the EMBL database address some of the problems of traditional flat-file formats and provide an efficient means for accessing and distributing EMBL data. CORBA also provides a flexible environment for users to develop their applications by building clients to our CORBA servers, which can be integrated into existing systems.

RHdb: the radiation hybrid database

Since July 1995, the European Bioinformatics Institute (EBI) has maintained RHdb (http://www.ebi.ac.uk/RHdb ), a public database for radiation hybrid data. Radiation hybrid mapping is an important technique for determining high resolution maps. RHdb is also served by CORBA servers. The EBI is an Outstation of the European Molecular Biology Laboratory (EMBL).

A proposal for a standard CORBA interface for genome maps

MOTIVATION

The scientific community urgently needs to standardize the exchange of biological data. This is helped by the use of a common protocol and the definition of shared data structures. We have based our standardization work on CORBA, a technology that has become a standard in the past years and allows interoperability between distributed objects.

RESULTS

We have defined an IDL specification for genome maps and present it to the scientific community. We have implemented CORBA servers based on this IDL to distribute RHdb and HuGeMap maps. The IDL will co-evolve with the needs of the mapping community.

AVAILABILITY

The standard IDL for genome maps is available at http:// corba.ebi.ac.uk/RHdb/EUCORBA/MapIDL.htm l. The IORs to browse maps from Infobiogen and EBI are at http://www.infobiogen.fr/services/Hugemap/IOR and http://corba.ebi.ac.uk/RHdb/EUCORBA/IOR CONTACT: manu@infobiogen.fr, tome@ebi.ac.uk

The radiation hybrid database

Since July 1995, the European Bioinformatics Institute (EBI) has maintained the Radiation Hybrid database (RHdb; http://www.ebi.ac. uk/RHdb ), a public database for radiation hybrid data. Radiation hybrid mapping is an important technique for determining high resolution maps. Recently, CORBA access has been added to RHdb. The EBI is an Outstation of the European Molecular Biology Laboratory (EMBL).

A physical map of 30,000 human genes

A map of 30,181 human gene-based markers was assembled and integrated with the current genetic map by radiation hybrid mapping. The new gene map contains nearly twice as many genes as the previous release, includes most genes that encode proteins of known function, and is twofold to threefold more accurate than the previous version. A redesigned, more informative and functional World Wide Web site (www.ncbi.nlm.nih.gov/genemap) provides the mapping information and associated data and annotations. This resource constitutes an important infrastructure and tool for the study of complex genetic traits, the positional cloning of disease genes, the cross-referencing of mammalian genomes, and validated human transcribed sequences for large-scale studies of gene expression.

The Radiation Hybrid Database

Since July 1995, the European Bioinformatics Institute (EBI) has maintained RHdb (http://www.ebi.ac.uk/RHdb/RHdb.html ), a public database for radiation hybrid data. Radiation hybrid mapping is an important technique for determining high resolution maps. Recently, CORBA access has been added to Rhdb. The EBI is an Outstation of the European Molecular Biology Laboratory (EMBL).

A CORBA server for the Radiation Hybrid DataBase

Modern biology depends on a wide range of software interacting with a large number of data sources, varying both in size, complexity and structure. The range of important databases in molecular biology and genetics makes it crucial to overcome the problems which this multiplicity presents. At EMBL-EBI we have started to use CORBA technology to support interoperability between a variety of databases, as well as to facilitate the integration of tools that access these databases. Within the Radiation Hybrid DataBase project we are confronted daily with the interoperation and linking issues. In this paper we present a CORBA infrastructure implemented to access the Radiation Hybrid DataBase.

Hydrophobic patches on protein subunit interfaces: characteristics and prediction

Hydrophobic patches, defined as clusters of neighboring apolar atoms deemed accessible on a given protein surface, have been investigated on protein subunit interfaces. The data were taken from known tertiary structures of multimeric protein complexes. Amino acid composition and preference, patch size distribution, and patch contact complementarity across associating subunits were examined and compared with hydrophobic patches found on the solvent-accessible surface of the multimeric complexes. The largest or second largest patch on the accessible surface of the entire subunit was involved in multimeric interfaces in 90% of the cases. These results should prove useful for subunit design and engineering as well as for prediction of subunit interface regions.

The Radiation Hybrid Database

Since July 1995, the European Bioinformatics Institute (EBI) has maintained RHdb, a public database for radiation hybrid data. Radiation hybrid data are used in the generation of alternative genetic maps as they can include non-polymorphic markers and are also powerful enough to order unresolved genetic clusters of polymorphic STSs. The EBI is an Outstation of the European Molecular Biology Laboratory (EMBL).

A method for detecting hydrophobic patches on protein surfaces

A method for the detection of hydrophobic patches on the surfaces of protein tertiary structures is presented. It delineates explicit contiguous pieces of surface of arbitrary size and shape that consist solely of carbon and sulphur atoms using a dot representation of the solvent-accessible surface. The technique is also useful in detecting surface segments with other characteristics, such as polar patches. Its potential as a tool in the study of protein-protein interactions and substrate recognition is demonstrated by applying the method to myoglobin, Leu/IIe/Val-binding protein, lipase, lysozyme, azurin, triose phosphate isomerase, carbonic anhydrase, and phosphoglycerate kinase. Only the largest patches, having sizes exceeding random expectation, are deemed meaningful. In addition to well-known hydrophobic patches on these proteins, a number of other patches are found, and their significance is discussed. The method is simple, fast, and robust. The program text is obtainable by anonymous ftp.

Hydrophobic patches on the surfaces of protein structures

A survey of hydrophobic patches on the surface of 112 soluble, monomeric proteins is presented. The largest patch on each individual protein averages around 400 A2 but can range from 200 to 1,200 A2. These areas are not correlated to the sizes of the proteins and only weakly to their apolar surface fraction. Ala, Lys, and Pro have dominating contributions to the apolar surface for smaller patches, while those of the hydrophobic amino acids become more important as the patch size increases. The hydrophilic amino acids expose an approximately constant fraction of their apolar area independent of patch size; the hydrophobic residue types reach similar exposure only in the larger patches. Though the mobility of residues on the surface is generally higher, it decreases for hydrophilic residues with increasing patch size. Several characteristics of hydrophobic patches catalogued here should prove useful in the design and engineering of proteins.

Tertiary structure of two-electron reduced Megasphaera elsdenii flavodoxin and some implications, as determined by two-dimensional 1H-NMR and restrained molecular dynamics

The tertiary structure of the non-crystallizable two-electron-reduced Megasphaera elsdenii flavodoxin (15 kDa, 137 amino acid residues) has been determined using nuclear Overhauser enhancement restraints extracted from two-dimensional 1H-NMR spectra. A tertiary structure satisfying the experimental restraints very well (maximum NOE violation of 66 pm) was obtained with use of restrained molecular dynamics, using 509 distance restraints (including one non-NOE) on a starting structure modeled from the crystal structure of one-electron-reduced Clostridium MP flavodoxin. The protein consists of a central parallel beta-sheet surrounded on both sides by two alpha-helices. The flavin is positioned at the periphery of the molecule. The tertiary structure of the protein is highly defined with the exception of the flavin. The latter is expected to result from performing the restrained molecular dynamics simulation without water molecules and without proper charges on the flavin. The flavin, including the phosphate, the ribityl side chain and the isoalloxazine ring, is solvent accessible under the experimental conditions used and evidenced by a two-dimensional amide exchange experiment. This accessibility is expected to be important in the redox potential regulation of the semiquinone/hydroquinone couple of the protein. The amide exchange against deuterons and several typical line shapes in the two-dimensional NMR spectra are consistent with the structure generated. The structure is discussed in detail.