Search

RNA sequencing (RNA-Seq) is a popular method for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can measure gene expression variation among wild-caught individuals and can yield important biological insights into organism function, sampling methods may also influence gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3’ RNA-Seq and whole mRNA-Seq. We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 minutes after euthanasia on estimated gene expression in blood, gill, and muscle. We detected 30% more genes with whole mRNA-Seq than with 3’ RNA-Seq and found that 58% of genes were significantly differently expressed between 3’ RNA-Seq and whole mRNA-Seq. Our findings indicate that 3’ RNA-Seq and whole mRNA-Seq are robust to the technical variables related to the field sampling approaches tested here with a lack of differential gene expression among sampling methods and tissue collection time after euthanasia. However, we found that gene expression varied based on which RNA-Seq library type was used on the same set of samples. Our study suggests researchers could safely rely on different fish sampling strategies in the field and save money and analyze more individuals using 3’ RNA-Seq, but should use whole mRNA-Seq when working with a species without good genomic resources, and when maximizing the number of genes identified and detecting alternative splicing are important.

Body coloration in ectotherms serves multiple biological functions, including avoiding predators, communicating with conspecific individuals, and involvement in thermoregulation. As ectotherms rely on environmental sources of heat to regulate their internal body temperature, stable melanistic body coloration or color change can be used to increase or decrease heat absorption and heat exchange with the environment. While the function of melanistic coloration for thermoregulation has been found to increase solar radiation absorption for heating in many diurnal ectotherms, research on crepuscular and nocturnal ectotherms is lacking. Since crepuscular and nocturnal ectotherms generally absorb heat from the substrate, coloration is likely under different selective pressures than in diurnal ectotherms. We tested if the proportion of dorsal melanistic body coloration is related to differences in body temperature heating and cooling rates in the crepuscular geckoEublepharis maculariusand whether changes in environmental temperature trigger color changes in this species. Temperature measurements of the geckos and of the environment were taken using infrared thermography and temperature loggers. Color data were obtained using objective photography and a newly developed custom software package. We found that body temperature reflected substrate temperatures, and that the proportion of melanistic coloration has no influence on heating or cooling rates or on color changes. These findings suggest that, inE. macularius, melanistic coloration may not be used for thermoregulation. Future research should further test the function of melanistic coloration in other crepuscular and nocturnal vertebrates to understand the evolution of melanistic pattern in animals active in low light conditions.

Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest quality squamate genomes to date for the leopard gecko,Eublepharis macularius(Eublepharidae). We compared this assembly to the previous, short-read only,E. maculariusreference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previousE. maculariusreference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified that 9 of the 19 chromosomes were assembled as single contigs, while the other 10 chromosomes were each scaffolded together from two or more contigs. We qualitatively identified that percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction previous cost estimates. This newE. maculariusreference assembly is available on NCBI at JAOPLA010000000. The genome version and its associated annotations are also available via this Figshare repositoryhttps://doi.org/10.6084/m9.figshare.20069273.

The Galapagos Archipelago is recognized as a natural laboratory for studying evolutionary processes. San Cristóbal was one of the first islands colonized by tortoises, which radiated from there across the archipelago to inhabit 10 islands. Here, we sequenced the mitochondrial control region from six historical giant tortoises from San Cristóbal (five long deceased individuals found in a cave and one found alive during an expedition in 1906) and discovered that the five from the cave are from a clade that is distinct among known Galapagos giant tortoises but closely related to the species from Española and Pinta Islands. The haplotype of the individual collected alive in 1906 is in the same clade as the haplotype in the contemporary population. To search for traces of a second lineage in the contemporary population on San Cristóbal, we closely examined the population by sequencing the mitochondrial control region for 129 individuals and genotyping 70 of these for both 21 microsatellite loci and >12,000 genome-wide single nucleotide polymorphisms [SNPs]. Only a single mitochondrial haplotype was found, with no evidence to suggest substructure based on the nuclear markers. Given the geographic and temporal proximity of the two deeply divergent mitochondrial lineages in the historical samples, they were likely sympatric, raising the possibility that the lineages coexisted. Without the museum samples, this important discovery of an additional lineage of Galapagos giant tortoise would not have been possible, underscoring the value of such collections and providing insights into the early evolution of this iconic radiation.

Identifying species with lower-than-expected cancer prevalence can help establish new models for understanding cancer resistance. Most studies of cancer prevalence have focused on mammals. Yet, other vertebrate groups vary tremendously in genetics, physiology, and ecology, which can all influence mechanisms of cancer resistance and may hold new keys to understanding cancer biology. Here, we present data on cancer prevalence in tetrapods, which includes all major vertebrate groups except fish, using necropsies from over a thousand different species. We investigated cancer prevalence within and among amphibians, birds, crocodilians, mammals, squamates, and turtles in relationship to body mass and lifespan. We are the first to (1) analyze non-avian reptile groups separately, (2) conduct statistical analyses appropriate for this type of discontinuous data, (3) more accurately use raw cancer occurrence data instead of species averages, and (4) look at how data resampling influences the robustness of results. We found remarkably low cancer prevalence in turtles and high prevalence in squamates and mammals. In contrast to previous studies, lifespan in mammals was negatively associated with neoplasia but positively associated with malignancy. We recovered the same results by reanalyzing data from these previous studies with our statistical approach. We also found that neoplasia prevalence was positively associated with body mass in amphibians and squamates. Overall, our results demonstrate a clear relationship between life history and neoplasia in most tetrapod groups. Our findings also indicate that choice of analytical methods is critical when examining cancer prevalence data. Finally, the exceptionally low cancer prevalence in turtles and extensive variation in cancer prevalence among tetrapod families hold particular promise for identifying species with novel mechanisms of cancer resistance.

