Don't put all your eggs in one basket: a cost‐effective and powerful method to optimize primer choice for rRNA environmental community analyses using the Fluidigm Access Array

With the increasing democratization of high‐throughput sequencing (HTS) technologies, along with the concomitant increase in sequence yield per dollar, many researchers are exploring HTS for microbial community ecology. Many elements of experimental design can drastically affect the final observed community structure, notably the choice of primers for amplification prior to sequencing. Some targeted microbes can fail to amplify due to primer‐targeted sequence divergence and be omitted from obtained sequences, leading to differences among primer pairs in the sequenced organisms even when targeting the same community. This potential source of taxonomic bias in HTS makes it prudent to investigate how primer choice will affect the sequenced community prior to investing in a costly community‐wide sequencing effort. Here, we use Fluidigm's microfluidic Access Arrays (IFC) followed by Illumina^® MiSeq Nano sequencing on a culture‐derived local mock community to demonstrate how this approach allows for a low‐cost combinatorial investigation of primer pairs and experimental samples (up to 48 primer pairs and 48 samples) to determine the most effective primers that maximize obtained communities whilst minimizing taxonomic biases.     

An improved method for utilizing high‐throughput amplicon sequencing to determine the diets of insectivorous animals

DNA analysis of predator faeces using high‐throughput amplicon sequencing (HTS) enhances our understanding of predator–prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (COI). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair “ZBJ” to results using the novel primer pair “ANML.” To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single‐copy (cloned) COI sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre‐ and post‐PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%), whereas ZBJ amplified 24–40 of 59 taxa (41%–68%). Furthermore, in an HTS comparison of field‐collected samples, the ANML primers detected nearly fourfold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.     

Amplicon –Based Metagenomic Analysis of Mixed Fungal Samples Using Proton Release Amplicon Sequencing

Next generation sequencing technology has revolutionised microbiology by allowing concurrent analysis of whole microbial communities. Here we developed and verified similar methods for the analysis of fungal communities using a proton release sequencing platform with the ability to sequence reads of up to 400 bp in length at significant depth. This read length permits the sequencing of amplicons from commonly used fungal identification regions and thereby taxonomic classification. Using the 400 bp sequencing capability, we have sequenced amplicons from the ITS1, ITS2 and LSU fungal regions to a depth of approximately 700,000 raw reads per sample. Representative operational taxonomic units (OTUs) were chosen by the USEARCH algorithm, and identified taxonomically through nucleotide blast (BLASTn). Combination of this sequencing technology with the bioinformatics pipeline allowed species recognition in two controlled fungal spore populations containing members of known identity and concentration. Each species included within the two controlled populations was found to correspond to a representative OTU, and these OTUs were found to be highly accurate representations of true biological sequences. However, the absolute number of reads attributed to each OTU differed among species. The majority of species were represented by an OTU derived from all three genomic regions although in some cases, species were only represented in two of the regions due to the absence of conserved primer binding sites or due to sequence composition. It is apparent from our data that proton release sequencing technologies can deliver a qualitative assessment of the fungal members comprising a sample. The fact that some fungi cannot be amplified by specific “conserved” primer pairs confirms our recommendation that a multi-region approach be taken for other amplicon-based metagenomic studies.     

Modular tagging of amplicons using a single PCR for high‐throughput sequencing

High‐throughput sequencing (HTS) of PCR amplicons is becoming the method of choice to sequence one or several targeted loci for phylogenetic and DNA barcoding studies. Although the development of HTS has allowed rapid generation of massive amounts of DNA sequence data, preparing amplicons for HTS remains a rate‐limiting step. For example, HTS platforms require platform‐specific adapter sequences to be present at the 5′ and 3′ end of the DNA fragment to be sequenced. In addition, short multiplex identifier (MID) tags are typically added to allow multiple samples to be pooled in a single HTS run. Existing methods to incorporate HTS adapters and MID tags into PCR amplicons are either inefficient, requiring multiple enzymatic reactions and clean‐up steps, or costly when applied to multiple samples or loci (fusion primers). We describe a method to amplify a target locus and add HTS adapters and MID tags via a linker sequence using a single PCR. We demonstrate our approach by generating reference sequence data for two mitochondrial loci (COI and 16S) for a diverse suite of insect taxa. Our approach provides a flexible, cost‐effective and efficient method to prepare amplicons for HTS.     

