Multiplexed microsatellite recovery using massively parallel sequencing

Conservation and management of natural populations requires accurate and inexpensive genotyping methods. Traditional microsatellite, or simple sequence repeat (SSR), marker analysis remains a popular genotyping method because of the comparatively low cost of marker development, ease of analysis and high power of genotype discrimination. With the availability of massively parallel sequencing (MPS), it is now possible to sequence microsatellite-enriched genomic libraries in multiplex pools. To test this approach, we prepared seven microsatellite-enriched, barcoded genomic libraries from diverse taxa (two conifer trees, five birds) and sequenced these on one lane of the Illumina Genome Analyzer using paired-end 80-bp reads. In this experiment, we screened 6.1 million sequences and identified 356,958 unique microreads that contained di- or trinucleotide microsatellites. Examination of four species shows that our conversion rate from raw sequences to polymorphic markers compares favourably to Sanger- and 454-based methods. The advantage of multiplexed MPS is that the staggering capacity of modern microread sequencing is spread across many libraries; this reduces sample preparation and sequencing costs to less than $400 (USD) per species. This price is sufficiently low that microsatellite libraries could be prepared and sequenced for all 1373 organisms listed as 'threatened' and 'endangered' in the United States for under $0.5 M (USD).     

Optimization of virus detection in cells using massively parallel sequencing

Massively parallel sequencing (MPS)-based virus detection has potential regulatory applications. We studied the ability of one of these approaches, based on degenerate oligonucleotide primer (DOP)-polymerase chain reaction (PCR), to detect viral sequences in cell lines known to express viral genes or particles. DOP-PCR was highly sensitive for the detection of small quantities of isolated viral sequences. Detected viral sequences included nodavirus, bracovirus, and endogenous retroviruses in High Five cells, porcine circovirus type 1 and porcine endogenous retrovirus in PK15 cells, human T-cell leukemia virus 1 in MJ cells, human papillomavirus 18 in HeLa cells, human herpesvirus 8 in BCBL-1 cells, and Epstein–Barr Virus in Raji cells. Illumina sequencing (for which primers were most efficiently added using PCR) provided greater sensitivity for virus detection than Roche 454 sequencing. Analyzing nucleic acids extracted both directly from samples and from capsid-enriched preparations provided useful information. Although there are limitations of these methods, these results indicate significant promise for the combination of nonspecific PCR and MPS in identifying contaminants in clinical and biological samples, including cell lines and reagents used to produce vaccines and therapeutic products.     

BreaKmer: detection of structural variation in targeted massively parallel sequencing data using kmers

Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.     

Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes

Background

Mycoheterotrophic plants are considered to associate very specifically with fungi. Mycoheterotrophic orchids are mostly associated with ectomycorrhizal fungi in temperate regions, or with saprobes or parasites in tropical regions. Although most mycoheterotrophic orchids occur in the tropics, few studies have been devoted to them, and the main conclusions about their specificity have hitherto been drawn from their association with ectomycorrhizal fungi in temperate regions.

Results

We investigated three Asiatic Neottieae species from ectomycorrhizal forests in Thailand. We found that all were associated with ectomycorrhizal fungi, such as Thelephoraceae, Russulaceae and Sebacinales. Based on ¹³C enrichment of their biomass, they probably received their organic carbon from these fungi, as do mycoheterotrophic Neottieae from temperate regions. Moreover, ¹³C enrichment suggested that some nearby green orchids received part of their carbon from fungi too. Nevertheless, two of the three orchids presented a unique feature for mycoheterotrophic plants: they were not specifically associated with a narrow clade of fungi. Some orchid individuals were even associated with up to nine different fungi.

Conclusion

Our results demonstrate that some green and mycoheterotrophic orchids in tropical regions can receive carbon from ectomycorrhizal fungi, and thus from trees. Our results reveal the absence of specificity in two mycoheterotrophic orchid-fungus associations in tropical regions, in contrast to most previous studies of mycoheterotrophic plants, which have been mainly focused on temperate orchids.     

Pseudo-Sanger sequencing: massively parallel production of long and near error-free reads using NGS technology

Background

Usually, next generation sequencing (NGS) technology has the property of ultra-high throughput but the read length is remarkably short compared to conventional Sanger sequencing. Paired-end NGS could computationally extend the read length but with a lot of practical inconvenience because of the inherent gaps. Now that Illumina paired-end sequencing has the ability of read both ends from 600 bp or even 800 bp DNA fragments, how to fill in the gaps between paired ends to produce accurate long reads is intriguing but challenging.

Results

We have developed a new technology, referred to as pseudo-Sanger (PS) sequencing. It tries to fill in the gaps between paired ends and could generate near error-free sequences equivalent to the conventional Sanger reads in length but with the high throughput of the Next Generation Sequencing. The major novelty of PS method lies on that the gap filling is based on local assembly of paired-end reads which have overlaps with at either end. Thus, we are able to fill in the gaps in repetitive genomic region correctly. The PS sequencing starts with short reads from NGS platforms, using a series of paired-end libraries of stepwise decreasing insert sizes. A computational method is introduced to transform these special paired-end reads into long and near error-free PS sequences, which correspond in length to those with the largest insert sizes. The PS construction has 3 advantages over untransformed reads: gap filling, error correction and heterozygote tolerance. Among the many applications of the PS construction is de novo genome assembly, which we tested in this study. Assembly of PS reads from a non-isogenic strain of Drosophila melanogaster yields an N50 contig of 190 kb, a 5 fold improvement over the existing de novo assembly methods and a 3 fold advantage over the assembly of long reads from 454 sequencing.

