Biogeographic origins of goannas (Varanidae): a molecular perspective

Fossil-informed biogeographic analysis suggests Eurasian regionalization in crown Squamata during the early Jurassic

Background

Squamata (lizards, snakes, and amphisbaenians) is a Triassic lineage with an extensive and complex biogeographic history, yet no large-scale study has reconstructed the ancestral range of early squamate lineages. The fossil record indicates a broadly Pangaean distribution by the end- Cretaceous, though many lineages (e.g., Paramacellodidae, Mosasauria, Polyglyphanodontia) subsequently went extinct. Thus, the origin and occupancy of extant radiations is unclear and may have been localized within Pangaea to specific plates, with potential regionalization to distinct Laurasian and Gondwanan landmasses during the Mesozoic in some groups.

Methods

We used recent tectonic models to code extant and fossil squamate distributions occurring on nine discrete plates for 9,755 species, with Jurassic and Cretaceous fossil constraints from three extinct lineages. We modeled ancestral ranges for crown Squamata from an extant-only molecular phylogeny using a suite of biogeographic models accommodating different evolutionary processes and fossil-based node constraints from known Jurassic and Cretaceous localities. We hypothesized that the best-fit models would not support a full Pangaean distribution (i.e., including all areas) for the origin of crown Squamata, but would instead show regionalization to specific areas within the fragmenting supercontinent, likely in the Northern Hemisphere where most early squamate fossils have been found.

Results

Incorporating fossil data reconstructs a localized origin within Pangaea, with early regionalization of extant lineages to Eurasia and Laurasia, while Gondwanan regionalization did not occur until the middle Cretaceous for Alethinophidia, Scolecophidia, and some crown Gekkotan lineages. While the Mesozoic history of extant squamate biogeography can be summarized as a Eurasian origin with dispersal out of Laurasia into Gondwana, their Cenozoic history is complex with multiple events (including secondary and tertiary recolonizations) in several directions. As noted by previous authors, squamates have likely utilized over-land range expansion, land-bridge colonization, and trans-oceanic dispersal. Tropical Gondwana and Eurasia hold more ancient lineages than the Holarctic (Rhineuridae being a major exception), and some asymmetries in colonization (e.g., to North America from Eurasia during the Cenozoic through Beringia) deserve additional study. Future studies that incorporate fossil branches, rather than as node constraints, into the reconstruction can be used to explore this history further.

Genetic variations and phylogenetic relationship of genus Uromastyx from Punjab Pakistan

During the present study, specimens were collected from selected sites of Cholistan desert and Kalabagh Game Reserve, Punjab province, Pakistan. Each captured specimen was tagged with voucher number and morphometric measurements were taken. The average snout to vent length was 172.559±1.40 mm and average weight was 92.1±1.30 g. The DNA of Uromastyx hardwickii was amplified and sequenced using 16S rRNA primer set. The obtained DNA sequence has shown reliable and clear species identification. After trimming ambiguous bases, the obtained 16S rRNA fragment was 520 bp while 16S rRNA fragments aligned with closely matched sequence from NCBI comprised of 510 bp. Closely matched sequences of genus Uromastyx were retrieved from NCBI in blast searches. Neighbour-joining tree of genus Uromastyx was constructed based on p-distance using MEGA X. The mean intraspecific variation was 0.095±0.01 while intraspecific variation was ranging from 0-1%. Similarly, interspecific variation of Uromastyx hardwikii with Saara asmussi, Uromastyx alfredschmidti, Uromastyx geyri, Uromastyx thomasi, Uromastyx alfredschmidti was 0-12%, 0-19%, 0-19%, 0-20%, 12-19% respectively. The newly produced DNA was submitted to NCBI and accession number was obtained (MW052563.1). Results of current study provided information about the molecular and morphological identification of Genus Uromastyx. In our recommendation, comprehensive molecular based identification of Pakistan's reptiles is required to report any new or subspecies from country.

Phylogenomics of Monitor Lizards and the Role of Competition in Dictating Body Size Disparity

