The origin of multiple clones in the parthenogenetic lizard species Darevskia rostombekowi

The all-female Caucasian rock lizard Darevskia rostombekowi and other unisexual species of this genus reproduce normally via true parthenogenesis. Typically, diploid parthenogenetic reptiles exhibit some amount of clonal diversity. However, allozyme data from D. rostombekowi have suggested that this species consists of a single clone. Herein, we test this hypothesis by evaluating variation at three variable microsatellite loci for 42 specimens of D. rostombekowi from four populations in Armenia. Analyses based on single nucleotide polymorphisms of each locus reveal five genotypes or presumptive clones in this species. All individuals are heterozygous at the loci. The major clone occurs in 24 individuals and involves three populations. Four rare clones involve one or several individuals from one or two populations. Most variation owes to parent-specific single nucleotide polymorphisms, which occur as heterozygotes. This result fails to reject the hypothesis of a single hybridization founder event that resulted in the initial formation of one major clone. The other clones appear to have originated via post-formation microsatellite mutations of the major clone.

Formal recognition of the species of Oreosaurus (Reptilia, Squamata, Gymnophthalmidae) from the Sierra Nevada de Santa Marta, Colombia

Oreosaurus is one of the two genera extracted from the former Riamasensu lato, which was recently recognized as polyphyletic. Oreosaurus is a small clade (five named and two undescribed species) of montane gymnophthalmid lizards and exhibits an exceptional distributional pattern. Its nominal and undescribed species are discontinuously distributed on the Cordillera de la Costa of Venezuela, the tepuis from the Chimantá massif in Venezuela, the highlands of the island of Trinidad, and the Sierra Nevada de Santa Marta in Colombia (SNSM). Herein, we describe the species of Oreosaurus that is endemic to the SNSM. Historically, this species associates with two names that are currently nomina nuda: Proctoporus serranus and P. specularis. Formal nomenclatural recognition of Oreosaurus serranussp. n. renders specularis a permanently unavailable name for this taxon. Oreosaurus serranussp. n. is the sister of all remaining congeners, and differs primarily from them in having only one pair of genial scales, as well as a unique pattern of scutellation. We provide an identification key to the species of Oreosaurus.

Analysis of the genetic variation in mitochondrial DNA, Y-chromosome sequences, and MC1R sheds light on the ancestry of Nigerian indigenous pigs

Background

The history of pig populations in Africa remains controversial due to insufficient evidence from archaeological and genetic data. Previously, a Western ancestry for West African pigs was reported based on loci that are involved in the determination of coat color. We investigated the genetic diversity of Nigerian indigenous pigs (NIP) by simultaneously analyzing variation in mitochondrial DNA (mtDNA), Y-chromosome sequence and the melanocortin receptor 1 (MC1R) gene.

Results

Median-joining network analysis of mtDNA D-loop sequences from 201 NIP and previously characterized loci clustered NIP with populations from the West (Europe/North Africa) and East/Southeast Asia. Analysis of partial sequences of the Y-chromosome in 57 Nigerian boars clustered NIP into lineage HY1. Finally, analysis of MC1R in 90 NIP resulted in seven haplotypes, among which the European wild boar haplotype was carried by one individual and the European dominant black by most of the other individuals (93%). The five remaining unique haplotypes differed by a single synonymous substitution from European wild type, European dominant black and Asian dominant black haplotypes.

Conclusions

Our results demonstrate a European and East/Southeast Asian ancestry for NIP. Analyses of MC1R provide further evidence. Additional genetic analyses and archaeological studies may provide further insights into the history of African pig breeds. Our findings provide a valuable resource for future studies on whole-genome analyses of African pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0326-1) contains supplementary material, which is available to authorized users.

The Agassiz's desert tortoise genome provides a resource for the conservation of a threatened species

Agassiz's desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.

A new color pattern of the Bungarus candidus complex (Squamata: Elapidae) from Vietnam based on morphological and molecular data

Kraits with black and white bands from Nui Chua National Park, central Vietnam are morphologically similar to the Burmese Krait, Bungarus magnimaculatus, however, analysis of molecular data finds them to be nested within the B. candidus complex.

