Evaluation of genome skimming to detect and characterise human and livestock helminths

The identification of gastrointestinal helminth infections of humans and livestock almost exclusively relies on the detection of eggs or larvae in faeces, followed by manual counting and morphological characterisation to differentiate species using microscopy-based techniques. However, molecular approaches based on the detection and quantification of parasite DNA are becoming more prevalent, increasing the sensitivity, specificity and throughput of diagnostic assays. High-throughput sequencing, from single PCR targets through to the analysis of whole genomes, offers significant promise towards providing information-rich data that may add value beyond traditional and conventional molecular approaches; however, thus far, its utility has not been fully explored to detect helminths in faecal samples. In this study, low-depth whole genome sequencing, i.e. genome skimming, has been applied to detect and characterise helminth diversity in a set of helminth-infected human and livestock faecal material. The strengths and limitations of this approach are evaluated using three methods to characterise and differentiate metagenomic sequencing data based on (i) mapping to whole mitochondrial genomes, (ii) whole genome assemblies, and (iii) a comprehensive internal transcribed spacer 2 (ITS2) database, together with validation using quantitative PCR (qPCR). Our analyses suggest that genome skimming can successfully identify most single and multi-species infections reported by qPCR and can provide sufficient coverage within some samples to resolve consensus mitochondrial genomes, thus facilitating phylogenetic analyses of selected genera, e.g. Ascaris spp. Key to this approach is both the availability and integrity of helminth reference genomes, some of which are currently contaminated with bacterial and host sequences. The success of genome skimming of faecal DNA is dependent on the availability of vouchered sequences of helminths spanning both taxonomic and geographic diversity, together with methods to detect or amplify minute quantities of parasite nucleic acids in mixed samples.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

Mitochondrial genomes of parasitic nematodes--progress and perspectives

Mitochondria are subcellular organelles in which oxidative phosphorylation and other important biochemical functions take place within the cell. Within these organelles is a mitochondrial (mt) genome, which is distinct from, but cooperates with, the nuclear genome of the cell. Studying mt genomes has implications for various fundamental areas, including mt biochemistry, physiology and molecular biology. Importantly, the mt genome is a rich source of markers for population genetic and systematic studies. To date, more than 696 mt genomes have been sequenced for a range of metazoan organisms. However, few of these are from parasitic nematodes, despite their socioeconomic importance and the need for fundamental investigations into areas such as nematode genetics, systematics and ecology. In this article, we review knowledge and recent progress in mt genomics of parasitic nematodes, summarize applications of mt gene markers to the study of population genetics, systematics, epidemiology and evolution of key nematodes, and highlight some prospects and opportunities for future research.     

Conservation and versatility of a new set of primers for long‐PCR amplification of complete insect mitochondrial genomes based on Haematobia irritans mtDNA sequences

The amplification of complete mitochondrial genomes by long PCR (polymerase chain reaction) has been a major contribution to the large‐scale sequencing of arthropodan mitochondrial genomes. In this work, we designed six conserved long‐PCR primers to successfully recover the entire mitochondrial genome of the horn fly Haematobia irritans (Diptera: Muscidae) in two overlapping fragments. The conservation and versatility of these primers were tested for 17 other species from four major insect orders: Diptera (14), Coleoptera (1), Lepidoptera (1) and Hymenoptera (1). The amplification of complete mitochondrial genomes in orders other than Diptera suggested an even broader application of these primers, especially within the Hexapoda.     

GPMiner: an integrated system for mining combinatorial cis-regulatory elements in mammalian gene group

Background

Massively parallel sequencing technology is revolutionizing approaches to genomic and genetic research. Since its advent, the scale and efficiency of Next-Generation Sequencing (NGS) has rapidly improved. In spite of this success, sequencing genomes or genomic regions with extremely biased base composition is still a great challenge to the currently available NGS platforms. The genomes of some important pathogenic organisms like Plasmodium falciparum (high AT content) and Mycobacterium tuberculosis (high GC content) display extremes of base composition. The standard library preparation procedures that employ PCR amplification have been shown to cause uneven read coverage particularly across AT and GC rich regions, leading to problems in genome assembly and variation analyses. Alternative library-preparation approaches that omit PCR amplification require large quantities of starting material and hence are not suitable for small amounts of DNA/RNA such as those from clinical isolates. We have developed and optimized library-preparation procedures suitable for low quantity starting material and tolerant to extremely high AT content sequences.

Results

