The Ocean Gene Atlas v2.0: online exploration of the biogeography and phylogeny of plankton genes

Testing hypothesis about the biogeography of genes using large data resources such as Tara Oceans marine metagenomes and metatranscriptomes requires significant hardware resources and programming skills. The new release of the ‘Ocean Gene Atlas’ (OGA2) is a freely available intuitive online service to mine large and complex marine environmental genomic databases. OGA2 datasets available have been extended and now include, from the Tara Oceans portfolio: (i) eukaryotic Metagenome-Assembled-Genomes (MAGs) and Single-cell Assembled Genomes (SAGs) (10.2E+6 coding genes), (ii) version 2 of Ocean Microbial Reference Gene Catalogue (46.8E+6 non-redundant genes), (iii) 924 MetaGenomic Transcriptomes (7E+6 unigenes), (iv) 530 MAGs from an Arctic MAG catalogue (1E+6 genes) and (v) 1888 Bacterial and Archaeal Genomes (4.5E+6 genes), and an additional dataset from the Malaspina 2010 global circumnavigation: (vi) 317 Malaspina Deep Metagenome Assembled Genomes (0.9E+6 genes). Novel analyses enabled by OGA2 include phylogenetic tree inference to visualize user queries within their context of sequence homologues from both the marine environmental dataset and the RefSeq database. An Application Programming Interface (API) now allows users to query OGA2 using command-line tools, hence providing local workflow integration. Finally, gene abundance can be interactively filtered directly on map displays using any of the available environmental variables. Ocean Gene Atlas v2.0 is freely-available at: https://tara-oceans.mio.osupytheas.fr/ocean-gene-atlas/.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

The Ocean barcode atlas: A web service to explore the biodiversity and biogeography of marine organisms

The Ocean Barcode Atlas (OBA) is a user friendly web service designed for biologists who wish to explore the biodiversity and biogeography of marine organisms locked in otherwise difficult to mine planetary scale DNA metabarcode data sets. Using just a web browser, a comprehensive picture of the diversity of a taxon or a barcode sequence is visualized graphically on world maps and interactive charts. Interactive results panels allow dynamic threshold adjustments and the display of diversity results in their environmental context measured at the time of sampling (temperature, oxygen, latitude, etc). Ecological analyses such as alpha and beta-diversity plots are produced via publication quality vector graphics representations. Currently, the Ocean Barcode Altas is deployed online with the (i) Tara Oceans eukaryotic 18S-V9 rDNA metabarcodes; (ii) Tara Oceans 16S/18S rRNA _mi Tags; and (iii) 16S-V4 V5 metabarcodes collected during the Malaspina-2010 expedition. Additional prokaryotic or eukaryotic plankton barcode data sets will be added upon availability, given they provide the required complement of barcodes (including raw reads to compute barcode abundance) associated with their contextual environmental variables. Ocean Barcode Atlas is a freely-available web service at: http://oba.mio.osupytheas.fr/ocean-atlas/.

Diversity and evolution of bacterial bioluminescence genes in the global ocean

Although bioluminescent bacteria are the most abundant and widely distributed of all light-emitting organisms, the biological role and evolutionary history of bacterial luminescence are still shrouded in mystery. Bioluminescence has so far been observed in the genomes of three families of Gammaproteobacteria in the form of canonical lux operons that adopt the CDAB(F)E(G) gene order. LuxA and luxB encode the two subunits of bacterial luciferase responsible for light-emission. Our deep exploration of public marine environmental databases considerably expands this view by providing a catalog of new lux homolog sequences, including 401 previously unknown luciferase-related genes. It also reveals a broader diversity of the lux operon organization, which we observed in previously undescribed configurations such as CEDA, CAED and AxxCE. This expanded operon diversity provides clues for deciphering lux operon evolution and propagation within the bacterial domain. Leveraging quantitative tracking of marine bacterial genes afforded by planetary scale metagenomic sampling, our study also reveals that the novel lux genes and operons described herein are more abundant in the global ocean than the canonical CDAB(F)E(G) operon.

Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms.

Video Abstract

•

A catalog of 47 million genes was generated from 370 globally distributed metagenomes

•

Meta-omics data integration disentangled the mechanisms of changes in transcript pools

•

Transcript pool changes of metabolic marker genes show distinct mechanistic patterns

•

Community turnover as a response to ocean warming may be strongest in polar regions