Multiplex PCR with the blunt hairpin primers for next generation sequencing

The multiplex PCR is one of the important methods to enrich the target DNAs for next generation sequencing. The non-specific amplification and interaction between the primers are the pivotal challenges of multiplex PCR. Here, we introduce the novel blunt hairpin primers for effective reducing the primer dimers and mispriming events. We also used a pair of auxiliary primers to enhance PCR efficiency. We simultaneously amplified 89 target regions from 44 samples and sequenced all amplicons on ion torrent PGM platform. Among all the filtrated amplicons (3438 different amplicons), 99.7, 97.6, 90.1 and 72.8% had sequencing depths fell within 200, 100, 50 and 25-fold range. The sequencing depth variations among all the samples were less than 27-fold. We also amplified multiplex regions with blunt hairpin, stick hairpin and normal linear primers, and the blunt hairpin primers could significantly reduce the amount of primer dimers and unspecific products.These results show that multiplex PCR with the blunt hairpin primers is a flexible, specific and economical target-region captured approach for the next generation sequencing.     

Target region amplification polymorphism: A novel marker technique for plant genotyping

The advent of large-scale DNA sequencing technology has generated a tremendous amount of sequence information for many important organisms. We have developed a rapid and efficient PCR-based technique, which uses bioinformatics tools and expressed sequence tag (EST) database information to generate polymorphic markers around targeted candidate gene sequences. This target region amplification polymorphism (TRAP) technique uses 2 primers of 18 nucleotides to generate markers. One of the primers, the fixed primer, is designed from the targeted EST sequence in the database; the second primer, the arbitrary primer, is an arbitrary sequence with either an AT-or GC-rich core to anneal with an intron or exon, respectively. PCR amplification is run for the first 5 cycles with an annealing temperature of 35°C, followed by 35 cycles with an annealing temperature of 50°C. For different plant species, each PCR reaction can generate as many as 50 scorable fragments with sizes ranging from 50–900 bp when separated on a 6.5% polyacrylamide sequencing gel. The TRAP technique should be useful in genotyping germplasm collections and in tagging genes governing desirable agronomic traits of crop plants.     

Thermodynamically balanced inside-out (TBIO) PCR-based gene synthesis: a novel method of primer design for high-fidelity assembly of longer gene sequences

A novel thermodynamically-balanced inside-out (TBIO) method of primer design was developed and compared with a thermodynamically-balanced conventional (TBC) method of primer design for PCR-based gene synthesis of codon-optimized gene sequences for the human protein kinase B-2 (PKB2; 1494 bp), p70 ribosomal S6 subunit protein kinase-1 (S6K1; 1622 bp) and phosphoinositide-dependent protein kinase-1 (PDK1; 1712 bp). Each of the 60mer TBIO primers coded for identical nucleotide regions that the 60mer TBC primers covered, except that half of the TBIO primers were reverse complement sequences. In addition, the TBIO and TBC primers contained identical regions of temperature- optimized primer overlaps. The TBC method was optimized to generate sequential overlapping fragments (~0.4–0.5 kb) for each of the gene sequences, and simultaneous and sequential combinations of overlapping fragments were tested for their ability to be assembled under an array of PCR conditions. However, no fully synthesized gene sequences could be obtained by this approach. In contrast, the TBIO method generated an initial central fragment (~0.4–0.5 kb), which could be gel purified and used for further inside-out bidirectional elongation by additional increments of 0.4–0.5 kb. By using the newly developed TBIO method of PCR-based gene synthesis, error-free synthetic genes for the human protein kinases PKB2, S6K1 and PDK1 were obtained with little or no corrective mutagenesis.     

Development of a genotype‐by‐sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.     

High-throughput siRNA-based functional target validation

