DePIE: Designing Primers for Protein Interaction Experiments

Several primer prediction and analysis programs have been developed for diverse applications. However, none of these existing programs can be directly used for the design of primers in protein interaction experiments, since proteins may have transmembrane domains (TMDs) and/or a signal peptide that must be excluded from experiments. Furthermore, it is frequently the case that a short restriction sequences must be added to each primer in order to clone PCR products into a given destination vectors for expression. DePIE, a web-based primer design tool, was developed to address these deficiencies. The program takes as input NCBI protein accession numbers and returns primer information including nucleotide sequences, thermodynamic melting temperature of the nucleotide sequences and the target positions. DePIE is implemented in JAVA, PERL and PHP and has proven to be very efficient in designing primers for our interaction experiments. DePIE services can be accessed at the web site: http://biocore.unl.edu/primer/primerPI.html.     

New Molecular Screening Tools for Analysis of Free-Living Diazotrophs in Soil

Free-living nitrogen-fixing prokaryotes (diazotrophs) are ubiquitous in soil and are phylogenetically and physiologically highly diverse. Molecular methods based on universal PCR detection of the nifH marker gene have been successfully applied to describe diazotroph populations in the environment. However, the use of highly degenerate primers and low-stringency amplification conditions render these methods prone to amplification bias, while less degenerate primer sets will not amplify all nifH genes. We have developed a fixed-primer-site approach with six PCR protocols using less degenerate to nondegenerate primer sets that all amplify the same nifH fragment as a previously published PCR protocol for universal amplification. These protocols target different groups of diazotrophs and allowed for direct comparison of the PCR products by use of restriction fragment length polymorphism fingerprinting. The new protocols were optimized on DNA from 14 reference strains and were subsequently tested with bulk DNA extracts from six soils. These analyses revealed that the new PCR primer sets amplified nifH sequences that were not detected by the universal primer set. Furthermore, they were better suited to distinguish between diazotroph populations in the different soils. Because the novel primer sets were not specific for monophyletic groups of diazotrophs, they do not serve as an identification tool; however, they proved powerful as fingerprinting tools for subsets of soil diazotroph communities.     

The development of highly transferable intron-spanning markers for temperate forage grasses

DNA markers have the potential to be a powerful tool for the molecular study and breeding of agronomic traits of temperate forage grasses, but some of these grasses have only limited sequence information available. We aimed to design highly transferable primer sets by using the abundant sequence information available for related crop species such as wheat and rice. The degree of similarity between the primer sequences of the wheat PCR-based landmark unique gene (PLUG) primer set and the corresponding sequences of rice orthologs was designated as the ??universality index?? (UI). We classified 359 PLUG primer pairs based on their UI values and found that primer pairs with high UI values showed higher successful amplification rates in Festuca and Lolium species than those with low UI values. Based on these results, we designed new primer sets, designated Conserved Three-prime-End Region (COTER) primers, with complete similarity to rice orthologs for eight bases at the 3?? end of each primer. COTER primer sets developed from both tall fescue and wheat showed high transferability in six temperate grasses (mean amplification rates of 95% for tall fescue primers; 79% for wheat primers). This strategy and primer information could also be useful for developing DNA markers for other grass species with little genome information available.     

Influence of Commonly Used Primer Systems on Automated Ribosomal Intergenic Spacer Analysis of Bacterial Communities in Environmental Samples

Due to the high diversity of bacteria in many ecosystems, their slow generation times, specific but mostly unknown nutrient requirements and syntrophic interactions, isolation based approaches in microbial ecology mostly fail to describe microbial community structure. Thus, cultivation independent techniques, which rely on directly extracted nucleic acids from the environment, are a well-used alternative. For example, bacterial automated ribosomal intergenic spacer analysis (B-ARISA) is one of the widely used methods for fingerprinting bacterial communities after PCR-based amplification of selected regions of the operon coding for rRNA genes using community DNA. However, B-ARISA alone does not provide any taxonomic information and the results may be severely biased in relation to the primer set selection. Furthermore, amplified DNA stemming from mitochondrial or chloroplast templates might strongly bias the obtained fingerprints. In this study, we determined the applicability of three different B-ARISA primer sets to the study of bacterial communities. The results from in silico analysis harnessing publicly available sequence databases showed that all three primer sets tested are specific to bacteria but only two primers sets assure high bacterial taxa coverage (1406f/23Sr and ITSF/ITSReub). Considering the study of bacteria in a plant interface, the primer set ITSF/ITSReub was found to amplify (in silico) sequences of some important crop species such as Sorghum bicolor and Zea mays. Bacterial genera and plant species potentially amplified by different primer sets are given. These data were confirmed when DNA extracted from soil and plant samples were analyzed. The presented information could be useful when interpreting existing B-ARISA results and planning B-ARISA experiments, especially when plant DNA can be expected.     

