Genome-wide parallelism underlies contemporary adaptation in urban lizards

Alternative splicing in parallel evolution and the evolutionary potential in sticklebacks

Repeatability of adaptation to similar environments provides opportunity to evaluate the predictability of natural selection. While many studies have investigated gene expression differences between populations adapted to contrasting environments, the role of post-transcriptional processes such as alternative splicing has rarely been evaluated in the context of parallel adaptation. To address the aforementioned knowledge gap, we reanalysed transcriptomic data from three pairs of threespine stickleback (Gasterosteus aculeatus) ecotypes adapted to marine or freshwater environment. First, we identified genes with repeated expression or splicing divergence across ecotype pairs, and compared the genetic architecture and biological processes between parallelly expressed and parallelly spliced loci. Second, we analysed the extent to which parallel adaptation was reflected at gene expression and alternative splicing levels. Finally, we tested how the two axes of transcriptional variation differed in their potential for evolutionary change. Although both repeated differential splicing and differential expression across ecotype pairs showed tendency for parallel divergence, the degree of parallelism was lower for splicing than expression. Furthermore, parallel divergences in splicing and expression were likely to be associated with distinct cis-regulatory genetic variants and functionally unique set of genes. Finally, we found that parallelly spliced genes showed higher nucleotide diversity than parallelly expressed genes, indicating splicing is less susceptible to genetic variation erosion during parallel adaptation. Our results provide novel insight into the role of splicing in parallel adaptation, and underscore the contribution of splicing to the evolutionary potential of wild populations under environmental change.

The skull of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae)

We provide a detailed and first description of the skull, hyoid apparatus, and trachea of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae). Cyclura is a radiation of iguanas restricted to islands of the Caribbean Sea. Species of Cyclura have high rates of endemism, and all species are severely threatened with extinction. Our anatomical description of this threatened iguana is based on high-resolution computed tomography scans of one adult, one putative adult or near adult, and one juvenile specimen, and includes three-dimensional segmented renderings and visualizations. We discuss some observations of intraspecific and ontogenetic variation, and provide a brief comparison with specimens of another species of Cyclura and published descriptions of other iguanas. Our study provides a cranial osteological framework for Cyclura and augments the body of knowledge on iguana anatomy generally. Finally, we posit that our description and future studies may facilitate identification of fossil Cyclura, which could help understand the paleobiogeography of the genus.

Online toolkits for collaborative and inclusive global research in urban evolutionary ecology

Urban evolutionary ecology is inherently interdisciplinary. Moreover, it is a field with global significance. However, bringing researchers and resources together across fields and countries is challenging. Therefore, an online collaborative research hub, where common methods and best practices are shared among scientists from diverse geographic, ethnic, and career backgrounds would make research focused on urban evolutionary ecology more inclusive. Here, we describe a freely available online research hub for toolkits that facilitate global research in urban evolutionary ecology. We provide rationales and descriptions of toolkits for: (1) decolonizing urban evolutionary ecology; (2) identifying and fostering international collaborative partnerships; (3) common methods and freely-available datasets for trait mapping across cities; (4) common methods and freely-available datasets for cross-city evolutionary ecology experiments; and (5) best practices and freely available resources for public outreach and communication of research findings in urban evolutionary ecology. We outline how the toolkits can be accessed, archived, and modified over time in order to sustain long-term global research that will advance our understanding of urban evolutionary ecology.

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

Novel Nuclear Roles for Testis-Specific ACTL7A and ACTL7B Supported by In Vivo Characterizations and AI Facilitated In Silico Mechanistic Modeling with Implications for Epigenetic Regulation in Spermiogenesis