Significance Beyond sugars, many types of nectar solutes play important ecological roles; however, the molecular basis for the diversity of nectar composition across species is less explored. One rare trait among flowering plants is the production of colored nectar, which may function to attract and guide prospective pollinators. Our findings indicate convergent evolution of a red-colored nectar has occurred across two distantly related plant species. Behavioral data show that the red pigment attracts diurnal geckos, the likely pollinator of one of these plants. These findings join a growing list of examples of distinct biochemical and molecular mechanisms underlying evolutionary convergence and provide a fascinating system for testing how interactions across species drive the evolution of novel pigments in an understudied context.

Pseudemys alabamensis is one of the most endangered turtle species in the United States due to its small population size and restricted geographic distribution in coastal Alabama and Mississippi. Increased urbanization and climate change impacts in the region further threaten this species. Populations of P. alabamensis are geographically isolated from one another by land and salt water, which could act as barriers to intraspecific gene flow. It is currently unknown how differentiated these isolated populations are from one another or whether they have experienced reductions in population size. Previous work found morphological differences between Alabama and Mississippi populations, suggesting that they may be evolutionarily distinct. Other Pseudemys turtles such as P. concinna and P. floridana occur within the same geographic area as P. alabamensis and are known to hybridize with each other. These more abundant species could further threaten the unique genetic identity of P. alabamensis through introgression. In order to evaluate the endangered status of P. alabamensis and the level of hybridization with other species, we used the mitochondrial (mtDNA) control region and nuclear microsatellite markers to assess genetic variation within and among populations of this species throughout its range and estimate admixture with co-occurring Pseudemys species. Genetic diversity of P. alabamensis was lower than expected at both markers (no variation in mtDNA and excess of homozygosity in microsatellites). We found evidence of genetic differentiation between Alabama and Mississippi populations as well as two populations (Fowl River, Alabama and Biloxi, Mississippi) with low estimated breeding sizes and signs of inbreeding. Finally, we found evidence of admixture of P. alabamensis with P. concinna/P. floridana and Pseudemys peninsularis (a species not native to Alabama or Mississippi). Our results indicate that P. alabamensis is highly endangered throughout its range and threatened by both low population sizes and hybridization. In order to improve the species’ chances of survival, focus should be placed on habitat preservation, maintenance of genetic diversity within both Mississippi and Alabama populations, and regular population monitoring activities such as nest surveillance and estimates of recruitment.

There are many costs associated with increased body size and longevity in animals, including the accumulation of genotoxic and cytotoxic damage that comes with having more cells and living longer. Yet, some species have overcome these barriers and have evolved remarkably large body sizes and long lifespans, sometimes within a narrow window of evolutionary time. Here, we demonstrate through phylogenetic comparative analysis that multiple turtle lineages, including Galapagos giant tortoises, concurrently evolved large bodies, long lifespans, and reduced cancer risk. We also show through comparative genomic analysis that Galapagos giant tortoises have gene duplications related to longevity and tumor suppression. To examine the molecular basis underlying increased body size and lifespan in turtles, we treated cell lines from multiple species, including Galapagos giant tortoises, with drugs that induce different types of cytotoxic stress. Our results indicate that turtle cells, in general, are resistant to oxidative stress related to aging, whereas Galapagos giant tortoise cells, specifically, are sensitive to endoplasmic reticulum stress, which may give this species an ability to mitigate the effects of cellular stress associated with increased body size and longevity.

There are many costs associated with increased body size and longevity in animals, including the accumulation of genotoxic and cytotoxic damage that comes with having more cells and living longer. Yet, some species have overcome these barriers and have evolved remarkably large body sizes and long lifespans, sometimes within a narrow window of evolutionary time. Here, we demonstrate through phylogenetic comparative analysis that multiple turtle lineages, including Galapagos giant tortoises, concurrently evolved large bodies, long lifespans, and reduced cancer risk. We also show through comparative genomic analysis that Galapagos giant tortoises have gene duplications related to longevity and tumor suppression. To examine the molecular basis underlying increased body size and lifespan in turtles, we treated cell lines from multiple species, including Galapagos giant tortoises, with drugs that induce different types of cytotoxic stress. Our results indicate that turtle cells, in general, are resistant to oxidative stress related to aging, while Galapagos giant tortoise cells, specifically, are sensitive to endoplasmic reticulum stress, which may give this species an ability to mitigate the effects of cellular stress associated with increased body size and longevity.

