Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Complexity in the timing of the first postnatal ecdysis in snakes.

Lepidosaurian reptiles, particularly snakes, periodically shed the outer epidermal layers of their skin (ecdysis) to restore or enhance vital functions such as regulating water and gaseous exchange, growth, and protection against insult, infection or physical injury. Although many studies have focused on the nature and mechanisms of skin shedding, little attention has been paid to the timing of the first ecdysis in neonates following birth or hatching. A recent study investigated patterns of the time to first postnatal ecdysis in snakes based on a large dataset taken from the literature. The analysis demonstrated patterns in the time to first postnatal ecdysis related to phylogeny as well as several life history traits. While this assessment provides important advances in our knowledge of this topic, data on known biophysical drivers of ecdysis - temperature and humidity - were largely unavailable and were not evaluated. The first postnatal ecdysis of neonatal snakes can be viewed as an adaptive adjustment to the transition from the aqueous environment of the embryo to the aerial environment of the newborn. Hence, the timing of the first postnatal ecdysis is logically influenced by the aerial environment into which a newborn snake or hatchling finds itself. Therefore, in this Commentary, we first emphasize the putative plasticity of ecdysis with respect to epidermal lipids that structure the water permeability barrier and are established or renewed during ecdysis to reduce transepidermal evaporative water loss. We then discuss the likely importance of biophysical variables as influential covariates that need future investigation as potential co-determinants of the timing of first postnatal ecdysis.

Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity.

Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and Phrynocephalus vlangalii) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous-viviparous pairs of lizards in different genera (Phrynocephalus, Eremias, Scincella, and Sphenomorphus) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity-viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.

Intestinal upregulation and specific dynamic action in snakes - Implications for the 'pay before pumping' hypothesis.

Animals which feed infrequently and on large prey, like many snake species, are characterized by a high magnitude of gut upregulation upon ingesting a meal. The intensity of intestinal upregulation was hypothesized to be proportional to the time and energy required for food processing (Specific-Dynamic-Action; SDA); hence, a positive correlation between the scope of intestinal growth and SDA response can be deduced. Such a correlation would support the so far not well established link between the intestinal and metabolic consequences of digestion. In this study I tested this prediction using an interspecific dataset on snakes gleaned from published sources. I found that SDAduration and SDAscope were positively correlated with post-feeding factorial increase in small intestine mass, but not with microvillar elongation. This indicates that a wide range of whole intestine remodelling (up- but potentially also downregulation) may temporarily prolong meal processing and that a greater magnitude of intestinal growth requires a stronger metabolic elevation. However, these effects do not seem large enough to drive the variation in the entire energetic costs of digestion, because SDAexpenditure was not affected either by intestinal or microvillar growth. I therefore propose that intestinal upregulation elicits non-negligible costs, but that these costs are a fairly small component of the whole SDAexpenditure.

Sex-specific growth is mirrored in feeding rate but not moulting frequency in a sexually dimorphic snake.

Sexual size dimorphism (SSD), commonly observed in snakes, may arise from a different growth rate between the sexes. This indicates a sex-specific resource intake that is in fact observable in free-living snakes. It is not so well known whether the sexes can express differential feeding rates under conditions unconstrained by spatial accessibility, competition, etc. Here, I studied sex-specific variation in growth, its correlate-moulting frequency, and feeding rate in a captive group of sexually dimorphic banded water snakes (Nerodia fasciata) with access to food unconstrained by predation, competition or space. I showed that the sexes did indeed differ in relative mass growth in that females grew faster than males (p = 0.02), but such differences were not apparent in the moulting rate (p = 0.19). Such differential growth was mirrored in the sex-specific feeding rate, with females ingesting a larger number of meals than males (p = 0.004). Such variation in feeding rate may be governed by an individual's energy expenditure and can be interpreted as a behavioural tendency that contributes to SSD development, independently of other behavioural characteristics. Sex-specific resource demands may drive the differential effects of increasing resource scarcity on both sexes.

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.

Optimal husbandry of hatchling Eastern Indigo Snakes (Drymarchon couperi) during a captive head-start program.

