A versatile and highly efficient toolkit including 102 nuclear markers for vertebrate phylogenomics, tested by resolving the higher level relationships of the caudata

A time-calibrated salamander phylogeny including 765 species and 503 genes

Recent time-calibrated amphibian phylogenies agree on the family-level relationships among extant salamanders but had disparate sampling regimes and inferred very different divergence times. For example, a recent phylogenomic study based on 220 nuclear loci had limited taxon sampling (41 species) and estimated relatively young divergence dates, whereas a more extensive supermatrix study based on 15 genes and 481 species estimated dates that were 22-45 million years older for major clades. Here, we combined phylogenomic and supermatrix approaches to estimate the largest salamander phylogeny to date based on molecular markers. Our matrix contained 765 salamander species and 503 genes (with 92.3% missing data overall). We included 284 more species than the previous largest salamander phylogeny (59% increase) and sampled approximately 93% of all currently described salamander species. Our dating analyses incorporated more than twice as many fossil calibration points within salamanders as previous studies. Maximum-likelihood estimates of tree topology yielded family-level relationships that were consistent with earlier studies. Nearly all species were placed in the expected genera, despite extensive missing data in many species. Bootstrap support was generally high across the tree but was poor in some clades where sampling of genes was limited (e.g., among some bolitoglossine salamanders). The dating analyses yielded age estimates for major clades that were generally intermediate between those from the previous phylogenomic and supermatrix analyses. We also provide a set of 200 time-calibrated trees for use in comparative analyses.

A comprehensive molecular phylogeny of the genus Sylvirana (Anura: Ranidae) highlights unrecognized diversity, revised classification and historical biogeography

The genus Sylvirana includes 12 species widely distributed in South China and Southeast Asia. The phylogenetic relationships and species diversity for Sylvirana and allied genera remain unresolved and controversial due to insufficient data and incomplete taxon sampling. Using a combined dataset of mitochondrial genes (16S and COI) and 101 nuclear genes obtained through the amplicon sequence capture approach, we generated the most comprehensive phylogenetic analysis for the genus Sylvirana to date, inferring diversity, phylogenetic relationships, and historical biogeography with unprecedented levels of taxon and geographic sampling. Our results conservatively reveal six undescribed species, mostly distributed in peninsular Indochina. Phylogenetic analyses strongly support the non-monophyly of Sylvirana with respect to Pterorana. Additionally, phylogenetic results place Sylvirana guentheri and Pelophylax lateralis into genus Humerana, supporting the inclusion of Hylarana latouchii, Papurana milleti, and Hylarana attigua within Pterorana + Sylvirana. The long-disputed species of Hylarana bannanica (previously Sylvirana) cluster with genus Papurana. Because the results of multiple non-monophyletic genera create taxonomic confusion, we suggest relegating all genera to subgenus rank of Hylarana. Sylvirana is a junior synonym of the Pterorana. Biogeographically, we trace the origin of Pterorana to Southeast Asia during the early Miocene, with subsequent dispersal thereafter. Our study shows that climatic changes may have profoundly influenced the diversification of Pterorana during the Miocene.

Incorporating New Datatypes to Enhance Species Delimitation: A Case Study in Rice Paddy Snakes (Homalopsidae: Hypsiscopus)

Homalopsids (Old World Mud Snakes) include 59 semiaquatic species in Asia and Australasia that display an array of morphological adaptations, behaviors, and microhabitat preferences. These attributes make homalopsids an ideal model system for broader questions in evolutionary biology, but the diversity of this understudied group of snakes is still being described. Recognized species diversity in rice paddy snakes (Hypsiscopus) has recently doubled after nearly 200 years of taxonomic stability. However, the evolutionary distinctiveness of some populations remains in question. In this study, we compare mainland Southeast Asian populations of Hypsiscopus east and west of the Red River Basin in Vietnam, a known biogeographic barrier in Asia, using an iterative approach with molecular phylogenetic reconstruction, machine-learning morphological quantitative statistics, and ecological niche modeling. Our analyses show that populations west of the Red River Basin represent an independent evolutionary lineage that is distinct in genetics, morphospace, and habitat suitability, and so warrants species recognition. The holotype of H. wettsteini, a species originally described in error from Costa Rica, grouped morphometrically with the population at the Red River Basin and eastward, and those west of the Red River Basin are referred to the recently described H. murphyi. The two species may have diversified due to a variety of geological and environmental factors, and their recognition exemplifies the importance of multifaceted approaches in taxonomy for downstream biogeographic studies on speciation scenarios.

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

UPrimer: A Clade-Specific Primer Design Program Based on Nested-PCR Strategy and Its Applications in Amplicon Capture Phylogenomics

