Robotic data acquisition with deep learning enables cell image-based prediction of transcriptomic phenotypes.

Single-cell whole-transcriptome analysis is the gold standard approach to identifying molecularly defined cell phenotypes. However, this approach cannot be used for dynamics measurements such as live-cell imaging. Here, we developed a multifunctional robot, the automated live imaging and cell picking system (ALPS) and used it to perform single-cell RNA sequencing for microscopically observed cells with multiple imaging modes. Using robotically obtained data that linked cell images and the whole transcriptome, we successfully predicted transcriptome-defined cell phenotypes in a noninvasive manner using cell image-based deep learning. This noninvasive approach opens a window to determine the live-cell whole transcriptome in real time. Moreover, this work, which is based on a data-driven approach, is a proof of concept for determining the transcriptome-defined phenotypes (i.e., not relying on specific genes) of any cell from cell images using a model trained on linked datasets.

Diverse mosquito-specific flaviviruses in the Bolivian Amazon basin.

The genus Flavivirus includes a range of mosquito-specific viruses in addition to well-known medically important arboviruses. Isolation and comprehensive genomic analyses of viruses in mosquitoes collected in Bolivia resulted in the identification of three novel flavivirus species. Psorophora flavivirus (PSFV) was isolated from Psorophora albigenu. The coding sequence of the PSFV polyprotein shares 60 % identity with that of the Aedes-associated lineage II insect-specific flavivirus (ISF), Marisma virus. Isolated PSFV replicates in both Aedes albopictus- and Aedes aegypti-derived cells, but not in mammalian Vero or BHK-21 cell lines. Two other flaviviruses, Ochlerotatus scapularis flavivirus (OSFV) and Mansonia flavivirus (MAFV), which were identified from Ochlerotatus scapularis and Mansonia titillans, respectively, group with the classical lineage I ISFs. The protein coding sequences of these viruses share only 60 and 40 % identity with the most closely related of known lineage I ISFs, including Xishuangbanna aedes flavivirus and Sabethes flavivirus, respectively. Phylogenetic analysis suggests that MAFV is clearly distinct from the groups of the current known Culicinae-associated lineage I ISFs. Interestingly, the predicted amino acid sequence of the MAFV capsid protein is approximately two times longer than that of any of the other known flaviviruses. Our results indicate that flaviviruses with distinct features can be found at the edge of the Bolivian Amazon basin at sites that are also home to dense populations of human-biting mosquitoes.

Full-length 16S rRNA gene amplicon analysis of human gut microbiota using MinION™ nanopore sequencing confers species-level resolution.

BACKGROUND: Species-level genetic characterization of complex bacterial communities has important clinical applications in both diagnosis and treatment. Amplicon sequencing of the 16S ribosomal RNA (rRNA) gene has proven to be a powerful strategy for the taxonomic classification of bacteria. This study aims to improve the method for full-length 16S rRNA gene analysis using the nanopore long-read sequencer MinION™. We compared it to the conventional short-read sequencing method in both a mock bacterial community and human fecal samples. RESULTS: We modified our existing protocol for full-length 16S rRNA gene amplicon sequencing by MinION™. A new strategy for library construction with an optimized primer set overcame PCR-associated bias and enabled taxonomic classification across a broad range of bacterial species. We compared the performance of full-length and short-read 16S rRNA gene amplicon sequencing for the characterization of human gut microbiota with a complex bacterial composition. The relative abundance of dominant bacterial genera was highly similar between full-length and short-read sequencing. At the species level, MinION™ long-read sequencing had better resolution for discriminating between members of particular taxa such as Bifidobacterium, allowing an accurate representation of the sample bacterial composition. CONCLUSIONS: Our present microbiome study, comparing the discriminatory power of full-length and short-read sequencing, clearly illustrated the analytical advantage of sequencing the full-length 16S rRNA gene.

Nucleotide Archival Format (NAF) enables efficient lossless reference-free compression of DNA sequences.

SUMMARY: DNA sequence databases use compression such as gzip to reduce the required storage space and network transmission time. We describe Nucleotide Archival Format (NAF)-a new file format for lossless reference-free compression of FASTA and FASTQ-formatted nucleotide sequences. Nucleotide Archival Format compression ratio is comparable to the best DNA compressors, while providing dramatically faster decompression. We compared our format with DNA compressors: DELIMINATE and MFCompress, and with general purpose compressors: gzip, bzip2, xz, brotli and zstd. AVAILABILITY AND IMPLEMENTATION: NAF compressor and decompressor, as well as format specification are available at https://github.com/KirillKryukov/naf. Format specification is in public domain. Compressor and decompressor are open source under the zlib/libpng license, free for nearly any use. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of a distinct lineage of aviadenovirus from crane feces.

