The development of three long universal nuclear protein-coding locus markers and their application to osteichthyan phylogenetics with nested PCR

Unexpected diversity and co-occurrence of phytotelmic frogs (Guibemantis) around Andasibe, one of the most intensively surveyed amphibian hotspots of Madagascar, and descriptions of three new species

The area around the Malagasy village of Andasibe, which includes Analamazaotra-Mantadia National Park as well as other protected areas, is characterized by very species-rich and well-studied communities of animals and plants, but new species are still regularly discovered. Three species of phytotelmic frogs of the subgenus Pandanusicola in the genus Guibemantis are known from this area, G. flavobrunneus, G. liber, and G. pulcher. Further Pandanusicola frogs from this area have been provisionally assigned to G. bicalcaratus or G. albolineatus, pending detailed taxonomic review. During preliminary exploration of the ecology of these specialized frogs that live and reproduce in the leaf axils of Pandanus screw pines, we noticed the syntopic presence of two differently colored and differently sized Pandanusicola in Andasibe that could not be unambiguously assigned to any known species. A genetic screening revealed that these correspond to yet two further species in the area. Based on our data, seven species of Pandanusicola occur in Andasibe and nearby forests: (1) G. liber, the only non-phytotelmic species of the subgenus in the region; (2) G. flavobrunneus which is the largest species and characterized by a diagnostic yellowish brown dorsal pattern; (3) G. pulcher, characterized by translucent-green color with purplish brown spotting not observed in any other species in the area; (4) G. methueni, a brownish species usually lacking contrasted dorsolateral bands that differs from the other species in the area by emitting a characteristic trill-like advertisement call series (rather than clicks or chirps) and according to our data is widespread along Madagascars east coast; as well as three new species: (5) G. ambakoana sp. nov., a brownish species, typically with contrasted incomplete light dorsolateral bands and with single click-like advertisement calls; (6) G. vakoa sp. nov., a species that is equally brownish but lacks contrasted light dorsolateral bands and that has single click-like advertisement calls of very short duration; and (7) G. rianasoa sp. nov., a species that is smaller sized and has less distinct femoral glands than all the others, and emits a short series of soft chirp-like advertisement calls. All these species are genetically highly distinct, with >5% uncorrected pairwise distances in the mitochondrial 16S rRNA gene, and lack of haplotype sharing in two nuclear-encoded genes. The co-occurrence of seven Pandanusicola frogs in a relatively small geographic area is unprecedented in Madagascar and calls for in-depth studies of a possible differentiation in habitat use and life history.

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.

Another step through the crux: a new microendemic rock-dwelling Paroedura (Squamata, Gekkonidae) from south-central Madagascar

Using an integrative taxonomic approach including genetic and morphological data, we formally describe a new microendemic gecko species belonging to the Paroedurabastardi clade, previously referred to as P.bastardi D. We name this taxon currently known from Anja Reserve and Tsaranoro Valley Forest (south-central Madagascar), as P.manongavatosp. nov. The new species differs from other species of the P.bastardi clade by ≥ 12.4% uncorrected p-distance at the mitochondrial 16S rRNA gene and it forms a monophyletic group in the COI mtDNA phylogenetic tree. It lacks haplotype sharing at the nuclear KIAA1239 and CMOS genes with the other species of the same complex, including the syntopic P.rennerae. Given its limited extent of occurrence and high levels of habitat fragmentation linked to forest clearances and fires, we propose the IUCN Red List Category of Critically Endangered, based on the B1ab(iii) criterion. The conservation value of Anja Reserve and Tsaranoro Valley Forest is remarkable. Preserving the remaining deciduous forest habitat is of paramount importance to protect these narrow-range reptile species.

Integrative revision of the Blommersia wittei complex, with description of a new species of frog from western and north-western Madagascar

