The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Ontogenetic change in the venom composition of one Mexican black-tailed rattlesnake (Crotalus molossus nigrescens) from Durango, Mexico

To corroborate the ontogenetic shift in the venom composition of the Mexican Black-tailed Rattlesnake (Crotalus molossus nigrescens) previously reported through the census approach, we evaluated the shift in the protein profile, lethality, and proteolytic and phospholipase activities of four venom samples obtained in 2015, 2018, 2019, and 2021 from one C. m. nigrescens individual (CMN06) collected in Durango, Mexico. We demonstrated that the venom of C. m. nigrescens changed from a myotoxin-rich venom to a phospholipase A2 and snake venom metalloproteinase-rich venom. Additionally, the proteolytic and phospholipase activities increased with age, but the lethality decreased approximately three times.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

Large-scale snake genome analyses provide insights into vertebrate development

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Venom phenotype conservation suggests integrated specialization in a lizard-eating snake

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.

Evolutionary allometry and ecological correlates of fang length evolution in vipers

Traits for prey acquisition form the phenotypic interface of predator-prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.

Intraspecific variation in male mating strategies in an African ground squirrel (Xerus inauris)

Male mating strategies respond to female availability such that variation in resources that affect spatial distribution can also alter cost–benefit tradeoffs within a population. In arid‐adapted species, rainfall alters reproduction, behavior, morphology, and population density such that populations differing in resource availability may also differ in successful reproductive strategies. Here, we compare two populations of Cape ground squirrels (Xerus inauris), a sub‐Saharan species with year‐round breeding and intense mating competition. Unlike most mammals where males resort to aggressive interactions over females, male X. inauris are tolerant of one another, relying instead on other nonaggressive pre‐ and postcopulatory strategies to determine reproductive success. Our findings suggest that differences in resource availability affect female distribution, which ultimately leads to intraspecific variation in male reproductive tactics and sexual morphology. Sperm competition, assessed by reproductive morphometrics, was more pronounced in our high resource site where females were distributed evenly across the landscape, whereas dominance seemed to be an important determinant of success in our low resource site where females were more aggregated. Both sites had similar mating intensities, and most males did not sire any offspring. However, our low resource site had a higher variance in fertilization success with fewer males siring multiple offspring compared with our high resource site where more individuals were successful. Our results lend support to resource models where variations in female spatial distribution attributed to environmental resources ultimately impact male reproductive behaviors and morphology.

De Novo Genome Assembly Highlights the Role of Lineage-Specific Duplications in the Evolution of Venom in Fea’s Viper

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A₂ (PLA₂) proteins expressed by genes that likely arose from lineage-specific PLA₂ gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

SARS-CoV-2 variants of concern Alpha and Delta show increased viral load in saliva

BACKGROUND

Higher viral loads in SARS-CoV-2 infections may be linked to more rapid spread of emerging variants of concern (VOC). Rapid detection and isolation of cases with highest viral loads, even in pre- or asymptomatic individuals, is essential for the mitigation of community outbreaks.

METHODS AND FINDINGS

In this study, we analyze Ct values from 1297 SARS-CoV-2 positive patient saliva samples collected at the Clemson University testing lab in upstate South Carolina. Samples were identified as positive using RT-qPCR, and clade information was determined via whole genome sequencing at nearby commercial labs. We also obtained patient-reported information on symptoms and exposures at the time of testing. The lowest Ct values were observed among those infected with Delta (median: 22.61, IQR: 16.72-28.51), followed by Alpha (23.93, 18.36-28.49), Gamma (24.74, 18.84-30.64), and the more historic clade 20G (25.21, 20.50-29.916). There was a statistically significant difference in Ct value between Delta and all other clades (all p.adj<0.01), as well as between Alpha and 20G (p.adj<0.05). Additionally, pre- or asymptomatic patients (n = 1093) showed the same statistical differences between Delta and all other clades (all p.adj<0.01); however, symptomatic patients (n = 167) did not show any significant differences between clades. Our weekly testing strategy ensures that cases are caught earlier in the infection cycle, often before symptoms are present, reducing this sample size in our population.

