Rapid introgression of non-native alleles following hybridization between a native Anolis lizard species and a cryptic invader across an urban landscape.

Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced-representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non-native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non-native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non-native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non-native genetic material even in the absence of ongoing immigration, indicating that selection favouring non-native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non-native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long-term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.

Maturation and refinement of the maculae and foveae in the Anolis sagrei lizard.

The fovea is a pit in the center of the macula, which is a region of the retina with a high concentration of photoreceptor cells, which accounts for a large degree of visual acuity in primates. The maturation of this primate visual acuity area is characterized by the shallowing and widening of the foveal pit, a decrease in the diameter of the rod-free zone, and an increase in photoreceptor cells packing after birth. Maturation occurs concurrently with progressing age, increasing eye size, and retinal length/area. These observations have led to the hypothesis that the maturation of the fovea might be a function of mechanical variables that remodel the retina. However, this has never been explored outside of primates. Here, we take advantage of the Anolis sagrei lizard, which has a bifoveated retina, to study maturation of the fovea and macula. Eyes were collected from male and female lizards-hatchling, 2-month, 4-month, 6-month, and adult. We found that Anolis maculae undergo a maturation process somewhat different than what has been observed in primates. Anole macular diameters actually increase in size and undergo minimal photoreceptor cell packing, possessing a near complete complement of these cells at the time of hatching. As the anole eye expands, foveal centers experience little change in overall retina cell density with most cell redistribution occurring at macular borders and peripheral retina areas. Gene editing technology has recently been developed in lizards; this study provides a baseline of normal retina maturation for future genetic manipulation studies in anoles.

Wind Resistance Mechanism of an Anole Lizard-Inspired Climbing Robot.

The stable operation of climbing robots exposed to high winds is of great significance for the health-monitoring of structures. This study proposes an anole lizard-like climbing robot inspired by its superior wind resistance. First, the stability mechanism of the anole lizard body in adhesion and desorption is investigated by developing adhesion and desorption models, respectively. Then, the hypothesis that the anole lizard improves its adhesion and stability performance through abdominal adjustment and trunk swing is tested by developing a simplified body model and kinematic model. After that, the structures of the toe, limb, and multi-stage flexible torso of the anole lizard-like climbing robot are designed. Subsequently, the aerodynamic behavior of the proposed robot under high-speed airflow are investigated using finite element analysis. The results show that when there is no obstacle, the climbing robot generates the normal force to enhance toepad friction and adhesion by tuning the abdomen's shape to create an air pressure difference between the back and abdomen. When there is an obstacle, a component force is obtained through periodic oscillation of the spine and tail to resist the frontal winds resulting from the vortex paths generated by the airflow behind the obstacle. These results confirm that the proposed hypothesis is correct. Finally, the adhesion and wind resistance performance of the anole lizard-like climbing robot is tested through the developed experimental platform. It is found that the adhesion force is equal to 50 N when the pre-pressure is 20 N. Further, it is shown that the normal pressure of the proposed robot can reach 76.6% of its weight in a high wind of 14 m/s.

Ocular elongation and retraction in foveated reptiles.

BACKGROUND: Pronounced asymmetric changes in ocular globe size during eye development have been observed in a number of species ranging from humans to lizards. In contrast, largely symmetric changes in globe size have been described for other species like rodents. We propose that asymmetric changes in the three-dimensional structure of the developing eye correlate with the types of retinal remodeling needed to produce areas of high photoreceptor density. To test this idea, we systematically examined three-dimensional aspects of globe size as a function of eye development in the bifoveated brown anole, Anolis sagrei. RESULTS: During embryonic development, the anole eye undergoes dynamic changes in ocular shape. Initially spherical, the eye elongates in the presumptive foveal regions of the retina and then proceeds through a period of retraction that returns the eye to its spherical shape. During this period of retraction, pit formation and photoreceptor cell packing are observed. We found a similar pattern of elongation and retraction associated with the single fovea of the veiled chameleon, Chamaeleo calyptratus. CONCLUSIONS: These results, together with those reported for other foveated species, support the idea that areas of high photoreceptor packing occur in regions where the ocular globe asymmetrically elongates and retracts during development.

Genome-wide epigenetic isolation by environment in a widespread Anolis lizard.