Conclusions

Our method generated near error-free long reads from NGS paired-end sequencing. We demonstrated that de novo assembly could benefit a lot from these Sanger-like reads. Besides, the characteristic of the long reads could be applied to such applications as structural variations detection and metagenomics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-711) contains supplementary material, which is available to authorized users.     

Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes

Background  

Molecular evolutionary studies share the common goal of elucidating historical relationships, and the common challenge of adequately sampling taxa and characters. Particularly at low taxonomic levels, recent divergence, rapid radiations, and conservative genome evolution yield limited sequence variation, and dense taxon sampling is often desirable. Recent advances in massively parallel sequencing make it possible to rapidly obtain large amounts of sequence data, and multiplexing makes extensive sampling of megabase sequences feasible. Is it possible to efficiently apply massively parallel sequencing to increase phylogenetic resolution at low taxonomic levels?     

16S rRNA sequences of uncultivated hot spring cyanobacterial mat inhabitants retrieved as randomly primed cDNA.

Cloning and analysis of cDNAs synthesized from rRNAs is one approach to assess the species composition of natural microbial communities. In some earlier attempts to synthesize cDNA from 16S rRNA (16S rcDNA) from the Octopus Spring cyanobacterial mat, a dominance of short 16S rcDNAs was observed, which appear to have originated only from certain organisms. Priming of cDNA synthesis from small ribosomal subunit RNA with random deoxyhexanucleotides can retrieve longer sequences, more suitable for phylogenetic analysis. Here we report the retrieval of 16S rRNA sequences from three formerly uncultured community members. One sequence type, which was retrieved three times from a total of five sequences analyzed, can be placed in the cyanobacterial phylum. A second sequence type is related to 16S rRNAs from green nonsulfur bacteria. The third sequence type may represent a novel phylogenetic type.     

16S rRNA sequences of uncultivated hot spring cyanobacterial mat inhabitants retrieved as randomly primed cDNA.

Cloning and analysis of cDNAs synthesized from rRNAs is one approach to assess the species composition of natural microbial communities. In some earlier attempts to synthesize cDNA from 16S rRNA (16S rcDNA) from the Octopus Spring cyanobacterial mat, a dominance of short 16S rcDNAs was observed, which appear to have originated only from certain organisms. Priming of cDNA synthesis from small ribosomal subunit RNA with random deoxyhexanucleotides can retrieve longer sequences, more suitable for phylogenetic analysis. Here we report the retrieval of 16S rRNA sequences from three formerly uncultured community members. One sequence type, which was retrieved three times from a total of five sequences analyzed, can be placed in the cyanobacterial phylum. A second sequence type is related to 16S rRNAs from green nonsulfur bacteria. The third sequence type may represent a novel phylogenetic type.     

Flow cytometry for enrichment and titration in massively parallel DNA sequencing

Massively parallel DNA sequencing is revolutionizing genomics research throughout the life sciences. However, the reagent costs and labor requirements in current sequencing protocols are still substantial, although improvements are continuously being made. Here, we demonstrate an effective alternative to existing sample titration protocols for the Roche/454 system using Fluorescence Activated Cell Sorting (FACS) technology to determine the optimal DNA-to-bead ratio prior to large-scale sequencing. Our method, which eliminates the need for the costly pilot sequencing of samples during titration is capable of rapidly providing accurate DNA-to-bead ratios that are not biased by the quantification and sedimentation steps included in current protocols. Moreover, we demonstrate that FACS sorting can be readily used to highly enrich fractions of beads carrying template DNA, with near total elimination of empty beads and no downstream sacrifice of DNA sequencing quality. Automated enrichment by FACS is a simple approach to obtain pure samples for bead-based sequencing systems, and offers an efficient, low-cost alternative to current enrichment protocols.     

Prenatal detection of aneuploidy and imbalanced chromosomal arrangements by massively parallel sequencing

Fetal chromosomal abnormalities are the most common reasons for invasive prenatal testing. Currently, G-band karyotyping and several molecular genetic methods have been established for diagnosis of chromosomal abnormalities. Although these testing methods are highly reliable, the major limitation remains restricted resolutions or can only achieve limited coverage on the human genome at one time. The massively parallel sequencing (MPS) technologies which can reach single base pair resolution allows detection of genome-wide intragenic deletions and duplication challenging karyotyping and microarrays as the tool for prenatal diagnosis. Here we reported a novel and robust MPS-based method to detect aneuploidy and imbalanced chromosomal arrangements in amniotic fluid (AF) samples. We sequenced 62 AF samples on Illumina GAIIx platform and with averagely 0.01× whole genome sequencing data we detected 13 samples with numerical chromosomal abnormalities by z-test. With up to 2× whole genome sequencing data we were able to detect microdeletion/microduplication (ranged from 1.4 Mb to 37.3 Mb of 5 samples from chorionic villus sampling (CVS) using SeqSeq algorithm. Our work demonstrated MPS is a robust and accurate approach to detect aneuploidy and imbalanced chromosomal arrangements in prenatal samples.