CONCLUSIONS

COVID-19 variants Alpha and Delta have substantially higher viral loads in saliva compared to more historic clades. This trend is especially observed in individuals who are pre- or asymptomatic, which provides evidence supporting higher transmissibility and more rapid spread of emerging variants. Understanding the viral load of variants spreading within a community can inform public policy and clinical decision making.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

SARS-CoV-2 variants of concern Alpha and Delta show increased viral load in saliva.. [PMID: 35194615 PMCID: PMC8863157 DOI: 10.1101/2022.02.10.22270797]

Background

Higher viral loads in SARS-CoV-2 infections may be linked to more rapid spread of emerging variants of concern (VOC). Rapid detection and isolation of cases with highest viral loads, even in pre- or asymptomatic individuals, is essential for the mitigation of community outbreaks.

Methods and Findings

In this study, we analyze Ct values from 1297 SARS-CoV-2 positive patient saliva samples collected at the Clemson University testing lab in upstate South Carolina. Samples were identified as positive using RT-qPCR, and clade information was determined via whole genome sequencing at nearby commercial labs. We also obtained patient-reported information on symptoms and exposures at the time of testing. The lowest Ct values were observed among those infected with Delta (median: 22.61, IQR: 16.72–28.51), followed by Alpha (23.93, 18.36–28.49), Gamma (24.74, 18.84–30.64), and the more historic clade 20G (25.21, 20.50–29.916). There was a statistically significant difference in Ct value between Delta and all other clades (all p.adj<0.01), as well as between Alpha and 20G (p.adj<0.05). Additionally, pre- or asymptomatic patients (n=1093) showed the same statistical differences between Delta and all other clades (all p.adj<0.01); however, symptomatic patients (n=167) did not show any significant differences between clades. Our weekly testing strategy ensures that cases are caught earlier in the infection cycle, often before symptoms are present, reducing this sample size in our population.

Conclusions

COVID-19 variants Alpha and Delta have substantially higher viral loads in saliva compared to more historic clades. This trend is especially observed in individuals who are pre- or asymptomatic, which provides evidence supporting higher transmissibility and more rapid spread of emerging variants. Understanding the viral load of variants spreading within a community can inform public policy and clinical decision making.

An integrative view of the toxic potential of Conophis lineatus (Dipsadidae: Xenodontinae), a medically relevant rear-fanged snake

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using μCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of 'harmless' non-front-fanged snakes.

Intraspecific venom variation of Mexican West Coast Rattlesnakes (Crotalus basiliscus) and its implications for antivenom production

Intraspecific variation in snake venoms has been widely documented worldwide. However, there are few studies on this subject in Mexico. Venom characterization studies provide important data used to predict clinical syndromes, to evaluate the efficacy of antivenoms and, in some cases, to improve immunogenic mixtures in the production of antivenoms. In the present work, we evaluated the intraspecific venom variation of Crotalus basiliscus, a rattlesnake of medical importance and whose venom is used in the immunization of horses to produce one of the Mexican antivenoms. Our results demonstrate that there is variation in biological and biochemical activities among adult venoms and that there is an ontogenetic change from juvenile to adult venoms. Juvenile venoms were more lethal and had higher percentages of crotamine and crotoxin, while adult venoms had higher percentages of snake venom metalloproteases (SVMPs). Additionally, we documented crotoxin-like PLA₂ variation in which specimens from Zacatecas, Sinaloa and Michoacán (except 1) lacked the neurotoxin, while the rest of the venoms had it. Finally, we evaluated the efficacy of three lots of Birmex antivenom and all three were able to neutralize the lethality of four representative venoms but were not able to neutralize crotamine. We also observed significant differences in the LD₅₀ values neutralized per vial among the different lots. Based on these results, we recommend including venoms containing crotamine in the production of antivenom for a better immunogenic mixture and to improve the homogeneity of lots.

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.

Duvernoy's Gland Transcriptomics of the Plains Black-Headed Snake, Tantilla nigriceps (Squamata, Colubridae): Unearthing the Venom of Small Rear-Fanged Snakes