Amplicon capture is a promising target sequence capture approach for phylogenomic analyses, and the design of clade-specific nuclear protein-coding locus (NPCL) amplification primers is crucial for its successful application. In this study, we developed a primer design program called UPrimer that can quickly design clade-specific NPCL amplification primers based on genome data, without requiring manual intervention. Unlike other available primer design programs, UPrimer uses a nested-PCR strategy that greatly improves the amplification success rate of the designed primers. We examined all available metazoan genome data deposited in NCBI and developed NPCL primer sets for 21 metazoan groups with UPrimer, covering a wide range of taxa, including arthropods, mollusks, cnidarians, echinoderms, and vertebrates. On average, each clade-specific NPCL primer set comprises ∼1,000 NPCLs. PCR amplification tests were performed in 6 metazoan groups, and the developed primers showed a PCR success rate exceeding 95%. Furthermore, we demonstrated a phylogenetic case study in Lepidoptera, showing how NPCL primers can be used for phylogenomic analyses with amplicon capture. Our results indicated that using 100 NPCL probes recovered robust high-level phylogenetic relationships among butterflies, highlighting the utility of the newly designed NPCL primer sets for phylogenetic studies. We anticipate that the automated tool UPrimer and the developed NPCL primer sets for 21 metazoan groups will enable researchers to obtain phylogenomic data more efficiently and cost-effectively and accelerate the resolution of various parts of the Tree of Life.

Frog phylogeny: A time-calibrated, species-level tree based on hundreds of loci and 5,242 species

Large-scale, time-calibrated phylogenies from supermatrix studies have become crucial for evolutionary and ecological studies in many groups of organisms. However, in frogs (anuran amphibians), there is a serious problem with existing supermatrix estimates. Specifically, these trees are based on a limited number of loci (15 or fewer), and the higher-level relationships estimated are discordant with recent phylogenomic estimates based on much larger numbers of loci. Here, we attempted to rectify this problem by generating an expanded supermatrix and combining this with data from phylogenomic studies. To assist in aligning ribosomal sequences for this supermatrix, we developed a new program (TaxonomyAlign) to help perform taxonomy-guided alignments. The new combined matrix contained 5,242 anuran species with data from 307 markers, but with 95% missing data overall. This dataset represented a 71% increase in species sampled relative to the previous largest supermatrix analysis of anurans (adding 2,175 species). Maximum-likelihood analyses generated a tree in which higher-level relationships (and estimated clade ages) were generally concordant with those from phylogenomic analyses but were more discordant with the previous largest supermatrix analysis. We found few obvious problems arising from the extensive missing data in most species. We also generated a set of 100 time-calibrated trees for use in comparative analyses. Overall, we provide an improved estimate of anuran phylogeny based on the largest number of combined taxa and markers to date. More broadly, we demonstrate the potential to combine phylogenomic and supermatrix analyses in other groups of organisms.

The salamander blastema within the broader context of metazoan regeneration

Throughout the animal kingdom regenerative ability varies greatly from species to species, and even tissue to tissue within the same organism. The sheer diversity of structures and mechanisms renders a thorough comparison of molecular processes truly daunting. Are "blastemas" found in organisms as distantly related as planarians and axolotls derived from the same ancestral process, or did they arise convergently and independently? Is a mouse digit tip blastema orthologous to a salamander limb blastema? In other fields, the thorough characterization of a reference model has greatly facilitated these comparisons. For example, the amphibian Spemann-Mangold organizer has served as an amazingly useful comparative template within the field of developmental biology, allowing researchers to draw analogies between distantly related species, and developmental processes which are superficially quite different. The salamander limb blastema may serve as the best starting point for a comparative analysis of regeneration, as it has been characterized by over 200 years of research and is supported by a growing arsenal of molecular tools. The anatomical and evolutionary closeness of the salamander and human limb also add value from a translational and therapeutic standpoint. Tracing the evolutionary origins of the salamander blastema, and its relatedness to other regenerative processes throughout the animal kingdom, will both enhance our basic biological understanding of regeneration and inform our selection of regenerative model systems.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

FishPIE: a universal phylogenetically informative exon markers set for ray-finned fishes

Understanding the evolutionary history of the highly diverse ray-finned fishes has been challenging, and the development of more universal primers for phylogenetic analyses may help overcoming these challenges. We developed FishPIE, a nested PCR primer set of 82 phylogenetically informative exon markers, and tested it on 203 species from 31 orders of Actinopterygii. We combined orthologous sequences of the FishPIE markers obtained from published genomes and transcriptomes and constructed the phylogeny of 710 species belonging to 190 families and 60 orders. The resulting phylogenies had topologies comparable to previous phylogenomic studies. We demonstrated that the FishPIE markers could address phylogenetic questions across broad taxonomic levels. By incorporating the newly sequenced taxa, we were able to shed new light on the phylogeny of the highly diverse Cypriniformes. Thus, FishPIE holds great promise for generating genetic data for broad taxonomic groups and accelerating our understanding of the fish tree of life.

•

FishPIE is a nested PCR primer set of 82 markers for fish phylogenetic analysis

•

The markers can be broadly applied to all orders of ray-finned fishes

•

Their phylogenetic performance is comparable to that of genomic analyses

Systematic position of the Clicking Frog (Kassinula Laurent, 1940), the problem of chimeric sequences and the revised classification of the family Hyperoliidae

