Novel herpesvirus in the critically endangered Galapagos pink land iguana

Virus surveillance in wildlife is important to understanding ecosystem health, taxonomy, and evolution. Nevertheless, viruses in reptiles, and specifically in squamates, continue to be understudied. Herein, we conducted a health assessment on the critically endangered Galapagos pink land iguana (Conolophus marthae) and the vulnerable Galapagos land iguana (Conolophus subcristatus). We collected oral and/or cloacal swabs from 47 clinically healthy iguanas and tested for adenovirus (cloacal swabs, n = 47) and herpesvirus (oral swabs, n = 45) using broad-spectrum PCRs. Two out of 38 (5.3 %) Galapagos pink land iguanas tested positive for herpesvirus, while no herpesvirus was detected in all Galapagos land iguanas (n = 7). Both herpesviral sequences were identical between them and divergent (61.9 % amino acid identity) when compared to the closest herpesvirus sequences available in GenBank/EMBL/DDBJ. The genetic distance between this and other herpesviruses is consistent with its classification as a novel virus species. All iguanas were negative for adenovirus. This is the first description of a herpesvirus in iguanas of the Galapagos islands, and the first report of a potential pathogen for the iconic Galapagos pink land iguana. Further research is needed to understand the implications of this virus in the conservation and management of one of the most endangered iguana species in the world.

Dispersal, not vicariance, explains the biogeographic origin of iguanas on Madagascar

Lizards of the clade Iguanidae (sensu lato) are primarily a New World group. Thus, the remarkable presence of an endemic lineage of iguanas (family Opluridae) on the isolated Indian Ocean island of Madagascar has long been considered a biogeographic anomaly. Previous work attributed this disjunct extant distribution to: (1) vicariance at about 140-165 Ma, caused by the breakup of Gondwana and the separation of South America, Africa, and Madagascar (with subsequent extinction of iguanas on Africa, and potentially other Gondwanan landmasses), (2) vicariance at about 80-90 Ma, caused by the sundering of hypothesized land-bridge connections between South America, Antarctica, India, and Madagascar, or (3) long-distance overwater dispersal from South America to Madagascar. Each hypothesis has been supported with molecular divergence dating analyses, and thus the biogeographic origin of the Opluridae is not yet well resolved. Here we utilize genetic sequences of ultraconserved elements for all Iguania families and the majority of Iguanidae (s.l.) genera, and morphological data for extant and fossil taxa (used for divergence dating analyses), to produce the most comprehensive dataset applied to date to test these origin hypotheses. We find strong support for a sister relationship between the Opluridae (Madagascar) and Leiosauridae (South America). Divergence of the Opluridae from Leiosauridae is dated to between the late Cretaceous and mid-Paleogene, at a time when Madagascar was already an island and was isolated from all other Gondwanan landmasses. Consequently, our results support a hypothesis of long-distance overwater dispersal of the Opluridae lineage, either directly from South America to Madagascar or potentially via Antarctica or Africa, leading to this radiation of iguanas in the Indian Ocean.

Defining intraspecific conservation units in the endemic Cuban Rock Iguanas (Cyclura nubila nubila)

Defining conservation units is an important step in species management and requires interpretation of the genetic diversity and ecological function of the taxon being considered. We used the endemic Cuban Rock Iguanas (Cyclura nubila nubila) as a model to highlight this challenge and examined patterns of its intraspecific genetic diversity across Cuba. We evaluated nuclear (microsatellite loci) and mitochondrial diversity across eight populations from the island and its off-shore cays, and applied the population genetics results for assignment of Management Unit (MU) status and Evolutionary Significant Units (ESUs) based on phylogeographic and time of divergence information. We identified at least six distinct Cuban Rock Iguana MUs, encompassing demographically isolated and genetically differentiated populations across Cuba, most with low effective population size, declining populations, and with high risk of inbreeding and genetic drift. Hence, each MU should be considered of urgent conservation priority. Given the key ecological seed dispersal role of C. n. nubila, the disappearance of any MU could trigger the loss of local ecological functional diversity and major negative impacts on their ecosystems. Two divergent ESUs were also identified, exhibiting an historical east-west geographic separation on Cuba. Based on a Caribbean phylogeographic assessment, our findings strengthen the conclusion that all geographically and evolutionarily differentiated Cyclura species and subspecies across the archipelago warrant ESU distinction.

Characterization and distribution of GHRH, PACAP, TRH, SST and IGF1 mRNAs in the green iguana

The somatotropic axis (SA) regulates numerous aspects of vertebrate physiology such as development, growth, and metabolism and has influence on several tissues including neural, immune, reproductive and gastric tract. Growth hormone (GH) is a key component of SA, it is synthesized and released mainly by pituitary somatotrophs, although now it is known that virtually all tissues can express GH, which, in addition to its well-described endocrine roles, also has autocrine/paracrine/intracrine actions. In the pituitary, GH expression is regulated by several hypothalamic neuropeptides including GHRH, PACAP, TRH and SST. GH, in turn, regulates IGF1 synthesis in several target tissues, adding complexity to the system since GH effects can be exerted either directly or mediated by IGF1. In reptiles, little is known about the SA components and their functional interactions. The aim of this work was to characterize the mRNAs of the principal SA components in the green iguana and to develop the tools that allow the study of the structural and functional evolution of this system in reptiles. By employing RT-PCR and RACE, the cDNAs encoding for GHRH, PACAP, TRH, SST and IGF1 were amplified and sequenced. Results showed that these cDNAs coded for the corresponding protein precursors of 154, 170, 243, 113, and 131 amino acids, respectively. Of these, GHRH, PACAP, SST and IGF1 precursors exhibited a high structural conservation with respect to its counterparts in other vertebrates. On the other hand, iguana's TRH precursor showed 7 functional copies of mature TRH (pyr-QHP-NH₂), as compared to 4 and 6 copies of TRH in avian and mammalian proTRH sequences, respectively. It was found that in addition to its primary production site (brain for GHRH, PACAP, TRH and SST, and liver for IGF1), they were also expressed in other peripheral tissues, i.e. testes and ovaries expressed all the studied mRNAs, whereas TRH and IGF1 mRNAs were observed ubiquitously in all tissues considered. These results show that the main SA components in reptiles of the Squamata Order maintain a good structural conservation among vertebrate phylogeny, and suggest important physiological interactions (endocrine, autocrine and/or paracrine) between them due to their wide peripheral tissue expression.