Organismal interactions drive the accumulation of diversity by influencing species ranges, morphology, and behavior. Interactions vary from agonistic to cooperative and should result in predictable patterns in trait and range evolution. However, despite a conceptual understanding of these processes, they have been difficult to model, particularly on macroevolutionary timescales and across broad geographic spaces. Here, we investigate the influence of biotic interactions on trait evolution and community assembly in monitor lizards (Varanus). Monitors are an iconic radiation with a cosmopolitan distribution and the greatest size disparity of any living terrestrial vertebrate genus. Between the colossal Komodo dragon Varanus komodoensis and the smallest Australian dwarf goannas, Varanus length and mass vary by multiple orders of magnitude. To test the hypothesis that size variation in this genus was driven by character displacement, we extended existing phylogenetic comparative methods which consider lineage interactions to account for dynamic biogeographic history and apply these methods to Australian monitors and marsupial predators. Incorporating both exon-capture molecular and morphological data sets we use a combined evidence approach to estimate the relationships among living and extinct varaniform lizards. Our results suggest that communities of Australian Varanus show high functional diversity as a result of continent-wide interspecific competition among monitors but not with faunivorous marsupials. We demonstrate that patterns of trait evolution resulting from character displacement on continental scales are recoverable from comparative data and highlight that these macroevolutionary patterns may develop in parallel across widely distributed sympatric groups.[Character displacement; comparative methods; phylogenetics; trait evolution; Varanus.].

A molecular approach to identifying the relationship between resource use and availability in Eurasian otters (Lutra lutra)

In South Korea, the Eurasian otter (Lutra lutra (Linnaeus, 1758)), a semi-aquatic carnivore, is found mainly in lower order streams that tend to have a low abundance of preferred prey fish species. To investigate the relationship between resource use and availability, we used DNA barcoding to identify otter diet items in 24 otter spraints (faeces) from 16 sites along the Nakdong River basin from 4 to 6 June 2014. At these sites fish availability was assessed using scoop nets and casting nets. Fish formed the bulk of otter diet, which included also frogs, mammals, and reptiles. By DNA barcoding (success rate: 72.38%), we identified 79 prey items from 105 bone remains. The diet comprised mostly fish, but frogs, mammals, and reptiles were also identified. The fish fauna and otter diet composition differed significantly. Across the study sites, members of the Cyprinidae dominated in netted samples, but occurred less frequently in otter diet. Because most Cyprinidae are fast swimmers, otters also fed on benthic fishes and frogs, suggesting limited foraging flexibility in otters and specialization on more slowly moving prey.

Revision of Varanus marathonensis (Squamata, Varanidae) based on historical and new material: morphology, systematics, and paleobiogeography of the European monitor lizards

Monitor lizards (genus Varanus) inhabited Europe at least from the early Miocene to the Pleistocene. Their fossil record is limited to about 40 localities that have provided mostly isolated vertebrae. Due to the poor diagnostic value of these fossils, it was recently claimed that all the European species described prior to the 21st century are not taxonomically valid and a new species, Varanus amnhophilis, was erected on the basis of fragmentary material including cranial elements, from the late Miocene of Samos (Greece). We re-examined the type material of Varanus marathonensis Weithofer, 1888, based on material from the late Miocene of Pikermi (Greece), and concluded that it is a valid, diagnosable species. Previously unpublished Iberian material from the Aragonian (middle Miocene) of Abocador de Can Mata (Vallès-Penedès Basin, Barcelona) and the Vallesian (late Miocene) of Batallones (Madrid Basin) is clearly referable to the same species on a morphological basis, further enabling to provide an emended diagnosis for this species. Varanus amnhophilis appears to be a junior subjective synonym of V. marathonensis. On the basis of the most complete fossil Varanus skeleton ever described, it has been possible to further resolve the internal phylogeny of this genus by cladistically analyzing 80 taxa coded for 495 morphological and 5729 molecular characters. Varanus marathonensis was a large-sized species distributed at relatively low latitudes in both southwestern and southeastern Europe from at least MN7+8 to MN12. Our cladistic analysis nests V. marathonensis into an eastern clade of Varanus instead of the African clade comprising Varanus griseus, to which it had been related in the past. At least two different Varanus lineages were present in Europe during the Neogene, represented by Varanus mokrensis (early Miocene) and V. marathonensis (middle to late Miocene), respectively.

