CiteULike: hpaces's library [111 articles]

A Metagenomic Survey of Microbes in Honey Bee Colony Collapse Disorder.

Diana L Cox-Foster — 2007-09-09T13:26:15-00:00

Science (6 September 2007)

In colony collapse disorder (CCD), honey bee colonies inexplicably lose their workers. CCD has resulted in a loss of 50 to 90% of colonies in beekeeping operations across the United States. The observation that irradiated combs from affected colonies can be repopulated with naive bees suggests that infection may contribute to CCD. We used an unbiased metagenomic approach to survey microflora in CCD hives, normal hives, and imported royal jelly. Candidate pathogens were screened for significance of association with CCD by examination of samples collected from several sites over a period of 3 years. One organism, Israeli acute paralysis virus of bees (IAPV), was strongly correlated with CCD.

MEGAN analysis of metagenomic data.

DH Huson — 2008-02-14T19:21:34-00:00

Genome Res, Vol. 17, No. 3. (March 2007), pp. 377-386.

Metagenomics is the study of the genomic content of a sample of organisms obtained from a common habitat using targeted or random sequencing. Goals include understanding the extent and role of microbial diversity. The taxonomical content of such a sample is usually estimated by comparison against sequence databases of known sequences. Most published studies use the analysis of paired-end reads, complete sequences of environmental fosmid and BAC clones, or environmental assemblies. Emerging sequencing-by-synthesis technologies with very high throughput are paving the way to low-cost random "shotgun" approaches. This paper introduces MEGAN, a new computer program that allows laptop analysis of large metagenomic data sets. In a preprocessing step, the set of DNA sequences is compared against databases of known sequences using BLAST or another comparison tool. MEGAN is then used to compute and explore the taxonomical content of the data set, employing the NCBI taxonomy to summarize and order the results. A simple lowest common ancestor algorithm assigns reads to taxa such that the taxonomical level of the assigned taxon reflects the level of conservation of the sequence. The software allows large data sets to be dissected without the need for assembly or the targeting of specific phylogenetic markers. It provides graphical and statistical output for comparing different data sets. The approach is applied to several data sets, including the Sargasso Sea data set, a recently published metagenomic data set sampled from a mammoth bone, and several complete microbial genomes. Also, simulations that evaluate the performance of the approach for different read lengths are presented.

Microbial Population Structures in the Deep Marine Biosphere

Julie Huber — 2007-10-05T07:57:32-00:00

Science, Vol. 318, No. 5847. (5 October 2007), pp. 97-100.

The analytical power of environmental DNA sequences for modeling microbial ecosystems depends on accurate assessments of population structure, including diversity (richness) and relative abundance (evenness). We investigated both aspects of population structure for microbial communities at two neighboring hydrothermal vents by examining the sequences of more than 900,000 microbial small-subunit ribosomal RNA amplicons. The two vent communities have different population structures that reflect local geochemical regimes. Descriptions of archaeal diversity were nearly exhaustive, but despite collecting an unparalleled number of sequences, statistical analyses indicated additional bacterial diversity at every taxonomic level. We predict that hundreds of thousands of sequences will be necessary to capture the vast diversity of microbial communities, and that different patterns of evenness for both high- and low-abundance taxa may be important in defining microbial ecosystem dynamics. 10.1126/science.1146689

Metagenomic analysis of the microbial community associated with the coral Porites astreoides

Linda Wegley — 2007-07-22T08:39:16-00:00

Environmental Microbiology, Vol. 0, No. 0. (0000), pp. ???-???.

Summary The coral holobiont is a dynamic assemblage of the coral animal, zooxanthellae, endolithic algae and fungi, Bacteria,Archaea and viruses. Zooxanthellae and some Bacteria form relatively stable and species-specific associations with corals. Other associations are less specific; coral-associated Archaea differ from those in the water column, but the same archaeal species may be found on different coral species. It has been hypothesized that the coral animal can adapt to differing ecological niches by switching its microbial associates. In the case of corals and zooxanthellae, this has been termed adaptive bleaching and it has important implications for carbon cycling within the coral holobiont and ultimately the survival of coral reefs. However, the roles of other components of the coral holobiont are essentially unknown. To better understand these other coral associates, a fractionation procedure was used to separate the microbes, mitochondria and viruses from the coral animal cells and zooxanthellae. The resulting metagenomic DNA was sequenced using pyrosequencing. Fungi, Bacteria and phage were the most commonly identified organisms in the metagenome. Three of the four fungal phyla were represented, including a wide diversity of fungal genes involved in carbon and nitrogen metabolism, suggesting that the endolithic community is more important than previously appreciated. In particular, the data suggested that endolithic fungi could be converting nitrate and nitrite to ammonia, which would enable fixed nitrogen to cycle within the coral holobiont. The most prominent bacterial groups were Proteobacteria (68%), Firmicutes (10%), Cyanobacteria (7%) and Actinobacteria (6%). Functionally, the bacterial community was primarily heterotrophic and included a number of pathways for the degradation of aromatic compounds, the most abundant being the homogentisate pathway. The most abundant phage family was the ssDNA Microphage and most of the eukaryotic viruses were most closely related to those known to infect aquatic organisms. This study provides a metabolic and taxonomic snapshot of microbes associated with the reef-building coral Porites astreoides and presents a basis for understanding how coral-microbial interactions structure the holobiont and coral reefs.

Structural divergence between the human and chimpanzee genomes

Kehrer-Sawatzki — 2007-01-23T22:56:06-00:00

Human Genetics, Vol. 120, No. 6. (January 2007), pp. 759-778.

Chromosomal speciation of humans and chimpanzees revisited: studies of DNA divergence within inverted regions.

JM Szamalek — 2008-07-17T07:44:29-00:00

Cytogenetic and genome research, Vol. 116, No. 1-2. (2007), pp. 53-60.