The Ocean Gene Atlas: exploring the biogeography of plankton genes online

The Ocean Gene Atlas is a web service to explore the biogeography of genes from marine planktonic organisms. It allows users to query protein or nucleotide sequences against global ocean reference gene catalogs. With just one click, the abundance and location of target sequences are visualized on world maps as well as their taxonomic distribution. Interactive results panels allow for adjusting cutoffs for alignment quality and displaying the abundances of genes in the context of environmental features (temperature, nutrients, etc.) measured at the time of sampling. The ease of use enables non-bioinformaticians to explore quantitative and contextualized information on genes of interest in the global ocean ecosystem. Currently the Ocean Gene Atlas is deployed with (i) the Ocean Microbial Reference Gene Catalog (OM-RGC) comprising 40 million non-redundant mostly prokaryotic gene sequences associated with both Tara Oceans and Global Ocean Sampling (GOS) gene abundances and (ii) the Marine Atlas of Tara Ocean Unigenes (MATOU) composed of >116 million eukaryote unigenes. Additional datasets will be added upon availability of further marine environmental datasets that provide the required complement of sequence assemblies, raw reads and contextual environmental parameters. Ocean Gene Atlas is a freely-available web service at: http://tara-oceans.mio.osupytheas.fr/ocean-gene-atlas/.

Microbial abundance, activity and population genomic profiling with mOTUs2

Metagenomic sequencing has greatly improved our ability to profile the composition of environmental and host-associated microbial communities. However, the dependency of most methods on reference genomes, which are currently unavailable for a substantial fraction of microbial species, introduces estimation biases. We present an updated and functionally extended tool based on universal (i.e., reference-independent), phylogenetic marker gene (MG)-based operational taxonomic units (mOTUs) enabling the profiling of >7700 microbial species. As more than 30% of them could not previously be quantified at this taxonomic resolution, relative abundance estimates based on mOTUs are more accurate compared to other methods. As a new feature, we show that mOTUs, which are based on essential housekeeping genes, are demonstrably well-suited for quantification of basal transcriptional activity of community members. Furthermore, single nucleotide variation profiles estimated using mOTUs reflect those from whole genomes, which allows for comparing microbial strain populations (e.g., across different human body sites).

Taxon Richness of "Megaviridae" Exceeds those of Bacteria and Archaea in the Ocean

Since the discovery of the giant mimivirus, evolutionarily related viruses have been isolated or identified from various environments. Phylogenetic analyses of this group of viruses, tentatively referred to as the family “Megaviridae”, suggest that it has an ancient origin that may predate the emergence of major eukaryotic lineages. Environmental genomics has since revealed that Megaviridae represents one of the most abundant and diverse groups of viruses in the ocean. In the present study, we compared the taxon richness and phylogenetic diversity of Megaviridae, Bacteria, and Archaea using DNA-dependent RNA polymerase as a common marker gene. By leveraging existing microbial metagenomic data, we found higher richness and phylogenetic diversity in this single viral family than in the two prokaryotic domains. We also obtained results showing that the evolutionary rate alone cannot account for the observed high diversity of Megaviridae lineages. These results suggest that the Megaviridae family has a deep co-evolutionary history with diverse marine protists since the early “Big-Bang” radiation of the eukaryotic tree of life.

Plankton networks driving carbon export in the oligotrophic ocean

The biological carbon pump is the process by which CO₂ is transformed to organic carbon via photosynthesis, exported through sinking particles, and finally sequestered in the deep ocean. While the intensity of the pump correlates with plankton community composition, the underlying ecosystem structure driving the process remains largely uncharacterised. Here we use environmental and metagenomic data gathered during the Tara Oceans expedition to improve our understanding of carbon export in the oligotrophic ocean. We show that specific plankton communities, from the surface and deep chlorophyll maximum, correlate with carbon export at 150 m and highlight unexpected taxa such as Radiolaria, alveolate parasites, as well as Synechococcus and their phages, as lineages most strongly associated with carbon export in the subtropical, nutrient-depleted, oligotrophic ocean. Additionally, we show that the relative abundance of just a few bacterial and viral genes can predict most of the variability in carbon export in these regions.

Deep sequencing of amplified Prasinovirus and host green algal genes from an Indian Ocean transect reveals interacting trophic dependencies and new genotypes