Optimal husbandry techniques are desirable for any headstart program, but frequently are unknown for rare species. Here we describe key reproductive variables and determine optimal incubation temperature and diet diversity for Eastern Indigo Snakes (Drymarchon couperi) grown in laboratory settings. Optimal incubation temperature was estimated from two variables dependent on temperature, shell dimpling, a surrogate for death from fungal infection, and deviation of an egg from an ovoid shape, a surrogate for death from developmental anomalies. Based on these relationships and size at hatching we determined optimal incubation temperature to be 26°C. Additionally, we used incubation data to assess the effect of temperature on duration of incubation and size of hatchlings. We also examined hatchling diets necessary to achieve optimal growth over a 21-month period. These snakes exhibited a positive linear relationship between total mass eaten and growth rate, when individuals were fed less than 1711 g of prey, and displayed constant growth for individuals exceeding 1711 g of prey. Similarly, growth rate increased linearly with increasing diet diversity up to a moderately diverse diet, followed by constant growth for higher levels of diet diversity. Of the two components of diet diversity, diet evenness played a stronger role than diet richness in explaining variance in hatchling growth. These patterns document that our goal of satiating snakes was achieved for some individuals but not others and that diets in which total grams consumed over the first 21 months of life is distributed equivalently among at least three prey genera yielded the fastest growth rates for hatchling snakes.

Three new endemic species of Epictia Gray, 1845 (Serpentes: Leptotyphlopidae) from the dry forest of northwestern Peru.

Three new blind snake species of the genus Epictia are described based on material collected in the Peruvian Regions Amazonas, Cajamarca and La Libertad. All three species are well differentiated from all congeners based on characteristics of their morphology and coloration. They share 10 scale rows around the middle of the tail and possess two supralabials with the anterior one in broad contact with the supraocular. Epictia septemlineata sp. nov. has 16 subcaudal scales, 257 mid-dorsal scale rows, a yellowish-white rostral, and a black terminal spine. Epictia vanwallachi sp. nov. exhibits 16 subcaudals, 188 mid-dorsal scale rows, a grayish-brown rostral, and a yellow terminal spine. Epictia antoniogarciai sp. nov. features 14-18 subcaudals, 195-208 mid-dorsal scale rows, a bright yellow or yellowish-white rostral, and the terminal spine and terminal portion of the tail yellow. All three species were collected in the interandean dry forest valleys of the Marañón River and its tributaries. This region is an area of endemism and warrants further attention from systematic and conservation biologists.

Fidelity of northern pine snakes (Pituophis m. melanoleucus) to natural and artificial hibernation sites in the New Jersey Pine Barrens.

Environmental managers require information on whether human-made hibernacula are used by rare snakes before constructing large numbers of them as mitigation measures. Fidelity of northern pine snakes (Pituophis m. melanoleucus) was examined in a 6-year study in the New Jersey Pine Barrens to determine whether they used natural and artificial hibernacula equally. Pine snakes used both artificial (human-made) and natural (snake-adapted) hibernacula. Most natural hibernacula were in abandoned burrows of large mammals. Occupancy rates were similar between natural and artificial hibernacula. Only 6 of 27 radio-tracked snakes did not shift hibernacula between years, whereas 78% shifted sites at least once, and fidelity from one year to the next was 42%. For snakes that switched hibernacula (n = 21), one switched among artificial hibernacula, 14 (65%) switched among natural hibernacula, and 6 (29%) switched from artificial to natural hibernacula. Data indicate that most pine snakes switch among hibernacula, mainly selecting natural hibernacula, suggesting that artificial dens are used, but protecting natural hibernacula should be a higher conservation priority.

A new blind snake (Typhlopidae) from northeastern Cuba.

Typhlops leptolepis sp. nov. is a new blind snake from Holguín Province, northeastern region of Cuba. The new species is characterized by its small and thin body size, as well as a narrow rostral in dorsal view, longer than broad, with parallel to curved sides, tapering toward anterior tip. The preocular is in contact with third supralabial only. It has 20 anterior scale rows reducing to 18 posteriorly at around midbody or posterior to midbody, moderate to high middorsal scale counts (250-308), and a peculiar coloration (head, neck, and tail whitish spotted in ventral view). Based on its morphology, the new species can be placed within the Typhlops lumbricalis species group and a key to the species belonging to this group is presented.

Fitness costs may explain the post-colonisation erosion of phenotypic plasticity.

Many organisms can adjust their phenotypes to match local environmental conditions via shifts in developmental trajectories, rather than relying on changes in gene frequencies wrought by natural selection. Adaptive developmental plasticity confers obvious benefits in terms of rapid response and higher mean fitness, so why is it not more common? Plausibly, adaptive plasticity also confers a cost; reshaping the phenotype takes time and energy, so that canalised control of trait values enhances fitness if the optimal phenotype remains the same from one generation to the next. Although this idea is central to interpreting the fitness consequences of adaptive plasticity, empirical data on costs of plasticity are scarce. In Australian tiger snakes, larger relative head size enhances maximal ingestible prey size on islands containing large prey. The trait arises via adaptive plasticity in snake populations on newly colonised islands but becomes genetically canalised on islands where snakes have been present for much longer periods. We experimentally manipulated relative head size in captive neonatal snakes to quantify the costs of adaptive plasticity. Although small-headed snakes were able to increase their head sizes when offered large prey, the delay in doing so, and their inability to consume large prey at the outset, significantly reduced their growth rates relative to conspecifics with larger heads at the beginning of the experiment. This study describes a proximate cause to the post-colonisation erosion of developmental plasticity recorded in tiger snake populations.

