Diversity of zooplankton assemblages at Station L4, Western English Channel, measured using next generation DNA sequencing

Background: Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel Next Generation Sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify the richness and diversity of a mixed zooplankton assemblage from a productive monitoring site in the Western English Channel. Methodology/Principle Findings: Plankton WP2 replicate net hauls (200 µm) were taken at the Western Channel Observatory long-term monitoring station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,042 sequences were obtained for all samples. The sequences clustered in to 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 138 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 75 taxonomic groups. Conclusions: The percentage of OTUs assigned to major eukaryotic taxonomic groups broadly aligns between the metagenetic and morphological analysis and are dominated by Copepoda. However, the metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for estimating diversity and species richness of zooplankton communities.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Automated parallel DNA sequencing on multiple channel microchips

We report automated DNA sequencing in 16-channel microchips. A microchip prefilled with sieving matrix is aligned on a heating plate affixed to a movable platform. Samples are loaded into sample reservoirs by using an eight-tip pipetting device, and the chip is docked with an array of electrodes in the focal plane of a four-color scanning detection system. Under computer control, high voltage is applied to the appropriate reservoirs in a programmed sequence that injects and separates the DNA samples. An integrated four-color confocal fluorescent detector automatically scans all 16 channels. The system routinely yields more than 450 bases in 15 min in all 16 channels. In the best case using an automated base-calling program, 543 bases have been called at an accuracy of >99%. Separations, including automated chip loading and sample injection, normally are completed in less than 18 min. The advantages of DNA sequencing on capillary electrophoresis chips include uniform signal intensity and tolerance of high DNA template concentration. To understand the fundamentals of these unique features we developed a theoretical treatment of cross-channel chip injection that we call the differential concentration effect. We present experimental evidence consistent with the predictions of the theory.

DNA sequencing and genotyping by transcriptional synthesis of chain-terminated RNA ladders and MALDI-TOF mass spectrometry

Sets of RNA ladders can be synthesized by transcription of a bacteriophage-encoded RNA polymerase using 3′-deoxynucleotides as chain terminators. These ladders can be used for sequencing of DNA. Using a nicked form of phage SP6 RNA polymerase in this study substantially enhanced yields of transcriptional sequencing ladders. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) of chain-terminated RNA ladders allowed DNA sequence determination of up to 56 nt. It is also demonstrated that A→G and C→T variations in heterozygous and homozygous samples can be unambiguously identified by the mass spectrometric analysis. As a step towards single-tube sequencing reactions, α-thiotriphosphate nucleotide analogs were used to overcome problems caused by chain terminator-independent, premature termination and by the small mass difference between natural pyrimidine nucleotides.

DNA sequencing by delayed extraction-matrix-assisted laser desorption/ionization time of flight mass spectrometry.

Matrix-assisted laser desorption/ionization (MALDI) time of flight mass spectrometry was used to detect and order DNA fragments generated by Sanger dideoxy cycle sequencing. This was accomplished by improving the sensitivity and resolution of the MALDI method using a delayed ion extraction technique (DE-MALDI). The cycle sequencing chemistry was optimized to produce as much as 100 fmol of each specific dideoxy terminated fragment, generated from extension of a 13-base primer annealed on 40- and 50-base templates. Analysis of the resultant sequencing mixture by DE-MALDI identified the appropriate termination products. The technique provides a new non-gel-based method to sequence DNA which may ultimately have considerable speed advantages over traditional methodologies.

A DNA sequencing strategy that requires only five bases of known terminal sequence for priming.

We have previously reported an enhanced version of sequencing by hybridization (SBH), termed positional SBH (PSBH). PSBH uses partially duplex probes containing single-stranded 3' overhangs, instead of simple single-stranded probes. Stacking interactions between the duplex probe and a single-stranded target allow us to reduce the probe sizes required to 5-base single-stranded overhangs. Here we demonstrate the use of PSBH to capture relatively long single-stranded DNA targets and perform standard solid-state Sanger sequencing on these primer-template complexes without ligation. Our results indicate that only 5 bases of known terminal sequence are required for priming. In addition, the partially duplex probes have the ability to capture their specific target from a mixture of five single-stranded targets with different 3'-terminal sequences. This indicates the potential utility of the PSBH approach to sequence mixtures of DNA targets without prior purification.

Fluorescence energy transfer dye-labeled primers for DNA sequencing and analysis.