Animals are exposed to different visual stimuli that influence how they perceive and interact with their environment. Visual information such as shape and colour can help the animal detect, discriminate and make appropriate behavioural decisions for mate selection, communication, camouflage, and foraging. In all major vertebrate groups, it has been shown that certain species can discriminate and prefer certain colours and that colours may increase the response to a stimulus. However, because colour is often studied together with other potentially confounding factors, it is still unclear to what extent colour discrimination plays a crucial role in the perception of and attention towards biologically relevant and irrelevant stimuli. To address these questions in reptiles, we assessed the response of three gecko speciesCorrelophus ciliatus, Eublepharis macularius, andPhelsuma laticaudato familiar and novel 2D images in colour or grayscale. We found that while all species responded more often to the novel than to the familiar images, colour information did not influence object discrimination. We also found that the duration of interaction with images was significantly longer for the diurnal species,P. laticauda, than for the two nocturnal species, but this was independent from colouration. Finally, no differences among sexes were observed within or across species. Our results indicate that geckos discriminate between 2D images of different content independent of colouration, suggesting that colouration does not increase detectability or intensity of the response. These results are essential for uncovering which visual stimuli produce a response in animals and furthering our understanding of how animals use colouration and colour vision.

RNA sequencing (RNA-Seq) is becoming a popular method for measuring gene expression in non-model organisms, including wild populations sampled in the field. While RNA-Seq can be used to measure gene expression variation among wild-caught individuals and can yield important biological insights into organismal function, technical variables may also influence gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish species, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq methods: 3’ RNA-Seq and whole mRNA-Seq. We evaluated the effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 minutes after euthanasia on estimated gene expression in blood, gill, muscle, and liver. We found higher RNA degradation in the liver compared to the other tissues. There were fewer expressed genes in blood compared to gill and muscle. We found no difference in gene expression among sampling methods or due to a delay in tissue collection. However, we detected fewer genes with 3’ RNA-Seq than with whole mRNA-Seq and found statistically significant differences in gene expression between 3’ RNA-Seq and whole mRNA-Seq. The magnitude and direction of these differences does not appear to be dependent on gene type or length. Our findings indicate that RNA-Seq is robust to the technical variables related to the field sampling techniques tested here but varies based on the tissue sampled and the RNA-Seq library used. This study advances understanding of usefulness of RNA-Seq to study gene expression variation in evolution, ecology, and conservation.

Cadmium (Cd) accumulates with aging and is elevated in long-lived species. Metallothioneins (MTs), small cysteine-rich proteins involved in metal homeostasis and Cd detoxification, are known to be related to longevity. However, the relationship between Cd accumulation, the role of MTs, and aging is currently unclear. Specifically, we do not know if long-lived species evolved an efficient metal stress response by upregulating their MT levels to reduce the toxic effects of environmental pollutants, such as Cd, that accumulate over their longer life span. It is also unknown if the number of MT genes, their expression, or both protect the organisms from potentially damaging effects during aging. To address these questions, we reanalyzed several cross-species studies and obtained data on MT expression and Cd accumulation in long-lived mouse models. We confirmed a relationship between species maximum life span in captive mammals and their Cd content in liver and kidney. We found that although the number of MT genes does not affect longevity, gene expression and protein amount of specific MT paralogs are strongly related to life span in mammals. MT expression rather than gene number may influence the high Cd levels and longevity of some species. In support of this, we found that overexpression of MT-1 accelerated Cd accumulation in mice and that tissue Cd was higher in long-lived mouse strains with high MT expression. We conclude that long-lived species have evolved a more efficient stress response by upregulating the expression of MT genes in presence of Cd, which contributes to elevated tissue Cd levels. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11357-021-00393-3.

The Galapagos Archipelago is recognized as a natural laboratory for studying evolutionary processes. San Cristóbal was one of the first islands colonized by tortoises, which radiated from there across the archipelago to inhabit 10 islands. Here, we sequenced the mitochondrial control region from six historical giant tortoises from San Cristóbal (five long deceased individuals found in a cave and one found alive during an expedition in 1906) and discovered that the five from the cave are from a clade that is distinct among known Galapagos giant tortoises but closely related to the species from Española and Pinta Islands. The haplotype of the individual collected alive in 1906 is in the same clade as the haplotype in the contemporary population. To search for traces of a second lineage in the contemporary population on San Cristóbal, we closely examined the population by sequencing the mitochondrial control region for 129 individuals and genotyping 70 of these for both 21 microsatellite loci and12,000 genome-wide single nucleotide polymorphisms [SNPs]. Only a single mitochondrial haplotype was found, with no evidence to suggest substructure based on the nuclear markers. Given the geographic and temporal proximity of the two deeply divergent mitochondrial lineages in the historical samples, they were likely sympatric, raising the possibility that the lineages coexisted. Without the museum samples, this important discovery of an additional lineage of Galapagos giant tortoise would not have been possible, underscoring the value of such collections and providing insights into the early evolution of this iconic radiation.

The dorsal surfaces of many taxonomic groups often feature repetitive pattern elements consisting of stripes, spots, or bands. Here, we investigate how distinct categories of camouflage pattern work by relating them to ecological and behavioral traits in 439 species of gecko. We use phylogenetic comparative methods to test outstanding hypotheses based on camouflage theory and research in other taxa. We found that bands are associated with nocturnal activity, suggesting bands provide effective camouflage for motionless geckos resting in refugia during the day. A predicted association between stripes and diurnal activity was not supported, suggesting that stripes do not work via dazzle camouflage mechanisms in geckos. This, along with a lack of support for our prediction that plain patterning should be associated with open habitats, suggests that similar camouflage patterns do not work in consistent ways across taxa. We also found that plain and striped lineages frequently switched between using open or closed habitats, whereas spotted lineages rarely transitioned. This suggests that pattern categories differ in how specialized or generalized their camouflage is. This result has ramifications for theory on how camouflage compromises to background heterogeneity and how camouflage pattern might influence evolutionary trajectories.