High-throughput sequencing of Prasinovirus DNA polymerase and host green algal (Mamiellophyceae) ribosomal RNA genes was used to analyse the diversity and distribution of these taxa over a ∼10 000 km latitudinal section of the Indian Ocean. New viral and host groups were identified among the different trophic conditions observed, and highlighted that although unknown prasinoviruses are diverse, the cosmopolitan algal genera Bathycoccus, Micromonas and Ostreococcus represent a large proportion of the host diversity. While Prasinovirus communities were correlated to both the geography and the environment, host communities were not, perhaps because the genetic marker used lacked sufficient resolution. Nevertheless, analysis of single environmental variables showed that eutrophic conditions strongly influence the distributions of both hosts and viruses. Moreover, these communities were not correlated, in their composition or specific richness. These observations could result from antagonistic dynamics, such as that illustrated in a prey-predator model, and/or because hosts might be under a complex set of selective pressures. Both of these reasons must be considered to interpret environmental surveys of viruses and hosts, because covariation does not always imply interaction.

Ocean plankton. Environmental characteristics of Agulhas rings affect interocean plankton transport

Agulhas rings provide the principal route for ocean waters to circulate from the Indo-Pacific to the Atlantic basin. Their influence on global ocean circulation is well known, but their role in plankton transport is largely unexplored. We show that, although the coarse taxonomic structure of plankton communities is continuous across the Agulhas choke point, South Atlantic plankton diversity is altered compared with Indian Ocean source populations. Modeling and in situ sampling of a young Agulhas ring indicate that strong vertical mixing drives complex nitrogen cycling, shaping community metabolism and biogeochemical signatures as the ring and associated plankton transit westward. The peculiar local environment inside Agulhas rings may provide a selective mechanism contributing to the limited dispersal of Indian Ocean plankton populations into the Atlantic.

Ocean plankton. Structure and function of the global ocean microbiome

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Ocean plankton. Patterns and ecological drivers of ocean viral communities

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

QDD version 3.1: a user-friendly computer program for microsatellite selection and primer design revisited: experimental validation of variables determining genotyping success rate

Microsatellite marker development has been greatly simplified by the use of high-throughput sequencing followed by in silico microsatellite detection and primer design. However, the selection of markers designed by the existing pipelines depends either on arbitrary criteria, or older studies on PCR success. Based on wet laboratory experiments, we have identified the following factors that are most likely to influence genotyping success rate: alignment score between the primers and the amplicon; the distance between primers and microsatellites; the length of the PCR product; target region complexity and the number of reads underlying the sequence. The QDD pipeline has been modified to include these most pertinent factors in the output to help the selection of markers. Furthermore, new features are also included in the present version: (i) not only raw sequencing reads are accepted as input, but also contigs, allowing the analysis of assembled high-coverage data; (ii) input data can be both in fasta and fastq format to facilitate the use of Illumina and IonTorrent reads; (iii) A comparison to known transposable elements allows their detection; (iv) A contamination check can be carried out by BLASTing potential markers against the nucleotide (nt) database of NCBI; (v) QDD3 is now also available imbedded into a virtual machine making installation easier and operating system independent. It can be used both on command-line version as well as integrated into a Galaxy server, providing a user-friendly interface, as well as the possibility to utilize a large variety of NGS tools.

Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities

Sequencing of 16S rDNA polymerase chain reaction (PCR) amplicons is the most common approach for investigating environmental prokaryotic diversity, despite the known biases introduced during PCR. Here we show that 16S rDNA fragments derived from Illumina-sequenced environmental metagenomes (mi tags) are a powerful alternative to 16S rDNA amplicons for investigating the taxonomic diversity and structure of prokaryotic communities. As part of the Tara Oceans global expedition, marine plankton was sampled in three locations, resulting in 29 subsamples for which metagenomes were produced by shotgun Illumina sequencing (ca. 700 Gb). For comparative analyses, a subset of samples was also selected for Roche-454 sequencing using both shotgun (m454 tags; 13 metagenomes, ca. 2.4 Gb) and 16S rDNA amplicon (454 tags; ca. 0.075 Gb) approaches. Our results indicate that by overcoming PCR biases related to amplification and primer mismatch, mi tags may provide more realistic estimates of community richness and evenness than amplicon 454 tags. In addition, mi tags can capture expected beta diversity patterns. Using mi tags is now economically feasible given the dramatic reduction in high-throughput sequencing costs, having the advantage of retrieving simultaneously both taxonomic (Bacteria, Archaea and Eukarya) and functional information from the same microbial community.