Fluorescent dye-labeled DNA primers have been developed that exploit fluorescence energy transfer (ET) to optimize the absorption and emission properties of the label. These primers carry a fluorescein derivative at the 5' end as a common donor and other fluorescein and rhodamine derivatives attached to a modified thymidine residue within the primer sequence as acceptors. Adjustment of the donor-acceptor spacing through the placement of the modified thymidine in the primer sequence allowed generation of four primers, all having strong absorption at a common excitation wavelength (488 nm) and fluorescence emission maxima of 525, 555, 580, and 605 nm. The ET efficiency of these primers ranges from 65% to 97%, and they exhibit similar electrophoretic mobilities by gel electrophoresis. With argon-ion laser excitation, the fluorescence of the ET primers and of the DNA sequencing fragments generated with ET primers is 2- to 6-fold greater than that of the corresponding primers or fragments labeled with single dyes. The higher fluorescence intensity of the ET primers allows DNA sequencing with one-fourth of the DNA template typically required when using T7 DNA polymerase. With single-stranded M13mp18 DNA as the template, a typical sequencing reaction with ET primers on a commercial sequencer provided DNA sequences with 99.8% accuracy in the first 500 bases. ET primers should be generally useful in the development of other multiplex DNA sequencing and analysis methods.

Highly accurate RNA and DNA sequencing: Application to longstanding questions in aging and cancer

Thesis (Ph.D.)--University of Washington, 2016-06

The use and implications of ribosomal DNA sequencing for the discrimination of digenean species

In just over a decade, the use of molecular approaches for the recognition of parasites has become commonplace. For trematodes, the internal transcribed spacer region of ribosomal DNA (ITS rDNA) has become the default region of choice. Here, we review the findings of 63 studies that report ITS rDNA sequence data for about 155 digenean species from 19 families, and then review the levels of variation that have been reported and how the variation has been interpreted. Overall, complete ITS sequences (or ITS1 or ITS2 regions alone) usually distinguish trematode species clearly, including combinations for which morphology gives ambiguous results. Closely related species may have few base differences and in at least one convincing case the ITS2 sequences of two good species are identical. In some cases, the ITS1 region gives greater resolution than the ITS2 because of the presence of variable repeat units that are generally lacking in the ITS2. Intraspecific variation is usually low and frequently apparently absent. Information on geographical variation of digeneans is limited but at least some of the reported variation probably reflects the presence of multiple species. Despite the accepted dogma that concerted evolution makes the individual representative of the entire species, a significant number of studies have reported at least some intraspecific variation. The significance of such variation is difficult to assess a posteriori, but it seems likely that identification and sequencing errors account for some of it and failure to recognise separate species may also be significant. Some reported variation clearly requires further analysis. The use of a yardstick to determine when separate species should be recognised is flawed. Instead, we argue that consistent genetic differences that are associated with consistent morphological or biological traits should be considered the marker for separate species. We propose a generalised approach to the use of rDNA to distinguish trematode species.

Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era

The class Kinetoplastea encompasses both free-living and parasitic species from a wide range of hosts. Several representatives of this group are responsible for severe human diseases and for economic losses in agriculture and livestock. While this group encompasses over 30 genera, most of the available information has been derived from the vertebrate pathogenic genera Leishmania and Trypanosoma. Recent studies of the previously neglected groups of Kinetoplastea indicated that the actual diversity is much higher than previously thought. This article discusses the known segment of kinetoplastid diversity and how gene-directed Sanger sequencing and next-generation sequencing methods can help to deepen our knowledge of these interesting protists.

The impact of whey protein consumption and exercise on the composition and diversity of the gut microbiota: a high through-put DNA sequencing approach

Advances in culture independent technologies over the last decade have highlighted the pivotal role which the gut microbiota plays in maintaining human health. Conversely, perturbations to the composition or actions of the ‘normal/functioning’ microbiota have been frequently associated with the pathogenesis of several disease states. Therefore the selective modulation of enteric microbial communities represents a viable target for the development of novel treatments for such diseases. Notably, while bovine whey proteins and exercise have been shown to positively influence several physiological processes, such as energy balance, their effect on the composition or functionality of the gut microbiota remains largely unknown. In this thesis, a variety of ex vivo, murine and human models are used in conjunction with high-throughput DNA sequencing-based analysis to provide valuable and novel insights into the impact of both whey proteins and exercise on enteric microbial communities. Overall the results presented in this thesis highlight that the consumption both whey protein isolate (WPI), and individual component proteins of whey such as bovine serum albumin (BSA) and lactoferrin, reduce high fat diet associated body weight gain and are associated with beneficial alterations within the murine gut microbiota. Although the impact of exercise on enteric microbial communities remains less clear, it may be that longer term investigations are required for the true effect of exercise on the gut microbiota to be fully elucidated.

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

DNA methylation analysis of phenotype specific stratified Indian population