Forensically informative nucleotide sequencing (FINS) for the first time authentication of Indian Varanus species: implication in wildlife forensics and conservation

Monitor lizards are Varanus species widely distributed, endangered reptile in the IUCN red data list. In India, based on the morphological and ecological characteristic, it is divided into four species viz. Bengal monitor lizard, Yellow monitor lizard, Desert monitor lizard and Water monitor lizard. These four species listed as Schedule I species in Indian Wildlife (Protection) Act 1972. This paper first attempt to present Forensically Informative Nucleotide Sequencing (FINS) for the Indian Varanus based on three mitochondrial genes. The molecular framework will be useful for the identification of Indian Varanus species and trade products derived from monitors and as such, have important applications for wildlife management and conservation. Here, we used known 14 individual skin pieces of four species of monitor lizards; the partial fragment of three mitochondrial genes (Cyt b, 12S rRNA, and 16S rRNA) were amplified for genetic study. In Cyt b, 12S rRNA and 16s rRNA, we observed, 5, 5 and 4 Haplotypes; 71, 69, and 43 Variables sites; 90, 89, and 50 Parsimony Informative sites within four species of Indian monitor lizards, respectively. Despite it, the nucleotide composition was T 26.4, C 32.8, A 29.2 and G11.6; T 18.8, C 29.7, A 34.0 and G 17.5; T 21.7, C 27.3, A 32.5 and G 18.5 in Cyt b, 12S rRNA and 16S rRNA, respectively. The neighbor joining phylogenetic tree and maximum parsimony tree of three mitochondrial genes, showed similar results and reveal that, there are two major clades are present in Indian monitor lizards.

Evaluating the effects of laboratory protocols on eDNA detection probability for an endangered freshwater fish

The effectiveness and accuracy of detection using environmental DNA (eDNA) is dependent on understanding the influence laboratory methods such as DNA extraction and PCR strategies have on detection probability. Ideally choice of sampling and extraction method will maximize eDNA yield and detection probability. Determining the survey effort required to reach a satisfactory detection probability (via increased PCR replicates or more sampling) could compensate for a lower eDNA yield if the sampling and extraction method has other advantages for a study, species or system. I analysed the effect of three different sampling and extraction methods on eDNA yield, detection probability and PCR replication for detecting the endangered freshwater fish Macquaria australasica from water samples. The impact of eDNA concentration, PCR strategy, target amplicon size and two marker regions: 12S (a mitochondrial gene) and 18S (a nuclear gene) was also assessed. The choice of sampling and extraction method and PCR strategy, rather than amplicon size and marker region, had the biggest effect on detection probability and PCR replication. The PCR replication effort required to achieve a detection probability of 0.95, ranged from 2 to 6 PCR replicates depending on the laboratory method used. As all methods yielded eDNA from which M. australasica was detected using the three target amplicons, differences in eDNA yield and detection probability between the three methods could be mitigated by determining the appropriate PCR replication effort. Evaluating the effect sampling and extraction methods will have on the detection probability and determining the laboratory protocols and PCR replication required to maximize detection and minimize false positives and negatives is a useful first step for eDNA occupancy studies.

Highly species-specific centromeric repetitive DNA sequences in lizards: molecular cytogenetic characterization of a novel family of satellite DNA sequences isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota)

Two novel repetitive DNA sequences, VSAREP1 and VSAREP2, were isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota) and characterized using molecular cytogenetics. The respective lengths and guanine-cytosine (GC) contents of the sequences were 190 bp and 57.5% for VSAREP1 and 185 bp and 59.7% for VSAREP2, and both elements were tandemly arrayed as satellite DNA in the genome. VSAREP1 and VSAREP2 were each located at the C-positive heterochromatin in the pericentromeric region of chromosome 2q, the centromeric region of chromosome 5, and 3 pairs of microchromosomes. This suggests that genomic compartmentalization between macro- and microchromosomes might not have occurred in the centromeric repetitive sequences of V. salvator macromaculatus. These 2 sequences did only hybridize to genomic DNA of V. salvator macromaculatus, but no signal was observed even for other squamate reptiles, including Varanus exanthematicus, which is a closely related species of V. salvator macromaculatus. These results suggest that these sequences were differentiated rapidly or were specifically amplified in the V. salvator macromaculatus genome.

