Pyrosequencing of the northern red oak (Quercus rubra L.) chloroplast genome reveals high quality polymorphisms for population management

Given the low intraspecific chloroplast diversity detected in northern red oak (Quercus rubra L.), more powerful genetic tools are necessary to accurately characterize Q. rubra chloroplast diversity and structure. We report the sequencing, assembly, and annotation of the chloroplast genome of northern red oak via pyrosequencing and a combination of de novo and reference-guided assembly (RGA). Chloroplast DNA from 16 individuals was separated into four MID-tagged pools for a Genome Sequencer 20 quarter-run (Roche Life Sciences, Indianapolis, IN, USA). A four-step assembly method was used to generate the Q. rubra chloroplast consensus sequence: (1) reads were assembled de novo into contigs, (2) de novo contigs were aligned to a reference genome and merged to produce a consensus sequence, (3) the consensus sequence was aligned to the reference sequence and gaps between contigs were filled with reference sequence to generate a "pseudoreference", and (4) reads were mapped to the pseudoreference using RGA to generate the draft chloroplast genome. One hundred percent of the pseudoreference sequence was covered with a minimum coverage of 2× and an average coverage of 43.75×. The 161,304-bp Q. rubra chloroplast genome draft sequence contained 137 genes and one rps19 pseudogene. The sequence was compared to that of Quercus robur and Q. nigra with 951 and 186 insertion/deletion or SNP polymorphisms detected, respectively. A total of 51 intraspecific polymorphisms were detected among four northern red oak individuals. The fully sequenced and annotated Q. rubra chloroplast genome containing locations of interspecific and intraspecific polymorphisms will be essential for studying population differentiation, phylogeography, and evolutionary history of this species as well as meeting management goals such as monitoring reintroduced populations, tracking wood products, and certifying seed lots and forests.     

Genetic structure of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill populations at gene-based EST-SSR and nuclear SSR markers

Sympatric hybridizing oak species provide a model system for studying local adaptation. Disjunct populations of Quercus rubra L. and Quercus ellipsoidalis E. J. Hill at the northern edge of their distribution may harbor important reservoirs of adaptive genetic variation. Genic (expressed sequence tag- simple sequence repeat?=?EST-SSR) and non-genic nuclear microsatellite (nuclear SSR?=?nSSR) markers were used to estimate neutral and potentially adaptive genetic variation in these two supposedly interfertile oak species showing different adaptations to drought. Eleven populations of putative Q. rubra and Q. ellipsoidalis located in the Western Upper Peninsula of Michigan were characterized using seven EST-SSRs and eight nSSRs. Bayesian cluster analysis revealed two distinct groups corresponding to each species with evidence of low levels of potential introgression. A comparison of the genetic structure of adult trees and seedlings revealed no evidence for selection against hybrids. Overall, similar levels of genetic variation and differentiation between populations and species were found at both EST-SSRs and nSSRs indicating that most EST-SSRs chosen reflect neutral variation. Two loci, 3A05 (nSSR) and GOT021 (EST-SSR, putative histidine kinase 4-like), were identified as putative outlier loci between species showing largely reduced variation in Q. ellipsoidalis. Future analyses of an increased number of EST-SSRs located in functional genes will allow the identification of genes involved in the reproductive isolation between both species.     

The development and evaluation of consensus chloroplast primer pairs that possess highly variable sequence regions in a diverse array of plant taxa