Background: DNA methylation and its perturbations are an established attribute to a wide spectrum of phenotypic variations and disease conditions. Indian traditional system practices personalized medicine through indigenous concept of distinctly descriptive physiological, psychological and anatomical features known as prakriti. Here we attempted to establish DNA methylation differences in these three prakriti phenotypes. Methods: Following structured and objective measurement of 3416 subjects, whole blood DNA of 147 healthy male individuals belonging to defined prakriti (Vata, Pitta and Kapha) between the age group of 20-30years were subjected to methylated DNA immunoprecipitation (MeDIP) and microarray analysis. After data analysis, prakriti specific signatures were validated through bisulfite DNA sequencing. Results: Differentially methylated regions in CpG islands and shores were significantly enriched in promoters/UTRs and gene body regions. Phenotypes characterized by higher metabolism (Pitta prakriti) in individuals showed distinct promoter (34) and gene body methylation (204), followed by Vata prakriti which correlates to motion showed DNA methylation in 52 promoters and 139 CpG islands and finally individuals with structural attributes (Kapha prakriti) with 23 and 19 promoters and CpG islands respectively. Bisulfite DNA sequencing of prakriti specific multiple CpG sites in promoters and 5'-UTR such as; LHX1 (Vata prakriti), SOX11 (Pitta prakriti) and CDH22 (Kapha prakriti) were validated. Kapha prakriti specific CDH22 5'-UTR CpG methylation was also found to be associated with higher body mass index (BMI). Conclusion: Differential DNA methylation signatures in three distinct prakriti phenotypes demonstrate the epigenetic basis of Indian traditional human classification which may have relevance to personalized medicine.

Caracterização de Bacilos Gram-Negativos Não Fermentadores não usuais em bacteremias pelas técnicas de Matrix-Assisted Laser Desorption IonizationTime of Flight Mass Spectrometry, sequenciamento de DNA e método fenotípico convencional

Alguns Bastonetes Gram-negativos não fermentadores (BGNNF) costumam ser considerados clinicamente pouco significantes e a sua implicação em infecções é subestimada. Devido à similaridade fenotípica, mudanças taxonômicas, baixa reatividade bioquímica e limitações nos bancos de dados em sistemas comerciais, a identificação de BGNNF é frequentemente equivocada, culminando com a denominação de diferentes micro-organismos apenas como BGNNF, por falta de melhor diferenciação. O objetivo desse estudo foi avaliar, por métodos fenotípico convencional, proteômico e molecular, a identificação de BGNNF incomuns isolados em hemoculturas de pacientes atendidos em um hospital universitário no Rio de Janeiro. Foram selecionadas 78 amostras isoladas de hemoculturas caracterizadas no laboratório clinico como BGNNF para a identificação por sequenciamento dos genes 16S RNA e recA, por um conjunto amplo de testes fenotípicos manuais e por MALDI-TOF MS. Os micro-organismos predominantes na amostragem foram genotipados pela técnica de eletroforese em gel de campo pulsado (PFGE). Pelo sequenciamento do gene 16S rRNA, a maioria das amostras (n=31; 40%) foi incluída no gênero Burkholderia, seguido de Pseudomonas stutzeri (10%) e Delftia acidovorans (4%). Os demais isolados foram agrupados em 27 diferentes espécies. O sequencimento do gene recA identificou a maioria das espécies de Burkholderia como Burkholderia contaminans (n=19; 24%). Os testes fenotípicos incluíram as 31 amostras apenas no CBc e para as outras 47 amostras, a concordância com o sequenciamento do gene 16S rRNA em nível de espécie foi de 64% (n=30) e apenas em gênero a concordância foi de 17% (n=8). A análise comparativa geral da identificação por MALDI-TOF MS com o sequenciamento do gene16S rRNA mostrou que 42% (n=33) das 78 amostras foram concordantes em nível de espécie e 45% (n=35) apenas em gênero. Excluindo as amostras do CBc, houve um aumento da concordância em nível de espécie para 60%. As discordâncias parecem ser devido às diferenças nos perfis proteicos das amostras em relação às amostras-referência do banco de dados do equipamento e podem ser aprimorados com a atualização de perfis no sistema. A análise do polimorfismo genético de B. contaminans mostrou a ausência de um clone disseminado causando surto, além da provável origem ambiental das infecções. Os setores de nefrologia e hemodiálise contribuíram com maior número de pacientes com amostras positivas (5 pacientes e 9 amostras). Os grupos clonais BcoD e BcoE foram encontrados em pacientes assistidos no mesmo setor com diferença de quatro meses (BcoD, nefrologia) e 1,5 ano (BcoE, hemodilálise), entre as culturas, respectivamente. As discordâncias entre as técnicas ocorreram principalmente devido a dificuldade de identificação das espécies do CBc. Os BGNNF incomuns são de difícil caracterização independente da metodologia usada e nenhum método por si só foi capaz de identificar todas as amostras.