Extreme mito-nuclear discordance in a peninsular lizard: the role of drift, selection, and climate

Nuclear and mitochondrial genomes coexist within cells but are subject to different tempos and modes of evolution. Evolutionary forces such as drift, mutation, selection, and migration are expected to play fundamental roles in the origin and maintenance of diverged populations; however, divergence may lag between genomes subject to different modes of inheritance and functional specialization. Herein, we explore whole mitochondrial genome data and thousands of nuclear single nucleotide polymorphisms to evidence extreme mito-nuclear discordance in the small black-tailed brush lizard, Urosaurus nigricaudus, of the Peninsula of Baja California, Mexico and southern California, USA, and discuss potential drivers. Results show three deeply divergent mitochondrial lineages dating back to the later Miocene (ca. 5.5 Ma) and Pliocene (ca. 2.8 Ma) that likely followed geographic isolation due to trans-peninsular seaways. This contrasts with very low levels of genetic differentiation in nuclear loci (F_ST < 0.028) between mtDNA lineages. Analyses of protein-coding genes reveal substantial fixed variation between mitochondrial lineages, of which a significant portion comes from non-synonymous mutations. A mixture of drift and selection is likely responsible for the rise of these mtDNA groups, albeit with little evidence of marked differences in climatic niche space between them. Finally, future investigations can look further into the role that mito-nuclear incompatibilities and mating systems play in explaining contrasting nuclear gene flow.

A new species of the genus Microhyla Tschudi, 1838 (Amphibia: Anura: Microhylidae) from Tay Nguyen Plateau, Central Vietnam

We describe a new species of Microhyla from Tram Lap forest, Gia Lai Province, Central Vietnam based on morphological, molecular, and acoustic data. The new species resembles M. butleri morphologically, but differs from all congeners by a combination of the following morphological attributes: (1) medium-sized adult snout-vent length 25.2-27.0 mm in 15 males and 30.5 mm in a single female; (2) body habitus moderately stocky; (3) head flat, snout rounded, slightly prominent in ventral profile; (4) dorsum and flanks slightly shagreened with evenly scattered tiny tubercles, ventral skin smooth; (5) first finger well developed, more than one-half the length of the second finger; (6) tips of three outer fingers slightly enlarged, forming weak disks and tips of all toes distinctly dilated into wide disks with narrow peripheral grooves; (7) finger and toe disks with dorsal median longitudinal grooves; (8) three palmar tubercles and two metatarsal tubercles; (9) tibiotarsal articulation of adpressed limb reaching slightly beyond the orbit; (10) webbing formula: I 1¾-2 II 1½-2¾ III 2-31/3 IV 3¼-1½ V; (11) in life, chin and throat yellowish to bright-orange with tiny dark brown speckling laterally; and (12) a call consisting of 15-26 pulses with a dominant frequency of 1.8-2.2 kHz (recorded at 18.5ºC). We also provide a preliminary genealogy of Microhyla based on analysis of a 2644 bp fragment of 12S-16S rRNA mitochondrial DNA. Based on the examed data, the new species and M. butleri are sister-species (genetic p-distance: 9.0%) and it can be distinguished from M. butleri by its morphology (size, webbing on toes, color) and advertisement call. Interspecific genetic p-distances between the new species and its congeners vary from 9.0% to 14.8%. Microhyla aurantiventris sp. nov. occurs in evergreen montane tropical forests at elevations around 1200 m a.s.l. and is known only from the type locality. The new species appears to be threatened due to intensive logging and agriculture plantation.

Convergent Evolution of Mitochondrial Genes in Deep-Sea Fishes

Deep seas have extremely harsh conditions including high hydrostatic pressure, total darkness, cold, and little food and oxygen. The adaptations of fishes to deep-sea environment apparently have occurred independently many times. The genetic basis of adaptation for obtaining their energy remains unknown. Mitochondria play a central role in aerobic respiration. Analyses of the available 2,161 complete mitochondrial genomes of 1,042 fishes, including 115 deep-sea species, detect signals of positive selection in mitochondrial genes in nine branches of deep-sea fishes. Aerobic metabolism yields much more energy per unit of source material than anaerobic metabolism. The adaptive evolution of the mtDNA may reflect that aerobic metabolism plays a more important role than anaerobic metabolism in deep-sea fishes, whose energy sources (food) are extremely limited. This strategy maximizes the usage of energy sources. Eleven mitochondrial genes have convergent/parallel amino acid changes between branches of deep-sea fishes. Thus, these amino acid sites may be functionally important in the acquisition of energy, and reflect convergent evolution during their independent invasion of the harsh deep-sea ecological niche.

