Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia

BACKGROUND

Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework.

RESULTS

Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions.

CONCLUSIONS

Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

Draft genome of six Cuban Anolis lizards and insights into genetic changes during their diversification

Background

Detecting genomic variants and their accumulation processes during species diversification and adaptive radiation is important for understanding the molecular and genetic basis of evolution. Anolis lizards in the West Indies are good models for studying evolutionary mechanisms because of the repeated evolution of their morphology and the ecology. We performed de novo genome assembly of six Cuban Anolis lizards with different ecomorphs and thermal habitats (Anolis isolepis, Anolis allisoni, Anolis porcatus, Anolis allogus, Anolis homolechis, and Anolis sagrei). We carried out a comparative analysis of these genome assemblies to investigate the genetic changes that occurred during their diversification.

Results

We reconstructed novel draft genomes with relatively long scaffolds and high gene completeness, with the scaffold N50 ranging from 5.56 to 39.79 Mb and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Comparing the repeat element compositions and landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species and separate expansions of several families of LINE in each Cuban trunk-ground species. Duplicated gene analysis suggested that the proportional differences in duplicated gene numbers among Cuban Anolis lizards may be associated with differences in their habitat ranges. Additionally, Pairwise Sequentially Markovian Coalescent analysis suggested that the effective population sizes of each species may have been affected by Cuba’s geohistory.

Conclusions

We provide draft genomes of six Cuban Anolis lizards and detected species and lineage-specific transposon accumulation and gene copy number changes that may be involved in adaptive evolution. The change processes in the past effective population size was also estimated, and the factors involved were inferred. These results provide new insights into the genetic basis of Anolis lizard diversification and are expected to serve as a stepping stone for the further elucidation of their diversification mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02086-7.

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

Phylogenomics and historical biogeography of West Indian Rock Iguanas (genus Cyclura)

The genus Cyclura includes nine extant species and six subspecies of West Indian Rock Iguanas and is one of the most imperiled genera of squamate reptiles globally. An understanding of species diversity, evolutionary relationships, diversification, and historical biogeography in this group is crucial for implementing sound long-term conservation strategies. We collected DNA samples from 1 to 10 individuals per taxon from all Cyclura taxa (n = 70 ingroup individuals), focusing where possible on incorporating individuals from different populations of each species. We also collected 1-2 individuals from each of seven outgroup species of iguanas (Iguana delicatissima; five Ctenosaura species) and Anolis sagrei (n = 12 total outgroup individuals). We used targeted genomic sequence capture to isolate and to sequence 1,872 loci comprising of 687,308 base pairs (bp) from each of the 82 individuals from across the nuclear genome. We extracted mitochondrial reads and assembled and annotated mitogenomes for all Cyclura taxa plus outgroup species. We present well-supported phylogenomic gene tree/species tree analyses for all extant species of Cyclura using ASTRAL-III, SVDQuartets, and StarBEAST2 methods, and discuss the taxonomic, biogeographic, and conservation implications of these data. We find a most recent common ancestor of the genus 9.91 million years ago. The earliest divergence within Cyclura separates C. pinguis from a clade comprising all other Cyclura. Within the latter group, a clade comprising C. carinata from the southern Lucayan Islands and C. ricordii from Hispaniola is the sister taxon to a clade comprising the other Cyclura. Among the other Cyclura, the species C. cornuta and C. stejnegeri (from Hispaniola and Isla Mona) form the sister taxon to a clade of species from Jamaica (C. collei), Cuba and Cayman Islands (C. nubila and C. lewisi), and the eastern (C. rileyi) and western (C. cychlura) Lucayan Islands. Cyclura cychlura and C. rileyi form a clade whose sister taxa are C. nubila and C. lewisi. Cyclura collei is the sister taxon to these four species combined.

Dewlap colour variation in Anolis sagrei is maintained among habitats within islands of the West Indies

Animal signals evolve in an ecological context. Locally adapting animal sexual signals can be especially important for initiating or reinforcing reproductive isolation during the early stages of speciation. Previous studies have demonstrated that dewlap colour in Anolis lizards can be highly variable between populations in relation to both biotic and abiotic adaptive drivers at relatively large geographical scales. Here, we investigated differentiation of dewlap coloration among habitat types at a small spatial scale, within multiple islands of the West Indies, to test the hypothesis that similar local adaptive processes occur over smaller spatial scales. We explored variation in dewlap coloration in the most widespread species of anole, Anolis sagrei, across three characteristic habitats spanning the Bahamas and the Cayman Islands, namely beach scrub, primary coppice forest and mangrove forest. Using reflectance spectrometry paired with supervised machine learning, we found significant differences in spectral properties of the dewlap between habitats within small islands, sometimes over very short distances. Passive divergence in dewlap phenotype associated with isolation-by-distance did not seem to explain our results. On the other hand, these habitat-specific dewlap differences varied in magnitude and direction across islands, and thus, our primary test for adaptation-parallel responses across islands-was not supported. We suggest that neutral processes or selection could be involved in several ways, including sexual selection. Our results shed new light on the scale at which signal colour polymorphism can be maintained in the presence of gene flow, and the relative role of local adaptation and other processes in driving these patterns of dewlap colour variation across islands.

DNA barcoding of the National Museum of Natural History reptile tissue holdings raises concerns about the use of natural history collections and the responsibilities of scientists in the molecular age