Viruses are believed to be ubiquitous; however, the diversity of viruses is largely unknown because of the bias of previous research toward pathogenic viruses. Deep sequencing is a promising and unbiased approach to detect viruses from animal-derived materials. Although cranes are known to be infected by several viruses such as influenza A viruses, previous studies targeted limited species of viruses, and thus viruses that infect cranes have not been extensively studied. In this study, we collected crane fecal samples in the Izumi plain in Japan, which is an overwintering site for cranes, and performed metagenomic shotgun sequencing analyses. We detected aviadenovirus-like sequences in the fecal samples and tentatively named the discovered virus crane-associated adenovirus 1 (CrAdV-1). We determined that our sequence accounted for approximately three-fourths of the estimated CrAdV-1 genome size (33,245 bp). The GC content of CrAdV-1 genome is 34.1%, which is considerably lower than that of other aviadenoviruses. Phylogenetic analyses revealed that CrAdV-1 clusters with members of the genus Aviadenovirus, but is distantly related to the previously identified aviadenoviruses. The protein sequence divergence between the DNA polymerase of CrAdV-1 and those of other aviadenoviruses is 45.2-46.8%. Based on these results and the species demarcation for the family Adenoviridae, we propose that CrAdV-1 be classified as a new species in the genus Aviadenovirus. Results of this study contribute to a deeper understanding of the diversity and evolution of viruses and provide additional information on viruses that infect cranes, which might lead to protection of the endangered species of cranes.

A partial nuclear genome of the Jomons who lived 3000 years ago in Fukushima, Japan.

The Jomon period of the Japanese Archipelago, characterized by cord-marked 'jomon' potteries, has yielded abundant human skeletal remains. However, the genetic origins of the Jomon people and their relationships with modern populations have not been clarified. We determined a total of 115 million base pair nuclear genome sequences from two Jomon individuals (male and female each) from the Sanganji Shell Mound (dated 3000 years before present) with the Jomon-characteristic mitochondrial DNA haplogroup N9b, and compared these nuclear genome sequences with those of worldwide populations. We found that the Jomon population lineage is best considered to have diverged before diversification of present-day East Eurasian populations, with no evidence of gene flow events between the Jomon and other continental populations. This suggests that the Sanganji Jomon people descended from an early phase of population dispersals in East Asia. We also estimated that the modern mainland Japanese inherited <20% of Jomon peoples' genomes. Our findings, based on the first analysis of Jomon nuclear genome sequence data, firmly demonstrate that the modern mainland Japanese resulted from genetic admixture of the indigenous Jomon people and later migrants.

Functional genomic analysis of the isolated potential probiotic Lactobacillus delbrueckii subsp. indicus TY-11 and its comparison with other Lactobacillus delbrueckii strains.