The human and chimpanzee karyotypes are distinguishable in terms of nine pericentric inversions. According to the recombination suppression model of speciation, these inversions could have promoted the process of parapatric speciation between hominoid populations ancestral to chimpanzees and humans. Were recombination suppression to have occurred in inversion heterozygotes, gene flow would have been reduced, resulting in the accumulation of genetic incompatibilities leading to reproductive isolation and eventual speciation. In an attempt to detect the molecular signature of such events, the sequence divergence of non-coding DNA was compared between humans and chimpanzees. Precise knowledge of the locations of the inversion breakpoints permitted accurate discrimination between inverted and non-inverted regions. Contrary to the predictions of the recombination suppression model, sequence divergence was found to be lower in inverted chromosomal regions as compared to non-inverted regions, albeit with borderline statistical significance. Thus, no signature of recombination suppression resulting from inversion heterozygosity appears to be detectable by analysis of extant human and chimpanzee non-coding DNA. The precise delineation of the inversion breakpoints may nevertheless still prove helpful in identifying potential speciation-relevant genes within the inverted regions.

Molecular mechanisms of chromosomal rearrangement during primate evolution.

H Kehrer-Sawatzki — 2008-07-17T07:40:01-00:00

Chromosome research : an international journal on the molecular, supramolecular and evolutionary aspects of chromosome biology, Vol. 16, No. 1. (2008), pp. 41-56.

Breakpoint analysis of the large chromosomal rearrangements which have occurred during primate evolution promises to yield new insights into the underlying mechanisms of mutagenesis. Comparison of these evolutionary breakpoints with those that are disease-associated in humans, and which occur during either meiotic or mitotic cell division, should help to identify basic mechanistic similarities as well as differences. It has recently become clear that segmental duplications (SDs) have had a very significant impact on genome plasticity during primate evolution. In comparisons of the human and chimpanzee genomes, SDs have been found in flanking regions of 70-80% of inversions and approximately 40% of deletions/duplications. A strong spatial association between primate-specific breakpoints and SDs has also become evident from comparisons of human with other mammalian genomes. The lineage-specific hyperexpansion of certain SDs observed in the genomes of human, chimpanzee, gorilla and gibbon is indicative of the intrinsic instability of some SDs in primates. However, since many primate-specific breakpoints map to regions lacking SDs, but containing interspersed high-copy repetitive sequence elements such as SINEs, LINEs, LTRs, alpha-satellites and (AT)( n ) repeats, we may infer that a range of different molecular mechanisms have probably been involved in promoting chromosomal breakage during the evolution of primate genomes.

DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis

Brad Sherman — 2007-11-03T09:38:04-00:00

BMC Bioinformatics, Vol. 8 (02 November 2007), 426.

BACKGROUND: Due to the complex and distributed nature of biological research, our current biological knowledge is spread over many redundant annotation databases maintained by many independent groups. Analysts usually need to visit many of these bioinformatics databases in order to integrate comprehensive annotation information for their genes, which becomes one of the bottlenecks, particularly for the analytic task associated with a large gene list. Thus, a highly centralized and ready-to-use gene-annotation knowledgebase is in demand for high throughput gene functional analysis. DESCRIPTION: The DAVID Knowledgebase is built around the DAVID Gene Concept, a single-linkage method to agglomerate tens of millions of gene/protein identifiers from a variety of public genomic resources into DAVID gene clusters. The grouping of such identifiers improves the cross-reference capability, particularly across NCBI and UniProt systems, enabling more than 40 publicly available functional annotation sources to be comprehensively integrated and centralized by the DAVID gene clusters. The simple, pair-wise, text format files which make up the DAVID Knowledgebase are freely downloadable for various data analysis uses. In addition, a well organized web interface allows users to query different types of heterogeneous annotations in a high-throughput manner. CONCLUSION: The DAVID Knowledgebase is designed to facilitate high throughput gene functional analysis. For a given gene list, it not only provides the quick accessibility to a wide range of heterogeneous annotation data in a centralized location, but also enriches the level of biological information for an individual gene. Moreover, the entire DAVID Knowledgebase is freely downloadable or searchable at http://david.abcc.ncifcrf.gov/knowledgebase/.

Recent acceleration of human adaptive evolution

John Hawks — 2007-12-26T21:54:38-00:00

Proceedings of the National Academy of Sciences, Vol. 104, No. 52. (26 December 2007), pp. 20753-20758.

Genomic surveys in humans identify a large amount of recent positive selection. Using the 3.9-million HapMap SNP dataset, we found that selection has accelerated greatly during the last 40,000 years. We tested the null hypothesis that the observed age distribution of recent positively selected linkage blocks is consistent with a constant rate of adaptive substitution during human evolution. We show that a constant rate high enough to explain the number of recently selected variants would predict (i) site heterozygosity at least 10-fold lower than is observed in humans, (ii) a strong relationship of heterozygosity and local recombination rate, which is not observed in humans, (iii) an implausibly high number of adaptive substitutions between humans and chimpanzees, and (iv) nearly 100 times the observed number of high-frequency linkage disequilibrium blocks. Larger populations generate more new selected mutations, and we show the consistency of the observed data with the historical pattern of human population growth. We consider human demographic growth to be linked with past changes in human cultures and ecologies. Both processes have contributed to the extraordinarily rapid recent genetic evolution of our species. 10.1073/pnas.0707650104

Beyond species: Il'ya Ivanov and his experiments on cross-breeding humans and anthropoid apes.

K Rossiianov — 2008-07-17T07:24:18-00:00

Science in context, Vol. 15, No. 2. (June 2002), pp. 277-316.

The possibility of crossing humans with other anthropoid species has been discussed in fiction as well as in scientific literature during the twentieth century. Professor Il'ya Ivanov's attempt to achieve this was critical for the beginning of organized primate research in the Soviet Union, and remains one of the most interesting and controversial experiments that was ever done on non-human primates. The possibility of removing the boundary that separates humans from other animal species, apes in particular, is loaded with important political meaning and violates cultural and ethical taboos. The history of Ivanov's scientific experiment thus helps to reveal some of the twentieth-century's important cultural conventions and hidden assumptions about human nature, species, and social hierarchy.

Beyond eugenics: the forgotten scandal of hybridizing humans and apes.

A Etkind — 2008-07-16T10:08:09-00:00

Studies in history and philosophy of biological and biomedical sciences, Vol. 39, No. 2. (June 2008), pp. 205-210.