Epigenetic changes can provide a pathway for organisms to respond to local environmental conditions by influencing gene expression. However, we still know little about the spatial distribution of epigenetic variation in natural systems, how it relates to the distribution of genetic variation and the environmental structure of the landscape, and the processes that generate and maintain it. Studies examining spatial patterns of genetic and epigenetic variation can provide valuable insights into how ecological and population processes contribute to epigenetic divergence across heterogeneous landscapes. Here, we perform a comparative analysis of spatial genetic and epigenetic variation based on 8,459 single nucleotide polymorphisms (SNPs) and 8,580 single methylation variants (SMVs) from eight populations of the Puerto Rican crested anole, Anolis cristatellus, an abundant lizard in the adaptive radiations of anoles on the Greater Antilles that occupies a diverse range of habitats. Using generalized dissimilarity modelling and multiple matrix regression, we found that genome-wide epigenetic differentiation is strongly correlated with environmental divergence, even after controlling for the underlying genetic structure. We also detected significant associations between key environmental variables and 96 SMVs, including 42 located in promoter regions or gene bodies. Our results suggest an environmental basis for population-level epigenetic differentiation in this system and contribute to better understanding how environmental gradients structure epigenetic variation in nature.

Squamate egg tooth development revisited using three-dimensional reconstructions of brown anole (Anolis sagrei, Squamata, Dactyloidae) dentition.

The egg tooth is a hatching adaptation, characteristic of all squamates. In brown anole embryos, the first tooth that starts differentiating is the egg tooth. It develops from a single tooth germ and, similar to the regular dentition of all the other vertebrates, the differentiating egg tooth of the brown anole passes through classic morphological and developmental stages named according to the shape of the dental epithelium: epithelial thickening, dental lamina, tooth bud, cap and bell stages. The differentiating egg tooth consists of three parts: the enamel organ, hard tissues and dental pulp. Shortly before hatching, the egg tooth connects with the premaxilla. Attachment tissue of the egg tooth does not undergo mineralization, which makes it different from the other teeth of most squamates. After hatching, odontoclasts are involved in resorption of the egg tooth's remains. This study shows that the brown anole egg tooth does not completely conform to previous reports describing iguanomorph egg teeth and reveals a need to investigate its development in the context of squamate phylogeny.

Ectoparasite extinction in simplified lizard assemblages during experimental island invasion.

Introduced species can become invasive, damaging ecosystems and disrupting economies through explosive population growth. One mechanism underlying population expansion in invasive populations is 'enemy release', whereby the invader experiences relaxation of agonistic interactions with other species, including parasites. However, direct observational evidence of release from parasitism during invasion is rare. We mimicked the early stages of invasion by experimentally translocating populations of mite-parasitized slender anole lizards (Anolis apletophallus) to islands that varied in the number of native anoles. Two islands were anole-free prior to the introduction, whereas a third island had a resident population of Gaige's anole (Anolis gaigei). We then characterized changes in trombiculid mite parasitism over multiple generations post-introduction. We found that mites rapidly went extinct on one-species islands, but that lizards introduced to the two-species island retained mites. After three generations, the two-species island had the highest total density and biomass of lizards, but the lowest density of the introduced species, implying that the 'invasion' had been less successful. This field-transplant study suggests that native species can be 'enemy reservoirs' that facilitate co-colonization of ectoparasites with the invasive host. Broadly, these results indicate that the presence of intact and diverse native communities may help to curb invasiveness.

First report on B chromosome content in a reptilian species: the case of Anolis carolinensis.

Supernumerary elements of the genome are often called B chromosomes. They usually consist of various autosomal sequences and, because of low selective pressure, are mostly pseudogenized and contain many repeats. There are numerous reports on B chromosomes in mammals, fish, invertebrates, plants, and fungi, but only a few of them have been studied using sequencing techniques. However, reptilian supernumerary chromosomes have been detected only cytogenetically and never sequenced or analyzed at the molecular level. One model squamate species with available genome sequence is Anolis carolinensis. The scope of the present article is to describe the genetic content of A. carolinensis supernumerary chromosomes. In this article, we confirm the presence of B chromosomes in this species by reverse painting and synaptonemal complex analysis. We applied low-pass high-throughput sequencing to analyze flow-sorted B chromosomes. Anole B chromosomes exhibit similar traits to other supernumerary chromosomes from different taxons: they contain two genes related to cell division control (INCENP and SPIRE2), are enriched in specific repeats, and show a high degree of pseudogenization. Therefore, the present study confirms that reptilian B chromosomes resemble supernumerary chromosomes of other taxons.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer.

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including FST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Estimating encounter rates as the first step of sexual selection in the lizard Anolis sagrei.