Design and Validation of Four New Primers for Next-Generation Sequencing To Target the 18S rRNA Genes of Gastrointestinal Ciliate Protozoa

Four new primers and one published primer were used to PCR amplify hypervariable regions within the protozoal 18S rRNA gene to determine which primer pair provided the best identification and statistical analysis. PCR amplicons of 394 to 498 bases were generated from three primer sets, sequenced using Roche 454 pyrosequencing with Titanium, and analyzed using the BLAST database (NCBI) and MOTHUR version 1.29. The protozoal diversity of rumen contents from moose in Alaska was assessed. In the present study, primer set 1, P-SSU-316F and GIC758R (amplicon of 482 bases), gave the best representation of diversity using BLAST classification, and the set amplified Entodinium simplex and Ostracodinium spp., which were not amplified by the other two primer sets. Primer set 2, GIC1080F and GIC1578R (amplicon of 498 bases), had similar BLAST results and a slightly higher percentage of sequences that were identified with a higher sequence identity. Primer sets 1 and 2 are recommended for use in ruminants. However, primer set 1 may be inadequate to determine protozoal diversity in nonruminants. The amplicons created by primer set 1 were indistinguishable for certain species within the genera Bandia, Blepharocorys, Polycosta, and Tetratoxum and between Hemiprorodon gymnoprosthium and Prorodonopsis coli, none of which are normally found in the rumen.     

Java web tools for PCR, in silico PCR, and oligonucleotide assembly and analysis

The polymerase chain reaction is fundamental to molecular biology and is the most important practical molecular technique for the research laboratory. We have developed and tested efficient tools for PCR primer and probe design, which also predict oligonucleotide properties based on experimental studies of PCR efficiency. The tools provide comprehensive facilities for designing primers for most PCR applications and their combinations, including standard, multiplex, long-distance, inverse, real-time, unique, group-specific, bisulphite modification assays, Overlap-Extension PCR Multi-Fragment Assembly, as well as a programme to design oligonucleotide sets for long sequence assembly by ligase chain reaction. The in silico PCR primer or probe search includes comprehensive analyses of individual primers and primer pairs. It calculates the melting temperature for standard and degenerate oligonucleotides including LNA and other modifications, provides analyses for a set of primers with prediction of oligonucleotide properties, dimer and G-quadruplex detection, linguistic complexity, and provides a dilution and resuspension calculator.     

Metabolic Primers for Detection of (Per)chlorate-Reducing Bacteria in the Environment and Phylogenetic Analysis of cld Gene Sequences

Natural attenuation of the environmental contaminant perchlorate is a cost-effective alternative to current removal methods. The success of natural perchlorate remediation is dependent on the presence and activity of dissimilatory (per)chlorate-reducing bacteria (DPRB) within a target site. To detect DPRB in the environment, two degenerate primer sets targeting the chlorite dismutase (cld) gene were developed and optimized. A nested PCR approach was used in conjunction with these primer sets to increase the sensitivity of the molecular detection method. Screening of environmental samples indicated that all products amplified by this method were cld gene sequences. These sequences were obtained from pristine sites as well as contaminated sites from which DPRB were isolated. More than one cld phylotype was also identified from some samples, indicating the presence of more than one DPRB strain at those sites. The use of these primer sets represents a direct and sensitive molecular method for the qualitative detection of (per)chlorate-reducing bacteria in the environment, thus offering another tool for monitoring natural attenuation. Sequences of cld genes isolated in the course of this project were also generated from various DPRB and provided the first opportunity for a phylogenetic treatment of this metabolic gene. Comparisons of the cld and 16S ribosomal DNA (rDNA) gene trees indicated that the cld gene does not track 16S rDNA phylogeny, further implicating the possible role of horizontal transfer in the evolution of (per)chlorate respiration.     

Flexgrepps--flexible greedy peptide pool search: computation of near-optimal sets of degenerate polypeptides for antigenic screening