Chemical use in society is growing rapidly and is one of the five major pressures on biodiversity worldwide. Since empirical toxicity studies of pollutants generally focus on a handful of model organisms, reliable approaches are needed to assess sensitivity to chemicals across the wide variety of species in the environment. Phylogenetic comparative methods (PCM) offer a promising approach for toxicity extrapolation incorporating known evolutionary relationships among species. If phylogenetic signal in toxicity data is high – i.e., closely related species are more similarly sensitive as compared to distantly related species – PCM could ultimately help predict species sensitivity when toxicity data are lacking. Here, we present the largest ever test of phylogenetic signal in toxicity data by combining phylogenetic data from fish with acute mortality data for 42 chemicals spanning 10 different chemical classes. Phylogenetic signal is high for some chemicals, particularly organophosphate pesticides, but not necessarily for many chemicals in other classes (e.g., metals, organochlorines). These results demonstrate that PCM may be useful for toxicity extrapolation in untested species for those chemicals with clear phylogenetic signal. This study provides a framework for using PCM to understand the patterns and causes of variation in species sensitivity to pollutants.

DNA damage has been linked to genomic instability and the progressive breakdown of cellular and organismal homeostasis, leading to the onset of disease and reduced longevity. Insults to DNA from endogenous sources include base deamination, base hydrolysis, base alkylation, and metabolism-induced oxidative damage that can lead to single-strand and double-strand DNA breaks. Alternatively, exposure to environmental pollutants, radiation or ultra-violet light, can also contribute to exogenously derived DNA damage. We previously validated a novel, high through-put approach to measure levels of DNA damage in cultured mammalian cells. This new CometChip Platform builds on the classical single cell gel electrophoresis or comet methodology used extensively in environmental toxicology and molecular biology. We asked whether the CometChip Platform could be used to measure DNA damage in samples derived from environmental field studies. To this end, we determined that nucleated erythrocytes from multiple species of turtle could be successfully evaluated in the CometChip Platform to quantify levels of DNA damage. In total, we compared levels of DNA damage in 40 animals from two species: the box turtle (Terrapene carolina) and the red-eared slider (Trachemys scripta elegans). Endogenous levels of DNA damage were identical between the two species, yet we did discover some sex-linked differences and changes in DNA damage accumulation. Based on these results, we confirm that the CometChip Platform allows for the measurement of DNA damage in a large number of samples quickly and accurately, and is particularly adaptable to environmental studies using field-collected samples. Environ. Mol. Mutagen. 59:322-333, 2018. © 2018 Wiley Periodicals, Inc.

Giant tortoises are among the longest-lived vertebrate animals and, as such, provide an excellent model to study traits like longevity and age-related diseases. However, genomic and molecular evolutionary information on giant tortoises is scarce. Here, we describe a global analysis of the genomes of Lonesome George—the iconic last member of Chelonoidis abingdonii—and the Aldabra giant tortoise (Aldabrachelys gigantea). Comparison of these genomes with those of related species, using both unsupervised and supervised analyses, led us to detect lineage-specific variants affecting DNA repair genes, inflammatory mediators and genes related to cancer development. Our study also hints at specific evolutionary strategies linked to increased lifespan, and expands our understanding of the genomic determinants of ageing. These new genome sequences also provide important resources to help the efforts for restoration of giant tortoise populations., We thank Banco Santander for funding a short stay of S.F.-R. and D.C.-I. at Yale University. V.Q. is supported by grants from the Principado de Asturias and Ministerio de Economía y Competitividad, including FEDER funding. L.F.K.K. is supported by an FPI fellowship associated with BFU2014-55090-P (FEDER). T.M.-B. is supported by MINECO BFU2017-86471-P (MINECO/FEDER, UE), an NIH U01 MH106874 grant, the Howard Hughes International Early Career programme, Obra Social ‘La Caixa’ and Secretaria d’Universitats i Recerca, and CERCA Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya. C.L.-O. is supported by grants from the European Research Council (DeAge; ERC Advanced Grant), Ministerio de Economía y Competitividad, Instituto de Salud Carlos III (RTICC) and Progeria Research Foundation. The Instituto Universitario de Oncología is supported by Fundación Bancaria Caja de Ahorros de Asturias. We also thank staff at the Galapagos National Park and Galapagos Conservancy for logistic and financial support.