Highly differentiated ZW sex microchromosomes in the Australian Varanus species evolved through rapid amplification of repetitive sequences

Transitions between sex determination systems have occurred in many lineages of squamates and it follows that novel sex chromosomes will also have arisen multiple times. The formation of sex chromosomes may be reinforced by inhibition of recombination and the accumulation of repetitive DNA sequences. The karyotypes of monitor lizards are known to be highly conserved yet the sex chromosomes in this family have not been fully investigated. Here, we compare male and female karyotypes of three Australian monitor lizards, Varanus acanthurus, V. gouldii and V. rosenbergi, from two different clades. V. acanthurus belongs to the acanthurus clade and the other two belong to the gouldii clade. We applied C-banding and comparative genomic hybridization to reveal that these species have ZZ/ZW sex micro-chromosomes in which the W chromosome is highly differentiated from the Z chromosome. In combination with previous reports, all six Varanus species in which sex chromosomes have been identified have ZZ/ZW sex chromosomes, spanning several clades on the varanid phylogeny, making it likely that the ZZ/ZW sex chromosome is ancestral for this family. However, repetitive sequences of these ZW chromosome pairs differed among species. In particular, an (AAT)n microsatellite repeat motif mapped by fluorescence in situ hybridization on part of W chromosome in V. acanthurus only, whereas a (CGG)n motif mapped onto the W chromosomes of V. gouldii and V. rosenbergi. Furthermore, the W chromosome probe for V. acanthurus produced hybridization signals only on the centromeric regions of W chromosomes of the other two species. These results suggest that the W chromosome sequences were not conserved between gouldii and acanthurus clades and that these repetitive sequences have been amplified rapidly and independently on the W chromosome of the two clades after their divergence.

Fossil-calibrated phylogeny and historical biogeography of Southeast Asian water monitors (Varanus salvator Complex)

We utilize robust geographical genetic sampling, a multilocus dataset, a new synthesis of numerous fossil calibration points, a time-calibrated phylogeny, and the Dispersal-Extinction-Cladogenesis model to test the prediction that widespread Southeast Asian water monitor species initially diversified on the Asian mainland and subsequently invaded the island archipelagos of the Philippines, Sundaland, and Wallacea. Our results strongly contradict these expectations and instead infer an initial water monitor radiation of range-restricted but highly divergent evolutionary lineages (now recognized as endemic species) in one archipelago around 3.6 mya, followed by an out-of-the-Philippines reinvasion of the mainland (2.2 mya), resulting in a few, widespread species that now inhabit most the islands of the Sunda Shelf and the Southeast Asian mainland as far north as Myanmar, as well as an out-of-the-Philippines invasion of Sulawesi (2.1 mya). Our analyses both confirm the importance of island archipelagos as drivers of diversification for mainland biodiversity and emphasize the global evolutionary significance and conservation priority of the Philippines for understanding processes of diversification in island archipelagos.

Molecular evidence for an Asian origin of monitor lizards followed by Tertiary dispersals to Africa and Australasia

Monitor lizards are emblematic reptiles that are widely distributed in the Old World. Although relatively well studied in vertebrate research, their biogeographic history is still controversial. We constructed a molecular dataset for 54 anguimorph species, including representatives of all families with detailed sampling of the Varanidae (38 species). Our results are consistent with an Asian origin of the Varanidae followed by a dispersal to Africa 41 (49-33) Ma, possibly via an Iranian route. Another major event was the dispersal of monitors to Australia in the Late Eocene-Oligocene 32 (39-26) Ma. This divergence estimate adds to the suggestion that Australia was colonized by several squamate lineages prior to the collision of the Australian plate with the Asian plate starting 25 Ma.

Earliest example of a giant monitor lizard (Varanus, Varanidae, Squamata)

BACKGROUND