We have used our optimized conditions in parallel with standard methods to prepare Illumina sequencing libraries from a non-clinical and a clinical isolate (containing ~53% host contamination). By analyzing and comparing the quality of sequence data generated, we show that our optimized conditions that involve a PCR additive (TMAC), produces amplified libraries with improved coverage of extremely AT-rich regions and reduced bias toward GC neutral templates.

Conclusion

We have developed a robust and optimized Next-Generation Sequencing library amplification method suitable for extremely AT-rich genomes. The new amplification conditions significantly reduce bias and retain the complexity of either extremes of base composition. This development will greatly benefit sequencing clinical samples that often require amplification due to low mass of DNA starting material.     

Toward next-generation sequencing of mitochondrial genomes — Focus on parasitic worms of animals and biotechnological implications

Helminths (worms) include parasitic nematodes (roundworms) and platyhelminths (flatworms). These worms are abundant, and many of them are of agricultural, aquacultural, veterinary and medical importance and cause substantial socioeconomic losses worldwide. The genetic characterization of parasitic nematodes using advanced molecular tools is central to the diagnosis of infections and the control of parasitism. The accurate analysis of genetic variation also underpins studies of their taxonomy, epidemiology and evolutionary history. Although the nuclear genome contains suitable genetic markers (e.g., in ribosomal DNA) for the identification of many species, the large size and high variability of the mt genome consistently provides a rich source of such markers for informative systematic and epidemiological studies both within and among species. There is significant value in establishing a practical platform for the rapid sequencing, annotation and analysis of mt genomic datasets to underpin such fundamental and applied studies of parasitic worms (= helminths). In the last decade, there have been some important advances in the mt genomics of helminths, but next-generation sequencing (NGS) technologies now provide opportunities for high throughput sequencing, assembly and annotation. In this article, we provide a background on mt genomics, cover technological challenges and recent advances, and provide a perspective on future mt genome research of parasitic helminths and its fundamental scientific and biotechnological implications.     

Amplification of mitochondrial small subunit ribosomal DNA of polypores and its potential for phylogenetic analysis

There has been a systematic need to seek adequate phylogenetic markers that can be applied in phylogenetic analyses of fungal taxa at various levels. The mitochondrial small subunit ribosomal DNA (mt SSU rDNA) is generally considered to be one of the molecules that are appropriate for phylogenetic analyses at a family level. In order to obtain universal primers for polypores of Hymenomycetes, mt SSU rRNA genes were cloned from Bjerkandera adusta, Ganoderma lucidum, Phlebiopsis gigantea, and Phellinus laevigatus and their sequences were determined. Based on the conserved sequences of cloned genes from polypores and Agrocybe aegerita, PCR primers were designed for amplification and sequencing of mt SSU rDNAs. New primers allowed effective amplification and sequencing of almost full-sized genes from representative species of polypores and related species. Phylogenetic relationships were resolved quite efficiently by mt SSU rDNA sequences, and they proved to be more useful in phylogenetic reconstruction of Ganoderma than nuclear internal transcribed spacer (ITS) rDNA sequences.     

Uniparental inheritance of cpDNA and the genetic control of sexual differentiation in Chlamydomonas reinhardtii

An intriguing feature of most eukaryotes is that chloroplast (cp) and mitochondrial (mt) genomes are inherited almost exclusively from one parent. Uniparental inheritance of cp/mt genomes was long thought to be a passive outcome, based on the fact that eggs contain multiple numbers of organelles, while male gametes contribute, at best, only a few cp/mtDNA. However, the process is likely to be more dynamic because uniparental inheritance occurs in organisms that produce gametes of identical sizes (isogamous). In Chlamydomonas reinhardtii, the uniparental inheritance of cp/mt genomes is achieved by a series of mating type-controlled events that actively eliminate the mating type minus (mt−) cpDNA. The method by which Chlamydomonas selectively degrades mt− cpDNA has long fascinated researchers, and is the subject of this review.     

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: An ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8–10 length variants. The motifs of the HTR arrays are about 353–362 bp and the number of repeats ranges from 1 to 11.     

Whole genome amplification and sequencing of a Daphnia resting egg