The venoms of small rear-fanged snakes (RFS) remain largely unexplored, despite increased recognition of their importance in understanding venom evolution more broadly. Sequencing the transcriptome of venom-producing glands has greatly increased the ability of researchers to examine and characterize the toxin repertoire of small taxa with low venom yields. Here, we use RNA-seq to characterize the Duvernoy’s gland transcriptome of the Plains Black-headed Snake, Tantilla nigriceps, a small, semi-fossorial colubrid that feeds on a variety of potentially dangerous arthropods including centipedes and spiders. We generated transcriptomes of six individuals from three localities in order to both characterize the toxin expression of this species for the first time, and to look for initial evidence of venom variation in the species. Three toxin families—three-finger neurotoxins (3FTxs), cysteine-rich secretory proteins (CRISPs), and snake venom metalloproteinases (SVMPIIIs)—dominated the transcriptome of T. nigriceps; 3FTx themselves were the dominant toxin family in most individuals, accounting for as much as 86.4% of an individual’s toxin expression. Variation in toxin expression between individuals was also noted, with two specimens exhibiting higher relative expression of c-type lectins than any other sample (8.7–11.9% compared to <1%), and another expressed CRISPs higher than any other toxin. This study provides the first Duvernoy’s gland transcriptomes of any species of Tantilla, and one of the few transcriptomic studies of RFS not predicated on a single individual. This initial characterization demonstrates the need for further study of toxin expression variation in this species, as well as the need for further exploration of small RFS venoms.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics

MOTIVATION

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide.

RESULTS

Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Genome-wide SNPs clarify lineage diversity confused by coloration in coralsnakes of the Micrurus diastema species complex (Serpentes: Elapidae)

New world coralsnakes of the genus Micrurus are a diverse radiation of highly venomous and brightly colored snakes that range from North Carolina to Argentina. Species in this group have played central roles in developing and testing hypotheses about the evolution of mimicry and aposematism. Despite their diversity and prominence as model systems, surprisingly little is known about species boundaries and phylogenetic relationships within Micrurus, which has substantially hindered meaningful analyses of their evolutionary history. Here we use mitochondrial genes together with thousands of nuclear genomic loci obtained via ddRADseq to study the phylogenetic relationships and population genomics of a subclade of the genus Micrurus: The M. diastema species complex. Our results indicate that prior species and species-group inferences based on morphology and color pattern have grossly misguided taxonomy, and that the M. diastema complex is not monophyletic. Based on our analyses of molecular data, we infer the phylogenetic relationships among species and populations, and provide a revised taxonomy for the group. Two non-sister species-complexes with similar color patterns are recognized, the M. distans and the M. diastema complexes, the first being basal to the monadal Micrurus and the second encompassing most North American monadal taxa. We examined all 13 species, and their respective subspecies, for a total of 24 recognized taxa in the M. diastema species complex. Our analyses suggest a reduction to 10 species, with no subspecific designations warranted, to be a more likely estimate of species diversity, namely, M. apiatus, M. browni, M. diastema, M. distans, M. ephippifer, M. fulvius, M. michoacanensis, M. oliveri, M. tener, and one undescribed species.

Trait differentiation and modular toxin expression in palm-pitvipers

Background

Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A₂ (PLA₂) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris).

Results

We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA₂ complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification.

Conclusion

Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

Intraspecific sequence and gene expression variation contribute little to venom diversity in sidewinder rattlesnakes ( Crotalus cerastes)