How individuals move through their environment dictates which other individuals they encounter, determining their social and reproductive interactions and the extent to which they experience sexual selection. Specifically, females rarely have the option of mating with all males in a population-they can only choose among the males they encounter. Further, quantifying phenotypic differences between the males that females encounter and those that sire females' offspring lends insight into how social and reproductive interactions shape male phenotypes. We used an explicitly spatio-temporal Markov chain model to estimate the number of potential mates of Anolis sagrei lizards from their movement behaviour, and used genetic paternity assignments to quantify sexual selection on males. Females frequently encountered and mated with multiple males, offering ample opportunity for female mate choice. Sexual selection favoured males that were bigger and moved over larger areas, though the effect of body size cannot be disentangled from last-male precedence. Our approach corroborates some patterns of sexual selection previously hypothesized in anoles based on describing them as territorial, whereas other results, including female multiple mating itself, are at odds with territorial polygyny, offering insight into discrepancies in other taxa between behavioural and genetic descriptions of mating systems.

When did anoles diverge? An analysis of multiple dating strategies.

Whereas most of the studies that discuss the evolutionary divergence of Anolis lizards have dated the clade's crown group in between 31 and 64 Ma, a single study has recovered a significantly older age for the same node (87 Ma). These differences also entail notable consequences on the preferred biogeographical hypothesis for the whole clade. Here we analyze a total of seven dating strategies by combining three calibration sources in independent BEAST runs to infer the most probable divergence timing for anole lizards (a mitochondrial rate for ND2 gene, the Anolis dominicanus fossil, and a group of fossils assigned to the Priscagamines, Iguanines, and Idontosaurus clades). Based on the estimated timing, we also addressed whether chronograms differ the most in deeper or shallower nodes by exploring the trend in the standard deviation of mean ages between chronograms across time. Next, we focus on the pattern for a single shallow node by hypothesizing the biogeography of the island-endemic Malpelo anole (Anolis agassizi), and evaluating the temporal congruence between the species' divergence and the island geology. The estimated set of ages suggests that anoles most likely diverged 72 Ma (71-73 Ma), with the crown group established around 58 Ma (51-65 Ma). Dispersal is therefore supported as the major driver in the biogeography of the group (and in Caribbean lineages in particular). Our analyses also indicated that (1) rate-based analyses pulled dates toward younger ages, (2) the differences in node ages between chronograms decrease towards the tips regardless of the position of the constrained node, and that (3) the estimated age for deep nodes (e.g. Anolis stem) is highly influenced when deep nodes are also constrained. The latter two results imply that the estimated age for shallower nodes is largely unaffected by the used temporal constraint. The congruence of all chronograms for the Malpelo anole also supports this finding. Anolis agassizi was found to have diverged before the emergence of Malpelo island in each analysis (anole: 19-31 Ma vs. Malpelo island: 16-17 Ma). Finally, we recommend when performing absolute dating analyses to first test for sequence saturation in the analyzed dataset (especially when calibrations are based on molecular rates). Our study also points out the importance of using multiple node constraints, especially when placed deeply in the tree, for fossil-based divergence dating analyses.

Anolis carolinensis as a model to understand the molecular and cellular basis of foveal development.

The fovea is an anatomical specialization of the central retina containing closely packed cone-photoreceptors providing an area of high acuity vision in humans and primates. Despite its key role in the clarity of vision, little is known about the molecular and cellular basis of foveal development, due to the absence of a foveal structure in commonly used laboratory animal models. Of the amniotes the retina in birds of prey and some reptiles do exhibit a typical foveal structure, but they have not been studied in the context of foveal development due to lack of availability of embryonic tissue, lack of captive breeding programs, and limited genomic information. However, the genome for the diurnal bifoveate reptile species Anolis carolinensis (green anole) was recently published and it is possible to collect embryos from this species in captivity. Here, we tested the feasibility of using the anole as a model to study foveal development. Eyes were collected at various stages of development for histological analysis, immunofluorescence, and apoptosis. We show that at embryonic stage (ES) 10 there is peak ganglion cell density at the incipient central foveal region and a single row of cone photoreceptor nuclei. At ES17 the foveal pit begins to form and at this stage there are 3-4 rows of cone nuclei. Post-hatching a further increase in cone density and lengthening of inner and outer segments is observed. A yellowish pigment was seen in the adult central foveal region, but not in the temporal fovea. At ES14 Pax6 was localized across the entire retina, but was more prominent in the ganglion cell layer (GCL) and the part of the inner nuclear layer (INL) containing amacrine cell bodies. However, at ES17 Pax6 expression in the ganglion cells of the central retina was markedly reduced. Bioinformatic analysis revealed that 86% of human candidate foveal hypoplasia genes had an orthologous gene or DNA sequence in the green anole. These findings provide the first insight into foveal morphogenesis in the green anole and suggest that it could be a very useful model for investigating the molecular signals driving foveal development, and thus inform on human foveal development and disease.

