On the origin of appetite: GLWamide in jellyfish represents an ancestral satiety neuropeptide.

Food intake is regulated by internal state. This function is mediated by hormones and neuropeptides, which are best characterized in popular model species. However, the evolutionary origins of such feeding-regulating neuropeptides are poorly understood. We used the jellyfish Cladonema to address this question. Our combined transcriptomic, behavioral, and anatomical approaches identified GLWamide as a feeding-suppressing peptide that selectively inhibits tentacle contraction in this jellyfish. In the fruit fly Drosophila, myoinhibitory peptide (MIP) is a related satiety peptide. Surprisingly, we found that GLWamide and MIP were fully interchangeable in these evolutionarily distant species for feeding suppression. Our results suggest that the satiety signaling systems of diverse animals share an ancient origin.

Single-cell genomics unveiled a cryptic cyanobacterial lineage with a worldwide distribution hidden by a dinoflagellate host.

Cyanobacteria are one of the most important contributors to oceanic primary production and survive in a wide range of marine habitats. Much effort has been made to understand their ecological features, diversity, and evolution, based mainly on data from free-living cyanobacterial species. In addition, symbiosis has emerged as an important lifestyle of oceanic microbes and increasing knowledge of cyanobacteria in symbiotic relationships with unicellular eukaryotes suggests their significance in understanding the global oceanic ecosystem. However, detailed characteristics of these cyanobacteria remain poorly described. To gain better insight into marine cyanobacteria in symbiosis, we sequenced the genome of cyanobacteria collected from a cell of a pelagic dinoflagellate that is known to host cyanobacterial symbionts within a specialized chamber. Phylogenetic analyses using the genome sequence revealed that the cyanobacterium represents an underdescribed lineage within an extensively studied, ecologically important group of marine cyanobacteria. Metagenomic analyses demonstrated that this cyanobacterial lineage is globally distributed and strictly coexists with its host dinoflagellates, suggesting that the intimate symbiotic association allowed the cyanobacteria to escape from previous metagenomic studies. Furthermore, a comparative analysis of the protein repertoire with related species indicated that the lineage has independently undergone reductive genome evolution to a similar extent as Prochlorococcus, which has the most reduced genomes among free-living cyanobacteria. Discovery of this cyanobacterial lineage, hidden by its symbiotic lifestyle, provides crucial insights into the diversity, ecology, and evolution of marine cyanobacteria and suggests the existence of other undiscovered cryptic cyanobacterial lineages.

The effect of the doublesex gene in body colour masculinization of the damselfly Ischnura senegalensis.

Odonata species display a remarkable diversity of colour patterns, including intrasexual polymorphisms. In the damselfly (Ischnura senegalensis), the expression of a sex-determining transcription factor, the doublesex (Isdsx) gene is reportedly associated with female colour polymorphism (CP) (gynomorph for female-specific colour and andromorph for male-mimicking colour). Here, the function of Isdsx in thoracic coloration was investigated by electroporation-mediated RNA interference (RNAi). RNAi of the Isdsx common region in males and andromorphic females reduced melanization and thus changed the colour pattern into that of gynomorphic females, while the gynomorphic colour pattern was not affected. By contrast, RNAi against the Isdsx long isoform produced no changes, suggesting that the Isdsx short isoform is important for body colour masculinization in both males and andromorphic females. When examining the expression levels of five genes with differences between sexes and female morphs, two melanin-suppressing genes, black and ebony, were expressed at higher levels in the Isdsx RNAi body area than a control area. Therefore, the Isdsx short isoform may induce thoracic colour differentiation by suppressing black and ebony, thereby generating female CP in I. senegalensis. These findings contribute to the understanding of the molecular and evolutionary mechanisms underlying female CP in Odonata.

Identification of a reptile lyssavirus in Anolis allogus provided novel insights into lyssavirus evolution.