Frogs of the Blommersia wittei complex are widespread in western and northern Madagascar, and are one of two clades of the family Mantellidae that have colonized the Comoran island of Mayotte. Based on a comprehensive set of DNA sequences of the mitochondrial 16S rRNA gene and the nuclear-encoded RAG1 and SACS genes, integrated with morphological and bioacoustic data, we here analyze the genetic differentiation of populations of this complex across Madagascar. We confirm that a candidate species named B. sp. Ca5 in previous studies represents a genetically well-defined evolutionary lineage distributed over much of western Madagascar, which we describe herein as Blommeria bara sp. nov. based on its molecular and bioacoustic differentiation. Blommersia wittei occurs across northern Madagascar but its type locality Ambanja, at the lower Sambirano river, is very close to the range of another, newly discovered microendemic lineage that was only found at two sites along the upper Sambirano river (here named as candidate species B. sp. Ca12). The B. wittei complex thus provides an example of a clade of closely related Malagasy frogs that contains species widespread over hundreds of kilometers, as well as extreme microendemics. For a full resolution of this species complex, more data need to be collected on the geographical contact among these two lineages, on the morphology and bioacoustics of B. sp. Ca12, and on the north-eastern populations of B. wittei at Sambava, which are weakly differentiated in mitochondrial genes but differ in bioacoustics and possibly in the extent of foot webbing.

A new species of Pandanus-dwelling frog from northern Madagascar related to Guibemantis pulcher

Populations of phytotelmic frogs from northern Madagascar assigned to Guibemantis (Pandanusicola) pulcher are known to differ genetically from populations further south in the eastern rainforest belt of the island, but to date, their status has not been analyzed in depth. We combined molecular genetic data with an examination of color pattern to clarify the taxonomy of these frogs. DNA sequences of both mitochondrial and nuclear-encoded genes were consistently differentiated between the northern populations and those occurring further south. Uncorrected pairwise distance in the 16S rRNA gene was 3.7‒4.3% and thus at a level usually characterizing distinct frog species in Madagascar. Furthermore, the northern specimens were characterized by more and smaller purplish-brown spots on their green dorsal surface, and a less distinct brown patch on the flanks. Although fully conclusive evidence for the species status of the northern lineage from bioacoustic differences, sympatric occurrence or narrow hybrid zone is currently lacking, such species-level distinctness is currently the most likely hypothesis. We therefore name the northern populations as Guibemantis (Pandanusicola) pulcherrimus sp. nov. The new species is known from Makira (type locality) and Bemanevika, and specimens morphologically assignable to this taxon have also been recorded from Masoala, Marojejy and Anjanaharibe-Sud.

The Andaman day gecko paradox: an ancient endemic without pronounced phylogeographic structure

The Andaman day gecko (Phelsuma andamanensis) is endemic to the Andaman Archipelago, located ~ 6000 km away from Madagascar where the genus Phelsuma likely evolved. We complemented existing phylogenetic data with additional markers to show that this species consistently branches off early in the evolution of the genus Phelsuma, and this early origin led us to hypothesize that island populations within the Andaman Archipelago could have further diversified. We sampled the Andaman day gecko from all major islands in the Andamans, developed new microsatellite markers and amplified mitochondrial markers to study population diversification. We detected high allelic diversity in microsatellites, but surprisingly poor geographical structuring. This study demonstrates that the Andaman day gecko has a panmictic population (K = 1), but with weak signals for two clusters that we name ‘North’ (North Andaman, Middle Andaman, Interview, Baratang, Neil, and Long Islands) and ‘South’ (Havelock, South Andaman, Little Andaman Islands). The mitochondrial COI gene uncovered wide haplotype sharing across islands with the presence of several private haplotypes (except for the Little Andaman Island, which only had an exclusive private haplotype) signalling ongoing admixture. This species differs from two other Andaman endemic geckos for which we provide comparative mitochondrial data, where haplotypes show a distinct phylogeographic structure. Testing population history scenarios for the Andaman day gecko using Approximate Bayesian Computation (ABC) supports two possible scenarios but fails to tease apart whether admixture or divergence produced the two weak clusters. Both scenarios agree that admixture and/or divergence prior to the onset of the last glacial maximum shaped the genetic diversity and structure detected in this study. ABC supports population expansion, possibly explained by anthropogenic food subsidies via plantations of cash crops, potentially coupled with human mediated dispersal resulting in the observed panmictic population. The Andaman day gecko may thus be a rare example of an island endemic reptile benefiting from habitat modification and increased movement in its native range.

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.

An integrative DNA barcoding framework of ladybird beetles (Coleoptera: Coccinellidae)