Resting eggs banks are unique windows that allow us to directly observe shifts in population genetics, and phenotypes over time as natural populations evolve. Though a variety of planktonic organisms also produce resting stages, the keystone freshwater consumer, Daphnia, is a well‐known model for paleogenetics and resurrection ecology. Nevertheless, paleogenomic investigations are limited largely because resting eggs do not contain enough DNA for genomic sequencing. In fact, genomic studies even on extant populations include a laborious preparatory phase of batch culturing dozens of individuals to generate sufficient genomic DNA. Here, we furnish a protocol to generate whole genomes of single ephippial (resting) eggs and single daphniids. Whole genomes of single ephippial eggs and single adults were amplified using Qiagen REPLI‐g Single Cell kit reaction, followed by NEBNext Ultra DNA Library Prep Kit for library construction and Illumina sequencing. We compared the quality of the single‐egg and single‐individual amplified genomes to the standard batch genomic DNA extraction in the absence of genome amplification. At mean 20× depth, coverage was essentially identical for the amplified single individual relative to the unamplified batch extracted genome (>90% of the genome was covered and callable). Finally, while amplification resulted in the slight loss of heterozygosity for the amplified genomes, estimates were largely comparable and illustrate the utility and limitations of this approach in estimating population genetic parameters over long periods of time in natural populations of Daphnia and also other small species known to produce resting stages.     

Whole Genome Amplification for Sequencing and Applications in Conservation Genetics

Abstract: We describe a method for rapidly amplifying whole genomes via a Phi29 DNA polymerase-mediated strand displacement reaction (SDR). Genomic amplification products derived from the SDR reaction resulted in high quantities of DNA suitable for polymerase chain reaction (PCR) amplification and sequencing of mitochondrial genomes. Control region sequences of DNA derived directly from PCR amplicons of extracted DNA were identical to those derived from PCR amplification of SDR genomic DNA. Effective SDR amplification and subsequent sequencing was successful across tissues sources ranging in age from 1 year to 19 years. Strand replacement reaction genomic amplification offers a means of obtaining large quantities of DNA from small amounts of tissue.     

Automated sequencing of complete mitochondrial genomes from laser-capture microdissected samples

Mitochondrial DNA mutations have been related to both aging and a variety of diseases such as cancer. Due to the relatively small size of the genome (16 kb) and with the use of automated DNA sequencing, the entire genome can be sequenced from clinical specimens in days. We present a reliable approach to complete mitochondrial genome sequencing from laser-capture microdissected human clinical cancer specimens that overcome the inherent limitations of relatively small tissue samples and partial DNA degradation, which are unavoidable when laser-capture microdissection is used to attain pure populations of cells from heterogeneous tissues obtained from surgical procedures. The acquisition of sufficient template combined with a standard set of 18 pairs of PCR primers allows for the efficient amplification of the genome. Subsequent single-stranded amplification is performed using 36 sequencing primers, and samples are run on an ABI PRISM 3100 Genetic Analyzer. The use of this procedure should allow even investigators with little experience sequencing from clinical specimens success in complete mitochondrial genome sequencing.     

Sequencing intractable DNA to close microbial genomes

Advancement in high throughput DNA sequencing technologies has supported a rapid proliferation of microbial genome sequencing projects, providing the genetic blueprint for in-depth studies. Oftentimes, difficult to sequence regions in microbial genomes are ruled "intractable" resulting in a growing number of genomes with sequence gaps deposited in databases. A procedure was developed to sequence such problematic regions in the "non-contiguous finished" Desulfovibrio desulfuricans ND132 genome (6 intractable gaps) and the Desulfovibrio africanus genome (1 intractable gap). The polynucleotides surrounding each gap formed GC rich secondary structures making the regions refractory to amplification and sequencing. Strand-displacing DNA polymerases used in concert with a novel ramped PCR extension cycle supported amplification and closure of all gap regions in both genomes. The developed procedures support accurate gene annotation, and provide a step-wise method that reduces the effort required for genome finishing.     

Genomic sequencing of single microbial cells from environmental samples

Recently developed techniques allow genomic DNA sequencing from single microbial cells [Lasken RS: Single-cell genomic sequencing using multiple displacement amplification. Curr Opin Microbiol 2007, 10:510-516]. Here, we focus on research strategies for putting these methods into practice in the laboratory setting. An immediate consequence of single-cell sequencing is that it provides an alternative to culturing organisms as a prerequisite for genomic sequencing. The microgram amounts of DNA required as template are amplified from a single bacterium by a method called multiple displacement amplification (MDA) avoiding the need to grow cells. The ability to sequence DNA from individual cells will likely have an immense impact on microbiology considering the vast numbers of novel organisms, which have been inaccessible unless culture-independent methods could be used. However, special approaches have been necessary to work with amplified DNA. MDA may not recover the entire genome from the single copy present in most bacteria. Also, some sequence rearrangements can occur during the DNA amplification reaction. Over the past two years many research groups have begun to use MDA, and some practical approaches to single-cell sequencing have been developed. We review the consensus that is emerging on optimum methods, reliability of amplified template, and the proper interpretation of 'composite' genomes which result from the necessity of combining data from several single-cell MDA reactions in order to complete the assembly. Preferred laboratory methods are considered on the basis of experience at several large sequencing centers where >70% of genomes are now often recovered from single cells. Methods are reviewed for preparation of bacterial fractions from environmental samples, single-cell isolation, DNA amplification by MDA, and DNA sequencing.     