Lyssaviruses (genus Lyssavirus) are negative-strand RNA viruses belonging to the family Rhabdoviridae. Although a lyssa-like virus (frog lyssa-like virus 1 [FLLV-1]), which is distantly related to lyssaviruses, was recently identified in frogs, a large phylogenetic gap exists between those viruses, and thus the evolution of lyssaviruses is unclear. In this study, we detected a lyssa-like virus from publicly available RNA-seq data obtained using the brain and skin of Anolis allogus (Spanish flag anole), which was designated anole lyssa-like virus 1 (ALLV-1), and determined its complete coding sequence. Via mapping analysis, we demonstrated that ALLV-1 was actively replicating in the original brain and skin samples. Phylogenetic analyses revealed that ALLV-1 is more closely related to lyssaviruses than FLLV-1. Overall, the topology of the tree is compatible with that of hosts, suggesting the long-term co-divergence of lyssa-like and lyssaviruses and vertebrates. The ψ region, which is a long 3' untranslated region of unknown origin present in the G mRNA of lyssaviruses (approximately 400-700 nucleotides), is also present in the genome of ALLV-1, but it is much shorter (approximately 180 nucleotides) than those of lyssaviruses. Interestingly, FLLV-1 lacks the ψ region, suggesting that the ψ region was acquired after the divergence of the FLLV-1 and ALLV-1/lyssavirus lineages. To the best of our knowledge, this is the first report to identify a lyssa-like virus in reptiles, and thus, our findings provide novel insights into the evolution of lyssaviruses.

Genetic and plastic variation in opsin gene expression, light sensitivity, and female response to visual signals in the guppy.

According to the sensory drive model, variation in visual properties can lead to diverse female preferences, which in turn results in a range of male nuptial colors by way of sexual selection. However, the cause of variation in visual properties and the mechanism by which variation drives female response to visual signals remain unclear. Here, we demonstrate that both differences in the long-wavelength-sensitive 1 (LWS-1) opsin genotype and the light environment during rearing lead to variation in opsin gene expression. Opsin expression variation affects the visual sensitivity threshold to long wavelengths of light. Moreover, a behavioral assay using digitally modified video images showed that the expression of multiple opsin genes is positively correlated with the female responsiveness to images of males with luminous orange spots. The findings suggest that genetic polymorphisms and light environment in habitats induce variations in opsin gene expression levels. The variations may facilitate variations in visual sensitivity and female responsiveness to male body colors within and among populations.

Environmental factors explain spawning day deviation from full moon in the scleractinian coral Acropora.

Broadcast-spawning scleractinian corals annually release their gametes with high levels of synchrony, both within and among species. However, the timing of spawning can vary inter-annually. In particular, the night of spawning relative to the full moon phase can vary considerably among years at some locations. Although multiple environmental factors can affect the night of spawning, their effects have not been quantitatively assessed at the multi-regional level. In this study, we analysed environmental factors that are potentially correlated with spawning day deviation, in relation to the full moon phase, in Acropora corals inhabiting seven reefs in Australia and Japan. We accordingly found that sea surface temperature and wind speed within one to two months prior to the full moon of the spawning month were strongly correlated with spawning day deviations. In addition, solar flux had a weak effect on the night of spawning. These findings indicate that Acropora have the capacity to adjust their development and physiology in response to environmental factors for fine-tuning the timing of synchronous spawning, thereby maximizing reproductive success and post-fertilization survival.

Human-specific mutations in VMAT1 confer functional changes and multi-directional evolution in the regulation of monoamine circuits.

BACKGROUND: Neurochemicals like serotonin and dopamine play crucial roles in human cognitive and emotional functions. Vesicular monoamine transporter 1 (VMAT1) transports monoamine neurotransmitters, and its variant (136Thr) is associated with various psychopathological symptoms and reduced monoamine uptake relative to 136Ile. We previously showed that two human-specific amino acid substitutions (Glu130Gly and Asn136Thr/Ile) of VMAT1 were subject to positive natural selection. However, the potential functional alterations caused by these substitutions (Glu130Gly and Asn136Thr) remain unclear. To assess functional changes in VMAT1 from an evolutionary perspective, we reconstructed ancestral residues and examined the role of these substitutions in monoamine uptake in vitro using fluorescent false neurotransmitters (FFN), which are newly developed substances used to quantitatively assay VMATs. RESULTS: Immunoblotting confirmed that all the transfected YFP-VMAT1 variants are properly expressed in HEK293T cells at comparable levels, and no significant difference was seen in the density and the size of vesicles among them. Our fluorescent assays revealed a significant difference in FFN206 uptake among VMAT1 variants: 130Glu/136Asn, 130Glu/136Thr, and 130Gly/136Ile showed significantly higher levels of FFN206 uptake than 130Gly/136Asn and 130Gly/136Thr, indicating that both 130Glu and 136Ile led to increased neurotransmitter uptake, for which 136Thr and 136Asn were comparable by contrast. CONCLUSIONS: These findings suggest that monoamine uptake by VMAT1 initially declined (from 130Glu/136Asn to 130Gly/136Thr) in human evolution, possibly resulting in higher susceptibility to the external environment of our ancestors.

