Optimizing read mapping to reference genomes to determine composition and species prevalence in microbial communities

The Human Microbiome Project (HMP) aims to characterize the microbial communities of 18 body sites from healthy individuals. To accomplish this, the HMP generated two types of shotgun data: reference shotgun sequences isolated from different anatomical sites on the human body and shotgun metagenomic sequences from the microbial communities of each site. The alignment strategy for characterizing these metagenomic communities using available reference sequence is important to the success of HMP data analysis. Six next-generation aligners were used to align a community of known composition against a database comprising reference organisms known to be present in that community. All aligners report nearly complete genome coverage (>97%) for strains with over 6X depth of coverage, however they differ in speed, memory requirement and ease of use issues such as database size limitations and supported mapping strategies. The selected aligner was tested across a range of parameters to maximize sensitivity while maintaining a low false positive rate. We found that constraining alignment length had more impact on sensitivity than does constraining similarity in all cases tested. However, when reference species were replaced with phylogenetic neighbors, similarity begins to play a larger role in detection. We also show that choosing the top hit randomly when multiple, equally strong mappings are available increases overall sensitivity at the expense of taxonomic resolution. The results of this study identified a strategy that was used to map over 3 tera-bases of microbial sequence against a database of more than 5,000 reference genomes in just over a month.     

Learning about human population history from ancient and modern genomes

Genome-wide data, both from SNP arrays and from complete genome sequencing, are becoming increasingly abundant and are now even available from extinct hominins. These data are providing new insights into population history; in particular, when combined with model-based analytical approaches, genome-wide data allow direct testing of hypotheses about population history. For example, genome-wide data from both contemporary populations and extinct hominins strongly support a single dispersal of modern humans from Africa, followed by two archaic admixture events: one with Neanderthals somewhere outside Africa and a second with Denisovans that (so far) has only been detected in New Guinea. These new developments promise to reveal new stories about human population history, without having to resort to storytelling.     

Optical mapping of DNA: single-molecule-based methods for mapping genomes

The technologies associated with DNA sequencing are rapidly evolving. Indeed, single-molecule DNA sequencing strategies are cheaper and faster than ever before. Despite this progress, every sequencing platform to date relies on reading the genome in small, abstract fragments, typically of less than 1000 bases in length. The overarching aim of the optical map is to complement the information derived from DNA sequencing by providing long-range context on which these short sequence reads can be built. This is typically done using an enzyme to target and modify at short DNA sequences of, say, six bases in length throughout the genome. By accurately placing these short pieces of sequence on long genomic DNA fragments, up to several millions of bases in length, a scaffold for sequence assembly can be obtained. This review focuses on three enzymatic approaches to optical mapping. Optical mapping was first developed using restriction enzymes to sequence-specifically cleave DNA that is immobilized on a surface. More recently, nicking enzymes have found application in the sequence-specific fluorescent labeling of DNA for optical mapping. Such covalent modification allows the DNA to be imaged in solution, and this, in combination with developing nanofluidic technologies, is enabling new high-throughput approaches to mapping. And, finally, this review will discuss the recent development of mapping with subdiffraction-limit precision using methyltransferase enzymes to label the DNA with an ultrahigh density.     

Molecular phylogenetics employing modern and ancient DNA

Comparative studies of DNA in recent populations and characterisation of ancient hereditary material have contributed very interesting facts to our understanding of evolution of modern mankind. Analysis of DNA homology in related species, assessment of mutations and polymorphisms in various populations and new DNA sequence data from prehistoric finds allowed - via sophisticated DNA extraction techniques, PCR, sequencing and digitalised processing of genetic information - insights into possible roots of Homo sapiens and related species, migration patterns and ancient cultural habits, thus enrhing the palaeoanthropological discipline. However, a presentation of this development would not be complete without pointing towards the methodological limitations and manifold presentations burdened with artifacts, data misinterpretation and unjustified conclusions. Presently, this modern field of research is in its consolidation phase and new parameters for quality control and authentication are being implemented to avoid spectacular but unfounded reports. It is expected that most of the problems connected to old biomolecules may be closely related to fossilisation parameters. The future challenge will be the full understanding of the complex and multi-faceted processes underlying diagenesis, including the elucidation of nucleic acid postmortem damage".     

Pathocoenoses ancient and modern