The drug discovery process pursued by major pharmaceutical companies for many years starts with target identification followed by high-throughput screening (HTS) with the goal of identifying lead compounds. To accomplish this goal, significant resources are invested into automation of the screening process or HTS. Robotic systems capable of handling thousands of data points per day are implemented across the pharmaceutical sector. Many of these systems are amenable to handling cell-based screening protocols as well. On the other hand, as companies strive to develop innovative products based on novel mechanisms of action(s), one of the current bottlenecks of the industry is the target validation process. Traditionally, bioinformatics and HTS groups operate separately at different stages of the drug discovery process. The authors describe the convergence and integration of HTS and bioinformatics to perform high-throughput target functional identification and validation. As an example of this approach, they initiated a project with a functional cell-based screen for a biological process of interest using libraries of small interfering RNA (siRNA) molecules. In this protocol, siRNAs function as potent gene-specific inhibitors. siRNA-mediated knockdown of the target genes is confirmed by TaqMan analysis, and genes with impacts on biological functions of interest are selected for further analysis. Once the genes are confirmed and further validated, they may be used for HTS to yield lead compounds.     

Expressed sequence tags for the chicken genome from a normalized 10-day-old White Leghorn whole embryo cDNA library: 1. DNA sequence characterization and linkage analysis

Expressed sequence tags (ESTs) provide a rapid and reliable method for gene discovery as well as a resource for the large-scale analysis of gene expression of known and unknown genes. Here we describe a normalized cDNA library developed from a 10-day-old White Leghorn chicken whole embryo. The utility of the library was evaluated by partial sequencing of 99 randomly selected insert-containing clones and the analysis of EST-targeted genomic regions for single nucleotide polymorphisms (SNPs) in the East Lansing chicken reference DNA mapping panel. Using stringent match criteria of percent identity of 80 or higher across a length of 50 or more bases, 46 ESTs matched database sequences including previously reported Gallus gallus genes. Thirty-seven of the 50 primer pairs developed from 50 unique ESTs amplified a single fragment. The size of the 37 amplicons ranged from 276 to 693 bp for a total of 17,508 and an average of 473. About 70% of the SNPs detected were either G-->A or C-->T transition. The number of SNPs detected within the amplicons from EST-targeted genomic regions ranged from 0 to 4 for a total of 65 and a frequency of about 1 every 470 bases. About 35% of the amplicons contained only 1 SNP, while 19% had 4 SNPs. Using the SNPs that were informative in the East Lansing reference panel, 17 ESTs were mapped on the East Lansing chicken genetic map. The ESTs described, as well as the nucleotide variants identified within the EST-targeted genomic regions, represent significant resources for genome analysis in the chicken.     

Molecular Profiling of the Clostridium leptum Subgroup in Human Fecal Microflora by PCR-Denaturing Gradient Gel Electrophoresis and Clone Library Analysis

A group-specific PCR-based denaturing gradient gel electrophoresis (DGGE) method was developed and combined with group-specific clone library analysis to investigate the diversity of the Clostridium leptum subgroup in human feces. PCR products (length, 239 bp) were amplified using C. leptum cluster-specific primers and were well separated by DGGE. The DGGE patterns of fecal amplicons from 11 human individuals revealed host-specific profiles; the patterns for fecal samples collected from a child for 3 years demonstrated the structural succession of the population in the first 2 years and its stability in the third year. A clone library was constructed with 100 clones consisting of 1,143-bp inserts of 16S rRNA gene fragments that were amplified from one adult fecal DNA with one forward universal bacterial primer and one reverse group-specific primer. Eighty-six of the clones produced the 239-bp C. leptum cluster-specific amplicons, and the remaining 14 clones did not produce these amplicons but still phylogenetically belong to the subgroup. Sixty-four percent of the clones were related to Faecalibacterium prausnitzii (similarity, 97 to 99%), 6% were related to Subdoligranulum variabile (similarity, ~99%), 2% were related to butyrate-producing bacterium A2-207 (similarity, 99%), and 28% were not identified at the species level. The identities of most bands in the DGGE profiles for the same adult were determined by comigration analysis with the 86 clones that harbored the 239-bp group-specific fragments. Our results suggest that DGGE combined with clone library analysis is an effective technique for monitoring and analyzing the composition of this important population in the human gut flora.     

Capture of assay template by multiplex PCR of long amplicons for genotyping SNPs and InDels with MALDI-TOF mass spectrometry

