Multi-locus phylogeographic and population genetic analysis of Anolis carolinensis: historical demography of a genomic model species

High-Density Geometric Morphometric Analysis of Intraspecific Cranial Integration in the Barred Grass Snake (Natrix helvetica) and Green Anole (Anolis carolinensis)

How do phenotypic associations intrinsic to an organism, such as developmental and mechanical processes, direct morphological evolution? Comparisons of intraspecific and clade-wide patterns of phenotypic covariation could inform how population-level trends ultimately dictate macroevolutionary changes. However, most studies have focused on analyzing integration and modularity either at macroevolutionary or intraspecific levels, without a shared analytical framework unifying these temporal scales. In this study, we investigate the intraspecific patterns of cranial integration in two squamate species: Natrix helvetica and Anolis carolinensis. We analyze their cranial integration patterns using the same high-density three-dimensional geometric morphometric approach used in a prior squamate-wide evolutionary study. Our results indicate that Natrix and Anolis exhibit shared intraspecific cranial integration patterns, with some differences, including a more integrated rostrum in the latter. Notably, these differences in intraspecific patterns correspond to their respective interspecific patterns in snakes and lizards, with few exceptions. These results suggest that interspecific patterns of cranial integration reflect intraspecific patterns. Hence, our study suggests that the phenotypic associations that direct morphological variation within species extend across micro- and macroevolutionary levels, bridging these two scales.

Climate shapes patterns of sexual size and shape dimorphism across the native range of the green anole lizard, Anolis carolinensis (Squamata: Dactyloidae)

Geographical variation in sexual size dimorphism (SSD) can result from the combined effects of environmental and sexual selection. To understand the determinants of SSD across geographical landscapes, we tested for relationships between SSD and climatic variables in the widespread lizard Anolis carolinensis. To distinguish alternative hypotheses for observed patterns of variation in SSD, we also examined sex-specific patterns of body size evolution and asked whether SSD was associated with certain patterns of sexual shape dimorphism. We found strong evidence for Rensch’s rule (an increase in male-biased SSD with average body size) in A. carolinensis and evidence for the reversed version of Bergmann’s rule (an increase in body size towards warmer environments) in males. Across populations, SSD was positively related to temperature; however, female body size was not related to any climatic variable, suggesting that the latitudinal gradient of SSD might be driven by a gradient in the intensity of sexual selection acting on males. Sexual size dimorphism was positively correlated with sexual dimorphism in head shape and negatively correlated with limb length dimorphism, suggesting that sexual selection in males might drive the evolution of SSD and that differences in size and limb shape between sexes might represent alternative strategies to avoid competition for the same resources.

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.

A highly-complete chromosome-scale genome assembly of the brown anole, Anolis sagrei, provides insight into the evolution of sex chromosomes and is a crucial resource for this model lizard species.

Phenotypic Convergence Is Not Mirrored at the Protein Level in a Lizard Adaptive Radiation

There are many compelling examples of molecular convergence at individual genes. However, the prevalence and the relative importance of adaptive genome-wide convergence remain largely unknown. Many recent works have reported striking examples of excess genome-wide convergence, but some of these studies have been called into question because of the use of inappropriate null models. Here, we sequenced and compared the genomes of 12 species of anole lizards that have independently converged on suites of adaptive behavioral and morphological traits. Despite extensive searches for a genome-wide signature of molecular convergence, we found no evidence supporting molecular convergence at specific amino acids either at individual genes or at genome-wide comparisons; we also uncovered no evidence supporting an excess of adaptive convergence in the rates of amino acid substitutions within genes. Our findings indicate that comprehensive phenotypic convergence is not mirrored at genome-wide protein-coding levels in anoles, and therefore, that adaptive phenotypic convergence is likely not constrained by the evolution of many specific protein sequences or structures.

Markers for genetic change

Background and aims

Wildlife conservation has focused primarily on species for the last decades. Recently, popular perception and laws have begun to recognize the central importance of genetic diversity in the conservation of biodiversity. How to incorporate genetic diversity in ongoing monitoring and management of wildlife is still an open question.

Methods

