Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species.

Studies of population genetics increasingly use next-generation DNA sequencing to identify microsatellite loci in nonmodel organisms. There are, however, relatively few studies that validate the feasibility of transitioning from marker development to experimental application across populations and species. North American coralsnakes of the Micrurus fulvius species complex occur in the United States and Mexico, and little is known about their population structure and phylogenetic relationships. This absence of information and population genetics markers is particularly concerning because they are highly venomous and have important implications on human health. To alleviate this problem in coralsnakes, we investigated the feasibility of using 454 shotgun sequences for microsatellite marker development. First, a genomic shotgun library from a single individual was sequenced (approximately 7.74 megabases; 26,831 reads) to identify potentially amplifiable microsatellite loci (PALs). We then hierarchically sampled 76 individuals from throughout the geographic distribution of the species complex and examined whether PALs were amplifiable and polymorphic. Approximately half of the loci tested were readily amplifiable from all individuals, and 80% of the loci tested for variation were variable and thus informative as population genetic markers. To evaluate the repetitive landscape characteristics across multiple snakes, we also compared microsatellite content between the coralsnake and two other previously sampled snakes, the venomous copperhead (Agkistrodon contortrix) and Burmese python (Python molurus).

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

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.

A multi-organ transcriptome resource for the Burmese Python (Python molurus bivittatus).

BACKGROUND: Snakes provide a unique vertebrate system for studying a diversity of extreme adaptations, including those related to development, metabolism, physiology, and venom. Despite their importance as research models, genomic resources for snakes are few. Among snakes, the Burmese python is the premier model for studying extremes of metabolic fluctuation and physiological remodelling. In this species, the consumption of large infrequent meals can induce a 40-fold increase in metabolic rate and more than a doubling in size of some organs. To provide a foundation for research utilizing the python, our aim was to assemble and annotate a transcriptome reference from the heart and liver. To accomplish this aim, we used the 454-FLX sequencing platform to collect sequence data from multiple cDNA libraries. RESULTS: We collected nearly 1 million 454 sequence reads, and assembled these into 37,245 contigs with a combined length of 13,409,006 bp. To identify known genes, these contigs were compared to chicken and lizard gene sets, and to all Genbank sequences. A total of 13,286 of these contigs were annotated based on similarity to known genes or Genbank sequences. We used gene ontology (GO) assignments to characterize the types of genes in this transcriptome resource. The raw data, transcript contig assembly, and transcript annotations are made available online for use by the broader research community. CONCLUSION: These data should facilitate future studies using pythons and snakes in general, helping to further contribute to the utilization of snakes as a model evolutionary and physiological system. This sequence collection represents a major genomic resource for the Burmese python, and the large number of transcript sequences characterized should contribute to future research in this and other snake species.

Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing.

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.

Rapid identification of thousands of copperhead snake (Agkistrodon contortrix) microsatellite loci from modest amounts of 454 shotgun genome sequence.

Optimal integration of next-generation sequencing into mainstream research requires re-evaluation 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 sample randomly ∼27 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, 4564 of which had flanking sequences suitable for polymerase chain reaction 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.

Effects of abiotic factors on the phylogenetic diversity of bacterial communities in acidic thermal springs.

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.

Protein activation of a ribozyme: the role of bacterial RNase P protein.

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.