Varanidae is a clade of tiny (<20 mm pre-caudal length [PCL]) to giant (>600 mm PCL) lizards first appearing in the Cretaceous. True monitor lizards (Varanus) are known from diagnostic remains beginning in the early Miocene (Varanus rusingensis), although extremely fragmentary remains have been suggested as indicating earlier Varanus. The paleobiogeographic history of Varanus and timing for origin of its gigantism remain uncertain.

METHODOLOGY/PRINCIPAL FINDINGS

A new Varanus from the Mytilini Formation (Turolian, Miocene) of Samos, Greece is described. The holotype consists of a partial skull roof, right side of a braincase, partial posterior mandible, fragment of clavicle, and parts of six vertebrae. A cladistic analysis including 83 taxa coded for 5733 molecular and 489 morphological characters (71 previously unincluded) demonstrates that the new fossil is a nested member of an otherwise exclusively East Asian Varanus clade. The new species is the earliest-known giant (>600 mm PCL) terrestrial lizard. Importantly, this species co-existed with a diverse continental mammalian fauna.

CONCLUSIONS/SIGNIFICANCE

The new monitor is larger (longer) than 99% of known fossil and living lizards. Varanus includes, by far, the largest limbed squamates today. The only extant non-snake squamates that approach monitors in maximum size are the glass-snake Pseudopus and the worm-lizard Amphisbaena. Mosasauroids were larger, but exclusively marine, and occurred only during the Late Cretaceous. Large, extant, non-Varanus, lizards are limbless and/or largely isolated from mammalian competitors. By contrast, our new Varanus achieved gigantism in a continental environment populated by diverse eutherian mammal competitors.

Vertebral morphometrics in Varanus

The forces acting on the vertebral column of varanid lizards differ greatly during swimming and walking. To examine the long-term impact of these forces, the dorsal vertebrae of terrestrial and aquatic species of Varanus were compared using 3-D laser scanning and morphometric analysis. There were significant differences between the two groups in vertebral anatomical features, particularly in the articular surfaces. Further analysis demonstrated that the dorsal vertebrae could be significantly divided into three groups of nearly equal size: the sternal group (dorsal vertebrae 1–7), the middle group (dorsal vertebrae 8–14), and the pelvic group (dorsal vertebrae 15–22). Within each of these groups there was significant differences between the terrestrial and aquatic species; these differences were least in the more conserved sternal region and greatest in the highly variable pelvic region. The results suggest that vertebral morphometrics can be used as a tool to further delineate the habitat preferences of extinct varanoid lizards, and could serve to highlight the functional transitions between terrestrial and aquatic species.

Rampant horizontal transfer of SPIN transposons in squamate reptiles

Transposable elements (TEs) are highly abundant in the genome and capable of mobility, two properties that make them particularly prone to transfer horizontally between organisms. Although the impact of horizontal transfer (HT) of TEs is well recognized in prokaryotes, the frequency of this phenomenon and its contribution to genome evolution in eukaryotes remain poorly appreciated. Here, we provide evidence that a DNA transposon called SPIN has colonized the genome of 17 species of reptiles representing nearly every major lineage of squamates, including 14 families of lizards, snakes, and amphisbaenians. Slot blot analyses indicate that SPIN has amplified to high copy numbers in most of these species, ranging from 2,000-28,000 copies per haploid genome. In contrast, we could not detect the presence of SPIN in any of the turtles (seven species from seven families) and crocodiles (four species) examined. Genetic distances between SPIN sequences from species belonging to different squamate families are consistently very low (average = 0.1), considering the deep evolutionary divergence of the families investigated (most are >100 My diverged). Furthermore, these distances fall below interfamilial distances calculated for two genes known to have evolved under strong functional constraint in vertebrates (RAG1, average = 0.24 and C-mos, average = 0.27). These data, combined with phylogenetic analyses, indicate that the widespread distribution of SPIN among squamates is the result of at least 13 independent events of HTs. Molecular dating and paleobiogeographical data suggest that these transfers took place during the last 50 My on at least three different continents (North America, South America and, Africa). Together, these results triple the number of known SPIN transfer events among tetrapods, provide evidence for a previously hypothesized transoceanic movement of SPIN transposons during the Cenozoic, and further underscore the role of HT in the evolution of vertebrate genomes.

