Museum Genomics.

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

Signals of differential introgression in the genome of natural hybrids of Caribbean anoles.

Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.

The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes.

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals.

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

The effects of climate and demographic history in shaping genomic variation across populations of the Desert Horned Lizard (Phrynosoma platyrhinos).

Species often experience spatial environmental heterogeneity across their range, and populations may exhibit signatures of adaptation to local environmental characteristics. Other population genetic processes, such as migration and genetic drift, can impede the effects of local adaptation. Genetic drift in particular can have a pronounced effect on population genetic structure during large-scale geographic expansions, where a series of founder effects leads to decreases in genetic variation in the direction of the expansion. Here, we explore the genetic diversity of a desert lizard that occupies a wide range of environmental conditions and that has experienced post-glacial expansion northwards along two colonization routes. Based on our analyses of a large SNP data set, we find evidence that both climate and demographic history have shaped the genetic structure of populations. Pronounced genetic differentiation was evident between populations occupying cold versus hot deserts, and we detected numerous loci with significant associations with climate. The genetic signal of founder effects, however, is still present in the genomes of the recently expanded populations, which comprise subsets of genetic variation found in the southern populations.

Discovery and characterization of New Delhi metallo-ß-lactamase-1 inhibitor peptides that potentiate meropenem-dependent killing of carbapenemase-producing Enterobacteriaceae.

OBJECTIVES: Metallo-ß-lactamases (MBLs) are an emerging class of antimicrobial resistance enzymes that degrade ß-lactam antibiotics, including last-resort carbapenems. Infections caused by carbapenemase-producing Enterobacteriaceae (CPE) are increasingly prevalent, but treatment options are limited. While several serine-dependent ß-lactamase inhibitors are formulated with commonly prescribed ß-lactams, no MBL inhibitors are currently approved for combinatorial therapies. New compounds that target MBLs to restore carbapenem activity against CPE are therefore urgently needed. Herein we identified and characterized novel synthetic peptide inhibitors that bound to and inhibited NDM-1, which is an emerging ß-lactam resistance mechanism in CPE. METHODS: We leveraged Surface Localized Antimicrobial displaY (SLAY) to identify and characterize peptides that inhibit NDM-1, which is a primary carbapenem resistance mechanism in CPE. Lead inhibitor sequences were chemically synthesized and MBCs and MICs were calculated in the presence/absence of carbapenems. Kinetic analysis with recombinant NDM-1 and select peptides tested direct binding and supported NDM-1 inhibitor mechanisms of action. Inhibitors were also tested for cytotoxicity. RESULTS: We identified approximately 1700 sequences that potentiated carbapenem-dependent killing against NDM-1 Escherichia coli. Several also enhanced meropenem-dependent killing of other CPE. Biochemical characterization of a subset indicated the peptides penetrated the bacterial periplasm and directly bound NDM-1 to inhibit enzymatic activity. Additionally, each demonstrated minimal haemolysis and cytotoxicity against mammalian cell lines. CONCLUSIONS: Our approach advances a molecular platform for antimicrobial discovery, which complements the growing need for alternative antimicrobials. We also discovered lead NDM-1 inhibitors, which serve as a starting point for further chemical optimization.

The transcriptome of the veiled chameleon (Chamaeleo calyptratus): A resource for studying the evolution and development of vertebrates.

PURPOSE: The veiled chameleon (Chamaeleo calyptratus) is an emerging model system for studying functional morphology and evolutionary developmental biology (evo-devo). Chameleons possess body plans that are highly adapted to an arboreal life style, featuring laterally compressed bodies, split hands/ft for grasping, a projectile tongue, turreted independently moving eyes, and a prehensile tail. Despite being one of the most phenotypically divergent clades of tetrapods, genomic resources for chameleons are severely lacking. METHODS: To address this lack of resources, we used RNAseq to generate 288 million raw Illumina sequence reads from four adult tissues (male and female eyes and gonads) and whole embryos at three distinct developmental stages. We used these data to assemble a largely complete de novo transcriptome consisting of only 82 952 transcripts. In addition, a majority of assembled transcripts (67%) were successfully annotated. RESULTS: We then demonstrated the utility of these data in the context of studying visual system evolution by examining the content of veiled chameleon opsin genes to show that chameleons possess all five ancestral tetrapod opsins. CONCLUSION: We present this de novo, annotated, multi-tissue transcriptome assembly for the Veiled Chameleon, Chamaeleo calyptratus, as a resource to address a range of evolutionary and developmental questions. The associated raw reads and final annotated transcriptome assembly are freely available for use on NCBI and Figshare, respectively.

