Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct “Seq-to-SSR” approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Optimal integration of next-generation sequencing into mainstream research requires reevaluation of how problems can be reasonably overcome and what questions can be asked. One potential application is the rapid acquisition of genomic information to identify microsatellite loci for evolutionary, population genetic and chromosome linkage mapping research on non-model and not previously sequenced organisms. Here, we report on results using high-throughput sequencing to obtain a large number of microsatellite loci from the venomous snake Agkistrodon contortrix, the copperhead. We used the 454 Genome Sequencer FLX next-generation sequencing platform to randomly sample approximately 26.8 Mbp (128,773 reads) of the copperhead genome, thus sampling about 2% of the genome of this species. We identified microsatellite loci in 11.3% of all reads obtained, with 14,612 microsatellite loci identified in total, 4,564 of which had flanking sequences suitable for PCR primer design. The random sequencing-based approach to identify microsatellites was rapid, cost-effective, and identified thousands of useful microsatellite loci in a previously unstudied species.

We conducted a comprehensive assessment of genomic repeat content in two snake genomes, the venomous copperhead (Agkistrodon contortrix) and the Burmese python (Python molurus bivittatus). These two genomes are both relatively small (∼1.4 Gb) but have surprisingly extensive differences in the abundance and expansion histories of their repeat elements. In the python, the readily identifiable repeat element content is low (21%), similar to bird genomes, whereas that of the copperhead is higher (45%), similar to mammalian genomes. The copperhead's greater repeat content arises from the recent expansion of many different microsatellites and transposable element (TE) families, and the copperhead had 23-fold greater levels of TE-related transcripts than the python. This suggests the possibility that greater TE activity in the copperhead is ongoing. Expansion of CR1 LINEs in the copperhead genome has resulted in TE-mediated microsatellite expansion (“microsatellite seeding”) at a scale several orders of magnitude greater than previously observed in vertebrates. Snakes also appear to be prone to horizontal transfer of TEs, particularly in the copperhead lineage. The reason that the copperhead has such a small genome in the face of so much recent expansion of repeat elements remains an open question, although selective pressure related to extreme metabolic performance is an obvious candidate. TE activity can affect gene regulation as well as rates of recombination and gene duplication, and it is therefore possible that TE activity played a role in the evolution of major adaptations in snakes; some evidence suggests this may include the evolution of venom repertoires.

Bacterial ribonuclease P (RNase P) belongs to a class of enzymes that utilize both RNAs and proteins to perform essential cellular functions. The bacterial RNase P protein is required to activate bacterial RNase P RNA in vivo, but previous studies have yielded contradictory conclusions regarding its specific functions. Here, we use biochemical and biophysical techniques to examine all of the proposed functions of the protein in both Escherichia coli and Bacillus subtilis RNase P. We demonstrate that the E. coli protein, but not the B. subtilis protein, stabilizes the global structure of RNase P RNA, although both proteins influence holoenzyme dimer formation and precursor tRNA recognition to different extents. By comparing each protein in complex with its cognate and noncognate RNA, we show that differences between the two types of holoenzymes reside primarily in the RNA and not the protein components of each. Our results reconcile previous contradictory conclusions regarding the role of the protein and support a model where the protein activates local RNA structures that manifest multiple holoenzyme properties.

Acidic thermal springs offer ideal environments for studying processes underlying extremophile microbial diversity. We used a carefully designed comparative analysis of acidic thermal springs in Yellowstone National Park to determine how abiotic factors (chemistry and temperature) shape acidophile microbial communities. Small-subunit rRNA gene sequences were PCR amplified, cloned, and sequenced, by using evolutionarily conserved bacterium-specific primers, directly from environmental DNA extracted from Amphitheater Springs and Roaring Mountain sediment samples. Energy-dispersive X-ray spectroscopy, X-ray diffraction, and colorimetric assays were used to analyze sediment chemistry, while an optical emission spectrometer was used to evaluate water chemistry and electronic probes were used to measure the pH, temperature, and E h of the spring waters. Phylogenetic-statistical analyses found exceptionally strong correlations between bacterial community composition and sediment mineral chemistry, followed by weaker but significant correlations with temperature gradients. For example, sulfur-rich sediment samples contained a high diversity of uncultured organisms related to Hydrogenobaculum spp., while iron-rich sediments were dominated by uncultured organisms related to a diverse array of gram-positive iron oxidizers. A detailed analysis of redox chemistry indicated that the available energy sources and electron acceptors were sufficient to support the metabolic potential of Hydrogenobaculum spp. and iron oxidizers, respectively. Principal-component analysis found that two factors explained 95% of the genetic diversity, with most of the variance attributable to mineral chemistry and a smaller fraction attributable to temperature.