Taxonomic position of Chinese voles of the tribe Arvicolini and the description of 2 new species from Xizang, China

China has 26 species in the tribe Arvicolini. The taxonomic status of these voles remains controversial despite much effort. Herein, we evaluate the taxonomic position of 22 species plus 2 unidentified taxa using mitochondrial DNA gene sequences (cytb + CO1). We also evaluate 18 species and 2 unidentified taxa using morphological data. Phylogenetic analyses of cytb resolve monophyly for the genera Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys with strong support. Stenocranius clusters with Chionomys but with very weak support. Analyses of concatenated cytb + CO1 resolve the same genera with strong support, but the topology of the tree differs from that of cytb in that Chionomys roots at the base of the tree independent of Stenocranius, which forms the sister-group of Lasiopodomys in a more terminal position. The matrilineal genealogy excludes the type species Arvicola amphibius from the rest of the Arvicolini. This species forms the sister-group of Ondatra with high support. Neodon includes N. irene, N. linzhiensis, N. fuscus, N. leucurus, N. sikimensis, Microtus clarkei, and 2 unidentified specimens. Alexandromys includes the former species Microtus oeconomus, M. kikuchii, M. limnophilus, M. fortis, and M. maximowiczii. Finally, Microtus has the subgenera Blanfordimys, Microtus, Mynomes, Pedomys, Pitymys, and Terricola, which includes the Chinese species M. agrestis, M. arvalis, and Blanfordimys juldaschi. General mixed Yule-coalescent species delimitation modeling demarcates 6 currently recognized species and 2 new species of Neodon. A principal component analysis of the morphological data among 7 matrilines shows that all variables have positive loadings of high magnitude on the 1st component. Canonical discriminant analysis for Neodon (including M. clarkei and 2 unidentified species) correctly classifies 93.0% of specimens. Overall, our analyses support the recognition of Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys as genera. Stenocranius includes Microtus gregalis, and the genealogical position of Stenocranius remains uncertain. The status of Arvicola requires further study. We assign M. clarkei to Neodon and describe 2 new species of Neodon.

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai-Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.

Intraspecific rearrangement of mitochondrial genome suggests the prevalence of the tandem duplication-random loss (TDLR) mechanism in Quasipaa boulengeri

BACKGROUND

Tandem duplication followed by random loss (TDRL) is the most frequently invoked model to explain the diversity of gene rearrangements in metazoan mitogenomes. The initial stages of gene rearrangement are difficult to observe in nature, which limits our understanding of incipient duplication events and the subsequent process of random loss. Intraspecific gene reorganizations may represent intermediate states, and if so they potentially shed light on the evolutionary dynamics of TDRL.

RESULTS

Nucleotide sequences in a hotspot of gene-rearrangement in 28 populations of a single species of frog, Quasipaa boulengeri, provide such predicted intermediate states. Gene order and phylogenetic analyses support a single tandem duplication event and a step-by-step process of random loss. Intraspecific gene rearrangements are not commonly found through comparison of all mitochondrial DNA records of amphibians and squamate reptiles in GenBank.

CONCLUSIONS

The intraspecific variation in Q. boulengeri provides insights into the rate of partial duplications and deletions within a mitogenome, and reveals that fixation and gene-distribution in mitogenomic reorganization is likely non-adaptive.

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis

Amphibian populations are experiencing catastrophic declines driven by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although horizontal gene transfer (HGT) facilitates the evolution and adaptation in many fungi by conferring novel function genes to the recipient fungi, inter-kingdom HGT in Bd remains largely unexplored. In this study, our investigation detects 19 bacterial genes transferred to Bd, including metallo-beta-lactamase and arsenate reductase that play important roles in the resistance to antibiotics and arsenates. Moreover, three probable HGT gene families in Bd are from plants and one gene family coding the ankyrin repeat-containing protein appears to come from oomycetes. The observed multi-copy gene families associated with HGT are probably due to the independent transfer events or gene duplications. Five HGT genes with extracellular locations may relate to infection, and some other genes may participate in a variety of metabolic pathways, and in doing so add important metabolic traits to the recipient. The evolutionary analysis indicates that all the transferred genes evolved under purifying selection, suggesting that their functions in Bd are similar to those of the donors. Collectively, our results indicate that HGT from diverse donors may be an important evolutionary driver of Bd, and improve its adaptations for infecting and colonizing host amphibians.