Predictability in evolution: Adaptation of the Bonaire anole (Anolis bonairensis) to an extreme environment

The extent to which evolution is deterministic (predictable), or random, is a fundamental question in evolution. This case study attempts to determine the extent to which interspecific divergence can be predicted from intraspecific trends related species. The mountainous Lesser Antilles are occupied by one or two anole species with very substantial intraspecific differences in the quantitative traits between xeric and rainforest habitats. These ecologically determined differences tend to be in parallel in each island species. A related species (Anolis bonairensis) lives on the far more xeric island of Bonaire, and this study tests the extent to which its interspecific divergence in hue and pattern traits can be predicted from the parallel intraspecific variation exhibited in Lesser Antillean anoles. Regression against a multivariate climate variable suggests that the hue and pattern of the Bonaire anole are consistently predicted from the ecologically determined intraspecific variation of its Lesser Antillean relatives. However, this predictability may be less consistent with other character systems, for example, scalation.

Differential responses of the somatotropic and thyroid axes to environmental temperature changes in the green iguana

Growth hormone (GH), together with thyroid hormones (TH), regulates growth and development, and has critical effects on vertebrate metabolism. In ectotherms, these physiological processes are strongly influenced by environmental temperature. In reptiles, however, little is known about the direct influences of this factor on the somatotropic and thyroid axes. Therefore, the aim of this study was to describe the effects of both acute (48h) and chronic (2weeks) exposure to sub-optimal temperatures (25 and 18°C) upon somatotropic and thyroid axis function of the green iguana, in comparison to the control temperature (30-35°C). We found a significant increase in GH release (2.0-fold at 25°C and 1.9-fold at 18°C) and GH mRNA expression (up to 3.7-fold), mainly under chronic exposure conditions. The serum concentration of insulin-like growth factor-I (IGF-I) was significantly greater after chronic exposure (18.5±2.3 at 25°C; 15.92±3.4 at 18°C; vs. 9.3±1.21ng/ml at 35°C), while hepatic IGF-I mRNA expression increased up to 6.8-fold. Somatotropic axis may be regulated, under acute conditions, by thyrotropin-releasing hormone (TRH) that significantly increased its hypothalamic concentration (1.45 times) and mRNA expression (0.9-fold above control), respectively; and somatostatin (mRNA expression increased 1.0-1.2 times above control); and under chronic treatment, by pituitary adenylate cyclase-activating peptide (PACAP mRNA expression was increased from 0.4 to 0.6 times). Also, it was shown that, under control conditions, injection of TRH stimulated a significant increase in circulating GH. On the other hand, while there was a significant rise in the hypothalamic content of TRH and its mRNA expression, this hormone did not appear to influence the thyroid axis activity, which showed a severe diminution in all conditions of cold exposure, as indicated by the decreases in thyrotropin (TSH) mRNA expression (up to one-eight of the control), serum T4 (from 11.6±1.09 to 5.3±0.58ng/ml, after 2weeks at 18°C) and T3 (from 0.87±0.09 to 0.05±0.01ng/ml, under chronic conditions at 25°C), and Type-2 deiodinase (D2) activity (from 992.5±224 to 213.6±26.4fmolI(125)T4/mgh). The reduction in thyroid activity correlates with the down-regulation of metabolism as suggested by the decrease in the serum glucose and free fatty acid levels. These changes apparently were independent of a possible stress response, at least under acute exposure to both temperatures and in chronic treatment to 25°C, since serum corticosterone had no significant changes in these conditions, while at chronic 18°C exposure, a slight increase (0.38 times above control) was found. Thus, these data suggest that the reptilian somatotropic and thyroid axes have differential responses to cold exposure, and that GH and TRH may play important roles associated to adaptation mechanisms that support temperature acclimation in the green iguana.

Genetic Evidence of Hybridization between the Endangered Native Species Iguana delicatissima and the Invasive Iguana iguana (Reptilia, Iguanidae) in the Lesser Antilles: Management Implications

The worldwide increase of hybridization in different groups is thought to have become more important with the loss of isolating barriers and the introduction of invasive species. This phenomenon could result in the extinction of endemic species. This study aims at investigating the hybridization dynamics between the endemic and threatened Lesser Antillean iguana (Iguana delicatissima) and the invasive common green iguana (Iguana iguana) in the Lesser Antilles, as well as assessing the impact of interspecific hybridization on the decline of I. delicatissima. 59 I. delicatissima (5 localities), 47 I. iguana (12 localities) and 27 hybrids (5 localities), who were all identified based on morphological characters, have been genotyped at 15 microsatellites markers. We also sequenced hybrids using ND4 mitochondrial loci to further investigate mitochondrial introgression. The genetic clustering of species and hybrid genetic assignment were performed using a comparative approach, through the implementation of a Discriminant Analysis of Principal Component (DAPC) based on statistics, as well as genetic clustering approaches based on the genetic models of several populations (Structure, NewHybrids and HIest), in order to get full characterization of hybridization patterns and introgression dynamics across the islands. The iguanas identified as hybrids in the wild, thanks to morphological analysis, were all genetically F1, F2, or backcrosses. A high proportion of individuals were also the result of a longer-term admixture. The absence of reproductive barriers between species leads to hybridization when species are in contact. Yet morphological and behavioral differences between species could explain why males I. iguana may dominate I. delicatissima, thus resulting in short-term species displacement and extinction by hybridization and recurrent introgression from I. iguana toward I. delicatissima. As a consequence, I. delicatissima gets eliminated through introgression, as observed in recent population history over several islands. These results have profound implications for species management of the endangered I. delicatissima and practical conservation recommendations are being discussed in the light of these findings.

Inferred vs realized patterns of gene flow: an analysis of population structure in the Andros Island Rock Iguana

