Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation.

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect 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 (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

Hybridization and low genetic diversity in the endangered Alabama red-bellied turtle (Pseudemys alabamensis).

Pseudemys alabamensis is one of the most endangered freshwater turtle species in the United States due to its restricted geographic distribution in coastal Alabama and Mississippi. Populations of P. alabamensis are geographically isolated from one another by land and saltwater, which could act as barriers to gene flow. It is currently unknown how differentiated these populations are from one another and 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 naturally within the same geographic area as P. alabamensis and are known to hybridize with each other. These more abundant species could 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 mitochondrial and nuclear microsatellite markers to assess genetic variation within and among populations of P. alabamensis throughout its range and estimate admixture with co-occurring Pseudemys species. In P. alabamensis, we found no variation in mitochondrial DNA and an excess of homozygosity in microsatellite data. Our results show genetic differentiation between Alabama and Mississippi populations of P. alabamensis, and low estimated breeding sizes and signs of inbreeding for two populations (Fowl River, Alabama and Biloxi, Mississippi). We also found evidence of admixture between P. alabamensis and P. concinna/P. floridana. Based on our results, 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 the Mississippi and Alabama populations, and routine population-monitoring activities such as nest surveillance and estimates of recruitment.

Convergent evolution of a blood-red nectar pigment in vertebrate-pollinated flowers.

Nearly 90% of flowering plants depend on animals for reproduction. One of the main rewards plants offer to pollinators for visitation is nectar. Nesocodon mauritianus (Campanulaceae) produces a blood-red nectar that has been proposed to serve as a visual attractant for pollinator visitation. Here, we show that the nectar's red color is derived from a previously undescribed alkaloid termed nesocodin. The first nectar produced is acidic and pale yellow in color, but slowly becomes alkaline before taking on its characteristic red color. Three enzymes secreted into the nectar are either necessary or sufficient for pigment production, including a carbonic anhydrase that increases nectar pH, an aryl-alcohol oxidase that produces a pigment precursor, and a ferritin-like catalase that protects the pigment from degradation by hydrogen peroxide. Our findings demonstrate how these three enzymatic activities allow for the condensation of sinapaldehyde and proline to form a pigment with a stable imine bond. We subsequently verified that synthetic nesocodin is indeed attractive to Phelsuma geckos, the most likely pollinators of Nesocodon We also identify nesocodin in the red nectar of the distantly related and hummingbird-visited Jaltomata herrerae and provide molecular evidence for convergent evolution of this trait. This work cumulatively identifies a convergently evolved trait in two vertebrate-pollinated species, suggesting that the red pigment is selectively favored and that only a limited number of compounds are likely to underlie this type of adaptation.

Capturing and analyzing pattern diversity: an example using the melanistic spotted patterns of leopard geckos.

Animal color patterns are widely studied in ecology, evolution, and through mathematical modeling. Patterns may vary among distinct body parts such as the head, trunk or tail. As large amounts of photographic data is becoming more easily available, there is a growing need for general quantitative methods for capturing and analyzing the full complexity and details of pattern variation. Detailed information on variation in color pattern elements is necessary to understand how patterns are produced and established during development, and which evolutionary forces may constrain such a variation. Here, we develop an approach to capture and analyze variation in melanistic color pattern elements in leopard geckos. We use this data to study the variation among different body parts of leopard geckos and to draw inferences about their development. We compare patterns using 14 different indices such as the ratio of melanistic versus total area, the ellipticity of spots, and the size of spots and use these to define a composite distance between two patterns. Pattern presence/absence among the different body parts indicates a clear pathway of pattern establishment from the head to the back legs. Together with weak within-individual correlation between leg patterns and main body patterns, this suggests that pattern establishment in the head and tail may be independent from the rest of the body. We found that patterns vary greatest in size and density of the spots among body parts and individuals, but little in their average shapes. We also found a correlation between the melanistic patterns of the two front legs, as well as the two back legs, and also between the head, tail and trunk, especially for the density and size of the spots, but not their shape or inter-spot distance. Our data collection and analysis approach can be applied to other organisms to study variation in color patterns between body parts and to address questions on pattern formation and establishment in animals.

Ecological, behavioral, and phylogenetic influences on the evolution of dorsal color pattern in geckos.

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.

Isolating and quantifying the role of developmental noise in generating phenotypic variation.

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.

Application of the CometChip platform to assess DNA damage in field-collected blood samples from turtles.

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.