A new species of kukri snake (Squamata: Colubridae: Oligodon Fitzinger, 1826) from Con Dao Islands, southern Vietnam

We describe a new kukri snake, Oligodon condaoensis sp. nov., from Con Dao Islands, southern Vietnam based on the morphological characters of four specimens. It differs from other congeners by a combination of the following characters: medium size in adults (total length up to 552 mm); 17-17-15 dorsal scale rows; deeply forked hemipenes without spines and papillae, extending to subcaudal 13 or 14; 11-13 maxillary teeth, the posterior three being enlarged; cloacal plate undivided; head scalation complete; nasal divided; presubocular absent; 168-176 ventrals; 33-37 subcaudals; overall dorsal coloration dark gray, faint body stripes present or absent; and ventral coloration cream to dark gray without rectangular blotches. We also provide a list of 43 species of amphibians and terrestrial reptiles recorded from Con Dao Islands.

VennPainter: A Tool for the Comparison and Identification of Candidate Genes Based on Venn Diagrams

VennPainter is a program for depicting unique and shared sets of genes lists and generating Venn diagrams, by using the Qt C++ framework. The software produces Classic Venn, Edwards’ Venn and Nested Venn diagrams and allows for eight sets in a graph mode and 31 sets in data processing mode only. In comparison, previous programs produce Classic Venn and Edwards’ Venn diagrams and allow for a maximum of six sets. The software incorporates user-friendly features and works in Windows, Linux and Mac OS. Its graphical interface does not require a user to have programing skills. Users can modify diagram content for up to eight datasets because of the Scalable Vector Graphics output. VennPainter can provide output results in vertical, horizontal and matrix formats, which facilitates sharing datasets as required for further identification of candidate genes. Users can obtain gene lists from shared sets by clicking the numbers on the diagram. Thus, VennPainter is an easy-to-use, highly efficient, cross-platform and powerful program that provides a more comprehensive tool for identifying candidate genes and visualizing the relationships among genes or gene families in comparative analysis.

Red River barrier and Pleistocene climatic fluctuations shaped the genetic structure of Microhyla fissipes complex (Anura: Microhylidae) in southern China and Indochina

South China and Indochina host striking species diversity and endemism. Complex tectonic and climatic evolutions appear to be the main drivers of the biogeographic patterns. In this study, based on the geologic history of this region, we test 2 hypotheses using the evolutionary history of Microhyla fissipes species complex. Using DNA sequence data from both mitochondrial and nuclear genes, we first test the hypothesis that the Red River is a barrier to gene flow and dispersal. Second, we test the hypothesis that Pleistocene climatic cycling affected the genetic structure and population history of these frogs. We detect 2 major genetic splits that associate with the Red River. Time estimation suggests that late Miocene tectonic movement associated with the Red River drove their diversification. Species distribution modeling (SDM) resolves significant ecological differences between sides of the Red River. Thus, ecological divergence also probably promoted and maintained the diversification. Genogeography, historical demography, and SDM associate patterns in southern China with climate changes of the last glacial maximum (LGM), but not Indochina. Differences in geography and climate between the 2 areas best explain the discovery. Responses to the Pleistocene glacial–interglacial cycling vary among species and regions.

Taxonomic revision of the Chinese Limnonectes (Anura, Dicroglossidae) with the description of a new species from China and Myanmar

Phylogenetic reconstructions derived from DNA sequence data play a central role in documenting the number of species in a complex. Such analyses are pointing to the existence of many cryptic species, especially in poorly understood groups such as the genus Limnonectes, and the L. kuhlii species complex in particular. To understand the Limnonectes frogs of China, we reconstruct the major matrilineal genealogy of Limnonectes from China and Southeast Asia based on 12S rRNA, tRNA^Val and 16S rRNA gene sequences. Based on new data we recognize five species of Limnonectes in China including L. bannaensis, L. fujianensis, L. fragilis, L. taylori (new record), and a new species from southern China and Myanmar. Phylogenetically, the new species is more closely related to the clade comprising L. taylori, L. megastomias, L. isanensis, L. nguyenorum, and L. jarujini from Thailand than to other Chinese species. This study supports previous findings of sympatric members of a species complex that are not each other's closest relatives.

