Selective Bacterial Community Enrichment between the Pitcher Plants Sarracenia minor and Sarracenia flava.

The interconnected and overlapping habitats present in natural ecosystems remain a challenge in determining the forces driving microbial community composition. The cuplike leaf structures of some carnivorous plants, including those of the family Sarraceniaceae, are self-contained ecological habitats that represent systems for exploring such microbial ecology questions. We investigated whether Sarracenia minor and Sarracenia flava cultivate distinct bacterial communities when sampled at the same geographic location and time. This sampling strategy eliminates many abiotic environmental variables present in other studies that compare samples harvested over time, and it could reveal biotic factors driving the selection of microbes. DNA extracted from the decomposing detritus trapped in each Sarracenia leaf pitcher was profiled using 16S rRNA amplicon sequencing. We identified a surprising amount of bacterial diversity within each pitcher, but we also discovered bacteria whose abundance was specifically enriched in one of the two Sarracenia species. These differences in bacterial community representation suggest some biotic influence of the Sarracenia plant on the bacterial composition of their pitchers. Overall, our results suggest that bacterial selection due to factors other than geographic location, weather, or prey availability is occurring within the pitchers of these two closely related plant species. This indicates that specific characteristics of S. minor and S. flava may play a role in fostering distinct bacterial communities. These confined, naturally occurring microbial ecosystems within Sarracenia pitchers may provide model systems to answer important questions about the drivers of microbial community composition, succession, and response to environmental perturbations. IMPORTANCE This study uses amplicon sequencing to compare the bacterial communities of environmental samples from the detritus of the leaf cavities of Sarracenia minor and Sarracenia flava pitcher plants. We sampled the detritus at the same time and in the same geographic location, eliminating many environmental variables present in other comparative studies. This study revealed that different species of Sarracenia contain distinct bacterial members within their pitchers, suggesting that these communities are not randomly established based on environmental factors and the prey pool but are potentially enriched for by the plants' chemical or physical environment. This study of these naturally occurring, confined microbial ecosystems will help further establish carnivorous pitcher plants as a model system for answering important questions about the development and succession of microbial communities.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

MT-Toolbox: improved amplicon sequencing using molecule tags.

BACKGROUND: Short oligonucleotides can be used as markers to tag and track DNA sequences. For example, barcoding techniques (i.e. Multiplex Identifiers or Indexing) use short oligonucleotides to distinguish between reads from different DNA samples pooled for high-throughput sequencing. A similar technique called molecule tagging uses the same principles but is applied to individual DNA template molecules. Each template molecule is tagged with a unique oligonucleotide prior to polymerase chain reaction. The resulting amplicon sequences can be traced back to their original templates by their oligonucleotide tag. Consensus building from sequences sharing the same tag enables inference of original template molecules thereby reducing effects of sequencing error and polymerase chain reaction bias. Several independent groups have developed similar protocols for molecule tagging; however, user-friendly software for build consensus sequences from molecule tagged reads is not readily available or is highly specific for a particular protocol. RESULTS: MT-Toolbox recognizes oligonucleotide tags in amplicons and infers the correct template sequence. On a set of molecule tagged test reads, MT-Toolbox generates sequences having on average 0.00047 errors per base. MT-Toolbox includes a graphical user interface, command line interface, and options for speed and accuracy maximization. It can be run in serial on a standard personal computer or in parallel on a Load Sharing Facility based cluster system. An optional plugin provides features for common 16S metagenome profiling analysis such as chimera filtering, building operational taxonomic units, contaminant removal, and taxonomy assignments. CONCLUSIONS: MT-Toolbox provides an accessible, user-friendly environment for analysis of molecule tagged reads thereby reducing technical errors and polymerase chain reaction bias. These improvements reduce noise and allow for greater precision in single amplicon sequencing experiments.

Development and mapping of SNP assays in allotetraploid cotton.

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin-streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N(50) of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F(2) population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome ('A' or 'D') each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.