SulfAtlas, the sulfatase database: state of the art and new developments.

SulfAtlas (https://sulfatlas.sb-roscoff.fr/) is a knowledge-based resource dedicated to a sequence-based classification of sulfatases. Currently four sulfatase families exist (S1-S4) and the largest family (S1, formylglycine-dependent sulfatases) is divided into subfamilies by a phylogenetic approach, each subfamily corresponding to either a single characterized specificity (or few specificities in some cases) or to unknown substrates. Sequences are linked to their biochemical and structural information according to an expert scrutiny of the available literature. Database browsing was initially made possible both through a keyword search engine and a specific sequence similarity (BLAST) server. In this article, we will briefly summarize the experimental progresses in the sulfatase field in the last 6 years. To improve and speed up the (sub)family assignment of sulfatases in (meta)genomic data, we have developed a new, freely-accessible search engine using Hidden Markov model (HMM) for each (sub)family. This new tool (SulfAtlas HMM) is also a key part of the internal pipeline used to regularly update the database. SulfAtlas resource has indeed significantly grown since its creation in 2016, from 4550 sequences to 162 430 sequences in August 2022.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes.

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.

Seasonal dynamics of marine protist communities in tidally mixed coastal waters.

Major seasonal community reorganizations and associated biomass variations are landmarks of plankton ecology. However, the processes of plankton community turnover rates have not been fully elucidated so far. Here, we analyse patterns of planktonic protist community succession in temperate latitudes, based on quantitative taxonomic data from both microscopy counts (cells >10 µm) and ribosomal DNA metabarcoding (size fraction >3 µm, 18S rRNA gene) from plankton samples collected bimonthly over 8 years (2009-2016) at the SOMLIT-Astan station (Roscoff, Western English Channel). Based on morphology, diatoms were clearly the dominating group all year round and over the study period. Metabarcoding uncovered a wider diversity spectrum and revealed the prevalence of Dinophyceae and diatoms but also of Cryptophyta, Chlorophyta, Cercozoa, Syndiniales and Ciliophora in terms of read counts and or richness. The use of morphological and molecular analyses in combination allowed improving the taxonomic resolution and to identify the sequence of the dominant species and OTUs (18S V4 rDNA-derived taxa) that drive annual plankton successions. We detected that some of these dominant OTUs were benthic as a result of the intense tidal mixing typical of the French coasts in the English Channel. Our analysis of the temporal structure of community changes point to a strong seasonality and resilience. The temporal structure of environmental variables (especially Photosynthetic Active Radiation, temperature and macronutrients) and temporal structures generated by species life cycles and or species interactions, are key drivers of the observed cyclic annual plankton turnover.

Insyght: navigating amongst abundant homologues, syntenies and gene functional annotations in bacteria, it's that symbol!

High-throughput techniques have considerably increased the potential of comparative genomics whilst simultaneously posing many new challenges. One of those challenges involves efficiently mining the large amount of data produced and exploring the landscape of both conserved and idiosyncratic genomic regions across multiple genomes. Domains of application of these analyses are diverse: identification of evolutionary events, inference of gene functions, detection of niche-specific genes or phylogenetic profiling. Insyght is a comparative genomic visualization tool that combines three complementary displays: (i) a table for thoroughly browsing amongst homologues, (ii) a comparator of orthologue functional annotations and (iii) a genomic organization view designed to improve the legibility of rearrangements and distinctive loci. The latter display combines symbolic and proportional graphical paradigms. Synchronized navigation across multiple species and interoperability between the views are core features of Insyght. A gene filter mechanism is provided that helps the user to build a biologically relevant gene set according to multiple criteria such as presence/absence of homologues and/or various annotations. We illustrate the use of Insyght with scenarios. Currently, only Bacteria and Archaea are supported. A public instance is available at http://genome.jouy.inra.fr/Insyght. The tool is freely downloadable for private data set analysis.

WiseScaffolder: an algorithm for the semi-automatic scaffolding of Next Generation Sequencing data.

BACKGROUND: The sequencing depth provided by high-throughput sequencing technologies has allowed a rise in the number of de novo sequenced genomes that could potentially be closed without further sequencing. However, genome scaffolding and closure require costly human supervision that often results in genomes being published as drafts. A number of automatic scaffolders were recently released, which improved the global quality of genomes published in the last few years. Yet, none of them reach the efficiency of manual scaffolding. RESULTS: Here, we present an innovative semi-automatic scaffolder that additionally helps with chimerae resolution and generates valuable contig maps and outputs for manual improvement of the automatic scaffolding. This software was tested on the newly sequenced marine cyanobacterium Synechococcus sp. WH8103 as well as two reference datasets used in previous studies, Rhodobacter sphaeroides and Homo sapiens chromosome 14 (http://gage.cbcb.umd.edu/). The quality of resulting scaffolds was compared to that of three other stand-alone scaffolders: SSPACE, SOPRA and SCARPA. For all three model organisms, WiseScaffolder produced better results than other scaffolders in terms of contiguity statistics (number of genome fragments, N50, LG50, etc.) and, in the case of WH8103, the reliability of the scaffolds was confirmed by whole genome alignment against a closely related reference genome. We also propose an efficient computer-assisted strategy for manual improvement of the scaffolding, using outputs generated by WiseScaffolder, as well as for genome finishing that in our hands led to the circularization of the WH8103 genome. CONCLUSION: Altogether, WiseScaffolder proved more efficient than three other scaffolders for both prokaryotic and eukaryotic genomes and is thus likely applicable to most genome projects. The scaffolding pipeline described here should be of particular interest to biologists wishing to take advantage of the high added value of complete genomes.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies.

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Joint segmentation, calling, and normalization of multiple CGH profiles.

The statistical analysis of array comparative genomic hybridization (CGH) data has now shifted to the joint assessment of copy number variations at the cohort level. Considering multiple profiles gives the opportunity to correct for systematic biases observed on single profiles, such as probe GC content or the so-called "wave effect." In this article, we extend the segmentation model developed in the univariate case to the joint analysis of multiple CGH profiles. Our contribution is multiple: we propose an integrated model to perform joint segmentation, normalization, and calling for multiple array CGH profiles. This model shows great flexibility, especially in the modeling of the wave effect that gives a likelihood framework to approaches proposed by others. We propose a new dynamic programming algorithm for break point positioning, as well as a model selection criterion based on a modified bayesian information criterion proposed in the univariate case. The performance of our method is assessed using simulated and real data sets. Our method is implemented in the R package cghseg.

Global Phylogeography of Marine Synechococcus in Coastal Areas Reveals Strong Community Shifts.

Marine Synechococcus comprise a numerically and ecologically prominent phytoplankton group, playing a major role in both carbon cycling and trophic networks in all oceanic regions except in the polar oceans. Despite their high abundance in coastal areas, our knowledge of Synechococcus communities in these environments is based on only a few local studies. Here, we use the global metagenome data set of the Ocean Sampling Day (June 21st, 2014) to get a snapshot of the taxonomic composition of coastal Synechococcus communities worldwide, by recruitment on a reference database of 141 picocyanobacterial genomes, representative of the whole Prochlorococcus, Synechococcus, and Cyanobium diversity. This allowed us to unravel drastic community shifts over small to medium scale gradients of environmental factors, in particular along European coasts. The combined analysis of the phylogeography of natural populations and the thermophysiological characterization of eight strains, representative of the four major Synechococcus lineages (clades I to IV), also brought novel insights about the differential niche partitioning of clades I and IV, which most often co-dominate the Synechococcus community in cold and temperate coastal areas. Altogether, this study reveals several important characteristics and specificities of the coastal communities of Synechococcus worldwide. IMPORTANCE Synechococcus is the second most abundant phytoplanktonic organism on Earth, and its wide genetic diversity allowed it to colonize all the oceans except for polar waters, with different clades colonizing distinct oceanic niches. In recent years, the use of global metagenomics data sets has greatly improved our knowledge of "who is where" by describing the distribution of Synechococcus clades or ecotypes in the open ocean. However, little is known about the global distribution of Synechococcus ecotypes in coastal areas, where Synechococcus is often the dominant phytoplanktonic organism. Here, we leverage the global Ocean Sampling Day metagenomics data set to describe Synechococcus community composition in coastal areas worldwide, revealing striking community shifts, in particular along the coasts of Europe. As temperature appears as an important driver of the community composition, we also characterize the thermal preferenda of 8 Synechococcus strains, bringing new insights into the adaptation to temperature of the dominant Synechococcus clades.

Comparative Thermophysiology of Marine Synechococcus CRD1 Strains Isolated From Different Thermal Niches in Iron-Depleted Areas.

Marine Synechococcus cyanobacteria are ubiquitous in the ocean, a feature likely related to their extensive genetic diversity. Amongst the major lineages, clades I and IV preferentially thrive in temperate and cold, nutrient-rich waters, whilst clades II and III prefer warm, nitrogen or phosphorus-depleted waters. The existence of such cold (I/IV) and warm (II/III) thermotypes is corroborated by physiological characterization of representative strains. A fifth clade, CRD1, was recently shown to dominate the Synechococcus community in iron-depleted areas of the world ocean and to encompass three distinct ecologically significant taxonomic units (ESTUs CRD1A-C) occupying different thermal niches, suggesting that distinct thermotypes could also occur within this clade. Here, using comparative thermophysiology of strains representative of these three CRD1 ESTUs we show that the CRD1A strain MITS9220 is a warm thermotype, the CRD1B strain BIOS-U3-1 a cold temperate thermotype, and the CRD1C strain BIOS-E4-1 a warm temperate stenotherm. Curiously, the CRD1B thermotype lacks traits and/or genomic features typical of cold thermotypes. In contrast, we found specific physiological traits of the CRD1 strains compared to their clade I, II, III, and IV counterparts, including a lower growth rate and photosystem II maximal quantum yield at most temperatures and a higher turnover rate of the D1 protein. Together, our data suggests that the CRD1 clade prioritizes adaptation to low-iron conditions over temperature adaptation, even though the occurrence of several CRD1 thermotypes likely explains why the CRD1 clade as a whole occupies most iron-limited waters.

Synergic Effects of Temperature and Irradiance on the Physiology of the Marine Synechococcus Strain WH7803.

Understanding how microorganisms adjust their metabolism to maintain their ability to cope with short-term environmental variations constitutes one of the major current challenges in microbial ecology. Here, the best physiologically characterized marine Synechococcus strain, WH7803, was exposed to modulated light/dark cycles or acclimated to continuous high-light (HL) or low-light (LL), then shifted to various stress conditions, including low (LT) or high temperature (HT), HL and ultraviolet (UV) radiations. Physiological responses were analyzed by measuring time courses of photosystem (PS) II quantum yield, PSII repair rate, pigment ratios and global changes in gene expression. Previously published membrane lipid composition were also used for correlation analyses. These data revealed that cells previously acclimated to HL are better prepared than LL-acclimated cells to sustain an additional light or UV stress, but not a LT stress. Indeed, LT seems to induce a synergic effect with the HL treatment, as previously observed with oxidative stress. While all tested shift conditions induced the downregulation of many photosynthetic genes, notably those encoding PSI, cytochrome b6/f and phycobilisomes, UV stress proved to be more deleterious for PSII than the other treatments, and full recovery of damaged PSII from UV stress seemed to involve the neo-synthesis of a fairly large number of PSII subunits and not just the reassembly of pre-existing subunits after D1 replacement. In contrast, genes involved in glycogen degradation and carotenoid biosynthesis pathways were more particularly upregulated in response to LT. Altogether, these experiments allowed us to identify responses common to all stresses and those more specific to a given stress, thus highlighting genes potentially involved in niche acclimation of a key member of marine ecosystems. Our data also revealed important specific features of the stress responses compared to model freshwater cyanobacteria.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria.

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

DomainSieve: a protein domain-based screen that led to the identification of dam-associated genes with potential link to DNA maintenance.

MOTIVATION: The Dam methyltransferase (DamMT) activity, broadly distributed in association with restriction endonucleases, as part of the restriction-modification defense systems, has evolved to become intimately associated with essential biological functions in a few organisms. In Escherichia coli, DamMT is involved in multiple aspects of DNA maintenance, replication initiation, daughter chromosome segregation, DNA mismatch repair, gene expression control, etc. The participation of DamMT in such a diverse set of functions required that other genes adapted, or emerged through evolution, in response to the DamMT-induced modification of the genomic environment. One example is SeqA, a protein that senses the methylation status of the origin of replication of the chromosome to control the timing of replication initiation. Interestingly, seqA is only present in a few DamMT-specifying proteobacteria. This observation led us to hypothesize that other genes, specifying related functions, might also be found in these organisms. To test this hypothesis, we implemented a large-scale comparative genomic screen meant to identify genes specifying DNA methylation sensing domains, probably involved in DNA maintenance functions. RESULTS: We carried out a phylogenetic analysis of DamMT, identifying two contrasting behaviors of the protein. Based on this phylogeny, we defined precisely a set of genomes, in which the protein activity is likely to be involved in DNA maintenance functions, the 'resident' dam genomes. We defined a second set of genomes, in which DamMT is not resident. We developped a new tool, 'DomainSieve', in order to screen these two sets for protein domains that are strictly associated with 'resident' dam genomes. This approach was rewarding and generated a list of genes, among which some, at least, specify activities with clear linkage to DamMT-dependent DNA methylation and DNA maintenance. AVAILABILITY: DomainSieve is implemented as a web resource and is accessible at http://stat.genopole.cnrs.fr/ds/.

Detecting local high-scoring segments: a first-stage approach for genome-wide association studies.

Genetic epidemiology aims at identifying biological mechanisms responsible for human diseases. Genome-wide association studies, made possible by recent improvements in genotyping technologies, are now promisingly investigated. In these studies, common first-stage strategies focus on marginal effects but lead to multiple-testing and are unable to capture the possibly complex interplay between genetic factors. We have adapted the use of the local score statistic, already successfully applied to analyse long molecular sequences. Via sum statistics, this method captures local and possible distant dependences between markers. Dedicated to genome-wide association studies, it is fast to compute, able to handle large datasets, circumvents the the multiple-testing problem and outlines a set of genomic regions (segments) for further analyses. Applied to simulated and real data, our approach outperforms classical Bonferroni and FDR corrections for multiple-testing. It is implemented in a software termed LHiSA for Local High-scoring Segments for Association and available at: http://stat.genopole.cnrs.fr/software/lhisa.

Mining Bacillus subtilis chromosome heterogeneities using hidden Markov models.

We present here the use of a new statistical segmentation method on the Bacillus subtilis chromosome sequence. Maximum likelihood parameter estimation of a hidden Markov model, based on the expectation-maximization algorithm, enables one to segment the DNA sequence according to its local composition. This approach is not based on sliding windows; it enables different compositional classes to be separated without prior knowledge of their content, size and localization. We compared these compositional classes, obtained from the sequence, with the annotated DNA physical map, sequence homologies and repeat regions. The first heterogeneity revealed discriminates between the two coding strands and the non-coding regions. Other main heterogeneities arise; some are related to horizontal gene transfer, some to t-enriched composition of hydrophobic protein coding strands, and others to the codon usage fitness of highly expressed genes. Concerning potential and established gene transfers, we found 9 of the 10 known prophages, plus 14 new regions of atypical composition. Some of them are surrounded by repeats, most of their genes have unknown function or possess homology to genes involved in secondary catabolism, metal and antibiotic resistance. Surprisingly, we notice that all of these detected regions are a + t-richer than the host genome, raising the question of their remote sources.

Matching the Diversity of Sulfated Biomolecules: Creation of a Classification Database for Sulfatases Reflecting Their Substrate Specificity.

Sulfatases cleave sulfate groups from various molecules and constitute a biologically and industrially important group of enzymes. However, the number of sulfatases whose substrate has been characterized is limited in comparison to the huge diversity of sulfated compounds, yielding functional annotations of sulfatases particularly prone to flaws and misinterpretations. In the context of the explosion of genomic data, a classification system allowing a better prediction of substrate specificity and for setting the limit of functional annotations is urgently needed for sulfatases. Here, after an overview on the diversity of sulfated compounds and on the known sulfatases, we propose a classification database, SulfAtlas (http://abims.sb-roscoff.fr/sulfatlas/), based on sequence homology and composed of four families of sulfatases. The formylglycine-dependent sulfatases, which constitute the largest family, are also divided by phylogenetic approach into 73 subfamilies, each subfamily corresponding to either a known specificity or to an uncharacterized substrate. SulfAtlas summarizes information about the different families of sulfatases. Within a family a web page displays the list of its subfamilies (when they exist) and the list of EC numbers. The family or subfamily page shows some descriptors and a table with all the UniProt accession numbers linked to the databases UniProt, ExplorEnz, and PDB.

Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis.

Bacteria adapt to environmental stimuli by adjusting their transcriptomes in a complex manner, the full potential of which has yet to be established for any individual bacterial species. Here, we report the transcriptomes of Bacillus subtilis exposed to a wide range of environmental and nutritional conditions that the organism might encounter in nature. We comprehensively mapped transcription units (TUs) and grouped 2935 promoters into regulons controlled by various RNA polymerase sigma factors, accounting for ~66% of the observed variance in transcriptional activity. This global classification of promoters and detailed description of TUs revealed that a large proportion of the detected antisense RNAs arose from potentially spurious transcription initiation by alternative sigma factors and from imperfect control of transcription termination.

MuGeN: simultaneous exploration of multiple genomes and computer analysis results.

MOTIVATION: The availability of increasing amounts of sequence data about completely sequenced genomes spurs the development of new methods in the fields of automated annotation, and of comparative genomics. Tools allowing the visualization of results produced by analysis methods, superimposed on possibly annotated sequence data, and enabling synchronized navigation in multiple genomes, provide new means for interactive genome exploration. This kind of visual inspection can be used as a basis to assess the quality of new analysis algorithms, or to discover genome portions to be subjected to in-depth studies. RESULTS: We propose a software package, MuGeN, built for navigating through multiple annotated genomes. It is capable of retrieving annotated sequences in several formats, stored in local files, or available in databases over the network. From these, it then generates an interactive display, or an image file, in most common formats suitable for printing, further editing or integrating in Web pages. Genome maps may be mixed with computer analysis results loaded from XML files, whose format is generic enough to be adapted to a majority of sequence oriented analysis methods. AVAILABILITY: MuGeN is available at http://www-mig.jouy.inra.fr/bdsi/MuGeN.