Pheromone Gene Diversification and the Evolution of Courtship Glands in Plethodontid Salamanders.

Proteinaceous pheromones that diversify through gene duplication can result in shifts in courtship cocktails that may serve as a mechanism for reproductive isolation. The molecular evolution of pheromones has been extensively studied in salamanders, but how these genes and associated novel courtship glands have codiversified has not been evaluated. In this study we used transcriptional analyses to examine the relationship between pheromone diversification and gland type in three divergent lineages of plethodontid salamanders. Our results revealed that plethodontid salamanders express up to eight divergent Sodefrin Precursor-like Factor genes (spf, representing both alpha and beta subfamilies) along with Plethodontid Modulating Factor (pmf) and Plethodontid Receptivity Factor (prf). Expression of pheromone genes is tissue specific with pmf, prf, and some spf genes restricted to the mental gland. In contrast, the caudal gland shows strong expression of the other spf genes. We found evidence for punctuated changes in pheromone cocktail composition related to the loss of metamorphosis, and subsequent extreme reduction of the mental gland, in a paedomorphic lineage. Our study provides insight into how pheromone diversification can be partitioned into unique glands, which may lead to cocktail specificity in behavioral modules during courtship.

Sociality sculpts similar patterns of molecular evolution in two independently evolved lineages of eusocial bees.

While it is well known that the genome can affect social behavior, recent models posit that social lifestyles can, in turn, influence genome evolution. Here, we perform the most phylogenetically comprehensive comparative analysis of 16 bee genomes to date: incorporating two published and four new carpenter bee genomes (Apidae: Xylocopinae) for a first-ever genomic comparison with a monophyletic clade containing solitary through advanced eusocial taxa. We find that eusocial lineages have undergone more gene family expansions, feature more signatures of positive selection, and have higher counts of taxonomically restricted genes than solitary and weakly social lineages. Transcriptomic data reveal that caste-affiliated genes are deeply-conserved; gene regulatory and functional elements are more closely tied to social phenotype than phylogenetic lineage; and regulatory complexity increases steadily with social complexity. Overall, our study provides robust empirical evidence that social evolution can act as a major and surprisingly consistent driver of macroevolutionary genomic change.

Genetic signatures of dominance hierarchies reveal conserved cis-regulatory and brain gene expression underlying aggression in a facultatively social bee.

Agonistic interactions among individuals can result in the formation of dominance hierarches that can reinforce individual behavior and social status. Such dominance hierarches precede the establishment of reproductive dominance, division of labor and caste formation in highly social insect taxa. As such, deciphering the molecular basis of aggression is fundamental in understanding the mechanisms of social evolution. Assessing the proximate mechanisms of aggression in incipiently social bees can provide insights into the foundations of genomic mechanisms of social behavior. Here, we measured the effects of aggression on brain gene expression in the incipiently social bee, Ceratina australensis. We examine the brain transcriptomic differences between individuals who have experienced recurrent winning, losing, or a change in rank during repeated encounters. Using comparative analyses across taxa, we identify deeply conserved candidate genes, pathways, and regulatory networks for the formation of social hierarchies.

Multiple stressors produce differential transcriptomic patterns in a stream-dwelling salamander.

BACKGROUND: Global biodiversity is decreasing at an alarming rate and amphibians are at the forefront of this crisis. Understanding the factors that negatively impact amphibian populations and effectively monitoring their health are fundamental to addressing this epidemic. Plasma glucocorticoids are often used to assess stress in amphibians and other vertebrates, but these hormones can be extremely dynamic and impractical to quantify in small organisms. Transcriptomic responses to stress hormones in amphibians have been largely limited to laboratory models, and there have been few studies on vertebrates that have evaluated the impact of multiple stressors on patterns of gene expression. Here we examined the gene expression patterns in tail tissues of stream-dwelling salamanders (Eurycea tynerensis) chronically exposed to the stress hormone corticosterone under different temperature regimes. RESULTS: We found unique transcriptional signatures for chronic corticosterone exposure that were independent of temperature variation. Several of the corticosterone responsive genes are known to be involved in immune system response (LY-6E), oxidative stress (GSTM2 and TRX), and tissue repair (A2M and FX). We also found many genes to be influenced by temperature (CIRBP, HSC71, HSP40, HSP90, HSP70, ZNF593). Furthermore, the expression patterns of some genes (GSTM2, LY-6E, UMOD, ZNF593, CIRBP, HSP90) show interactive effects of temperature and corticosterone exposure, compared to each treatment alone. Through a series of experiments we also showed that stressor induced patterns of expression were largely consistent across ages, life cycle modes, and tissue regeneration. CONCLUSIONS: Outside of thermal stressors, the application of transcriptomes to monitor the health of non-human vertebrate systems has been vastly underinvestigated. Our study suggests that transcriptomic patterns harbor stressor specific signatures that can be highly informative for monitoring the diverse stressors of amphibian populations.