We tested a panel of multiplexed, high-throughput sequenced introns in the small mammal communities of two UNESCO World Heritage Sites on different continents to assess their viability for large-scale monitoring of genetic variability in a spectrum of diverse species. To enhance applicability across other systems, the bioinformatic pipeline for primer design was outlined.

Results

The number of loci amplified and amplification evenness decreased as phylogenetic distance increased from the reference taxa, yet several loci were still variable across multiple mammal orders.

Conclusions

Genetic variability found is informative for population genetic analyses and for addressing phylogeographic and phylogenetic questions, illustrated by small mammal examples here.

Competing native and invasive Anolis lizards exhibit thermal preference plasticity in opposite directions

Invasive species have emerged as a significant problem in the age of anthropogenic change. Behavior can be key to invasive species success and is strongly affected by temperature. Therefore, knowledge of the temperature dependence of behavior is likely critical to understand invasive species dynamics and their interactions with native species. In this study, we tested for differences in thermal preference plasticity and temperature-dependent activity levels in a pair of congeneric lizards found in the United States: the invasive Anolis sagrei and the native A. carolinensis. We predicted that A. sagrei would demonstrate greater thermal preference plasticity and would utilize a higher and/or wider range of activity temperatures than A. carolinensis. Both would point to plasticity allowing A. sagrei to behaviorally exploiting thermal conditions that A. carolinensis cannot. We found that both species exhibited plasticity in thermal preference, but in opposite directions: preferred temperatures of A. carolinensis increased with acclimation temperature, while those of A. sagrei decreased. As a result, which species had a higher thermal preference changed with acclimation conditions. We saw no difference in overall field activity rates between the species, but that A. sagrei did tend to be active over a broader range of body temperatures. In sum, we found little evidence that differences in thermal preference plasticity between the species allow A. sagrei to remain active at a higher or broader temperature range than A. carolinensis. Nonetheless, the thermal preference data suggest complementary thermal preferences between the species that could promote microclimatic partitioning, though more work is required to test this idea.

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

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post-colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor-poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.

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 F_ST 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.

Recent Secondary Contacts, Linked Selection, and Variable Recombination Rates Shape Genomic Diversity in the Model Species Anolis carolinensis

Gaining a better understanding on how selection and neutral processes affect genomic diversity is essential to gain better insights into the mechanisms driving adaptation and speciation. However, the evolutionary processes affecting variation at a genomic scale have not been investigated in most vertebrate lineages. Here, we present the first population genomics survey using whole genome resequencing in the green anole (Anolis carolinensis). Anoles have been intensively studied to understand mechanisms underlying adaptation and speciation. The green anole in particular is an important model to study genome evolution. We quantified how demography, recombination, and selection have led to the current genetic diversity of the green anole by using whole-genome resequencing of five genetic clusters covering the entire species range. The differentiation of green anole's populations is consistent with a northward expansion from South Florida followed by genetic isolation and subsequent gene flow among adjacent genetic clusters. Dispersal out-of-Florida was accompanied by a drastic population bottleneck followed by a rapid population expansion. This event was accompanied by male-biased dispersal and/or selective sweeps on the X chromosome. We show that the interaction between linked selection and recombination is the main contributor to the genomic landscape of differentiation in the anole genome.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Local adaptation in mainland anole lizards: Integrating population history and genome-environment associations

Environmental gradients constrain physiological performance and thus species' ranges, suggesting that species occurrence in diverse environments may be associated with local adaptation. Genome-environment association analyses (GEAA) have become central for studies of local adaptation, yet they are sensitive to the spatial orientation of historical range expansions relative to landscape gradients. To test whether potentially adaptive genotypes occur in varied climates in wide-ranged species, we implemented GEAA on the basis of genomewide data from the anole lizards Anolis ortonii and Anolis punctatus, which expanded from Amazonia, presently dominated by warm and wet settings, into the cooler and less rainy Atlantic Forest. To examine whether local adaptation has been constrained by population structure and history, we estimated effective population sizes, divergence times, and gene flow under a coalescent framework. In both species, divergence between Amazonian and Atlantic Forest populations dates back to the mid-Pleistocene, with subsequent gene flow. We recovered eleven candidate genes involved with metabolism, immunity, development, and cell signaling in A. punctatus and found no loci whose frequency is associated with environmental gradients in A. ortonii. Distinct signatures of adaptation between these species are not associated with historical constraints or distinct climatic space occupancies. Similar patterns of spatial structure between selected and neutral SNPs along the climatic gradient, as supported by patterns of genetic clustering in A. punctatus, may have led to conservative GEAA performance. This study illustrates how tests of local adaptation can benefit from knowledge about species histories to support hypothesis formulation, sampling design, and landscape gradient characterization.