Carbon source accounting for fish using combined DNA and stable isotope analyses in a regulated lowland river weir pool

Determining the source and flow of carbon, energy and nutrients through food webs is essential for understanding ecological connectivity and thus determining the impact of management practices on biodiversity. We combined DNA sequencing, microarrays and stable isotope analyses to test whether this approach would allow us to resolve the carbon flows through food webs in a weir pool on the lower Murray River, a highly impacted, complex and regulated ecosystem in southern Australia. We demonstrate that small fish in the Murray River consume a wide range of food items, but that a significant component of carbon and nitrogen entering the food web during dry periods in summer, but not spring, is derived from nonconventional sources other than in-channel primary producers. This study also showed that isotopic analyses alone cannot distinguish food sources and that a combined approach is better able to elucidate food-consumer dynamics. Our results highlight that a major river ecosystem, stressed by reduced environmental flows, can rapidly undergo significant and previously undetected changes that impact on the ecology of the system as a whole.

Dragon's paradise lost: palaeobiogeography, evolution and extinction of the largest-ever terrestrial lizards (Varanidae)

BACKGROUND

The largest living lizard species, Varanus komodoensis Ouwens 1912, is vulnerable to extinction, being restricted to a few isolated islands in eastern Indonesia, between Java and Australia, where it is the dominant terrestrial carnivore. Understanding how large-bodied varanids responded to past environmental change underpins long-term management of V. komodoensis populations.

METHODOLOGY/PRINCIPAL FINDINGS

We reconstruct the palaeobiogeography of Neogene giant varanids and identify a new (unnamed) species from the island of Timor. Our data reject the long-held perception that V. komodoensis became a giant because of insular evolution or as a specialist hunter of pygmy Stegodon. Phyletic giantism, coupled with a westward dispersal from mainland Australia, provides the most parsimonious explanation for the palaeodistribution of V. komodoensis and the newly identified species of giant varanid from Timor. Pliocene giant varanid fossils from Australia are morphologically referable to V. komodoensis suggesting an ultimate origin for V. komodoensis on mainland Australia (>3.8 million years ago). Varanus komodoensis body size has remained stable over the last 900,000 years (ka) on Flores, a time marked by major faunal turnovers, extinction of the island's megafauna, the arrival of early hominids by 880 ka, co-existence with Homo floresiensis, and the arrival of modern humans by 10 ka. Within the last 2000 years their populations have contracted severely.

CONCLUSIONS/SIGNIFICANCE

Giant varanids were once a ubiquitous part of Subcontinental Eurasian and Australasian faunas during the Neogene. Extinction played a pivotal role in the reduction of their ranges and diversity throughout the late Quaternary, leaving only V. komodoensis as an isolated long-term survivor. The events over the last two millennia now threaten its future survival.

The Liang Bua faunal remains: a 95k.yr. sequence from Flores, East Indonesia

Excavations at Liang Bua, a limestone cave on the island of Flores, East Indonesia, have yielded a well-dated archaeological and faunal sequence spanning the last 95k.yr., major climatic fluctuations, and two human species -H. floresiensis from 95 to 17k.yr.(1), and modern humans from 11k.yr. to the present. The faunal assemblage comprises well-preserved mammal, bird, reptile and mollusc remains, including examples of island gigantism in small mammals and the dwarfing of large taxa. Together with evidence from Early-Middle Pleistocene sites in the Soa Basin, it confirms the long-term isolation, impoverishment, and phylogenetic continuity of the Flores faunal community. The accumulation of Stegodon and Komodo dragon remains at the site in the Pleistocene is attributed to Homo floresiensis, while predatory birds, including an extinct species of owl, were largely responsible for the accumulation of the small vertebrates. The disappearance from the sequence of the two large-bodied, endemic mammals, Stegodon florensis insularis and Homo floresiensis, was associated with a volcanic eruption at 17 ka and precedes the earliest evidence for modern humans, who initiated use of mollusc and shell working, and began to introduce a range of exotic animals to the island. Faunal introductions during the Holocene included the Sulawesi warty pig (Sus celebensis) at about 7ka, followed by the Eurasian pig (Sus scrofa), Long-tailed macaque, Javanese porcupine, and Masked palm civet at about 4ka, and cattle, deer, and horse - possibly by the Portuguese within historic times. The Holocene sequence at the site also documents local faunal extinctions - a result of accelerating human population growth, habitat loss, and over-exploitation.