How many species of Apodemus and Rattus occur in China? A survey based on mitochondrial cyt b and morphological analyses

Apodemus (mice) and Rattus (rats) are the top rodent reservoirs for zoonoses in China, yet little is known about their diversity. We reexamined the alpha diversity of these two genera based on a new collection of specimens from China and their cyt b sequences in GenBank. We also tested whether species could be identified using external and craniodental measurements exclusively. Measurements from 147 specimens of Apodemus and 236 specimens of Rattus were used for morphological comparisons. We analysed 74 cyt b sequences of Apodemus and 100 cyt b sequences of Rattus to facilitate phylogenetic estimations. Results demonstrated that nine species of Apodemus and seven species of Rattus, plus a new subspecies of Rattus nitidus, are distributed in China. Principal component analysis using external and craniodental measurements revealed that measurements alone could not separate the recognized species. The occurrence of Rattus pyctoris in China remains uncertain.

The adaptation of codon usage of +ssRNA viruses to their hosts

Viruses depend on their host's cellular structure to survive. Most of them do not have tRNAs, their translation relies on hosts' tRNA pools. Over the course of evolution, viruses needed to optimally exploit cellular processes of their host. Thus, codon usage of a virus should coevolve with its host to efficiently and rapidly replicate. Some viruses can invade a broad spectrum of hosts (BSTVs), while others can invade a narrow spectrum only (NSTVs). Consequently, we test the hypothesis that similarity of codon usage preference and the degree of matching between BSTVs and their hosts will be lower than that of NSTVs, which only need to coevolve with few hosts. We compare the patterns of codon usage in 255 virus genomes to test this hypothesis. Our results show that NSTVs have a higher degree of matching to their hosts' tRNA pools than BSTVs. Further, analysis of the effective number of codons (ENC) infers that codon usage bias of NSTVs is relatively stronger than that of BSTVs. Thus, codon usage of NSTVs tends to better match their host than that of BSTVs. This supports the hypothesis that viruses adapt to the expression system of their host(s).

Evolutionary dynamics of avian influenza A H7N9 virus across five waves in mainland China, 2013-2017

Since its emergence in March 2013, novel avian influenza A H7N9 virus has triggered five epidemics of human infections in China. This raises concerns about the pandemic threat of this quickly evolving H7N9 subtype for humans. In this study, we evaluated all available genomes for H7N9 and H9N2 influenza A viruses. Our assessment discovered that H7N9 of the 1st wave had the lowest nucleotide diversity, which then experienced substantial and rapid population expansion from a small founder population. From the 2nd wave, their nucleotide diversity increased quickly, indicating that H7N9 viruses had acquired larger populations and mutations after their initial emergence in 2013. After the phylogeographic divergence in the 2nd wave, although the HA and NA genes from different regions differed, compared to previous epidemics, the evolving H7N9 viruses in the 5th wave lost most of their previous clades. The highly pathogenic avian influenza (HPAI) H7N9 viruses in the 5th wave clustered together, and clustered close to the low pathogenic avian influenza (LPAI) virus isolated from the Pearl River Delta in the 3rd and 4th waves. This result supports the origin of HPAI H7N9 viruses was in the Pearl River Delta. In the 5th wave, although both HPAI and LPAI H7N9 viruses were isolated from the Pearl River Delta, their HA and NA genes were phylogenetically distinct.

Landscape limits gene flow and drives population structure in Agassiz's desert tortoise (Gopherus agassizii)