The systematics of the African frog family Hyperoliidae has undergone turbulent changes in last decades. Representatives of several genera have not been genetically investigated or with only limited data, and their phylogenetic positions are thus still not reliably known. This is the case of the De Witte's Clicking Frog (Kassinula wittei) which belongs to a monotypic genus. This miniature frog occurs in a poorly studied region, southeastern Democratic Republic of the Congo, northern Zambia, Angola. So far it is not settled whether this genus belongs to the subfamily Kassininae as a relative of the genus Kassina, or to the subfamily Hyperoliinae as a relative of the genus Afrixalus. Here we present for the first time a multilocus phylogenetic reconstruction (using five nuclear and one mitochondrial marker) of the family Hyperoliidae, including Kassinula. We demonstrate with high confidence that Kassinula is a member of Hyperoliinae belonging to a clade also containing Afrixalus (sub-Saharan Africa), Heterixalus (Madagascar) and Tachycnemis (Seychelles). We find that Kassinula represents a divergent lineage (17-25 Mya), which supports its separate genus-level status, but its exact systematic position remains uncertain. We propose to name the clade to which the above four genera belong as the tribe Tachycnemini Channing, 1989. A new taxonomy of the family Hyperoliidae was recently proposed by Dubois et al. (2021: Megataxa 5, 1-738). We demonstrate here that the new taxonomy was based on a partially erroneous phylogenetic reconstruction resulting from a supermatrix analysis of chimeric DNA sequences combining data from two families, Hyperoliidae and Arthroleptidae (the case of Cryptothylax). We therefore correct the erroneous part and propose a new, revised suprageneric taxonomy of the family Hyperoliidae. We also emphasize the importance of inspecting individual genetic markers before their concatenation or coalescent-based tree reconstructions to avoid analyses of chimeric DNA sequences producing incorrect phylogenetic reconstructions. Especially when phylogenetic reconstructions are used to propose taxonomies and systematic classifications.

Endemism, invasion, and overseas dispersal: the phylogeographic history of the Lesser Antillean frog, Eleutherodactylus johnstonei

Cryptogenic species are those whose native and introduced ranges are unknown. The extent and long history of human migration rendered numerous species cryptogenic. Incomplete knowledge regarding the origin and native habitat of a species poses problems for conservation management and may confound ecological and evolutionary studies. The Lesser Antilles pose a particular challenge with regard to cryptogenic species because these islands have been anthropogenically connected since before recorded history. Here, we use population genetic and phylogeographic tools in an attempt to determine the origin of Eleutherodactylus johnstonei, a frog species with a potentially widespread introduced range and whose native range within the Lesser Antilles is unknown. Based on elevated estimates of genetic diversity and within-island geographic structure not present elsewhere in the range, we identify Montserrat as the native island of E. johnstonei. We also document two major clades within E. johnstonei, only one of which is the primary source of introduced populations throughout the Americas. Our results demonstrate the utility of genetic tools for resolving cryptogenic species problems and highlight E. johnstonei as a potential system for understanding differences in invasive potential among sister lineages.

Congolius, a new genus of African reed frog endemic to the central Congo: A potential case of convergent evolution

The reed frog genus Hyperolius (Afrobatrachia, Hyperoliidae) is a speciose genus containing over 140 species of mostly small to medium-sized frogs distributed in sub-Saharan Africa. Its high level of colour polymorphism, together with in anurans relatively rare sexual dichromatism, make systematic studies more difficult. As a result, the knowledge of the diversity and taxonomy of this genus is still limited. Hyperolius robustus known only from a handful of localities in rain forests of the central Congo Basin is one of the least known species. Here, we have used molecular methods for the first time to study the phylogenetic position of this taxon, accompanied by an analysis of phenotype based on external (morphometric) and internal (osteological) morphological characters. Our phylogenetic results undoubtedly placed H. robustus out of Hyperolius into a common clade with sympatric Cryptothylax and West African Morerella. To prevent the uncovered paraphyly, we place H. robustus into a new genus, Congolius. The review of all available data suggests that the new genus is endemic to the central Congolian lowland rain forests. The analysis of phenotype underlined morphological similarity of the new genus to some Hyperolius species. This uniformity of body shape (including cranial shape) indicates that the two genera have either retained ancestral morphology or evolved through convergent evolution under similar ecological pressures in the African rain forests.

Optimizing the widely used nuclear protein-coding gene primers in beetle phylogenies and their application in the genus Sasajiscymnus Vandenberg (Coleoptera: Coccinellidae)

Advances in genomic biology and the increasing availability of genomic resources allow developing hundreds of nuclear protein-coding (NPC) markers, which can be used in phylogenetic research. However, for low taxonomic levels, it may be more practical to select a handful of suitable molecular loci for phylogenetic inference. Unfortunately, the presence of degenerate primers of NPC markers can be a major impediment, as the amplification success rate is low and they tend to amplify nontargeted regions. In this study, we optimized five NPC fragments widely used in beetle phylogenetics (i.e., two parts of carbamoyl-phosphate synthetase: CADXM and CADMC, Topoisomerase, Wingless and Pepck) by reducing the degenerate site of primers and the length of target genes slightly. These five NPC fragments and 6 other molecular loci were amplified to test the monophyly of the coccinellid genus Sasajiscymnus Vandenberg. The analysis of our molecular data set clearly supported the genus Sasajiscymnus may be monophyletic but confirmation with an extended sampling is required. A fossil-calibrated chronogram was generated by BEAST, indicating an origin of the genus at the end of the Cretaceous (77.87 Myr). Furthermore, a phylogenetic informativeness profile was generated to compare the phylogenetic properties of each gene more explicitly. The results showed that COI provides the strongest phylogenetic signal among all the genes, but Pepck, Topoisomerase, CADXM and CADMC are also relatively informative. Our results provide insight into the evolution of the genus Sasajiscymnus, and also enrich the molecular data resources for further study.