Barcoding blood meals: New vertebrate-specific primer sets for assigning taxonomic identities to host DNA from mosquito blood meals

The transmission dynamics of mosquito-vectored pathogens are, in part, mediated by mosquito host-feeding patterns. These patterns are elucidated using blood meal analysis, a collection of serological and molecular techniques that determine the taxonomic identities of the host animals from which blood meals are derived. Modern blood meal analyses rely on polymerase chain reaction (PCR), DNA sequencing, and bioinformatic comparisons of blood meal DNA sequences to reference databases. Ideally, primers used in blood meal analysis PCRs amplify templates from a taxonomically diverse range of vertebrates, produce a short amplicon, and avoid co-amplification of non-target templates. Few primer sets that fit these requirements are available for the cytochrome c oxidase subunit I (COI) gene, the species identification marker with the highest taxonomic coverage in reference databases. Here, we present new primer sets designed to amplify fragments of the DNA barcoding region of the vertebrate COI gene, while avoiding co-amplification of mosquito templates, without multiplexed or nested PCR. Primers were validated using host vertebrate DNA templates from mosquito blood meals of known origin, representing all terrestrial vertebrate classes, and field-collected mosquito blood meals of unknown origin. We found that the primers were generally effective in amplifying vertebrate host, but not mosquito DNA templates. Applied to the sample of unknown mosquito blood meals, > 98% (60/61) of blood meals samples were reliably identified, demonstrating the feasibility of identifying mosquito hosts with the new primers. These primers are beneficial in that they can be used to amplify COI templates from a diverse range of vertebrate hosts using standard PCR, thereby streamlining the process of identifying the hosts of mosquitoes, and could be applied to next generation DNA sequencing and metabarcoding approaches.

Inferring responses to climate dynamics from historical demography in neotropical forest lizards

We apply a comparative framework to test for concerted demographic changes in response to climate shifts in the neotropical lowland forests, learning from the past to inform projections of the future. Using reduced genomic (SNP) data from three lizard species codistributed in Amazonia and the Atlantic Forest (Anolis punctatus, Anolis ortonii, and Polychrus marmoratus), we first reconstruct former population history and test for assemblage-level responses to cycles of moisture transport recently implicated in changes of forest distribution during the Late Quaternary. We find support for population shifts within the time frame of inferred precipitation fluctuations (the last 250,000 y) but detect idiosyncratic responses across species and uniformity of within-species responses across forest regions. These results are incongruent with expectations of concerted population expansion in response to increased rainfall and fail to detect out-of-phase demographic syndromes (expansions vs. contractions) across forest regions. Using reduced genomic data to infer species-specific demographical parameters, we then model the plausible spatial distribution of genetic diversity in the Atlantic Forest into future climates (2080) under a medium carbon emission trajectory. The models forecast very distinct trajectories for the lizard species, reflecting unique estimated population densities and dispersal abilities. Ecological and demographic constraints seemingly lead to distinct and asynchronous responses to climatic regimes in the tropics, even among similarly distributed taxa. Incorporating such constraints is key to improve modeling of the distribution of biodiversity in the past and future.

Factors restricting the range expansion of the invasive green anole Anolis carolinensis on Okinawa Island, Japan