Cohort variation in offspring growth and survival: prenatal and postnatal factors in a late-maturing viviparous snake.

1. Recruitment to adulthood plays an important role in the population dynamics of late-maturing organisms as it is usually variable. Compared to birds and mammals, few studies assessing the contributions to this variation of environmental factors, offspring traits and maternal traits have been carried out for late-maturing snakes. 2. Cohort variation in recruitment through offspring growth and survival in the meadow viper (Vipera ursinii ursinii) was evaluated from 13 years of mark-recapture data collected at Mont Ventoux, France. In this species, females are mature at the age of 4-6 years and adult survival and fecundity rates are high and constant over time. 3. Offspring were difficult to catch during the first 3 years of their lives, but their mean annual probability of survival was reasonably high (0.48 +/- 0.11 SE). Mass and body condition at birth (mass residuals) varied significantly between years, decreased with litter size, and increased with maternal length. 4. Cohorts of offspring in better condition at birth grew faster, but offspring growth was not affected by sex, habitat or maternal traits. 5. Survival varied considerably between birth cohorts, some cohorts having a high-survival rate and others having essentially no survivors. No difference in mass or body condition at birth was found between cohorts with 'no survival' and 'good survival'. However, offspring survival in cohorts with good survival was positively correlated with mass at birth and negatively correlated with body condition at birth. 6. Thus, variation in offspring performance was influenced by direct environmental effects on survival and indirect environmental effects on growth, mediated by body condition at birth. Effects of maternal traits were entirely channelled through offspring traits.

Evolutionary biology: adaptive developmental plasticity in snakes.

The morphology of organisms is generally well matched to their environment, presumably because expression of their genes is tailored either at the population or the individual level to suit local conditions: for example, snake populations that persistently encounter large prey may accumulate gene mutations that specify a large head size, or head growth may be increased in individual snakes to meet local demands (adaptive developmental plasticity). Here we test the relative contributions of genetics and environment to the jaw sizes of two tiger snake populations: one that consumes small prey on the mainland, and an island population that relies on larger prey and has a larger jaw size. Although the idea of adaptive plasticity in response to environmental pressures is controversial, we find that both factors influence the difference in jaw size between the two populations, and the influence of developmental plasticity is greater in the island population.

Combining genetic and demographic monitoring better informs conservation of an endangered urban snake.

Conversion and fragmentation of wildlife habitat often leads to smaller and isolated populations and can reduce a species' ability to disperse across the landscape. As a consequence, genetic drift can quickly lower genetic variation and increase vulnerability to extirpation. For species of conservation concern, quantification of population size and connectivity can clarify the influence of genetic drift in local populations and provides important information for conservation management and recovery strategies. Here, we used genome-wide single nucleotide polymorphism (SNP) data and capture-mark-recapture methods to evaluate the genetic diversity and demography within seven focal sites of the endangered San Francisco gartersnake (Thamnophis sirtalis tetrataenia), a species affected by alteration and isolation of wetland habitats throughout its distribution. The primary goals were to determine the population structure and degree of genetic isolation among T. s. tetrataenia populations and estimate effective size and population abundance within sites to better understand the present and future importance of genetic drift. We also used temporally sampled datasets to examine the magnitude of genetic change over time. We found moderate population genetic structure throughout the San Francisco Peninsula that partitions sites into northern and southern regional clusters. Point estimates of both effective size and population abundance were generally small (≤ 100) for a majority of the sites, and estimates were particularly low in the northern populations. Genetic analyses of temporal datasets indicated an increase in genetic differentiation, especially for the most geographically isolated sites, and decreased genetic diversity over time in at least one site (Pacifica). Our results suggest that drift-mediated processes as a function of small population size and reduced connectivity from neighboring populations may decrease diversity and increase differentiation. Improving genetic diversity and connectivity among T. s. tetrataenia populations could promote persistence of this endangered snake.

Cross-Species BAC Mapping Highlights Conservation of Chromosome Synteny across Dragon Lizards (Squamata: Agamidae).