Sequencing of high G+C microbial genomes using the ultrafast pyrosequencing technology

Next generation pyrosequencing of high G + C content genomes still poses problems to automated sequencing and assembly processes which necessitates cost and time intensive manual work in order to finish such genomes completely. The sequencing of the high G + C actinomycete Actinoplanes sp. SE50/110 was performed with standard pyrosequencing technology (454 Life Sciences) and revealed a high number of gaps. The reasons for the introduction of gaps were analyzed on a previously known 41 kb long DNA reference sequence from Actinoplanes sp. SE50/110, hosting the acarbose biosynthesis gene cluster. Mapping of the sequencing results on the reference gene cluster sequence revealed a fragmentation into 30 contiguous sequences of different lengths. The gaps between these sequences were characterized by extremely low read coverage which strongly correlated with the G + C content in the gap regions in a negative manner. Furthermore, the gap-sequences contained strong stem-loop structures which hindered the amplification of these sequences during the emulsion PCR. Being significantly underrepresented or absent in the subsequent sequencing process, these sequences lead to weakly or uncovered genomic regions which forces the assembly algorithm to output multiple contiguous sequences instead of one finished genome. However, by applying a different pyrosequencing protocol, it was possible to sequence the complete acarbose biosynthesis gene cluster. The changes to the protocol include longer read length and addition of chemicals to the amplification chemistry, which reduces the self-annealing of DNA fragments during the amplification process and enables the complete reconstruction of high G + C content genomes without manual intervention.     

Conserved PCR primer set designing for closely-related species to complete mitochondrial genome sequencing using a sliding window-based PSO algorithm

Background

Complete mitochondrial (mt) genome sequencing is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. For long template sequencing, i.e., like the entire mtDNA, it is essential to design primers for Polymerase Chain Reaction (PCR) amplicons which are partly overlapping each other. The presented chromosome walking strategy provides the overlapping design to solve the problem for unreliable sequencing data at the 5′ end and provides the effective sequencing. However, current algorithms and tools are mostly focused on the primer design for a local region in the genomic sequence. Accordingly, it is still challenging to provide the primer sets for the entire mtDNA.

Methodology/Principal Findings

The purpose of this study is to develop an integrated primer design algorithm for entire mt genome in general, and for the common primer sets for closely-related species in particular. We introduce ClustalW to generate the multiple sequence alignment needed to find the conserved sequences in closely-related species. These conserved sequences are suitable for designing the common primers for the entire mtDNA. Using a heuristic algorithm particle swarm optimization (PSO), all the designed primers were computationally validated to fit the common primer design constraints, such as the melting temperature, primer length and GC content, PCR product length, secondary structure, specificity, and terminal limitation. The overlap requirement for PCR amplicons in the entire mtDNA is satisfied by defining the overlapping region with the sliding window technology. Finally, primer sets were designed within the overlapping region. The primer sets for the entire mtDNA sequences were successfully demonstrated in the example of two closely-related fish species. The pseudo code for the primer design algorithm is provided.

Conclusions/Significance

In conclusion, it can be said that our proposed sliding window-based PSO algorithm provides the necessary primer sets for the entire mt genome amplification and sequencing.     

Finding of Bov-B LINE retroelement in parthenogenetic and bisexual lizard species of the genus Darevskia (Lacertidae)

The Bov-B LINE retrotransposon was first discovered in Ruminantia and was long considered to be specific for this order. Later, this mobile element was described in snakes and some lizard species. Analysis of phylogenetic relationships of Bov-B LINE elements from different ruminants, snakes, and lizard species led to the suggestion on horizontal transfer of this retrotransposon from Squamata to Ruminantia. In the Squamata group, Bov-B LINE element was found in all snakes and some lizard species examined. The element was not detected in the genomes of some species of the genera Lacerta and Podarcis. In the present study, using PCR amplification and sequencing of PCR products, Bov-B LINE element was identified in the genomes of parthenogenetic and bisexual species of the genus Darevskia (Lacertidae), as well as in such species as Lacerta agilis and Zootoca vivipara, where this retrotransposon had not been not detected before.     

Beyond barcoding: A mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.