The green anole Anolis carolinensis invaded the Ogasawara Islands in Japan, drove various native species to extinction, and its distribution expanded 14 years after initial establishment. A. carolinensis invaded Okinawa Island, but it has not expanded its distribution in more than 25 years, although its density is extremely high in the southern region. To determine whether A. carolinensis has the potential to expand its distribution on Okinawa Island, we performed phylogenetic analysis of mitochondrial ND2 DNA sequences to study the origin of A. carolinensis that invaded Okinawa Island. We further used a species distribution model (MaxEnt) based on the distribution of native populations in North America to identify ecologically suitable areas on Okinawa Island. Nucleotide sequence analysis shows that the invader A. carolinensis originated in the western part of the Gulf Coast and inland areas of the United States and that a portion of the anoles on Okinawa was not introduced via the Ogasawara Islands. The MaxEnt predictions indicate that most areas in Okinawa Island are suitable for A. carolinensis. Therefore, A. carolinensis may have the potential to expand its distribution in Okinawa Island. The predictions indicate that habitat suitability is high in areas of high annual mean temperature and urbanized areas. The values of precipitation in summer in the northern region of Okinawa Island were higher compared with those of North America, which reduced the habitat suitability in Okinawa Island. Adaptation to low temperatures, an increase in the mean temperature through global warming, and an increase in open environments through land development will likely expand the distribution of A. carolinensis in Okinawa Island. Therefore, we must continue to monitor the introduced populations and be alert to the possibility that city planning that increases open environments may cause their range to expand.

LINE Insertion Polymorphisms are Abundant but at Low Frequencies across Populations of Anolis carolinensis

Vertebrate genomes differ considerably in size and structure. Among the features that show the most variation is the abundance of Long Interspersed Nuclear Elements (LINEs). Mammalian genomes contain 100,000s LINEs that belong to a single clade, L1, and in most species a single family is usually active at a time. In contrast, non-mammalian vertebrates (fish, amphibians and reptiles) contain multiple active families, belonging to several clades, but each of them is represented by a small number of recently inserted copies. It is unclear why vertebrate genomes harbor such drastic differences in LINE composition. To address this issue, we conducted whole genome resequencing to investigate the population genomics of LINEs across 13 genomes of the lizard Anolis carolinensis sampled from two geographically and genetically distinct populations in the Eastern Florida and the Gulf Atlantic regions of the United States. We used the Mobile Element Locator Tool to identify and genotype polymorphic insertions from five major clades of LINEs (CR1, L1, L2, RTE and R4) and the 41 subfamilies that constitute them. Across these groups we found large variation in the frequency of polymorphic insertions and the observed length distributions of these insertions, suggesting these groups vary in their activity and how frequently they successfully generate full-length, potentially active copies. Though we found an abundance of polymorphic insertions (over 45,000) most of these were observed at low frequencies and typically appeared as singletons. Site frequency spectra for most LINEs showed a significant shift toward low frequency alleles compared to the spectra observed for total genomic single nucleotide polymorphisms. Using Tajima's D, FST and the mean number of pairwise differences in LINE insertion polymorphisms, we found evidence that negative selection is acting on LINE families in a length-dependent manner, its effects being stronger in the larger Eastern Florida population. Our results suggest that a large effective population size and negative selection limit the expansion of polymorphic LINE insertions across these populations and that the probability of LINE polymorphisms reaching fixation is extremely low.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

Population Genetic Structure of the Dwarf Seahorse (Hippocampus zosterae) in Florida

The dwarf seahorse (Hippocampus zosterae) is widely distributed throughout near-shore habitats of the Gulf of Mexico and is of commercial significance in Florida, where it is harvested for the aquarium and curio trades. Despite its regional importance, the genetic structure of dwarf seahorse populations remains largely unknown. As an aid to ongoing conservation efforts, we employed three commonly applied mtDNA markers (ND4, DLoop and CO1) to investigate the genetic structuring of H. zosterae in Florida using samples collected throughout its range in the state. A total of 1450 bp provided sufficient resolution to delineate four populations of dwarf seahorses, as indicated by significant fixation indices. Despite an overall significant population structure, we observed evidence of interbreeding between individuals from geographically distant sites, supporting the hypothesis that rafting serves to maintain a degree of population connectivity. All individuals collected from Pensacola belong to a single distinct subpopulation, which is highly differentiated from the rest of Floridian dwarf seahorses sampled. Our findings highlight the utility of mtDNA markers in evaluating barriers to gene flow and identifying genetically distinct populations, which are vital to the development of comprehensive conservation strategies for exploited taxa.

