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

A comparison of gene organisations and phylogenetic relationships of all 22 squamate species listed in South Korea using complete mitochondrial DNA

Studies using complete mitochondrial genome data have the potential to increase our understanding on gene organisations and evolutionary species relationships. In this study, we compared complete mitochondrial genomes between all 22 squamate species listed in South Korea. In addition, we constructed Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI) phylogenetic trees using 13 mitochondrial protein-coding genes. The mitochondrial genes for all six species in the suborder Sauria followed the same organisation as the sequenced Testudines (turtle) outgroup. In contrast, 16 snake species in the suborder Serpentes contained some gene organisational variations. For example, all snake species contained a second control region (CR2), while three species in the family Viperidae had a translocated tRNA-Pro gene region. In addition, the snake species, Elapheschrenckii, carried a tRNA-Pro pseudogene. We were also able to identify a translocation of a tRNA-Asn gene within the five tRNA (WANCY gene region) gene clusters for two true sea snake species in the subfamily Hydrophiinae. Our BI phylogenetic tree was also well fitted against currently known Korean squamate phylogenetic trees, where each family and genus unit forms monophyletic clades and the suborder Sauria is paraphyletic to the suborder Serpentes. Our results may form the basis for future northeast Asian squamate phylogenetic studies.

Toxicological impact of organophosphorus Chlorpyrifos 48%EC pesticide on hemocytes, biochemical disruption, and molecular changes in Biomphalaria alexandrina snails

Organophosphorus pesticides like Chlorpyrifos 48%EC were widely used to control agricultural pests. The present study aimed to evaluate the toxic effects of Chlorpyrifos 48%EC on B. alexandrina snails, the intermediate host of Schistosoma mansoni. After exposure of snails to serial concentrations to determine the LC₅₀, thirty snails for each sublethal concentration (LC₁₀ 2.1 and LC₂₅ 5.6 mg/l) in each group were exposed for 24 h followed by another 24 h for recovery. After recovery random samples were collected from hemolymph and tissue to measure the impacts on Phagocytic index, histological, biochemical, and molecular parameters. The current results showed a toxic effect of Chlorpyrifos 48%EC on adult B. alexandrina snails after 24 h of exposure at LC₅₀ 9.6 mg/l. After exposure to the sub-lethal concentrations of this pesticide, it decreased the total number of hemocytes and the percentage of small cells, while increased the percentage of hyalinocytes. The granulocyte percentage was increased after exposure to LC₁₀, while after LC₂₅, it was decreased compared to the control group. Also, the light microscopical examination showed that some granulocytes have plenty of granules, vacuoles and filopodia. Some hyalinocytes were contained shrinked nuclei, incomplete cell division and forming pseudopodia. Besides, the phagocytic index of hemocytes was significantly increased than control in all treated groups. Also, these sub-lethal concentrations increased MDA and SOD activities, while, tissue NO, GST and TAC contents were significantly decreased after exposure. Levels of Testosterone (T) and Estradiol (E) were increased significantly after exposure compared with control group. The present results showed that the concentration of DNA and RNA was highly decreased after exposure to LC₁₀, ₂₅ than the control group. Therefore, B. alexandrina snails could be used as a bio monitor of the chemical pollution. Besides, this pesticide could reduce the transmission of schistosomiasis as it altered the biological system of these snails.

Dynamics of epigenetic modifiers and environmentally sensitive proteins in a reptile with temperature induced sex reversal

The mechanisms by which sex is determined, and how a sexual phenotype is stably maintained during adulthood, has been the focus of vigorous scientific inquiry. Resources common to the biomedical field (automated staining and imaging platforms) were leveraged to provide the first immunofluorescent data for a reptile species with temperature induced sex reversal. Two four-plex immunofluorescent panels were explored across three sex classes (sex reversed ZZf females, normal ZWf females, and normal ZZm males). One panel was stained for chromatin remodelling genes JARID2 and KDM6B, and methylation marks H3K27me3, and H3K4me3 (Jumonji Panel). The other CaRe panel stained for environmental response genes CIRBP and RelA, and H3K27me3 and H3K4me3. Our study characterised tissue specific expression and cellular localisation patterns of these proteins and histone marks, providing new insights to the molecular characteristics of adult gonads in a dragon lizard Pogona vitticeps. The confirmation that mammalian antibodies cross react in P. vitticeps paves the way for experiments that can take advantage of this new immunohistochemical resource to gain a new understanding of the role of these proteins during embryonic development, and most importantly for P. vitticeps, the molecular underpinnings of sex reversal.

The complete mitochondrial genome of the parthenogenetic Caucasian rock lizard Darevskia unisexualis (Squamata: lacertidae) contains long tandem repeat formed by 59 bp monomer

The first complete mitochondrial genome sequence of parthenogenetic Caucasian rock lizard Darevskia unisexualis (Lacertidae family) is determined by hybrid assembly with Illumina HiSeq and PacBio RS II platforms. The circular 21.4 kbp mitogenome contains 13 protein-coding genes, 12S and 16S rRNA genes, 20 tRNAs, two pseudogenized tRNAs, and one long tandem repeats with 4.1 kbp length formed by 59 bp monomer repeated x70.6 times located before control region. This finding represents a new example of mitogenome variation in lizards of hybrid origin, providing fundamental data for following study of a unique hybridization system formed by parthenogenetic and bisexual species in the mountain steppe of central Armenia.

The first iguanian lizard from the Mesozoic of Africa

The fossil record shows that iguanian lizards were widely distributed during the Late Cretaceous. However, the biogeographic history and early evolution of one of its most diverse and peculiar clades (acrodontans) remain poorly known. Here, we present the first Mesozoic acrodontan from Africa, which also represents the oldest iguanian lizard from that continent. The new taxon comes from the Kem Kem Beds in Morocco (Cenomanian, Late Cretaceous) and is based on a partial lower jaw. The new taxon presents a number of features that are found only among acrodontan lizards and shares greatest similarities with uromastycines, specifically. In a combined evidence phylogenetic dataset comprehensive of all major acrodontan lineages using multiple tree inference methods (traditional and implied weighting maximum-parsimony, and Bayesian inference), we found support for the placement of the new species within uromastycines, along with Gueragama sulamericana (Late Cretaceous of Brazil). The new fossil supports the previously hypothesized widespread geographical distribution of acrodontans in Gondwana during the Mesozoic. Additionally, it provides the first fossil evidence of uromastycines in the Cretaceous, and the ancestry of acrodontan iguanians in Africa.

Gene rearrangements in gekkonid mitochondrial genomes with shuffling, loss, and reassignment of tRNA genes

BACKGROUND

Vertebrate mitochondrial genomes (mitogenomes) are 16-18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

RESULTS

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNAArg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNAArg gene appears to have been converted to a tRNAGln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNALeu genes, as well as a translocation of the tRNAAla gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNAGlu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNAGlu gene adjacent to the 5' end of the major noncoding region.

CONCLUSIONS

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

A novel model of double replications and random loss accounts for rearrangements in the Mitogenome of Samariscus latus (Teleostei: Pleuronectiformes)

BACKGROUND

Although more than one thousand complete mitochondrial DNA (mtDNA) sequences have been determined in teleostean fishes, only a few gene rearrangements have been observed, and genome-scale rearrangements are even rarer. However, flatfishes (Pleuronectiformes) have been identified as having diverse types of mitochondrial gene rearrangements. It has been reported that tongue soles and the blue flounder mitogenomes exhibit different types of large-scale gene rearrangements.

RESULTS

In the present study, the complete mitochondrial genome of another flatfish, Samariscus latus, was sequenced, and genome-scale rearrangements were observed. The genomic features of this flounder are different from those of any other studied vertebrates, including flatfish species too. The mitogenome of S. latus is characterized by the duplication and translocation of the control region (CR). The genes located between the two CRs are divided into two clusters in which their relative orders are maintained.

CONCLUSIONS

We propose a "Double Replications and Random Loss" model to explain the rearrangement events in S. latus mitogenome. This model consists of the following steps. First, the CR was duplicated and translocated. Subsequently, double replications of the mitogenome were successively initiated from the two CRs, leading to the duplication of the genes between the two CRs. Finally, one of each pair of duplicated genes was lost in a random event.

Forearm posture and mobility in quadrupedal dinosaurs

Quadrupedality evolved four independent times in dinosaurs; however, the constraints associated with these transitions in limb anatomy and function remain poorly understood, in particular the evolution of forearm posture and rotational ability (i.e., active pronation and supination). Results of previous qualitative studies are inconsistent, likely due to an inability to quantitatively assess the likelihood of their conclusions. We attempt to quantify antebrachial posture and mobility using the radius bone because its morphology is distinct between extant sprawled taxa with a limited active pronation ability and parasagittal taxa that have an enhanced ability to actively pronate the manus. We used a sliding semi-landmark, outline-based geometric morphometric approach of the proximal radial head and a measurement of the angle of curvature of the radius in a sample of 189 mammals, 49 dinosaurs, 35 squamates, 16 birds, and 5 crocodilians. Our results of radial head morphology showed that quadrupedal ceratopsians, bipedal non-hadrosaurid ornithopods, and theropods had limited pronation/supination ability, and sauropodomorphs have unique radial head morphology that likely allowed limited rotational ability. However, the curvature of the radius showed that no dinosaurian clade had the ability to cross the radius about the ulna, suggesting parallel antebrachial elements for all quadrupedal dinosaurs. We conclude that the bipedal origins of all quadrupedal dinosaur clades could have allowed for greater disparity in forelimb posture than previously appreciated, and future studies on dinosaur posture should not limit their classifications to the overly simplistic extant dichotomy.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Mosaic gene conversion after a tandem duplication of mtDNA sequence in Diomedeidae (albatrosses)

Although the tandem duplication of mitochondrial (mt) sequences, especially those of the control region (CR), has been detected in metazoan species, few studies have focused on the features of the duplicated sequence itself, such as the gene conversion rate, distribution patterns of the variation, and relative rates of evolution between the copies. To investigate the features of duplicated mt sequences, we partially sequenced the mt genome of 16 Phoebastria albatrosses belonging to three species (P. albatrus, P. nigripes, and P. immutabilis). More than 2,300 base pairs of tandemly-duplicated sequence were shared by all three species. The observed gene arrangement was shared in the three Phoebastria albatrosses and suggests that the duplication event occurred in the common ancestor of the three species. Most of the copies in each individual were identical or nearly identical, and were maintained through frequent gene conversions. By contrast, portions of CR domains I and III had different phylogenetic signals, suggesting that gene conversion had not occurred in those sections after the speciation of the three species. Several lines of data, including the heterogeneity of the rate of molecular evolution, nucleotide differences, and putative secondary structures, suggests that the two sequences in CR domain I are maintained through selection; however, additional studies into the mechanisms of gene conversion and mtDNA synthesis are required to confirm this hypothesis.

Mitochondrial genomes of acrodont lizards: timing of gene rearrangements and phylogenetic and biogeographic implications

Background

Acrodonta consists of Agamidae and Chamaeleonidae that have the characteristic acrodont dentition. These two families and Iguanidae sensu lato are members of infraorder Iguania. Phylogenetic relationships and historical biogeography of iguanian lizards still remain to be elucidated in spite of a number of morphological and molecular studies. This issue was addressed by sequencing complete mitochondrial genomes from 10 species that represent major lineages of acrodont lizards. This study also provided a good opportunity to compare molecular evolutionary modes of mitogenomes among different iguanian lineages.

Results

Acrodontan mitogenomes were found to be less conservative than iguanid counterparts with respect to gene arrangement features and rates of sequence evolution. Phylogenetic relationships were constructed with the mitogenomic sequence data and timing of gene rearrangements was inferred on it. The result suggested highly lineage-specific occurrence of several gene rearrangements, except for the translocation of the tRNA^Pro gene from the 5' to 3' side of the control region, which likely occurred independently in both agamine and chamaeleonid lineages. Phylogenetic analyses strongly suggested the monophyly of Agamidae in relation to Chamaeleonidae and the non-monophyly of traditional genus Chamaeleo within Chamaeleonidae. Uromastyx and Brookesia were suggested to be the earliest shoot-off of Agamidae and Chamaeleonidae, respectively. Together with the results of relaxed-clock dating analyses, our molecular phylogeny was used to infer the origin of Acrodonta and historical biogeography of its descendant lineages. Our molecular data favored Gondwanan origin of Acrodonta, vicariant divergence of Agamidae and Chamaeleonidae in the drifting India-Madagascar landmass, and migration of the Agamidae to Eurasia with the Indian subcontinent, although Laurasian origin of Acrodonta was not strictly ruled out.

Conclusions

We detected distinct modes of mitogenomic evolution among iguanian families. Agamidae was highlighted in including a number of lineage-specific mitochondrial gene rearrangements. The mitogenomic data provided a certain level of resolution in reconstructing acrodontan phylogeny, although there still remain ambiguous relationships. Our biogeographic implications shed a light on the previous hypothesis of Gondwanan origin of Acrodonta by adding some new evidence and concreteness.

The structure of the Australian and South American marsupial mitochondrial control region

BACKGROUND AND AIMS

The mitochondrial control region (CR) was studied across five marsupialian orders, in order to give a detailed overview of its features.

RESULTS

The CR is organised into three domains similar to the CR of placental mammals. However, the conservation of different features among the marsupial orders is in general more strict. In the first domain, two conserved blocks extended termination-associated sequences (ETAS 1 and ETAS 2) are present in all marsupial orders. In the third domain, the three conserved sequence blocks (CSB 1, CSB 2 and CSB 3) are present and complete, with CSB 1 being duplicated in four of five marsupial orders.

CONCLUSIONS

The nucleotide frequency and secondary structures of the repeats were typical for marsupial species. The repeats are generally AT-rich except in Dasyuridae and Paucituberculata, which show a significant increase in GC content.

Tandem duplication of mitochondrial DNA in the black-faced spoonbill, Platalea minor

Mitochondrial (mt) heteroplasmy in the control region (CR) of the black-faced spoonbill was investigated using LA-PCR. To avoid amplification of transpositioned nuclear genome fragment from mtDNA (numt), PCR product of the almost-complete mitochondrial genome was amplified using primers designed to anneal on the COIII gene. Then nested LA-PCR product was amplified between the cyt b and 12S rRNA genes using the almost-complete mitochondrial genome PCR product as a template. Nucleotide sequencing revealed tandem duplication composed of two units. The first contains cyt b-1, tRNA(Thr)-1, tRNA(Pro)-1, ND6-1, tRNA(Glu)-1 and CR1, and the second consists of cyt b-2, tRNA(Thr)-2, tRNA(Pro)-2, ND6-2, tRNA(Glu)-2 and CR2, followed by tRNA(Phe) and 12S rRNA. The duplicated cyt b-2 sequence coincided with 499 bp at the 3' end of cyt b-1. With the exception of the CR, the other genes in the duplicated sequence were identical to the original corresponding gene. Even though both CR1 and CR2 contain functional blocks, such as a poly-C site, a goose hairpin and a TAS structure in Domain I, the 3' end of CR1 was followed by a 112 bp sequence (non-coding region) that was not found in CR2 or in sequence homology analysis of similar genes. Meanwhile, CR2 ended in a complicated repeat sequence. The 5' franking region in the Domain I (Region A) and the 3' franking region in the Domain I (Region B) of the two CRs evolve in quite different manners: Region A was highly variable between CR1 and CR2 in the same individuals, while Region B was almost identical between them, which indicates concerted evolution.

Complete mitochondrial genome of the frillneck lizard (Chlamydosaurus kingii, Reptilia; Agamidae), another squamate with two control regions

Using PCR, the complete mitochondrial genome was sequenced in three frillneck lizards (Chlamydosaurus kingii). The mitochondria spanned over 16,761bp. As in other vertebrates, two rRNA genes, 22 tRNA genes and 13 protein coding genes were identified. However, similar to some other squamate reptiles, two control regions (CRI and CRII) were identified, spanning 801 and 812 bp, respectively. Our results were compared with another Australian member of the family Agamidae, the bearded dragon (Pogana vitticeps). The overall base composition of the light-strand sequence largely mirrored that observed in P vitticeps. Furthermore, similar to P. vitticeps, we observed an insertion 801 bp long between the ND5 and ND6 genes. However, in contrast to P vitticeps we did not observe a conserved sequence block III region. Based on a comparison among the three frillneck lizards, we also present data on the proportion of variable sites within the major mitochondrial regions.

Population genetic structure, gene flow and sex-biased dispersal in frillneck lizards (Chlamydosaurus kingii)

By using both mitochondrial and nuclear multiloci markers, we explored population genetic structure, gene flow and sex-specific dispersal of frillneck lizards (Chlamydosaurus kingii) sampled at three locations, separated by 10 to 50 km, in a homogenous savannah woodland in tropical Australia. Apart from a recombinant lizard, the mitochondrial analyses revealed two nonoverlapping haplotypes/populations, while the nuclear markers showed that the frillneck lizards represented three separate clusters/populations. Due to the small population size of the mtDNA, fixation may occur via founder effects and/or drift. We therefore suggest that either of these two processes, or a combination of the two, are the most likely causes of the discordant results obtained from the mitochondrial and the nuclear markers. In contrast to the nonoverlapping mitochondrial haplotypes, in 12 out of 74 lizards, mixed nuclear genotypes were observed, hence revealing a limited nuclear gene flow. Although gene flow should ultimately result in a blending of the populations, we propose that the distinct nuclear population structure is maintained by frequent fires resulting in local bottlenecks, and concomitant spatial separation of the frillneck lizard populations. Limited mark-recapture data and the difference in distribution of the mitochondrial and nuclear markers suggest that the mixed nuclear genotypes were caused by juvenile male-biased dispersal.

The molecular evolutionary tree of lizards, snakes, and amphisbaenians

Squamate reptiles (lizards, snakes, amphisbaenians) number approximately 8200 living species and are a major component of the world's terrestrial vertebrate diversity. Recent molecular phylogenies based on protein-coding nuclear genes have challenged the classical, morphology-based concept of squamate relationships, requiring new classifications, and drawing new evolutionary and biogeographic hypotheses. Even the key and long-held concept of a dichotomy between iguanians (approximately 1470 sp.) and scleroglossans (all other squamates) has been refuted because molecular trees place iguanians in a highly nested position. Together with snakes and anguimorphs, iguanians form a clade--Toxicofera--characterized by the presence of toxin secreting oral glands and demonstrating a single early origin of venom in squamates. Consequently, neither the varanid lizards nor burrowing lineages such as amphisbaenians or dibamid lizards are the closest relative of snakes. The squamate timetree shows that most major groups diversified in the Jurassic and Cretaceous, 200-66 million years (Myr) ago. In contrast, five of the six families of amphisbaenians arose during the early Cenozoic, approximately 60-40 Myr ago, and oceanic dispersal on floating islands apparently played a significant role in their distribution on both sides of the Atlantic Ocean. Among snakes, molecular data support the basic division between the small fossorial scolecophidians (approximately 370 sp.) and the alethinophidians (all other snakes, approximately 2700 sp.). They show that the alethinophidians were primitively macrostomatan and that this condition was secondarily lost by burrowing lineages. The diversification of alethinophidians resulted from a mid-Cretaceous vicariant event, the separation of South America from Africa, giving rise to Amerophidia (aniliids and tropidophiids) and Afrophidia (all other alethinophidians). Finally, molecular phylogenies have made it possible to draw a detailed evolutionary history of venom among advanced snakes (Caenophidia), a key functional innovation underlying their radiation (approximately 2500 sp.).

Mitochondrial DNA recombination in a free-ranging Australian lizard

Mitochondrial DNA (mtDNA) is the traditional workhorse for reconstructing evolutionary events. The frequent use of mtDNA in such analyses derives from the apparent simplicity of its inheritance: maternal and lacking bi-parental recombination. However, in hybrid zones, the reproductive barriers are often not completely developed, resulting in the breakdown of male mitochondrial elimination mechanisms, leading to leakage of paternal mitochondria and transient heteroplasmy, resulting in an increased possibility of recombination. Despite the widespread occurrence of heteroplasmy and the presence of the molecular machinery necessary for recombination, we know of no documented example of recombination of mtDNA in any terrestrial wild vertebrate population. By sequencing the entire mitochondrial genome (16761bp), we present evidence for mitochondrial recombination in the hybrid zone of two mitochondrial haplotypes in the Australian frillneck lizard (Chlamydosaurus kingii).

Effect of taxon sampling on recovering the phylogeny of squamate reptiles based on complete mitochondrial genome and nuclear gene sequence data

The complete nucleotide sequences of the mitochondrial (mt) genomes of three species of squamate lizards: Blanus cinereus (Amphisbaenidae), Anguis fragilis (Anguidae), and Tarentola mauritanica (Geckkonidae) were determined anew. The deduced amino acid sequences of all 13 mt protein-coding genes were combined into a single data set and phylogenetic relationships among main squamate lineages were analyzed under maximum likelihood (ML) and Bayesian Inference (BI). Within Squamata, the monophyly of Iguanidae, Anguimorpha, Amphisbaenia, Gekkota, Serpentes, and Acrodonta received high statistical support with both methods. It is particularly striking that this is the first molecular analysis that recovers the monophyly of Scincomorpha (including Scincidae, Xantusiidae, Cordylidae, and Lacertidae), although it is only supported in the Bayesian analysis, and it is sensitive to changes in the outgroup (see below). Phylogenetic relationships among the main squamate lineages could not be resolved with ML but received strong support with BI (above 95%). The newly reconstructed phylogeny of squamates does not support the Iguania-Scleroglossa split. Acrodonta and Serpentes form a clade, which is the sister group of the remaining squamate lineages. Within these, Gekkota were the first branching out, followed by Amphisbaenia, and a clade including Anguimorpha as sister group of Scincomorpha + Iguanidae. The recovered topology differed substantially from previously reported hypotheses on squamate relationships, and the relative effect of using different outgroups, genes, and taxon samplings were explored. The sister group relationship of Serpentes + Acrodonta, and their relative basal position within Squamata could be due to a long-branch attraction artifact. Phylogenetic relationships among Scincomorpha, Amphisbaenia, and Anguimorpha were found to be rather unresolved. Future improving of squamate phylogenetic relationships would rely on finding snake and acrodont species with slower mt evolutionary rates, ensuring thorough taxon coverage of squamate diversity, and incorporating more nuclear genes with appropriate evolutionary rates.

Mitogenomic perspectives into iguanid phylogeny and biogeography: Gondwanan vicariance for the origin of Madagascan oplurines

Complete or nearly complete nucleotide sequences of mitochondrial genomes (mtDNAs) were determined from eight species which, together with previous mtDNA data for two other taxa, cover most subfamilies of Iguanidae sensu lato. These iguanid mtDNAs were found to be rather conservative with respect to gene arrangements and molecular evolutionary rates, which contrasts with mtDNAs of Acrodonta (Agamidae and Chamaeleonidae) in which several gene rearrangements and highly accelerated molecular evolutionary rates have been known. Phylogenetic analyses consistently suggested the earliest shoot-off of a Malagasy subfamily Oplurinae and an affinity of Polychrotinae and Tropidurinae sensu stricto. However, even with the ample molecular characters derived from complete mtDNA sequences, phylogenetic relationships between iguanid subfamilies were poorly resolved in general, presumably due to the rapid ancient cladogenesis. Divergence time estimation without assuming the molecular clock suggested the Late Triassic/Early Jurassic divergence of Iguanidae from acrodonts and the Middle/Late Jurassic divergence of Oplurinae from the other iguanids. Together with geological and paleontological evidence, these results led us to propose Gondwanan vicariance for the origin of Malagasy oplurines without invoking a land bridge connection between South America/Antarctica and drifting Madagascar/India.

Complete mitochondrial genomes of three neobatrachian anurans: a case study of divergence time estimation using different data and calibration settings

We sequenced the whole mitochondrial (mt) genomes of three neobatrachian species: Japanese tree frog Hyla japonica, Japanese common toad Bufo japonicus, and narrow-mouthed toad Microhyla okinavensis. The gene arrangements of these genomes diverged from that of basal anurans (suborder Archaeobatrachia), but are the same as that of the members of derived frogs (i.e., superfamily Hyloidae and Ranoidae in suborder Neobatrachia), suggesting the one-time occurrence of a gene rearrangement event in an ancestral lineage of derived anurans. Furthermore, several distinct repeat motifs including putative termination-associated sequences (TASs) and conserved sequence blocks (CSBs) were observed in the control regions (CRs) of B. japonicus and H. japonica, while no repeat motifs were found in that of M. okinavensis. Phylogenetic analyses using both nucleotide and amino acid data of mt genes support monophyly of neobatrachians. The estimated divergence time based on amino acid data with multiple reference points suggests that the three living amphibian orders may have originated in the Carboniferous period, and that the divergences of anurans had occurred between the Permian and Tertiary periods. We also checked the influence of the data types and the settings of reference times on divergence time estimation. The resultant divergence times estimated from several datasets and reference time settings suggest that the substitution saturation of nucleotide data may lead to overestimated (i.e., older) branching times, especially for early divergent taxa. We also found a highly accelerated substitution rate in neobatrachian mt genes, and fast substitution possibly resulted in overestimation. To correct this erroneous estimation, it is efficient to apply several reference points among neobatrachians.

Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region

BACKGROUND

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

RESULTS

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

CONCLUSION

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.

Partitioned Bayesian analyses, dispersal-vicariance analysis, and the biogeography of Chinese toad-headed lizards (Agamidae: Phrynocephalus): a re-evaluation

The toad-headed lizards of genus Phrynocephalus are distributed from northwestern China to Turkey and are one of the major components of the central Asian desert fauna. To date, published morphological and molecular phylogenetic hypotheses of Phrynocephalus are only partially congruent, and the relationships within the genus are still far from clear. We re-analyzed published mitochondrial gene sequence data (12S, 16S, cyt b, ND4-tRNA(Leu)) by employing partition-specific modeling in a combined DNA analysis to clarify existing gaps in the phylogeny of Chinese Phrynocephalus. Using this phylogenetic framework, we inferred the genus' historical biogeography by using weighted ancestral-area analysis and dispersal-vicariance analysis in combination with a Bayesian relaxed molecular-clock approach and paleogeographical data. The partitioned Bayesian analyses support the monophyly of Phrynocephalus and its sister-group relationship with Laudakia. An earlier finding demonstrating the monophyly of the viviparous group is corroborated. However, our hypothesis of internal relationships of the oviparous group differs from a previous hypothesis as our results do not support monophyly of the oviparous taxa. Instead, the viviparous taxa form a clade with many oviparous taxa exclusive of P. helioscopus and P. mystaceus. Our results also suggest that: (1) P. putjatia is a valid species, comprising populations from Guide, Qinghai Province and Tianzhu, Gansu Province; (2) P. hongyuanensis is not a valid species, synonymized instead with P. vlangalii; (3) P. zetangensis is not a valid species and should be included in P. theobaldi; (4) the population occurring in Kuytun, Xinjiang Uygur Autonomous Region is recognized as P. guttatus instead of P. versicolor; and (5) the Lanzhou population of P. frontalis is part of P. przewalskii. Congruent with previous hypotheses, the uplift of the Tibetan Plateau played a fundamental role in the diversification of Phrynocephalus. An evolutionary scenario combining aspects of vicariance and dispersal is necessary to explain the distribution of Phrynocephalus. Bayesian divergence-time estimation suggests that Phrynocephalus originated at the Middle-Late Miocene boundary (15.16-10.4 Ma), and diversified from Late Miocene to Pleistocene from a center of origin in Central Asia, Tarim Basin, and Junggar Basin temperate desert, followed by several rapid speciation events in a relatively short time. The proposed biogeographic scenarios also indicate that the Tarim Basin desert may be the secondary diversification center, followed by Junggar Basin temperate desert and Alashan Plateau temperate desert. In the viviparous group, the allopatric speciation of P. theobaldi and P. vlangalii may have been caused by the uplifting of Tanggula Mountain Ranges. In addition, the results of this study make an important contribution to understanding the uplift of the Tibetan Plateau and Tian Shan Mountains and the biogeography of the entire region.

The complete mitochondrial genome of the Green Lizard Lacerta viridis viridis (Reptilia: Lacertidae) and its phylogenetic position within squamate reptiles

For the first time the complete mitochondrial genome was sequenced for a member of Lacertidae. Lacerta viridis viridis was sequenced in order to compare the phylogenetic relationships of this family to other reptilian lineages. Using the long-polymerase chain reaction (long PCR) we characterized a mitochondrial genome, 17,156 bp long showing a typical vertebrate pattern with 13 protein coding genes, 22 transfer RNAs (tRNA), two ribosomal RNAs (rRNA) and one major noncoding region. The noncoding region of L. v. viridis was characterized by a conspicuous 35 bp tandem repeat at its 5' terminus. A phylogenetic study including all currently available squamate mitochondrial sequences demonstrates the position of Lacertidae within a monophyletic squamate group. We obtained a narrow relationship of Lacertidae to Scincidae, Iguanidae, Varanidae, Anguidae, and Cordylidae. Although, the internal relationships within this group yielded only a weak resolution and low bootstrap support, the revealed relationships were more congruent with morphological studies than with recent molecular analyses.

Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent Anomalurus

Background

Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious.

Results

We have sequenced and performed extensive evolutionary analyses on the mitochondrial genome of the scaly-tailed flying squirrel Anomalurus sp., an enigmatic rodent whose phylogenetic affinities have been obscure and extensively debated. Our phylogenetic analyses of the coding regions of available complete mitochondrial genome sequences from Euarchontoglires suggest that Anomalurus is a sister taxon to the Hystricognathi, and that this clade represents the most basal divergence among sampled Rodentia. Bayesian dating methods incorporating a relaxed molecular clock provide divergence-time estimates which are consistently in agreement with the fossil record and which indicate a rapid radiation within Glires around 60 million years ago.

Conclusion

Taken together, the data presented provide a working hypothesis as to the phylogenetic placement of Anomalurus, underline the utility of mitochondrial sequences in the resolution of even relatively deep divergences and go some way to explaining the difficulty of conclusively resolving higher-level relationships within Glires with available data and methodologies.

Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.