Repeated evolution of circadian clock dysregulation in cavefish populations.

Circadian rhythms are nearly ubiquitous throughout nature, suggesting they are critical for survival in diverse environments. Organisms inhabiting largely arrhythmic environments, such as caves, offer a unique opportunity to study the evolution of circadian rhythms in response to changing ecological pressures. Populations of the Mexican tetra, Astyanax mexicanus, have repeatedly invaded caves from surface rivers, where individuals must contend with perpetual darkness, reduced food availability, and limited fluctuations in daily environmental cues. To investigate the molecular basis for evolved changes in circadian rhythms, we investigated rhythmic transcription across multiple independently-evolved cavefish populations. Our findings reveal that evolution in a cave environment has led to the repeated disruption of the endogenous biological clock, and its entrainment by light. The circadian transcriptome shows widespread reductions and losses of rhythmic transcription and changes to the timing of the activation/repression of core-transcriptional clock. In addition to dysregulation of the core clock, we find that rhythmic transcription of the melatonin regulator aanat2 and melatonin rhythms are disrupted in cavefish under darkness. Mutants of aanat2 and core clock gene rorca disrupt diurnal regulation of sleep in A. mexicanus, phenocopying circadian modulation of sleep and activity phenotypes of cave populations. Together, these findings reveal multiple independent mechanisms for loss of circadian rhythms in cavefish populations and provide a platform for studying how evolved changes in the biological clock can contribute to variation in sleep and circadian behavior.

Performance of gene expression analyses using de novo assembled transcripts in polyploid species.

MOTIVATION: Quality of gene expression analyses using de novo assembled transcripts in species that experienced recent polyploidization remains unexplored. RESULTS: Differential gene expression (DGE) analyses using putative genes inferred by Trinity, Corset and Grouper performed slightly differently across five plant species that experienced various polyploidy histories. In species that lack recent polyploidy events that occurred in the past several millions of years, DGE analyses using de novo assembled transcriptomes identified 54-82% of the differentially expressed genes recovered by mapping reads to the reference genes. However, in species that experienced more recent polyploidy events, the percentage decreased to 21-65%. Gene co-expression network analyses using de novo assemblies versus mapping to the reference genes recovered the same module that significantly correlated with treatment in one species that lacks recent polyploidization. AVAILABILITY AND IMPLEMENTATION: Commands and scripts used in this study are available at https://bitbucket.org/lychen83/chen_et_al_2018_benchmark_dge/; Analysis files are available at Dryad doi: 10.5061/dryad.4p6n481. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Unique transcriptional signatures of sleep loss across independently evolved cavefish populations.

Animals respond to sleep loss with compensatory rebound sleep, and this is thought to be critical for the maintenance of physiological homeostasis. Sleep duration varies dramatically across animal species, but it is not known whether evolutionary differences in sleep duration are associated with differences in sleep homeostasis. The Mexican cavefish, Astyanax mexicanus, has emerged as a powerful model for studying the evolution of sleep. While eyed surface populations of A. mexicanus sleep approximately 8 hr each day, multiple blind cavefish populations have converged on sleep patterns that total as little as 2 hr each day, providing the opportunity to examine whether the evolution of sleep loss is accompanied by changes in sleep homeostasis. Here, we examine the behavioral and molecular response to sleep deprivation across four independent populations of A. mexicanus. Our behavioral analysis indicates that surface fish and all three cavefish populations display robust recovery sleep during the day following nighttime sleep deprivation, suggesting sleep homeostasis remains intact in cavefish. We profiled transcriptome-wide changes associated with sleep deprivation in surface fish and cavefish. While the total number of differentially expressed genes was not greater for the surface population, the surface population exhibited the highest number of uniquely differentially expressed genes than any other population. Strikingly, a majority of the differentially expressed genes are unique to individual cave populations, suggesting unique expression responses are exhibited across independently evolved cavefish populations. Together, these findings suggest sleep homeostasis is intact in cavefish despite a dramatic reduction in overall sleep duration.