Probiotics refer to living microorganisms that exert a variety of beneficial effects on human health. On the contrary, they also can cause infection, produce toxins within the body, and transfer antibiotic-resistant genes to the other microorganisms in the digestive tract necessitating a comprehensive safety assessment. This study aimed to conduct functional genomic analysis and some relevant biochemical tests to uncover the probiotic potentials of Lactobacillus delbrueckii subsp. indicus TY-11 isolated from native yogurt in Bangladesh. We also performed transmission electron microscopic (TEM) analysis, comparative genomic study as well as phylogenetic tree construction with 332 core genes from 262 genomes. The strain TY-11 was identified as Lactobacillus delbrueckii subsp. indicus, whose genome (1,916,674 bp) contained 1911 CDS, and no gene was identified for either antibiotic resistance or toxic metabolites. It carried genes for the degradation of toxic metabolites, treatment of lactose intolerance, toll-like receptor 2-dependent innate immune response, heat and cold shock, bile salts tolerance, and acidic pH tolerance. Genes were annotated for inhibiting pathogenic bacteria by inhibitory substances [bacteriocin: Helveticin-J (331 bp) and Enterolysin-A (275 bp), hydrogen peroxide, and acid]; blockage of adhesion sites; and competition for nutrients. The genes involved in its metabolic pathway were detected as suitable for digesting indigestible nutrients in the human gut. The TY-11 genome possessed an additional 37 core genes of subspecies indicus which were deficient in the core genome of the most popular subsp. bulgaricus. During the phenotypic testing, the isolate TY-11 demonstrated high antagonistic activity (inhibition zone of 21.33 ± 1.53 mm) against Escherichia coli ATCC 8739 and was not sensitive to any of the 10 tested antibiotics. This study was the first study to explore the molecular insights into probiotic roles, including antimicrobial activities and antibiotic sensitivity, of a representative strain (TY-11) of Lactobacillus delbrueckii subsp. indicus. IMPORTANCE: This study aimed to conduct functional genomic analysis to uncover the probiotic potential of Lactobacillus delbrueckii subsp. indicus TY-11 isolated from native yogurt in Bangladesh. We also performed transmission electron microscopic (TEM) analysis, comparative genomic study as well as phylogenetic tree construction with 332 core genes from 262 genomes. In our current investigation, we revealed a number of common and unique excellences of the probiotic Lactobacillus delbrueckii subsp. indicus TY-11 that are likely to be important to illustrate its intestinal residence and probiotic roles. This is the first study to explore the molecular insights into intestinal residence and probiotic roles, including antimicrobial activities and antibiotic sensitivity, of a representative strain (TY-11) of Lactobacillus delbrueckii subsp. indicus.

Nanopore Sequencing Data Analysis of 16S rRNA Genes Using the GenomeSync-GSTK System.

With the development of nanopore sequencing technology, long reads of DNA sequences can now be determined rapidly from various samples. This protocol introduces the GenomeSync-GSTK system for bacterial species identification in a given sample using nanopore sequencing data of 16S rRNA genes as an example. GenomeSync is a collection of genome sequences designed to provide easy access to genomic data of the species as demanded. GSTK (genome search toolkit) is a set of scripts for managing local homology searches using genomes obtained from the GenomeSync database. Based on this protocol, nanopore sequencing data analyses of metagenomes and amplicons could be efficiently performed. We also noted reanalysis in conjunction with future developments in nanopore sequencing technology and the accumulation of genome sequencing data.

Global Antimicrobial Resistance Gene Study of Helicobacter pylori: Comparison of Detection Tools, ARG and Efflux Pump Gene Analysis, Worldwide Epidemiological Distribution, and Information Related to the Antimicrobial-Resistant Phenotype.

We conducted a global-scale study to identify H. pylori antimicrobial-resistant genes (ARG), address their global distribution, and understand their effect on the antimicrobial resistance (AMR) phenotypes of the clinical isolates. We identified ARG using several well-known tools against extensive bacterial ARG databases, then analyzed their correlation with clinical antibiogram data from dozens of patients across countries. This revealed that combining multiple tools and databases, followed by manual selection of ARG from the annotation results, produces more conclusive results than using a single tool or database alone. After curation, the results showed that H. pylori has 42 ARG against 11 different antibiotic classes (16 genes related to single antibiotic class resistance and 26 genes related to multidrug resistance). Further analysis revealed that H. pylori naturally harbors ARG in the core genome, called the 'Set of ARG commonly found in the Core Genome of H. pylori (ARG-CORE)', while ARG-ACC-the ARG in the accessory genome-are exclusive to particular strains. In addition, we detected 29 genes of potential efflux pump-related AMR that were mostly categorized as ARG-CORE. The ARG distribution appears to be almost similar either by geographical or H. pylori populations perspective; however, some ARG had a unique distribution since they tend to be found only in a particular region or population. Finally, we demonstrated that the presence of ARG may not directly correlate with the sensitive/resistance phenotype of clinical patient isolates but may influence the minimum inhibitory concentration phenotype.

SARS-CoV-2 HaploGraph: visualization of SARS-CoV-2 haplotype spread in Japan.

Since the early phase of the coronavirus disease 2019 (COVID-19) pandemic, a number of research institutes have been sequencing and sharing high-quality severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes to trace the route of infection in Japan. To provide insight into the spread of COVID-19, we developed a web platform named SARS-CoV-2 HaploGraph to visualize the emergence timing and geographical transmission of SARS-CoV-2 haplotypes. Using data from the GISAID EpiCoV database as of June 4, 2022, we created a haplotype naming system by determining the ancestral haplotype for each epidemic wave and showed prefecture- or region-specific haplotypes in each of four waves in Japan. The SARS-CoV-2 HaploGraph allows for interactive tracking of virus evolution and of geographical prevalence of haplotypes, and aids in developing effective public health control strategies during the global pandemic. The code and the data used for this study are publicly available at: https://github.com/ktym/covid19/.