Problems concerning the concept of biocoenosis in ecology (the antecedent of the pathocoenosis concept) are discussed first of all. Six main problems are identified: the problem of emergent properties of ecological communities; the problem of ambiguity; the problem of heterogeneity; the boundary problem; the problem of retrospective differential diagnosis; and the problem of explaining change over time. The rest of the paper gives illustrations of these problems in relation to human pathogens drawn mainly from the interactions of malaria with other diseases, particularly but not exclusively in the Mediterranean world, from antiquity through to modern times.     

GRIM-Filter: Fast seed location filtering in DNA read mapping using processing-in-memory technologies

Background

Seed location filtering is critical in DNA read mapping, a process where billions of DNA fragments (reads) sampled from a donor are mapped onto a reference genome to identify genomic variants of the donor. State-of-the-art read mappers 1) quickly generate possible mapping locations for seeds (i.e., smaller segments) within each read, 2) extract reference sequences at each of the mapping locations, and 3) check similarity between each read and its associated reference sequences with a computationally-expensive algorithm (i.e., sequence alignment) to determine the origin of the read. A seed location filter comes into play before alignment, discarding seed locations that alignment would deem a poor match. The ideal seed location filter would discard all poor match locations prior to alignment such that there is no wasted computation on unnecessary alignments.

Results

We propose a novel seed location filtering algorithm, GRIM-Filter, optimized to exploit 3D-stacked memory systems that integrate computation within a logic layer stacked under memory layers, to perform processing-in-memory (PIM). GRIM-Filter quickly filters seed locations by 1) introducing a new representation of coarse-grained segments of the reference genome, and 2) using massively-parallel in-memory operations to identify read presence within each coarse-grained segment. Our evaluations show that for a sequence alignment error tolerance of 0.05, GRIM-Filter 1) reduces the false negative rate of filtering by 5.59x–6.41x, and 2) provides an end-to-end read mapper speedup of 1.81x–3.65x, compared to a state-of-the-art read mapper employing the best previous seed location filtering algorithm.

Conclusion

GRIM-Filter exploits 3D-stacked memory, which enables the efficient use of processing-in-memory, to overcome the memory bandwidth bottleneck in seed location filtering. We show that GRIM-Filter significantly improves the performance of a state-of-the-art read mapper. GRIM-Filter is a universal seed location filter that can be applied to any read mapper. We hope that our results provide inspiration for new works to design other bioinformatics algorithms that take advantage of emerging technologies and new processing paradigms, such as processing-in-memory using 3D-stacked memory devices.

     

Accessory DNA in the genomes of representatives of the Escherichia coli reference collection

Different strains of the Escherichia coli reference collection (ECOR) differ widely in chromosomal size. To analyze the nature of the differential gene pool carried by different strains, we have followed an approach in which random amplified polymorphic DNA (RAPD) was used to generate several PCR fragments. Those present in some but not all the strains were screened by hybridization to assess their distribution throughout the ECOR collection. Thirteen fragments with various degrees of occurrence were sequenced. Three of them corresponded to RAPD markers of widespread distribution. Of these, two were housekeeping genes shown by hybridization to be present in all the E. coli strains and in Salmonella enterica LT2; the third fragment contained a paralogous copy of dnaK with widespread, but not global, distribution. The other 10 RAPD markers were found in only a few strains. However, hybridization results demonstrated that four of them were actually present in a large selection of the ECOR collection (between 42 and 97% of the strains); three of these fragments contained open reading frames associated with phages or plasmids known in E. coli K-12. The remaining six fragments were present in only between one and four strains; of these, four fragments showed no similarity to any sequence in the databases, and the other two had low but significant similarity to a protein involved in the Klebsiella capsule synthesis and to RNA helicases of archaeal genomes, respectively. Their percent GC, dinucleotide content, and codon adaptation index suggested an exogenous origin by horizontal transfer. These results can be interpreted as reflecting the presence of a large pool of strain-specific genes, whose origin could be outside the species boundaries.     

Evolution of ancient satellite DNAs in sturgeon genomes