Exploring nucleo-cytoplasmic large DNA viruses in Tara Oceans microbial metagenomes

Nucleo-cytoplasmic large DNA viruses (NCLDVs) constitute a group of eukaryotic viruses that can have crucial ecological roles in the sea by accelerating the turnover of their unicellular hosts or by causing diseases in animals. To better characterize the diversity, abundance and biogeography of marine NCLDVs, we analyzed 17 metagenomes derived from microbial samples (0.2–1.6 μm size range) collected during the Tara Oceans Expedition. The sample set includes ecosystems under-represented in previous studies, such as the Arabian Sea oxygen minimum zone (OMZ) and Indian Ocean lagoons. By combining computationally derived relative abundance and direct prokaryote cell counts, the abundance of NCLDVs was found to be in the order of 10⁴–10⁵ genomes ml⁻¹ for the samples from the photic zone and 10²–10³ genomes ml⁻¹ for the OMZ. The Megaviridae and Phycodnaviridae dominated the NCLDV populations in the metagenomes, although most of the reads classified in these families showed large divergence from known viral genomes. Our taxon co-occurrence analysis revealed a potential association between viruses of the Megaviridae family and eukaryotes related to oomycetes. In support of this predicted association, we identified six cases of lateral gene transfer between Megaviridae and oomycetes. Our results suggest that marine NCLDVs probably outnumber eukaryotic organisms in the photic layer (per given water mass) and that metagenomic sequence analyses promise to shed new light on the biodiversity of marine viruses and their interactions with potential hosts.

mRNA deep sequencing reveals 75 new genes and a complex transcriptional landscape in Mimivirus

Mimivirus, a virus infecting Acanthamoeba, is the prototype of the Mimiviridae, the latest addition to the nucleocytoplasmic large DNA viruses. The Mimivirus genome encodes close to 1000 proteins, many of them never before encountered in a virus, such as four amino-acyl tRNA synthetases. To explore the physiology of this exceptional virus and identify the genes involved in the building of its characteristic intracytoplasmic "virion factory," we coupled electron microscopy observations with the massively parallel pyrosequencing of the polyadenylated RNA fractions of Acanthamoeba castellanii cells at various time post-infection. We generated 633,346 reads, of which 322,904 correspond to Mimivirus transcripts. This first application of deep mRNA sequencing (454 Life Sciences [Roche] FLX) to a large DNA virus allowed the precise delineation of the 5' and 3' extremities of Mimivirus mRNAs and revealed 75 new transcripts including several noncoding RNAs. Mimivirus genes are expressed across a wide dynamic range, in a finely regulated manner broadly described by three main temporal classes: early, intermediate, and late. This RNA-seq study confirmed the AAAATTGA sequence as an early promoter element, as well as the presence of palindromes at most of the polyadenylation sites. It also revealed a new promoter element correlating with late gene expression, which is also prominent in Sputnik, the recently described Mimivirus "virophage." These results-validated genome-wide by the hybridization of total RNA extracted from infected Acanthamoeba cells on a tiling array (Agilent)--will constitute the foundation on which to build subsequent functional studies of the Mimivirus/Acanthamoeba system.

Metagenome annotation using a distributed grid of undergraduate students

The Annotathon is a novel bioinformatics teaching environment, where undergraduate students join in a community annotation effort. Besides being a rewarding educational tool, it holds the added promise of potentially useful scientific findings.

MicroArray Facility: a laboratory information management system with extended support for Nylon based technologies

Background

High throughput gene expression profiling (GEP) is becoming a routine technique in life science laboratories. With experimental designs that repeatedly span thousands of genes and hundreds of samples, relying on a dedicated database infrastructure is no longer an option.

GEP technology is a fast moving target, with new approaches constantly broadening the field diversity. This technology heterogeneity, compounded by the informatics complexity of GEP databases, means that software developments have so far focused on mainstream techniques, leaving less typical yet established techniques such as Nylon microarrays at best partially supported.

Results

MAF (MicroArray Facility) is the laboratory database system we have developed for managing the design, production and hybridization of spotted microarrays. Although it can support the widely used glass microarrays and oligo-chips, MAF was designed with the specific idiosyncrasies of Nylon based microarrays in mind. Notably single channel radioactive probes, microarray stripping and reuse, vector control hybridizations and spike-in controls are all natively supported by the software suite. MicroArray Facility is MIAME supportive and dynamically provides feedback on missing annotations to help users estimate effective MIAME compliance. Genomic data such as clone identifiers and gene symbols are also directly annotated by MAF software using standard public resources. The MAGE-ML data format is implemented for full data export. Journalized database operations (audit tracking), data anonymization, material traceability and user/project level confidentiality policies are also managed by MAF.