Although universal or consensus chloroplast primers are available, they are limited by their number and genomic distribution. Therefore, a set of consensus chloroplast primer pairs for simple sequence repeats (ccSSRs) analysis was constructed from tobacco (Nicotiana tabacum L.) chloroplast sequences. These were then tested for their general utility in the genetic analysis of a diverse array of plant taxa. In order to increase the number of ccSSRs beyond that previously reported, the target sequences for SSR motifs was set at A or T (n 7) mononucleotide repeats. Each SSR sequence motif, along with ±200-bp flanking sequences from the first of each mononucleotide base repeat, was screened for homologies with chloroplast DNA sequences of other plant species in GenBank databases using BLAST search procedures. Twenty three putative marker loci that possessed conserved flanking sequence spans were selected for consensus primer pair construction using commercially available computer algorithms. All primer pairs produced amplicons after PCR employing genomic DNA from members of the Cucurbitaceae (six species) and Solanaceae (four species). Sixteen, 22 and 19 of the initial 23 primer pairs were successively amplified by PCR using template DNA from species of the Apiaceae (two species), Brassicaceae (one species) and Fabaceae (two species), respectively. Twenty of 23 primer pairs were also functional in three monocot species of the Liliaceae [onion (Allium cepa L.) and garlic (Allium sativum L.)], and the Poaceae [oat (Avena sativa L.)]. Sequence analysis of selected ccSSR fragments suggests that ccSSR length and sequence variation could be useful as a tool for investigating the genetic relationships within a genus or closely related taxa (i.e., tribal level). In order to provide for a marker system having significant coverage of the cucumber chloroplast genome, ccSSR primers were strategically "recombined" and named recombined consensus chloroplast primers (RCCP) for PCR analysis. Successful amplification after extended-length PCR of 16 RCCP primer pairs from cucumber (Cucumis sativus L.) DNA suggested that the amplicons detected are representative of the cucumber chloroplast genome. These RCCP pairs, therefore, could be useful as an initial molecular tool for investigation of traits related to a chloroplast gene(s) in cucumber, and other closely related species.Communicated by C. Möllers     

Development and Characterization of Genomic and Gene-Based Microsatellite Markers in North American Red Oak Species

Oaks (Quercus: Fagaceae) are ecological and economic keystones of many forested ecosystems but effective genetic management strategies are hindered by high levels of phenotypic plasticity within species and frequent hybridization among them. These same features, however, make oak communities suited for the study of speciation, hybridization, and genetic adaptation. Efforts to develop new and to adapt existing genomic resources to less-studied members of this genus should not only improve oak conservation and management but also aid the study of fundamental evolutionary processes. Here, we present a suite of 27 highly polymorphic simple sequence repeat (SSR) markers tested in four North American red oak (Quercus section Lobatae) species: Q. rubra, Q. ellipsoidalis, Q. coccinea, and Q. velutina. Five markers are genomic SSRs (gSSRs) — four novel and one previously transferred from Q. petraea — and 22 are gene-based SSRs derived from Q. robur and Q. petraea expressed sequence tags (EST-SSRs). Overall, levels of polymorphism detected with these primer pairs were high, with gene diversity (H _e) averaging 0.66 across all loci in natural populations. In addition, we show that EST-SSR markers may have the potential to detect divergent selection at stress-resistance candidate genes among closely related oak species.     

Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies

16S ribosomal RNA gene (rDNA) amplicon analysis remains the standard approach for the cultivation-independent investigation of microbial diversity. The accuracy of these analyses depends strongly on the choice of primers. The overall coverage and phylum spectrum of 175 primers and 512 primer pairs were evaluated in silico with respect to the SILVA 16S/18S rDNA non-redundant reference dataset (SSURef 108 NR). Based on this evaluation a selection of ‘best available’ primer pairs for Bacteria and Archaea for three amplicon size classes (100–400, 400–1000, ≥1000 bp) is provided. The most promising bacterial primer pair (S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21), with an amplicon size of 464 bp, was experimentally evaluated by comparing the taxonomic distribution of the 16S rDNA amplicons with 16S rDNA fragments from directly sequenced metagenomes. The results of this study may be used as a guideline for selecting primer pairs with the best overall coverage and phylum spectrum for specific applications, therefore reducing the bias in PCR-based microbial diversity studies.     

Evidence for selection on a CONSTANS-like gene between two red oak species

Background and Aims

Hybridizing species such as oaks may provide a model to study the role of selection in speciation with gene flow. Discrete species'' identities and different adaptations are maintained among closely related oak species despite recurrent gene flow. This is probably due to ecologically mediated selection at a few key genes or genomic regions. Neutrality tests can be applied to identify so-called outlier loci, which demonstrate locus-specific signatures of divergent selection and are candidate genes for further study.

Methods

Thirty-six genic microsatellite markers, some with putative functions in flowering time and drought tolerance, and eight non-genic microsatellite markers were screened in two population pairs (n = 160) of the interfertile species Quercus rubra and Q. ellipsoidalis, which are characterized by contrasting adaptations to drought. Putative outliers were then tested in additional population pairs from two different geographic regions (n = 159) to support further their potential role in adaptive divergence.

