Juan M. Daza

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.

The Burmese python genome reveals the molecular basis for extreme adaptation in snakes

Significance The molecular basis of morphological and physiological adaptations in snakes is largely unknown. Here, we study these phenotypes using the genome of the Burmese python (Python molurus bivittatus), a model for extreme phenotypic plasticity and metabolic adaptation. We discovered massive rapid changes in gene expression that coordinate major changes in organ size and function after feeding. Many significantly responsive genes are associated with metabolism, development, and mammalian diseases. A striking number of genes experienced positive selection in ancestral snakes. Such genes were related to metabolism, development, lungs, eyes, heart, kidney, and skeletal structure—all highly modified features in snakes. Snake phenotypic novelty seems to be driven by the system-wide coordination of protein adaptation, gene expression, and changes in genome structure. Snakes possess many extreme morphological and physiological adaptations. Identification of the molecular basis of these traits can provide novel understanding for vertebrate biology and medicine. Here, we study snake biology using the genome sequence of the Burmese python (Python molurus bivittatus), a model of extreme physiological and metabolic adaptation. We compare the python and king cobra genomes along with genomic samples from other snakes and perform transcriptome analysis to gain insights into the extreme phenotypes of the python. We discovered rapid and massive transcriptional responses in multiple organ systems that occur on feeding and coordinate major changes in organ size and function. Intriguingly, the homologs of these genes in humans are associated with metabolism, development, and pathology. We also found that many snake metabolic genes have undergone positive selection, which together with the rapid evolution of mitochondrial proteins, provides evidence for extensive adaptive redesign of snake metabolic pathways. Additional evidence for molecular adaptation and gene family expansions and contractions is associated with major physiological and phenotypic adaptations in snakes; genes involved are related to cell cycle, development, lungs, eyes, heart, intestine, and skeletal structure, including GRB2-associated binding protein 1, SSH, WNT16, and bone morphogenetic protein 7. Finally, changes in repetitive DNA content, guanine-cytosine isochore structure, and nucleotide substitution rates indicate major shifts in the structure and evolution of snake genomes compared with other amniotes. Phenotypic and physiological novelty in snakes seems to be driven by system-wide coordination of protein adaptation, gene expression, and changes in the structure of the genome.

Comparative phylogeographic analyses illustrate the complex evolutionary history of threatened cloud forests of northern Mesoamerica.

Comparative phylogeography can elucidate the influence of historical events on current patterns of biodiversity and can identify patterns of co-vicariance among unrelated taxa that span the same geographic areas. Here we analyze temporal and spatial divergence patterns of cloud forest plant and animal species and relate them to the evolutionary history of naturally fragmented cloud forests–among the most threatened vegetation types in northern Mesoamerica. We used comparative phylogeographic analyses to identify patterns of co-vicariance in taxa that share geographic ranges across cloud forest habitats and to elucidate the influence of historical events on current patterns of biodiversity. We document temporal and spatial genetic divergence of 15 species (including seed plants, birds and rodents), and relate them to the evolutionary history of the naturally fragmented cloud forests. We used fossil-calibrated genealogies, coalescent-based divergence time inference, and estimates of gene flow to assess the permeability of putative barriers to gene flow. We also used the hierarchical Approximate Bayesian Computation (HABC) method implemented in the program msBayes to test simultaneous versus non-simultaneous divergence of the cloud forest lineages. Our results show shared phylogeographic breaks that correspond to the Isthmus of Tehuantepec, Los Tuxtlas, and the Chiapas Central Depression, with the Isthmus representing the most frequently shared break among taxa. However, dating analyses suggest that the phylogeographic breaks corresponding to the Isthmus occurred at different times in different taxa. Current divergence patterns are therefore consistent with the hypothesis of broad vicariance across the Isthmus of Tehuantepec derived from different mechanisms operating at different times. This study, coupled with existing data on divergence cloud forest species, indicates that the evolutionary history of contemporary cloud forest lineages is complex and often lineage-specific, and thus difficult to capture in a simple conservation strategy.

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 copperheads 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.

Complex evolution in the Neotropics: The origin and diversification of the widespread genus Leptodeira (Serpentes: Colubridae)

Lineage diversification in the Neotropics is an interesting topic in evolutionary biology but is also one of the least understood. The abiotic and biotic complexity of the region precludes generalizations that can be drawn regarding the historical evolutionary processes responsible for the diversity observed. The snake genus Leptodeira provides an excellent opportunity to investigate such processes because it spans the entire Neotropical region. In this study, we infer the phylogenetic position of Leptodeira within Dipsadinae, estimate evolutionary relationships among and within Leptodeira species, and estimate the diversification time and biogeography of the genus. Three mitochondrial gene regions were sequenced for individuals representing all the Leptodeira species and most subspecies currently recognized. Additionally, two nuclear protein-coding gene regions were sequenced for representatives of each species and several genera within the Dipsadinae. We infer that several Leptodeira species are either paraphyletic or polyphyletic as currently recognized, and that most recognized subspecies are not monophyletic lineages. Despite the taxonomic discordance with evolutionary relationships, clades appear to correspond very well to major biogeographic regions of Mexico, Central America and South America. Our results thus highlight the important role of the Miocene and Pliocene for lineage diversification in the Neotropics. Additionally, our time estimates suggest that recent intraspecific phylogeographic structure is likely the result of habitat and climatic fluctuations during the Pleistocene. Cumulatively, our inferences of lineage diversification within Leptodeira suggest a complex evolutionary scenario in the Mexican transition zone and a north to south expansion with a final colonization of the tropics in South America.