Is body shape of varanid lizards linked with retreat choice?

In our earlier analysis of Varanus body shape, size was a dominating factor with some qualitative phylogenetic patterns and grouping of species into ecological categories. With a phylogeny and an improved capacity to account for the effects of size, we have reanalysed our morphometric data for male Australian goannas (Varanus spp.) using an increased number of specimens and species to examine whether variations in body shape can be accounted for by retreat choice, as it can for Western Australian Ctenophorus dragon lizards. After accounting for body size in the current analysis, four ecotypes based on retreat choice (i.e. those that retreat to oblique crevices between large rocks or rock faces, those that retreat to burrows dug into the ground, those that retreat to spaces under rocks or in tree hollows, and those that retreat to trees but not tree hollows) accounted for much of the variation in body shape. There is a phylogenetic pattern to the ecotypes, but accounting for phylogenetic effects did not weaken the link between body shape and ecotype based on retreat choice. This suggests that there are large differences in body shape among ecotypes, and shape is relatively independent of phylogeny. The strong link between shape and choice of retreat site in Varanus spp. is consistent with that for Ctenophorus spp. We speculate on why there might be a strong link between retreat choice and body shape for both Varanus and Ctenophorus.

Sauria SINEs: Novel short interspersed retroposable elements that are widespread in reptile genomes

SINEs are short interspersed retrotransposable elements that invade new genomic sites. Their retrotransposition depends on reverse transcriptase and endonuclease activities encoded by partner LINEs (long interspersed elements). Recent genomic research has demonstrated that retroposons account for at least 40% of the human genome. Hitherto, more than 30 families of SINEs have been characterized in mammalian genomes, comprising approximately 4600 extant species; the distribution and extent of SINEs in reptilian genomes, however, are poorly documented. With more than 7400 species of lizards and snakes, Squamata constitutes the largest and most diverse group of living reptiles. We have discovered and characterized a novel SINE family, Sauria SINEs, whose members are widely distributed among genomes of lizards, snakes, and tuataras. Sauria SINEs comprise a 5' tRNA-related region, a tRNA-unrelated region, and a 3' tail region (containing short tandem repeats) derived from LINEs. We distinguished eight Sauria SINE subfamilies in genomes of four major squamate lineages and investigated their evolutionary relationships. Our data illustrate the overall efficacy of Sauria SINEs as novel retrotransposable markers for elucidation of squamate evolutionary history. We show that all Sauria SINEs share an identical 3' sequence with Bov-B LINEs and propose that they utilize the enzymatic machinery of Bov-B LINEs for their own retrotransposition. This finding, along with the ubiquity of Bov-B LINEs previously demonstrated in squamate genomes, suggests that these LINEs have been an active partner of Sauria SINEs since this SINE family was generated more than 200 million years ago.

A molecular phylogeny of the Australian monitor lizards (Squamata:Varanidae) inferred from mitochondrial DNA sequences

To date no complete phylogeny of all of the currently recognised Indo-Australian varanid species and subspecies has been published. This paper presents a comprehensive mitochondrial gene phylogeny of these lizards. A portion of the mitochondrial genome comprising part of the ND4 gene and three adjacent tRNA genes (hereafter referred to as ND4) was analysed alone and, for a subset of the taxa, combined with previously published mitochondrial data. Similar tree topologies were produced by both datasets although combining the data helped resolve some of the unresolved or weakly supported nodes in the ND4 analyses. The monophyly of the Indo-Australian group was strongly supported in all analyses. This group comprised three major lineages: the gouldii group, the Odatria group and the varius group. Mitochondrial ND4 nucleotide sequences were successfully amplified from all of the Indo-Australian monitor species and subspecies currently recognised and, as such, is the first comprehensive phylogenetic study of the Australian monitor lizards published. Analysis of the tempo of diversification and evolution of preferred habitat use identified six episodes of increased net speciation rate, with two closely adjacent episodes showing the highest rates of diversification and correlating with the appearance of all preferred habitat types. The comprehensive molecular phylogenetic framework will also be useful for the identification of varanid species and traded products derived from monitors and, as such, has important applications for wildlife management and conservation.