Seagrass meadows form ecologically and economically valuable coastal habitat on every continental margin except the Antarctic, but their areal extent is declining by approximately 2–5 % per year. Seagrass wasting disease is a contributing factor in these declines, with the protist Labyrinthula identified as the etiologic agent. To help elucidate the role of Labyrinthula spp. in global seagrass declines, we surveyed roughly one fourth of all seagrass species to identify Labyrinthula diversity at the strain and/or species level, combining results from culturing methods and two common nuclear DNA markers: the ITS and 18S regions of the ribosomal RNA gene complex. After assaying a subset of the resulting isolates (of which 170 were newly sequenced), we produced a cladogenic context for putative seagrass-pathogenic versus non-pathogenic Labyrinthula while also defining host and geographic ranges. Assays also suggest that pathogenicity is consistently high (when present; and, even when comparing susceptibility of US East- versus West Coast Zostera marina hosts) while virulence is variable, that some isolate-host combinations have the potential for host cross-infection, and that several modes of transmission can be effective. Taken together, these data provide additional means for delimiting putative species of Labyrinthula, suggesting at least five seagrass-pathogenic and perhaps ten or more non-pathogenic marine “species”, yielding a working definition for ecologists and epidemiologists attempting to reconcile the sundry data related to seagrass wasting disease.

Oceanic islands are often inhabited by endemic species that have undergone substantial morphological evolutionary change due to processes of multiple colonizations from various source populations, dispersal, and local adaptation. Galápagos marine iguanas are an example of an island endemic exhibiting high morphological diversity, including substantial body size variation among populations and sexes, but the causes and magnitude of this variation are not well understood. We obtained morphological measurements from marine iguanas throughout their distribution range. These data were combined with genetic and local environmental data from each population to investigate the effects of evolutionary history and environmental conditions on body size and shape variation and sexual dimorphism. Our results indicate that body size and shape are highly variable among populations. Sea surface temperature and island perimeter, but not evolutionary history as depicted by phylogeographic patterns in this species, explain variation in body size among populations. Conversely, evolutionary history, but not environmental parameters or island size, was found to influence variation in body shape among populations. Finally, in all populations except one, we found strong sexual dimorphism in body size and shape in which males are larger, with higher heads than females, while females have longer heads than males. Differences among populations suggest that plasticity and/or genetic adaptation may shape body size and shape variation in marine iguanas. This study will help target future investigations to address the contribution of plasticity versus genetic adaptation on size and shape variation in marine iguanas.

Genotypic variation, environmental variation, and their interaction may produce variation in the developmental process and cause phenotypic differences among individuals. Developmental noise, which arises during development from stochasticity in cellular and molecular processes when genotype and environment are fixed, also contributes to phenotypic variation. While evolutionary biology has long focused on teasing apart the relative contribution of genes and environment to phenotypic variation, our understanding of the role of developmental noise has lagged due to technical difficulties in directly measuring the contribution of developmental noise. The influence of developmental noise is likely underestimated in studies of phenotypic variation due to intrinsic mechanisms within organisms that stabilize phenotypes and decrease variation. Since we are just beginning to appreciate the extent to which phenotypic variation due to stochasticity is potentially adaptive, the contribution of developmental noise to phenotypic variation must be separated and measured to fully understand its role in evolution. Here, we show that variation in the component of the developmental process corresponding to environmental and genetic factors (here treated together as a unit called the LALI-type) versus the contribution of developmental noise, can be distinguished for leopard gecko (Eublepharis macularius) head color patterns using mathematical simulations that model the role of random variation (corresponding to developmental noise) in patterning. Specifically, we modified the parameters of simulations corresponding to variation in the LALI-type to generate the full range of phenotypic variation in color pattern seen on the heads of eight leopard geckos. We observed that over the range of these parameters, variation in color pattern due to LALI-type variation exceeds that due to developmental noise in the studied gecko cohort. However, the effect of developmental noise on patterning is also substantial. Our approach addresses one of the major goals of evolutionary biology: to quantify the role of stochasticity in shaping phenotypic variation., Author summary The observable characteristics of an organism make up its phenotype. Variation among phenotypes is due to genetic differences, environmental factors and developmental noise (effects due to inherent stochasticity) during development. We used mathematical models to investigate the contributions of variation of the developmental process due to genetic and environmental factors (treated in this work as a single unit) versus developmental noise (unavoidable variation within the developmental program) to the development of pigment patterns on gecko heads. We found that for our cohort, the proportion of phenotypic variation due to variation in the unit composed of genotypic and environmental variation is larger than that due to developmental noise. Furthermore, by allowing the parameters of the mathematical model to vary, we generated the full extent of potential phenotypic pattern variation that could occur on the head of geckos. This serves to further study the influence of the buffering mechanisms (canalization, selection, and developmental stability) limiting phenotypic variation. This approach can be applied to any regular morphological trait that results from self-organized processes such as reaction-diffusion mechanisms, including the frequently found striped and spotted patterns of animal pigmentation patterning, patterning of bones in vertebrate limbs, and body segmentation in segmented animals.

Self-righting, the capacity of an animal to self-turn after falling on its back, is a fitness-related trait. Delayed self-righting can result in loss of mating opportunities or death. Traits involved in self-righting may therefore be under selection. Galápagos giant tortoises have two main shell morphologies - saddleback and domed – that have been proposed to be adaptive. The more sloped shape on the sides of the shell and the longer extension of neck and legs of the saddlebacks could have evolved to optimize self-righting. The drier environments with more uneven surfaces where the saddleback tortoises occur increases their risk to fall on their back while walking. The ability to fast overturn could reduce the danger of dying. To test this hypothesis, we used 3D shell reconstructions of 89 Galápagos giant tortoises from three domed and two saddleback species to compare self-righting potential of the two shell morphotypes. Our results indicate that saddleback shells require higher energy input to self-right than domed ones. This suggests that several traits associated with the saddleback shell morphology could have evolved to facilitate self-righting. Studying the functional performances of fitness-related traits, as in this work, could provide important insight into the adaptive value of traits.