Expression analyses of cave mollies (Poecilia mexicana) reveal key genes involved in the early evolution of eye regression.

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced-albeit functional-eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana, shedding further light on regressive evolution in cavefish.

The evolution of a series of behavioral traits is associated with autism-risk genes in cavefish.

BACKGROUND: An essential question in evolutionary biology is whether shifts in a set of polygenic behaviors share a genetic basis across species. Such a behavioral shift is seen in the cave-dwelling Mexican tetra, Astyanax mexicanus. Relative to surface-dwelling conspecifics, cavefish do not school (asocial), are hyperactive and sleepless, adhere to a particular vibration stimulus (imbalanced attention), behave repetitively, and show elevated stress hormone levels. Interestingly, these traits largely overlap with the core symptoms of human autism spectrum disorder (ASD), raising the possibility that these behavioral traits are underpinned by a similar set of genes (i.e. a repeatedly used suite of genes). RESULT: Here, we explored whether modification of ASD-risk genes underlies cavefish evolution. Transcriptomic analyses revealed that > 58.5% of 3152 cavefish orthologs to ASD-risk genes are significantly up- or down-regulated in the same direction as genes in postmortem brains from ASD patients. Enrichment tests suggest that ASD-risk gene orthologs in A. mexicanus have experienced more positive selection than other genes across the genome. Notably, these positively selected cavefish ASD-risk genes are enriched for pathways involved in gut function, inflammatory diseases, and lipid/energy metabolism, similar to symptoms that frequently coexist in ASD patients. Lastly, ASD drugs mitigated cavefish's ASD-like behaviors, implying shared aspects of neural processing. CONCLUSION: Overall, our study indicates that ASD-risk genes and associated pathways (especially digestive, immune and metabolic pathways) may be repeatedly used for shifts in polygenic behaviors across evolutionary time.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs.

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S-along with related lineages from non-sulfidic environments-to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae).

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish.

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H2 S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H2 S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H2 S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H2 S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H2 S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Multiple mating and reproductive skew in parental and introgressed females of the live-bearing fish Xiphophorus birchmanni.

Just as mating patterns can promote speciation or hybridization, the presence of hybridization can shape mating patterns within a population. In this study, we characterized patterns of multiple mating and reproductive skew in a naturally hybridizing swordtail fish species, Xiphophorus birchmanni. We quantified multiple mating using microsatellite markers to genotype embryos from 43 females collected from 2 wild populations. We also used a suite of single-nucleotide polymorphism markers to categorize females and their inferred mates as either parental X. birchmanni or as introgressed individuals, which carried alleles from a sister species, X. malinche. We found that parental and introgressed X. birchmanni females mated multiply with both parental and introgressed males. We found no difference in mating patterns or reproductive skew between parental and introgressed X. birchmanni females. However, nonintrogressed X. birchmanni males mated more often with large, fecund females. These females also had the greatest levels of skew in fertilization success of males. Thus, our results show that X. birchmanni has a polygynandrous mating system and that introgression of X. malinche alleles has only subtle effects on mating patterns in this species.

Co-evolution of AR gene copy number and structural complexity in endocrine therapy resistant prostate cancer.

Androgen receptor (AR) inhibition is standard of care for advanced prostate cancer (PC). However, efficacy is limited by progression to castration-resistant PC (CRPC), usually due to AR re-activation via mechanisms that include AR amplification and structural rearrangement. These two classes of AR alterations often co-occur in CRPC tumors, but it is unclear whether this reflects intercellular or intracellular heterogeneity of AR. Resolving this is important for developing new therapies and predictive biomarkers. Here, we analyzed 41 CRPC tumors and 6 patient-derived xenografts (PDXs) using linked-read DNA-sequencing, and identified 7 tumors that developed complex, multiply-rearranged AR gene structures in conjunction with very high AR copy number. Analysis of PDX models by optical genome mapping and fluorescence in situ hybridization showed that AR residing on extrachromosomal DNA (ecDNA) was an underlying mechanism, and was associated with elevated levels and diversity of AR expression. This study identifies co-evolution of AR gene copy number and structural complexity via ecDNA as a mechanism associated with endocrine therapy resistance.