Traits can evolve rapidly through changes in gene expression or protein-coding sequences. However, these forms of genetic variation can be correlated and changes to one can influence the other. As a result, we might expect traits lacking differential expression to preferentially evolve through changes in protein sequences or morphological adaptation. Given the lack of differential expression across the distribution of sidewinder rattlesnakes ( Crotalus cerastes), we tested this hypothesis by comparing the coding regions of genes expressed in the venom gland transcriptomes and fang morphology. We calculated Tajima's D and F_ST across four populations comparing toxin and nontoxin loci. Overall, we found little evidence of directional selection or differentiation between populations, suggesting that changes to protein sequences do not underlie the evolution of sidewinder venom or that toxins are under extremely variant selection pressures. Although low-expression toxins do not have higher sequence divergence between populations, they do have more standing variation on which selection can act. Additionally, we found significant differences in fang length among populations. The lack of differential expression and sequence divergence suggests sidewinders-given their generalist diet, moderate gene flow and environmental variation-are under stabilizing selection which functions to maintain a generalist phenotype. Overall, we demonstrate the importance of examining the relationship between gene expression and protein-coding changes to understand the evolution of complex traits.

Ontogenetic Change in the Venom of Mexican Black-Tailed Rattlesnakes (Crotalus molossus nigrescens)

Ontogenetic changes in venom composition have important ecological implications due the relevance of venom in prey acquisition and defense. Additionally, intraspecific venom variation has direct medical consequences for the treatment of snakebite. However, ontogenetic changes are not well documented in most species. The Mexican Black-tailed Rattlesnake (Crotalus molossus nigrescens) is large-bodied and broadly distributed in Mexico. To document venom variation and test for ontogenetic changes in venom composition, we obtained venom samples from twenty-seven C. m. nigrescens with different total body lengths (TBL) from eight states in Mexico. The primary components in the venom were detected by reverse-phase HPLC, western blot, and mass spectrometry. In addition, we evaluated the biochemical (proteolytic, coagulant and fibrinogenolytic activities) and biological (LD₅₀ and hemorrhagic activity) activities of the venoms. Finally, we tested for recognition and neutralization of Mexican antivenoms against venoms of juvenile and adult snakes. We detected clear ontogenetic venom variation in C. m. nigrescens. Venoms from younger snakes contained more crotamine-like myotoxins and snake venom serine proteinases than venoms from older snakes; however, an increase of snake venom metalloproteinases was detected in venoms of larger snakes. Venoms from juvenile snakes were, in general, more toxic and procoagulant than venoms from adults; however, adult venoms were more proteolytic. Most of the venoms analyzed were hemorrhagic. Importantly, Mexican antivenoms had difficulties recognizing low molecular mass proteins (<12 kDa) of venoms from both juvenile and adult snakes. The antivenoms did not neutralize the crotamine effect caused by the venom of juveniles. Thus, we suggest that Mexican antivenoms would have difficulty neutralizing some human envenomations and, therefore, it may be necessary improve the immunization mixture in Mexican antivenoms to account for low molecular mass proteins, like myotoxins.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

Genetic assessment of Staphylococcus aureus in an underreported locality: Ambulatory care clinic

BACKGROUND

Staphylococcus aureus has strong association with anthropogenic environments. This association has not been well supported by use of genetic tools. The aim of this study was to phylogenetically relate numerous isolates from three environments - NCBI samples from hospitals, a community, and a previously unexplored healthcare environment: an ambulatory care clinic (ACC).

METHODS

This study incorporated hospital samples from NCBI, a community database from the University of Central Florida (UCF), and newly added samples taken from employees of an ambulatory care clinic located at UCF. Samples were collected from nasal swabs of employees, and positive samples were cultured, extracted, and sequenced at seven MLST loci and one virulence locus (spa). MLST sequences were used in eBURST and TCS population structure analyses and all sequences were incorporated into a phylogenetic reconstruction of relationships.

RESULTS

A total of 185 samples were incorporated in this study (15 NCBI sequences from hospital infections, 29 from the ACC, and 141 from the community). In both phylogenetic and population genetics analyses, samples proved to be panmixic, with samples not segregating monophyletically based on sample origin.

CONCLUSION

Samples isolated from ambulatory care clinics are not significantly differentiated from either community or hospital samples at the representative loci chosen. These results strengthen previous conclusions that S. aureus may exhibit high genetic similarity across anthropogenic environments.

Phenotypic Variation in Mojave Rattlesnake (Crotalus scutulatus) Venom Is Driven by Four Toxin Families