Efficient compression of SARS-CoV-2 genome data using Nucleotide Archival Format.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome data are essential for epidemiology, vaccine development, and tracking emerging variants. Millions of SARS-CoV-2 genomes have been sequenced during the pandemic. However, downloading SARS-CoV-2 genomes from databases is slow and unreliable, largely due to suboptimal choice of compression method. We evaluated the available compressors and found that Nucleotide Archival Format (NAF) would provide a drastic improvement compared with current methods. For Global Initiative on Sharing Avian Flu Data's (GISAID) pre-compressed datasets, NAF would increase efficiency 52.2 times for gzip-compressed data and 3.7 times for xz-compressed data. For DNA DataBank of Japan (DDBJ), NAF would improve throughput 40 times for gzip-compressed data. For GenBank and European Nucleotide Archive (ENA), NAF would accelerate data distribution by a factor of 29.3 times compared with uncompressed FASTA. This article provides a tutorial for installing and using NAF. Offering a NAF download option in sequence databases would provide a significant saving of time, bandwidth, and disk space and accelerate biological and medical research worldwide.

MinION, a portable long-read sequencer, enables rapid vaginal microbiota analysis in a clinical setting.

BACKGROUND: It has been suggested that the local microbiota in the reproductive organs is relevant to women's health and may also affect pregnancy outcomes. Analysis of partial 16S ribosomal RNA (rRNA) gene sequences generated by short-read sequencers has been used to identify vaginal and endometrial microbiota, but it requires a long time to obtain the results, making it unsuitable for rapid bacterial identification from a small specimen amount in a clinical context. METHODS: We developed a simple workflow using the nanopore sequencer MinION that allows high-resolution and rapid differentiation of vaginal microbiota. Vaginal samples collected from 18 participants were subjected to DNA extraction and full-length 16S rRNA gene sequencing with MinION. RESULTS: The principal coordinate analysis showed no differences in the bacterial compositions regardless of the sample collection method. The analysis of vaginal microbiota could be completed with a total analysis time of approximately four hours, allowing same-day results. Taxonomic profiling by MinION sequencing revealed relatively low diversity of the vaginal bacterial community, identifying the prevailing Lactobacillus species and several causative agents of bacterial vaginosis. CONCLUSIONS: Full-length 16S rRNA gene sequencing analysis with MinION provides a rapid means for identifying vaginal bacteria with higher resolution. Species-level profiling of human vaginal microbiota by MinION sequencing can allow the analysis of associations with conditions such as genital infections, endometritis, and threatened miscarriage.

Helicobacter pylori genomes reveal Paleolithic human migration to the east end of Asia.

A virulence bacterium, Helicobacter pylori, evolved parallel to its host human, therefore, can work as a marker for tracing the human migration. We found H. pylori strains indigenous in the southernmost islands of Japanese Archipelago, Okinawa, and defined them as hspOkinawa and hpRyukyu. Genome data of the strains revealed that hspOkinawa diverged from other East Asian strains about 20,000 years ago, and that hpRyukyu diverged about 45,000 years ago. The closest strains of hpRyukyu were found from Afghanistan, Punjab, and Nepal, which suggest this strain originated in the central Asia and traveled across the Eurasian continent during Paleolithic era. The divergence date of hpRyukyu corresponds with human fossil records in Okinawa. Although it is controversial from human DNA analyses whether descendants of the Paleolithic migrants remain in the modern Japanese population, this study reveals that the bacterium of Paleolithic origin remains in the stomachs of current Japanese.

A circulating subset of iNKT cells mediates antitumor and antiviral immunity.