Mis-priming associated with uncharacterised single nucleotide polymorphisms (SNPs) may lead to failure of PCR for genotyping. This is particularly troublesome in high-throughput SNP genotyping applications relying on multiplex PCR (2–40-plex) generating many short amplicons (80–120 bp) of similar size, an approach best suited for whole genome scans. However, if the target SNPs are clustered within a few target genes one option to ameliorate this is to increase the amplicon length, effectively reducing the potential for primer/template interactions and mis-priming. We tested this approach in a diverse population of 372 Eucalyptus pilularis individuals (π = 8.11 × 10⁻³, H _e = 0.75) using a modified Sequenom iPLEX gold assay. Four candidate genes (MYB1, MYB2, CAD and CCR) were amplified in a single long range multiplex capture PCR generating 6 long amplicons ranging in size from 907 to 2,225 bp. This contrasts with the standard approach which would have required the amplification of 98 short amplicons in 4 multiplex reactions. These 6 long amplicons provided the assay template for 98 assays (87 SNP and 11 InDel) within the 4 candidate genes. Reaction results indicated that longer amplicons could provide a suitable template for genotyping assays, with 90.8% of assays functional and 84.3% of assays suitable for downstream analysis. Additional advantages of this approach were the capacity for troubleshooting using gel electrophoresis and savings of 94% in capture primer synthesis costs. This approach will have the greatest relevance for candidate gene approaches for association testing in uncharacterised populations of organisms with high sequence diversity.     

Performance of amplicon and shotgun sequencing for accurate biomass estimation in invertebrate community samples

New applications of DNA and RNA sequencing are expanding the field of biodiversity discovery and ecological monitoring, yet questions remain regarding precision and efficiency. Due to primer bias, the ability of metabarcoding to accurately depict biomass of different taxa from bulk communities remains unclear, while PCR‐free whole mitochondrial genome (mitogenome) sequencing may provide a more reliable alternative. Here, we used a set of documented mock communities comprising 13 species of freshwater macroinvertebrates of estimated individual biomass, to compare the detection efficiency of COI metabarcoding (three different amplicons) and shotgun mitogenome sequencing. Additionally, we used individual COI barcoding and de novo mitochondrial genome sequencing, to provide reference sequences for OTU assignment and metagenome mapping (mitogenome skimming), respectively. We found that, even though both methods occasionally failed to recover very low abundance species, metabarcoding was less consistent, by failing to recover some species with higher abundances, probably due to primer bias. Shotgun sequencing results provided highly significant correlations between read number and biomass in all but one species. Conversely, the read–biomass relationships obtained from metabarcoding varied across amplicons. Specifically, we found significant relationships for eight of 13 (amplicons B1FR‐450 bp, FF130R‐130 bp) or four of 13 (amplicon FFFR, 658 bp) species. Combining the results of all three COI amplicons (multiamplicon approach) improved the read–biomass correlations for some of the species. Overall, mitogenomic sequencing yielded more informative predictions of biomass content from bulk macroinvertebrate communities than metabarcoding. However, for large‐scale ecological studies, metabarcoding currently remains the most commonly used approach for diversity assessment.     

ARDRA对植物根瘤内共生放线菌Frankia多样性的研究

应用原核生物16S rDNA特异性引物rD1和fD1,通过ARDRA(amplified ribosomal DNA restrictionanalysis)法直接扩增自中国云南、东北地区赤杨属3种植物和沙棘属1种植物根瘤内Frankia DNA,得到一长约1500bp的扩增产物.选用两种内切酶HaeⅢ、AfaⅠ联合对扩增产物进行酶切,得到稳定的酶切图谱.将所测8个感染赤杨的Frankia样本区分为3个不同的组,所测3个感染沙棘的Frankia样本区分为3个不同的组.显示根瘤内Frankia存在丰富的遗传多样性.     

A targeted capture approach to generating reference sequence databases for chloroplast gene regions