Genomic resources for a model in adaptation and speciation research: characterization of the Poecilia mexicana transcriptome.

BACKGROUND: Elucidating the genomic basis of adaptation and speciation is a major challenge in natural systems with large quantities of environmental and phenotypic data, mostly because of the scarcity of genomic resources for non-model organisms. The Atlantic molly (Poecilia mexicana, Poeciliidae) is a small livebearing fish that has been extensively studied for evolutionary ecology research, particularly because this species has repeatedly colonized extreme environments in the form of caves and toxic hydrogen sulfide containing springs. In such extreme environments, populations show strong patterns of adaptive trait divergence and the emergence of reproductive isolation. Here, we used RNA-sequencing to assemble and annotate the first transcriptome of P. mexicana to facilitate ecological genomics studies in the future and aid the identification of genes underlying adaptation and speciation in the system. DESCRIPTION: We provide the first annotated reference transcriptome of P. mexicana. Our transcriptome shows high congruence with other published fish transcriptomes, including that of the guppy, medaka, zebrafish, and stickleback. Transcriptome annotation uncovered the presence of candidate genes relevant in the study of adaptation to extreme environments. We describe general and oxidative stress response genes as well as genes involved in pathways induced by hypoxia or involved in sulfide metabolism. To facilitate future comparative analyses, we also conducted quantitative comparisons between P. mexicana from different river drainages. 106,524 single nucleotide polymorphisms were detected in our dataset, including potential markers that are putatively fixed across drainages. Furthermore, specimens from different drainages exhibited some consistent differences in gene regulation. CONCLUSIONS: Our study provides a valuable genomic resource to study the molecular underpinnings of adaptation to extreme environments in replicated sulfide spring and cave environments. In addition, this study adds to the increasing number of genomic resources in the family Poeciliidae, which are widely used in comparative analyses of behavior, ecology, evolution, and medical genetics.

Contrasting Patterns of Rapid Molecular Evolution within the p53 Network across Mammal and Sauropsid Lineages.

Cancer is a threat to multicellular organisms, yet the molecular evolution of pathways that prevent the accumulation of genetic damage has been largely unexplored. The p53 network regulates how cells respond to DNA-damaging stressors. We know little about p53 network molecular evolution as a whole. In this study, we performed comparative genetic analyses of the p53 network to quantify the number of genes within the network that are rapidly evolving and constrained, and the association between lifespan and the patterns of evolution. Based on our previous published data set, we used genomes and transcriptomes of 34 sauropsids and 32 mammals to analyze the molecular evolution of 45 genes within the p53 network. We found that genes in the network exhibited evidence of positive selection and divergent molecular evolution in mammals and sauropsids. Specifically, we found more evidence of positive selection in sauropsids than mammals, indicating that sauropsids have different targets of selection. In sauropsids, more genes upstream in the network exhibited positive selection, and this observation is driven by positive selection in squamates, which is consistent with previous work showing rapid divergence and adaptation of metabolic and stress pathways in this group. Finally, we identified a negative correlation between maximum lifespan and the number of genes with evidence of divergent molecular evolution, indicating that species with longer lifespans likely experienced less variation in selection across the network. In summary, our study offers evidence that comparative genomic approaches can provide insights into how molecular networks have evolved across diverse species.

Nonrandom RNAseq gene expression associated with RNAlater and flash freezing storage methods.

RNA sequencing is a popular next-generation sequencing technique for assaying genome-wide gene expression profiles. Nonetheless, it is susceptible to biases that are introduced by sample handling prior gene expression measurements. Two of the most common methods for preserving samples in both field-based and laboratory conditions are submersion in RNAlater and flash freezing in liquid nitrogen. Flash freezing in liquid nitrogen can be impractical, particularly for field collections. RNAlater is a solution for stabilizing tissue for longer-term storage as it rapidly permeates tissue to protect cellular RNA. In this study, we assessed genome-wide expression patterns in 30-day-old fry collected from the same brood at the same time point that were flash-frozen in liquid nitrogen and stored at -80°C or submerged and stored in RNAlater at room temperature, simulating conditions of fieldwork. We show that sample storage is a significant factor influencing observed differential gene expression. In particular, genes with elevated GC content exhibit higher observed expression levels in liquid nitrogen flash-freezing relative to RNAlater storage. Further, genes with higher expression in RNAlater relative to liquid nitrogen experience disproportionate enrichment for functional categories, many of which are involved in RNA processing. This suggests that RNAlater may elicit a physiological response that has the potential to bias biological interpretations of expression studies. The biases introduced to observed gene expression arising from mimicking many field-based studies are substantial and should not be ignored.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome.

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