A new species of the horned toad Megophrys Kuhl & Van Hasselt, 1822 (Anura, Megophryidae) from southwest China

A new species of the genus Megophrys is described from Guizhou Province, China. Molecular phylogenetic analyses based on mitochondrial DNA and nuclear DNA sequences all strongly supported the new species as an independent clade sister to M.minor and M.jiangi. The new species could be distinguished from its congeners by a combination of the following characters: body size moderate (SVL 43.4–44.1 mm in males, and 44.8–49.8 mm in females; vomerine teeth absent; tongue not notched behind; a small horn-like tubercle at the edge of each upper eyelid; tympanum distinctly visible, rounded; two metacarpal tubercles on palm; relative finger lengths II < I < V < III; toes without webbing; heels overlapping when thighs are positioned at right angles to the body; tibiotarsal articulation reaching the level between tympanum and eye when leg stretched forward; in breeding males, an internal single subgular vocal sac in male, and the nuptial pads with black spines on dorsal surface of bases of the first two fingers.

A large-scale systematic framework of Chinese snakes based on a unified multilocus marker system

Snakes are one of the most diverse groups of terrestrial vertebrates, with approximately 3500 extant species. A robust phylogeny and taxonomy of snakes is crucial for us to know, study and protect them. For a large group such as snakes, broad-scale phylogenetic reconstructions largely rely on data integration. Increasing the compatibility of the data from different researches is thus important, which can be facilitated by standardization of the loci used in systematic analyses. In this study, we proposed a unified multilocus marker system for snake systematics by conflating 5 mitochondrial markers, 19 vertebrate-universal nuclear protein coding (NPC) markers and 72 snake-specific noncoding intron markers. This marker system is an addition to the large squamate conserved locus set (SqCL) for studies preferring a medium-scale data set. We applied this marker system to over 440 snake samples and constructed the currently most comprehensive systematic framework of the snakes in China. Robust snake phylogenetic relationships were recovered at both deep and shallow evolutionary depths, demonstrating the usefulness of this multilocus marker system. Discordance was revealed by a parallel comparison between the snake tree based on the multilocus marker system and that based on only the mitochondrial loci, highlighting the necessity of using multiple types of markers to better understand the snake evolutionary histories. The divergence times of different snake groups were estimated with the nuclear data set. Our comprehensive snake tree not only confirms many important nodes inferred in previous studies but also contributes new insights into many snake phylogenetic relationships. Suggestions are made for the current Chinese snake taxonomy.

Two novel cricetine mitogenomes: Insight into the mitogenomic characteristics and phylogeny in Cricetinae (Rodentia: Cricetidae)

Both Cricetus cricetus and Phodopus sungorus mitochondrial genomes (mitogenomes) were sequenced and elaborated for the first time in the present study. Their mitogenomes contained 37 genes and showed typical characteristics of the vertebrate mitogenome. Comparative analysis of 10 cricetine mitogenomes indicated that they shared similar characteristics with those of other cricetines in terms of genes arrangement, nucleotide composition, codon usage, tRNA structure, nucleotide skew and the origin of replication of light strand. Phylogenetic relationship of the subfamily Cricetinae was reconstructed using mitogenomes data with the methods of Bayesian Inference and Maximum Likelihood. Phylogenetic analysis indicated that Cricetulus kamensis was at basal position and phylogenetically distant from all other Cricetulus species but had a close relationship with the group of Phodopus, and supported that the genus Urocricetus deserved as a separate genus rank. The phylogenetic status of Tscherskia triton represented a separate clade corresponding to a diversified cricetine lineage (Cricetulus, Allocricetulus, and Cricetus).

Sequence capture across large phylogenetic scales by using pooled PCR-generated baits: A case study of Lepidoptera

Sequence capture across large phylogenetic scales is not easy because hybridization capture is only effective when the genetic distance between the bait and target is small. Here, we propose a simple but effective strategy to tackle this issue: pooling DNA from a number of selected representative species of different clades to prepare PCR-generated baits to minimize the genetic distance between the bait and target. To demonstrate the utility of this strategy, we newly developed a set of universal nuclear markers (including 94 nuclear protein-coding genes) for Lepidoptera, a superdiverse insect group. We used a DNA pool from six lepidopteran species (representing six superfamilies) to prepare PCR baits for the 94 markers. These homemade PCR baits were used to capture sequence data from 43 species of 17 lepidopteran families, and 94% of the target loci were recovered. We constructed two data sets from the obtained data (one containing ~90 kb target coding sequences and the other containing ~120 kb target + flanking coding sequences). Both data sets yielded highly similar and well-resolved trees with 90% of nodes having >95% bootstrap support. Our capture experiment indicated that using DNA mixtures pooled from different clade-representative species of Lepidoptera to prepare PCR baits can reliably capture a large number of targeted nuclear markers across different Lepidoptera lineages. We hope that this newly developed nuclear marker set will serve as a new phylogenetic tool for Lepidoptera phylogenetics, and the PCR bait preparation strategy can facilitate the application of sequence capture techniques by researchers to accelerate data collection.

Phylogeography of the Tibetan hamster Cricetulus kamensis in response to uplift and environmental change in the Qinghai-Tibet Plateau

AIM

The evolutionary process of an organism provides valuable data toward an understanding of the Earth evolution history. To investigate the relationship between the uplift of the Qinghai-Tibet Plateau (QTP) and mammalian evolution since the late Cenozoic, the geographic distribution of genetic variations in the Tibetan hamster Cricetulus kamensis was investigated using phylogeographical methods. In particular, population divergence, demographic history, genetic variation, and the prediction of species distribution area were investigated.

LOCATION

The Qinghai-Tibet Plateau.

METHODS

A total of 53 specimens, representing 13 geographic populations, were collected from the QTP. The phylogeographical pattern and demographic history of C. kamensis were analyzed, and the probable factors in the QTP uplift and the Quaternary glacial periods were inferred from one nuclear and four mitochondrial genes. Furthermore, the species distribution model (SDM) was used to predict changes in potentially suitable habitats since the last Interglacial.

RESULTS

Phylogenetic analysis demonstrated that two major genetic differentiations of the C. kamensis population occurred during the Early Pleistocene that were influenced by the Qing-Zang tectonic movement from the Middle Pliocene to the Early Pleistocene. Genetic distance between two major clades indicated low genetic divergence. Demographic history analysis showed that the C. kamensis population was affected by the Quaternary glacial period. SDM analysis indicated that C. kamensis was endemic to the QTP and the suitable habitat was affected by climate change, especially during the Last Glacial Maximum (LGM).

MAIN CONCLUSION

Our results indicated that the QTP uplift led to the population divergence of C. kamensis, and vicariance well accounted for the geographic distribution of genetic variation in C. kamensis as a result of genetic divergence and lack of gene flow. The genetic distance shows that C. alticola may be a subspecies of C. kamensis. Demographic history analysis suggests that the QTP was affected by the last glacial period. SDM analysis supports that almost the entire QTP is covered by a huge ice sheet during the LGM.

The first set of universal nuclear protein-coding loci markers for avian phylogenetic and population genetic studies

Multiple nuclear markers provide genetic polymorphism data for molecular systematics and population genetic studies. They are especially required for the coalescent-based analyses that can be used to accurately estimate species trees and infer population demographic histories. However, in avian evolutionary studies, these powerful coalescent-based methods are hindered by the lack of a sufficient number of markers. In this study, we designed PCR primers to amplify 136 nuclear protein-coding loci (NPCLs) by scanning the published Red Junglefowl (Gallus gallus) and Zebra Finch (Taeniopygia guttata) genomes. To test their utility, we amplified these loci in 41 bird species representing 23 Aves orders. The sixty-three best-performing NPCLs, based on high PCR success rates, were selected which had various mutation rates and were evenly distributed across 17 avian autosomal chromosomes and the Z chromosome. To test phylogenetic resolving power of these markers, we conducted a Neoavian phylogenies analysis using 63 concatenated NPCL markers derived from 48 whole genomes of birds. The resulting phylogenetic topology, to a large extent, is congruence with results resolved by previous whole genome data. To test the level of intraspecific polymorphism in these makers, we examined the genetic diversity in four populations of the Kentish Plover (Charadrius alexandrinus) at 17 of NPCL markers chosen at random. Our results showed that these NPCL markers exhibited a level of polymorphism comparable with mitochondrial loci. Therefore, this set of pan-avian nuclear protein-coding loci has great potential to facilitate studies in avian phylogenetics and population genetics.

A new species of the odorous frog genus Odorrana (Amphibia, Anura, Ranidae) from southwestern China

The genus Odorrana is widely distributed in the mountains of East and Southeastern Asia. An increasing number of new species in the genus have been recognized especially in the last decade. Phylogenetic studies of the O. schmackeri species complex with wide distributional range also revealed several cryptic species. Here, we describe a new species in the species complex from Guizhou Province of China. Phylogenetic analyses based on mitochondrial DNA indicated the new species as a monophyly clustered into the Odorrana clade and sister to O. schmackeri, and nuclear DNA also indicated it as an independent lineage separated from its related species. Morphologically, the new species can be distinguished from its congeners based on a combination of the following characters: (1) having smaller body size in males (snout-vent length (SVL) <43.3 mm); (2) head longer than wide; (3) dorsolateral folds absent; (4) tympanum of males large and distinct, tympanum diameter twice as long as width of distal phalanx of finger III; (5) two metacarpal tubercles; (6) relative finger lengths: II < I < IV < III; (7) tibiotarsal articulation reaching to the level between eye to nostril when leg stretched forward; (8) disks on digits with circum-marginal grooves; (9) toes fully webbed to disks; (10) the first subarticular tubercle on fingers weak; (11) having white pectoral spinules, paired subgular vocal sacs located at corners of throat, light yellow nuptial pad on the first finger in males.

Prevalence of cryptic species in morphologically uniform taxa - Fast speciation and evolutionary radiation in Asian frogs

Diversity and distributions of cryptic species have long been a vexing issue. Identification of species boundaries is made difficult by the lack of obvious morphological differences. Here, we investigate the cryptic diversity and evolutionary history of an underappreciated group of Asian frog species (Megophrys) to explore the pattern and dynamic of amphibian cryptic species. We sequenced four mitochondrial genes and five nuclear genes and delineated species using multiple approaches, combining DNA and mating-call data. A Bayesian species tree was generated to estimate divergence times and to reconstruct ancestral ranges. Macroevolutionary analyses and hybridization tests were conducted to explore the evolutionary dynamics of this cryptic group. Our phylogenies support the current subgenera. We revealed 43 cryptic species, 158% higher than previously thought. The species-delimitation results were further confirmed by mating-call data and morphological divergence. We found that these Asian frogs entered China from the Sunda Shelf 48 Mya, followed by an ancient radiation event during middle Miocene. We confirmed the efficiency of the multispecies coalescent model for delimitation of species with low morphological diversity. Species diversity of Megophrys is severely underappreciated, and species distributions have been misestimated as a result.

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein-coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single-copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next-generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p-distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.

Phylogenomics reveals rapid, simultaneous diversification of three major clades of Gondwanan frogs at the Cretaceous-Paleogene boundary

Frogs (Anura) are one of the most diverse groups of vertebrates and comprise nearly 90% of living amphibian species. Their worldwide distribution and diverse biology make them well-suited for assessing fundamental questions in evolution, ecology, and conservation. However, despite their scientific importance, the evolutionary history and tempo of frog diversification remain poorly understood. By using a molecular dataset of unprecedented size, including 88-kb characters from 95 nuclear genes of 156 frog species, in conjunction with 20 fossil-based calibrations, our analyses result in the most strongly supported phylogeny of all major frog lineages and provide a timescale of frog evolution that suggests much younger divergence times than suggested by earlier studies. Unexpectedly, our divergence-time analyses show that three species-rich clades (Hyloidea, Microhylidae, and Natatanura), which together comprise ∼88% of extant anuran species, simultaneously underwent rapid diversification at the Cretaceous-Paleogene (K-Pg) boundary (KPB). Moreover, anuran families and subfamilies containing arboreal species originated near or after the KPB. These results suggest that the K-Pg mass extinction may have triggered explosive radiations of frogs by creating new ecological opportunities. This phylogeny also reveals relationships such as Microhylidae being sister to all other ranoid frogs and African continental lineages of Natatanura forming a clade that is sister to a clade of Eurasian, Indian, Melanesian, and Malagasy lineages. Biogeographical analyses suggest that the ancestral area of modern frogs was Africa, and their current distribution is largely associated with the breakup of Pangaea and subsequent Gondwanan fragmentation.

Phylogeny and biogeography of South Chinese brown frogs (Ranidae, Anura)

Few studies have explored the role of Cenozoic tectonic evolution in shaping the patterns and processes of extant animal distributions in and around East Asia. In this study, we selected South Chinese brown frogs as a model to examine the phylogenetic and biogeographical consequences of Miocene tectonic events within South China and its margins. We used mitochondrial and nuclear molecular data to reconstruct phylogenetic interrelationships among Chinese brown frogs using Bayesian and maximum likelihood analyses. The phylogeny results show that there are four main clades of Chinese brown frogs. Excepting the three commonly known Chinese brown frog species groups, R. maoershanensis forms an independent clade nearest to the R. japonica group. Phylogeny and P-distance analyses confirmed R. maoershanensis as a valid species. Among South Chinese brown frogs, there are four subclades associated with four geographical areas: (I) R. maoershanensis; (II) R. japonica; (III) R. chaochiaoensis; and (IV) other species of the R. longicrus species group. Divergence times, estimated using mitochondrial sequences, place the vicariance events among the four subclades in the middle to late Miocene epoch. Our results suggest that (1) South Chinese brown frogs originated due to a vicariance event separating them from the R. chensinensis species group at the time of the Geological movement (~18 million years ago, Ma) in southern Tibet and the Himalayan region; (2) the separation and speciation of R. maoershanensis from the R. japonica group occurred due to the dry climate at approximately 16 Ma; (3) South Chinese brown frogs migrated from South China to Japan at the time (~10.8 Ma) that the global sea-level fell and the East China Sea Shelf Basin was swamp facies, when a land gallery may have formed across the sea to connect the two areas; and (4) R. chaochiaoensis separated from other species of the R. longicrus species group during the uplift of the Tibetan Plateau at approximately 9.5 Ma.

Comprehensive species set revealing the phylogeny and biogeography of Feliformia (Mammalia, Carnivora) based on mitochondrial DNA

Extant Feliformia species are one of the most diverse radiations of Carnivora (~123 species). Despite substantial recent interest in their conservation, diversification, and systematic study, no previous phylogeny contains a comprehensive species set, and no biogeography of this group is available. Here, we present a phylogenetic estimate for Feliformia with a comprehensive species set and establish a historical biogeography based on mitochondrial DNA. Both the Bayesian and maximum likelihood phylogeny for Feliformia are elucidated in our analyses and are strongly consistent with many groups recognized in previous studies. The mitochondrial phylogenetic relationships of Felidae were for the first time successfully reconstructed in our analyses with strong supported. When divergence times and dispersal/vicariance histories were compared with historical sea level changes, four dispersal and six vicariance events were identified. These vicariance events were closely related with global sea level changes. The transgression of sea into the lowland plains between Eurasia and Africa may have caused the vicariance in these regions. A fall in the sea level during late Miocene to Pliocene produced the Bering strait land bridge, which assisted the migration of American Feliformia ancestors from Asia to North America. In contrast with the ‘sweepstakes hypothesis’, our results suggest that the climate cooling during 30–27 Ma assisted Feliformia migration from the African mainland to Madagascar by creating a short-lived ice bridge across the Mozambique Channel. Lineages-through-time plots revealed a large increase in lineages since the Mid-Miocene. During the Mid-Miocene Climatic Optimum, the ecosystems and population of Feliformia rapidly expanded. Subsequent climate cooling catalyzed immigration, speciation, and the extinction of Feliformia.

Simulation-Based Evaluation of Hybridization Network Reconstruction Methods in the Presence of Incomplete Lineage Sorting

Hybridization events generate reticulate species relationships, giving rise to species networks rather than species trees. We report a comparative study of consensus, maximum parsimony, and maximum likelihood methods of species network reconstruction using gene trees simulated assuming a known species history. We evaluate the role of the divergence time between species involved in a hybridization event, the relative contributions of the hybridizing species, and the error in gene tree estimation. When gene tree discordance is mostly due to hybridization and not due to incomplete lineage sorting (ILS), most of the methods can detect even highly skewed hybridization events between highly divergent species. For recent divergences between hybridizing species, when the influence of ILS is sufficiently high, likelihood methods outperform parsimony and consensus methods, which erroneously identify extra hybridizations. The more sophisticated likelihood methods, however, are affected by gene tree errors to a greater extent than are consensus and parsimony.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

A Genome-Scale Investigation of How Sequence, Function, and Tree-Based Gene Properties Influence Phylogenetic Inference

Molecular phylogenetic inference is inherently dependent on choices in both methodology and data. Many insightful studies have shown how choices in methodology, such as the model of sequence evolution or optimality criterion used, can strongly influence inference. In contrast, much less is known about the impact of choices in the properties of the data, typically genes, on phylogenetic inference. We investigated the relationships between 52 gene properties (24 sequence-based, 19 function-based, and 9 tree-based) with each other and with three measures of phylogenetic signal in two assembled data sets of 2,832 yeast and 2,002 mammalian genes. We found that most gene properties, such as evolutionary rate (measured through the percent average of pairwise identity across taxa) and total tree length, were highly correlated with each other. Similarly, several gene properties, such as gene alignment length, Guanine-Cytosine content, and the proportion of tree distance on internal branches divided by relative composition variability (treeness/RCV), were strongly correlated with phylogenetic signal. Analysis of partial correlations between gene properties and phylogenetic signal in which gene evolutionary rate and alignment length were simultaneously controlled, showed similar patterns of correlations, albeit weaker in strength. Examination of the relative importance of each gene property on phylogenetic signal identified gene alignment length, alongside with number of parsimony-informative sites and variable sites, as the most important predictors. Interestingly, the subsets of gene properties that optimally predicted phylogenetic signal differed considerably across our three phylogenetic measures and two data sets; however, gene alignment length and RCV were consistently included as predictors of all three phylogenetic measures in both yeasts and mammals. These results suggest that a handful of sequence-based gene properties are reliable predictors of phylogenetic signal and could be useful in guiding the choice of phylogenetic markers.

The evolution of reproductive diversity in Afrobatrachia: A phylogenetic comparative analysis of an extensive radiation of African frogs

The reproductive modes of anurans (frogs and toads) are the most diverse of terrestrial vertebrates, and a major challenge is identifying selective factors that promote the evolution or retention of reproductive modes across clades. Terrestrialized anuran breeding strategies have evolved repeatedly from the plesiomorphic fully aquatic reproductive mode, a process thought to occur through intermediate reproductive stages. Several selective forces have been proposed for the evolution of terrestrialized reproductive traits, but factors such as water systems and co-evolution with ecomorphologies have not been investigated. We examined these topics in a comparative phylogenetic framework using Afrobatrachian frogs, an ecologically and reproductively diverse clade representing more than half of the total frog diversity found in Africa (∼400 species). We infer direct development has evolved twice independently from terrestrialized reproductive modes involving subterranean or terrestrial oviposition, supporting evolution through intermediate stages. We also detect associations between specific ecomorphologies and oviposition sites, and demonstrate arboreal species exhibit an overall shift toward using lentic water systems for breeding. These results indicate that changes in microhabitat use associated with ecomorphology, which allow access to novel sites for reproductive behavior, oviposition, or larval development, may also promote reproductive mode diversity in anurans.

Tandem amino acid repeats in the green anole (Anolis carolinensis) and other squamates may have a role in increasing genetic variability

Background

Tandem amino acid repeats are characterised by the consecutive recurrence of a single amino acid. They exhibit high rates of length mutations in addition to point mutations and have been proposed to be involved in genetic plasticity. Squamate reptiles (lizards and snakes) diversify in both morphology and physiology. The underlying mechanism is yet to be understood. In a previous phylogenomic analysis of reptiles, the density of tandem repeats in an anole lizard diverged heavily from that of the other reptiles. To gain further insight into the tandem amino acid repeats in squamates, we analysed the repeat content in the green anole (Anolis carolinensis) proteome and compared the amino acid repeats in a large orthologous protein data set from six vertebrates (the Western clawed frog, the green anole, the Chinese softshell turtle, the zebra finch, mouse and human).

Results

Our results revealed that the number of amino acid repeats in the green anole exceeded those found in the other five species studied. Species-only repeats were found in high proportion in the green anole but not in the other five species, suggesting that the green anole had gained many amino acid repeats in either the Anolis or the squamate lineage. Since the amino acid repeat containing genes in the green anole were highly enriched in genes related to transcription and development, an important family of developmental genes, i.e., the Hox family, was further studied in a wide collection of squamates. Abundant amino acid repeats were also observed, implying the general high tolerance of amino acid repeats in squamates. A particular enrichment of amino acid repeats was observed in the central class Hox genes that are known to be responsible for defining cervical to lumbar regions.

Conclusions

Our study suggests that the abundant amino acid repeats in the green anole, and possibly in other squamates, may play a role in increasing the genetic variability, and contribute to the evolutionary diversity of this clade.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2430-y) contains supplementary material, which is available to authorized users.

