SIS: a program to generate draft genome sequence scaffolds for prokaryotes

Background: Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Results: We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of SIS, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that SIS has overall better performance. Conclusions: SIS is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of SIS in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Giant magnetofossils and hyperthermal events

Magnetotactic bacteria biomineralize magnetic minerals with precisely controlled size, morphology, and stoichiometry. These cosmopolitan bacteria are widely observed in aquatic environments. If preserved after burial, the inorganic remains of magnetotactic bacteria act as magnetofossils that record ancient geomagnetic field variations. They also have potential to provide paleoenvironmental information. In contrast to conventional magnetofossils, giant magnetofossils (most likely produced by eukaryotic organisms) have only been reported once before from Paleocene-Eocene Thermal Maximum (PETM; 55.8 Ma) sediments on the New Jersey coastal plain. Here, using transmission electron microscopic observations, we present evidence for abundant giant magnetofossils, including previously reported elongated prisms and spindles, and new giant bullet-shaped magnetite crystals, in the Southern Ocean near Antarctica, not only during the PETM, but also shortly before and after the PETM. Moreover, we have discovered giant bullet-shaped magnetite crystals from the equatorial Indian Ocean during the Mid-Eocene Climatic Optimum (similar to 40 Ma). Our results indicate a more widespread geographic, environmental, and temporal distribution of giant magnetofossils in the geological record with a link to "hyperthermal" events. Enhanced global weathering during hyperthermals, and expanded suboxic diagenetic environments, probably provided more bioavailable iron that enabled biomineralization of giant magnetofossils. Our micromagnetic modelling indicates the presence of magnetic multi-domain (i.e., not ideal for navigation) and single domain (i.e., ideal for navigation) structures in the giant magnetite particles depending on their size, morphology and spatial arrangement. Different giant magnetite crystal morphologies appear to have had different biological functions, including magnetotaxis and other non-navigational purposes. Our observations suggest that hyperthermals provided ideal conditions for giant magnetofossils, and that these organisms were globally distributed. Much more work is needed to understand the interplay between magnetofossil morphology, climate, nutrient availability, and environmental variability.