Distance, environmental heterogeneity and local adaptation can strongly influence population structure and connectivity. Understanding how these factors shape the genomic landscape of threatened species is a major goal in conservation genomics and wildlife management. Herein, we use thousands (6,859) of single nucleotide polymorphism markers and spatial data from hundreds of individuals (n = 646) to re-evaluate the population structure of Agassiz's desert tortoise (Gopherus agassizii). Analyses resolve from 4 to 8 spatially well-defined clusters across the range. Western, central, and southern populations within the Western Mojave recovery unit are consistent throughout, while analyses sometimes merge other recovery units depending on the level of clustering. Causal modeling consistently associates genetic connectivity with least-cost distance, based on multiple landscape features associated with tortoise habitat, better than geographic distance. Some features include elevation, soil depth, rock volume, precipitation, and vegetation coverage, suggesting that physical, climatic, and biotic landscape features have played a strong evolutionary role restricting gene flow between populations. Further, 12 highly differentiated outlier loci have associated functions that may be involved with neurogenesis, wound healing, lipid metabolism, and possibly vitellogenesis. Together, these findings have important implications for recovery programs, such as translocations, population augmentation, reproduction in captivity and the identification of ecologically important genes, opening new venues for conservation genomics in desert tortoises.

CRISPR/Cas9 Application for Gene Copy Fate Survey of Polyploid Vertebrates

Polyploidization occurs widely in eukaryotes, and especially in plants. Polyploid plants and some fishes have been commercialized. Typically, severe genomic perturbations immediately follow polyploidization and little is known about how polyploid offspring survives the genetic and epigenetic changes. Investigations into this require the identification of genes related to polyploidization and the discrimination of dosage-balance from paternal and maternal copies, and regardless of the mechanism being either autopolyploidization or allopolyploidization. New approaches and technologies may discern the mosaic of novel gene functions gained through the recombination of paternal and maternal genes in allopolyploidization. Modifications of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) with CRISPR-associated system (Cas) protein 9 (CRISPR/Cas9) have been employed in studies of polyploidization of plants. However, the approach has seldom been applied to polyploidization in vertebrates. Herein, we use CRISPR/Cas9 to trace gene-fate in tetraploid goldfish, and specifically to identify the functional differentiation of two divergent copies of fgf20a, which are expressed differently throughout embryonic development. We expect this gene editing system will be applicable to studies of polyploids and the genetic improvement of polyploid livestock.

An Evolutionary Genomic Perspective on the Breeding of Dwarf Chickens

The evolutionary history for dwarfism in chickens remains an enigma. Herein, we explore the evolution of the Serama, the smallest breed of chicken. Leveraging comparative population genomics, analyses identify several genes that are potentially associated with the growth and development of bones and muscles. These genes, and in particular both POU1F1 and IGF1, are under strong positive selection. Three allopatric dwarf bantams (Serama, Yuanbao, and Daweishan) with different breeding-histories, form distinct clusters and exhibit unique population structures. Parallel genetic mechanisms underlay their variation in body size. These findings provide insights into the multiple and complex pathways, depending on genomic variation, that chicken can take in response to aviculture selection for dwarfism.

A new skink of the genus Sphenomorphus Fitzinger, 1843 (Squamata: Scincidae) from Hon Ba Nature Reserve, southern Vietnam

A new forest skink of the genus Sphenomorphus Fitzinger, 1843 is described from Khanh Hoa Province, southern Vietnam based on morphological characters of four specimens and a fragment of 653 nucleotides of the gene COI. Sphenomorphus yersini sp. nov. is characterized by the following morphological characters: medium size in adults (snout-vent length up to 55 mm); tail length/snout-vent length ratio 1.81; toes reach to fingers when limbs adpressed; midbody scale rows 32-34, smooth; paravertebral scales 61-69; ventral scale rows 58-67; subcaudal scales 112; supraoculars four, rarely five; prefrontals in broad contact with one another; loreal scales two; tympanum deeply sunk; smooth lamellae beneath finger and toe IV 10-12 and 18-20 respectively; a pair of enlarged precloacal scales; hemipenis deeply forked and asymmetrical with two differently sized smooth lobes. The new species differs from its most similar congener, Sphenomorphus buenloicus Darevsky Nguyen, 1983, by 16.4-16.7% uncorrected p-distance in COI sequences.

Reassessment of the taxonomic status of Craseomys and three controversial species of Myodes and Alticola (Rodentia: Arvicolinae)