Shifts in Selective Pressures on Snake Phototransduction Genes Associated with Photoreceptor Transmutation and Dim-Light Ancestry.

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole-genome and targeted capture sequencing data. We used these data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures.

Multi-species comparisons of snakes identify coordinated signalling networks underlying post-feeding intestinal regeneration.

Several snake species that feed infrequently in nature have evolved the ability to massively upregulate intestinal form and function with each meal. While fasting, these snakes downregulate intestinal form and function, and upon feeding restore intestinal structure and function through major increases in cell growth and proliferation, metabolism and upregulation of digestive function. Previous studies have identified changes in gene expression that underlie this regenerative growth of the python intestine, but the unique features that differentiate this extreme regenerative growth from non-regenerative post-feeding responses exhibited by snakes that feed more frequently remain unclear. Here, we leveraged variation in regenerative capacity across three snake species-two distantly related lineages ( Crotalus and Python) that experience regenerative growth, and one ( Nerodia) that does not-to infer molecular mechanisms underlying intestinal regeneration using transcriptomic and proteomic approaches. Using a comparative approach, we identify a suite of growth, stress response and DNA damage response signalling pathways with inferred activity specifically in regenerating species, and propose a hypothesis model of interactivity between these pathways that may drive regenerative intestinal growth in snakes.

ThetaMater: Bayesian estimation of population size parameter θ from genomic data.