DNA damage has been linked to genomic instability and the progressive breakdown of cellular and organismal homeostasis, leading to the onset of disease and reduced longevity. Insults to DNA from endogenous sources include base deamination, base hydrolysis, base alkylation, and metabolism-induced oxidative damage that can lead to single-strand and double-strand DNA breaks. Alternatively, exposure to environmental pollutants, radiation or ultra-violet light, can also contribute to exogenously derived DNA damage. We previously validated a novel, high through-put approach to measure levels of DNA damage in cultured mammalian cells. This new CometChip Platform builds on the classical single cell gel electrophoresis or comet methodology used extensively in environmental toxicology and molecular biology. We asked whether the CometChip Platform could be used to measure DNA damage in samples derived from environmental field studies. To this end, we determined that nucleated erythrocytes from multiple species of turtle could be successfully evaluated in the CometChip Platform to quantify levels of DNA damage. In total, we compared levels of DNA damage in 40 animals from two species: the box turtle (Terrapene carolina) and the red-eared slider (Trachemys scripta elegans). Endogenous levels of DNA damage were identical between the two species, yet we did discover some sex-linked differences and changes in DNA damage accumulation. Based on these results, we confirm that the CometChip Platform allows for the measurement of DNA damage in a large number of samples quickly and accurately, and is particularly adaptable to environmental studies using field-collected samples. Environ. Mol. Mutagen. 59:322-333, 2018. © 2018 Wiley Periodicals, Inc.

The metallothionein (MT) gene superfamily consists of metal-binding proteins involved in various metal detoxification and storage mechanisms. The evolution of this gene family in vertebrates has mostly been studied in mammals using sparse taxon or gene sampling. Genomic databases and available data on MT protein function and expression allow a better understanding of the evolution and functional divergence of the different MT types. We recovered 77 MT coding sequences from 20 representative vertebrates with annotated complete genomes. We found multiple MT genes, also in reptiles, which were thought to have only one MT type. Phylogenetic and synteny analyses indicate the existence of a eutherian MT1 and MT2, a tetrapod MT3, an amniote MT4, and fish MT. The optimal gene-tree/species-tree reconciliation analyses identified the best root in the fish clade. Functional analyses reveal variation in hydropathic index among protein domains, likely correlated with their distinct flexibility and metal affinity. Analyses of functional divergence identified amino acid sites correlated with functional divergence among MT types. Uncovering the number of genes and sites possibly correlated with functional divergence will help to design cost-effective MT functional and gene expression studies. This will permit further understanding of the distinct roles and specificity of these proteins and to properly target specific MT for different types of functional studies. Therefore, this work presents a critical background on the molecular evolution and functional divergence of vertebrate MTs to carry out further detailed studies on the relationship between heavy metal metabolism and tolerances among vertebrates.

Amphibian declines have been documented worldwide. Chytridiomycosis, a disease caused by the aquatic fungus Batrachochytrium dendrobatidis (Bd), is one of the causes associated with these declines. Occurrence, rate of infection and mortality due to chytridiomycosis in amphibians depend on multiple factors including habitat and life-style (aquatic/terrestrial). Bd infection is lower in terrestrial than in aquatic species, but a fully terrestrial life-style alone may not explain the absence of Bd in some species. Low individual dispersal, decreasing the occurrence of contact with infected organisms, could also favour lower Bd infection. To survey the occurrence of Bd infection in fully terrestrial salamanders with low dispersal, we sampled the Sardinian Hydromantes species to measure their level of infection. Bd was not detected and likely absent in Sardinian Hydromantes. This phenomenon could be explained by a combination of terrestrial habitat, low dispersal, and occurrence mostly in habitats where other amphibians do not occur.

Neutral rates of molecular evolution vary across species, and this variation has been shown to be related to biological traits. One of the first patterns to be observed in vertebrates has been an inverse relationship between body mass (BM) and substitution rates. The effects of three major life-history traits (LHT) that covary with BM - metabolic rate, generation time and longevity (LON) - have been invoked to explain this relationship. However, most of the theoretical and empirical evidence supporting this relationship comes from endothermic vertebrates, that is, mammals and birds, in which the environmental conditions, especially temperature, do not have a direct impact on cellular and molecular biology. We analysed the variations in mitochondrial and nuclear rates of synonymous substitution across 224 turtle species and examined their correlation with two LHT (LON and BM) and two environmental variables [latitude (LAT) and habitat]. Our analyses indicate that in turtles, neutral rates of molecular evolution are hardly correlated with LON or BM. Rather, both the mitochondrial and nuclear substitution rates are significantly correlated with LAT - faster evolution in the tropics - and especially so for aquatic species. These results question the generality of the relationships reported in mammals and birds and suggest that environmental factors might be the strongest determinants of the mutation rate in ectotherms.