The genera Myodes (red-backed voles) and Alticola (mountain voles) appear to be sister taxa based on morphological similarities, but molecular analyses fail to resolve them as monophyletic genera owing to the uncertain taxonomic status of Craseomys and Phaulomys. As a result of incomplete sampling of related specimens, ongoing controversies on the taxonomic positions of several generic and specific taxa necessitate further clarifications. Herein, we combined molecular, morphometric, and geometric morphometric approaches to analyze 217 specimens of 10 taxa of Myodes and Alticola systematically. We sequenced three genes (Cytb, COI, GHR) de novo from specimens with fresh tissues, and published sequences for M. shanseius and A. stoliczkanus for the first time. Based on this new molecular dataset, we produced phylogenetic trees using Bayesian inference, maximum likelihood, and maximum parsimony approaches. Our molecular and morphological analyses both identified three primary clades within Myodes and Alticola. The Craseomys-Phaulomys clade consistently separated from Myodes sensu stricto (s. str.) and Alticola s. str.-Platycranius. Our results support the resurrection of the genus Craseomys and the treatment of Phaulomys as its junior synonym. As Craseomys shanseius clustered with C. rufocanus in three gene phylogenies and this assessment was congruent with morphological results, we assigned C. shanseius to a subspecies of C. rufocanus. Specimens from one sampling site in Pulan County of Tibet possess M3 patterns typical of A. stoliczkanus and A. stracheyi, despite clustering together in matrilineal genealogy. Thus, we tentatively assigned A. stracheyi as a junior synonym of A. stoliczkanus. Our analyses confirmed the validity of A. semicanus and unambiguously distinguished it from A. argentatus by the ratio of tail length to head-body length, color of tail and feet, M3 pattern, and distribution.

Phylogeny of Riama (Squamata: Gymnophthalmidae), impact of phenotypic evidence on molecular datasets, and the origin of the Sierra Nevada de Santa Marta endemic fauna

Riama is the most speciose genus of the Neotropical lizard family Gymnophthalmidae. Its more than 30 montane species occur throughout the northern Andes, the Cordillera de la Costa (CC) in Venezuela, and Trinidad. We present the most comprehensive phylogenetic analysis of Riama to date based on a total evidence (TE) approach and direct optimization of molecular and morphological evidence. Analyses use DNA sequences from four loci and 35 phenotypic characters. The dataset consists of 55 ingroup terminals representing 25 of the 30 currently recognized species of Riama plus five undescribed taxa, including an endemic species from the Sierra Nevada de Santa Marta (SNSM) in Colombia, and 66 outgroup terminals of 47 species. Analysis results in a well-supported hypothesis in which Riama is polyphyletic, with its species falling into three clades. The Tepuian Anadia mcdiarmidi nests within one clade of Riama, and the recently resurrected Pantodactylus nests within Cercosaura. Accordingly, we propose a monophyletic taxonomy that reflects historical relationships. Analysis of character evolution indicates that the presence/absence of prefrontals-a cornerstone of the early genus-level taxonomy of cercosaurines-is optimally explained as having been plesiomorphically present in the most recent common ancestor of Cercosaurinae and lost in that of the immediately less inclusive clade. Multiple independent reversals to present and subsequent returns to absent occur within this clade. To evaluate the impact of phenotypic evidence on our results, we compare our TE results with results obtained from analyses using only molecular data. Although phenotypic evidence comprises only 1.2% of the TE matrix, its inclusion alters both the topology and support values of the clades that do not differ. Finally, current phylogenetic evidence reveals a SNSM-CC-Trinidad-tepuis biogeographical link. We hypothesize that an ancient connection facilitated the exchange of species between the SNSM and the CC.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Cenozoic aridization in Central Eurasia shaped diversification of toad-headed agamas (Phrynocephalus; Agamidae, Reptilia)