Phenotypic diversity generated through altered gene expression is a primary mechanism facilitating evolutionary response in natural systems. By linking the phenotype to genotype through transcriptomics, it is possible to determine what changes are occurring at the molecular level. High phenotypic diversity has been documented in rattlesnake venom, which is under strong selection due to its role in prey acquisition and defense. Rattlesnake venom can be characterized by the presence (Type A) or absence (Type B) of a type of neurotoxic phospholipase A2 (PLA2), such as Mojave toxin, that increases venom toxicity. Mojave rattlesnakes (Crotalus scutulatus), represent this diversity as both venom types are found within this species and within a single panmictic population in the Sonoran Desert. We used comparative venom gland transcriptomics of nine specimens of C. scutulatus from this region to test whether expression differences explain diversity within and between venom types. Type A individuals expressed significantly fewer toxins than Type B individuals owing to the diversity of C-type lectins (CTLs) and snake venom metalloproteinases (SVMPs) found in Type B animals. As expected, both subunits of Mojave toxin were exclusively found in Type A individuals but we found high diversity in four additional PLA2s that was not associated with a venom type. Myotoxin a expression and toxin number variation was not associated with venom type, and myotoxin a had the highest range of expression of any toxin class. Our study represents the most comprehensive transcriptomic profile of the venom type dichotomy in rattlesnakes and C. scutulatus. Even intra-specifically, Mojave rattlesnakes showcase the diversity of snake venoms and illustrate that variation within venom types blurs the distinction of the venom dichotomy.

Phylogeography of the Central American lancehead Bothrops asper (SERPENTES: VIPERIDAE)

The uplift and final connection of the Central American land bridge is considered the major event that allowed biotic exchange between vertebrate lineages of northern and southern origin in the New World. However, given the complex tectonics that shaped Middle America, there is still substantial controversy over details of this geographical reconnection, and its role in determining biogeographic patterns in the region. Here, we examine the phylogeography of Bothrops asper, a widely distributed pitviper in Middle America and northwestern South America, in an attempt to evaluate how the final Isthmian uplift and other biogeographical boundaries in the region influenced genealogical lineage divergence in this species. We examined sequence data from two mitochondrial genes (MT-CYB and MT-ND4) from 111 specimens of B. asper, representing 70 localities throughout the species’ distribution. We reconstructed phylogeographic patterns using maximum likelihood and Bayesian methods and estimated divergence time using the Bayesian relaxed clock method. Within the nominal species, an early split led to two divergent lineages of B. asper: one includes five phylogroups distributed in Caribbean Middle America and southwestern Ecuador, and the other comprises five other groups scattered in the Pacific slope of Isthmian Central America and northwestern South America. Our results provide evidence of a complex transition that involves at least two dispersal events into Middle America during the final closure of the Isthmus.

A cryptic palm-pitviper species (Squamata: Viperidae: Bothriechis) from the Costa Rican highlands, with notes on the variation within B. nigroviridis

Middle America is one of the most biodiverse regions in the world, harboring an exceptional number of rare and endemic species. This is especially true of Middle American cloud forests, where montane specialists occupy restricted, high-elevation ranges making them attractive candidates for investigating historical biogeography and speciation. One such highland-restricted species, the black speckled palm-pitviper (Bothriechis nigroviridis), occupies the Central, Tilarán, and Talamanca Cordilleras in Costa Rica and Panama. In this study, we investigate the genetic and morphological variation among populations of B. nigroviridis by inferring a multilocus phylogeny (21 individuals) and analyzing meristic scale characters with a principal component analysis (64 individuals). We find B. nigroviridis sensu stricto to be composed of two deeply divergent lineages, one with a restricted range in the northern and central Cordillera Talamanca and the other ranging throughout the Central, Tilarán, and Talamanca Cordilleras. Furthermore, these two lineages are morphologically distinct, with previously unrecognized differences in several characters allowing us to name and diagnose a new species B. nubestris sp. nov. We also examine the genetic and morphological variation within B. nigroviridis and discuss biogeographic hypotheses that may have led to the diversification of Bothriechis lineages.