Invariant natural killer T (iNKT) cells are a group of innate-like T lymphocytes that recognize lipid antigens. They are supposed to be tissue resident and important for systemic and local immune regulation. To investigate the heterogeneity of iNKT cells, we recharacterized iNKT cells in the thymus and peripheral tissues. iNKT cells in the thymus were divided into three subpopulations by the expression of the natural killer cell receptor CD244 and the chemokine receptor CXCR6 and designated as C0 (CD244-CXCR6-), C1 (CD244-CXCR6+), or C2 (CD244+CXCR6+) iNKT cells. The development and maturation of C2 iNKT cells from C0 iNKT cells strictly depended on IL-15 produced by thymic epithelial cells. C2 iNKT cells expressed high levels of IFN-γ and granzymes and exhibited more NK cell-like features, whereas C1 iNKT cells showed more T cell-like characteristics. C2 iNKT cells were influenced by the microbiome and aging and suppressed the expression of the autoimmune regulator AIRE in the thymus. In peripheral tissues, C2 iNKT cells were circulating that were distinct from conventional tissue-resident C1 iNKT cells. Functionally, C2 iNKT cells protected mice from the tumor metastasis of melanoma cells by enhancing antitumor immunity and promoted antiviral immune responses against influenza virus infection. Furthermore, we identified human CD244+CXCR6+ iNKT cells with high cytotoxic properties as a counterpart of mouse C2 iNKT cells. Thus, this study reveals a circulating subset of iNKT cells with NK cell-like properties distinct from conventional tissue-resident iNKT cells.

Rapid profiling of drug-resistant bacteria using DNA-binding dyes and a nanopore-based DNA sequencer.

Spread of drug-resistant bacteria is a serious problem worldwide. We thus designed a new sequence-based protocol that can quickly identify bacterial compositions of clinical samples and their drug-resistance profiles simultaneously. Here we utilized propidium monoazide (PMA) that prohibits DNA amplifications from dead bacteria, and subjected the original and antibiotics-treated samples to 16S rRNA metagenome sequencing. We tested our protocol on bacterial mixtures, and observed that sequencing reads derived from drug-resistant bacteria were significantly increased compared with those from drug-sensitive bacteria when samples were treated by antibiotics. Our protocol is scalable and will be useful for quickly profiling drug-resistant bacteria.

Diagnosis of pleural empyema/parapneumonic effusion by next-generation sequencing.

BACKGROUND: Although a microbiological diagnosis of pleural infection is clinically important, it is often complicated by prior antibiotic treatment and/or difficulties with culturing some bacterial species. Therefore, we aimed to identify probable causative bacteria in pleural empyema/parapneumonic effusions by combining 16S ribosomal RNA (rRNA) gene amplification and next-generation sequencing (NGS). METHODS: Pleural fluids were collected from 19 patients with infectious effusions and nine patients with non-infectious malignant effusions. We analysed DNA extracted from the pleural fluid supernatant by NGS using the Genome Search Toolkit and GenomeSync database, either directly or after PCR amplification of the 16S rRNA gene. Infectious and non-infectious effusions were distinguished by semi-quantitative PCR of the 16S rRNA gene. RESULTS: Only 8 (42%) effusions were culture-positive, however, NGS of the 16S rRNA gene amplicon identified 14 anaerobes and 7 aerobes/facultative anaerobes in all patients, including Streptococcus sp. (n = 6), Fusobacterium sp. (n = 5), Porphyromonas sp. (n = 5), and Prevotella sp. (n = 4), accounting for >10% of the total genomes. The culture and NGS results were discordant for 3 out of 8 patients, all of whom had previously been treated with antibiotics. Total (2ΔCT value in semi-quantitative PCR of the 16S rRNA gene) and specific (total bacterial load multiplied by the proportion of primary bacteria in NGS) bacterial loads could efficiently distinguish empyema/parapneumonic effusion from non-infectious effusion. CONCLUSION: Combining NGS with semi-quantitative PCR can facilitate the diagnosis of pleural empyema/parapneumonic effusion and its causal bacteria.

Small fraction of testicular cancer cases may be causatively related to CHEK2 inactivating germ-line mutations: evidence for somatic loss of the remaining CHEK2 allele in the tumor tissue.