Ecological data, the primary source of information on patterns and rates of migration, can be integrated with genetic data to more accurately describe the realized connectivity between geographically isolated demes. In this paper we implement this approach and discuss its implications for managing populations of the endangered Andros Island Rock Iguana, Cyclura cychlura cychlura. This iguana is endemic to Andros, a highly fragmented landmass of large islands and smaller cays. Field observations suggest that geographically isolated demes were panmictic due to high, inferred rates of gene flow. We expand on these observations using 16 polymorphic microsatellites to investigate the genetic structure and rates of gene flow from 188 Andros Iguanas collected across 23 island sites. Bayesian clustering of specimens assigned individuals to three distinct genotypic clusters. An analysis of molecular variance (AMOVA) indicates that allele frequency differences are responsible for a significant portion of the genetic variance across the three defined clusters (F_st =  0.117, p0.01). These clusters are associated with larger islands and satellite cays isolated by broad water channels with strong currents. These findings imply that broad water channels present greater obstacles to gene flow than was inferred from field observation alone. Additionally, rates of gene flow were indirectly estimated using BAYESASS 3.0. The proportion of individuals originating from within each identified cluster varied from 94.5 to 98.7%, providing further support for local isolation. Our assessment reveals a major disparity between inferred and realized gene flow. We discuss our results in a conservation perspective for species inhabiting highly fragmented landscapes.

Characterization of pituitary growth hormone and its receptor in the green iguana (Iguana iguana)