This paper examines the available evidence on one of the most radical ideas in the history of eugenics and utopianism. In the mid-1920s, the zoology professor Ilia Ivanov submitted to the Soviet government a project for hybridizing humans and apes by means of artificial insemination. He received substantial financing and organized expeditions to Africa to catch apes for his experiments. His project caused an international sensation. The American Association for the Advancement of Atheism announced its fund-raising campaign to support Ivanov's project but gave it a scandalously racist interpretation. Ivanov's own motivation remained unclear, as did the motivation of those in the Bolshevik government who supported Ivanov until his arrest in 1930. This paper discusses three hypothetical reasons for Ivanov's adventure: first, hybridization between humans and apes, should it be successful, would support the atheist propaganda of the Bolsheviks; second, regardless of the success of hybridization, Ivanov would catch and bring to Russia apes, which were necessary for the rejuvenation programs that were fashionable among the Bolshevik elite; and third, hybridization, should it be successful, would pave the way to the New Socialist Man whose 'construction by scientific means' was the official purpose of the Bolsheviks. Ivanov's ideas were arguably important for the American proponent of reform eugenics, Herman Muller, and for the Soviet anthropologist Boris Porshnev.

Diversity of microRNAs in human and chimpanzee brain.

Eugene Berezikov — 2006-11-01T08:27:30-00:00

Nat Genet (29 October 2006)

We used massively parallel sequencing to compare the microRNA (miRNA) content of human and chimpanzee brains, and we identified 447 new miRNA genes. Many of the new miRNAs are not conserved beyond primates, indicating their recent origin, and some miRNAs seem species specific, whereas others are expanded in one species through duplication events. These data suggest that evolution of miRNAs is an ongoing process and that along with ancient, highly conserved miRNAs, there are a number of emerging miRNAs.

Mammalian mirtron genes.

E Berezikov — 2007-11-02T06:59:17-00:00

Mol Cell, Vol. 28, No. 2. (26 October 2007), pp. 328-336.

Mirtrons are alternative precursors for microRNA biogenesis that were recently described in invertebrates. These short hairpin introns use splicing to bypass Drosha cleavage, which is otherwise essential for the generation of canonical animal microRNAs. Using computational and experimental strategies, we now establish that mammals have mirtrons as well. We identified 3 mirtrons that are well conserved and expressed in diverse mammals, 16 primate-specific mirtrons, and 46 candidates supported by limited cloning evidence in primates. As with some fly and worm mirtrons, the existence of well-conserved mammalian mirtrons indicates their relatively ancient incorporation into endogenous regulatory pathways. However, as worms, flies, and mammals each have different sets of mirtrons, we hypothesize that different animals may have independently evolved the capacity for this hybrid small RNA pathway. This notion is supported by our observation of several clade-specific features of mammalian and invertebrate mirtrons.

Inter- and intra-combinatorial regulation by transcription factors and microRNAs

Yiming Zhou — 2007-10-31T09:34:56-00:00

BMC Genomics, Vol. 8 (30 October 2007), 396.

Identification and Classification of Conserved RNA Secondary Structures in the Human Genome.

Jakob Skou Pedersen — 2006-04-24T15:27:06-00:00

PLoS Comput Biol, Vol. 2, No. 4. (21 April 2006)

The discoveries of microRNAs and riboswitches, among others, have shown functional RNAs to be biologically more important and genomically more prevalent than previously anticipated. We have developed a general comparative genomics method based on phylogenetic stochastic context-free grammars for identifying functional RNAs encoded in the human genome and used it to survey an eight-way genome-wide alignment of the human, chimpanzee, mouse, rat, dog, chicken, zebra-fish, and puffer-fish genomes for deeply conserved functional RNAs. At a loose threshold for acceptance, this search resulted in a set of 48,479 candidate RNA structures. This screen finds a large number of known functional RNAs, including 195 miRNAs, 62 histone 3'UTR stem loops, and various types of known genetic recoding elements. Among the highest-scoring new predictions are 169 new miRNA candidates, as well as new candidate selenocysteine insertion sites, RNA editing hairpins, RNAs involved in transcript auto regulation, and many folds that form singletons or small functional RNA families of completely unknown function. While the rate of false positives in the overall set is difficult to estimate and is likely to be substantial, the results nevertheless provide evidence for many new human functional RNAs and present specific predictions to facilitate their further characterization.

Primate MicroRNAs miR-220 and miR-492 Lie within Processed Pseudogenes

Devor — 2006-03-09T15:58:57-00:00

Journal of Heredity, Vol. 97, No. 2. (15 March 2006), pp. 186-190.

MicroRNAs (miRNAs) are a new and abundant class of small, noncoding RNAs. To date, the evolutionary history of most of these loci appears to be marked by duplication and divergence. The ultimate origin of miRNAs remains an open question. A survey of the genomic context of more than 300 human miRNA loci revealed that two primate-specific miRNAs, miR-220 and miR-492, each lie within a processed pseudogene. In silico and in vitro examinations of these two loci suggest that this is a rare phenomenon requiring the juxtaposition of a specific combination of factors. Thus it appears that, while processed pseudogenes are good candidates for miRNA incubators, it is unlikely that more than a very small percentage of new miRNAs arise this way.

Initial sequence and comparative analysis of the cat genome

Joan Pontius — 2007-11-02T00:25:10-00:00

Genome Res., Vol. 17, No. 11. (1 November 2007), pp. 1675-1689.

The genome sequence (1.9-fold coverage) of an inbred Abyssinian domestic cat was assembled, mapped, and annotated with a comparative approach that involved cross-reference to annotated genome assemblies of six mammals (human, chimpanzee, mouse, rat, dog, and cow). The results resolved chromosomal positions for 663,480 contigs, 20,285 putative feline gene orthologs, and 133,499 conserved sequence blocks (CSBs). Additional annotated features include repetitive elements, endogenous retroviral sequences, nuclear mitochondrial (numt) sequences, micro-RNAs, and evolutionary breakpoints that suggest historic balancing of translocation and inversion incidences in distinct mammalian lineages. Large numbers of single nucleotide polymorphisms (SNPs), deletion insertion polymorphisms (DIPs), and short tandem repeats (STRs), suitable for linkage or association studies were characterized in the context of long stretches of chromosome homozygosity. In spite of the light coverage capturing [~]65% of euchromatin sequence from the cat genome, these comparative insights shed new light on the tempo and mode of gene/genome evolution in mammals, promise several research applications for the cat, and also illustrate that a comparative approach using more deeply covered mammals provides an informative, preliminary annotation of a light (1.9-fold) coverage mammal genome sequence. 10.1101/gr.6380007