Invasive invertebrates associated with highly duplicated gene content.

Invasion of alien species has led to serious problems, including the destruction of native ecosystems. In general, invasive species adapt to new environments rapidly, suggesting that they have high genetic diversity that can directly influence environmental adaptability. However, it is not known how genomic architecture causes genetic diversity that leads to invasiveness. Recent studies have showed that the proportion of duplicated genes (PD ) in whole animal genomes correlate with environmental variability within a habitat. Here, we show that PD and propagule size significantly explain the differences in species categories (invasive species, noninvasive species, and parasites). PD correlated negatively with the propagule size. The residual values of regression of PD on propagule size revealed that the invasive species had higher PD values and larger propagule size than those of the noninvasive species, whereas the parasites had lower PD values and smaller propagule size than those of others. There were no correlations between the invasive species and other genomic factors including the genome size, number of genes, and certain gene families. Our results suggest that the PD values of a genome might be a potential genomic source causing genetic variations for adaptation to diverse environments. The results also showed that the invasiveness status of a species would be predicted by the residual values of regression of PD on propagule size. Our innovative approach provides a measure to estimate the environmental adaptability of organisms based on genomic data.

Candidate genes associated with color morphs of female-limited polymorphisms of the damselfly Ischnura senegalensis.

Many Odonata species exhibit female-limited polymorphisms, where one morph is similar to the conspecific male in body color and other traits (andromorph), whereas one or more other morphs differ from the male (gynomorphs). Here we investigated the differentially expressed transcripts (DETs) among males and two female morph groups (gynomorphs and andromorphs) using RNA-seq to identify candidate transcripts encoding female-limited polymorphisms in the damselfly Ischnura senegalensis. Seven DETs that had significantly different expression levels between males and gynomorphs, but not between males and andromorphs, were identified. The expression levels of four of these candidate genes, doublesex (dsx), black, ebony, and chaoptin (chp), were selected for further analysis using qRT-PCR. Sequence analysis of the dsx amplicons revealed that this gene produced at least three transcripts. Two short transcripts were mainly expressed in males and andromorphs, whereas the long transcript was specifically expressed in both morph female groups; that is, the expression pattern of the dsx splice variants in andromorphs was an intermediate between that of males and gynomorphs. Because the dsx gene functions as a transcription factor that regulates the sex-specific expression of multiple genes, its splice variants in I. senegalensis may explain why the andromorph is female but exhibits some masculinized traits. Because we did not detect different coding sequences of the candidate genes among the different morphs, a diallelic genomic region controlling alternative splicing of dsx, thus determining female-limited polymorphism in I. senegalensis most likely lies in a non-coding region of the dsx gene or in a gene upstream of it.

Balanced genetic diversity improves population fitness.

Although genetic diversity within a population is suggested to improve population-level fitness and productivity, the existence of these effects is controversial because empirical evidence for an ecological effect of genetic diversity and the underlying mechanisms is scarce and incomplete. Here, we show that the natural single-gene behavioural polymorphism (Rover and sitter) in Drosophila melanogaster has a positive effect on population fitness. Our simple numerical model predicted that the fitness of a polymorphic population would be higher than that expected with two monomorphic populations, but only under balancing selection. Moreover, this positive diversity effect of genetic polymorphism was attributable to a complementarity effect, rather than to a selection effect. Our empirical tests using the behavioural polymorphism in D. melanogaster clearly supported the model predictions. These results provide direct evidence for an ecological effect of genetic diversity on population fitness and its condition dependence.