Convergent evolution of reduced energy demands in extremophile fish.

Convergent evolution in organismal function can arise from nonconvergent changes in traits that contribute to that function. Theory predicts that low resource availability and high maintenance costs in extreme environments select for reductions in organismal energy demands, which could be attained through modifications of body size or metabolic rate. We tested for convergence in energy demands and underlying traits by investigating livebearing fish (genus Poecilia) that have repeatedly colonized toxic, hydrogen sulphide-rich springs. We quantified variation in body size and routine metabolism across replicated sulphidic and non-sulphidic populations in nature, modelled total organismal energy demands, and conducted a common-garden experiment to test whether population differences had a genetic basis. Sulphidic populations generally exhibited smaller body sizes and lower routine metabolic rates compared to non-sulphidic populations, which together caused significant reductions in total organismal energy demands in extremophile populations. Although both mechanisms contributed to variation in organismal energy demands, variance partitioning indicated reductions of body size overall had a greater effect than reductions of routine metabolism. Finally, population differences in routine metabolism documented in natural populations were maintained in common-garden reared individuals, indicating evolved differences. In combination with other studies, these results suggest that reductions in energy demands may represent a common theme in adaptation to physiochemical stressors. Selection for reduced energy demand may particularly affect body size, which has implications for life history evolution in extreme environments.

Reduction of Energetic Demands through Modification of Body Size and Routine Metabolic Rates in Extremophile Fish.

Variation in energy availability or maintenance costs in extreme environments can exert selection for efficient energy use, and reductions in organismal energy demand can be achieved in two ways: reducing body mass or metabolic suppression. Whether long-term exposure to extreme environmental conditions drives adaptive shifts in body mass or metabolic rates remains an open question. We studied body size variation and variation in routine metabolic rates in locally adapted populations of extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulfide-rich springs and caves. We quantified size distributions and routine metabolic rates in wild-caught individuals from four habitat types. Compared with ancestral populations in nonsulfidic surface habitats, extremophile populations were characterized by significant reductions in body size. Despite elevated metabolic rates in cave fish, the body size reduction precipitated in significantly reduced energy demands in all extremophile populations. Laboratory experiments on common garden-raised fish indicated that elevated routine metabolic rates in cave fish likely have a genetic basis. The results of this study indicate that adaptation to extreme environments directly impacts energy metabolism, with fish living in cave and sulfide spring environments expending less energy overall during routine metabolism.

Sexual ornaments, body morphology, and swimming performance in naturally hybridizing swordtails (teleostei: xiphophorus).

Determining the costs of sexual ornaments is complicated by the fact that ornaments are often integrated with other, non-sexual traits, making it difficult to dissect the effect of ornaments independent of other aspects of the phenotype. Hybridization can produce reduced phenotypic integration, allowing one to evaluate performance across a broad range of multivariate trait values. Here we assess the relationship between morphology and performance in the swordtails Xiphophorus malinche and X. birchmanni, two naturally-hybridizing fish species that differ extensively in non-sexual as well as sexual traits. We took advantage of novel trait variation in hybrids to determine if sexual ornaments incur a cost in terms of locomotor ability. For both fast-start and endurance swimming, hybrids performed at least as well as the two parental species. The sexually-dimorphic sword did not impair swimming performance per se. Rather, the sword negatively affected performance only when paired with a sub-optimal body shape. Studies seeking to quantify the costs of ornaments should consider that covariance with non-sexual traits may create the spurious appearance of costs.