Key Results

A marker located in the coding sequence of a putative CONSTANS-like (COL) gene was repeatedly identified as under strong divergent selection across all three geographically disjunct population pairs. COL genes are involved in the photoperiodic control of growth and development and are implicated in the regulation of flowering time.

Conclusions

The location of the polymorphism in the Quercus COL gene and given the potential role of COL genes in adaptive divergence and reproductive isolation makes this a promising candidate speciation gene. Further investigation of the phenological characteristics of both species and flowering time pathway genes is suggested in order to elucidate the importance of phenology genes for the maintenance of species integrity. Next-generation sequencing in multiple population pairs in combination with high-density genetic linkage maps could reveal the genome-wide distribution of outlier genes and their potential role in reproductive isolation between these species.     

Single nucleotide polymorphism-mismatch primer development for rapid molecular identification of selected microalgal species

Microalgae have been a great source for food, cosmetic, pharmacological, and biofuel production. The adoption of effective diagnostic assays for monitoring all stages of algal cultivation has become essential. In addition to microscopy identification, molecular assays can aid greatly in the identification and monitoring of algal species of interest. In this study the 18S ribosomal RNA (rRNA) sequences of 12 microalgal species and/or strains were used to design algal identification primers. Sequence alignment revealed five highly variable regions and multiple unique single nucleotide polymorphisms (SNPs). To design target algae specific primers, a SNP identified as unique to each microalgal species was incorporated into the 3’-terminus of forward and reverse primer pairs, respectively. To further enhance primer specificity, transverse mutation was introduced into each primer at the third base upstream of the respective SNP. The SNP-mismatch primer pairs yield size-specific amplicons, enabling the rapid molecular detection of 12 microalgae by circumventing cloning and sequencing. To verify the primer specificity, two SNP-mismatch primer pairs designed for Chlorella sorokiniana DOE1412 and wildtype species of Scenedesmus were tested in the outdoor reactor run inoculated with C. sorokiniana DOE1412. The primer pairs were able to identify C. sorokiniana DOE1412 as well as the environmental invader Scenedesmus sp. Furthermore, the “relative concentration” of two microalgae was accessed throughout the entire cultivation run. The use of SNPs primers designed in this study offers a cost-effective, easy to use alternative for routine monitoring of microalgal cultures in laboratories, in scale-ups, and in cultivation reactors, independent of the production platform.

     

Taxonomic individuality of Leonurus cardiaca and Leonurus quinquelobatus in view of morphological and molecular studies

The main goal of this study was to determine the number and taxonomic rank of taxa belonging to the complex Leonurus cardiaca agg. in Poland. Based on statistical analysis of selected features, two morphological forms of this plant were distinguished. In order to determine their genetic polymorphism and the relationships between them, the nuclear, mitochondrial and chloroplast genomes were analysed with the use of RAPD and PCR–RFLP markers. 39 RAPD primers produced a total of 234 nuclear DNA fragments, of which 128 were polymorphic and distributed almost equally between two forms. It was found that 87 % of the compared pairs of RAPD profiles differ from each other. Five chloroplast and two mitochondrial primer pairs were used to amplify non-coding regions of organelle genomes. Restriction analysis revealed uniformity of mtDNA and occurrence of two cpDNA haplotypes, corresponding to naked and hairy forms of L. cardiaca agg. The obtained results justifies the recognition of these forms as separate species L. cardiaca s. s. L. and L. quinquelobatus Gilib., respectively. The distribution of both species in Poland is given in the paper.     

New Strategy for Identification of Novel cry-Type Genes from Bacillus thuringiensis Strains

We designed five degenerate primers for detection of novel cry genes from Bacillus thuringiensis strains. An efficient strategy was developed based on a two-step PCR approach with these primers in five pair combinations. In the first step, only one of the primer pairs is used in the PCR, which allows amplification of DNA fragments encoding protein regions that include consensus domains of representative proteins belonging to different Cry groups. A second PCR is performed by using the first-step amplification products as DNA templates and the set of five primer combinations. Cloning and sequencing of the last-step amplicons allow both the identification of known cry genes encoding Cry proteins covering a wide phylogenetic distance and the detection and characterization of cry-related sequences from novel B. thuringiensis isolates.     