A mechanistic role for nuclear function of testis-specific actin related proteins (ARPs) is proposed here through contributions of ARP subunit swapping in canonical chromatin regulatory complexes. This is significant to our understanding of both mechanisms controlling regulation of spermiogenesis, and the expanding functional roles of the ARPs in cell biology. Among these roles, actins and ARPs are pivotal not only in cytoskeletal regulation, but also in intranuclear chromatin organization, influencing gene regulation and nucleosome remodeling. This study focuses on two testis-specific ARPs, ACTL7A and ACTL7B, exploring their intranuclear activities and broader implications utilizing combined in vivo, in vitro, and in silico approaches. ACTL7A and ACTL7B, previously associated with structural roles, are hypothesized here to serve in chromatin regulation during germline development. This study confirms the intranuclear presence of ACTL7B in spermatocytes and round spermatids, revealing a potential role in intranuclear processes, and identifies a putative nuclear localization sequence conserved across mammalian ACTL7B, indicating a potentially unique mode of nuclear transport which differs from conventional actin. Ablation of ACTL7B leads to varied transcriptional changes reported here. Additionally, in the absence of ACTL7A or ACTL7B there is a loss of intranuclear localization of HDAC1 and HDAC3, which are known regulators of epigenetic associated acetylation changes that in turn regulate gene expression. Thus, these HDACs are implicated as contributors to the aberrant gene expression observed in the KO mouse testis transcriptomic analysis. Furthermore, this study employed and confirmed the accuracy of in silico models to predict ARP interactions with Helicase-SANT-associated (HSA) domains, uncovering putative roles for testis-specific ARPs in nucleosome remodeling complexes. In these models, ACTL7A and ACTL7B were found capable of binding to INO80 and SWI/SNF nucleosome remodeler family members in a manner akin to nuclear actin and ACTL6A. These models thus implicate germline-specific ARP subunit swapping within chromatin regulatory complexes as a potential regulatory mechanism for chromatin and associated molecular machinery adaptations in nuclear reorganizations required during spermiogenesis. These results hold implications for male fertility and epigenetic programing in the male-germline that warrant significant future investigation. In summary, this study reveals that ACTL7A and ACTL7B play intranuclear gene regulation roles in male gametogenesis, adding to the multifaceted roles identified also spanning structural, acrosomal, and flagellar stability. ACTL7A and ACTL7B unique nuclear transport, impact on HDAC nuclear associations, impact on transcriptional processes, and proposed mechanism for involvement in nucleosome remodeling complexes supported by AI facilitated in silico modeling contribute to a more comprehensive understanding of the indispensable functions of ARPs broadly in cell biology, and specifically in male fertility.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Genetic depletion does not prevent rapid evolution in island-introduced lizards

Experimental introductions of species have provided some of the most tractable examples of rapid phenotypic changes, which may reflect plasticity, the impact of stochastic processes, or the action of natural selection. Yet to date, very few studies have investigated the neutral and potentially adaptive genetic impacts of experimental introductions. We dissect the role of these processes in shaping the population differentiation of wall lizards in three Croatian islands (Sušac, Pod Kopište, and Pod Mrčaru), including the islet of Pod Mrčaru, where experimentally introduced lizards underwent rapid (~30 generations) phenotypic changes associated with a shift from an insectivorous to a plant-based diet. Using a genomic approach (~82,000 ddRAD loci), we confirmed a founder effect during introduction and very low neutral genetic differentiation between the introduced population and its source. However, genetic depletion did not prevent rapid population growth, as the introduced lizards exhibited population genetic signals of expansion and are known to have reached a high density. Our genome-scan analysis identified just a handful of loci showing large allelic shifts between ecologically divergent populations. This low overall signal of selection suggests that the extreme phenotypic differences observed among populations are determined by a small number of large-effect loci and/or that phenotypic plasticity plays a major role in phenotypic changes. Nonetheless, functional annotation of the outlier loci revealed some candidate genes relevant to diet-induced adaptation, in agreement with the hypothesis of directional selection. Our study provides important insights on the evolutionary potential of bottlenecked populations in response to new selective pressures on short ecological timescales.

The genomic response to urbanization in the damselfly Ischnura elegans

The complex and rapid environmental changes brought about by urbanization pose significant challenges to organisms. The multifaceted effects of urbanization often make it difficult to define and pinpoint the very nature of adaptive urban phenotypes. In such situations, scanning genomes for regions differentiated between urban and non-urban populations may be an attractive approach. Here, we investigated the genomic signatures of adaptation to urbanization in the damselfly Ischnura elegans sampled from 31 rural and urban localities in three geographic regions: southern and northern Poland, and southern Sweden. Genome-wide variation was assessed using more than 370,000 single nucleotide polymorphisms (SNPs) genotyped by ddRADseq. Associations between SNPs and the level of urbanization were tested using two genetic environment association methods: Latent Factors Mixed Models and BayPass. While we found numerous candidate SNPs and a highly significant overlap between candidates identified by the two methods within the geographic regions, there was a distinctive lack of repeatability between the geographic regions both at the level of individual SNPs and of genomic regions. However, we found "synapse organization" at the top of the functional categories enriched among the genes located in the proximity of the candidate urbanization SNPs. Interestingly, the overall significance of "synapse organization" was built up by the accretion of different genes associated with candidate SNPs in different geographic regions. This finding is consistent with the highly polygenic nature of adaptation, where the response may be achieved through a subtle adjustment of allele frequencies in different genes that contribute to adaptive phenotypes. Taken together, our results point to a polygenic adaptive response in the nervous system, specifically implicating genes involved in synapse organization, which mirrors the findings from several genomic and behavioral studies of adaptation to urbanization in other taxa.

Habitat-associated Genomic Variation in a Wall Lizard from an Oceanic Island