The desert tortoise trichotomy: Mexico hosts a third, new sister-species of tortoise in the Gopherus morafkai-G. agassizii group

Desert tortoises (Testudines; Testudinidae; Gopherusagassizii group) have an extensive distribution throughout the Mojave, Colorado, and Sonoran desert regions. Not surprisingly, they exhibit a tremendous amount of ecological, behavioral, morphological and genetic variation. Gopherusagassizii was considered a single species for almost 150 years but recently the species was split into the nominate form and Morafka’s desert tortoise, Gopherusmorafkai, the latter occurring south and east of the Colorado River. Whereas a large body of literature focuses on tortoises in the United States, a dearth of investigations exists for Mexican animals. Notwithstanding, Mexican populations of desert tortoises in the southern part of the range of Gopherusmorafkai are distinct, particularly where the tortoises occur in tropical thornscrub and tropical deciduous forest. Recent studies have shed light on the ecology, morphology and genetics of these southern ‘desert’ tortoises. All evidence warrants recognition of this clade as a distinctive taxon and herein we describe it as Gopherusevgoodeisp. n. The description of the new species significantly reduces and limits the distribution of Gopherusmorafkai to desertscrub habitat only. By contrast, Gopherusevgoodeisp. n. occurs in thornscrub and tropical deciduous forests only and this leaves it with the smallest range of the three sister species. We present conservation implications for the newly described Gopherusevgoodei, which already faces impending threats.

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

Genomic incompatibilities in the diploid and tetraploid offspring of the goldfish × common carp cross

Polyploidy is much rarer in animals than in plants but it is not known why. The outcome of combining two genomes in vertebrates remains unpredictable, especially because polyploidization seldom shows positive effects and more often results in lethal consequences because viable gametes fail to form during meiosis. Fortunately, the goldfish (maternal) × common carp (paternal) hybrids have reproduced successfully up to generation 22, and this hybrid lineage permits an investigation into the genomics of hybridization and tetraploidization. The first two generations of these hybrids are diploids, and subsequent generations are tetraploids. Liver transcriptomes from four generations and their progenitors reveal chimeric genes (>9%) and mutations of orthologous genes. Characterizations of 18 randomly chosen genes from genomic DNA and cDNA confirm the chimera. Some of the chimeric and differentially expressed genes relate to mutagenesis, repair, and cancer-related pathways in 2nF1. Erroneous DNA excision between homologous parental genes may drive the high percentage of chimeric genes, or even more potential mechanisms may result in this phenomenon. Meanwhile, diploid offspring show paternal-biased expression, yet tetraploids show maternal-biased expression. These discoveries reveal that fast and unstable changes are mainly deleterious at the level of transcriptomes although some offspring still survive their genomic abnormalities. In addition, the synthetic effect of genome shock might have resulted in greatly reduced viability of 2nF2 hybrid offspring. The goldfish × common carp hybrids constitute an ideal system for unveiling the consequences of intergenomic interactions in hybrid vertebrate genomes and their fertility.

Phylogenetic Relationships of Five Asian Schilbid Genera Including Clupisoma (Siluriformes: Schilbeidae)