Multilocus phylogenetic analyses of Hispaniolan and Bahamian trunk anoles (distichus species group)

The distichus species group includes six species and 21 subspecies of trunk ecomorph anoles distributed across Hispaniola and its satellite islands as well as the northern Bahamas. Although this group has long served as a model system for studies of reproductive character displacement, adaptation, behavior and speciation, it has never been the subject of a comprehensive phylogenetic analysis. Our goal here is to generate a multilocus phylogenetic dataset (one mitochondrial and seven nuclear loci) and to use this dataset to infer phylogenetic relationships among the majority of the taxa assigned to the distichus species group. We use these phylogenetic trees to address three topics about the group's evolution. First, we consider longstanding taxonomic controversies about the status of several species and subspecies assigned to the distichus species group. Second, we investigate the biogeographic history of the group and specifically test the hypotheses that historical division of Hispaniola into two paleo-islands contributed to the group's diversification and that Bahamian and Hispaniolan satellite island populations are derived from colonists from the main Hispaniolan landmass. Finally, third, we use comparative phylogenetic analyses to test the hypothesis that divergence between pale yellow and darkly pigmented orange or red dewlap coloration has occurred repeatedly across the distichus species group.

Lizards and LINEs: selection and demography affect the fate of L1 retrotransposons in the genome of the green anole (Anolis carolinensis)

Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila, applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard (Anolis carolinensis). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis. We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

Genetic variation in the green anole lizard (Anolis carolinensis) reveals island refugia and a fragmented Florida during the quaternary

The green anole lizard (Anolis carolinensis) is a model organism for behavior and genomics that is native to the southeastern United States. It is currently thought that the ancestors of modern green anoles dispersed to peninsular Florida from Cuba. However, the climatic changes and geological features responsible for the early diversification of A. carolinensis in North America have remained largely unexplored. This is because previous studies (1) differ in their estimates of the divergence times of populations, (2) are based on a single genetic locus or (3) did not test specific hypotheses regarding the geologic and topographic history of Florida. Here we provide a multi-locus study of green anole genetic diversity and find that the Florida peninsula contains a larger number of genetically distinct populations that are more diverse than those on the continental mainland. As a test of the island refugia hypothesis in Pleistocene Florida, we use a coalescent approach to estimate the divergence times of modern green anole lineages. We find that all demographic events occurred during or after the Upper Pliocene and suggest that green anole diversification was driven by population divergence on interglacial island refugia in Florida during the Lower Pleistocene, while the region was often separated from continental North America. When Florida reconnected to the mainland, two separate dispersal events led to the expansion of green anole populations across the Atlantic Seaboard and Gulf Coastal Plain.

Classification and evolutionary analysis of the basic helix-loop-helix gene family in the green anole lizard, Anolis carolinensis

Helix-loop-helix (bHLH) proteins play essential regulatory roles in a variety of biological processes. These highly conserved proteins form a large transcription factor superfamily, and are commonly identified in large numbers within animal, plant, and fungal genomes. The bHLH domain has been well studied in many animal species, but has not yet been characterized in non-avian reptiles. In this study, we identified 102 putative bHLH genes in the genome of the green anole lizard, Anolis carolinensis. Based on phylogenetic analysis, these genes were classified into 43 families, with 43, 24, 16, 3, 10, and 3 members assigned into groups A, B, C, D, E, and F, respectively, and 3 members categorized as "orphans". Within-group evolutionary relationships inferred from the phylogenetic analysis were consistent with highly conserved patterns observed for introns and additional domains. Results from phylogenetic analysis of the H/E(spl) family suggest that genome and tandem gene duplications have contributed to this family's expansion. Our classification and evolutionary analysis has provided insights into the evolutionary diversification of animal bHLH genes, and should aid future studies on bHLH protein regulation of key growth and developmental processes.