The superfamily Testudinoidea is the most diverse and widely distributed clade of extant turtles. Surprisingly, despite an extensive fossil record, and increasing amount of molecular data available, the temporal origin of this group is still largely unknown. To address this issue, we used a comprehensive molecular dataset to perform phylogenetic and molecular dating analyses, as well as seven fossil constraints to calibrate the ages of the nodes in the phylogeny. The molecular dataset includes the complete mitochondrial genomes of 37 turtle species, including newly sequenced mitochondrial genomes of Phrynops hilarii, Emys orbicularis, Rhinoclemmys punctularia, and Chelonoidis nigra, and four nuclear markers. Our results revealed that the earliest divergences within crown testudinoids occurred around 95.0 Mya, in the early Late Cretaceous, earlier than previously reported, raising new questions about the historical biogeography of this group.

Next Generation Sequencing technologies (NGS) are rapidly invading many evolutionary and ecological fields, such as phylogenomics, molecular evolution, population genomics and molecular ecology. Among the potential targets of NGS is transcriptome sequencing, a fast and relatively cheap way to generate massive amounts of coding sequence data, offering promising perspectives for the analysis of molecular diversity in the wild. A number of molecular ecology research groups therefore may switch from DNA-based to RNA-based typing in the near future. Sample preparation from natural populations, however, requires specific care and protocols when RNA is the target. Furthermore, NGS sequencing of transcriptome requires high amount of good-quality RNA. Here we present the results of RNA extraction experiments from various samples of 39 animal species caught in the wild. We compared tissue preparation and storage conditions, evaluated and improved standard RNA extraction protocols, and achieved RNA yield and quality suitable for NGS in all cases. We derive general guidelines for the production of ready-to-sequence RNA in nonmodel animals sampled in the field.

Galapagos tortoises possess two main shell forms, domed and saddleback, that correlate with the biogeographic history of this species group. However, the lack of description of morphological shell variation within and among populations has prevented the under- standing of the contribution of evolutionary forces and the potential role of ontogeny in shaping morphological shell differences. Here, we analyze two lineages of Galapagos tortoises inhabiting Santa Cruz Island by applying geometric morphometrics in combination with a photogrammetry 3D reconstruction method on a set of tortoises of different ages (from juvenile to adult). The aim of this study is to describe morphological features on the carapace that could be used for taxonomic recog- nition by taking into account confounding factors, such as the morphological changes occurring during growth. Our results indicate that despite the shared similarities of growth patterns and of morphological changes observed during growth, the two lineages and the differ- ent sexes can be distinguished on the basis of distinct carapace features. Lineages differ by the shape of the vertebral (especially concerning their width) and pleural scutes, with one lineage having a more compressed carapace shape, whereas the other possesses a carapace that is more elongated and expanded toward the sides as well as an higher positioning of the first vertebral scute. Furthermore, females have a more elongated and wider carapace shape than males. Finally, carapace shape changes with growth, with vertebral scutes becoming narrower and pleural scutes becoming larger during late ontogeny. J. Morphol. 272:379-386, 2011. 2011 Wiley-Liss, Inc.

Cancer researchers have traditionally used the mouse and the rat as staple model organisms. These animals are very short-lived, reproduce rapidly, and are highly prone to cancer. They have been very useful for modeling some human cancer types and testing experimental treatments; however, these cancer-prone species offer little for understanding the mechanisms of cancer resistance. Recent technological advances have expanded research bestiary to non-standard model organisms that possess unique traits of very high value to humans such as cancer resistance and longevity. In recent years, several discoveries have been made in non-standard mammalian species providing new insights on the natural mechanisms of cancer resistance. These include mechanisms of cancer resistance in the naked mole rat, blind mole rat and elephant. In each of these species, evolution took a different path leading to novel mechanisms. Many other long-lived mammalian species display cancer resistance, including whales, grey squirrels, microbats, cows and horses. Understanding the molecular mechanisms of cancer resistance in all of these species is important and timely as ultimately, these mechanisms could be harnessed for the development of human cancer therapies.

We examined two image-based methods, photogrammetry and stereo vision, used for reconstructing the threedimensional form of biological organisms under field conditions. We also developed and tested a third ‘hybrid’ method, which combines the other two techniques. We tested these three methodologies using two different cameras to obtain digital images of museum and field sampled specimens of giant tortoises. Both the precision and repeatability of the methods were assessed statistically on the same specimens by comparing geodesic and Euclidean measurements made on the digital models with linear measurements obtained with caliper and flexible tape. We found no substantial difference between the three methods in measuring the Euclidean distance between landmarks, but spatially denser models (stereo vision and ‘hybrid’) were more accurate for geodesic distances. The use of different digital cameras did not influence the results. Image-based methods require only inexpensive instruments and appropriate software, and allow reconstruction of the three-dimensional forms (including their curved surfaces) of organisms sampled in the field. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95, 425–436.

Malagasy poison frogs of the genus Mantella with its 16 species are currently sub-divided into 5 major groups. Of these, the Mantella betsileo group is traditionally understood as containing four species, Mantella betsileo, M. expectata, M. viridis and M. manery, while the M. laevigata group is considered to be monospecific. A phylogenetic analysis of samples from multiple localities of all species in these two groups, based on sequences of the mitochondrial cytochrome b gene, shows the existence of several well-distinct clades in what is currently considered to be Mantella betsileo: (1) central-western populations from Kirindy, Isalo, and near Antsirabe close to the Betsileo region, to which the name M. betsileo is to be applied, (2) populations of the north-east and north-west, which are closely related to M. viridis and to which the name M. ebenaui is to be applied, and (3) a clade from southernmost Madagascar and from the Tsingy de Bemaraha, which is sister to M. expectata and furthermore includes important intra-clade variation, therefore probably representing one or two undescribed species. Our data also support a large genetic distance of M. manery to all other species and its probable sister-group relationship to the sympatric M. laevigata; M. manery is consequently transferred from the M. betsileo group to the M. laevigata group.