A new set of mono‐ and dinucleotide chloroplast microsatellites in Fagaceae

A new set of 14 chloroplast microsatellites, represented by mono‐ and dinucleotide repeats, was optimized in the three main species of the Fagaceae (Castanea sativa, Fagus sylvatica and Quercus petraea). The intraspecific variation was tested in some natural populations. The polymorphic microsatellites displayed two or three variants. Conservation of the primer pairs was checked on an additional set of six species of the Fagaceae and on Fraxinus excelsior. All the primer pairs produced a fragment of the expected size in the Fagaceae species while no amplification was obtained with 36% of the primers in F. excelsior.     

The chloroplast genome sequence of Syzygium cumini (L.) and its relationship with other angiosperms

Only a few studies to date have conducted comparative genomics in the Myrtaceae family. Here, we report the complete sequence and bioinformatics analysis of the chloroplast genome of Syzygium cumini (L.), one of the family members. The size of S. cumini cp genome was within the range of reported angiosperm chloroplast genomes. Comparison of S. cumini cpDNA sequence with previously reported partial sequences of S. cumini revealed several SNPs that resulted in non-synonymous mutations in maturase K and NADH-plastoquinone oxidoreductase subunit-5. These polymorphic characters might serve as intra-specific markers to address whether lineage sorting from polymorphic ancestry has occurred. Comparison of the S. cumini chloroplast genome with related dicots revealed an expansion in the intergenic spacer located between IRA/large single copy (LSC) border and the first gene of LSC region, driven by sequence of 54 bp. This type of variation in the intergenic regions can be utilized in the development of species-specific vectors for chloroplast genetic engineering. Several of the longer (30–40 bp) repeats were found to be conserved in other dicot species, suggesting that they might be widespread in angiosperm chloroplast genomes.     

Complete Chloroplast Genome of Tanaecium tetragonolobum: The First Bignoniaceae Plastome

Bignoniaceae is a Pantropical plant family that is especially abundant in the Neotropics. Members of the Bignoniaceae are diverse in many ecosystems and represent key components of the Tropical flora. Despite the ecological importance of the Bignoniaceae and all the efforts to reconstruct the phylogeny of this group, whole chloroplast genome information has not yet been reported for any members of the family. Here, we report the complete chloroplast genome sequence of Tanaecium tetragonolobum (Jacq.) L.G. Lohmann, which was reconstructed using de novo and referenced-based assembly of single-end reads generated by shotgun sequencing of total genomic DNA in an Illumina platform. The gene order and organization of the chloroplast genome of T. tetragonolobum exhibits the general structure of flowering plants, and is similar to other Lamiales chloroplast genomes. The chloroplast genome of T. tetragonolobum is a circular molecule of 153,776 base pairs (bp) with a quadripartite structure containing two single copy regions, a large single copy region (LSC, 84,612 bp) and a small single copy region (SSC, 17,586 bp) separated by inverted repeat regions (IRs, 25,789 bp). In addition, the chloroplast genome of T. tetragonolobum has 38.3% GC content and includes 121 genes, of which 86 are protein-coding, 31 are transfer RNA, and four are ribosomal RNA. The chloroplast genome of T. tetragonolobum presents a total of 47 tandem repeats and 347 simple sequence repeats (SSRs) with mononucleotides being the most common and di-, tri-, tetra-, and hexanucleotides occurring with less frequency. The results obtained here were compared to other chloroplast genomes of Lamiales available to date, providing new insight into the evolution of chloroplast genomes within Lamiales. Overall, the evolutionary rates of genes in Lamiales are lineage-, locus-, and region-specific, indicating that the evolutionary pattern of nucleotide substitution in chloroplast genomes of flowering plants is complex. The discovery of tandem repeats within T. tetragonolobum and the presence of divergent regions between chloroplast genomes of Lamiales provides the basis for the development of markers at various taxonomic levels. The newly developed markers have the potential to greatly improve the resolution of molecular phylogenies.     