Conserved Genes Underlie Phenotypic Plasticity in an Incipiently Social Bee.

Despite a strong history of theoretical work on the mechanisms of social evolution, relatively little is known of the molecular genetic changes that accompany transitions from solitary to eusocial forms. Here, we provide the first genome of an incipiently social bee that shows both solitary and social colony organization in sympatry, the Australian carpenter bee Ceratina australensis. Through comparative analysis, we provide support for the role of conserved genes and cis-regulation of gene expression in the phenotypic plasticity observed in nest-sharing, a rudimentary form of sociality. Additionally, we find that these conserved genes are associated with caste differences in advanced eusocial species, suggesting these types of mechanisms could pave the molecular pathway from solitary to eusocial living. Genes associated with social nesting in this species show signatures of being deeply conserved, in contrast to previous studies in other bees showing novel and faster-evolving genes are associated with derived sociality. Our data provide support for the idea that the earliest social transitions are driven by changes in gene regulation of deeply conserved genes.

Heterochrony repolarized: a phylogenetic analysis of developmental timing in plethodontid salamanders.

BACKGROUND: Disentangling evolutionary shifts in developmental timing (heterochony) is dependent upon accurate estimates of ancestral patterns. However, many classic assessments of heterochronic patterns predate robust phylogenetic hypotheses and methods for trait reconstruction, and therefore may have been polarized with untested 'primitive' conditions. Here we revisit the heterochronic modes of development that underlie the evolution of metamorphosis, maturation, and paedomorphosis in plethodontid salamanders. We focus on the tribe Spelerpini, which is a diverse clade that exhibits tremendous variation in timing of metamorphosis and maturation, as well as multiple independent instances of larval form paedomorphosis. Based on morphology and biogeography, early investigators concluded that the most recent common ancestors of plethodontids, and also spelerpines, were large salamanders, with very long larval periods and late maturation times. This prevailing assumption influenced subsequent heterochronic assessments, which concluded that most modern spelerpines (with shorter larval periods) were derived through multiple independent accelerations in larval development. It was also concluded that most occurrences of larval form paedomorphosis in this clade resulted from progenesis (acceleration of gonadal development relative to metamorphosis). RESULTS: By reconstructing the time to metamorphosis on a molecular-based phylogeny of plethodontids, we find that ancestral spelerpines likely had relatively shorter larval periods than previously proposed. Taken together with the credibility interval from our ancestral state estimation we show that very long larval periods are likely derived decelerations, only a few lineages have undergone appreciable accelerations in metamorphic timing, and the remaining taxa have lower probabilities of being different than the ancestral condition (possibly due to stasis). Reconstructing maturation age across nodes concomitant with the evolution of larval form paedomorphosis in one large radiation does not show clear evidence of progenesis, but more likely indicates a case of neoteny (delayed metamorphosis). CONCLUSIONS: This study demonstrates cases in plethodontid salamanders where phylogenetic-based character reconstructions reject previously hypothesized ancestral life history conditions. As a result, several prior hypotheses of heterochronic evolution in this family are reversed.

Reduced effects of thyroid hormone on gene expression and metamorphosis in a paedomorphic plethodontid salamander.

It has been over a century since Gudernatsch (1912, Wilhelm Roux Arch Entwickl Mech Org 35:457-483) demonstrated that mammalian thyroid gland extracts can stimulate tadpole metamorphosis. Despite the tremendous developmental diversity of amphibians, mechanisms of metamorphosis have mostly been studied in a few model systems. This limits our understanding of the processes that influence the evolution of developmental aberrations. Here we isolated thyroid hormone receptors alpha (TRα) and beta (TRß) from Oklahoma salamanders (Eurycea tynerensis), which exhibit permanently aquatic (paedomorphic) or biphasic (metamorphic) developmental modes in different populations. We found that TRα and TRß were upregulated by thyroid hormone (T3 ) in tail tissues of larvae from metamorphic populations, but basal levels of TR expression and T3 responsiveness were reduced in larvae from paedomorphic populations. Likewise, we found that T3 treatment resulted in complete loss of larval epibranchials in larvae from metamorphic populations, but little to no epibranchial remodeling occurred in larvae from paedomorphic populations over the same duration. This is the first study to directly demonstrate reduced gene expression and metamorphic responses to T3 in a paedomorphic plethodontid compared to metamorphic conspecifics, and the first salamander system to show differential expression of thyroid hormone receptors associated with alternative developmental patterns.

Larval masquerade: a new species of paedomorphic salamander (Caudata: Plethodontidae: Eurycea) from the Ouachita Mountains of North America.