Longitudinal genetic analyses of Staphylococcus aureus nasal carriage dynamics in a diverse population

Background

Staphylococcus aureus (SA) nasal colonization plays a critical role in the pathogenesis of staphylococcal infections and SA eradication from the nares has proven to be effective in reducing endogenous infections. To understand SA nasal colonization and its relation with consequent disease, assessment of nasal carriage dynamics and genotypic diversity among a diverse population is a necessity.

Results

We have performed extensive longitudinal monitoring of SA nasal carriage isolates in 109 healthy individuals over a period of up to three years. Longitudinal sampling revealed that 24% of the individuals were persistent SA nasal carriers while 32% were intermittent. To assess the genetic relatedness between different SA isolates within our cohort, multi locus sequence typing (MLST) was performed. MLST revealed that not only were strains colonizing intermittent and persistent nasal carriers genetically similar, belonging to the same clonal complexes, but strain changes within the same host were also observed over time for both types of carriers. More highly discriminating genetic analyses using the hypervariable regions of staphylococcal protein A and clumping factor B virulence genes revealed no preferential colonization of specific SA strains in persistent or intermittent carriers. Moreover, we observed that a subset of persistent and intermittent carriers retained clinically relevant community-acquired methicillin-resistant SA (CA-MRSA) strains in their nares over time.

Conclusions

The findings of this study provides added perspective on the nasal carriage dynamics between strains colonizing persistent and intermittent carriers; an area currently in need of assessment given that persistent carriers are at greater risk of autoinfection than intermittent carriers.

Phylogenetic relationships among Staphylococcus species and refinement of cluster groups based on multilocus data

Background

Estimates of relationships among Staphylococcus species have been hampered by poor and inconsistent resolution of phylogenies based largely on single gene analyses incorporating only a limited taxon sample. As such, the evolutionary relationships and hierarchical classification schemes among species have not been confidently established. Here, we address these points through analyses of DNA sequence data from multiple loci (16S rRNA gene, dnaJ, rpoB, and tuf gene fragments) using multiple Bayesian and maximum likelihood phylogenetic approaches that incorporate nearly all recognized Staphylococcus taxa.

Results

We estimated the phylogeny of fifty-seven Staphylococcus taxa using partitioned-model Bayesian and maximum likelihood analysis, as well as Bayesian gene-tree species-tree methods. Regardless of methodology, we found broad agreement among methods that the current cluster groups require revision, although there was some disagreement among methods in resolution of higher order relationships. Based on our phylogenetic estimates, we propose a refined classification for Staphylococcus with species being classified into 15 cluster groups (based on molecular data) that adhere to six species groups (based on phenotypic properties).

Conclusions

Our findings are in general agreement with gene tree-based reports of the staphylococcal phylogeny, although we identify multiple previously unreported relationships among species. Our results support the general importance of such multilocus assessments as a standard in microbial studies to more robustly infer relationships among recognized and newly discovered lineages.

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

Evolutionary analyses of Staphylococcus aureus identify genetic relationships between nasal carriage and clinical isolates

Nasal carriage of Staphylococcus aureus has long been hypothesized to be a major vector for the transmission of virulent strains throughout the community. To address this hypothesis, we have analyzed the relatedness between a cohort of nasal carriage strains and clinical isolates to understand better the genetic conformity therein. To assess the relatedness between nasal carriage and clinical isolates of S. aureus, a genetic association study was conducted using multilocus sequence typing (MLST) and typing of the hypervariable regions of clumping factor and fibronectin binding protein genes. At all loci analyzed, genetic associations between both nasal carriage and clinical isolates were observed. Computational analyses of MLST data indicate that nasal carriage and clinical isolates belong to the same genetic clusters (clades), despite differences in sequence type assignments. Genetic analyses of the hypervariable regions from the clumping factor and fibronectin binding protein genes revealed that not only do clinically relevant strains belong to identical genetic lineages as the nasal carriage isolates within our cohort, but they also exhibit 100% sequence similarity within these regions. The findings of this report indicate that strains of S. aureus being carried asymptomatically throughout the community via nasal colonization are genetically related to those responsible for high levels of morbidity and mortality.