Pituitary growth hormone (GH) has been studied in most vertebrate groups; however, only a few studies have been carried out in reptiles. Little is known about pituitary hormones in the order Squamata, to which the green iguana (gi) belongs. In this work, we characterized the hypophysis of Iguana iguana morphologically. The somatotrophs (round cells of 7.6-10 μm containing 250- to 300-nm secretory granules where the giGH is stored) were found, by immunohistochemistry and in situ hybridization, exclusively in the caudal lobe of the pars distalis, whereas the lactotrophs were distributed only in the rostral lobe. A pituitary giGH-like protein was obtained by immuno-affinity chromatography employing a heterologous antibody against chicken GH. giGH showed molecular heterogeneity (22, 44, and 88 kDa by SDS-PAGE/Western blot under non-reducing conditions and at least four charge variants (pIs 6.2, 6.5, 6.9, 7.4) by isoelectric focusing. The pituitary giGH cDNA (1016 bp), amplified by PCR and RACE, encodes a pre-hormone of 218 aa, of which 190 aa correspond to the mature protein and 28 aa to the signal peptide. The giGH receptor cDNA was also partially sequenced. Phylogenetic analyses of the amino acid sequences of giGH and giGHR homologs in vertebrates suggest a parallel evolution and functional relationship between the GH and its receptor.

Identification of a novel splicing form of amelogenin gene in a reptile, Ctenosaura similis

Amelogenin, the major enamel matrix protein in tooth development, has been demonstrated to play a significant role in tooth enamel formation. Previous studies have identified the alternative splicing of amelogenin in many mammalian vertebrates as one mechanism for amelogenin heterogeneous expression in teeth. While amelogenin and its splicing forms in mammalian vertebrates have been cloned and sequenced, the amelogenin gene, especially its splicing forms in non-mammalian species, remains largely unknown. To better understand the mechanism underlying amelogenin evolution, we previously cloned and characterized an amelogenin gene sequence from a squamate, the green iguana. In this study, we employed RT-PCR to amplify the amelogenin gene from the black spiny-tailed iguana Ctenosaura similis teeth, and discovered a novel splicing form of the amelogenin gene. The transcript of the newly identified iguana amelogenin gene (named C. Similis-T2L) is 873 nucleotides long encoding an expected polypeptide of 206 amino acids. The C. Similis-T2L contains a unique exon denominated exon X, which is located between exon 5 and exon 6. The C. Similis-T2L contains 7 exons including exon 1, 2, 3, 5, X, 6, and 7. Analysis of the secondary and tertiary structures of T2L amelogenin protein demonstrated that exon X has a dramatic effect on the amelogenin structures. This is the first report to provide definitive evidence for the amelogenin alternative splicing in non-mammalian vertebrates, revealing a unique exon X and the splicing form of the amelogenin gene transcript in Ctenosaura similis.

Genotoxic effects in blood cells of Mus musculus and Iguana iguana living near coal mining areas in Colombia

Coal is a mixture of chemicals with the capacity of promoting biochemical changes that may lead to DNA damage. In this study, the comet assay in peripheral blood cells, and the micronucleus test in blood smears were used to evaluate potential genotoxic effects derived from exposure to coal mining activities on wild populations of Mus musculus and Iguana iguana. Four locations from Colombia were evaluated: La Loma and La Jagua de Ibirico, two municipalities located near coal mining fields at the Department of Cesar; and Valledupar and Arjona, cities used as reference sites, both localized at least 100 and 200km far from the mines, respectively. Compared to Valledupar and Arjona, animals collected in close proximity to coal mining areas showed highest percentages of DNA damage for both species, evidencing that living around coal mining fields may result in an increase of DNA lesions in blood cells of rodents and reptiles. The results for micronucleus test were conflicting. Mice from Arjona had greater number of cells with micronucleus than those from the other studied locations, probably as a result of infection found by blood parasites. In summary, it was demonstrated that animals living around coal mining areas have a greater chance of having DNA damage, as measured by the comet assay, than those from sites far from the coal dust source.

Multi-locus estimates of population structure and migration in a fence lizard hybrid zone

A hybrid zone between two species of lizards in the genus Sceloporus (S. cowlesi and S. tristichus) on the Mogollon Rim in Arizona provides a unique opportunity to study the processes of lineage divergence and merging. This hybrid zone involves complex interactions between 2 morphologically and ecologically divergent subspecies, 3 chromosomal groups, and 4 mitochondrial DNA (mtDNA) clades. The spatial patterns of divergence between morphology, chromosomes and mtDNA are discordant, and determining which of these character types (if any) reflects the underlying population-level lineages that are of interest has remained impeded by character conflict. The focus of this study is to estimate the number of populations interacting in the hybrid zone using multi-locus nuclear data, and to then estimate the migration rates and divergence time between the inferred populations. Multi-locus estimates of population structure and gene flow were obtained from 12 anonymous nuclear loci sequenced for 93 specimens of Sceloporus. Population structure estimates support two populations, and this result is robust to changes to the prior probability distribution used in the Bayesian analysis and the use of spatially-explicit or non-spatial models. A coalescent analysis of population divergence suggests that gene flow is high between the two populations, and that the timing of divergence is restricted to the Pleistocene. The hybrid zone is more accurately described as involving two populations belonging to S. tristichus, and the presence of S. cowlesi mtDNA haplotypes in the hybrid zone is an anomaly resulting from mitochondrial introgression.

Phylogenetic evidence for lateral gene transfer in the intestine of marine iguanas

Background

Lateral gene transfer (LGT) appears to promote genotypic and phenotypic variation in microbial communities in a range of environments, including the mammalian intestine. However, the extent and mechanisms of LGT in intestinal microbial communities of non-mammalian hosts remains poorly understood.

Methodology/Principal Findings

We sequenced two fosmid inserts obtained from a genomic DNA library derived from an agar-degrading enrichment culture of marine iguana fecal material. The inserts harbored 16S rRNA genes that place the organism from which they originated within Clostridium cluster IV, a well documented group that habitats the mammalian intestinal tract. However, sequence analysis indicates that 52% of the protein-coding genes on the fosmids have top BLASTX hits to bacterial species that are not members of Clostridium cluster IV, and phylogenetic analysis suggests that at least 10 of 44 coding genes on the fosmids may have been transferred from Clostridium cluster XIVa to cluster IV. The fosmids encoded four transposase-encoding genes and an integrase-encoding gene, suggesting their involvement in LGT. In addition, several coding genes likely involved in sugar transport were probably acquired through LGT.

Conclusion

Our phylogenetic evidence suggests that LGT may be common among phylogenetically distinct members of the phylum Firmicutes inhabiting the intestinal tract of marine iguanas.

Characterization and evolution of MHC class II B genes in Galápagos marine iguanas (Amblyrhynchus cristatus)

Major histocompatibility complex (MHC) class II molecules play a key role in the adaptive immune system of vertebrates. Class II B genes appear to evolve in a very different manner in mammals and birds. Orthology is commonly observed among mammal loci, while genes tend to cluster phylogenetically within bird species. Here we present class II B data from a representative of another major group of amniotes, the squamates (i.e. lizards, snakes, amphisbaenians), with the ultimate goal of placing mammalian and avian MHC evolution into a broader context. In this study, eight class II B cDNA sequences were obtained from the Galápagos marine iguana (Amblyrhynchus cristatus) which were divided into five locus groups, Amcr-DAB1 through -DAB5, based on similarities along most of the coding and noncoding portions of the transcribed gene. All marine iguana sequences were monophyletic with respect to class II genes from other vertebrates indicating that they originated from a common ancestral locus after squamates split from other reptiles. The beta-1 domain, which is involved in antigen binding, exhibited signatures of positive selection as well as interlocus gene conversion in both long and short tracts-a pattern also observed in birds and fish, but not in mammals. On the other hand, the beta-2 domain was divergent between gene groups, which is characteristic of mammals. Based on these results, we preliminarily show that squamate class II B genes have been shaped by a unique blend of evolutionary forces that have been observed in differing degrees in other vertebrates.

[Specific organization of ribosomal DNA arrays in Squamata]

A first report on structural organization of ribosomal DNA arrays in some members of the order Squamata is presented. The data obtained point to unusually small (for vertebrates) size of the rDNA repetitive unit (approximately, 10 to 15 kb) in the lizard species examined. Analysis of BAC library of Uta stansburiana (Iguania) showed that haploid genome of this lizard contained a single cluster, consisting of about ten rDNA repeats. Determination of the extent of rDNA unit repetition in some other representatives of the order Squamata, using the method of comparative real-time PCR, showed that the number of rDNA units varied from one or several dozens in Iguanina to several hundred repeats in Scincomorpha and Varonoidea. The results are discussed in terms of an ambiguous position of the family Iguania on the evolutionary trees constructed based on morphological and molecular data.

Molecular and morphological analysis of the critically endangered Fijian iguanas reveals cryptic diversity and a complex biogeographic history

The Pacific iguanas of the Fijian and Tongan archipelagos are a biogeographic enigma in that their closest relatives are found only in the New World. They currently comprise two genera and four species of extinct and extant taxa. The two extant species, Brachylophus fasciatus from Fiji, Tonga, and Vanuatu and Brachylophus vitiensis from western Fiji, are of considerable conservation concern with B. vitiensis listed as critically endangered. A recent molecular study has shown that Brachylophus comprised three evolutionarily significant units. To test these conclusions and to reevaluate the phylogenetic and biogeographic relationships within Brachylophus, we generated an mtDNA dataset consisting of 1462 base pairs for 61 individuals from 13 islands, representing both Brachylophus species. Unweighted parsimony analyses and Bayesian analyses produced a well-resolved phylogenetic hypothesis supported by high bootstrap values and posterior probabilities within Brachylophus. Our data reject the monophyly of specimens previously believed to comprise B. fasciatus. Instead, our data demonstrate that living Brachylophus comprise three robust and well-supported clades that do not correspond to current taxonomy. One of these clades comprises B. fasciatus from the Lau group of Fiji and Tonga (type locality for B. fasciatus), while a second comprises putative B. fasciatus from the central regions of Fiji, which we refer to here as B. n. sp. Animals in this clade form the sister group to B. vitiensis rather than other B. fasciatus. We herein describe this clade as a new species of Brachylophus based on molecular and morphological data. With only one exception, every island is home to one or more unique haplotypes. We discuss alternative biogeographic hypotheses to explain their distribution in the Pacific and the difficulties of distinguishing these. Together, our molecular and taxonomic results have important implications for future conservation initiatives for the Pacific iguanas.

Genetic impact of a severe El Niño event on Galápagos marine iguanas (Amblyrhynchus cristatus)

The El Niño-Southern Oscillation (ENSO) is a major source of climatic disturbance, impacting the dynamics of ecosystems worldwide. Recent models predict that human-generated rises in green-house gas levels will cause an increase in the strength and frequency of El Niño warming events in the next several decades, highlighting the need to understand the potential biological consequences of increased ENSO activity. Studies have focused on the ecological and demographic implications of El Niño in a range of organisms, but there have been few systematic attempts to measure the impact of these processes on genetic diversity in populations. Here, we evaluate whether the 1997–1998 El Niño altered the genetic composition of Galápagos marine iguana populations from eleven islands, some of which experienced mortality rates of up to 90% as a result of El Niño warming. Specifically, we measured the temporal variation in microsatellite allele frequencies and mitochondrial DNA diversity (mtDNA) in samples collected before (1991/1993) and after (2004) the El Niño event. Based on microsatellite data, only one island (Marchena) showed signatures of a genetic bottleneck, where the harmonic mean of the effective population size (N_e) was estimated to be less than 50 individuals during the period between samplings. Substantial decreases in mtDNA variation between time points were observed in populations from just two islands (Marchena and Genovesa). Our results suggests that, for the majority of islands, a single, intense El Niño event did not reduce marine iguana populations to the point where substantial neutral genetic diversity was lost. In the case of Marchena, simultaneous changes to both nuclear and mitochondrial DNA variation may also be the result of a volcanic eruption on the island in 1991. Therefore, studies that seek to evaluate the genetic impact of El Niño must also consider the confounding or potentially synergistic effect of other environmental and biological forces shaping populations.

Gene expression and tissue distribution of cytoglobin and myoglobin in the Amphibia and Reptilia: possible compensation of myoglobin with cytoglobin in skeletal muscle cells of anurans that lack the myoglobin gene

Cytoglobin (Cygb), a recently discovered vertebrate cytoplasmic heme-binding globin, is considered to be in a clade with vertebrate myoglobin (Mb), which is exclusively distributed in the cytoplasm of cardiac and skeletal muscles as an oxygen storage protein. GenBank databases (NCBI and JGI) and gene synteny analyses showed the absence of the Mb gene (mb) in two anuran amphibians, Xenopus laevis and X. tropicalis. Here we conducted comparative studies on the gene expression and tissue distribution of Cygb and Mb in anuran and reptilian tissues. Cygb and Mb genes were cloned from a reptile, iguana (Iguana iguana). Two types of cygb (cygb-1 and -2) were cloned, with lengths of 1066 and 1034 bp, and 196 and 193 amino acid residues, respectively. Their nucleotide and amino acid sequence identities were 90 and 87%, respectively. The Mb gene covered 1416 bp with an open reading frame of 465 bp, giving rise to a 154 amino acid protein. The distal ligand-binding histidine at E7, the proximal heme-binding histidine at F8, and the phenylalanine residue at CD1 were conserved in Mb and Cygb. The nucleotide and amino acid sequence identity of I. iguana cygb-1 and cygb-2 against X. laevis cygb were approximately 67% and 65%, respectively. RT-PCR demonstrated that X. laevis cygb was uniquely expressed in the heart and skeletal muscles, and faintly in the liver and spleen, which was quite contrasted with Iguana and the other vertebrates, where mb is exclusively expressed in the heart and skeletal muscles. Immunohistochemical analyses showed the distribution of Cygb in the cytoplasm of skeletal muscle cells. Interestingly, Cygb in the heart was localized in the nuclei. Considering the absence of mb in the Anura, we hypothesize that Cygb in muscle cells of anurans compensates for the lack of Mb for the storage and intracellular transportation of oxygen.

Identification and characterization of a squamate reptilian amelogenin gene: Iguana iguana

As the principal components of the developing tooth enamel matrix, amelogenins play a significant role in tooth enamel formation and organization. In order to elucidate the structure and function of amelogenins in the evolution of enamel, we have selected the Iguana iguana as a squamate model organism. Here we report the first complete squamate amelogenin sequence available as of yet and document unique features of Iguana amelogenins and enamel. Transmission electron microscopy documented randomly oriented Iguana enamel crystals during the elongation phase compared with organized enamel crystal patterns at comparable stages in mammals. Sequencing of PCR amplified products revealed a full-length I. iguana amelogenin cDNA containing 877 nucleotides with a 564 nucleotide coding sequence encoding 187 amino acids. The homologies of the newly discovered I. iguana amelogenin amino acid sequence with the published mouse, caiman (Palaeosuchus), and snake (Elaphe) amelogenin were 41.3%, 53.5%, and 55.5%, respectively. On Western blots one major protein with a molecular weight of 24 kDa, and two minor proteins with molecular weights of 28 and 13.5 kDa, respectively, were detected based on the cross-reactivity of antisera against recombinant Rana pipiens amelogenin proteins. Sequence analysis revealed a moderate sequence homology between mammalian and reptilian amelogenin genes. A significant alteration was the deletion of the hydrophilic GSP sequence from exon 3 in the mouse sequence resulting in a conversion to a hydrophobic region in Iguana. Together, these findings identified a novel amelogenin cDNA sequence in the squamate reptilian I. iguana and functional implications for the evolution of amelogenins and enamel in squamates.

Evolution of body size in Galapagos marine iguanas

Body size is one of the most important traits of organisms and allows predictions of an individual's morphology, physiology, behaviour and life history. However, explaining the evolution of complex traits such as body size is difficult because a plethora of other traits influence body size. Here I review what we know about the evolution of body size in a group of island reptiles and try to generalize about the mechanisms that shape body size. Galapagos marine iguanas occupy all 13 larger islands in this Pacific archipelago and have maximum island body weights between 900 and 12 000g. The distribution of body sizes does not match mitochondrial clades, indicating that body size evolves independently of genetic relatedness. Marine iguanas lack intra- and inter-specific food competition and predators are not size-specific, discounting these factors as selective agents influencing body size. Instead I hypothesize that body size reflects the trade-offs between sexual and natural selection. We found that sexual selection continuously favours larger body sizes. Large males establish display territories and some gain over-proportional reproductive success in the iguanas' mating aggregations. Females select males based on size and activity and are thus responsible for the observed mating skew. However, large individuals are strongly selected against during El Niño-related famines when dietary algae disappear from the intertidal foraging areas. We showed that differences in algae sward ('pasture') heights and thermal constraints on large size are causally responsible for differences in maximum body size among populations. I hypothesize that body size in many animal species reflects a trade-off between foraging constraints and sexual selection and suggest that future research could focus on physiological and genetic mechanisms determining body size in wild animals. Furthermore, evolutionary stable body size distributions within populations should be analysed to better understand selection pressures on individual body size.

Intercontinental community convergence of ecology and morphology in desert lizards

Evolutionary ecologists have long debated the extent to which communities in similar environments but different geographic regions exhibit convergence. On the one hand, if species' adaptations and community structure are determined by environmental features, convergence would be expected. However, if historical contingencies have long-lasting effects convergence would be unlikely. Most studies to date have emphasized the differences between communities in similar environments and little quantitative evidence for convergence exists. The application of comparative phylogenetic methods to ecological studies provides an opportunity to further investigate hypotheses of convergence. We compared the evolutionary patterns of structural ecology and morphology of 42 species of iguanian lizards from deserts of Australia and North America. Using a comparative approach, we found that evolutionary convergence of ecology and morphology occurs both in overall, community-wide patterns and in terms of pairs of highly similar intercontinental pairs of species. This result indicates that in these desert lizards, deterministic adaptive evolution shapes community patterns and overrides the historical contingencies unique to particular lineages.

Inter-population variation of carotenoids in Galápagos land iguanas (Conolophus subcristatus)

Carotenoids have received much attention from biologists because of their ecological and evolutionary implications in vertebrate biology. We sampled Galápagos land iguanas (Conolophus subcristatus) to investigate the types and levels of blood carotenoids and the possible factors affecting inter-population variation. Blood samples were collected from populations from three islands within the species natural range (Santa Cruz, Isabela, and Fernandina) and one translocated population (Venecia). Lutein and zeaxanthin were the predominant carotenoids found in the serum. In addition, two metabolically modified carotenoids (anhydrolutein and 3'-dehydrolutein) were also identified. Differences in the carotenoid types were not related to sex or locality. Instead, carotenoid concentration varied across the localities, it was higher in females, and it was positively correlated to an index of body condition. Our results suggest a possible sex-related physiological role of xanthophylls in land iguanas. The variation in the overall carotenoid concentration between populations seems to be related to the differences in local abundance and type of food within and between islands.

Characterization of a ribonuclease gene and encoded protein from the reptile, Iguana iguana

In this work we identify an intronless open reading frame encoding an RNase A ribonuclease from genomic DNA from the Iguana iguana IgH2 cell line. The iguana RNase is expressed primarily in pancreas, and represents the majority of the specific enzymatic activity in this tissue. The encoded sequence shares many features with its better-known mammalian counterparts including the crucial His12, Lys40 and His114 catalytic residues and efficient hydrolytic activity against yeast tRNA substrate (k(cat)/K(m)=6 x 10(4) M(-1) s(-1)), albeit at a reduced pH optimum (pH 6.0). Although the catalytic activity of the iguana RNase is not diminished by human placental RI, iguana RNase is not bactericidal nor is it cytotoxic even at micromolar concentrations. Phylogenetic analysis indicates moderate (46%) amino acid sequence similarity to a pancreatic RNase isolated from Chelydra serpentina (snapping turtle) although no specific relationship could be determined between these RNases and the pancreatic ribonucleases characterized among mammalian species. Further analysis of ribonucleases from non-mammalian vertebrate species is needed in order to define relationships and lineages within the larger RNase A gene superfamily.

Tempo and mode of evolutionary radiation in iguanian lizards

Identification of general properties of evolutionary radiations has been hindered by the lack of a general statistical and phylogenetic approach applicable across diverse taxa. We present a comparative analytical framework for examining phylogenetic patterns of diversification and morphological disparity with data from four iguanian-lizard taxa that exhibit substantially different patterns of evolution. Taxa whose diversification occurred disproportionately early in their evolutionary history partition more of their morphological disparity among, rather than within, subclades. This inverse relationship between timing of diversification and morphological disparity within subclades may be a general feature that transcends the historically contingent properties of different evolutionary radiations.

Sampling strategies for delimiting species: genes, individuals, and populations in the Liolaemus elongatus-kriegi complex (Squamata: Liolaemidae) in Andean-Patagonian South America

Recovery of evolutionary history and delimiting species boundaries in widely distributed, poorly known groups requires extensive geographic sampling, but sampling regimes are difficult to design a priori because evolutionary diversity is often "hidden" by inadequate taxonomy. Large data sets are needed, and these provide unique challenges for analysis when they span intra- and interspecific levels of divergence. However, protocols have been designed to combine methods of analysis for DNA sequences that exhibit both very shallow and relatively deeper divergences. In this study, we combined several tree-based phylogeny reconstruction methods with nested-clade analysis to extract maximum historical signal at various levels in the poorly known Liolaemus elongatus-kriegi lizard complex in temperate South America. We implemented a recently descrirbed tree-based protocol for DNA sequences to test for species boundaries, and we propose modifications to accommodate large data sets and gene regions with heterogeneous substitution rates. Combining haplotype trees with nested-clade analyses allowed testing of species boundaries on the basis of a priori defined criteria. The results obtained suggest that the number of putative species in the L. elongatus-kriegi complex could be doubled. We discuss these findings in the context of the advantages and limitations of a combined approach for retrieval of maximum historical information in large data sets and with reference to the yet formidable unresolved issues of sampling strategies.

Comparative cytogenetics and supernumerary chromosomes in the Brazilian lizard genus Enyalius (Squamata, Polychrotidae)

Cytogenetical analyses based on conventional and differential staining were performed for the first time on five species of the Brazilian lizard genus Enyalius: E. bibronii, E. bilineatus, E. iheringii, E. leechii, and E. perditus. The species share a similar 2n = 36 (12M + 24m) karyotype, comprised of 12 metacentric or submetacentric macrochromosomes, except for an acrocentric pair 6 that characterizes E. bibronii. The 24 microchromosomes were acrocentrics, but in E. perditus two meta/submetacentric microchromosome pairs were unambiguously identified. Karyotypes with 2n = 37 and 2n = 37/38 chromosomes were also observed in some specimens of E. bilineatus as a result of the presence of supernumerary chromosomes (Bs). Ag-NORs were always located at the distal region of the long arm of the submetacentric pair 2. The constitutive heterochromatin was mostly restricted to the pericentromeric regions of some macrochromosomes and microchromosomes. A XX:XY mechanism of sex determination with a dot-like Y microchromosome occurs in E. bilineatus, E. leechii, and E. perditus.

Geographic variation in sexual selection among populations of an iguanid lizard, Sauromalus obesus (=ater)

Geographic variation in selection pressures may result in population divergence and speciation, especially if sexual selection varies among populations. Yet spatial variation in targets and intensity of sexual selection is well studied in only a few species. Even more rare are simultaneous studies of multiple populations combining observations from natural settings with controlled behavioral experiments. We investigated how sexual selection varies among populations of the chuckwalla, Sauromalus obesus. Chuckwallas are sexually dimorphic in color, and males vary in coloration among populations. Using field observations and multiple regression techniques, we investigated how sexual selection acts on various male traits in three populations in which males differed in coloration. The influence of sexual selection on male coloration was then investigated in more detail using controlled experiments. Results from field observations indicate that phenotypic selection was acting on territory quality in all three populations. In two populations, selection was also acting either directly or indirectly on male coloration. Male color likely functions as an indicator of food resources to females because male color is based partly on carotenoid pigments. In controlled experiments, significantly more females from these two populations chose males with brighter colors over dull males, a result consistent with studies on carotenoid pigments in other taxa. In a third population, no evidence of sexual selection on male coloration was found in either the field study or controlled experiment. Lack of female preferences for male color in this population, in which chuckwalla densities are low and home ranges are large, may result from searching costs to females.

The mitochondrial genomes of the iguana (Iguana iguana) and the caiman (Caiman crocodylus): implications for amniote phylogeny

The complete mitochondrial genomes of two reptiles, the common iguana (Iguana iguana) and the caiman (Caiman crocodylus), were sequenced in order to investigate phylogenetic questions of tetrapod evolution. The addition of the two species allows analysis of reptilian relationships using data sets other than those including only fast-evolving species. The crocodilian mitochondrial genomes seem to have evolved generally at a higher rate than those of other vertebrates. Phylogenetic analyses of 2889 amino-acid sites from 35 mitochondrial genomes supported the bird-crocodile relationship, lending no support to the Haematotherma hypothesis (with birds and mammals representing sister groups). The analyses corroborated the view that turtles are at the base of the bird-crocodile branch. This position of the turtles makes Diapsida paraphyletic. The origin of the squamates was estimated at 294 million years (Myr) ago and that of the turtles at 278 Myr ago. Phylogenetic analysis of mammalian relationships using the additional outgroups corroborated the Marsupionta hypothesis, which joins the monotremes and the marsupials to the exclusion of the eutherians.

Phylogeography of the Caribbean rock iguana (Cyclura): implications for conservation and insights on the biogeographic history of the West Indies

The Caribbean rock iguana, Cyclura, has had an unstable intrageneric taxonomy and an unclear phylogenetic position within the family Iguanidae. We use mtDNA sequence data to address these issues and explore the phylogeographic history of the genus. ND4 to leucine tRNA sequence data were collected from multiple individuals of each of the eight species of Cyclura (including 15 of 16 subspecies) and from four localities of Iguana iguana (representative of this species' broad geographic range). This data set was combined with sequence data from Sites et al. (1996, Mol. Biol. Evol. 13, 1087-1105) and analyzed under maximum-parsimony and maximum-likelihood optimization criteria. The ND4 region provided good resolution for the majority of nodes, as indicated by high bootstrap support. In agreement with several recent molecular studies, Cyclura is recovered as monophyletic and is not closely related to any other genus, whereas Iguana is strongly supported as the sister taxon to Sauromalus. This result is statistically more likely than other published hypotheses of Iguanid relationships. Cyclura shows a southeast to northwest speciation sequence in the Caribbean, with the most ancient lineage on the Puerto Rican Bank. The amount of interspecific sequence divergence within Cyclura (maximum 11.4%) is very high in comparison to data from other iguanid taxa at this locus, suggesting that this group either has been in the Caribbean for a very long time or has gone through a very rapid rate of evolution at this locus. Using dates from other published studies, we calculate a molecular clock that suggests that Cyclura colonized the Caribbean between 15 and 35 mya. Several questions regarding subspecific taxonomy are raised in the analysis and await further investigation using a more rapidly evolving marker such as nuclear microsatellites.

War of the Iguanas: conflicting molecular and morphological phylogenies and long-branch attraction in iguanid lizards

Recent studies based on different types of data (i.e., morphology, molecules) have found strongly conflicting phylogenies for the genera of iguanid lizards but have been unable to explain the basis for this incongruence. We reanalyze published data from morphology and from the mitochondrial ND4, cytochrome b, 12S, and 16S genes to explore the sources of incongruence and resolve these conflicts. Much of the incongruence centers on the genus Cyclura, which is the sister taxon of Iguana, according to parsimony analyses of the morphology and the ribosomal genes, but is the sister taxon of all other Iguanini, according to the protein-coding genes. Maximum likelihood analyses show that there has been an increase in the rate of nucleotide substitution in Cyclura in the two protein-coding genes (ND4 and cytochrome b), although this increase is not as clear when parsimony is used to estimate branch lengths. Parametric simulations suggest that Cyclura may be misplaced by the protein-coding genes as a result of long-branch attraction; even when Cyclura and Iguana are sister taxa in a simulated phylogeny, Cyclura is still placed as the basal member of the Iguanini by parsimony analysis in 55% of the replicates. A similar long-branch attraction problem may also exist in the morphological data with regard to the placement of Sauromalus with the Galápagos iguanas (Amblyrhynchus and Conolophus). The results have many implications for the analysis of diverse data sets, the impact of long branches on parsimony and likelihood methods, and the use of certain protein-coding genes in phylogeny reconstruction.

Molecular tests of phylogenetic taxonomies: a general procedure and example using four subfamilies of the lizard family Iguanidae

A general procedure is described for examining when results of molecular phylogenetic analyses warrant formal revision of taxonomies constructed using morphological characters. We illustrate this procedure with tests of monophyly for four subfamilies in the lizard family Iguanidae using 1561 aligned base positions (838 phylogenetically informative) of mitochondrial DNA sequences, representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Ten new sequences ranging in length from 1732 to 1751 bases are compared with 12 previously reported sequences and 67 morphological characters (54 phylogenetically informative) from the literature. New morphological character states are provided for Sator. Phylogenies derived from the molecular and combined data are in agreement but both conflict with phylogenetic inferences from the morphological data alone. Strong support is found for the monophyly of the subfamilies Crotaphytinae and Phrynosomatinae. Monophyly of the Iguaninae is weakly supported in each analysis. All analyses suggest that the Tropidurinae is not monophyletic but the hypothesis of monophyly cannot be rejected. A phylogenetic taxonomy is proposed in which the Tropidurinae* is maintained as a metataxon (denoted with an asterisk), for which monophyly has not been demonstrated. Within the Phrynosomatinae, the close relationship of Sator and Sceloporus is questioned and an alternative hypothesis in which Sator is the sister taxon to a clade comprising Petrosaurus, Sceloporus, and Urosaurus is presented. Statistical tests of monophyly provide a powerful way to evaluate support for taxonomic groupings. Use of the metataxon prevents premature taxonomic rearrangements where support is lacking.

A new derived and highly polymorphic chromosomal race of Liolaemus monticola (Iguanidae) from the 'Norte Chico' of Chile

A multiple Robertsonian fission chromosomal race of the Liolaemus monticola complex in Chile is described and is shown to be the most derived and the most complex among the Liolaemus examined thus far. The 29 karyotyped lizards analysed from the locality of Mina Hierro Viejo, Petorca, Provincia de Valparaiso, Chile, exhibited a diploid chromosomal number ranging from 42 to 44, and several polymorphisms. The polymorphisms included: a pair 1 fission; a pair 2 fission plus a pericentric inversion in one of the fission products, which moved the NOR and satellite from the tip of the long arm of the metacentric 2 to the short arm of the fission product; a fission in pair 3; a polymorphism for an enlarged chromosome pair 6; and a polymorphism for a pericentric inversion in pair 7. This population is fixed for a fission of chromosome pair 4. A total of 76% of the lizards analysed were polymorphic for one or more pairs of chromosomes. We have compared these data with other Liolaemus monticola chromosomal races and calculated the Hardy-Weinberg ratios for the polymorphic chromosome pairs in this Multiple-Fission race. Karyotypic differences between the Northern (2n = 38-40) and the Multiple-Fission (2n = 42-44) races were attributed mainly to Robertsonian fissions, an enlarged chromosome and pericentric inversions involving the macrochromosomes and one microchromosome pair.

Evolutionary age of the Galápagos iguanas predates the age of the present Galápagos islands

New geological findings suggest that volcanoes existed over the Galápagos hotspot long before today's islands emerged less than 5 million years ago. The evolution of some of Galápagos' biota might have taken place on these former islands. This study investigates the evolutionary history of two of the archipelagos' older vertebrate taxa, the endemic Galápagos marine and land iguana (genera Amblyrhynchus and Conolophus). Mitochondrial rDNA sequences (in total about one kilobase of the 12S and 16S genes) were obtained from all extant genera of the family Iguanidae and the outgroup Oplurus. The phylogenetic analyses suggest that the Galápagos iguanas are sister taxa. Rate comparisons between the iguanid sequences and a corresponding set of sequences from ungulates with known fossil ages date their separation time at 10 million years, or more. The results strengthen the hypothesis that extended speciation times in the Galápagos are possible and provide an estimate of the minimum time inhabited islands of the archipelago may have existed.

Character congruence and phylogenetic signal in molecular and morphological data sets: a case study in the living Iguanas (Squamata, Iguanidae)

The lizard family Iguanidae comprises eight living genera distributed throughout the New and Old World, and includes several island endemics. We reconstruct phylogenetic relationships among these genera using 90 previously published morphological characters, to which we add a molecular (mtDNA sequence) data set that includes 742 nucleotides of the ND4 gene and the complete sequences of the histidine, serine, and leucine tRNAs (217 nucleotides). Trees were initially constructed separately from these three data sets, and then tested for significant conflict in topologies that would suggest the influence of different evolutionary processes. The three data sets were then combined, and a single tree was obtained from the total evidence that permitted identification of potential sources of character incongruence. Several additional analyses of the combined data sets were repeated with sequential deletion of successive classes of homoplastic characters, and we show that the same single tree topology is recovered in most cases. However, part of the tree structure collapses when the matrix of combined characters is completely purged of all homoplastic characters. We argue that this extreme results in an unacceptable loss of phylogenetic information, and we present a single phylogenetic hypothesis for all living genera of iguanas. We show that this hypothesis is significantly more parsimonious than either of two previously published trees, and we discuss the evolution and biogeography of the Iguanidae based on the preferred hypothesis.

Y3 is the most conserved small RNA component of Ro ribonucleoprotein complexes in vertebrate species

YRNAs are small cytoplasmic RNAs that are components of the Ro ribonucleoprotein complex. This complex, which also includes the 60-kDa Ro protein, is a human autoantigen which is conserved among vertebrates, and is of unknown function. Multiple sequences with YRNA homology, known as YRNA-like sequences, have been detected in rabbit, mouse, duck, iguana and frog genomes with human Y cDNA probes. As judged by Northern blots of total RNA, however, not all of these genomic YRNA-like sequences are expressed. Complementary DNA and putative gene sequences for iguana Y3 (iY3) and iguana Y4 (iY4) Ro RNAs have been determined and used, along with previously sequenced human and frog Ro YRNA sequences, to construct the most likely Y3 and Y4 RNA secondary structures. The data presented indicate that Y3 is the most conserved Ro RNA, not only by its more consistent presence in other species, but also at the levels of sequence divergence and secondary structure similarity. The differences observed between the secondary structure solutions for the Y3 and Y4 Ro RNAs are consistent with the possibility that these RNAs perform different cellular functions.