Identification and characterization of tyrosine kinases in anole lizard indicate the conserved tyrosine kinase repertoire in vertebrates.

The tyrosine kinases (TKs) play principal roles in regulation of multicellular aspects of the organism and are implicated in many cancer types and congenital disorders. The anole lizard has recently been introduced as a model organism for laboratory-based studies of organismal function and field studies of ecology and evolution. However, the TK family of anole lizard has not been systematically identified and characterized yet. In this study, we identified 82 TK-encoding genes in the anole lizard genome and classified them into 28 subfamilies through phylogenetic analysis, with no member from ROS and STYK1 subfamilies identified. Although TK domain sequences and domain organization in each subfamily were conserved, the total number of TKs in different species was much variable. In addition, extensive evolutionary analysis in metazoans indicated that TK repertoire in vertebrates tends to be remarkably stable. Phylogenetic analysis of Eph subfamily indicated that the divergence of EphA and EphB occurred prior to the whole genome duplication (WGD) but after the split of Urochordates and vertebrates. Moreover, the expression pattern analysis of lizard TK genes among 9 different tissues showed that 14 TK genes exhibited tissue-specific expression and 6 TK genes were widely expressed. Comparative analysis of TK expression suggested that the tissue specifically expressed genes showed different expression pattern but the widely expressed genes showed similar pattern between anole lizard and human. These results may provide insights into the evolutionary diversification of animal TK genes and would aid future studies on TK protein regulation of key growth and developmental processes.

Climate-related environmental variation in a visual signalling device: the male and female dewlap in Anolis sagrei lizards.

Animals communicate using a variety of signals that differ dramatically among and within species. The astonishing dewlap diversity in anoles has attracted considerable attention in this respect. Yet, the evolutionary processes behind it remain elusive and have mostly been explored for males only. Here, we considered Anolis sagrei males and females to study signal divergence among populations. First, we assessed the degree of variation in dewlap design (size, pattern and colour) and displays by comparing 17 populations distributed across the Caribbean. Second, we assessed whether the observed dewlap diversity is associated with variation in climate-related environmental conditions. Results showed that populations differed in all dewlap characteristics, with the exception of display rate in females. We further found that males and females occurring in 'xeric' environments had a higher proportion of solid dewlaps with higher UV reflectance. In addition, lizards inhabiting 'mesic' environments had primarily marginal dewlaps showing high reflectance in red. For dewlap display, a correlation with environment was only observed in males. Our study provides evidence for a strong relationship between signal design and prevailing environmental conditions, which may result from differential selection on signal efficacy. Moreover, our study highlights the importance of including females when studying dewlaps in an evolutionary context.

Exercise training reveals trade-offs between endurance performance and immune function, but does not influence growth, in juvenile lizards.

Acquired energetic resources allocated to a particular trait cannot then be re-allocated to a different trait. This often results in a trade-off between survival and reproduction for the adults of many species, but such a trade-off may be manifested differently in juveniles not yet capable of reproduction. Whereas adults may allocate resources to current and/or future reproduction, juveniles can only allocate to future reproduction. Thus, juveniles should allocate resources toward traits that increase survival and their chances of future reproductive success. We manipulated allocation of resources to performance, via endurance exercise training, to examine trade-offs among endurance capacity, immune function and growth in juvenile green anole lizards. We trained male and female captive anoles on a treadmill for 8 weeks, with increasing intensity, and compared traits with those of untrained individuals. Our results show that training enhanced endurance capacity equally in both sexes, but immune function was suppressed only in females. Training had no effect on growth, but males had higher growth rates than females. Previous work showed that trained adults have enhanced growth, so juvenile growth is either insensitive to stimulation with exercise, or they are already growing at maximal rates. Our results add to a growing body of literature indicating that locomotor performance is an important part of life-history trade-offs that are sex and age specific.

Microsatellite Analysis of the Population Genetic Structure of Anolis carolinensis Introduced to the Ogasawara Islands.