Mitochondrial genome of Pogona vitticepes (Reptilia; Agamidae): control region duplication and the origin of Australasian agamids

The complete mitochondrial DNA sequence for an Australian agamid Pogona vitticepes was determined. Twenty-two tRNA genes, two rRNA genes, thirteen protein-coding genes, and two control regions were identified in this mitochondrial genome. The second control region was inserted between NADH dehydrogenase subunits 5 and 6 genes. The duplication of the control region was found in all Australasian agamids examined and was not found in other Asian or African taxa. The two control regions had nearly identical sequences within species but they were divergent among species, suggesting their concerted sequence evolution. Phylogenetic analyses including divergence time estimation without assuming the molecular clock suggested that the duplication of the control region occurred on a lineage leading to the Australasian agamids 25-45 million years ago after their divergence from a Southeast Asian Physignathus cocincinus. Our finding thus supports the recent dispersal origin of Australasian agamids in connection with plate tectonic movement of Australia to the proximity of Southeast Asia.

A comparative test of adaptive hypotheses for sexual size dimorphism in lizards

It is commonly argued that sexual size dimorphism (SSD) in lizards has evolved in response to two primary, nonexclusive processes: (1) sexual selection for large male size, which confers an advantage in intrasexual mate competition (intrasexual selection hypothesis), and (2) natural selection for large female size, which confers a fecundity advantage (fecundity advantage hypothesis). However, outside of several well-studied lizard genera, the empirical support for these hypotheses has not been examined with appropriate phylogenetic control. We conducted a comparative phylogenetic analysis to test these hypotheses using literature data from 497 lizard populations representing 302 species and 18 families. As predicted by the intrasexual selection hypothesis, male aggression and territoriality are correlated with SSD, but evolutionary shifts in these categorical variables each explain less than 2% of the inferred evolutionary change in SSD. We found stronger correlations between SSD and continuous estimates of intrasexual selection such as male to female home range ratio and female home range size. These results are consistent with the criticism that categorical variables may obscure much of the actual variation in intrasexual selection intensity needed to explain patterns in SSD. In accordance with the fecundity advantage hypothesis, SSD is correlated with clutch size, reproductive frequency, and reproductive mode (but not fecundity slope, reduced major axis estimator of fecundity slope, length of reproductive season, or latitude). However, evolutionary shifts in clutch size explain less than 8% of the associated change in SSD, which also varies significantly in the absence of evolutionary shifts in reproductive frequency and mode. A multiple regression model retained territoriality and clutch size as significant predictors of SSD, but only 16% of the variation in SSD is explained using these variables. Intrasexual selection for large male size and fecundity selection for large female size have undoubtedly helped to shape patterns of SSD across lizards, but the comparative data at present provide only weak support for these hypotheses as general explanations for SSD in this group. Future work would benefit from the consideration of alternatives to these traditional evolutionary hypotheses, and the elucidation of proximate mechanisms influencing growth and SSD within populations.

Phylogeny of the lizard subfamily Lygosominae (Reptilia: Scincidae), with special reference to the origin of the new world taxa

Phylogenetic relationships of the three lygosomine skink genera occurring both in the Old World and the New World (Mabuya, Scincella and Sphenomorphus) were inferred from mitochondrial DNA sequence of 12S and 16S rRNA genes. Results strongly suggested the non-monophyly for any of these three genera. Within the Mabuya group, Asian members appear to have diverged first, leaving the Neotropical and the Afro-Malagasy Mabuya as sister groups. These relationships, together with the absence of extant or fossil representatives of the Mabuya group from North America, strongly suggest the trans-Atlantic dispersals of Mabuya from Africa to Neotropics. Our results also indicated a closer affinity of the New World Scincella with the New World Sphenomorphus than with the Old World Scincella. Such relationships suggest the trans-Beringian dispersal of the common ancestor from Asia and its subsequent divergence into the North American Scincella and the Neotropical Sphenomorphus.