Even though ladybirds are well known as economically important biological control agents, an integrative framework of DNA barcoding research was not available for the family so far. We designed and present a set of efficient mini-barcoding primers to recover full DNA barcoding sequences for Coccinellidae, even for specimens collected 40 years ago. Based on these mini-barcoding primers, we obtained 104 full DNA barcode sequences for 104 species of Coccinellidae, in which 101 barcodes were newly reported for the first time. We also downloaded 870 COI barcode sequences (658 bp) from GenBank and BOLD database, belonging to 108 species within 46 genera, to assess the optimum genetic distance threshold and compare four methods of species delimitation (GMYC, bPTP, BIN and ABGD) to determine the most accurate approach for the family. The results suggested the existence of a 'barcode gap' and that 3% is likely an appropriate genetic distance threshold to delimit species of Coccinellidae using DNA barcodes. Species delimitation analyses confirm ABGD as an accurate and efficient approach, more suitable than the other three methods. Our research provides an integrative framework for DNA barcoding and descriptions of new taxa in Coccinellidae. Our results enrich DNA barcoding public reference libraries, including data for Chinese coccinellids. This will facilitate taxonomic identification and biodiversity monitoring of ladybirds using metabarcoding.

The First Transcriptome Assembly of Yenyuan Stream Salamander (Batrachuperus yenyuanensis) Provides Novel Insights into Its Molecular Evolution

The Yenyuan stream salamander (Batrachuperus yenyuanensis) has been previously evaluated with regards to phylogeny, population genetics, and hematology, but genomic information is sparse due to the giant genome size of salamanders which contain highly repetitive sequences, thus resulting in the lack of a complete reference genome. This study evaluates the encoding genetic sequences and provides the first transcriptome assembly of Yenyuan stream salamander based on mixed samples from the liver, spermary, muscle and spleen tissues. Using this transcriptome assembly and available encoding sequences from other vertebrates, the gene families, phylogenetic status, and species divergence time were compared or estimated. A total of 13,750 encoding sequences were successfully obtained from the transcriptome assembly of Yenyuan stream salamander, estimated to contain 40.1% of the unigenes represented in tetrapod databases. A total of 88.79% of these genes could be annotated to a biological function by current databases. Through gene family clustering, we found multiple possible isoforms of the Scribble gene-whose function is related to regeneration-based on sequence similarity. Meanwhile, we constructed a robust phylogenetic tree based on 56 single-copy orthologues, which indicates that based on phylogenetic position, the Yenyuan stream salamander presents the closest relationship with the Chinese giant salamander (Andrias davidianus) of the investigated vertebrates. Based on the fossil-calibrated phylogeny, we estimated that the lineage divergence between the ancestral Yenyuan stream salamander and the Chinese giant salamander may have occurred during the Cretaceous period (~78.4 million years ago). In conclusion, this study not only provides a candidate gene that is valuable for exploring the remarkable capacity of regeneration in the future, but also gives an interesting insight into the understanding of Yenyuan stream salamander by this first transcriptome assembly.

Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Reconstructing reliable timescales for species evolution is an important and indispensable goal of modern biogeography. However, many factors influence the estimation of divergence times, and uncertainty in the inferred time trees remains a major issue that is often insufficiently acknowledged. We here focus on a fundamental problem of time tree analysis: the combination of slow-evolving (nuclear DNA) and fast-evolving (mitochondrial DNA) markers in a single time tree. Both markers differ in their suitability to infer divergences at different time scales (the 'genome-timescale-dilemma'). However, strategies to infer shallow and deep divergences in a single time tree have rarely been compared empirically. Using Mediterranean amphibians as model system that is exceptional in its geographic and taxonomic completeness of available genetic information, we analyze 202 lineages of western Palearctic amphibians across the entire Mediterranean region. We compiled data of four nuclear and five mitochondrial genes and used twelve fossil calibration points widely acknowledged for amphibian evolution. We reconstruct time trees for an extensive lineage-level data set and compare the performances of the different trees: the first tree is based on primary fossil calibration and mitochondrial DNA, while the second tree is based on a combination of primary fossil and on secondary calibrations taken from a nuclear tree using mitochondrial DNA (two-step protocol). Focusing on a set of nodes that are most likely explained by vicariance, we statistically compare the reconstructed alternative time trees by applying a biogeographical plausibility test. Our two-step protocol outperformed the alternative approach in terms of spatial and temporal plausibility. It allows us to infer scenarios for Mediterranean amphibian evolution in eight geographic provinces. We identified several tectonic and climatic events explaining the majority of Mediterranean amphibian divergences, with Plio-Pleistocene climatic fluctuations being the dominant driver for intrageneric evolution. However, often more than one event could be invoked for a specific split. We give recommendations for the use of secondary calibrations in future molecular clock analyses at the community level.