Seven at one blow: the origin of major lineages of the viviparous Lycian salamanders (Lyciasalamandra Veith and Steinfartz, 2004) was triggered by a single paleo-historic event

The number of tectonic and climatic events that are used to explain speciation processes in the eastern Mediterranean region is low compared to the western Mediterranean. Among them, the emergence of the mid-Aegean trench and the Messinian Salinity Crisis (MSC) often concurred with speciation time estimates that were inferred from molecular data. We here present a dated molecular phylogeny of Lyciasalamandra from Turkey and Greece based on ca. 4500 bp of the mitochondrial genome (3000 bp of three nuclear genes appeared to be completely inconclusive due to their extremely low degree of variation among taxa). Seven major lineages emerged simultaneously from a basal hard polytomy. A scenario that dates this polytomy to 12.3 and 10.2 million years ago, around the final emergence of the mid-Aegean trench, appears to be most plausible. The MSC can be made responsible for first intraspecific divergence events within L. luschani, L. fazilae and L. flavimembris. Further diversification can be explained by Pliocene and Pleistocene glaciations. Based on levels of molecular differentiation we suggest the recently described species L. arikani, L. irfani and L. yehudahi to be treated as subspecies of L. billae.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

An integrative approach to phylogeography: investigating the effects of ancient seaways, climate, and historical geology on multi-locus phylogeographic boundaries of the Arboreal Salamander (Aneides lugubris)