Comparisons of behavioural and TRPA1 heat sensitivities in three sympatric Cuban Anolis lizards.

Thermal tolerances of organisms play a role in defining geographic ranges and occurrence of species. In Cuba, three sympatric species of Anolis lizards (Anolis allogus, Anolis homolechis and Anolis sagrei) inhabit different thermal microhabitats. A previous study found that these species showed distinct gene expression patterns in response to temperature stimuli, suggesting the genetically distinct thermal physiology among species. To investigate whether the Anolis species inhabiting locally distinct thermal habitats diverge their thermal tolerances, we first conducted behavioural experiments to analyse the temperatures at which the three Anolis species escape from heat source. Then, for each of the three species, we isolated cDNA encoding a putative molecular heat sensor, transient receptor potential ion channel ankyrin 1 (TRPA1), which has been suggested to play a role on eliciting behavioural responses to heat stimuli. We performed electrophysiological analysis to quantify activation temperature of Anolis TRPA1 to see whether the pattern of divergence in TRPA1 responses is congruent with that of divergence in behavioural responses. We found that temperatures triggering behavioural and TRPA1 responses were significantly lower for shade-dwelling species (A. allogus) than for sun-dwelling species (A. homolechis and A. sagrei). The ambient temperature of shade habitats where A. allogus occurs stays relatively cool compared to that of open habitats where A. homolechis and A. sagrei occur and bask. The high temperature thresholds of A. homolechis and A. sagrei may reflect their heat tolerances that would benefit these species to inhabit the open habitats.

Sex chromosome turnover contributes to genomic divergence between incipient stickleback species.

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration.

Natural selection in a population of Drosophila melanogaster explained by changes in gene expression caused by sequence variation in core promoter regions.

BACKGROUND: Understanding the evolutionary forces that influence variation in gene regulatory regions in natural populations is an important challenge for evolutionary biology because natural selection for such variations could promote adaptive phenotypic evolution. Recently, whole-genome sequence analyses have identified regulatory regions subject to natural selection. However, these studies could not identify the relationship between sequence variation in the detected regions and change in gene expression levels. We analyzed sequence variations in core promoter regions, which are critical regions for gene regulation in higher eukaryotes, in a natural population of Drosophila melanogaster, and identified core promoter sequence variations associated with differences in gene expression levels subjected to natural selection. RESULTS: Among the core promoter regions whose sequence variation could change transcription factor binding sites and explain differences in expression levels, three core promoter regions were detected as candidates associated with purifying selection or selective sweep and seven as candidates associated with balancing selection, excluding the possibility of linkage between these regions and core promoter regions. CHKov1, which confers resistance to the sigma virus and related insecticides, was identified as core promoter regions that has been subject to selective sweep, although it could not be denied that selection for variation in core promoter regions was due to linked single nucleotide polymorphisms in the regulatory region outside core promoter regions. Nucleotide changes in core promoter regions of CHKov1 caused the loss of two basal transcription factor binding sites and acquisition of one transcription factor binding site, resulting in decreased gene expression levels. Of nine core promoter regions regions associated with balancing selection, brat, and CG9044 are associated with neuromuscular junction development, and Nmda1 are associated with learning, behavioral plasticity, and memory. Diversity of neural and behavioral traits may have been maintained by balancing selection. CONCLUSIONS: Our results revealed the evolutionary process occurring by natural selection for differences in gene expression levels caused by sequence variation in core promoter regions in a natural population. The sequences of core promoter regions were diverse even within the population, possibly providing a source for natural selection.

Effects of light environment during growth on the expression of cone opsin genes and behavioral spectral sensitivities in guppies (Poecilia reticulata).

