Species of Anolis lizards of the West Indies that naturally inhabit hot and open areas also tend to thrive in urban areas. In this study, transcriptome was sequenced for nine species of Cuban Anolis lizards that are closely related to each other, but inhabit different thermal microhabitats. Using PAML and HyPhy software, we attempted to identify genes and amino acid sites under positive selection in the common ancestral branch of A. porcatus and A. allisoni, and the branch of A. sagrei, which inhabit hot and open areas, and thrive in urban areas. Although there were no genes where positive selection was commonly detected on both of the tested branches, positive selection was detected in genes involved in the stress response (e.g., DNA damage and oxidative stress) and cardiac function, which could be related to adaptive evolution of tolerance to heat or ultraviolet radiation, on both branches. These findings suggest that adaptive evolution of the response to stress caused by heat or ultraviolet radiation might have occurred in ancestors of Anolis species inhabiting hot and open areas and might be related to the current thriving in urban areas of them.

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.

The detection of various type of 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 the mechanism of the evolution because of the repeated evolution of their morphology and the ecology. In this study, 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). As a result, we obtained six novel draft genomes with relatively long and high gene completeness, with scaffold N50 ranging from 5.56–39.79 Mb, and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Subsequently, we performed comparative analysis of genomic contents including those of mainland Anolis lizards to estimate genetic variations that had emerged and accumulated during the diversification of Anolis lizards. Comparing the repeat element compositions and repeat landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species, LTR accumulation observed only in A. carolinensis, and separate expansions of several families of LINE in each of Cuban trunk-ground species. The analysis of duplicated genes suggested that the proportional difference of duplicated gene number among Cuban Anolis lizards may be associated to the difference of their habitat range. Furthermore, Pairwise Sequentially Markovian Coalescent analysis proposed that the effective population sizes of each species might have been affected by Cuba’s geohistory. Hence, these six novel draft genome assemblies and detected genetic variations can be a springboard for the further genetic elucidation of the Anolis lizard’s diversification.SignificanceAnolis lizard in the West Indies is excellent model for studying the mechanisms of speciation and adaptive evolution. Still, due to a lack of genome assemblies, genetic variations and accumulation process of them involved in the diversification remain largely unexplored. In this study, we reported the novel genome assemblies of six Cuban Anolis lizards and analyzed evolution of genome contents. From comparative genomic analysis and inferences of genetic variation accumulation process, we detected species- and lineage-specific transposon accumulation processes and gene copy number evolution, considered to be associated with the adaptation to their habitats. Additionally, we estimated past effective population sizes and the results suggested its relationship to Cuba’s geohistory.