Automatic recognition of anuran species based on syllable identification

Abstract Monitoring of biological populations is well known for being a complex task that involves high operational costs, unknown reproductive intervals of the studied species, and difficult visualization of isolated individuals (due to their mimetic and cryptic capabilities). Therefore, the development of new methodologies able to measure quantities of individuals in specific biological populations without direct contact is desired. Species and individual recognition, based on acoustic analysis of their calls (Bioacoustics), is possible for many animals and has proven to be a useful tool in the study and monitoring of animal species. In this paper, an unsupervised methodology for anuran automatic identification is proposed; it is based on the use of a fuzzy classifier and Mel Frequency Cepstral Coefficients. This methodology is able to detect species not presented in the training stage, although they belong to different populations. Additionally, correlations among species of the same genus can be determined through the similarities of their calls. For testing the proposed method, two different datasets with species from the northeastern Colombia (Choco and Antioquia departments with 103 and 813 mating calls respectively) were used. In validation tests performed, accuracies between 99.38% and 100% were achieved in all species by applying the proposed methodology to both datasets. Thirteen different species of anurans in both datasets were correctly identified.

MORPHOMETRIC VARIATION AND ITS EFFECT ON REPRODUCTIVE POTENTIAL IN FEMALE COLOMBIAN SLIDER TURTLES (TRACHEMYS CALLIROSTRIS CALLIROSTRIS)

We evaluated morphometric variation and reproductive effort in Trachemys callirostris callirostris females from nine populations in the Mompos Depression of northern Colombia. Two hypotheses were evaluated: (1) phenotypes of nesting females co-vary according to geographic distances among sites and/or specific environmental conditions at each site, and (2) female body size influences their reproductive potential and offspring characteristics. Populations were shown to differ in both mean female body size and body shape, with differences more related to local habitat characteristics than to geographic distances among sites. Larger females occurred in sites with higher annual precipitation and lower hunting pressure. Female size was highly correlated with clutch size and mass. Although habitat characteristics affected female size, habitat alone was not a predictor of reproductive potential. Specific habitat conditions and hunting both influence female phenotypes in this region and thereby also influence key demographic parameters, which have clear conservation implications for this highly exploited species.

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

BackgroundSnakes 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.ResultsWe 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.ConclusionThese 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.

Dynamic nucleotide mutation gradients and control region usage in squamate reptile mitochondrial genomes

Gradients of nucleotide bias and substitution rates occur in vertebrate mitochondrial genomes due to the asymmetric nature of the replication process. The evolution of these gradients has previously been studied in detail in primates, but not in other vertebrate groups. From the primate study, the strengths of these gradients are known to evolve in ways that can substantially alter the substitution process, but it is unclear how rapidly they evolve over evolutionary time or how different they may be in different lineages or groups of vertebrates. Given the importance of mitochondrial genomes in phylogenetics and molecular evolutionary research, a better understanding of how asymmetric mitochondrial substitution gradients evolve would contribute key insights into how this gradient evolution may mislead evolutionary inferences, and how it may also be incorporated into new evolutionary models. Most snake mitochondrial genomes have an additional interesting feature, 2 nearly identical control regions, which vary among different species in the extent that they are used as origins of replication. Given the expanded sampling of complete snake genomes currently available, together with 2 additional snakes sequenced in this study, we reexamined gradient strength and CR usage in alethinophidian snakes as well as several lizards that possess dual CRs. Our results suggest that nucleotide substitution gradients (and corresponding nucleotide bias) and CR usage is highly labile over the ∼200 m.y. of squamate evolution, and demonstrates greater overall variability than previously shown in primates. The evidence for the existence of such gradients, and their ability to evolve rapidly and converge among unrelated species suggests that gradient dynamics could easily mislead phylogenetic and molecular evolutionary inferences, and argues strongly that these dynamics should be incorporated into phylogenetic models.

Fat frogs, mobile genes: unexpected phylogeographic patterns for the ornate chorus frog (Pseudacris ornata)

The southeastern coastal plain of the United States is a region marked by extraordinary phylogeographic congruence that is frequently attributed to the changing sea levels that occurred during the glacial‐interglacial cycles of the Pleistocene epoch. A phylogeographic break corresponding to the Apalachicola River has been suggested for many species studied to date that are endemic to this region. Here, we used this pattern of phylogeographic congruence to develop and test explicit hypotheses about the genetic structure in the ornate chorus frog (Pseudacris ornata). Using 1299 bp of mtDNA sequence and seven nuclear microsatellite markers in 13 natural populations of P. ornata, we found three clades corresponding to geographically distinct regions; one spans the Apalachicola River (Southern Clade), one encompasses Georgia and South Carolina (Central Clade) and a third comprises more northerly individuals (Northern Clade). However, it does not appear that typical phylogeographic barriers demarcate these clades. Instead, isolation by distance across the range of the entire species explained the pattern of genetic variation that we observed. We propose that P. ornata was historically widespread in the southeastern United States, and that a balance between genetic drift and migration was the root of the genetic divergence among populations. Additionally, we investigated fine‐scale patterns of genetic structure and found the spatial scale at which there was significant genetic structure varied among the regions studied. Furthermore, we discuss our results in light of other phylogeographic studies of southeastern coastal plain organisms and in relation to amphibian conservation and management.