BACKGROUND: The visual system is important for animals for mate choice, food acquisition, and predator avoidance. Animals possessing a visual system can sense particular wavelengths of light emanating from objects and their surroundings and perceive their environments by processing information contained in these visual perceptions of light. Visual perception in individuals varies with the absorption spectra of visual pigments and the expression levels of opsin genes, which may be altered according to the light environments. However, which light environments and the mechanism by which they change opsin expression profiles and whether these changes in opsin gene expression can affect light sensitivities are largely unknown. This study determined whether the light environment during growth induced plastic changes in opsin gene expression and behavioral sensitivity to particular wavelengths of light in guppies (Poecilia reticulata). RESULTS: Individuals grown under orange light exhibited a higher expression of long wavelength-sensitive (LWS) opsin genes and a higher sensitivity to 600-nm light than those grown under green light. In addition, we confirmed that variations in the expression levels of LWS opsin genes were related to the behavioral sensitivities to long wavelengths of light. CONCLUSIONS: The light environment during the growth stage alters the expression levels of LWS opsin genes and behavioral sensitivities to long wavelengths of light in guppies. The plastically enhanced sensitivity to background light due to changes in opsin gene expression can enhance the detection and visibility of predators and foods, thereby affecting survival. Moreover, changes in sensitivities to orange light may lead to changes in the discrimination of orange/red colors of male guppies and might alter female preferences for male color patterns.

Differentially expressed genes associated with adaptation to different thermal environments in three sympatric Cuban Anolis lizards.

How animals achieve evolutionary adaptation to different thermal environments is an important issue for evolutionary biology as well as for biodiversity conservation in the context of recent global warming. In Cuba, three sympatric species of Anolis lizards (Anolis allogus, A. homolechis and A. sagrei) inhabit different thermal microhabitats, thereby providing an excellent opportunity to examine how they have adapted to different environmental temperatures. Here, we performed RNA-seq on the brain, liver and skin tissues from these three species to analyse their transcriptional responses at two different temperatures. In total, we identified 400, 816 and 781 differentially expressed genes (DEGs) between the two temperatures in A. allogus, A. homolechis and A. sagrei, respectively. Only 62 of these DEGs were shared across the three species, indicating that global transcriptional responses have diverged among these species. Gene ontology (GO) analysis showed that large numbers of ribosomal protein genes were DEGs in the warm-adapted A. homolechis, suggesting that the upregulation of protein synthesis is an important physiological mechanism in the adaptation of this species to hotter environments. GO analysis also showed that GO terms associated with circadian regulation were enriched in all three species. A gene associated with circadian regulation, Nr1d1, was detected as a DEG with opposite expression patterns between the cool-adapted A. allogus and the hot-adapted A. sagrei. Because the environmental temperature fluctuates more widely in open habitats than in forests throughout the day, the circadian thermoregulation could also be important for adaptation to distinct thermal habitats.

Lack of genetic variation prevents adaptation at the geographic range margin in a damselfly.

What limits a species' distribution in the absence of physical barriers? Genetic load due to asymmetric gene flow and the absence of genetic variation due to lack of gene flow are hypothesized to constrain adaptation to novel environments in marginal populations, preventing range expansion. Here, we examined the genetic structure and geographic variation in morphological traits in two damselflies (Ischnura asiatica and I. senegalensis) along a latitudinal gradient in Japan, which is the distribution centre of I. asiatica and the northern limit of I. senegalensis. Genomewide genetic analyses found a loss of genetic diversity at the edge of distribution in I. senegalensis but consistently high diversity in I. asiatica. Gene flow was asymmetric in a south-north direction in both species. Although body size and wing loading showed decreasing latitudinal clines (smaller in north) in I. asiatica in Japan, increasing latitudinal clines (larger in north) in these phenotypic markers were observed in I. senegalensis, particularly near the northern boundary, which coincided well with the location where genetic diversity began a sharp decline. In ectothermic animals, increasing latitudinal cline in these traits was suggested to be established when they failed to adapt to thermal gradient. Therefore, our findings support the possibility that a lack of genetic variation rather than geneflow swamping is responsible for the constraint of adaptation at the margin of geographic distribution.

Contribution of nonohnologous duplicated genes to high habitat variability in mammals.