The phylogenetic relationships of Asian schilbid catfishes of the genera Clupisoma, Ailia, Horabagrus, Laides and Pseudeutropius are poorly understood, especially those of Clupisoma. Herein, we reconstruct the phylogeny of 38 species of catfishes belonging to 28 genera and 14 families using the concatenated mitochondrial genes COI, cytb, and 16S rRNA, as well as the nuclear genes RAG1 and RAG2. The resulting phylogenetic trees consistently place Clupisoma as the sister taxon of Laides, and the five representative Asian schilbid genera form two monophyletic groups with the relationships (Ailia (Laides, Clupisoma)) and (Horabagrus, Pseudeutropius). The so-called "Big Asia" lineage relates distantly to African schilbids. Independent analyses of the mitochondrial and nuclear DNA data yield differing trees for the two Asian schilbid groups. Analyses of the mitochondrial gene data support a sister-group relationship for (Ailia (Laides, Clupisoma)) and the Sisoroidea and a sister-taxon association of (Horabagrus, Pseudeutropius) and the Bagridae. In contrast, analyses of the combined nuclear data indicate (Ailia (Laides, Clupisoma)) to be the sister group to (Horabagrus, Pseudeutropius). Our results indicate that the Horabagridae, recognized by some authors as consisting of Horabagrus, Pseudeutropius and Clupisoma does not include the latter genus. We formally erect a new family, Ailiidae fam. nov. for a monophyletic Asian group comprised of the genera Ailia, Laides and Clupisoma.

Horizontal functional gene transfer from bacteria to fishes

Invertebrates can acquire functional genes via horizontal gene transfer (HGT) from bacteria but fishes are not known to do so. We provide the first reliable evidence of one HGT event from marine bacteria to fishes. The HGT appears to have occurred after emergence of the teleosts. The transferred gene is expressed and regulated developmentally. Its successful integration and expression may change the genetic and metabolic repertoire of fishes. In addition, this gene contains conserved domains and similar tertiary structures in fishes and their putative donor bacteria. Thus, it may function similarly in both groups. Evolutionary analyses indicate that it evolved under purifying selection, further indicating its conserved function. We document the first likely case of HGT of functional gene from prokaryote to fishes. This discovery certifies that HGT can influence vertebrate evolution.

Genetic adaptations of the plateau zokor in high-elevation burrows

The plateau zokor (Myospalax baileyi) spends its entire life underground in sealed burrows. Confronting limited oxygen and high carbon dioxide concentrations, and complete darkness, they epitomize a successful physiological adaptation. Here, we employ transcriptome sequencing to explore the genetic underpinnings of their adaptations to this unique habitat. Compared to Rattus norvegicus, genes belonging to GO categories related to energy metabolism (e.g. mitochondrion and fatty acid beta-oxidation) underwent accelerated evolution in the plateau zokor. Furthermore, the numbers of positively selected genes were significantly enriched in the gene categories involved in ATPase activity, blood vessel development and respiratory gaseous exchange, functional categories that are relevant to adaptation to high altitudes. Among the 787 genes with evidence of parallel evolution, and thus identified as candidate genes, several GO categories (e.g. response to hypoxia, oxygen homeostasis and erythrocyte homeostasis) are significantly enriched, are two genes, EPAS1 and AJUBA, involved in the response to hypoxia, where the parallel evolved sites are at positions that are highly conserved in sequence alignments from multiple species. Thus, accelerated evolution of GO categories, positive selection and parallel evolution at the molecular level provide evidences to parse the genetic adaptations of the plateau zokor for living in high-elevation burrows.

Two New Species of Liuixalus (Rhacophoridae, Anura): Evidence from Morphological and Molecular Analyses

Due to small body sizes, superficial similarities in morphologies, and obscure activity behaviors, the phylogeny and taxonomy of species in the genus Liuixalus were very troublesome. Some species might comprise a complex of cryptic species. To investigate the species of group, we constructed the matrilineal genealogy of the genus using 16s rRNA mitochondrial DNA sequences. Analyses recovered six well supported matrilines that involved L. romeri, L. ocellatus, L. hainanus, L. calcarius, Liuixalus shiwandashansp. nov. and Liuixalus jinxiuensissp. nov., though the historical relationships among them remained unresolved. Currently, Liuixalus included 4 species, distributed eastwards from northern Vietnam to Hong Kong, China. Based on genealogical and morphological distinctiveness, we described Liuixalus jinxiuensissp. nov. from the type locality Mt. Dayao, Jinxiu, China and Liuixalus shiwandashansp. nov. from the type locality Mt. Shiwanda, China. A combination of morphological measurements, genetic, bioacoustic and osteological analysis was served to diagnose the new taxa.