Conclusion

MicroArray Facility is a complete data management system for microarray producers and end-users. Particular care has been devoted to adequately model Nylon based microarrays. The MAF system, developed and implemented in both private and academic environments, has proved a robust solution for shared facilities and industry service providers alike.

Mice lacking the UBC4-testis gene have a delay in postnatal testis development but normal spermatogenesis and fertility

Activation of ubiquitination occurs during spermatogenesis and is dependent on the induction of isoforms of the UBC4 family of ubiquitin-conjugating enzymes. The UBC4-testis isoform is testis specific, is induced in round spermatids, and demonstrates biochemical functions distinct from a ubiquitously expressed isoform UBC4-1. To explore further the function of UBC4-testis, mice bearing inactivation of this gene were produced. Homozygous (-/-) mice showed normal body growth and fertility. Although testis weight and morphology were normal in testes from adult mice, examination of young mice during the first wave of spermatogenesis revealed that testes were approximately 10% smaller in weight at 40 and 45 days of age but had become normal at 65 days of age. Overall protein content, levels of ubiquitinated proteins, and ubiquitin-conjugating activity did not differ between wild-type and homozygous (-/-) mice. Spermatid number, as well as the motility of spermatozoa isolated from the epididymis, was also normal in homozygous (-/-) mice. To determine whether the germ cells lacking UBC4-testis might be more sensitive to stress, testes from wild-type and knockout mice were exposed to heat stress by implantation in the abdominal cavity. Testes from both strains of mice showed similar rates of degeneration in response to heat. The lack of an obvious phenotype did not appear to be due to induction of other UBC4 isoforms, as shown by two-dimensional gel immunoblotting. Our data indicate that UBC4-testis plays a role in early maturation of the testis and suggest that the many UBC4 isoforms have mixed redundant and specific functions.

Submission of microarray data to public repositories

The Microarray Gene Expression Data Society believe that the time is right for journals to require that microarray data be deposited in public repositories, as a condition for publication

Feature extraction and signal processing for nylon DNA microarrays

Background

High-density DNA microarrays require automatic feature extraction methodologies and softwares. These can be a potential source of non-reproducibility of gene expression measurements. Variation in feature location or in signal integration methodology may be a significant contribution to the observed variance in gene expression levels.

Results

We explore sources of variability in feature extraction from DNA microarrays on Nylon membrane with radioactive detection. We introduce a mathematical model of the signal emission and derive methods for correcting biases such as overshining, saturation or variation in probe amount. We also provide a quality metric which can be used qualitatively to flag weak or untrusted signals or quantitatively to modulate the weight of each experiment or gene in higher level analyses (clustering or discriminant analysis).

Conclusions

Our novel feature extraction methodology, based on a mathematical model of the radioactive emission, reduces variability due to saturation, neighbourhood effects and variable probe amount. Furthermore, we provide a fully automatic feature extraction software, BZScan, which implements the algorithms described in this paper.

DNA microarray data and contextual analysis of correlation graphs

BACKGROUND

DNA microarrays are used to produce large sets of expression measurements from which specific biological information is sought. Their analysis requires efficient and reliable algorithms for dimensional reduction, classification and annotation.

RESULTS

We study networks of co-expressed genes obtained from DNA microarray experiments. The mathematical concept of curvature on graphs is used to group genes or samples into clusters to which relevant gene or sample annotations are automatically assigned. Application to publicly available yeast and human lymphoma data demonstrates the reliability of the method in spite of its simplicity, especially with respect to the small number of parameters involved.

CONCLUSIONS

We provide a method for automatically determining relevant gene clusters among the many genes monitored with microarrays. The automatic annotations and the graphical interface improve the readability of the data. A C++ implementation, called Trixy, is available from http://tagc.univ-mrs.fr/bioinformatics/trixy.html.

Breast cancer revisited using DNA array-based gene expression profiling