DNA analysis can reveal the origins and dispersal patterns of invasive species. The green anole Anolis carolinensis is one such alien animal, which has been dispersed widely by humans from its native North America to many Pacific Ocean islands. In the Ogasawara (Bonin) Islands, this anole was recorded from Chichi-jima at the end of the 1960s, and then from Haha-jima in the early 1980s. These two islands are inhabited. In 2013, it was also found on the uninhabited Ani-jima, close to Chichi-jima. Humans are thought to have introduced the anole to Haha-jima, while the mode of introduction to Ani-jima is unknown. To clarify its dispersal patterns within and among these three islands, we assessed the fine-scale population genetic structure using five microsatellite loci. The results show a homogeneous genetic structure within islands, but different genetic structures among islands, suggesting that limited gene flow occurs between islands. The recently established Ani-jima population may have originated from several individuals simultaneously, or by repeated immigration from Chichi-jima. We must consider frequent incursions among these islands to control these invasive lizard populations and prevent their negative impact on native biodiversity.

Involvement of different mesotocin (oxytocin homologue) populations in sexual and aggressive behaviours of the brown anole.

The oxytocin (OT) family of neuropeptides are known to modulate social behaviours and anxiety in mammals and birds. We investigated cell numbers and neural activity, assessed as Fos induction, within magnocellular and parvocellular populations of neurons producing the OT homologue mesotocin (MT, Ile(8)-oxytocin). This was conducted within the male brown anole lizard, Anolis sagrei, following agonistic or courtship encounters with a conspecific. Both neurons colocalizing and not colocalizing corticotropin-releasing factor (CRF) were examined. Parvocellular neurons of the paraventricular nucleus exhibited a positive correlation between courtship frequency and Fos colocalization, regardless of whether they produce just MT or MT + CRF. Magnocellular populations showed only trends towards positive relationships with courtship and no cell populations showed aggression-related Fos induction. These findings are novel because they demonstrate the involvement of MT neurons in male social behaviour, especially in reptiles for whom the involvement of MT in social behaviour was previously unknown.

Comparative proteomic analysis of tail regeneration in the green anole lizard, Anolis carolinensis.

As amniote vertebrates, lizards are the most closely related organisms to humans capable of appendage regeneration. Lizards can autotomize, or release their tails as a means of predator evasion, and subsequently regenerate a functional replacement. Green anoles (Anolis carolinensis) can regenerate their tails through a process that involves differential expression of hundreds of genes, which has previously been analyzed by transcriptomic and microRNA analysis. To investigate protein expression in regenerating tissue, we performed whole proteomic analysis of regenerating tail tip and base. This is the first proteomic data set available for any anole lizard. We identified a total of 2,646 proteins - 976 proteins only in the regenerating tail base, 796 only in the tail tip, and 874 in both tip and base. For over 90% of these proteins in these tissues, we were able to assign a clear orthology to gene models in either the Ensembl or NCBI databases. For 13 proteins in the tail base, 9 proteins in the tail tip, and 10 proteins in both regions, the gene model in Ensembl and NCBI matched an uncharacterized protein, confirming that these predictions are present in the proteome. Ontology and pathways analysis of proteins expressed in the regenerating tail base identified categories including actin filament-based process, ncRNA metabolism, regulation of phosphatase activity, small GTPase mediated signal transduction, and cellular component organization or biogenesis. Analysis of proteins expressed in the tail tip identified categories including regulation of organelle organization, regulation of protein localization, ubiquitin-dependent protein catabolism, small GTPase mediated signal transduction, morphogenesis of epithelium, and regulation of biological quality. These proteomic findings confirm pathways and gene families activated in tail regeneration in the green anole as well as identify uncharacterized proteins whose role in regrowth remains to be revealed. This study demonstrates the insights that are possible from the integration of proteomic and transcriptomic data in tail regrowth in the green anole, with potentially broader application to studies in other regenerative models.

On the role of TFEC in reptilian coloration.

Reptilian species, particularly snakes and lizards, are emerging models of animal coloration. Here, I focus on the role of the TFEC transcription factor in snake and lizard coloration based on a study on wild-type and piebald ball pythons. Genomic mapping previously identified a TFEC mutation linked to the piebald ball python phenotype. The association of TFEC with skin coloration was further supported by gene-editing experiments in the brown anole lizard. However, novel histological analyses presented here reveal discrepancies between the ball python and the anole TFEC mutants phenotype, cautioning against broad generalizations. Indeed, both wild-type and piebald ball pythons completely lack iridophores, whereas the TFEC anole lizard mutants lose their iridophores compared to the wild-type anole. Based on these findings, I discuss the potential role of the MiT/TFE family in skin pigmentation across vertebrate lineages and advocate the need for developmental analyses and additional gene-editing experiments to explore the reptilian coloration diversity.

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.