Summary: We describe ThetaMater, an open source R package comprising a suite of functions for efficient and scalable Bayesian estimation of the population size parameter θ from genomic data. Availability and implementation: ThetaMater is available at GitHub (https://github.com/radamsRHA/ThetaMater). Contact: todd.castoe@uta.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Assessing the Impacts of Positive Selection on Coalescent-Based Species Tree Estimation and Species Delimitation.

The assumption of strictly neutral evolution is fundamental to the multispecies coalescent model and permits the derivation of gene tree distributions and coalescent times conditioned on a given species tree. In this study, we conduct computer simulations to explore the effects of violating this assumption in the form of species-specific positive selection when estimating species trees, species delimitations, and coalescent parameters under the model. We simulated data sets under an array of evolutionary scenarios that differ in both speciation parameters (i.e., divergence times, strength of selection) and experimental design (i.e., number of loci sampled) and incorporated species-specific positive selection occurring within branches of a species tree to identify the effects of selection on multispecies coalescent inferences. Our results highlight particular evolutionary scenarios and parameter combinations in which inferences may be more, or less, susceptible to the effects of positive selection. In some extreme cases, selection can decrease error in species delimitation and increase error in species tree estimation, yet these inferences appear to be largely robust to the effects of positive selection under many conditions likely to be encountered in empirical data sets.

GppFst: genomic posterior predictive simulations of FST and dXY for identifying outlier loci from population genomic data.

Summary: We introduce GppFst, an open source R package that generates posterior predictive distributions of FST and dx under a neutral coalescent model to identify putative targets of selection from genomic data. Availability and Implementation: GppFst is available at ( https://github.com/radamsRHA/GppFst ). Contact: todd.castoe@uta.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Plasticity and local adaptation explain lizard cold tolerance.

How does climate variation limit the range of species and what does it take for species to colonize new regions? In this issue of Molecular Ecology, Campbell-Staton et al. () address these broad questions by investigating cold tolerance adaptation in the green anole lizard (Anolis carolinensis) across a latitudinal transect. By integrating physiological data, gene expression data and acclimation experiments, the authors disentangle the mechanisms underlying cold adaptation. They first establish that cold tolerance adaptation in Anolis lizards follows the predictions of the oxygen- and capacity-limited thermal tolerance hypothesis, which states that organisms are limited by temperature thresholds at which oxygen supply cannot meet demand. They then explore the drivers of cold tolerance at a finer scale, finding evidence that northern populations are adapted to cooler thermal regimes and that both phenotypic plasticity and heritable genetic variation contribute to cold tolerance. The integration of physiological and gene expression data further highlights the varied mechanisms that drive cold tolerance adaptation in Anolis lizards, including both supply-side and demand-side adaptations that improve oxygen economy. Altogether, their work provides new insight into the physiological and genetic mechanisms underlying adaptation to new climatic niches and demonstrates that cold tolerance in northern lizard populations is achieved through the synergy of physiological plasticity and local genetic adaptation for thermal performance.

Novel ecological and climatic conditions drive rapid adaptation in invasive Florida Burmese pythons.

Invasive species provide powerful in situ experimental systems for studying evolution in response to selective pressures in novel habitats. While research has shown that phenotypic evolution can occur rapidly in nature, few examples exist of genomewide adaptation on short "ecological" timescales. Burmese pythons (Python molurus bivittatus) have become a successful and impactful invasive species in Florida over the last 30 years despite major freeze events that caused high python mortality. We sampled Florida Burmese pythons before and after a major freeze event in 2010 and found evidence for directional selection in genomic regions enriched for genes associated with thermosensation, behaviour and physiology. Several of these genes are linked to regenerative organ growth, an adaptive response that modulates organ size and function with feeding and fasting in pythons. Independent histological and functional genomic data sets provide additional layers of support for a contemporary shift in invasive Burmese python physiology. In the Florida population, a shift towards maintaining an active digestive system may be driven by the fitness benefits of maintaining higher metabolic rates and body temperature during freeze events. Our results suggest that a synergistic interaction between ecological and climatic selection pressures has driven adaptation in Florida Burmese pythons, demonstrating the often-overlooked potential of rapid adaptation to influence the success of invasive species.

Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus).

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene.

Growth and stress response mechanisms underlying post-feeding regenerative organ growth in the Burmese python.

BACKGROUND: Previous studies examining post-feeding organ regeneration in the Burmese python (Python molurus bivittatus) have identified thousands of genes that are significantly differentially regulated during this process. However, substantial gaps remain in our understanding of coherent mechanisms and specific growth pathways that underlie these rapid and extensive shifts in organ form and function. Here we addressed these gaps by comparing gene expression in the Burmese python heart, liver, kidney, and small intestine across pre- and post-feeding time points (fasted, one day post-feeding, and four days post-feeding), and by conducting detailed analyses of molecular pathways and predictions of upstream regulatory molecules across these organ systems. RESULTS: Identified enriched canonical pathways and upstream regulators indicate that while downstream transcriptional responses are fairly tissue specific, a suite of core pathways and upstream regulator molecules are shared among responsive tissues. Pathways such as mTOR signaling, PPAR/LXR/RXR signaling, and NRF2-mediated oxidative stress response are significantly differentially regulated in multiple tissues, indicative of cell growth and proliferation along with coordinated cell-protective stress responses. Upstream regulatory molecule analyses identify multiple growth factors, kinase receptors, and transmembrane receptors, both within individual organs and across separate tissues. Downstream transcription factors MYC and SREBF are induced in all tissues. CONCLUSIONS: These results suggest that largely divergent patterns of post-feeding gene regulation across tissues are mediated by a core set of higher-level signaling molecules. Consistent enrichment of the NRF2-mediated oxidative stress response indicates this pathway may be particularly important in mediating cellular stress during such extreme regenerative growth.

Hormonally Mediated Increases in Sex-Biased Gene Expression Accompany the Breakdown of Between-Sex Genetic Correlations in a Sexually Dimorphic Lizard.

The evolution of sexual dimorphism is predicted to occur through reductions in between-sex genetic correlations (rmf) for shared traits, but the physiological and genetic mechanisms that facilitate these reductions remain largely speculative. Here, we use a paternal half-sibling breeding design in captive brown anole lizards (Anolis sagrei) to show that the development of sexual size dimorphism is mirrored by the ontogenetic breakdown of rmf for body size and growth rate. Using transcriptome data from the liver (which integrates growth and metabolism), we show that sex-biased gene expression also increases dramatically between ontogenetic stages bracketing this breakdown of rmf. Ontogenetic increases in sex-biased expression are particularly evident for genes involved in growth, metabolism, and cell proliferation, suggesting that they contribute to both the development of sexual dimorphism and the breakdown of rmf. Mechanistically, we show that treatment of females with testosterone stimulates the expression of male-biased genes while inhibiting the expression of female-biased genes, thereby inducing male-like phenotypes at both organismal and transcriptomic levels. Collectively, our results suggest that sex-specific modifiers such as testosterone can orchestrate sex-biased gene expression to facilitate the phenotypic development of sexual dimorphism while simultaneously reducing genetic correlations that would otherwise constrain the independent evolution of the sexes.

Contrasting gene expression programs correspond with predator-induced phenotypic plasticity within and across generations in Daphnia.

Research has shown that a change in environmental conditions can alter the expression of traits during development (i.e., "within-generation phenotypic plasticity") as well as induce heritable phenotypic responses that persist for multiple generations (i.e., "transgenerational plasticity", TGP). It has long been assumed that shifts in gene expression are tightly linked to observed trait responses at the phenotypic level. Yet, the manner in which organisms couple within- and TGP at the molecular level is unclear. Here we tested the influence of fish predator chemical cues on patterns of gene expression within- and across generations using a clone of Daphnia ambigua that is known to exhibit strong TGP but weak within-generation plasticity. Daphnia were reared in the presence of predator cues in generation 1, and shifts in gene expression were tracked across two additional asexual experimental generations that lacked exposure to predator cues. Initial exposure to predator cues in generation 1 was linked to ~50 responsive genes, but such shifts were 3-4× larger in later generations. Differentially expressed genes included those involved in reproduction, exoskeleton structure and digestion; major shifts in expression of genes encoding ribosomal proteins were also identified. Furthermore, shifts within the first-generation and transgenerational shifts in gene expression were largely distinct in terms of the genes that were differentially expressed. Such results argue that the gene expression programmes involved in within- vs. transgeneration plasticity are fundamentally different. Our study provides new key insights into the plasticity of gene expression and how it relates to phenotypic plasticity in nature.

Expression of venom gene homologs in diverse python tissues suggests a new model for the evolution of snake venom.

Snake venom gene evolution has been studied intensively over the past several decades, yet most previous studies have lacked the context of complete snake genomes and the full context of gene expression across diverse snake tissues. We took a novel approach to studying snake venom evolution by leveraging the complete genome of the Burmese python, including information from tissue-specific patterns of gene expression. We identified the orthologs of snake venom genes in the python genome, and conducted detailed analysis of gene expression of these venom homologs to identify patterns that differ between snake venom gene families and all other genes. We found that venom gene homologs in the python are expressed in many different tissues outside of oral glands, which illustrates the pitfalls of using transcriptomic data alone to define "venom toxins." We hypothesize that the python may represent an ancestral state prior to major venom development, which is supported by our finding that the expansion of venom gene families is largely restricted to highly venomous caenophidian snakes. Therefore, the python provides insight into biases in which genes were recruited for snake venom systems. Python venom homologs are generally expressed at lower levels, have higher variance among tissues, and are expressed in fewer organs compared with all other python genes. We propose a model for the evolution of snake venoms in which venom genes are recruited preferentially from genes with particular expression profile characteristics, which facilitate a nearly neutral transition toward specialized venom system expression.

Phylogeographic and population genetic analyses reveal multiple species of Boa and independent origins of insular dwarfism.

Boa is a Neotropical genus of snakes historically recognized as monotypic despite its expansive distribution. The distinct morphological traits and color patterns exhibited by these snakes, together with the wide diversity of ecosystems they inhabit, collectively suggest that the genus may represent multiple species. Morphological variation within Boa also includes instances of dwarfism observed in multiple offshore island populations. Despite this substantial diversity, the systematics of the genus Boa has received little attention until very recently. In this study we examined the genetic structure and phylogenetic relationships of Boa populations using mitochondrial sequences and genome-wide SNP data obtained from RADseq. We analyzed these data at multiple geographic scales using a combination of phylogenetic inference (including coalescent-based species delimitation) and population genetic analyses. We identified extensive population structure across the range of the genus Boa and multiple lines of evidence for three widely-distributed clades roughly corresponding with the three primary land masses of the Western Hemisphere. We also find both mitochondrial and nuclear support for independent origins and parallel evolution of dwarfism on offshore island clusters in Belize and Cayos Cochinos Menor, Honduras.