Breast cancer is a complex genetic disease characterized by the accumulation of multiple molecular alterations. The resulting clinical heterogeneity makes current diagnostic and therapeutic strategies less than perfectly adapted to each patient. Pathological and clinical factors are insufficient to capture the complex cascade of events that drive the clinical behavior of tumors. High-throughput molecular technologies provide novel tools to tackle this complexity. In particular, DNA arrays allow the simultaneous and quantitative analysis of the mRNA expression levels of thousands of genes in a single assay. Potential applications are multiple in the cancer field and the first research results are promising; comprehensive gene expression profiles of breast tumors are providing insights into mammary oncogenesis and are revealing new tumor subgroups previously indistinguishable. Significant advances will be the identification of new diagnostic, prognostic and predictive biomarkers as well as the discovery of new potential therapeutic targets. This review presents recent applications of DNA arrays in breast cancer research and discusses some issues to address in the near future to allow the technology to reach its full potential.

Prognosis of breast cancer and gene expression profiling using DNA arrays

Breast cancer is a complex genetic disease characterized by the accumulation of multiple molecular alterations. The resulting clinical heterogeneity makes current therapeutic strategies-based on clinicopathlogical factors-less than perfectly adapted to each patient. Today, DNA arrays, by allowing the simultaneous and quantitative analysis of the mRNA expression levels of thousands of genes in a single assay, provide novel tools to tackle this complexity. Potential applications are multiple in the cancer field and the first research results are promising. Using home-made DNA arrays in an approach easily compatible with academic research-nylon support and radioactive detection-we identified a predictor set of 23 genes whose expression patterns differentiated two groups of breast cancer patients with different survival after adjuvant chemotherapy. We then validated and further extended these results in a larger, independent and homogeneous series of poor prognosis primary breast cancers treated with adjuvant anthracyclin-based chemotherapy. We confirmed the prognostic classification provided by the 23-gene set predictor. We then improved the predictor set and refined the classification by sorting the tumors into three classes with significantly different long-term survival. These results show the potential of the technology with an accessible approach for academic research teams. We also showed that nylon DNA arrays with radioactive detection are associated with excellent sensitivity, an advantage in clinical situations where the amount of available material is limited.

EMBL-Align: a new public nucleotide and amino acid multiple sequence alignment database

The submission of multiple sequence alignment data to EMBL has grown 30-fold in the past 10 years, creating a problem of archiving them. The EBI has developed a new public database of multiple sequence alignments called EMBL-Align. It has a dedicated web-based submission tool, Webin-Align. Together they represent a comprehensive data management solution for alignment data. Webin-Align accepts all the common alignment formats and can display data in CLUSTALW format as well as a new standard EMBL-Align flat file format. The alignments are stored in the EMBL-Align database and can be queried from the EBI SRS (Sequence Retrieval System) server.

AVAILABILITY

Webin-Align: http://www.ebi.ac.uk/embl/Submission/align_top.html, EMBL-Align: ftp://ftp.ebi.ac.uk/pub/databases/embl/align, http://srs.ebi.ac.uk/

Minimum information about a microarray experiment (MIAME)-toward standards for microarray data

Microarray analysis has become a widely used tool for the generation of gene expression data on a genomic scale. Although many significant results have been derived from microarray studies, one limitation has been the lack of standards for presenting and exchanging such data. Here we present a proposal, the Minimum Information About a Microarray Experiment (MIAME), that describes the minimum information required to ensure that microarray data can be easily interpreted and that results derived from its analysis can be independently verified. The ultimate goal of this work is to establish a standard for recording and reporting microarray-based gene expression data, which will in turn facilitate the establishment of databases and public repositories and enable the development of data analysis tools. With respect to MIAME, we concentrate on defining the content and structure of the necessary information rather than the technical format for capturing it.

Divergent N-terminal sequences target an inducible testis deubiquitinating enzyme to distinct subcellular structures