Metabarcoding has improved the way we understand plants within our environment, from their ecology and conservation to invasive species management. The notion of identifying plant taxa within environmental samples relies on the ability to match unknown sequences to known reference libraries. Without comprehensive reference databases, species can go undetected or be incorrectly assigned, leading to false‐positive and false‐negative detections. To improve our ability to generate reference sequence databases, we developed a targeted capture approach using the OZBaits_CP V1.0 set, designed to capture chloroplast gene regions across the entirety of flowering plant diversity. We focused on generating a reference database for coastal temperate plant species given the lack of reference sequences for these taxa. Our approach was successful across all specimens with a target gene recovery rate of 92%, which was achieved in a single assay (i.e., samples were pooled), thus making this approach much faster and more efficient than standard barcoding. Further testing of this database highlighted 80% of all samples could be discriminated to family level across all gene regions with some genes achieving greater resolution than others—which was also dependent on the taxon of interest. Thus, we demonstrate the importance of generating reference sequences across multiple chloroplast gene regions as no single loci are sufficient to discriminate across all plant groups. The targeted capture approach outlined in this study provides a way forward to achieve this.     

Screening of virus specific primers against the Cassava mosaic virus

Abstract

Five different Indian cassava mosaic virus (ICMV) specific primers were used to screen the virus from CMD affected samples collected from the different parts of Tamil Nadu. Out of five specific primers, three were designed to amplify the specific viral genes of ICMV and two were used for detection of ICMV. All primers amplified specific regions of the virus in all samples. The specific primer for amplification of coat protein gene of ICMV amplified 800 bp of coat protein gene from both ICMV and Sri Lankan cassava mosaic virus (SLCMV) infected samples invariably. The specific primer for amplifying movement protein (MP) gene amplified about 900 bp of movement protein gene from all CMD infected cassava samples. Likewise, 800 bp of nuclear shuttle protein (NSP) gene was amplified from all the samples. The primer ICMV A amplified 700 bp of PCR product from mosaic diseased cassava samples. A 300 bp product from DNA A of the virus amplified in all samples using the primer ICMV A₁.     

In vivo phage display: identification of organ-specific peptides using deep sequencing and differential profiling across tissues

In vivo phage display is widely used for identification of organ- or disease-specific homing peptides. However, the current in vivo phage biopanning approaches fail to assess biodistribution of specific peptide phages across tissues during the screen, thus necessitating laborious and time-consuming post-screening validation studies on individual peptide phages. Here, we adopted bioinformatics tools used for RNA sequencing for analysis of high-throughput sequencing (HTS) data to estimate the representation of individual peptides during biopanning in vivo. The data from in vivo phage screen were analyzed using differential binding—relative representation of each peptide in the target organ versus in a panel of control organs. Application of this approach in a model study using low-diversity peptide T7 phage library with spiked-in brain homing phage demonstrated brain-specific differential binding of brain homing phage and resulted in identification of novel lung- and brain-specific homing peptides. Our study provides a broadly applicable approach to streamline in vivo peptide phage biopanning and to increase its reproducibility and success rate.     

MARA: a novel approach for highly multiplexed locus-specific SNP genotyping using high-density DNA oligonucleotide arrays

We have developed a locus-specific DNA target preparation method for highly multiplexed single nucleotide polymorphism (SNP) genotyping called MARA (Multiplexed Anchored Runoff Amplification). The approach uses a single primer per SNP in conjunction with restriction enzyme digested, adapter-ligated human genomic DNA. Each primer is composed of common sequence at the 5′ end followed by locus-specific sequence at the 3′ end. Following a primary reaction in which locus-specific products are generated, a secondary universal amplification is carried out using a generic primer pair corresponding to the oligonucleotide and genomic DNA adapter sequences. Allele discrimination is achieved by hybridization to high-density DNA oligonucleotide arrays. Initial multiplex reactions containing either 250 primers or 750 primers across nine DNA samples demonstrated an average sample call rate of ~95% for 250- and 750-plex MARA. We have also evaluated >1000- and 4000-primer plex MARA to genotype SNPs from human chromosome 21. We have identified a subset of SNPs corresponding to a primer conversion rate of ~75%, which show an average call rate over 95% and concordance >99% across seven DNA samples. Thus, MARA may potentially improve the throughput of SNP genotyping when coupled with allele discrimination on high-density arrays by allowing levels of multiplexing during target generation that far exceed the capacity of traditional multiplex PCR.