A recent study suggested a role of CHEK2 loss-of-function germ-line pathogenic variants in the predisposition to testicular cancer (TC) (AlDubayan et al. JAMA Oncol 5:514-522, 2019). We attempted to validate this finding relying on the high population frequency of recurrent CHEK2 pathogenic variants in Slavic populations. CHEK2 pathogenic alleles (c.1100delC (p.Thr367Metfs); del5395 [del ex9-10]; IVS2 + 1G > A [c.444 + 1G > A]) were detected in 7/280 (2.5%) TC patients vs. 3/424 (0.7%) healthy men and 6/1007 (0.6%) healthy women [OR 4.0 (95% CI 1.5-11), p = 0.009 for pooled control groups]. Somatic CHEK2 loss-of-heterozygosity (LOH) was detected in 4 out of 6 tumors available for analysis; strikingly all these instances of LOH involved inactivation of the wild-type allele. The CHEK2 c.470T > C (p.Ile157Thr) variant was detected in 21/280 (7.5%) affected vs. 22/424 (5.2%) non-affected men [OR 1.5 (95% CI 0.8-2.7), p = 0.3]. Somatic CHEK2 LOH was revealed only in 6 out of 21 tumors obtained from CHEK2 c.470T > C (p.Ile157Thr) carriers, with the C-allele lost in two cases and T-allele deleted in four tumors. The results of comparison of allele frequencies in TC patients versus population controls coupled with the data on CHEK2 LOH status in tumor tissues support the association of CHEK2 pathogenic variants with TC risk.

MISHIMA--a new method for high speed multiple alignment of nucleotide sequences of bacterial genome scale data.

BACKGROUND: Large nucleotide sequence datasets are becoming increasingly common objects of comparison. Complete bacterial genomes are reported almost everyday. This creates challenges for developing new multiple sequence alignment methods. Conventional multiple alignment methods are based on pairwise alignment and/or progressive alignment techniques. These approaches have performance problems when the number of sequences is large and when dealing with genome scale sequences. RESULTS: We present a new method of multiple sequence alignment, called MISHIMA (Method for Inferring Sequence History In terms of Multiple Alignment), that does not depend on pairwise sequence comparison. A new algorithm is used to quickly find rare oligonucleotide sequences shared by all sequences. Divide and conquer approach is then applied to break the sequences into fragments that can be aligned independently by an external alignment program. These partial alignments are assembled together to form a complete alignment of the original sequences. CONCLUSIONS: MISHIMA provides improved performance compared to the commonly used multiple alignment methods. As an example, six complete genome sequences of bacteria species Helicobacter pylori (about 1.7 Mb each) were successfully aligned in about 6 hours using a single PC.

Estimation of bacterial species phylogeny through oligonucleotide frequency distances.

Classification of bacteria is mainly based on sequence comparisons of certain homologous genes such as 16S rRNA. Recently there are challenges to classify bacteria using oligonucleotide frequency pattern of nonhomologous sequences. However, the evolutionary significance of oligonucleotides longer than tetra-nucleotide is not studied well. We performed phylogenetic analysis by using the Euclidean distances calculated from the di to deca-nucleotide frequencies in bacterial genomes, and compared these oligonucleotide frequency-based tree topologies with those for 16S rRNA gene and concatenated seven genes. When oligonucleotide frequency-based trees were constructed for bacterial species with similar GC content, their topologies at genus and family level were congruent with those based on homologous genes. Our results suggest that oligonucleotide frequency is useful not only for classification of bacteria, but also for estimation of their phylogenetic relationships for closely related species.

Sequence Compression Benchmark (SCB) database-A comprehensive evaluation of reference-free compressors for FASTA-formatted sequences.

BACKGROUND: Nearly all molecular sequence databases currently use gzip for data compression. Ongoing rapid accumulation of stored data calls for a more efficient compression tool. Although numerous compressors exist, both specialized and general-purpose, choosing one of them was difficult because no comprehensive analysis of their comparative advantages for sequence compression was available. FINDINGS: We systematically benchmarked 430 settings of 48 compressors (including 29 specialized sequence compressors and 19 general-purpose compressors) on representative FASTA-formatted datasets of DNA, RNA, and protein sequences. Each compressor was evaluated on 17 performance measures, including compression strength, as well as time and memory required for compression and decompression. We used 27 test datasets including individual genomes of various sizes, DNA and RNA datasets, and standard protein datasets. We summarized the results as the Sequence Compression Benchmark database (SCB database, http://kirr.dyndns.org/sequence-compression-benchmark/), which allows custom visualizations to be built for selected subsets of benchmark results. CONCLUSION: We found that modern compressors offer a large improvement in compactness and speed compared to gzip. Our benchmark allows compressors and their settings to be compared using a variety of performance measures, offering the opportunity to select the optimal compressor on the basis of the data type and usage scenario specific to a particular application.