Morphometric variation in oaks of the Apostle Islands in Wisconsin: evidence of hybridization between Quercus rubra and Q. ellipsoidalis (Fagaceae)

The Apostle Islands in Lake Superior are populated by trees that are clearly related to Quercus rubra L. However, several islands have trees with morphological characteristics suggestive of hybridization with Q. ellipsoidalis Hill. Leaf specimens were collected from trees in five locations: the outermost island, an intermediate island, the nearest-shore island, the northeast shoreline, and an inland forest about 24 km from the shoreline. Seventeen landmarks were digitized for two to five leaves per tree. These landmarks were used to generate nine linear characters and three angles. These characters, along with the number of bristle tips per leaf, were used in various combinations for several principal component analyses. In addition, the landmark configurations were examined using rotational-fit methods. The patterns observed in both types of analysis indicate phenotypic variation coincident with a line connecting the two most distant sample sites. The location nearest the geographic center of this line is also nearest the center of the two-dimensional view of phenotypic variation. Trees at each site illustrate a distinctive pattern in the rotational-fit analyses, and patterns of co-variation in the morphometric characters are different for each site. The observed morphometric variation is consistent with the hypothesis that there is hybridization between these two species, most likely in the form of introgression from Q. ellipsoidalis into Q. rubra.     

The Complete Chloroplast Genome of Banana (Musa acuminata,Zingiberales): Insight into Plastid Monocotyledon Evolution

Background

Banana (genus Musa) is a crop of major economic importance worldwide. It is a monocotyledonous member of the Zingiberales, a sister group of the widely studied Poales. Most cultivated bananas are natural Musa inter-(sub-)specific triploid hybrids. A Musa acuminata reference nuclear genome sequence was recently produced based on sequencing of genomic DNA enriched in nucleus.

Methodology/Principal Findings

The Musa acuminata chloroplast genome was assembled with chloroplast reads extracted from whole-genome-shotgun sequence data. The Musa chloroplast genome is a circular molecule of 169,972 bp with a quadripartite structure containing two single copy regions, a Large Single Copy region (LSC, 88,338 bp) and a Small Single Copy region (SSC, 10,768 bp) separated by Inverted Repeat regions (IRs, 35,433 bp). Two forms of the chloroplast genome relative to the orientation of SSC versus LSC were found. The Musa chloroplast genome shows an extreme IR expansion at the IR/SSC boundary relative to the most common structures found in angiosperms. This expansion consists of the integration of three additional complete genes (rps15, ndhH and ycf1) and part of the ndhA gene. No such expansion has been observed in monocots so far. Simple Sequence Repeats were identified in the Musa chloroplast genome and a new set of Musa chloroplastic markers was designed.

Conclusion

The complete sequence of M. acuminata ssp malaccensis chloroplast we reported here is the first one for the Zingiberales order. As such it provides new insight in the evolution of the chloroplast of monocotyledons. In particular, it reinforces that IR/SSC expansion has occurred independently several times within monocotyledons. The discovery of new polymorphic markers within Musa chloroplast opens new perspectives to better understand the origin of cultivated triploid bananas.     

A set of primers for the amplification of chloroplast microsatellites in Quercus

The increase in demand for the certification of oak seed lots, as well as control of the geographical origin of oak wood, has led us to develop powerful genetic markers permitting us to discriminate among provenance regions. With the aim of detecting new chloroplast variants, we have identified 17 potential cpSSRs motifs from available oak sequences and tested their variability among French oak populations. Six loci were polymorphic at the intraspecific level in Quercus petraea and Q. robur. Moreover, conservation of the primer pairs was checked on a set of 21 forest tree species and they were all shown to work well on several Quercus species, and even within Fagacaea.     

Development of high‐resolution DNA barcodes for Dioscorea species discrimination and phylogenetic analysis