Natural history collections are essential to a wide variety of studies in biology because they maintain large collections of specimens and associated data, including genetic material (e.g., tissues) for DNA sequence data, yet they are currently under-funded and collection staff have high workloads. With the advent of aggregate databases and advances in sequencing technologies, there is an increased demand on collection staff for access to tissue samples and associated data. Scientists are rapidly developing large DNA barcode libraries, DNA sequences of specific genes for species across the tree of life, in order to document and conserve biodiversity. In doing so, mistakes are made. For instance, inconsistent taxonomic information is commonly taken from different lending institutions and deposited in data repositories, such as the Barcode of Life Database (BOLD) and GenBank, despite explicit disclaimers regarding the need for taxonomic verification by the lending institutions. Such errors can have profound effects on subsequent research based on these mis-labelled sequences in data repositories. Here, we present the production of a large DNA barcode library of reptiles from the National Museum of Natural History tissue holdings. The library contains 2,758 sequences (2,205 COI and 553 16S) from 2260 specimens (four crocodilians, 37 turtles, and 2,219 lizards, including snakes), representing 583 named species, from 52 countries. In generating this library, we noticed several common mistakes made by scientists depositing DNA barcode data in public repositories (e.g., BOLD and GenBank). Our goal is to raise awareness of these concerns and offer advice to avoid such mistakes in the future to maintain accurate DNA barcode libraries to properly document Earth’s biodiversity.

Phylogenetic inference of where species spread or split across barriers

SignificanceGeography molds how species evolve in space. Strong geographical barriers to movement, for instance, both inhibit dispersal between regions and allow isolated populations to diverge as new species. Weak barriers, by contrast, permit species range expansion and persistence. These factors present a conundrum: How strong must a barrier be before between-region speciation outpaces dispersal? We designed a phylogenetic model of dispersal, extinction, and speciation that allows regional features to influence rates of biogeographic change and applied it to the neotropical radiation of Anolis lizards. Separation by water induces a threefold steeper barrier to movement than equivalent distances over land. Our model will help biologists detect relationships between evolutionary processes and the spatial contexts in which they operate.

Changes in selection pressure can facilitate hybridization during biological invasion in a Cuban lizard

Hybridization is among the evolutionary mechanisms most frequently hypothesized to drive the success of invasive species, in part because hybrids are common in invasive populations. One explanation for this pattern is that biological invasions coincide with a change in selection pressures that limit hybridization in the native range. To investigate this possibility, we studied the introduction of the brown anole (Anolis sagrei) in the southeastern United States. We find that native populations are highly genetically structured. In contrast, all invasive populations show evidence of hybridization among native-range lineages. Temporal sampling in the invasive range spanning 15 y showed that invasive genetic structure has stabilized, indicating that large-scale contemporary gene flow is limited among invasive populations and that hybrid ancestry is maintained. Additionally, our results are consistent with hybrid persistence in invasive populations resulting from changes in natural selection that occurred during invasion. Specifically, we identify a large-effect X chromosome locus associated with variation in limb length, a well-known adaptive trait in anoles, and show that this locus is often under selection in the native range, but rarely so in the invasive range. Moreover, we find that the effect size of alleles at this locus on limb length is much reduced in hybrids among divergent lineages, consistent with epistatic interactions. Thus, in the native range, epistasis manifested in hybrids can strengthen extrinsic postmating isolation. Together, our findings show how a change in natural selection can contribute to an increase in hybridization in invasive populations.

Evolution of island lizards remains a mystery

Lizards that live in the Greater Antilles exploit a large range of skeletal variations to adapt to similar habitats, in defiance of the theory of plasticity-led evolution.

Gene expression of the IGF hormones and IGF binding proteins across time and tissues in a model reptile

The insulin and insulin-like signaling (IIS) network regulates cellular processes including pre- and postnatal growth, cellular development, wound healing, reproduction, and longevity. Despite their importance in the physiology of vertebrates, the study of the specific functions of the top regulators of the IIS network, insulin-like growth factors (IGFs) and IGF binding proteins (IGFBPs), has been mostly limited to a few model organisms. To expand our understanding of this network, we performed quantitative gene expression of IGF hormones in liver and qualitative expression of IGFBPs across tissues and developmental stages in a model reptile, the brown anole lizard (Anolis sagrei). We found that lizards express IGF2 across all life stages (preoviposition embryos to adulthood) and at a higher level than IGF1, which is opposite to patterns seen in laboratory rodents but similar to those seen in humans and other vertebrate models. IGFBP expression was ubiquitous across tissues (brain, gonad, heart, liver, skeletal muscle, tail, and regenerating tail) in adults, apart from IGFBP5, which was variable. These findings provide an essential foundation for further developing the anole lizard as a physiological and biomedical reptile model, as well as expanding our understanding of the function of the IIS network across species.

Establishment of brown anoles (Anolis sagrei) across a southern California county and potential interactions with a native lizard species

The brown anole, Anolis sagrei, is a native species to the Caribbean; however, A. sagrei has invaded multiple parts of the USA, including Florida, Louisiana, Hawai'i and more recently California. The biological impacts of A. sagrei invading California are currently unknown. Evidence from the invasion in Taiwan shows that they spread quickly and when immediate action is not taken eradication stops being a viable option. In Orange County, California, five urban sites, each less than 100 ha, were surveyed for an average of 49.2 min. Approximately 200 A. sagrei were seen and verified across all survey sites. The paucity of native lizards encountered during the surveys within these sites suggests little to no overlap between the dominant diurnal western fence lizard, Sceloporus occidentalis, and A. sagrei. This notable lack of overlap could indicate a potentially disturbing reality that A. sagrei are driving local extirpations of S. occidentalis.