'Chumanzee' evolution: the urge to diverge and merge

Todd Disotell — 2006-11-24T22:00:28-00:00

Genome Biology, Vol. 7 (24 November 2006), 240.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

Hildegard Kehrer-Sawatzki — 2007-01-16T17:58:20-00:00

Human Mutation, Vol. 28, No. 2. (2007), pp. 99-130.

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome. Hum Mutat 28(2), 99-130, 2007. © 2006 Wiley-Liss, Inc.

The new paradigm of flow cell sequencing

Robert Holt — 2008-06-02T17:49:08-00:00

Genome Res., Vol. 18, No. 6. (1 June 2008), pp. 839-846.

DNA sequencing is in a period of rapid change, in which capillary sequencing is no longer the technology of choice for most ultra-high-throughput applications. A new generation of instruments that utilize primed synthesis in flow cells to obtain, simultaneously, the sequence of millions of different DNA templates has changed the field. We compare and contrast these new sequencing platforms in terms of stage of development, instrument configuration, template format, sequencing chemistry, throughput capability, operating cost, data handling issues, and error models. While these platforms outperform capillary instruments in terms of bases per day and cost per base, the short length of sequence reads obtained from most instruments and the limited number of samples that can be run simultaneously imposes some practical constraints on sequencing applications. However, recently developed methods for paired-end sequencing and for array-based direct selection of desired templates from complex mixtures extend the utility of these platforms for genome analysis. Given the ever increasing demand for DNA sequence information, we can expect continuous improvement of this new generation of instruments and their eventual replacement by even more powerful technology. 10.1101/gr.073262.107

The impact of next-generation sequencing technology on genetics

Elaine Mardis — 2008-02-13T12:35:19-00:00

Trends in Genetics, Vol. In Press, Corrected Proof

If one accepts that the fundamental pursuit of genetics is to determine the genotypes that explain phenotypes, the meteoric increase of DNA sequence information applied toward that pursuit has nowhere to go but up. The recent introduction of instruments capable of producing millions of DNA sequence reads in a single run is rapidly changing the landscape of genetics, providing the ability to answer questions with heretofore unimaginable speed. These technologies will provide an inexpensive, genome-wide sequence readout as an endpoint to applications ranging from chromatin immunoprecipitation, mutation mapping and polymorphism discovery to noncoding RNA discovery. Here I survey next-generation sequencing technologies and consider how they can provide a more complete picture of how the genome shapes the organism.

From The Cover: An algorithm for progressive multiple alignment of sequences with insertions

Ari Loytynoja — 2006-11-25T22:12:28-00:00

PNAS, Vol. 102, No. 30. (26 July 2005), pp. 10557-10562.

Dynamic programming algorithms guarantee to find the optimal alignment between two sequences. For more than a few sequences, exact algorithms become computationally impractical, and progressive algorithms iterating pairwise alignments are widely used. These heuristic methods have a serious drawback because pairwise algorithms do not differentiate insertions from deletions and end up penalizing single insertion events multiple times. Such an unrealistically high penalty for insertions typically results in overmatching of sequences and an underestimation of the number of insertion events. We describe a modification of the traditional alignment algorithm that can distinguish insertion from deletion and avoid repeated penalization of insertions and illustrate this method with a pair hidden Markov model that uses an evolutionary scoring function. In comparison with a traditional progressive alignment method, our algorithm infers a greater number of insertion events and creates gaps that are phylogenetically consistent but spatially less concentrated. Our results suggest that some insertion/deletion "hot spots" may actually be artifacts of traditional alignment algorithms. 10.1073/pnas.0409137102

Phylogeny-Aware Gap Placement Prevents Errors in Sequence Alignment and Evolutionary Analysis

Ari Loytynoja — 2008-06-20T03:07:52-00:00

Science, Vol. 320, No. 5883. (20 June 2008), pp. 1632-1635.

Genetic sequence alignment is the basis of many evolutionary and comparative studies, and errors in alignments lead to errors in the interpretation of evolutionary information in genomes. Traditional multiple sequence alignment methods disregard the phylogenetic implications of gap patterns that they create and infer systematically biased alignments with excess deletions and substitutions, too few insertions, and implausible insertion-deletion-event histories. We present a method that prevents these systematic errors by recognizing insertions and deletions as distinct evolutionary events. We show theoretically and practically that this improves the quality of sequence alignments and downstream analyses over a wide range of realistic alignment problems. These results suggest that insertions and sequence turnover are more common than is currently thought and challenge the conventional picture of sequence evolution and mechanisms of functional and structural changes. 10.1126/science.1158395

Millimeter-scale genetic gradients and community-level molecular convergence in a hypersaline microbial mat

Victor Kunin — 2008-06-05T08:57:43-00:00

Mol Syst Biol, Vol. 4 (3 June 2008)

To investigate the extent of genetic stratification in structured microbial communities, we compared the metagenomes of 10 successive layers of a phylogenetically complex hypersaline mat from Guerrero Negro, Mexico. We found pronounced millimeter-scale genetic gradients that were consistent with the physicochemical profile of the mat. Despite these gradients, all layers displayed near-identical and acid-shifted isoelectric point profiles due to a molecular convergence of amino-acid usage, indicating that hypersalinity enforces an overriding selective pressure on the mat community.

KEGG for linking genomes to life and the environment

Minoru Kanehisa — 2007-12-13T06:01:37-00:00

Nucl. Acids Res. (12 December 2007), gkm882.