This study characterizes a repetitive DNA family of sequences in sturgeon, the PstI satellite DNA. We have found a high degree of preservation for these sequences, which are present in all 13 species analyzed, including within the genera Acipenser, Huso, and Scaphirhynchus of the family Acipenseridae. This is one of the most ancient satellite DNAs found to date, because it has been estimated to be more than 100 million years old. Alternatively, to the current view that most satellite DNAs are species-specific or preserved in a few closely related species, the PstI family and other previously characterized sturgeon satellite DNA, the HindIII, represent the most fascinating exceptions to the rapid sequence change usually undergone by satellite DNAs. Here, we compare the evolutionary pattern of these two satellite DNA families, PstI and HindIII, which differ markedly in length, sequence, and nucleotide composition. We have found that, in contrast to the situation in most other living beings, a high degree of preservation, a slow sequence change rate and slowed concerted evolution, appears to be a general rule for sturgeon satellite DNAs. The possible causes for all these features are discussed in the light of the evolutionary specifics found within these ancient organisms.     

Classical catalase: ancient and modern

This review describes the historical difficulties in devising a kinetically satisfactory mechanism for the classical catalase after its identification as a unique catalytic entity in 1902 and prior to the breakthrough 1947 analysis by Chance and co-workers which led to the identification of peroxide compounds I and II. The role of protons in the formation of these two ferryl complexes is discussed and current problems of inhibitory ligand and hydrogen donor binding at the active site are outlined, especially the multiple roles involving formate or formic acid. A previous mechanism of NADPH-dependent catalase protection against substrate inhibition is defended. A revised model linking the catalytic ('catalatic') action and the one-electron side reactions involving compound II is suggested. And it is concluded that, contrary to an idea proposed in 1963 that eukaryotic catalase might be a 'fossil enzyme', current thinking gives it a central role in the redox protective processes of long term importance for human and other eukaryotic and prokaryotic life.     

Waging war on physical inactivity: using modern molecular ammunition against an ancient enemy.

A hypothesis is presented based on a coalescence of anthropological estimations of Homo sapiens' phenotypes in the Late Paleolithic era 10,000 years ago, with Darwinian natural selection synergized with Neel's idea of the so-called thrifty gene. It is proposed that humans inherited genes that were evolved to support a physically active lifestyle. It is further postulated that physical inactivity in sedentary societies directly contributes to multiple chronic health disorders. Therefore, it is imperative to identify the underlying genetic and cellular/biochemical bases of why sedentary living produces chronic health conditions. This will allow society to improve its ability to effect beneficial lifestyle changes and hence improve the overall quality of living. To win the war against physical inactivity and the myriad of chronic health conditions produced because of physical inactivity, a multifactorial approach is needed, which includes successful preventive medicine, drug development, optimal target selection, and efficacious clinical therapy. All of these approaches require a thorough understanding of fundamental biology and how the dysregulated molecular circuitry caused by physical inactivity produces clinically overt disease. The purpose of this review is to summarize the vast armamentarium at our disposal in the form of the extensive scientific basis underlying how physical inactivity affects at least 20 of the most deadly chronic disorders. We hope that this information will provide readers with a starting point for developing additional strategies of their own in the ongoing war against inactivity-induced chronic health conditions.     

Molecular genetic analysis of the polymorphic apolipoprotein E in modern and ancient human DNA samples

In this report, methodical bases for the molecular genetic analysis of the three common apolipoprotein E alleles APOE*2, APOE*3 and APOE*4 in DNA isolated from ancient human skeletal remains are described. Considering that ancient DNA target regions for amplification are generally quite small, the detection method is based on short amplification products in the range from 71 bp to 75 bp. The applicability of the modified method for APOE genotyping was examined in modern human DNA samples.     

Merging ancient and modern DNA: extinct seabird taxon rediscovered in the North Tasman Sea

Ancient DNA has revolutionized the way in which evolutionary biologists research both extinct and extant taxa, from the inference of evolutionary history to the resolution of taxonomy. Here, we present, to our knowledge, the first study to report the rediscovery of an ‘extinct’ avian taxon, the Tasman booby (Sula tasmani), using classical palaeontological data combined with ancient and modern DNA data. Contrary to earlier work, we show an overlap in size between fossil and modern birds in the North Tasman Sea (classified currently as S. tasmani and Sula dactylatra fullagari, respectively). In addition, we show that Holocene fossil birds have mitochondrial control region sequences that are identical to those found in modern birds. These results indicate that the Tasman booby is not an extinct taxon: S. dactylatra fullagari O''Brien & Davies, 1990 is therefore a junior synonym of Sula tasmani van Tets, Meredith, Fullagar & Davidson, 1988 and all North Tasman Sea boobies should be known as S. d. tasmani. In addition to reporting the rediscovery of an extinct avian taxon, our study highlights the need for researchers to be cognizant of multidisciplinary approaches to understanding taxonomy and past biodiversity.