Although synthesizing and utilizing individual peptides and DNA primers has become relatively inexpensive, massively parallel probing and next-generation sequencing approaches have dramatically increased the number of molecules that can be subjected to screening; this, in turn, requires vast numbers of peptides and therefore results in significant expenses. To alleviate this issue, pools of related molecules are often used to downselect prior to testing individual sequences. A computational selection process to create pools of related sequences at large scale has not been reported for peptides. In the case of PCR primers, there have been successful attempts to address this problem by designing degenerate primers that can be produced at the same cost as conventional, unique primers and then be used to amplify several different genomic regions. We present an algorithm, "FlexGrePPS" (Flexible Greedy Peptide Pool Search), that can create a near-optimal set of peptide pools. This approach is also applicable to nucleotide sequences and outperforms most DNA primer selection programs. For the proteomic compression with FlexGrePPS, the main body of our work presented here, we demonstrate the feasibility of the computation of an exhaustive cover of pathogenic proteomes with degenerate peptides that lend themselves to antigenic screening. Furthermore, we present preliminary data that demonstrate the experimental utility of highly degenerate peptides for antigenic screening. FlexGrePPS provides a near-optimal solution for proteomic compression and there are no programs available for comparison. We also demonstrate computational performance of our GreedyPrime implementation, which is a modified version of FlexGrePPS applicable to the design of degenerate primers and is comparable to existing programs for the design of degenerate primers. Specifically, we focus on the comparisons with PAMPS and DPS-DIP, software tools that have recently been shown to be superior to other methods. FlexGrePPS forms the foundation of a novel antigenic screening methodology that is based on the representation of an entire proteome by near-optimal degenerate peptide pools. Our preliminary wet lab data indicate that the approach will likely prove successful in comprehensive wet lab studies, and hence will dramatically reduce the expenses for antigenic screening and make whole proteome screening feasible. Although FlexGrePPS was designed for computational performance in order to handle vast data sets, there is the very surprising finding that even for small data sets the primer design version of FlexGrePPS, GreedyPrime, offers similar or even superior results for MP-DPD and most MDPD instances when compared to existing methods; despite the much longer run times, other approaches did not fare significantly better in reducing the original data sets to degenerate primers. The FlexGrePPS and GreedyPrime programs are available at no charge under the GNU LGPL license at http://sourceforge.net/projects/flexgrepps/.     

The translation of recombinant proteins in E. coli can be improved by in silico generating and screening random libraries of a −70/+96 mRNA region with respect to the translation initiation codon

Recombinant protein translation in Escherichia coli may be limited by stable (i.e. low free energy) secondary structures in the mRNA translation initiation region. To circumvent this issue, we have set-up a computer tool called ‘ExEnSo’ (Expression Enhancer Software) that generates a random library of 8192 sequences, calculates the free energy of secondary structures of each sequence in the −70/+96 region (base 1 is the translation initiation codon), and then selects the sequence having the highest free energy. The software uses this ‘optimized’ sequence to create a 5′ primer that can be used in PCR experiments to amplify the coding sequence of interest prior to sub-cloning into a prokaryotic expression vector. In this article, we report how ExEnSo was set-up and the results obtained with nine coding sequences with low expression levels in E. coli. The free energy of the −70/+96 region of all these coding sequences was increased compared to the non-optimized sequences. Moreover, the protein expression of eight out of nine of these coding sequences was increased in E. coli, indicating a good correlation between in silico and in vivo results. ExEnSo is available as a free online tool.     

Editor's choice: AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding—but rarely from non-coding—DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.     

pr2-primers: An 18S rRNA primer database for protists

