Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

Identification of an integrated stress and growth response signaling switch that directs vertebrate intestinal regeneration.

BACKGROUND: Snakes exhibit extreme intestinal regeneration following months-long fasts that involves unparalleled increases in metabolism, function, and tissue growth, but the specific molecular control of this process is unknown. Understanding the mechanisms that coordinate these regenerative phenotypes provides valuable opportunities to understand critical pathways that may control vertebrate regeneration and novel perspectives on vertebrate regenerative capacities. RESULTS: Here, we integrate a comprehensive set of phenotypic, transcriptomic, proteomic, and phosphoproteomic data from boa constrictors to identify the mechanisms that orchestrate shifts in metabolism, nutrient uptake, and cellular stress to direct phases of the regenerative response. We identify specific temporal patterns of metabolic, stress response, and growth pathway activation that direct regeneration and provide evidence for multiple key central regulatory molecules kinases that integrate these signals, including major conserved pathways like mTOR signaling and the unfolded protein response. CONCLUSION: Collectively, our results identify a novel switch-like role of stress responses in intestinal regeneration that forms a primary regulatory hub facilitating organ regeneration and could point to potential pathways to understand regenerative capacity in vertebrates.

Population genomic analyses of schistosome parasites highlight critical challenges facing endgame elimination efforts.

Schistosomiasis persists in Asian regions despite aggressive elimination measures. To identify factors enabling continued parasite transmission, we performed reduced representation genome sequencing on Schistosoma japonicum miracidia collected across multiple years from transmission hotspots in Sichuan, China. We discovered strong geographic structure, suggesting that local, rather than imported, reservoirs are key sources of persistent infections in the region. At the village level, parasites collected after referral for praziquantel treatment are closely related to local pre-treatment populations. Schistosomes within villages are also highly related, suggesting that only a few parasites from a limited number of hosts drive re-infection. The close familial relationships among miracidia from different human hosts also implicate short transmission routes among humans. At the individual host level, genetic evidence indicates that multiple humans retained infections following referral for treatment. Our findings suggest that end-game schistosomiasis control measures should focus on completely extirpating local parasite reservoirs and confirming successful treatment of infected human hosts.

Patterns of relatedness and genetic diversity inferred from whole genome sequencing of archival blood fluke miracidia (Schistosoma japonicum).

Genomic approaches hold great promise for resolving unanswered questions about transmission patterns and responses to control efforts for schistosomiasis and other neglected tropical diseases. However, the cost of generating genomic data and the challenges associated with obtaining sufficient DNA from individual schistosome larvae (miracidia) from mammalian hosts have limited the application of genomic data for studying schistosomes and other complex macroparasites. Here, we demonstrate the feasibility of utilizing whole genome amplification and sequencing (WGS) to analyze individual archival miracidia. As an example, we sequenced whole genomes of 22 miracidia from 11 human hosts representing two villages in rural Sichuan, China, and used these data to evaluate patterns of relatedness and genetic diversity. We also down-sampled our dataset to test how lower coverage sequencing could increase the cost effectiveness of WGS while maintaining power to accurately infer relatedness. Collectively, our results illustrate that population-level WGS datasets are attainable for individual miracidia and represent a powerful tool for ultimately providing insight into overall genetic diversity, parasite relatedness, and transmission patterns for better design and evaluation of disease control efforts.

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.

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.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome.

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

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.

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.