The changing views on the evolutionary relationships of extant Salamandridae (Amphibia: Urodela)

The phylogenetic relationships among members of the family Salamandridae have been repeatedly investigated over the last 90 years, with changing character and taxon sampling. We review the changing composition and the phylogenetic position of salamandrid genera and species groups and add a new phylogeny based exclusively on sequences of nuclear genes. Salamandrina often changed its position depending on the characters used. It was included several times in a clade together with the primitive newts (Echinotriton, Pleurodeles, Tylototriton) due to their seemingly ancestral morphology. The latter were often inferred as a monophyletic clade. Respective monophyly was almost consistently established in all molecular studies for true salamanders (Chioglossa, Lyciasalamandra, Mertensiella, Salamandra), modern Asian newts (Cynops, Laotriton, Pachytriton, Paramesotriton) and modern New World newts (Notophthalmus, Taricha). Reciprocal non-monophyly has been established through molecular studies for the European mountain newts (Calotriton, Euproctus) and the modern European newts (Ichthyosaura, Lissotriton, Neurergus, Ommatotriton, Triturus) since Calotriton was identified as the sister lineage of Triturus. In pre-molecular studies, their respective monophyly had almost always been assumed, mainly because a complex courtship behaviour shared by their respective members. Our nuclear tree is nearly identical to a mito-genomic tree, with all but one node being highly supported. The major difference concerns the position of Calotriton, which is no longer nested within the modern European newts. This has implications for the evolution of courtship behaviour of European newts. Within modern European newts, Ichthyosaura and Lissotriton changed their position compared to the mito-genomic tree. Previous molecular trees based on seemingly large nuclear data sets, but analysed together with mitochondrial data, did not reveal monophyly of modern European newts since taxon sampling and nuclear gene coverage was too poor to obtain conclusive results. We therefore conclude that mitochondrial and nuclear data should be analysed on their own.

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.

Expression and characterization of the zebrafish orthologue of the human FOLR1 gene during embryogenesis

It has been well established that many types of rapidly dividing normal and diseased cells require an increased amount of folate for DNA replication and repair as well as cellular metabolism. Thus one of folate's cognate receptors, Folate Receptor 1 (FOLR1) is usually up-regulated in rapidly dividing cells, including many types of cancerous tumors. Because zebrafish have become a model organism for understanding conserved vertebrate cellular pathways and human disease, there has been an increased need to identify and elucidate orthologous zebrafish genes that are central to known human maladies. The cells of all early animal embryos go through a phase of rapid division (cleavage) where particular cell cycle checkpoints are skipped until a specification event occurs directing these embryonic stem cells to their fated germ layer cell type. Interestingly, this rapid cell division that ignores cell cycle checkpoints is also observed in many cancers. Developing blastula and tumor cells both require folr1 expression to obtain folate. In this report we have identified the expression pattern of the zebrafish gene zgc:165502, located on chromosome 15. Using computational and comparative methods and molecular biology techniques such as reverse transcription polymerase chain reaction (RT-PCR) and whole mount in situ hybridization (WISH) during embryogenesis, we demonstrate that zgc:165502 is the zebrafish orthologue of the human FOLR1 gene. Understanding when and where FOLR1 orthologues are expressed in different biomedical model organisms such as the zebrafish will help researchers design better experiments to study the endogenous FOLR1 activity.

Stem caecilian from the Triassic of Colorado sheds light on the origins of Lissamphibia

The origin of the limbless caecilians remains a lasting question in vertebrate evolution. Molecular phylogenies and morphology support that caecilians are the sister taxon of batrachians (frogs and salamanders), from which they diverged no later than the early Permian. Although recent efforts have discovered new, early members of the batrachian lineage, the record of pre-Cretaceous caecilians is limited to a single species, Eocaecilia micropodia The position of Eocaecilia within tetrapod phylogeny is controversial, as it already acquired the specialized morphology that characterizes modern caecilians by the Jurassic. Here, we report on a small amphibian from the Upper Triassic of Colorado, United States, with a mélange of caecilian synapomorphies and general lissamphibian plesiomorphies. We evaluated its relationships by designing an inclusive phylogenetic analysis that broadly incorporates definitive members of the modern lissamphibian orders and a diversity of extinct temnospondyl amphibians, including stereospondyls. Our results place the taxon confidently within lissamphibians but demonstrate that the diversity of Permian and Triassic stereospondyls also falls within this group. This hypothesis of caecilian origins closes a substantial morphologic and temporal gap and explains the appeal of morphology-based polyphyly hypotheses for the origins of Lissamphibia while reconciling molecular support for the group's monophyly. Stem caecilian morphology reveals a previously unrecognized stepwise acquisition of typical caecilian cranial apomorphies during the Triassic. A major implication is that many Paleozoic total group lissamphibians (i.e., higher temnospondyls, including the stereospondyl subclade) fall within crown Lissamphibia, which must have originated before 315 million years ago.