The mechanism by which genetic systems affect environmental adaptation is a focus of considerable attention in the fields of ecology, evolution, and conservation. However, the genomic characteristics that constrain adaptive evolution have remained unknown. A recent study showed that the proportion of duplicated genes in whole Drosophila genomes correlated with environmental variability within habitat, but it remains unclear whether the correlation is observed even in vertebrates whose genomes including a large number of duplicated genes generated by whole-genome duplication (WGD). Here, we focus on fully sequenced mammalian genomes that experienced WGD in early vertebrate lineages and show that the proportion of small-scale duplication (SSD) genes in the genome, but not that of WGD genes, is significantly correlated with habitat variability. Moreover, species with low habitat variability have a higher proportion of lost duplicated genes, particularly SSD genes, than those with high habitat variability. These results indicate that species that inhabit variable environments may maintain more SSD genes in their genomes and suggest that SSD genes are important for adapting to novel environments and surviving environmental changes. These insights may be applied to predicting invasive and endangered species.

Developmental stages for the divergence of relative limb length between a twig and a trunk-ground Anolis lizard species.

The divergent evolution of niche-related traits can facilitate adaptive radiation, yet identification of the genetic or molecular mechanisms underlying such trait changes remains a major challenge in evolutionary biology. Conducting a detailed morphological comparison along growth trajectories is a powerful method for observing the formation of differences in niche-related traits. Here, we focused on hindlimb length of Anolis lizards, differences in which are related to adaptation for use of different microhabitats. We measured the length of hindlimb skeletons in different ecomorphs of anole lizards (A. sagrei, a trunk-ground ecomorph with long hindlimbs, and A. angusticeps, a twig ecomorph with short hindlimbs) from early embryonic stages to adulthood, to determine which hindlimb elements mainly differentiate the species and the timing of the formation of these differences. With respect to the digit, differences between the species mainly occurred during the embryonic stages of interdigit reduction, when the cartilage of the distal phalanges was simultaneously forming. In addition, we compared the relative length of developing autopods in early embryonic stages using whole-mount in situ hybridization before the formation of the cartilaginous bones, and the results showed that the relative growth rate of the Hoxa11-negative distal region in A. sagrei was greater than that in A. angusticeps. Our results show that there are several important developmental stages for hindlimb length differentiation between A. angusticeps and A. sagrei, depending on which hindlimb element is considered. In particular, the species differences were largely due to variations in digit length, which arose at early embryonic stages.

Flower color polymorphism maintained by overdominant selection in Sisyrinchium sp.

Negative frequency-dependent selection derived from positive frequency-dependent foraging is the best-known selection force maintaining genetic polymorphism within a population. However, in flowering plants, positive frequency-dependent foraging by pollinators is expected to accelerate the loss of low-frequency morphs by conferring a fitness advantage to the common morph, leading to monomorphism. In Japan, a non-native species, Sisyrinchium sp., exhibits conspicuous flower color polymorphism within a population comprising both purple morphs (homozygous recessive) and white morphs (heterozygous or homozygous dominant). Here we quantified genotype-specific reproductive success in order to reveal the contribution of overdominant selection on the maintenance of flower color polymorphism in this species. In artificial pollination experiments using individuals with identified genotypes, female reproductive success was higher in the heterozygote than in either homozygote. The frequency of purple morphs in natural populations (ca. 31%) is similar to the frequency predicted by overdominant selection (25%). Our results suggest that overdominant selection contributes to the maintenance of color morphs in the natural population of this species.

Functional morphology and comparative anatomy of appendicular musculature in Cuban Anolis lizards with different locomotor habits.

We examined the diversity of the musculoskeletal morphology in the limbs of Anolis lizards with different habitats and identified variations in functional and morphological adaptations to different ecologies or behaviors. Dissection and isolation of 40 muscles from the fore- and hindlimbs of five species of Anolis were performed, and the muscle mass and length of the moment arm were compared after body size effects were removed. Ecologically and behaviorally characteristic morphological differences were observed in several muscles. Well-developed hindlimb extensors were observed in ground-dwelling species, A. sagrei and A. bremeri, and were considered advantageous for running, whereas adept climber species possessed expanded femoral retractors for weight-bearing during climbing. Moreover, morphological variations were observed among arboreal species. Wider excursions of the forelimb joint characterized A. porcatus, presumably enabling branch-to-branch locomotion, while A. equestris and A. angusticeps possessed highly developed adductor muscles for grasping thick branches or twigs. These findings suggest divergent evolution of musculoskeletal characteristic in the limbs within the genus Anolis, with correlations observed among morphological traits, locomotor performance, and habitat uses.