Metabarcoding of microbial eukaryotes (collectively known as protists) has developed tremendously in the last decade, almost solely relying on the 18S rRNA gene. As microbial eukaryotes are extremely diverse, many primers and primer pairs have been developed. To cover a relevant and representative fraction of the protist community in a given study system, an informed primer choice is necessary, as no primer pair can target all protists equally well. As such, a smart primer choice is very difficult even for experts and there are very few online resources available to list existing primers. We built a database listing 285 primers and 83 unique primer pairs that have been used for eukaryotic 18S rRNA gene metabarcoding. In silico performance of primer pairs was tested against two sequence databases: PR² version 4.12.0 for eukaryotes and a subset of silva version 132 for bacteria and archaea. We developed an R -based web application enabling browsing of the database, visualization of the taxonomic distribution of the amplified sequences with the number of mismatches, and testing any user-defined primer or primer set ( https://app.pr2-primers.org ). Taxonomic specificity of primer pairs, amplicon size and location of mismatches can also be determined. We identified universal primer sets that matched the largest number of sequences and analysed the specificity of some primer sets designed to target certain groups. This tool enables guided primer choices that will help a wide range of researchers to include protists as part of their investigations.     

Novel PCR Primers for the Archaeal Phylum Thaumarchaeota Designed Based on the Comparative Analysis of 16S rRNA Gene Sequences

Based on comparative phylogenetic analysis of 16S rRNA gene sequences deposited in an RDP database, we constructed a local database of thaumarchaeotal 16S rRNA gene sequences and developed a novel PCR primer specific for the archaeal phylum Thaumarchaeota. Among 9,727 quality-filtered (chimeral-checked, size >1.2 kb) archaeal sequences downloaded from the RDP database, 1,549 thaumarchaeotal sequences were identified and included in our local database. In our study, Thaumarchaeota included archaeal groups MG-I, SAGMCG-I, SCG, FSCG, RC, and HWCG-III, forming a monophyletic group in the phylogenetic tree. Cluster analysis revealed 114 phylotypes for Thaumarchaeota. The majority of the phylotypes (66.7%) belonged to the MG-I and SCG, which together contained most (93.9%) of the thaumarchaeotal sequences in our local database. A phylum-directed primer was designed from a consensus sequence of the phylotype sequences, and the primer’s specificity was evaluated for coverage and tolerance both in silico and empirically. The phylum-directed primer, designated THAUM-494, showed >90% coverage for Thaumarchaeota and <1% tolerance to non-target taxa, indicating high specificity. To validate this result experimentally, PCRs were performed with THAUM-494 in combination with a universal archaeal primer (ARC917R or 1017FAR) and DNAs from five environmental samples to construct clone libraries. THAUM-494 showed a satisfactory specificity in empirical studies, as expected from the in silico results. Phylogenetic analysis of 859 cloned sequences obtained from 10 clone libraries revealed that >95% of the amplified sequences belonged to Thaumarchaeota. The most frequently sampled thaumarchaeotal subgroups in our samples were SCG, MG-I, and SAGMCG-I. To our knowledge, THAUM-494 is the first phylum-level primer for Thaumarchaeota. Furthermore, the high coverage and low tolerance of THAUM-494 will make it a potentially valuable tool in understanding the phylogenetic diversity and ecological niche of Thaumarchaeota.     

Validation and application of a PCR primer set to quantify fungal communities in the soil environment by real-time quantitative PCR

Fungi constitute an important group in soil biological diversity and functioning. However, characterization and knowledge of fungal communities is hampered because few primer sets are available to quantify fungal abundance by real-time quantitative PCR (real-time Q-PCR). The aim in this study was to quantify fungal abundance in soils by incorporating, into a real-time Q-PCR using the SYBRGreen® method, a primer set already used to study the genetic structure of soil fungal communities. To satisfy the real-time Q-PCR requirements to enhance the accuracy and reproducibility of the detection technique, this study focused on the 18S rRNA gene conserved regions. These regions are little affected by length polymorphism and may provide sufficiently small targets, a crucial criterion for enhancing accuracy and reproducibility of the detection technique. An in silico analysis of 33 primer sets targeting the 18S rRNA gene was performed to select the primer set with the best potential for real-time Q-PCR: short amplicon length; good fungal specificity and coverage. The best consensus between specificity, coverage and amplicon length among the 33 sets tested was the primer set FR1 / FF390. This in silico analysis of the specificity of FR1 / FF390 also provided additional information to the previously published analysis on this primer set. The specificity of the primer set FR1 / FF390 for Fungi was validated in vitro by cloning - sequencing the amplicons obtained from a real time Q-PCR assay performed on five independent soil samples. This assay was also used to evaluate the sensitivity and reproducibility of the method. Finally, fungal abundance in samples from 24 soils with contrasting physico-chemical and environmental characteristics was examined and ranked to determine the importance of soil texture, organic carbon content, C∶N ratio and land use in determining fungal abundance in soils.     

Selection and validation of reference genes for real-time RT-PCR studies in the non-model species Delomys sublineatus, an endemic Brazilian rodent

Quantitative real-time RT-PCR (qRT-PCR) is a sensitive technique for gene expression analysis. A critical factor for creating reliable data in relative quantification is the normalization of the expression data of genes of interest. Therefore the needed normalization factor is calculated out of the expression data of co-amplified genes that are stable expressed in the certain sample material, the so-called reference genes. In this study, we demonstrate the important process of validating potential reference genes using a non-model species. As there are almost no sequences known of the Pallid Atlantic Forest Rat (Delomys sublineatus), a rodent used as indicator species in conservation studies of the endangered Brazilian rainforest, suitable primer sets are more problematic to find than in model species. Out of nine tested primer sets designed for the fully sequenced Mus musculus, five could be used for the establishment of a proper running SYBR-Green assay and validation of their constant expression. qRT-PCR results of 12 cDNAs of Delomys livers were analyzed with three different validation software programs: BestKeeper, NormFinder and geNorm. Our approach showed that out of the five (Sdha, Canx, Pgk1, Actb and Actg1) potential reference genes, the first four should be used for accurate normalization in further relative quantification analyses. Transferring data from close-by model organisms makes high sensitive real-time RT-PCR applicable even to free-ranging non-model organisms. Our approach might be suitable for other non-model organisms.     

Novel primers for 16S rRNA-based archaeal community analyses in environmental samples

Next generation sequencing technologies for in depth analyses of complex microbial communities rely on rational primer design based on up-to-date reference databases. Most of the 16S rRNA-gene based analyses of environmental Archaea community composition use PCR primers developed from small data sets several years ago, making an update long overdue. Here we present a new set of archaeal primers targeting the 16S rRNA gene designed from 8500 aligned archaeal sequences in the SILVA database. The primers 340F-1000R showed a high archaeal specificity (< 1% bacteria amplification) covering 93 and 97% of available sequences for Crenarchaeota and Euryarchaeota respectively. In silico tests of the primers revealed at least 38% higher coverage for Archaea compared to other commonly used primers. Empirical tests with clone libraries confirmed the high specificity of the primer pair to Archaea in three biomes: surface waters in the Arctic Ocean, the pelagic zone of a temperate lake and a methanogenic bioreactor. The clone libraries featured both Euryarchaeota and Crenarchaeota in variable proportions and revealed dramatic differences in the archaeal community composition and minimal phylogenetic overlap between samples.     

BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes

Bisulfite genomic sequencing is the most widely used technique to analyze the 5-methylation of cytosines, the prevalent covalent DNA modification in mammals. The process is based on the selective transformation of unmethylated cytosines to uridines. Then, the investigated genomic regions are PCR amplified, subcloned and sequenced. During sequencing, the initially unmethylated cytosines are detected as thymines. The efficacy of bisulfite PCR is generally low; mispriming and non-specific amplification often occurs due to the T richness of the target sequences. In order to ameliorate the efficiency of PCR, we developed a new primer-design software called BiSearch, available on the World Wide Web. It has the unique property of analyzing the primer pairs for mispriming sites on the bisulfite-treated genome and determines potential non-specific amplification products with a new search algorithm. The options of primer-design and analysis for mispriming sites can be used sequentially or separately, both on bisulfite-treated and untreated sequences. In silico and in vitro tests of the software suggest that new PCR strategies may increase the efficiency of the amplification.     

MiCA: A Web-Based Tool for the Analysis of Microbial Communities Based on Terminal-Restriction Fragment Length Polymorphisms of 16S and 18S rRNA Genes

A web-based resource, Microbial Community Analysis (MiCA), has been developed to facilitate studies on microbial community ecology that use analyses of terminal-restriction fragment length polymorphisms (T-RFLP) of 16S and 18S rRNA genes. MiCA provides an intuitive web interface to access two specialized programs and a specially formatted database of 16S ribosomal RNA sequences. The first program performs virtual polymerase chain reaction (PCR) amplification of rRNA genes and restriction of the amplicons using primer sequences and restriction enzymes chosen by the user. This program, in silico PCR and Restriction (ISPaR), uses a binary encoding of DNA sequences to rapidly scan large numbers of sequences in databases searching for primer annealing and restriction sites while permitting the user to specify the number of mismatches in primer sequences. ISPaR supports multiple digests with up to three enzymes. The number of base pairs between the 5′ and 3′ primers and the proximal restriction sites can be reported, printed, or exported in various formats. The second program, APLAUS, infers a plausible community structure(s) based on T-RFLP data supplied by a user. APLAUS estimates the relative abundances of populations and reports a listing of phylotypes that are consistent with the empirical data. MiCA is accessible at .     

PlasmaDNA: a free, cross-platform plasmid manipulation program for molecular biology laboratories

Background

Most molecular biology experiments, and the techniques associated with this field of study, involve a great deal of engineering in the form of molecular cloning. Like all forms of engineering, perfect information about the starting material is crucial for successful completion of design and strategies.

Results

We have generated a program that allows complete in silico simulation of the cloning experiment. Starting with a primary DNA sequence, PlasmaDNA looks for restriction sites, open reading frames, primer annealing sequences, and various common domains. The databases are easily expandable by the user to fit his most common cloning needs. PlasmaDNA can manage and graphically represent multiple sequences at the same time, and keeps in memory the overhangs at the end of the sequences if any. This means that it is possible to virtually digest fragments, to add the digestion products to the project, and to ligate together fragments with compatible ends to generate the new sequences. Polymerase Chain Reaction (PCR) fragments can also be virtually generated using the primer database, automatically adding to the fragments any 5' extra sequences present in the primers.

Conclusion

PlasmaDNA is a program available both on Windows and Apple operating systems, designed to facilitate molecular cloning experiments by building a visual map of the DNA. It then allows the complete planning and simulation of the cloning experiment. It also automatically updates the new sequences generated in the process, which is an important help in practice. The capacity to maintain multiple sequences in the same file can also be used to archive the various steps and strategies involved in the cloning of each construct. The program is freely available for download without charge or restriction.     

Group-specific PCR primers for the phylum Acidobacteria designed based on the comparative analysis of 16S rRNA gene sequences

We performed a comprehensive phylogenetic analysis of the phylum Acidobacteria and developed novel, group-specific PCR primers for Acidobacteria and its class-level subgroups. Acidobacterial 16S rRNA gene sequences deposited in the RDP database were used to construct a local database then subsequently analyzed. A total of 556 phylotypes were observed and the majority of the phylotypes belonged to five major subgroups (subgroups 1, 2, 3, 4, and 6), which comprised > 80% of the acidobacterial sequences in the RDP database. Phylum-specific and subgroup-specific primers were designed from the consensus sequences of the phylotype sequences, and the specificities of the designed primers were evaluated both in silico and empirically for coverage and tolerance. The phylum-specific primer ACIDO, which was designed in this study, showed increased coverage for Acidobacteria, as compared to the previous phylum-specific primer 31F. However, the tolerance of the primer ACIDO for non-target sequences was slightly higher than that of the primer 31F. We also developed subgroup-specific PCR primers for the major subgroups of Acidobacteria, except for subgroup 4. Subgroup-specific primers S1, S2, and S3, which targeted subgroups 1, 2, and 3, respectively, showed high coverage for their target subgroups and low tolerance for non-target sequences. However, the primer S6 targeting subgroup 6 showed a lower specificity in its empirical evaluation than expected from the in silico results. The subgroup-specific primers, as well as the phylum-specific primer designed in this study, will be valuable tools in understanding the phylogenetic diversity and ecological niche of the phylum Acidobacteria and its subgroups.     

Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—a baiting and iterative mapping approach

We present an in silico approach for the reconstruction of complete mitochondrial genomes of non-model organisms directly from next-generation sequencing (NGS) data—mitochondrial baiting and iterative mapping (MITObim). The method is straightforward even if only (i) distantly related mitochondrial genomes or (ii) mitochondrial barcode sequences are available as starting-reference sequences or seeds, respectively. We demonstrate the efficiency of the approach in case studies using real NGS data sets of the two monogenean ectoparasites species Gyrodactylus thymalli and Gyrodactylus derjavinoides including their respective teleost hosts European grayling (Thymallus thymallus) and Rainbow trout (Oncorhynchus mykiss). MITObim appeared superior to existing tools in terms of accuracy, runtime and memory requirements and fully automatically recovered mitochondrial genomes exceeding 99.5% accuracy from total genomic DNA derived NGS data sets in <24 h using a standard desktop computer. The approach overcomes the limitations of traditional strategies for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information at hand and represents a fast and highly efficient in silico alternative to laborious conventional strategies relying on initial long-range PCR. We furthermore demonstrate the applicability of MITObim for metagenomic/pooled data sets using simulated data. MITObim is an easy to use tool even for biologists with modest bioinformatics experience. The software is made available as open source pipeline under the MIT license at https://github.com/chrishah/MITObim.