BACKGROUND

Phylogeography is an important tool that can be used to reveal cryptic biodiversity and to better understand the processes that promote lineage diversification. We studied the phylogeographic history of the Arboreal Salamander (Aneides lugubris), a wide-ranging species endemic to the California floristic province. We used multi-locus data to reconstruct the evolutionary history of A. lugubris and to discover the geographic location of major genetic breaks within the species. We also used species distribution modeling and comparative phylogeography to better understand the environmental factors that have shaped the genetic history of A. lugubris.

RESULTS

We found six major mitochondrial clades in A. lugubris. Nuclear loci supported the existence of at least three genetically distinct groups, corresponding to populations north of the San Francisco Bay and in the Sierra Nevada, in the Santa Cruz Mountains, and in the central coast and southern California. All of the genetic breaks in mitochondrial and nuclear loci corresponded to regions where historical barriers to dispersal have been observed in other species. Geologic or water barriers likely were the most important factors restricting gene flow among clades. Climatic unsuitability during glacial maximum may have contributed to the isolation of the mitochondrial clades in the central coast and southern California. A projection of our species distribution model to a future scenario with a moderate amount of climate change suggests that most of the range of A. lugubris will remain climatically suitable, but climatic conditions in the Sierra Nevada and low elevation areas in Southern California are likely to deteriorate.