Species with truncated developmental patterns may go undetected if they resemble the juveniles of their close relatives. Herein we present an example of this phenomenon with the description of a highly divergent, relict species of stream-dwelling plethodontid salamander from the Ouachita Mountains of North America. Both mitochondrial and nuclear sequence data show that this new species is most closely related to its syntopic relative, Eurycea multiplicata. Interestingly, E. multiplicata exhibits the ancestral biphasic (metamorphic) life cycle, whereas the new species maintains an aquatic larval form throughout life (paedomorphic) and superficially resembles larval E. multiplicata. The new species is the first known paedomorphic plethodontid from the Ouachita Mountains, and the most divergent paedomorphic salamander discovered in over seventy years. This species represents an independent instance of the evolution of paedomorphosis associated with a porous streambed, which may facilitate vertical seasonal movements. This new species currently has an extremely limited known distribution and is of immediate conservation concern.

Evolution of paedomorphosis in plethodontid salamanders: ecological correlates and re-evolution of metamorphosis.

Life-history modes can profoundly impact the biology of a species, and a classic example is the dichotomy between metamorphic (biphasic) and paedomorphic (permanently aquatic) life-history strategies in salamanders. However, despite centuries of research on this system, several basic questions about the evolution of paedomorphosis in salamanders have not been addressed. Here, we use a nearly comprehensive, time-calibrated phylogeny of spelerpine plethodontids to reconstruct the evolution of paedomorphosis and to test if paedomorphosis is (1) reversible; (2) associated with living in caves; (3) associated with relatively dry climatic conditions on the surface; and (4) correlated with limited range size and geographic dispersal. We find that paedomorphosis arose multiple times in spelerpines. We also find evidence for re-evolution of metamorphosis after several million years of paedomorphosis in a lineage of Eurycea from the Edwards Plateau region of Texas. We also show for the first time using phylogenetic comparative methods that paedomorphosis is highly correlated with cave-dwelling, arid surface environments, and small geographic range sizes, providing insights into both the causes and consequences of this major life history transition.

Distribution, abundance, and genetic diversity of Clinostomum spp. metacercariae (Trematoda:Digenea) in a modified Ozark stream system.

Land-use alterations can have profound influences on faunal distributions, including host-parasite relationships. Yellow grub trematodes ( Clinostomum spp.) have complex life cycles involving 3 hosts: a snail, a fish or amphibian, and a bird. Here, we analyze the distribution, prevalence, intensity, abundance, and genetic diversity of encysting metacercariae of Clinostomum spp. in salamanders and fishes throughout an aquatic system that includes a natural Ozark stream and man-made ponds. We found Clinostomum sp. infecting permanently aquatic Oklahoma salamanders ( Eurycea tynerensis ; 56% prevalence) and larval grotto salamanders ( Eurycea spelaea ) immediately downstream from a man-made pond. However, Clinostomum sp. did not infect any salamanders in the spring that supplies this pond, or in sections farther downstream (~0.5 and 2 km). Metacercariae of Clinostomum sp. were present in ~90% of introduced largemouth bass ( Micropterus salmoides ) in the man-made pond adjunct to the stream. Morphological examination and phylogenetic analyses based on the mitochondrial gene cytochrome oxidase 1 ( Co1 ) and the nuclear ribosomal gene 18S show that fishes and salamanders at this site are primarily infected with Clinostomum marginatum . There is a relatively high degree of mitochondrial haplotype diversity in C. marginatum at this site but no consistent genetic difference between parasites in largemouth bass from the man-made pond and those in salamanders from the stream. Based on the microgeographic distribution and relationships of metacercariae of C. marginatum at this site, we hypothesize that the adjunct man-made pond has created an ecological situation that brings the cercariae of this parasite into contact with novel stream salamander hosts.

Spontaneous activity of dopaminergic retinal neurons.

Dopaminergic local circuit neurons in the retina (DA cells) show robust, spontaneous, tetrodotoxin-sensitive pacemaking. To investigate the mechanism underlying this behavior, we characterized the sodium current and a subset of the potassium currents in the cells in voltage-clamp experiments. We found that there is a persistent component of the sodium current in DA cells which activates at more depolarized potentials than the transient component of the current. The transient component was completely inactivated at -50 mV, but DA cells remained able to fire spontaneous action potentials when potassium channels were partially blocked and the membrane potential remained above -40 mV. Based on these electrophysiological data, we developed a reduced computer model that reproduced the major features of DA cells. In simulations at the physiological resting potential, the persistent component of the sodium current was both necessary and sufficient to account for spontaneous activity, and the major contribution of the transient component of the sodium current was to initiate the depolarization of the model cell during the interspike interval. When tonic inhibition was simulated by lowering the input impedance of the model cell, the transient component played a larger role.