Dragon lizards (Squamata: Agamidae) comprise about 520 species in six subfamilies distributed across Asia, Australasia and Africa. Only five species are known to have sex chromosomes. All of them possess ZZ/ZW sex chromosomes, which are microchromosomes in four species from the subfamily Amphibolurinae, but much larger in Phrynocephalus vlangalii from the subfamily Agaminae. In most previous studies of these sex chromosomes, the focus has been on Australian species from the subfamily Amphibolurinae, but only the sex chromosomes of the Australian central bearded dragon (Pogona vitticeps) are well-characterized cytogenetically. To determine the level of synteny of the sex chromosomes of P. vitticeps across agamid subfamilies, we performed cross-species two-colour FISH using two bacterial artificial chromosome (BAC) clones from the pseudo-autosomal regions of P. vitticeps. We mapped these two BACs across representative species from all six subfamilies as well as two species of chameleons, the sister group to agamids. We found that one of these BAC sequences is conserved in macrochromosomes and the other in microchromosomes across the agamid lineages. However, within the Amphibolurinae, there is evidence of multiple chromosomal rearrangements with one of the BACs mapping to the second-largest chromosome pair and to the microchromosomes in multiple species including the sex chromosomes of P. vitticeps. Intriguingly, no hybridization signal was observed in chameleons for either of these BACs, suggesting a likely agamid origin of these sequences. Our study shows lineage-specific evolution of sequences/syntenic blocks and successive rearrangements and reveals a complex history of sequences leading to their association with important biological processes such as the evolution of sex chromosomes and sex determination.

Evolutionary Insights of the ZW Sex Chromosomesin Snakes: A New Chapter Added by the AmazonianPuffing Snakes of the Genus Spilotes.

Amazonian puffing snakes (Spilotes; Colubridae) are snakes widely distributed in the Neotropical region. However, chromosomal data are scarce in this group and, when available, are only limited to karyotype description using conventional staining. In this paper, we focused on the process of karyotype evolution and trends for sex chromosomes in two Amazonian Puffer Snakes (S.pulllatus and S.sulphureus). We performed an extensive karyotype characterization using conventional and molecular cytogenetic approaches. The karyotype of S.sulphureus (presented here for the first time) exhibits a 2n = 36, similar to that previously described in S.pullatus. Both species have highly differentiated ZZ/ZW sex chromosomes, where the W chromosome is highly heterochromatic in S.pullatus but euchromatic in S.sulphureus. Both W chromosomes are homologous between these species as revealed by cross-species comparative genomic hybridization, even with heterogeneous distributions of several repetitive sequences across their genomes, including on the Z and on the W chromosomes. Our study provides evidence that W chromosomes in these two species have shared ancestry.

Estimating survival rates of uncatchable animals: the myth of high juvenile mortality in reptiles.

Survival rates of juvenile reptiles are critical population parameters but are difficult to obtain through mark-recapture programs because these small, secretive animals are rarely caught. This scarcity has encouraged speculation that survival rates of juveniles are very low, and we test this prediction by estimating juvenile survival rates indirectly. A simple mathematical model calculates the annual juvenile survival rate needed to maintain a stable population size, using published data on adult survival rates, reproductive output, and ages at maturity in 109 reptile populations encompassing 57 species. Counter to prediction, estimated juvenile survival rates were relatively high (on average, only about 13% less than those of conspecific adults) and highly correlated with adult survival rates. Overall, survival rates during both juvenile and adult life were higher in turtles than in snakes, and higher in snakes than in lizards. As predicted from life history theory, rates of juvenile survival were higher in species that produce large offspring, and higher in viviparous squamates than in oviparous species. Our analyses challenge the widely held belief that juvenile reptiles have low rates of annual survival and suggest instead that sampling problems and the elusive biology of juvenile reptiles have misled researchers in this respect.

A mid-Cretaceous embryonic-to-neonate snake in amber from Myanmar.

We present the first known fossilized snake embryo/neonate preserved in early Late Cretaceous (Early Cenomanian) amber from Myanmar, which at the time, was an island arc including terranes from Austral Gondwana. This unique and very tiny snake fossil is an articulated postcranial skeleton, which includes posterior precloacal, cloacal, and caudal vertebrae, and details of squamation and body shape; a second specimen preserves a fragment of shed skin interpreted as a snake. Important details of skeletal ontogeny, including the stage at which snake zygosphene-zygantral joints began to form along with the neural arch lamina, are preserved. The vertebrae show similarities to those of fossil Gondwanan snakes, suggesting a dispersal route of Gondwanan faunas to Laurasia. Finally, the new species is the first Mesozoic snake to be found in a forested environment, indicating greater ecological diversity among early snakes than previously thought.