CONCLUSIONS

Aneides lugubris contains substantial cryptic genetic diversity as a result of historical isolation of populations. At least two (and perhaps three) evolutionarily significant units in A. lugubris merit protection; all six mitochondrial clades should be considered as management units within the species.

Enlarged Multilocus Data set Provides Surprisingly Younger Time of Origin for the Plethodontidae, the Largest Family of Salamanders

Deep phylogenetic relationships of the largest salamander family Plethodontidae have been difficult to resolve, probably reflecting a rapid diversification early in their evolutionary history. Here, data from 50 independent nuclear markers (total 48,582 bp) are used to reconstruct the phylogeny and divergence times for plethodontid salamanders, using both concatenation and coalescence-based species tree analyses. Our results robustly resolve the position of the enigmatic eastern North American four-toed salamander (Hemidactylium) as the sister taxon of Batrachoseps + Tribe Bolitoglossini, thus settling a long-standing question. Furthermore, we statistically reject sister taxon status of Karsenia and Hydromantes, the only plethodontids to occur outside the Americas, leading us to new biogeographic hypotheses. Contrary to previous long-standing arguments that plethodontid salamanders are an old lineage originating in the Cretaceous (more than 90 Ma), our analyses lead to the hypothesis that these salamanders are much younger, arising close to the K-T boundary (~66 Ma). These time estimates are highly stable using alternative calibration schemes and dating methods. Our data simulation highlights the potential risk of making strong arguments about phylogenetic timing based on inferences from a handful of nuclear genes, a common practice. Based on the newly obtained timetree and ancestral area reconstruction results, we argue that (i) the classic "Out of Appalachia" hypothesis of plethodontid origins is problematic; (ii) the common ancestor of extant plethodontids may have originated in northwestern North America in the early Paleocene; (iii) origins of Eurasian plethodontids likely result from two separate dispersal events from western North America via Beringia in the late Eocene (~42 Ma) and the early Miocene (~23 Ma), respectively.