Variability of antarctic sea ice: and changes in carbon dioxide

A definitive long-term decrease in the extent of antarctic sea ice is not detectable from 9 years (1973 to 1981) of year-round satellite observations and limited prior data. Regional interannual variability is large, with sea ice decreasing in some regions while increasing in others. A significant decrease in overall ice extent during the mid-1970's, previously suggested to reflect warming induced by carbon dioxide, has not been maintained. In particular, the extent of ice in the Weddell Sea region has rebounded after a large decrease concurrent with a major oceanographic anomaly, the Weddell polynya. Over the 9 years, the trends are nearly the same in all seasons, but for periods of 3 to 5 years, greater winter ice maxima are associated with lesser summer ice minima. The decrease of the mid-1970's was preceded by an increase in ice extent from 1966 to 1972, further indicating the presence of cyclical components of variation that obscure any long-term trends that might be caused by a warming induced by carbon dioxide.

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

Phylogeographic structure and historical demography of the western diamondback rattlesnake (Crotalus atrox): A perspective on North American desert biogeography

The western diamondback rattlesnake (Crotalus atrox) is a prominent member of North American desert and semi-arid ecosystems, and its importance extends from its impact on the region's ecology and imagery, to its medical relevance as a large deadly venomous snake. We used mtDNA sequences to identify population genetic structure and historical demographic patterns across the range of this species, and relate these to broader patterns of historical biogeography of desert and semi-arid regions of the southwestern USA and adjacent Mexico. We inferred a Late Pliocene divergence between peninsular and continental lineages of Crotalus, followed by an Early Mid Pleistocene divergence across the continental divide within C. atrox. Within desert regions (Sonoran and Chihuahuan Deserts, Southern Plains, and Tamaulipan Plain) we observed population structure indicating isolation of populations in multiple Pleistocene refugia on either side of the continental divide, which we attempt to identify. Evidence of post-glacial population growth and range expansion was inferred, particularly in populations east of the continental divide. We observed clear evidence of (probably recent) gene flow across the continental divide and secondary contact of haplotype lineages. This recent gene flow appears to be particularly strong in the West-to-East direction. Our results also suggest that Crotalus tortugensis (Tortuga Island rattlesnake) and a population of 'C. atrox' inhabiting Santa Cruz Island (in the Gulf of California) previously suggested to be an unnamed species, are in fact deeply phylogenetically nested within continental lineages of C. atrox. Accordingly, we suggest C. tortugensis and 'C. atrox' from Santa Cruz Island be placed in the synonymy of C. atrox.

Hylid frog phylogeny and sampling strategies for speciose clades

How should characters and taxa be sampled to resolve efficiently the phylogeny of ancient and highly speciose groups? We addressed this question empirically in the treefrog family Hylidae, which contains > 800 species and may be nonmonophyletic with respect to other anuran families. We sampled 81 species (54 hylids and 27 outgroups) for two mitochondrial genes (12S, ND1), two nuclear genes (POMC, c-myc), and morphology (144 characters) in an attempt to resolve higher-level relationships. We then added 117 taxa to the combined data set, many of which were sampled for only one gene (12S). Despite the relative incompleteness of the majority of taxa, the resulting trees placed all taxa in the expected higher-level clades with strong support, despite some taxa being > 90% incomplete. Furthermore, we found no relationship between the completeness of a taxon and the support (parsimony bootstrap or Bayesian posterior probabilities) for its localized placement on the tree. Separate analysis of the data set with the most taxa (12S) gives a somewhat problematic estimate of higher-level relationships, suggesting that data sets scored only for some taxa (ND1, nuclear genes, morphology) are important in determining the outcome of the combined analysis. The results show that hemiphractine hylids are not closely related to other hylids and should be recognized as a distinct family. They also show that the speciose genus Hyla is polyphyletic, but that its species can be arranged into three monophyletic genera. A new classification of hylid frogs is proposed. Several potentially misleading signals in the morphological data are discussed.

Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes)

