Signaling dynamics in coexisting monoclonal cell subpopulations unveil mechanisms of resistance to anti-cancer compounds.

BACKGROUND: Tumor heterogeneity is a main contributor of resistance to anti-cancer targeted agents though it has proven difficult to study. Unfortunately, model systems to functionally characterize and mechanistically study dynamic responses to treatment across coexisting subpopulations of cancer cells remain a missing need in oncology. METHODS: Using single cell cloning and expansion techniques, we established monoclonal cell subpopulations (MCPs) from a commercially available epidermal growth factor receptor (EGFR)-mutant non-small cell lung cancer cell line. We then used this model sensitivity to the EGFR inhibitor osimertinib across coexisting cell populations within the same tumor. Pathway-centered signaling dynamics associated with response to treatment and morphological characteristics of the MCPs were assessed using Reverse Phase Protein Microarray. Signaling nodes differentially activated in MCPs less sensitive to treatment were then pharmacologically inhibited to identify target signaling proteins putatively implicated in promoting drug resistance. RESULTS: MCPs demonstrated highly heterogeneous sensitivities to osimertinib. Cell viability after treatment increased > 20% compared to the parental line in selected MCPs, whereas viability decreased by 75% in other MCPs. Reduced treatment response was detected in MCPs with higher proliferation rates, EGFR L858R expression, activation of EGFR binding partners and downstream signaling molecules, and expression of epithelial-to-mesenchymal transition markers. Levels of activation of EGFR binding partners and MCPs' proliferation rates were also associated with response to c-MET and IGFR inhibitors. CONCLUSIONS: MCPs represent a suitable model system to characterize heterogeneous biomolecular behaviors in preclinical studies and identify and functionally test biological mechanisms associated with resistance to targeted therapeutics.

Thermal tolerance plasticity and dynamics of thermal tolerance in Eublepharis macularius: Implications for future climate-driven heat stress.

The intensity and duration of heat waves, as well as average global temperatures, are expected to increase due to climate change. Heat waves can cause physiological stress and reduce fitness in animals. Species can reduce overheating risk through phenotypic plasticity, which allows them to raise their thermal tolerance limits over time. This mechanism could be important for ectotherms whose body temperatures are directly influenced by available environmental temperatures. Geckos are a large, diverse group of ectotherms that vary in their thermal habitats and times of daily activity, which could affect how they physiologically adjust to heat waves. Data on thermal physiology are scarce for reptiles, with only one study in geckos. Understanding thermal tolerance and plasticity, and their relationship, is essential for understanding how some species are able to adjust or adapt to changing temperatures. In this study, we estimated thermal tolerance and plasticity, and their interaction, in the crepuscular gecko, Eublepharis macularius, a species that is emerging as a model for reptile biology. After estimating basal thermal tolerance for 28 geckos, thermal tolerance was measured for each individual a second time at several timepoints (3, 6, or 24 h) to determine thermal tolerance plasticity. We found that thermal tolerance plasticity (1) does not depend on the basal thermal tolerance of the organism, (2) was highest after 6 h from initial heat shock, and (3) was negatively influenced by individual body mass. Our findings contribute to the increasing body of work focused on understanding the influence of biological and environmental factors on thermal tolerance plasticity in organisms and provide phenotypic data to further investigate the molecular basis of thermal tolerance plasticity in organisms.

Melanistic coloration does not influence thermoregulation in the crepuscular gecko Eublepharis macularius.

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 melanistic coloration for thermoregulation functions to increase solar radiation absorption and consequently heating in many diurnal ectotherms, research on crepuscular and nocturnal ectotherms is lacking. Since crepuscular and nocturnal ectotherms generally absorb heat from the substrate, in these organisms melanistic coloration may have other primary functions beside thermoregulation. As such, in this work we hypothesized that the proportion of dorsal melanistic body coloration would not influence heating and cooling rates in the crepuscular gecko, Eublepharis macularius, and that changes in environmental temperature would not 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 support that melanistic coloration in E. macularius may not be used for thermoregulation and strengthen the hypothesis that in animals active in low light conditions, melanistic coloration may be used instead for camouflage or other functions.

The revised reference genome of the leopard gecko (Eublepharis macularius) provides insight into the considerations of genome phasing and assembly.

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. macularius reference 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 previous E. macularius reference 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 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the 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 of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000.

The revised reference genome of the leopard gecko ( Eublepharis macularius ) provides insight into the considerations of genome phasing and assembly.

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. macularius reference 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 previous E. macularius reference 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 new E. macularius reference assembly is available on NCBI at JAOPLA010000000. The genome version and its associated annotations are also available via this Figshare repository https://doi.org/10.6084/m9.figshare.20069273 .

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.

A new lineage of Galapagos giant tortoises identified from museum samples.

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.

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.

Concurrent Evolution of Antiaging Gene Duplications and Cellular Phenotypes in Long-Lived Turtles.

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.

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.

Elevated metallothionein expression in long-lived species mediates the influence of cadmium accumulation on aging.

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.

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.

Giant tortoise genomes provide insights into longevity and age-related disease.

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.

Mixed phylogenetic signal in fish toxicity data across chemical classes.

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.

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.

Self-righting potential and the evolution of shell shape in Galápagos tortoises.

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.