KEGG (http://www.genome.jp/kegg/) is a database of biological systems that integrates genomic, chemical and systemic functional information. KEGG provides a reference knowledge base for linking genomes to life through the process of PATHWAY mapping, which is to map, for example, a genomic or transcriptomic content of genes to KEGG reference pathways to infer systemic behaviors of the cell or the organism. In addition, KEGG provides a reference knowledge base for linking genomes to the environment, such as for the analysis of drug-target relationships, through the process of BRITE mapping. KEGG BRITE is an ontology database representing functional hierarchies of various biological objects, including molecules, cells, organisms, diseases and drugs, as well as relationships among them. KEGG PATHWAY is now supplemented with a new global map of metabolic pathways, which is essentially a combined map of about 120 existing pathway maps. In addition, smaller pathway modules are defined and stored in KEGG MODULE that also contains other functional units and complexes. The KEGG resource is being expanded to suit the needs for practical applications. KEGG DRUG contains all approved drugs in the US and Japan, and KEGG DISEASE is a new database linking disease genes, pathways, drugs and diagnostic markers. 10.1093/nar/gkm882

Vir-Mir db: prediction of viral microRNA candidate hairpins.

SC Li — 2008-05-11T15:08:43-00:00

Nucleic acids research, Vol. 36, No. Database issue. (January 2008)

MicroRNAs have been found in various organisms and play essential roles in gene expression regulation of many critical cellular processes. Large-scale computational prediction of miRNAs has been conducted for many organisms using known genomic sequences; however, there has been no such effort for the thousands of known viral genomes. Some viruses utilize existing host cellular pathways for their own benefit. Furthermore, viruses are capable of encoding miRNAs and using them to repress host genes. Thus, identifying potential miRNAs in all viral genomes would be valuable to virologists who study virus-host interactions. Based on our previously reported hairpin secondary structure and feature selection filters, we have examined the 2266 available viral genome sequences for putative miRNA hairpins and identified 33 691 hairpin candidates in 1491 genomes. Evaluation of the system performance indicated that our discovery pipeline exhibited 84.4% sensitivity. We established an interface for users to query the predicted viral miRNA hairpins based on taxonomic classification, and a host target gene prediction service based on the RNAhybrid program and the 3'-UTR gene sequences of human, mouse, rat, zebrafish, rice and Arabidopsis. The viral miRNA prediction database (Vir-Mir) can be accessed via http://alk.ibms.sinica.edu.tw.

What is a gene, post-ENCODE? History and updated definition.

MB Gerstein — 2007-06-15T21:49:06-00:00

Genome Res, Vol. 17, No. 6. (June 2007), pp. 669-681.

While sequencing of the human genome surprised us with how many protein-coding genes there are, it did not fundamentally change our perspective on what a gene is. In contrast, the complex patterns of dispersed regulation and pervasive transcription uncovered by the ENCODE project, together with non-genic conservation and the abundance of noncoding RNA genes, have challenged the notion of the gene. To illustrate this, we review the evolution of operational definitions of a gene over the past century-from the abstract elements of heredity of Mendel and Morgan to the present-day ORFs enumerated in the sequence databanks. We then summarize the current ENCODE findings and provide a computational metaphor for the complexity. Finally, we propose a tentative update to the definition of a gene: A gene is a union of genomic sequences encoding a coherent set of potentially overlapping functional products. Our definition sidesteps the complexities of regulation and transcription by removing the former altogether from the definition and arguing that final, functional gene products (rather than intermediate transcripts) should be used to group together entities associated with a single gene. It also manifests how integral the concept of biological function is in defining genes.

Genome analysis of the platypus reveals unique signatures of evolution

2008-05-07T23:43:22-00:00

Nature, Vol. 453, No. 7192. (May 2008), pp. 175-183.

Mapping and sequencing of structural variation from eight human genomes

Jeffrey Kidd — 2008-04-30T19:31:59-00:00

Nature, Vol. 453, No. 7191., pp. 56-64.

Poor prognosis in carcinoma is associated with a gene expression signature of aberrant PTEN tumor suppressor pathway activity.

LH Saal — 2007-11-20T16:05:43-00:00

Proc Natl Acad Sci U S A, Vol. 104, No. 18. (1 May 2007), pp. 7564-7569.

Pathway-specific therapy is the future of cancer management. The oncogenic phosphatidylinositol 3-kinase (PI3K) pathway is frequently activated in solid tumors; however, currently, no reliable test for PI3K pathway activation exists for human tumors. Taking advantage of the observation that loss of PTEN, the negative regulator of PI3K, results in robust activation of this pathway, we developed and validated a microarray gene expression signature for immunohistochemistry (IHC)-detectable PTEN loss in breast cancer (BC). The most significant signature gene was PTEN itself, indicating that PTEN mRNA levels are the primary determinant of PTEN protein levels in BC. Some PTEN IHC-positive BCs exhibited the signature of PTEN loss, which was associated to moderately reduced PTEN mRNA levels cooperating with specific types of PIK3CA mutations and/or amplification of HER2. This demonstrates that the signature is more sensitive than PTEN IHC for identifying tumors with pathway activation. In independent data sets of breast, prostate, and bladder carcinoma, prediction of pathway activity by the signature correlated significantly to poor patient outcome. Stathmin, encoded by the signature gene STMN1, was an accurate IHC marker of the signature and had prognostic significance in BC. Stathmin was also pathway-pharmacodynamic in vitro and in vivo. Thus, the signature or its components such as stathmin may be clinically useful tests for stratification of patients for anti-PI3K pathway therapy and monitoring therapeutic efficacy. This study indicates that aberrant PI3K pathway signaling is strongly associated with metastasis and poor survival across carcinoma types, highlighting the enormous potential impact on patient survival that pathway inhibition could achieve.

What is principal component analysis?

Markus Ringnér — 2008-03-09T04:13:08-00:00

Nature Biotechnology, Vol. 26, No. 3., pp. 303-304.

Sequencing and analysis of Neanderthal genomic DNA.

JP Noonan — 2006-11-23T23:51:10-00:00

Science, Vol. 314, No. 5802. (17 November 2006), pp. 1113-1118.

Our knowledge of Neanderthals is based on a limited number of remains and artifacts from which we must make inferences about their biology, behavior, and relationship to ourselves. Here, we describe the characterization of these extinct hominids from a new perspective, based on the development of a Neanderthal metagenomic library and its high-throughput sequencing and analysis. Several lines of evidence indicate that the 65,250 base pairs of hominid sequence so far identified in the library are of Neanderthal origin, the strongest being the ascertainment of sequence identities between Neanderthal and chimpanzee at sites where the human genomic sequence is different. These results enabled us to calculate the human-Neanderthal divergence time based on multiple randomly distributed autosomal loci. Our analyses suggest that on average the Neanderthal genomic sequence we obtained and the reference human genome sequence share a most recent common ancestor approximately 706,000 years ago, and that the human and Neanderthal ancestral populations split approximately 370,000 years ago, before the emergence of anatomically modern humans. Our finding that the Neanderthal and human genomes are at least 99.5% identical led us to develop and successfully implement a targeted method for recovering specific ancient DNA sequences from metagenomic libraries. This initial analysis of the Neanderthal genome advances our understanding of the evolutionary relationship of Homo sapiens and Homo neanderthalensis and signifies the dawn of Neanderthal genomics.

Reconstructing phylogenies and phenotypes: a molecular view of human evolution

Bradley — 2008-04-05T06:23:35-00:00

Journal of Anatomy, Vol. 212, No. 4. (April 2008), pp. 337-353.

ALLPATHS: De novo assembly of whole-genome shotgun microreads.

Jonathan Butler — 2008-03-18T00:16:38-00:00

Genome Res (13 March 2008)

New DNA sequencing technologies deliver data at dramatically lower costs but demand new analytical methods to take full advantage of the very short reads that they produce. We provide an initial, theoretical solution to the challenge of de novo assembly from whole-genome shotgun "microreads." For 11 genomes of sizes up to 39 Mb, we generated high-quality assemblies from 80x coverage by paired 30-base simulated reads modeled after real Illumina-Solexa reads. The bacterial genomes of Campylobacter jejuni and Escherichia coli assemble optimally, yielding single perfect contigs, and larger genomes yield assemblies that are highly connected and accurate. Assemblies are presented in a graph form that retains intrinsic ambiguities such as those arising from polymorphism, thereby providing information that has been absent from previous genome assemblies. For both C. jejuni and E. coli, this assembly graph is a single edge encompassing the entire genome. Larger genomes produce more complicated graphs, but the vast majority of the bases in their assemblies are present in long edges that are nearly always perfect. We describe a general method for genome assembly that can be applied to all types of DNA sequence data, not only short read data, but also conventional sequence reads.

An evaluation of human protein-protein interaction data in the public domain

Suresh Mathivanan — 2007-05-24T13:23:24-00:00

BMC Bioinformatics, Vol. 7, No. Suppl 5. (2006)

BACKGROUND:Protein-protein interaction (PPI) databases have become a major resource for investigating biological networks and pathways in cells. A number of publicly available repositories for human PPIs are currently available. Each of these databases has their own unique features with a large variation in the type and depth of their annotations.RESULTS:We analyzed the major publicly available primary databases that contain literature curated PPI information for human proteins. This included BIND, DIP, HPRD, IntAct, MINT, MIPS, PDZBase and Reactome databases. The number of binary non-redundant human PPIs ranged from 101 in PDZBase and 346 in MIPS to 11,367 in MINT and 36,617 in HPRD. The number of genes annotated with at least one interactor was 9,427 in HPRD, 4,975 in MINT, 4,614 in IntAct, 3,887 in BIND and <1,000 in the remaining databases. The number of literature citations for the PPIs included in the databases was 43,634 in HPRD, 11,480 in MINT, 10,331 in IntAct, 8,020 in BIND and <2,100 in the remaining databases.CONCLUSION:Given the importance of PPIs, we suggest that submission of PPIs to repositories be made mandatory by scientific journals at the time of manuscript submission as this will minimize annotation errors, promote standardization and help keep the information up to date. We hope that our analysis will help guide biomedical scientists in selecting the most appropriate database for their needs especially in light of the dramatic differences in their content.

The Diploid Genome Sequence of an Individual Human

Samuel Levy — 2007-09-04T11:15:47-00:00

PLoS Biology, Vol. 5, No. 10. (1 October 2007), e254.

Presented here is a genome sequence of an individual human. It was produced from ∼32 million random DNA fragments, sequenced by Sanger dideoxy technology and assembled into 4,528 scaffolds, comprising 2,810 million bases (Mb) of contiguous sequence with approximately 7.5-fold coverage for any given region. We developed a modified version of the Celera assembler to facilitate the identification and comparison of alternate alleles within this individual diploid genome. Comparison of this genome and the National Center for Biotechnology Information human reference assembly revealed more than 4.1 million DNA variants, encompassing 12.3 Mb. These variants (of which 1,288,319 were novel) included 3,213,401 single nucleotide polymorphisms (SNPs), 53,823 block substitutions (2–206 bp), 292,102 heterozygous insertion/deletion events (indels)(1–571 bp), 559,473 homozygous indels (1–82,711 bp), 90 inversions, as well as numerous segmental duplications and copy number variation regions. Non-SNP DNA variation accounts for 22% of all events identified in the donor, however they involve 74% of all variant bases. This suggests an important role for non-SNP genetic alterations in defining the diploid genome structure. Moreover, 44% of genes were heterozygous for one or more variants. Using a novel haplotype assembly strategy, we were able to span 1.5 Gb of genome sequence in segments >200 kb, providing further precision to the diploid nature of the genome. These data depict a definitive molecular portrait of a diploid human genome that provides a starting point for future genome comparisons and enables an era of individualized genomic information.

DNA sequencing: bench to bedside and beyond.

Clyde A Hutchison — 2007-09-19T06:05:24-00:00

Nucleic Acids Res (12 September 2007)

Fifteen years elapsed between the discovery of the double helix (1953) and the first DNA sequencing (1968). Modern DNA sequencing began in 1977, with development of the chemical method of Maxam and Gilbert and the dideoxy method of Sanger, Nicklen and Coulson, and with the first complete DNA sequence (phage varphiX174), which demonstrated that sequence could give profound insights into genetic organization. Incremental improvements allowed sequencing of molecules >200 kb (human cytomegalovirus) leading to an avalanche of data that demanded computational analysis and spawned the field of bioinformatics. The US Human Genome Project spurred sequencing activity. By 1992 the first 'sequencing factory' was established, and others soon followed. The first complete cellular genome sequences, from bacteria, appeared in 1995 and other eubacterial, archaebacterial and eukaryotic genomes were soon sequenced. Competition between the public Human Genome Project and Celera Genomics produced working drafts of the human genome sequence, published in 2001, but refinement and analysis of the human genome sequence will continue for the foreseeable future. New 'massively parallel' sequencing methods are greatly increasing sequencing capacity, but further innovations are needed to achieve the 'thousand dollar genome' that many feel is prerequisite to personalized genomic medicine. These advances will also allow new approaches to a variety of problems in biology, evolution and the environment.

Discovering Motifs in Ranked Lists of DNA Sequences

Eran Eden — 2007-03-30T07:44:02-00:00

PLoS Computational Biology, Vol. 3, No. 3. (1 March 2007), e39.

Computational methods for discovery of sequence elements that are enriched in a target set compared with a background set are fundamental in molecular biology research. One example is the discovery of transcription factor binding motifs that are inferred from ChIP–chip (chromatin immuno-precipitation on a microarray) measurements. Several major challenges in sequence motif discovery still require consideration: (i) the need for a principled approach to partitioning the data into target and background sets; (ii) the lack of rigorous models and of an exact p-value for measuring motif enrichment; (iii) the need for an appropriate framework for accounting for motif multiplicity; (iv) the tendency, in many of the existing methods, to report presumably significant motifs even when applied to randomly generated data. In this paper we present a statistical framework for discovering enriched sequence elements in ranked lists that resolves these four issues. We demonstrate the implementation of this framework in a software application, termed DRIM (discovery of rank imbalanced motifs), which identifies sequence motifs in lists of ranked DNA sequences. We applied DRIM to ChIP–chip and CpG methylation data and obtained the following results. (i) Identification of 50 novel putative transcription factor (TF) binding sites in yeast ChIP–chip data. The biological function of some of them was further investigated to gain new insights on transcription regulation networks in yeast. For example, our discoveries enable the elucidation of the network of the TF ARO80. Another finding concerns a systematic TF binding enhancement to sequences containing CA repeats. (ii) Discovery of novel motifs in human cancer CpG methylation data. Remarkably, most of these motifs are similar to DNA sequence elements bound by the Polycomb complex that promotes histone methylation. Our findings thus support a model in which histone methylation and CpG methylation are mechanistically linked. Overall, we demonstrate that the statistical framework embodied in the DRIM software tool is highly effective for identifying regulatory sequence elements in a variety of applications ranging from expression and ChIP–chip to CpG methylation data. DRIM is publicly available at http://bioinfo.cs.technion.ac.il/drim.

Ringo - an R/Bioconductor package for analyzing ChIP-chip readouts

Joern Toedling — 2007-06-27T13:35:34-00:00

BMC Bioinformatics, Vol. 8 (26 June 2007), 221.

Multiple alignment of protein sequences with repeats and rearrangements

Phuong — 2006-12-06T11:41:10-00:00

Nucleic Acids Research, Vol. 34, No. 20. (November 2006), pp. 5932-5942.

Automated generation of heuristics for biological sequence comparison

Guy Slater — 2007-05-29T13:09:02-00:00

BMC Bioinformatics, Vol. 6, No. 1. (2005)

BACKGROUND:Exhaustive methods of sequence alignment are accurate but slow, whereas heuristic approaches run quickly, but their complexity makes them more difficult to implement. We introduce bounded sparse dynamic programming (BSDP) to allow rapid approximation to exhaustive alignment. This is used within a framework whereby the alignment algorithms are described in terms of their underlying model, to allow automated development of efficient heuristic implementations which may be applied to a general set of sequence comparison problems.RESULTS:The speed and accuracy of this approach compares favourably with existing methods. Examples of its use in the context of genome annotation are given.CONCLUSIONS:This system allows rapid implementation of heuristics approximating to many complex alignment models, and has been incorporated into the freely available sequence alignment program, exonerate.

Chaos game representation for comparison of whole genomes

Jijoy Joseph — 2007-12-22T18:31:31-00:00

BMC Bioinformatics, Vol. 7, No. 1. (2006)

BACKGROUND:Chaos game representation of genome sequences has been used for visual representation of genome sequence patterns as well as alignment-free comparisons of sequences based on oligonucleotide frequencies. However the potential of this representation for making alignment-based comparisons of whole genome sequences has not been exploited.RESULTS:We present here a fast algorithm for identifying all local alignments between two long DNA sequences using the sequence information contained in CGR points. The local alignments can be depicted graphically in a dot-matrix plot or in text form, and the significant similarities and differences between the two sequences can be identified. We demonstrate the method through comparison of whole genomes of several microbial species. Given two closely related genomes we generate information on mismatches, insertions, deletions and shuffles that differentiate the two genomes.CONCLUSION:Addition of the possibility of large scale sequence alignment to the repertoire of alignment-free sequence analysis applications of chaos game representation, positions CGR as a powerful sequence analysis tool.

The ENCODE (ENCyclopedia Of DNA Elements) Project.

E Encode — 2005-04-04T06:53:10-00:00

Science, Vol. 306, No. 5696. (22 October 2004), pp. 636-640.

The ENCyclopedia Of DNA Elements (ENCODE) Project aims to identify all functional elements in the human genome sequence. The pilot phase of the Project is focused on a specified 30 megabases (approximately 1%) of the human genome sequence and is organized as an international consortium of computational and laboratory-based scientists working to develop and apply high-throughput approaches for detecting all sequence elements that confer biological function. The results of this pilot phase will guide future efforts to analyze the entire human genome.

Genomic plasticity in prokaryotes: the case of the square haloarchaeon.

S Cuadros-Orellana — 2007-12-02T15:44:29-00:00

ISME J, Vol. 1, No. 3. (July 2007), pp. 235-245.

The variability in genome content among closely related strains of prokaryotes has been one of the most remarkable discoveries of genomics. One way to approach the description of this so-called pan-genome is to compare one reference strain genome with metagenomic sequences from the environment. We have applied this approach to one extreme aquatic habitat, saturated brines in a solar saltern. The genome of Haloquadratum walsbyi strain DSM 16790 was compared to an environmental metagenome obtained from the exact site of its isolation. This approach revealed that some regions of the strain genome were scarcely represented in the metagenome. Here we have analyzed these genomic islands (GI) in the genome of DSM 16790 and compared them with the complete sequence of some fosmids from the environmental library. Two of the islands, GI 2 and GI 4, overlapped with two large guanine and cytosine (GC)-rich regions that showed evidence of high variability through mobile elements. GI 3 seemed to be a phage or phage-remnant acquired by the reference genome, but not present in most environmental lineages. Most differential gene content was related to small molecule transport and detection, probably reflecting adaptation to different pools of organic nutrients. GI 1 did not possess traces of mobile elements and had normal GC content. This island contained the main cluster of cell envelope glycoproteins and the variability found was different from the other GIs. Rather than containing different genes it consisted of homologs with low similarity. This variation might reflect a phage evasion strategy.The ISME Journal (2007) 1, 235-245; doi:10.1038/ismej.2007.35; published online 31 May 2007.

Long-term evolution of transposable elements

Arnaud Le Rouzic — 2007-11-27T03:46:29-00:00

Proceedings of the National Academy of Sciences (26 November 2007), 0705238104.

Transposable elements are often considered parasitic DNA sequences, able to invade the genome of their host thanks to their self-replicating ability. This colonization process has been extensively studied, both theoretically and experimentally, but their long-term coevolution with the genomes is still poorly understood. In this work, we aim to challenge previous population genetics models by considering features of transposable elements as quantitative, rather than discrete, variables. We also describe more realistic transposable element dynamics by accounting for the variability of the insertion effect, from deleterious to adaptive, as well as mutations leading to a loss of transposition activity and to nonautonomous copies. Individual-based simulations of the behavior of a transposable-element family over several thousand generations show different ways in which active or inactive copies can be maintained for a very long time. Results reveal an unexpected impact of genetic drift on the "junk DNA" content of the genome and strongly question the likelihood of the sustainable long-term stable transposition-selection equilibrium on which numerous previous works were based. 10.1073/pnas.0705238104

Herpesviral Protein Networks and Their Interaction with the Human Proteome

Peter Uetz — 2006-01-17T17:52:01-00:00

Science, Vol. 311, No. 5758. (13 January 2006), pp. 239-242.

The comprehensive yeast two-hybrid analysis of intraviral protein interactions in two members of the herpesvirus family, Kaposi sarcoma-associated herpesvirus (KSHV) and varicella-zoster virus (VZV), revealed 123 and 173 interactions, respectively. Viral protein interaction networks resemble single, highly coupled modules, whereas cellular networks are organized in separate functional submodules. Predicted and experimentally verified interactions between KSHV and human proteins were used to connect the viral interactome into a prototypical human interactome and to simulate infection. The analysis of the combined system showed that the viral network adopts cellular network features and that protein networks of herpesviruses and possibly other intracellular pathogens have distinguishing topologies.

Comparative analysis of the complete genome sequence of the plant growth–promoting bacterium Bacillus amyloliquefaciens FZB42

Xiao Chen — 2007-09-11T23:43:23-00:00

Nature Biotechnology, Vol. 25, No. 9. (19 August 2007), pp. 1007-1014.

Genome-wide transcriptional changes in Streptococcus gordonii in response to competence signaling peptide.

MM Vickerman — 2007-12-11T09:06:57-00:00

J Bacteriol, Vol. 189, No. 21. (November 2007), pp. 7799-7807.

Streptococcus gordonii is a primary colonizer of the multispecies biofilm on tooth surfaces forming dental plaque and a potential agent of endocarditis. The recent completion of the genome sequence of the naturally competent strain Challis allowed the design of a spotted oligonucleotide microarray to examine a genome-wide response of this organism to environmental stimuli such as signal peptides. Based on temporal responses to synthetic competence signaling peptide (CSP) as indicated by transformation frequencies, the S. gordonii transcriptome was analyzed at various time points after CSP exposure. Microarray analysis identified 35 candidate early genes and 127 candidate late genes that were up-regulated at 5 and 15 min, respectively; these genes were often grouped in clusters. Results supported published findings on S. gordonii competence, showing up-regulation of 12 of 16 genes that have been reported to affect transformation frequencies in this species. Comparison of CSP-induced S. gordonii transcriptomes to results published for Streptococcus pneumoniae strains identified both conserved and species-specific genes. Putative intergenic regulatory sites, such as the conserved combox sequence thought to be a binding site for competence sigma factor, were found preceding S. gordonii late responsive genes. In contrast, S. gordonii early CSP-responsive genes were not preceded by the direct repeats found in S. pneumoniae. These studies provide the first insights into a genome-wide transcriptional response of an oral commensal organism. They offer an extensive analysis of transcriptional changes that accompany competence in S. gordonii and form a basis for future intra- and interspecies comparative analyses of this ecologically important phenotype.

The complete genome sequence of Campylobacter jejuni strain 81116 (NCTC11828).

BM Pearson — 2007-12-11T09:04:30-00:00

J Bacteriol, Vol. 189, No. 22. (November 2007), pp. 8402-8403.

Campylobacter jejuni is a major human enteric pathogen that displays genetic variability via genomic reorganization and phase variation. This variability can adversely affect the outcomes and reproducibility of experiments. C. jejuni strain 81116 (NCTC11828) has been suggested to be a genetically stable strain (G. Manning, B. Duim, T. Wassenaar, J. A. Wagenaar, A. Ridley, and D. G. Newell, Appl. Environ. Microbiol. 67:1185-1189, 2001), is amenable to genetic manipulation, and is infective for chickens. Here we report the finished annotated genome sequence of C. jejuni strain 81116.

Paradoxical DNA Repair and Peroxide Resistance Gene Conservation in Bacillus pumilus SAFR-032.

J Gioia — 2007-12-11T09:03:25-00:00

PLoS ONE, Vol. 2, No. 9. (2007)

BACKGROUND: Bacillus spores are notoriously resistant to unfavorable conditions such as UV radiation, gamma-radiation, H(2)O(2), desiccation, chemical disinfection, or starvation. Bacillus pumilus SAFR-032 survives standard decontamination procedures of the Jet Propulsion Lab spacecraft assembly facility, and both spores and vegetative cells of this strain exhibit elevated resistance to UV radiation and H(2)O(2) compared to other Bacillus species. PRINCIPAL FINDINGS: The genome of B. pumilus SAFR-032 was sequenced and annotated. Lists of genes relevant to DNA repair and the oxidative stress response were generated and compared to B. subtilis and B. licheniformis. Differences in conservation of genes, gene order, and protein sequences are highlighted because they potentially explain the extreme resistance phenotype of B. pumilus. The B. pumilus genome includes genes not found in B. subtilis or B. licheniformis and conserved genes with sequence divergence, but paradoxically lacks several genes that function in UV or H(2)O(2) resistance in other Bacillus species. SIGNIFICANCE: This study identifies several candidate genes for further research into UV and H(2)O(2) resistance. These findings will help explain the resistance of B. pumilus and are applicable to understanding sterilization survival strategies of microbes.