We hypothesize the phylogenetic relationships of the agamid genus Phrynocephalus to assess how past environmental changes shaped the evolutionary and biogeographic history of these lizards and especially the impact of paleogeography and climatic factors. Phrynocephalus is one of the most diverse and taxonomically confusing lizard genera. As a key element of Palearctic deserts, it serves as a promising model for studies of historical biogeography and formation of arid habitats in Eurasia. We used 51 samples representing 33 of 40 recognized species of Phrynocephalus covering all major areas of the genus. Molecular data included four mtDNA (COI, ND2, ND4, Cytb; 2,703 bp) and four nuDNA protein-coding genes (RAG1, BDNF, AKAP9, NKTR; 4,188 bp). AU-tests were implemented to test for significant differences between mtDNA- and nuDNA-based topologies. A time-calibrated phylogeny was estimated using a Bayesian relaxed molecular clock with nine fossil calibrations. We reconstructed the ancestral area of origin, biogeographic scenarios, body size, and the evolution of habitat preference. Phylogenetic analyses of nuDNA genes recovered a well-resolved and supported topology. Analyses detected significant discordance with the less-supported mtDNA genealogy. The position of Phrynocephalus mystaceus conflicted greatly between the two datasets. MtDNA introgression due to ancient hybridization best explained this result. Monophyletic Phrynocephalus contained three main clades: (I) oviparous species from south-western and Middle Asia; (II) viviparous species of Qinghai–Tibetan Plateau (QTP); and (III) oviparous species of the Caspian Basin, Middle and Central Asia. Phrynocephalus originated in late Oligocene (26.9 Ma) and modern species diversified during the middle Miocene (14.8–13.5 Ma). The reconstruction of ancestral areas indicated that Phrynocephalus originated in Middle East–southern Middle Asia. Body size miniaturization likely occurred early in the history of Phrynocephalus. The common ancestor of Phrynocephalus probably preferred sandy substrates with the inclusion of clay or gravel. The time of Agaminae radiation and origin of Phrynocephalus in the late Oligocene significantly precedes the landbridge between Afro-Arabia and Eurasia in the Early Miocene. Diversification of Phrynocephalus coincides well with the mid-Miocene climatic transition when a rapid cooling of climate drove progressing aridification and the Paratethys salinity crisis. These factors likely triggered the spreading of desert habitats in Central Eurasia, which Phrynocephalus occupied. The origin of the viviparous Tibetan clade has been associated traditionally with uplifting of the QTP; however, further studies are needed to confirm this. Progressing late Miocene aridification, the decrease of the Paratethys Basin, orogenesis, and Plio–Pleistocene climate oscillations likely promoted further diversification within Phrynocephalus. We discuss Phrynocephalus taxonomy in scope of the new analyses.

A new species of Leptolalax (Anura: Megophryidae) from Son Tra Peninsula, central Vietnam

Morphological, acoustic and molecular analyses result in the description of Leptolalax rowleyae sp. nov., a new species of frog in the Megophryidae, belonging to the L. applebyi Rowley Cao species group from central Vietnam. It differs from its congeners by a combination of the following morphological attributes: (1) adult SVL 23.4-25.4 mm in males and 27-27.8 mm in females; (2) presence of distinct dark/brown dorsolateral markings, including black spots on flanks; (3) pinkish milk-white to light brown chest and belly with numerous white speckles; (4) tympanum distinct; (5) absence of webbing or lateral dermal fringes on fingers and toes; (6) pectoral glands comparatively small (3.3-4.7% of SVL); (7) ventrolateral glands indistinct; and (8) iris bicolored with copper tint in upper half fading to golden in lower third of iris. The male advertisement call of the new species consists of 4-6 notes, lacking a distinct introductory note, with an average dominant frequency of 3.2-3.5 kHz. The description of the tadpole constitutes the first description of larval morphology for a member of the L. appleybi species group. Genetically, an uncorrected sequence divergence of 7.4% for 16S rRNA separates the new species from its two closest relatives, L. ardens and L. melicus in the L. applebyi species group.

A new species of genus Fejervarya (Anura: Dicroglossidae) from northern Thailand

We describe a new species of frog in the dicroglossid genus Fejervarya from Ban Monjong, Omkoi District, Chiang Mai Province, northern Thailand. Analysis of DNA sequence data from the mitochondrial gene 16S, advertisement calls, and morphological distinctiveness support recognition of the new species. Matrilineal genealogy suggests that the new population from Chiang Mai is a sister taxon to the South Asian clade that includes F. syhadrensis, F. granosa, and F. pierrei. The new species, Fejervarya chiangmaiensissp. nov., differs morphologically from its congeners by its relatively small body size and proportions and the presence of dorsal warts and dermal ridges. Discovery of this new species indicates that the biodiversity of amphibians in this region remains underestimated.