The lizard Teira dugesii exhibits morphological divergence between beach and inland habitats in the face of gene flow, within the volcanic island of Madeira, Portugal. Here, we analyzed genomic data obtained by genotyping-by-sequencing, which provided 16,378 single nucleotide polymorphisms (SNPs) from 94 individuals sampled from 15 sites across Madeira. Ancient within-island divergence in allopatry appears to have mediated divergence in similar species within other Atlantic islands, but this hypothesis was not supported for T. dugesii. Across all samples, a total of 168 SNPs were classified as statistical outliers using pcadapt and OutFLANK. Redundancy analysis (RDA) revealed that 17 of these outliers were associated with beach/inland habitats. The SNPs were located within 16 sequence tags and 15 of these were homologous with sequences in a 31 Mb region on chromosome 3 of a reference wall lizard genome (the remaining tag could not be associated with any chromosome). We further investigated outliers through contingency analyses of allele frequencies at each of four pairs of adjacent beach-inland sites. The majority of the outliers detected by the RDA were confirmed at two pairs of these matched sites. These analyses also suggested some parallel divergence at different localities. Six other outliers were associated with site elevation, four of which were located on chromosome 5 of the reference genome. Our study lends support to a previous hypothesis that divergent selection between gray shingle beaches and inland regions overcomes gene flow and leads to the observed morphological divergence between populations in these adjacent habitats.

Variation in Pheidole nodus (Hymenoptera: Formicidae) functional morphology across urban parks

Background

Habitat fragmentation and consequent population isolation in urban areas can impose significant selection pressures on individuals and species confined to urban islands, such as parks. Despite many comparative studies on the diversity and structure of ant community living in urban areas, studies on ants' responses to these highly variable ecosystems are often based on assemblage composition and interspecific mean trait values, which ignore the potential for high intraspecific functional trait variation among individuals.

Methods

Here, we examined differences in functional traits among populations of the generalist ant Pheidole nodus fragmented between urban parks. We used pitfall trapping, which is more random and objective than sampling colonies directly, despite a trade-off against sample size. We then tested whether trait-filtering could explain phenotypic differences among urban park ant populations, and whether ant populations in different parks exhibited different phenotypic optima, leading to positional shifts in anatomical morphospace through the regional ant meta-population.

Results

Intraspecific morphological differentiation was evident across this urban region. Populations had different convex hull volumes, positioned differently over the morphospace.

Conclusions

Fragmentation and habitat degradation reduced phenotypic diversity and, ultimately, changed the morphological optima of populations in this urban landscape. Considering ants' broad taxonomic and functional diversity and their important role in ecosystems, further work over a variety of ant taxa is necessary to ascertain those varied morphological response pathways operating in response to population segregation in urban environments.

Biodiversity: What makes a city bird?

The impacts of urbanisation on biodiversity varies greatly across species. A new study shows how the intrinsic species properties underlying urban tolerance vary globally according to environmental context. This has important implications for conserving biodiversity in a rapidly urbanising world.

Urban evolutionary ecology brings exaptation back into focus

The contribution of pre-existing phenotypic variation to evolution in novel environments has long been appreciated. Nevertheless, evolutionary ecologists have struggled with communicating these aspects of the adaptive process. In 1982, Gould and Vrba proposed terminology to distinguish character states shaped via natural selection for the roles they currently serve ('adaptations') from those shaped under preceding selective regimes ('exaptations'), with the intention of replacing the inaccurate 'preadaptation'. Forty years later, we revisit Gould and Vrba's ideas which, while often controversial, continue to be widely debated and highly cited. We use the recent emergence of urban evolutionary ecology as a timely opportunity to reintroduce the ideas of Gould and Vrba as an integrated framework to understand contemporary evolution in novel environments.

Genome-wide parallelism underlies contemporary adaptation in urban lizards

Urbanization drastically transforms landscapes, resulting in fragmentation, degradation, and the loss of local biodiversity. Yet, urban environments also offer opportunities to observe rapid evolutionary change in wild populations that survive and even thrive in these novel habitats. In many ways, cities represent replicated "natural experiments" in which geographically separated populations adaptively respond to similar selection pressures over rapid evolutionary timescales. Little is known, however, about the genetic basis of adaptive phenotypic differentiation in urban populations nor the extent to which phenotypic parallelism is reflected at the genomic level with signatures of parallel selection. Here, we analyzed the genomic underpinnings of parallel urban-associated phenotypic change in Anolis cristatellus, a small-bodied neotropical lizard found abundantly in both urbanized and forested environments. We show that phenotypic parallelism in response to parallel urban environmental change is underlain by genomic parallelism and identify candidate loci across the Anolis genome associated with this adaptive morphological divergence. Our findings point to polygenic selection on standing genetic variation as a key process to effectuate rapid morphological adaptation. Identified candidate loci represent several functions associated with skeletomuscular development, morphology, and human disease. Taken together, these results shed light on the genomic basis of complex morphological adaptations, provide insight into the role of contingency and determinism in adaptation to novel environments, and underscore the value of urban environments to address fundamental evolutionary questions.