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 A2 (PLA2) 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 PLA2 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

Transcriptomics of Tasmanian devil (Sarcophilus harrisii) ear tissue reveals homogeneous gene expression patterns across a heterogeneous landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil’s geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex

The devil is in the details: Genomics of transmissible cancers in Tasmanian devils

Cancer poses one of the greatest human health threats of our time. Fortunately, aside from a few rare cases of cancer transmission in immune-suppressed organ transplant recipients or a small number of transmission events from mother to fetus, cancers are not spread from human to human. However, transmissible cancers have been detected in vertebrate and invertebrate animals, sometimes with devastating effects. Four examples of transmissible cancers are now known: 1) canine transmissible venereal tumor (CTVT) in dogs, 2) a tumor in a laboratory population of Syrian hamsters that is no longer cultured, 3) infectious neoplasias in at least four species of bivalve mollusks, and 4) two independently derived transmissible cancers (devil facial tumor disease [DFTD]) in Tasmanian devils (Fig 1A and 1B). The etiologic agents of CTVT, the bivalve cancers, and DFTDare the transplants (allografts) of the neoplastic cells themselves, but the etiologic agent is unknown for the hamster tumor.The effects of these transmissible cancers on their respective host populations vary. CTVT is spread in dogs through sexual contact and is at least 11,000 years old, placing the timing of its origin close to that of the domestication of dogs. Although genomic analyses of the tumor suggest evasion of multiple components of the dog immune system, dogs most commonly survive and often show evidence of spontaneous tumor regression within a year of initial diagnosis. For the infectious bivalve neoplasias, which have existed for at least 40 years, population effects vary from enzootic infections with no noticeable effects on population sizes to evidence of a catastrophic population decline. In Tasmanian devils (Fig 1A), the first infectious tumor discovered (DFT1; Fig 1B) has spread across approximately 95% of the geographic range of Tasmanian devils since 1996 (Fig 1C). DFTD is almost always fatal (Fig 1B), with >90% declines in infected localities and an overall species-wide decline exceeding 80%. Transmission dynamics appear consistent with frequency dependence, with DFTD spread by biting during social interactions, resulting in predictions of extinction from standard epidemiological models. Despite these predictions, long-infected devil populations persist at reduced densities, suggesting that individual-level variability in fecundity and tumor growth rate in infected individuals are key for understanding epidemiological dynamics. Additionally, the origin of the second, independent lineage of DFTD (i.e., DFT2) within 20 years of the discovery of DFT1 suggests that transmissible cancers may be a recurring part of the Tasmanian devils' evolutionary history, without causing extinction

Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations

Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes

The genomic basis of tumor regression in Tasmanian devils (Sarcophilus harrisii)

Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers