Xiaomao Zeng

Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: A case study of salamanders (Order Caudata)

In the practice of molecular dating, substitution saturation will bias the results if not properly modeled. Date estimates based on commonly used mitochondrial DNA sequences likely suffer from this problem because of their high substitution rate. Nevertheless, the patterns and extent of such expected bias remain unknown for many major evolutionary lineages, which often differ in ages, available calibrations, and substitution rates of their mitochondrial genome. In this case study of salamanders, we used estimates based on multiple nuclear exons to assess the effects of saturation on dating divergences using mitochondrial genome sequences on a timescale of ~200-300 My. The results indicated that, due to saturation for older divergences and in the absence of younger effective calibration points, dates derived from the mitochondrial data were considerably overestimated and systematically biased toward the calibration point for the ingroup root. The overestimate might be as great as 3-10 times (about 20 My) older than actual divergence dates for recent splitting events and 40 My older for events that are more ancient. For deep divergences, dates estimated were strongly compressed together. Furthermore, excluding the third codon positions of protein-coding genes or only using the RNA genes or second codon positions did not considerably improve the performance. In the order Caudata, slowly evolving markers such as nuclear exons are preferred for dating a phylogeny covering a relatively wide time span. Dates estimated from these markers can be used as secondary calibrations for dating recent events based on rapidly evolving markers for which mitochondrial DNA sequences are attractive candidates due to their short coalescent time. In other groups, similar evaluation should be performed to facilitate the choice of markers for molecular dating and making inferences from the results.

A phylogeny of Chinese species in the genus Phrynocephalus (Agamidae) inferred from mitochondrial DNA sequences

We investigated the phylogenetic relationships among most Chinese species of lizards in the genus Phrynocephalus (118 individuals collected from 56 populations of 14 well-defined species and several unidentified specimens) using four mitochondrial gene fragments (12S rRNA, 16S rRNA, cytochrome b, and ND4-tRNA(LEU)). The partition-homogeneity tests indicated that the combined dataset was homogeneous, and maximum-parsimony (MP), neighbor-joining (NJ), maximum-likelihood (ML) and Bayesian (BI) analyses were performed on this combined dataset (49 haplotypes including outgroups for 2058bp in total). The maximum-parsimony analysis resulted in 24 equally parsimonious trees, and their strict consensus tree shows that there are two major clades representing the Chinese Phrynocephalus species: the viviparous group (Clade A) and the oviparous group (Clade B). The trees derived from Bayesian, ML, and NJ analyses were topologically identical to the MP analysis except for the position of P. mystaceus. All analyses left the nodes for the oviparous group, the most basal clade within the oviparous group, and P. mystaceus unresolved. The phylogenies further suggest that the monophyly of the viviparous species may have resulted from vicariance, while recent dispersal may have been important in generating the pattern of variation among the oviparous species.

COI is better than 16S rRNA for DNA barcoding Asiatic salamanders (Amphibia: Caudata: Hynobiidae)

The 5′ region of the mitochondrial DNA (mtDNA) gene cytochrome c oxidase I (COI) is the standard marker for DNA barcoding. However, because COI tends to be highly variable in amphibians, sequencing is often challenging. Consequently, another mtDNA gene, 16S rRNA gene, is often advocated for amphibian barcoding. Herein, we directly compare the usefulness of COI and 16S in discriminating species of hynobiid salamanders using 130 individuals. Species identification and classification of these animals, which are endemic to Asia, are often based on morphology only. Analysis of Kimura 2‐parameter genetic distances (K2P) documents the mean intraspecific variation for COI and 16S rRNA genes to be 1.4% and 0.3%, respectively. Whereas COI can always identify species, sometimes 16S cannot. Intra‐ and interspecific genetic divergences occasionally overlap in both markers, thus reducing the value of a barcoding gap to identify genera. Regardless, COI is the better DNA barcoding marker for hynobiids. In addition to the comparison of two potential markers, high levels of intraspecific divergence in COI (>5%) suggest that both Onychodactylus fischeri and Salamandrella keyserlingii might be composites of cryptic species.

How many species are in the genus Batrachuperus? A phylogeographical analysis of the stream salamanders (family Hynobiidae) from southwestern China

Phylogeographical analysis of DNA sequence data has been routinely used to test species boundaries using the monophyly criterion; however, a complementary criterion, reproductive isolation, is often ignored. We used a combination of phylogenetic and population genetic approaches to determine species boundaries among stream salamanders in the genus Batrachuperus. First, cytochrome b sequence data from 174 Batrachuperus individuals, sampled from 78 populations, were used to reconstruct historical relationships within the genus. Second, allozyme data for 14 presumptive nuclear loci, from 463 individuals sampled from 60 populations, were collected and analysed to assess population similarity or disparity, as well as potential reproductive isolation. The DNA sequence data grouped all populations into seven major monophyletic groups, and the allozyme data provided evidence for reproductive isolation among four of the seven groups, thereby supporting the species status of these groups. The allozyme data suggested that two of the other groups share the same gene pool, and therefore belong to a single species. Finally, the allozyme data revealed two reproductively isolated units within the seventh group, which we suggest represents a case of ‘budding speciation’ based on the DNA gene tree. In total, seven species of the genus Batrachuperus were defined, two of which were previously unknown. The phylogeographical analysis also revealed that vicariance events might have dominated the evolutionary history of this group, but the speciation events might precede the formation of the existing mountain topology. This study demonstrates the importance of including frequency data from multiple nuclear gene loci in determining species boundaries.

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

BackgroundAlthough mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.ResultsGene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and dN/dS ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.ConclusionsNo evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Coalescence patterns of endemic Tibetan species of stream salamanders (Hynobiidae: Batrachuperus)

Orogenesis of topographically diverse montane regions often drives complex evolutionary histories of species. The extensive biodiversity of the eastern edge of the Tibetan Plateau, which gradually decreases eastwardly, facilitates a comparison of historical patterns. We use coalescence methods to compare species of stream salamanders (Batrachuperus) that occur at high and low elevations. Coalescent simulations reveal that closely related species are likely to have been influenced by different drivers of diversification. Species living in the western high‐elevation region with its northsouth extending mountains appear to have experienced colonization via dispersal followed by isolation and divergence. In contrast, species on the eastern low‐elevation region, which has many discontinuous mountain ranges, appear to have experienced fragmentation, sometimes staged, of wide‐ranging ancestral populations. The two groups of species appear to have been affected differently by glaciation. High‐elevation species, which are more resistant to cooler temperatures, appear to have experienced population declines as recently as the last glaciation (0.016–0.032 Ma). In contrast, salamanders dwelling in the warmer and wetter habitats at low‐elevation environs appear to have been affected less by the relatively recent, milder glaciation, and more so by harsher, extensive glaciations (0.5–0.175 Ma). Thus, elevation, topography and cold tolerance appear to drive evolutionary patterns of diversification and demography even among closely related taxa. The comparison of multiple species in genealogical analyses can lead to an understanding of the evolutionary drivers.

Rediscovery of Protohynobius puxiongensis (Caudata: Hynobiidae) and its phylogenetic position based on complete mitochondrial genomes

The mysterious Asian hynobiid salamander, Protohynobius puxiongensis, was described based on a single specimen collected in 1965 and never found again since then. Because the specimen had an internasal bone, Pr. puxiongensis was thought to retain a primitive character lost by a common ancestor of all other hynobiid salamanders, and it was thus considered to be not only a new genus, but also a new subfamily. This conclusion bothered herpetologists for decades because it was based on only one specimen and one character without other living specimens being rediscovered. After years of field effort, we rediscovered living individuals of Pr. puxiongensis at its type locality. All characters observed in rediscovered specimens are identical to the original description of the holotype except the internasal bone, implying that the internasal bone observed in the holotype may be just an individual variation. To examine the phylogenetic position of Pr. puxiongensis, we sequenced complete mitochondrial genomes for this species, together with two Pseudohynobius species. By combining 18 published hynobiid mitochondrial genomes and our new sequences, we reconstructed a comprehensive phylogenetic relationship of Hynobiidae at the genus level. Our results indicate that Pr. puxiongensis is deeply nested within the hynobiid phylogeny. It is the sister group of the Pseudohynobius species, and the validity of subfamily Protohynobiinae is not supported.