The genus Dioscorea is widely distributed in tropical and subtropical regions, and is economically important in terms of food supply and pharmaceutical applications. However, DNA barcodes are relatively unsuccessful in discriminating between Dioscorea species, with the highest discrimination rate (23.26%) derived from matK sequences. In this study, we compared genic and intergenic regions of three Dioscorea chloroplast genomes and found that the density of SNPs and indels in intergenic sites was about twice and seven times higher than that of SNPs and indels in the genic regions, respectively. A total of 52 primer pairs covering highly variable regions were designed and seven pairs of primers had 80%–100% PCR success rate. PCR amplicons of 73 Dioscorea individuals and assembled sequences of 47 Dioscorea SRAs were used for estimating intraspecific and interspecific divergence for the seven loci: The rpoB‐trnC locus had the highest interspecific divergence. Automatic barcoding gap discovery (ABGD), Poisson tree processes (PTP), and generalized mixed Yule coalescence (GMYC) analysis were applied for species delimitation based on the seven loci and successfully identified the majority of species, except for species in the Enantiophyllum section. Phylogenetic analysis of 51 Dioscorea individuals (28 species) showed that most individuals belonging to the same species tended to cluster in the same group. Our results suggest that the variable loci derived from comparative analysis of plastid genome sequences could be good DNA barcode candidates for taxonomic analysis and species delimitation.     

16S rRNA Gene-Based Detection of Tetrachloroethene-Dechlorinating Desulfuromonas and Dehalococcoides Species

Members of the genera Desulfuromonas and Dehalococcoides reductively dechlorinate tetrachloroethene (PCE) and trichloroethene. Two primer pairs specific to hypervariable regions of the 16S rRNA genes of the Dehalococcoides group (comprising Dehalococcoides ethenogenes and Dehalococcoides sp. strain FL2) and the acetate-oxidizing, PCE-dechlorinating Desulfuromonas group (comprising Desulfuromonas sp. strain BB1 and Desulfuromonas chloroethenica) were designed. The detection threshold of a nested PCR approach using universal bacterial primers followed by a second PCR with the Desulfuromonas dechlorinator-targeted primer pair was 1 × 10³ BB1 cells added per gram (wet weight) of sandy aquifer material. Total community DNA isolated from sediments of three Michigan rivers and six different chloroethene-contaminated aquifer samples was used as template in nested PCR. All river sediment samples yielded positive signals with the BB1- and the Dehalococcoides-targeted primers. One chloroethene-contaminated aquifer tested positive with the Dehalococcoides-targeted primers, and another contaminated aquifer tested positive with the Desulfuromonas dechlorinator-targeted primer pair. Restriction fragment analysis of the amplicons could discriminate strain BB1 from other known Desulfuromonas species. Microcosm studies confirmed the presence of PCE-dechlorinating, acetate-oxidizing Desulfuromonas and hydrogenotrophic Dehalococcoides species in samples yielding positive PCR signals with the specific primers.     

A PCR-based toolbox for the culture-independent quantification of total bacterial abundances in plant environments

A major obstacle in the culture-independent estimation of the abundance of bacteria associated with plants is contamination with plant organelles, which precludes the use of universal rRNA bacterial primers in quantitative PCR applications. We present here a PCR-based method that allows a priori determination of the degree of chloroplast and mitochondrial contamination in DNA samples from plant environments. It is based on differential digestibility of chloroplast, mitochondrial and bacterial small subunit rRNA gene amplicons with the restriction enzymes AfeI and BbvCI. Using this method, we demonstrated for field-grown lettuce plants that even a gentle washing protocol, designed to recover the microbial community and its metagenome from the leaf surface, resulted in substantial contamination with chloroplast DNA. This finding cautions against the use of universal primer pairs that do not exclude chloroplast DNA from amplification, because they risk overestimation of bacterial population sizes. In contrast, contamination with mitochondrial 18S rRNA was minor in the lettuce phyllosphere. These findings were confirmed by real-time PCR using primer sets specific for small subunit rRNA genes from bacteria, chloroplasts, and mitochondria. Based on these results, we propose two primer pairs (534f/783r and mito1345f/mito1430r) which between them offer an indirect means of faithfully estimating bacterial abundances on plants, by deduction of the mito1345f/mito1430r-based mitochondrial count from that obtained with 534f/783r, which amplifies both bacterial and mitochondrial DNA but excludes chloroplast. In this manner, we estimated the number of total bacteria on most leaves of field-grown lettuce to be between 10⁵ and 10⁶ g^− 1 of leaf, which was 1-3 orders of magnitudes higher than the number of colony-forming units that were retrieved from the same leaf surfaces on agar plates.     

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.