The subfamily Crotalinae (pitvipers) contains over 190 species of venomous snakes distributed in both the Old and New World. We incorporated an extensive sampling of taxa (including 28 of 29 genera), and sequences of four mitochondrial gene fragments (2.3kb) per individual, to estimate the phylogeny of pitvipers based on maximum parsimony and Bayesian phylogenetic methods. Our Bayesian analyses incorporated complex mixed models of nucleotide evolution that allocated independent models to various partitions of the dataset within combined analyses. We compared results of unpartitioned versus partitioned Bayesian analyses to investigate how much unpartitioned (versus partitioned) models were forced to compromise estimates of model parameters, and whether complex models substantially alter phylogenetic conclusions to the extent that they appear to extract more phylogenetic signal than simple models. Our results indicate that complex models do extract more phylogenetic signal from the data. We also address how differences in phylogenetic results (e.g., bipartition posterior probabilities) obtained from simple versus complex models may be interpreted in terms of relative credibility. Our estimates of pitviper phylogeny suggest that nearly all recently proposed generic reallocations appear valid, although certain Old and New World genera (Ovophis, Trimeresurus, and Bothrops) remain poly- or paraphyletic and require further taxonomic revision. While a majority of nodes were resolved, we could not confidently estimate the basal relationships among New World genera and which lineage of Old World species is most closely related to this New World group.

Multiple major increases and decreases in mitochondrial substitution rates in the plant family Geraniaceae

BACKGROUND

Rates of synonymous nucleotide substitutions are, in general, exceptionally low in plant mitochondrial genomes, several times lower than in chloroplast genomes, 10-20 times lower than in plant nuclear genomes, and 50-100 times lower than in many animal mitochondrial genomes. Several cases of moderate variation in mitochondrial substitution rates have been reported in plants, but these mostly involve correlated changes in chloroplast and/or nuclear substitution rates and are therefore thought to reflect whole-organism forces rather than ones impinging directly on the mitochondrial mutation rate. Only a single case of extensive, mitochondrial-specific rate changes has been described, in the angiosperm genus Plantago.

RESULTS

We explored a second potential case of highly accelerated mitochondrial sequence evolution in plants. This case was first suggested by relatively poor hybridization of mitochondrial gene probes to DNA of Pelargonium hortorum (the common geranium). We found that all eight mitochondrial genes sequenced from P. hortorum are exceptionally divergent, whereas chloroplast and nuclear divergence is unexceptional in P. hortorum. Two mitochondrial genes were sequenced from a broad range of taxa of variable relatedness to P. hortorum, and absolute rates of mitochondrial synonymous substitutions were calculated on each branch of a phylogenetic tree of these taxa. We infer one major, approximately 10-fold increase in the mitochondrial synonymous substitution rate at the base of the Pelargonium family Geraniaceae, and a subsequent approximately 10-fold rate increase early in the evolution of Pelargonium. We also infer several moderate to major rate decreases following these initial rate increases, such that the mitochondrial substitution rate has returned to normally low levels in many members of the Geraniaceae. Finally, we find unusually little RNA editing of Geraniaceae mitochondrial genes, suggesting high levels of retroprocessing in their history.

CONCLUSION

The existence of major, mitochondrial-specific changes in rates of synonymous substitutions in the Geraniaceae implies major and reversible underlying changes in the mitochondrial mutation rate in this family. Together with the recent report of a similar pattern of rate heterogeneity in Plantago, these findings indicate that the mitochondrial mutation rate is a more plastic character in plants than previously realized. Many molecular factors could be responsible for these dramatic changes in the mitochondrial mutation rate, including nuclear gene mutations affecting the fidelity and efficacy of mitochondrial DNA replication and/or repair and--consistent with the lack of RNA editing--exceptionally high levels of "mutagenic" retroprocessing. That the mitochondrial mutation rate has returned to normally low levels in many Geraniaceae raises the possibility that, akin to the ephemerality of mutator strains in bacteria, selection favors a low mutation rate in plant mitochondria.