Ubiquitin-specific processing proteases (UBPs) presently form the largest enzyme family in the ubiquitin system, characterized by a core region containing conserved motifs surrounded by divergent sequences, most commonly at the N-terminal end. The functions of these divergent sequences remain unclear. We identified two isoforms of a novel testis-specific UBP, UBP-t1 and UBP-t2, which contain identical core regions but distinct N termini, thereby permitting dissection of the functions of these two regions. Both isoforms were germ cell specific and developmentally regulated. Immunocytochemistry revealed that UBP-t1 was induced in step 16 to 19 spermatids while UBP-t2 was expressed in step 18 to 19 spermatids. Immunoelectron microscopy showed that UBP-t1 was found in the nucleus while UBP-t2 was extranuclear and was found in residual bodies. For the first time, we show that the differential subcellular localization was due to the distinct N-terminal sequences. When transfected into COS-7 cells, the core region was expressed throughout the cell but the UBP-t1 and UBP-t2 isoforms were concentrated in the nucleus and the perinuclear region, respectively. Fusions of each N-terminal end with green fluorescent protein yielded the same subcellular localization as the native proteins, indicating that the N-terminal ends were sufficient for determining differential localization. Interestingly, UBP-t2 colocalized with anti-gamma-tubulin immunoreactivity, indicating that like several other components of the ubiquitin system, a deubiquitinating enzyme is associated with the centrosome. Regulated expression and alternative N termini can confer specificity of UBP function by restricting its temporal and spatial loci of action.

The EMBL nucleotide sequence database

The European Molecular Biology Laboratory (EMBL) Nucleotide Sequence Database (http://www.ebi.ac. uk/embl/index.html ) is maintained at the European Bioinformatics Institute (EBI) in an international collaboration with the DNA Data Bank of Japan (DDBJ) and GenBank (USA). Data is exchanged amongst the collaborative databases on a daily basis. The major contributors to the EMBL database are individual authors and genome project groups. WEBIN is the preferred web-based submission system for individual submitters, whilst automatic procedures allow incorporation of sequence data from large-scale genome sequencing centres and from the European Patent Office (EPO). Database releases are produced quarterly. Network services allow free access to the most up-to-date data collection via Internet and WWW interfaces. EBI's Sequence Retrieval System (SRS) is a network browser for databanks in molecular biology, integrating and linking the main nucleotide and protein databases plus many specialised databases. For sequence similarity searching a variety of tools (e.g., BLITZ, FASTA, BLAST) are available which allow external users to compare their own sequences against the most currently available data in the EMBL Nucleotide Sequence Database and SWISS-PROT.

The EMBL Nucleotide Sequence Database. Contributing and accessing data

The European Molecular Biology Laboratory Nucleotide Sequence Database receives sequence and sequence annotation data from genome projects, sequencing centers, individual scientists, and patent offices. Data may be most efficiently submitted to the database using the Internet based submission tool WEBIN or via previously established genome project accounts. Biologist curators will review the data and provide accession numbers within two working days. Non-confidential data are exchanged daily in an international collaboration between EMBL. DDBJ (the DNA Databank of Japan) and GenBank (USA) and may be accessed and retrieved via the Internet with the Sequence Retrieval System (SRS). Sequence database searching algorithms (e.g., Blitz, Fasta, Blast) are available for comparison of query to database sequences.

A novel rat homolog of the Saccharomyces cerevisiae ubiquitin-conjugating enzymes UBC4 and UBC5 with distinct biochemical features is induced during spermatogenesis

The Saccharomyces cerevisiae ubiquitin-conjugating enzymes (E2s) UBC4 and UBC5 are essential for degradation of short-lived and abnormal proteins. We previously identified rat cDNAs encoding two E2s with strong sequence similarity to UBC4 and UBC5. These E2 isoforms are widely expressed in rat tissues, consistent with a fundamental cellular function for these E2s. We now report a new isoform, 8A, which despite having >91% amino acid identity with the other isoforms, shows several novel features. Expression of the 8A isoform appears restricted to the testis, is absent in early life, but is induced during puberty. Hypophysectomy reduced expression of the 8A isoform. In situ hybridization studies indicated that 8A mRNA is expressed mainly in round spermatids. Immunoblot analyses showed that 8A protein is found not only in subfractions of germ cells enriched in round spermatids but also in subfractions containing residual bodies extruded from more mature elongated spermatids, indicating that the protein possesses a longer half-life than the mRNA. Unlike all previously identified mammalian and plant homologs of S. cerevisiae UBC4, which possess a basic pI, the 8A isoform is unique in possessing an acidic pI. The small differences in sequence between the 8A isoform and other rat isoforms conferred differences in biochemical function. The 8A isoform was less effective than an isoform with a basic pI or ineffective in conjugating ubiquitin to certain fractions of testis proteins. Thus, although multiple isoforms of a specific E2 may exist to ensure performance of a critical cellular function, our data demonstrate, for the first time, that multiple genes also permit highly specialized regulation of expression of specific isoforms and that subtle differences in E2 primary structure can dictate conjugation of ubiquitin to different subsets of cellular proteins.