The Magdalena River Turtle (Podocnemis lewyana) is a Colombian endemic species, endangered due to human exploitation and habitat destruction. To date, this species is poorly known ecologically and data on its genetic diversity are lacking. Here we report on the first genetic survey of the species across its distribution range. We obtained mitochondrial DNA sequences (488 bp) of the cytochrome b gene from 109 individuals. Samples belong to populations located at several different localities, grouped in five regions, along the four main river basins: Magdalena, Cauca, San Jorge, and Sinú drainages. We found two haplotypes, which differ in only one nucleotide substitution and which are represented with different frequencies in the five geographic regions. These results suggest that P. lewyana harbors little genetic variation and is a genetically uniform species, but more variable markers (i.e., microsatellites) should be used to unravel fine-scale phylogeographic structures in this species.

The release of large quantities of chemicals into the environment represents a major source of environmental disturbance. In recent years, the focus of ecotoxicology has shifted from describing the effects of chemical contaminants on individual species to developing more integrated approaches for predicting and evaluating long term effects of chemicals across species and ecosystems. Traditional ecotoxicology is typically based on data of sensitivity to a contaminant of a few surrogate species and often considers little variability in chemical sensitivity within and among taxonomic groups. This approach assumes that evolutionary history and phylogenetic relatedness among species have little or no impact on species' sensitivity to chemical compounds. Few studies have tested this assumption. Using phylogenetic comparative methods and published data for amphibians, we show that sensitivity to copper sulfate, a commonly used pesticide, exhibits a strong phylogenetic signal when controlling for experimental temperature. Our results indicate that evolutionary history needs to be accounted for to make accurate predictions of amphibian sensitivity to this contaminant under different temperature scenarios. Since physiological and metabolic traits showing high phylogenetic signal likely underlie variation in species sensitivity to chemical stressors, future studies should evaluate and predict species vulnerability to pollutants using evolutionarily informed approaches.

Background Amelotin (AMTN) is an ameloblast-secreted protein that belongs to the secretory calcium-binding phosphoprotein (SCPP) family, which originated in early vertebrates. In rodents, AMTN is expressed during the maturation stage of amelogenesis only. This expression pattern strongly differs from the spatiotemporal expression of other ameloblast-secreted SCPPs, such as the enamel matrix proteins (EMPs). Furthermore, AMTN was characterized in rodents only. In this study, we applied various approaches, including in silico screening of databases, PCRs and transcriptome sequencing to characterize AMTN sequences in sauropsids and amphibians, and compared them to available mammalian and coelacanth sequences. Results We showed that (i) AMTN is tooth (enamel) specific and underwent pseudogenization in toothless turtles and birds, and (ii) the AMTN structure changed during tetrapod evolution. To infer AMTN function, we studied spatiotemporal expression of AMTN during amelogenesis in a salamander and a lizard, and compared the results with available expression data from mouse. We found that AMTN is expressed throughout amelogenesis in non-mammalian tetrapods, in contrast to its expression limited to enamel maturation in rodents. Conclusions Taken together our findings suggest that AMTN was primarily an EMP. Its functions were conserved in amphibians and sauropsids while a change occurred early in the mammalian lineage, modifying its expression pattern during amelogenesis and its gene structure. These changes likely led to a partial loss of AMTN function and could have a link with the emergence of prismatic enamel in mammals. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0329-x) contains supplementary material, which is available to authorized users.

Aim Amphibians are a model group for studies of the biogeographical origins of salt-intolerant taxa on oceanic islands. We used the Gulf of Guinea islands to explore the biogeographical origins of island endemism of one species of frog, and used this to gain insights into potential colonization mechanisms. Location Sao Tome and Principe, two of the four major islands in the Gulf of Guinea, West Africa, are truly oceanic and have an exceptionally high biodiversity. Methods Mitochondrial DNA is used to test the endemic status of a frog from Sao Tome and compare it with congeneric taxa from tropical Africa. Existing data on surface currents, surface salinity, atmospheric circulation and bird migration in the Gulf of Guinea are summarized to address hypotheses concerning colonization mechanisms. Results The endemic status of Ptychadena newtoni (Bocage) is supported here by mitochondrial DNA sequences, and analysis of this and other molecular data indicates that an East African species close to Ptychadena mascareniensis (Dumeril and Bibron) is its nearest relative. We refute the possibility that this population was anthropogenically introduced, in favour of a natural dispersal mechanism. Main conclusions With six endemic frogs and one caecilian, the Gulf of Guinea islands harbour a diverse amphibian fauna. Five of these species appear to have their closest relatives in East Africa. Insufficient evidence exists for transportation by storms, birds or rafts alone. However, we propose a synergy of rafting, favourable surface currents and a reduction in salinity of surface waters. Catastrophic events, or wet periods in climatic history, could allow freshwater paths to open far enough to enable continental flora and fauna to reach these and other isolated oceanic islands.