Expression of Immune-Related Genes during Loach (Misgurnus anguillicaudatus) Embryonic and Early Larval Development

Early life stage mortality in fish is one of the problems faced by loach aquaculture. However, our understanding of immune system in early life stage fish is still incomplete, and the information available is restricted to a few fish species. In the present work, we investigated the expression of immune-related transcripts in loach during early development. In fishes, recombination-activating gene 1 (RAG-1) and sacsin (SACS) have been considered as immunological function. In this study, the expression of the both genes was assessed throughout the early developmental stages of loach using real-time PCR method. maRAG-1 mRNA was first detected in 0 dph, observed the increased mostly until 40 dph. Significant expression of maRAG-1 was detected in 0 to 40 dph. These patterns of expression may suggest that the loach start to develop its function after hatching. On the other hand, maSACS was detected in unfertilized oocyte to molura stages and 0 to 40 dph. maSACS mRNA transcripts were detected in unfertilized oocytes, suggesting that they are maternally transferred.

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.

Parallel tagged amplicon sequencing of relatively long PCR products using the Illumina HiSeq platform and transcriptome assembly

In phylogenetics and population genetics, a large number of loci are often needed to accurately resolve species relationships. Normally, loci are enriched by PCR and sequenced by Sanger sequencing, which is expensive when the number of amplicons is large. Next-generation sequencing (NGS) techniques are increasingly used for parallel amplicon sequencing, which reduces sequencing costs tremendously, but has not reduced preparation costs very much. Moreover, for most current NGS methods, amplicons need to be purified and quantified before sequencing and their lengths are also restricted (normally <700 bp). Here, we describe an approach to sequence pooled amplicons of any length using the Illumina platform. Using this method, amplicons are pooled at equal volume rather than at equal concentration, thus eliminating the laborious purification and quantification steps. We then shear the pooled amplicons, repair the ends, add sample identifying linkers and pool multiple samples prior to Illumina library preparation. Data are then assembled using the transcriptome assembly program trinity, which is optimized to deal with templates of highly varying quantities. We demonstrated the utility of our approach by recovering 93.5% of the target amplicons (size up to 1650 bp) in full length for a 16 taxa × 101 loci project, using ~2.0 GB of Illumina HiSeq paired-end 90-bp data. Overall, we validate a rapid, cost-effective and scalable approach to sequence a large number of targeted loci from a large number of samples that is particularly suitable for both phylogenetics and population genetics studies that require a modest scale of data.

Impact of gene molecular evolution on phylogenetic reconstruction: a case study in the rosids (Superorder Rosanae, Angiosperms)

Rate of substitution of genomic regions is among the most debated intrinsic features that impact phylogenetic informativeness. However, this variable is also coupled with rates of nonsynonymous substitutions that underscore the nature and degree of selection on the selected genes. To empirically address these variables, we constructed four completely overlapping data sets of plastid matK, atpB, rbcL, and mitochondrial matR genes and used the rosid lineage (angiosperms) as a working platform. The genes differ in combinations of overall rates of nucleotide and amino acid substitutions. Tree robustness, homoplasy, accuracy in contrast to a reference tree, and phylogenetic informativeness are evaluated. The rapidly evolving/unconstrained matK faired best, whereas remaining genes varied in degrees of contribution to rosid phylogenetics across the lineage's 108 million years evolutionary history. Phylogenetic accuracy was low with the slowly evolving/unconstrained matR despite least amount of homoplasy. Third codon positions contributed the highest amount of parsimony informative sites, resolution and informativeness, but magnitude varied with gene mode of evolution. These findings are in clear contrast with the views that rapidly evolving regions and the 3rd codon position have inevitable negative impact on phylogenetic reconstruction at deep historic level due to accumulation of multiple hits and subsequent elevation in homoplasy and saturation. Relaxed evolutionary constraint in rapidly evolving genes distributes substitutions across codon positions, an evolutionary mode expected to reduce the frequency of multiple hits. These findings should be tested at deeper evolutionary histories.