Potential Pandemic of H7N9 Avian Influenza A Virus in Human

Since 2013, the H7N9 avian influenza A virus (AIV) has caused human infections and to the extent of now surpassing H5N1. This raises an alarm about the potential of H7N9 to become a pandemic problem. Our compilation of the amino acid changes required for AIVs to cross the species-barrier discovers 58 that have very high proportions in both the human- and avian-isolated H7N9 viruses. These changes correspond with sporadic human infections that continue to occur in regions of avian infections. Among the six internal viral genes, amino acid changes do not differ significantly between H9N2 and H7N9, except for V100A in PA, and K526R, D627K, and D701N in PB2. H9N2 AIVs provide internal genes to H7N9. Most of the amino acid changes in H7N9 appear to come directly from H9N2. Seventeen amino acid substitutions appear to have fixed quickly by the 5th wave. Among these, six amino acid sites in HA1 are receptor binding sites, and PB2-A588V was shown to promote the adaptation of AIVs to mammals. The accelerated fixation of mutations may promote the adaptation of H7N9 to human, but need further functional evidence. Although H7N9 AIVs still cannot efficiently transmit between humans, they have the genetic makeup associated with human infections. These viruses must be controlled in poultry to remove the threat of it becoming a human pandemic event.

Capture Hybridization of Long-Range DNA Fragments for High-Throughput Sequencing

Capture hybridization coupled with high-throughput sequencing (HTS) has become one of the most popular approaches to address some scientific problems not only for fundamental evolution but also for ecology and human disease in recent years. However, the technical problem of limited probe capture ability affects its widespread application. Here, we propose to capture hybridize long-range DNA fragments for HTS (termed LR-LCH). We provide a case of three amphibian samples to examine LR-LCH with 2 kb libraries and comparison of standard capture hybridization with 480 bp libraries. Capture sensitivity increased from an average 13.57% of standard capture hybridization to an average 19.80% of LR-LCH; capture efficiency also increased from an average 72.56% of standard capture hybridization to an average 97.71% of LR-LCH. These indicate that longer fragments in the library generally contain both relatively variable regions and relatively conservative regions. The divergent parts of target DNA are enriched along with conservative parts of DNA sequence that effectively captured during hybridization. We present a protocol that allows users to overcome the low capture sensitivity problem for high divergent regions.

Detection of Potential Problematic Cytb Gene Sequences of Fishes in GenBank

Fishes are, by far, the most diverse group of vertebrates. Their classification relies heavily on morphology. In practice, the correct morphological identification of species often depends on personal experience because many species vary in their body shape, color and other external characters. Thus, the identification of a species may be prone to errors. Due to the rapid development of molecular biology, the number of sequences of fishes deposited in GenBank has grown explosively. These published data likely contain errors owing to invalid or incorrectly identified species. The erroneous data can lead to downstream problems. Thus, it is critical that such errors get identified and corrected. A strategy based on DNA barcoding can detect potentially erroneous data, especially when intraspecific K2P variation exceeds interspecific K2P divergence. Analyses of the most used DNA marker for fishes (mitochondrial Cytb) discovers that intraspecific differences of fishes are generally less than 1%, while interspecific differences are generally higher than 10%. Based on this ruler, our analyses identify 1,303 potential problematic Cytb sequences of fishes in GenBank and point to taxonomic problems, errors in identification, genetic introgression and other concerns. Care must be taken to avoid the perpetuation of errors when using these available data.

A new species of the genus Rana from Henan, central China (Anura, Ranidae)

A new species of brown frog Rana luanchuanensis Zhao & Yuan, sp. n. is described from Luanchuan County, western Henan, central China. The mitochondrial genealogy suggests that the new species is the sister taxon to the clade including R. amurensis and R. coreana, and is separated by uncorrected pairwise distances more than 12.5%. Morphologically, this new species differs from its congeners by a suite of characters. Analyses of partial sequences of cytochrome oxidase subunit I (COI) resolve the new species as a single matriline.