Genetic Diversity of Eastern Batrachuperus (Caudata: Hynobiidae)

Abstract Using mitochondrial cytochrome b and 16s gene sequences, we examined the genetic diversity of eastern Batrachuperus distributed in China along the eastern slopes of the Tibet Plateau. The group showed a moderate pairwise divergence compared to other salamander groups. The highest cytochrome b divergence is 10.07%. A phylogenetic analysis of the mitochondrial haplotypes revealed several evolutionarily independent lineages, including the Hongya and Pengxian populations of Batrachuperus tibetanus that probably represent undescribed species. The genetic divergence of Batrachuperus is clearly much higher than currently perceived. The distribution patterns of the mitochondrial haplotypes of the northwestern populations indicated that the dispersal or divergence of Batrachuperus might not be confined to river valleys, although all species in the genus are fully aquatic.

A De Novo Case of Floating Chromosomal Polymorphisms by Translocation in Quasipaa boulengeri (Anura, Dicroglossidae)

Very few natural polymorphisms involving interchromosomal reciprocal translocations are known in amphibians even in vertebrates. In this study, thirty three populations, including 471 individuals of the spiny frog Quasipaa boulengeri, were karyotypically examined using Giemsa stain or FISH. Five different karyomorphs were observed. The observed heteromorphism was autosomal but not sex-related, as the same heteromorphic chromosomes were found both in males and females. Our results indicated that the variant karyotypes resulted from a mutual interchange occurring between chromosomes 1 and 6. The occurrence of a nearly whole-arm translocation between chromosome no. 1 and no. 6 gave rise to a high frequency of alternate segregation and probably resulted in the maintenance of the translocation polymorphisms in a few populations. The translocation polymorphism is explained by different frequencies of segregation modes of the translocation heterozygote during meiosis. Theoretically, nine karyomorphs should be investigated, however, four expected karyotypes were not found. The absent karyomorphs may result from recessive lethal mutations, position effects, duplications and deficiencies. The phylogenetic inference proved that all populations of Q. boulengeri grouped into a monophyletic clade. The mutual translocation likely evolved just once in this species and the dispersal of the one karyomorph (type IV) can explain the chromosomal variations among populations.

A reinvestigation of phylogeny and divergence times of Hynobiidae (Amphibia, Caudata) based on 29 nuclear genes.

Although several recent studies have investigated the major phylogenetic relationships within Hynobiidae, their evolutionary history remains partially resolved and the phylogenetic positions of some genera, particularly Pachyhynobius and Salamandrella are still disputed. Notably, previous studies relied primarily on mitochondrial DNA data and concatenated analyses; thus, a new investigation based on multiple nuclear genes and species-tree inference is needed. Here, we provide an in-depth phylogenetic analysis, based on 29 nuclear genes comprising 29,232bp of data from a comprehensive taxonomic sampling (24 hynobiids and 7 outgroups), using both concatenated and species-tree methods. Our results robustly resolved most genus-level relationships within Hynobiidae, including the placement of Salamandrella as the sister group to a clade containing Batrachuperus, Liua and Pseudohynobius, and the placement of Pachyhynobius as the sister group to a clade containing all hynobiids excluding Onychodactylus, Paradactylodon and Ranodon. Time estimates based on our data suggest that the major group of living hynobiids (excluding Onychodactylus) originated approximately 40Ma, ∼50% younger than estimates from mtDNA data (62.5Ma) but 10% older than estimates from three nuclear genes (36Ma). Our results highlight the benefits of using a large number of nuclear loci to infer both phylogeny and time for relatively old lineages.