Development, validation and high-throughput analysis of sequence markers in nonmodel species

DNA sequences derived from multiple regions of the nuclear genome are essential for historical inferences in the fields of phylogeography and phylogenetics. The appropriate markers should be single-copy, variable, easy to amplify from multiple samples and easy to sequence using high-throughput technologies. This may be difficult to achieve for species lacking sequenced genomes and particularly challenging for species possessing large genomes, which consist mostly of repetitive sequences. Here, we present a cost-effective, broadly applicable framework for designing, validating and high-throughput sequencing of multiple markers in nonmodel species without sequenced genomes. We demonstrate its utility in two closely related species of newts, representatives of urodeles, a vertebrate group characterized by large genomes. We show that over 80 markers, c. 600 bp each, developed mainly from 3' untranslated transcript regions (3'UTR) may be effectively multiplexed and sequenced. Data are further processed using standard, freely available bioinformatic tools, producing phase-resolved sequences. The approach does not require barcoded PCR primers, and the cost of library preparation is independent of the number of markers investigated. We hope that this approach will be of broad interest for researchers working at the interface of population genetics and phylogenetics, exploring deep intraspecific genetic structure, species boundaries and phylogeographies of closely related species.

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

Resolving difficult nodes for any part of the vertebrate tree of life often requires analyzing a large number of loci. Developing molecular markers that are workable for the groups of interest is often a bottleneck in phylogenetic research. Here, on the basis of a nested polymerase chain reaction (PCR) strategy, we present a universal toolkit including 102 nuclear protein-coding locus (NPCL) markers for vertebrate phylogenomics. The 102 NPCL markers have a broad range of evolutionary rates, which makes them useful for a wide range of time depths. The new NPCL toolkit has three important advantages compared with all previously developed NPCL sets: 1) the kit is universally applicable across vertebrates, with a PCR success rate of 94.6% in 16 widely divergent tested vertebrate species; 2) more than 90% of PCR reactions produce strong and single bands of the expected sizes that can be directly sequenced; and 3) all cleanup PCR reactions can be sequenced with only two specific universal primers. To test its actual phylogenetic utility, 30 NPCLs from this toolkit were used to address the higher level relationships of living salamanders. Of the 639 target PCR reactions performed on 19 salamanders and several outgroup species, 632 (98.9%) were successful, and 602 (94.1%) were directly sequenced. Concatenation and species-tree analyses on this 30-locus data set produced a fully resolved phylogeny and showed that Cryptobranchoidea (Cryptobranchidae + Hynobiidae) branches first within the salamander tree, followed by Sirenidae. Our experimental tests and our demonstration for a particular case show that our NPCL toolkit is a highly reliable, fast, and cost-effective approach for vertebrate phylogenomic studies and thus has the potential to accelerate the completion of many parts of the vertebrate tree of life.

Identifying localized biases in large datasets: a case study using the avian tree of life

Large-scale multi-locus studies have become common in molecular phylogenetics, with new studies continually adding to previous datasets in an effort to fully resolve the tree of life. Total evidence analyses that combine existing data with newly collected data are expected to increase the power of phylogenetic analyses to resolve difficult relationships. However, they might be subject to localized biases, with one or a few loci having a strong and potentially misleading influence upon the results. To examine this possibility we combined a newly collected 31-locus dataset that includes representatives of all major avian lineages with a published dataset of 19 loci that has a comparable number of sites (Hackett et al., 2008. Science 320, 1763-1768). This allowed us to explore the advantages of conducting total evidence analyses, and to determine whether it was also important to analyze new datasets independent of published ones. The total evidence analysis yielded results very similar to the published results, with only slightly increased support at a few nodes. However, analyzing the 31- and 19-locus datasets separately highlighted several differences. Two clades received strong support in the published dataset and total evidence analysis, but the support appeared to reflect bias at a single locus (β-fibrinogen [FGB]). The signal in FGB that supported these relationships was sufficient to result in their recovery with bootstrap support, even when combined with 49 loci lacking that signal. FGB did not appear to have a substantial impact upon the results of species tree methods, but another locus (brain-derived neurotrophic factor [BDNF]) did have an impact upon those analyses. These results demonstrated that localized biases can influence large-scale phylogenetic analyses but they also indicated that considering independent evidence and exploring multiple analytical approaches could reveal them.