Resolving arthropod phylogeny: exploring phylogenetic signal within 41 kb of protein-coding nuclear gene sequence

Taxonomic revision of the Taractrocera ilia (Waterhouse) complex (Lepidoptera : Hesperiidae) from north-western Australia and mainland New Guinea based on morphological and molecular data

The taxonomy of the Taractrocera ilia complex is revised based on evidence from multiple sources, including phenotypic characters (wing colour pattern elements, morphology of adult and immature stages), ecology and molecular data. Three species are recognised in the complex: T. psammopetra Braby, sp. nov. is restricted to north-western Australia where it has a disjunct distribution (Kimberley and the central Arnhem Land plateau); T. ilia (Waterhouse, 1932) sensu stricto is endemic to the central Arnhem Land plateau where it is narrowly sympatric with T. psammopetra; while T. beta Evans, 1934, stat. nov. is allopatric and endemic to mainland New Guinea. Each species is illustrated and diagnosed based on comparative morphology of the adult stage, including androconia and male and female genitalia, and information on type material, biology, distribution and habitat is summarised and discussed. A phylogenetic hypothesis of the Australian species of Taractrocera based on combined analysis of molecular data (2.6 kb: mitochondrial COI and nuclear CAD, MDH, wingless) revealed that T. psammopetra is monophyletic and sister to T. ilia, with an overall mean divergence of 3.4% between the two species; the genus Taractrocera, however, emerged as polyphyletic. For T. psammopetra and T. ilia, the life history, including details of larval food plant specialisation and morphology of the immature stages, is described for the first time. Larvae of T. psammopetra and T. ilia specialise on ‘resurrection’ grasses, Micraira spp. (Poaceae), which typically grow as moss-like, mat-forming pioneer plants on shallow rock surfaces comprising open sandstone pavements or on rock ledges and under rock overhangs of sandstone cliffs. These food plants, together with the butterflies that they support, are a characteristic element of the sandstone plateaux of the Kimberley and central Arnhem Land, and they exemplify patterns of diversification, endemism and adaptation to these unique environments within the Australian monsoon tropics. Facultative larval diapause is recorded for both T. psammopetra and T. ilia in which larvae (all instars, but particularly the final instar) may suspend feeding for up to six months during the dry season when the food plants typically dehydrate. It is hypothesised that larval diapause may be the most flexible strategy to ameliorate seasonal adversity of the dry season and the onset, duration and inter-annual cycle of monsoon rainfall during the wet season, which is highly variable and unpredictable, given that Micraira can revive rapidly from a desiccated state after rainfall. http://zoobank.org/urn:lsid:zoobank.org:pub:51438A5E-048C-4C72-AD22-3FD176ECBA55

Phylotranscriptomics: saturated third codon positions radically influence the estimation of trees based on next-gen data

Recent advancements in molecular sequencing techniques have led to a surge in the number of phylogenetic studies that incorporate large amounts of genetic data. We test the assumption that analyzing large number of genes will lead to improvements in tree resolution and branch support using moths in the superfamily Bombycoidea, a group with some interfamilial relationships that have been difficult to resolve. Specifically, we use a next-gen data set that included 19 taxa and 938 genes (∼1.2M bp) to examine how codon position and saturation might influence resolution and node support among three key families. Maximum likelihood, parsimony, and species tree analysis using gene tree parsimony, on different nucleotide and amino acid data sets, resulted in largely congruent topologies with high bootstrap support compared with prior studies that included fewer loci. However, for a few shallow nodes, nucleotide and amino acid data provided high support for conflicting relationships. The third codon position was saturated and phylogenetic analysis of this position alone supported a completely different, potentially misleading sister group relationship. We used the program RADICAL to assess the number of genes needed to fix some of these difficult nodes. One such node originally needed a total of 850 genes but only required 250 when synonymous signal was removed. Our study shows that, in order to effectively use next-gen data to correctly resolve difficult phylogenetic relationships, it is necessary to assess the effects of synonymous substitutions and third codon positions.

The ovary structure and oogenesis in the basal crustaceans and hexapods. Possible phylogenetic significance

Recent large-scale phylogenetic analyses of exclusively molecular or combined molecular and morphological characters support a close relationship between Crustacea and Hexapoda. The growing consensus on this phylogenetic link is reflected in uniting both taxa under the name Pancrustacea or Tetraconata. Several recent molecular phylogenies have also indicated that the monophyletic hexapods should be nested within paraphyletic crustaceans. However, it is still contentious exactly which crustacean taxon is the sister group to Hexapoda. Among the favored candidates are Branchiopoda, Malacostraca, Remipedia and Xenocarida (Remipedia + Cephalocarida). In this context, we review morphological and ultrastructural features of the ovary architecture and oogenesis in these crustacean groups in search of traits potentially suitable for phylogenetic considerations. We have identified a suite of morphological characters which may prove useful in further comparative studies.

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.

Phylogenetic position of Myriapoda revealed by 454 transcriptome sequencing

Myriapods had been considered closely allied to hexapods (insects and relatives). However, analyses of molecular sequence data have consistently placed Myriapoda either as a sister group of Pancrustacea, comprising crustaceans and hexapods, and thereby supporting the monophyly of Mandibulata, or retrieved Myriapoda as a sister group of Chelicerata (spiders, ticks, mites and allies). In addition, the relationships among the four myriapod groups (Pauropoda, Symphyla, Diplopoda, Chilopoda) are unclear. To resolve the phylogeny of myriapods and their relationship to other main arthropod groups, we collected transcriptome data from the symphylan Symphylella vulgaris, the centipedes Lithobius forficatus and Scolopendra dehaani, and the millipedes Polyxenus lagurus, Glomeris pustulata and Polydesmus angustus by 454 sequencing. We concatenated a multiple sequence alignment that contained 1550 orthologous single copy genes (1,109,847 amino acid positions) from 55 euarthropod and 14 outgroup taxa. The final selected alignment included 181 genes and 37,425 amino acid positions from 55 taxa, with eight myriapods and 33 other euarthropods. Bayesian analyses robustly recovered monophyletic Mandibulata, Pancrustacea and Myriapoda. Most analyses support a sister group relationship of Symphyla in respect to a clade comprising Chilopoda and Diplopoda. Inclusion of additional sequence data from nine myriapod species resulted in an alignment with poor data density, but broader taxon average. With this dataset we inferred Diplopoda+Pauropoda as closest relatives (i.e., Dignatha) and recovered monophyletic Helminthomorpha. Molecular clock calculations suggest an early Cambrian emergence of Myriapoda ∼513 million years ago and a late Cambrian divergence of myriapod classes. This implies a marine origin of the myriapods and independent terrestrialization events during myriapod evolution.

Evolutionary implications of dipluran hexamerins

Hexamerin, as a member of the highly conserved arthropod hemocyanin superfamily, has been shown to be a good marker for the phylogenetic study of insects. However, few studies have been conducted on hexamerins in basal hexapods. The first Diplura hexamerin CspHex1 was reported only recently (Pick and Burmester, 2009). Remarkably, CspHex1 was suggested to have evolved from hexapod hemocyanin subunit type 2, which is very different from all insect hexamerins originated from hexapod hemocyanin subunit type 1. Does this finding suggest double or even multiple origins of hexamerins in Hexapoda? To find more evidence on the evolution of dipluran hexamerins, eight putative hexamerin gene sequences were obtained from three dipluran species, as were three hemocyanin genes from two collembolan species. Unexpectedly, after adding the new sequences into the phylogenetic analyses, all dipluran hexamerins including CspHex1 grouped together and as sister to the insect hexamerins, with high likelihood and Bayesian support. Our analysis supports a single origin of the hexamerins in Hexapoda, and suggests the close relationship between Diplura and Insecta. In addition, our study indicates that a relatively comprehensive taxa sampling is essential to solve some problems in phylogenetic reconstruction.

Molecular phylogeny of Myriapoda provides insights into evolutionary patterns of the mode in post-embryonic development

Myriapoda, a subphylum of Arthropoda, comprises four classes, Chilopoda, Diplopoda, Pauropoda, and Symphyla. While recent molecular evidence has shown that Myriapoda is monophyletic, the internal phylogeny, which is pivotal for understanding the evolutionary history of myriapods, remains unresolved. Here we report the results of phylogenetic analyses and estimations of divergence time and ancestral state of myriapods. Phylogenetic analyses were performed based on three nuclear protein-coding genes determined from 19 myriapods representing the four classes (17 orders) and 11 outgroup species. The results revealed that Symphyla whose phylogenetic position has long been debated is the sister lineage to all other myriapods, and that the interordinal relationships within classes were consistent with traditional classifications. Ancestral state estimation based on the tree topology suggests that myriapods evolved from an ancestral state that was characterized by a hemianamorphic mode of post-embryonic development and had a relatively low number of body segments and legs.

Detection of implausible phylogenetic inferences using posterior predictive assessment of model fit

Systematic phylogenetic error caused by the simplifying assumptions made in models of molecular evolution may be impossible to avoid entirely when attempting to model evolution across massive, diverse data sets. However, not all deficiencies of inference models result in unreliable phylogenetic estimates. The field of phylogenetics lacks a direct method to identify cases where model specification adversely affects inferences. Posterior predictive simulation is a flexible and intuitive approach for assessing goodness-of-fit of the assumed model and priors in a Bayesian phylogenetic analysis. Here, I propose new test statistics for use in posterior predictive assessment of model fit. These test statistics compare phylogenetic inferences from posterior predictive data sets to inferences from the original data. A simulation study demonstrates the utility of these new statistics. The new tests reject the plausibility of inferred tree lengths or topologies more often when data/model combinations produce biased inferences. I also apply this approach to exemplar empirical data sets, highlighting the value of the novel assessments.

Decisive data sets in phylogenomics: lessons from studies on the phylogenetic relationships of primarily wingless insects

Phylogenetic relationships of the primarily wingless insects are still considered unresolved. Even the most comprehensive phylogenomic studies that addressed this question did not yield congruent results. To get a grip on these problems, we here analyzed the sources of incongruence in these phylogenomic studies by using an extended transcriptome data set. Our analyses showed that unevenly distributed missing data can be severely misleading by inflating node support despite the absence of phylogenetic signal. In consequence, only decisive data sets should be used which exclusively comprise data blocks containing all taxa whose relationships are addressed. Additionally, we used Four-cluster Likelihood Mapping (FcLM) to measure the degree of congruence among genes of a data set, as a measure of support alternative to bootstrap. FcLM showed incongruent signal among genes, which in our case is correlated neither with functional class assignment of these genes nor with model misspecification due to unpartitioned analyses. The herein analyzed data set is the currently largest data set covering primarily wingless insects, but failed to elucidate their interordinal phylogenetic relationships. Although this is unsatisfying from a phylogenetic perspective, we try to show that the analyses of structure and signal within phylogenomic data can protect us from biased phylogenetic inferences due to analytical artifacts.

Selecting informative subsets of sparse supermatrices increases the chance to find correct trees

Background

Character matrices with extensive missing data are frequently used in phylogenomics with potentially detrimental effects on the accuracy and robustness of tree inference. Therefore, many investigators select taxa and genes with high data coverage. Drawbacks of these selections are their exclusive reliance on data coverage without consideration of actual signal in the data which might, thus, not deliver optimal data matrices in terms of potential phylogenetic signal. In order to circumvent this problem, we have developed a heuristics implemented in a software called mare which (1) assesses information content of genes in supermatrices using a measure of potential signal combined with data coverage and (2) reduces supermatrices with a simple hill climbing procedure to submatrices with high total information content. We conducted simulation studies using matrices of 50 taxa × 50 genes with heterogeneous phylogenetic signal among genes and data coverage between 10–30%.

Results

With matrices of 50 taxa × 50 genes with heterogeneous phylogenetic signal among genes and data coverage between 10–30% Maximum Likelihood (ML) tree reconstructions failed to recover correct trees. A selection of a data subset with the herein proposed approach increased the chance to recover correct partial trees more than 10-fold. The selection of data subsets with the herein proposed simple hill climbing procedure performed well either considering the information content or just a simple presence/absence information of genes. We also applied our approach on an empirical data set, addressing questions of vertebrate systematics. With this empirical dataset selecting a data subset with high information content and supporting a tree with high average boostrap support was most successful if information content of genes was considered.

Conclusions

Our analyses of simulated and empirical data demonstrate that sparse supermatrices can be reduced on a formal basis outperforming the usually used simple selections of taxa and genes with high data coverage.

When integration fails: Prokaryote phylogeny and the tree of life

Much is being written these days about integration, its desirability and even its necessity when complex research problems are to be addressed. Seldom, however, do we hear much about the failure of such efforts. Because integration is an ongoing activity rather than a final achievement, and because today's literature about integration consists mostly of manifesto statements rather than precise descriptions, an examination of unsuccessful integration could be illuminating to understand better how it works. This paper will examine the case of prokaryote phylogeny and its apparent failure to achieve integration within broader tree-of-life accounts of evolutionary history (often called 'universal phylogeny'). Despite the fact that integrated databases exist of molecules pertinent to the phylogenetic reconstruction of all lineages of life, and even though the same methods can be used to construct phylogenies wherever the organisms fall on the tree of life, prokaryote phylogeny remains at best only partly integrated within tree-of-life efforts. I will examine why integration does not occur, compare it with integrative practices in animal and other eukaryote phylogeny, and reflect on whether there might be different expectations of what integration should achieve. Finally, I will draw some general conclusions about integration and its function as a 'meta-heuristic' in the normative commitments guiding scientific practice.

Embryonic neurogenesis in Pseudopallene sp. (Arthropoda, Pycnogonida) includes two subsequent phases with similarities to different arthropod groups

BACKGROUND

Studies on early neurogenesis have had considerable impact on the discussion of the phylogenetic relationships of arthropods, having revealed striking similarities and differences between the major lineages. In Hexapoda and crustaceans, neurogenesis involves the neuroblast, a type of neural stem cell. In each hemi-segment, a set of neuroblasts produces neural cells by repeated asymmetrical and interiorly directed divisions. In Euchelicerata and Myriapoda, neurogenesis lacks neural stem cells, featuring instead direct immigration of neural cell groups from fixed sites in the neuroectoderm. Accordingly, neural stem cells were hitherto assumed to be an evolutionary novelty of the Tetraconata (Hexapoda + crustaceans). To further test this hypothesis, we investigated neurogenesis in Pycnogonida, or sea spiders, a group of marine arthropods with close affinities to euchelicerates.

RESULTS

We studied neurogenesis during embryonic development of Pseudopallene sp. (Callipallenidae), using fluorescent histochemical staining and immunolabelling. Embryonic neurogenesis has two phases. The first phase shows notable similarities to euchelicerates and myriapods. These include i) the lack of morphologically different cell types in the neuroectoderm; ii) the formation of transiently identifiable, stereotypically arranged cell internalization sites; iii) immigration of predominantly post-mitotic ganglion cells; and iv) restriction of tangentially oriented cell proliferation to the apical cell layer. However, in the second phase, the formation of a central invagination in each hemi-neuromere is accompanied by the differentiation of apical neural stem cells. The latter grow in size, show high mitotic activity and an asymmetrical division mode. A marked increase of ganglion cell numbers follows their differentiation. Directly basal to the neural stem cells, an additional type of intermediate neural precursor is found.

CONCLUSIONS

Embryonic neurogenesis of Pseudopallene sp. combines features of central nervous system development that have been hitherto described separately in different arthropod taxa. The two-phase character of pycnogonid neurogenesis calls for a thorough reinvestigation of other non-model arthropods over the entire course of neurogenesis. With the currently available data, a common origin of pycnogonid neural stem cells and tetraconate neuroblasts remains unresolved. To acknowledge this, we present two possible scenarios on the evolution of arthropod neurogenesis, whereby Myriapoda play a key role in the resolution of this issue.

Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha

BACKGROUND

Molecular phylogenetic analyses have revealed that Hexapoda and Crustacea form a common clade (the Pancrustacea), which is now widely accepted among zoologists; however, the origin of Hexapoda remains unresolved. The main problems are the unclear relationships among the basal hexapod lineages, Protura (proturans), Collembola (springtails), Diplura (diplurans), and Ectognatha (bristletails, silverfishes, and all winged insects). Mitogenomic analyses have challenged hexapod monophyly and suggested the reciprocal paraphyly of Hexapoda and Crustacea, whereas studies based on nuclear molecular data support the monophyletic origin of hexapods. Additionally, there are significant discrepancies with respect to these issues between the results of morphological and molecular studies. To investigate these problems, we performed phylogenetic analyses of Pancrustacea based on the protein sequences of three orthologous nuclear genes encoding the catalytic subunit of DNA polymerase delta and the largest and second largest subunits of RNA polymerase II from 64 species of arthropods, including representatives of all hexapod orders.

RESULTS

Phylogenetic analyses were conducted based on the inferred amino acid (aa) sequences (~3400 aa in total) of the three genes using the maximum likelihood (ML) method and Bayesian inference. Analyses were also performed with additional datasets generated by excluding long-branch taxa or by using different outgroups. These analyses all yielded essentially the same results. All hexapods were clustered into a common clade, with Branchiopoda as its sister lineage, whereas Crustacea was paraphyletic. Within Hexapoda, the lineages Ectognatha, Palaeoptera, Neoptera, Polyneoptera, and Holometabola were each confirmed to be monophyletic with robust support, but monophyly was not supported for Entognatha (Protura + Collembola + Diplura), Ellipura (Protura + Collembola), or Nonoculata (Protura + Diplura). Instead, our results showed that Protura is the sister lineage to all other hexapods and that Diplura or Diplura + Collembola is closely related to Ectognatha.

CONCLUSION

This is the first study to include all hexapod orders in a phylogenetic analysis using multiple nuclear protein-coding genes to investigate the phylogeny of Hexapoda, with an emphasis on Entognatha. The results strongly support the monophyletic origin of hexapods but reject the monophyly of Entognatha, Ellipura, and Nonoculata. Our results provided the first molecular evidence in support of Protura as the sister group to other hexapods. These findings are expected to provide additional insights into the origin of hexapods and the processes involved in the adaptation of insects to life on land.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Arthropod phylogenetics in light of three novel millipede (myriapoda: diplopoda) mitochondrial genomes with comments on the appropriateness of mitochondrial genome sequence data for inferring deep level relationships

BACKGROUND

Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda.

RESULTS

The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic.

CONCLUSIONS

The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.

Phylogeny of the Paracalanidae Giesbrecht, 1888 (Crustacea: Copepoda: Calanoida)

The Paracalanidae are ecologically-important marine planktonic copepods that occur in the epipelagic zone in temperate and tropical waters. They are often the dominant taxon - in terms of biomass and abundance - in continental shelf regions. As primary consumers, they form a vital link in the pelagic food web between primary producers and higher trophic levels. Despite the ecological importance of the taxon, evolutionary and systematic relationships within the family remain largely unknown. A multigene phylogeny including 24 species, including representatives for all seven genera, was determined based on two nuclear genes, small-subunit (18S) ribosomal RNA and Histone 3 (H3) and one mitochondrial gene, cytochrome c oxidase subunit I (COI). The molecular phylogeny was well supported by Maximum likelihood and Bayesian inference analysis; all genera were found to be monophyletic, except for Paracalanus, which was separated into two distinct clades: the Paracalanus aculeatus group and Paracalanus parvus group. The molecular phylogeny also confirmed previous findings that Mecynocera and Calocalanus are genera of the family Paracalanidae. For comparison, a morphological phylogeny was created for 35 paracalanid species based on 54 morphological characters derived from published descriptions. The morphological phylogeny did not resolve all genera as monophyletic and bootstrap support was not strong. Molecular and morphological phylogenies were not congruent in the positioning of Bestiolina and the Paracalanus species groups, possibly due to the lack of sufficient phylogenetically-informative morphological characters.

Serotonin-immunoreactive neurons in the ventral nerve cord of Remipedia (Crustacea): support for a sister group relationship of Remipedia and Hexapoda?

Background

Remipedia were initially seen as a primitive taxon within Pancrustacea based on characters considered ancestral, such as the homonomously segmented trunk. Meanwhile, several morphological and molecular studies proposed a more derived position of Remipedia within Pancrustacea, including a sister group relationship to Hexapoda. Because of these conflicting hypotheses, fresh data are crucial to contribute new insights into euarthropod phylogeny. The architecture of individually identifiable serotonin-immunoreactive neurons has successfully been used for phylogenetic considerations in Euarthropoda. Here, we identified neurons in three species of Remipedia with an antiserum against serotonin and compared our findings to reconstructed ground patterns in other euarthropod taxa. Additionally, we traced neurite connectivity and neuropil outlines using antisera against acetylated α-tubulin and synapsin.

Results

The ventral nerve cord of Remipedia displays a typical rope-ladder-like arrangement of separate metameric ganglia linked by paired longitudinally projecting connectives. The peripheral projections comprise an intersegmental nerve, consisting of two branches that fuse shortly after exiting the connectives, and the segmental anterior and posterior nerve. The distribution and morphology of serotonin-immunoreactive interneurons in the trunk segments is highly conserved within the remipede species we analyzed, which allows for the reconstruction of a ground pattern: two posterior and one anterior pair of serotonin-immunoreactive neurons that possess a single contralateral projection. Additionally, three pairs of immunoreactive neurons are found in the medial part of each hemiganglion. In one species (Cryptocorynetes haptodiscus), the anterior pair of immunoreactive neurons is missing.

Conclusions

The anatomy of the remipede ventral nerve cord with its separate metameric ganglia mirrors the external morphology of the animal’s trunk. The rope-ladder-like structure and principal architecture of the segmental ganglia in Remipedia corresponds closely to that of other Euarthropoda. A comparison of the serotonin-immunoreactive cell arrangement of Remipedia to reconstructed ground patterns of major euarthropod taxa supports a homology of the anterior and posterior neurons in Pancrustacea. These neurons in Remipedia possess unbranched projections across the midline, pointing towards similarities to the hexapod pattern. Our findings are in line with a growing number of phylogenetic investigations proposing Remipedia to be a rather derived crustacean lineage that perhaps has close affinities to Hexapoda.

Elongation factor-1α, a putative single-copy nuclear gene, has divergent sets of paralogs in an arachnid

Identification of paralogy in candidate nuclear loci is an important prerequisite in phylogenetics and statistical phylogeography, but one that is often overlooked. One marker commonly assumed to be a single-copy gene and claimed to harbor great utility for inferring recent divergences is elongation factor-1α (EF-1α). To test this hypothesis, we systematically cloned EF-1α in three disjunct populations of the harvestman Metasiro americanus. Here we show that EF-1α has a large number of paralogs in this species. The paralogs do not evolve in a concerted manner, and the paralogs diverged prior to the population divergence. Moreover, the paralogs of M. americanus are not comparable to the highly divergent EF-1α paralogs found in bees and spiders, which are easily recognized and separated through the use of specific primers. We demonstrate statistically that our detection of paralogs cannot be attributed to amplification error. The presence of EF-1α paralogs in M. americanus prevents its use in statistical phylogeography, and the presence of out-paralogs argues against its use in phylogenetic inference among recently diverged clades. These data contradict the common assumption that EF-1α is for most or all taxa a single-copy gene, or that it has a small number of paralogs that are homogenized through gene conversion, unequal crossing over, or other processes.

A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies)

Background

Higher-level relationships within the Lepidoptera, and particularly within the species-rich subclade Ditrysia, are generally not well understood, although recent studies have yielded progress. We present the most comprehensive molecular analysis of lepidopteran phylogeny to date, focusing on relationships among superfamilies.

Methodology / Principal Findings

483 taxa spanning 115 of 124 families were sampled for 19 protein-coding nuclear genes, from which maximum likelihood tree estimates and bootstrap percentages were obtained using GARLI. Assessment of heuristic search effectiveness showed that better trees and higher bootstrap percentages probably remain to be discovered even after 1000 or more search replicates, but further search proved impractical even with grid computing. Other analyses explored the effects of sampling nonsynonymous change only versus partitioned and unpartitioned total nucleotide change; deletion of rogue taxa; and compositional heterogeneity. Relationships among the non-ditrysian lineages previously inferred from morphology were largely confirmed, plus some new ones, with strong support. Robust support was also found for divergences among non-apoditrysian lineages of Ditrysia, but only rarely so within Apoditrysia. Paraphyly for Tineoidea is strongly supported by analysis of nonsynonymous-only signal; conflicting, strong support for tineoid monophyly when synonymous signal was added back is shown to result from compositional heterogeneity.

Conclusions / Significance

Support for among-superfamily relationships outside the Apoditrysia is now generally strong. Comparable support is mostly lacking within Apoditrysia, but dramatically increased bootstrap percentages for some nodes after rogue taxon removal, and concordance with other evidence, strongly suggest that our picture of apoditrysian phylogeny is approximately correct. This study highlights the challenge of finding optimal topologies when analyzing hundreds of taxa. It also shows that some nodes get strong support only when analysis is restricted to nonsynonymous change, while total change is necessary for strong support of others. Thus, multiple types of analyses will be necessary to fully resolve lepidopteran phylogeny.

A molecular phylogeny for yponomeutoidea (insecta, Lepidoptera, ditrysia) and its implications for classification, biogeography and the evolution of host plant use

BACKGROUND

Yponomeutoidea, one of the early-diverging lineages of ditrysian Lepidoptera, comprise about 1,800 species worldwide, including notable pests and insect-plant interaction models. Yponomeutoids were one of the earliest lepidopteran clades to evolve external feeding and to extensively colonize herbaceous angiosperms. Despite the group's economic importance, and its value for tracing early lepidopteran evolution, the biodiversity and phylogeny of Yponomeutoidea have been relatively little studied.

METHODOLOGY/PRINCIPAL FINDINGS

Eight nuclear genes (8 kb) were initially sequenced for 86 putative yponomeutoid species, spanning all previously recognized suprageneric groups, and 53 outgroups representing 22 families and 12 superfamilies. Eleven to 19 additional genes, yielding a total of 14.8 to 18.9 kb, were then sampled for a subset of taxa, including 28 yponomeutoids and 43 outgroups. Maximum likelihood analyses were conducted on data sets differing in numbers of genes, matrix completeness, inclusion/weighting of synonymous substitutions, and inclusion/exclusion of "rogue" taxa. Monophyly for Yponomeutoidea was supported very strongly when the 18 "rogue" taxa were excluded, and moderately otherwise. Results from different analyses are highly congruent and relationships within Yponomeutoidea are well supported overall. There is strong support overall for monophyly of families previously recognized on morphological grounds, including Yponomeutidae, Ypsolophidae, Plutellidae, Glyphipterigidae, Argyresthiidae, Attevidae, Praydidae, Heliodinidae, and Bedelliidae. We also assign family rank to Scythropiinae (Scythropiidae stat. rev.), which in our trees are strongly grouped with Bedelliidae, in contrast to all previous proposals. We present a working hypothesis of among-family relationships, and an informal higher classification. Host plant family associations of yponomeutoid subfamilies and families are non-random, but show no trends suggesting parallel phylogenesis. Our analyses suggest that previous characterizations of yponomeutoids as predominantly Holarctic were based on insufficient sampling.

CONCLUSIONS/SIGNIFICANCE

We provide the first robust molecular phylogeny for Yponomeutoidea, together with a revised classification and new insights into their life history evolution and biogeography.

Resolving discrepancy between nucleotides and amino acids in deep-level arthropod phylogenomics: differentiating serine codons in 21-amino-acid models

BACKGROUND

In a previous study of higher-level arthropod phylogeny, analyses of nucleotide sequences from 62 protein-coding nuclear genes for 80 panarthopod species yielded significantly higher bootstrap support for selected nodes than did amino acids. This study investigates the cause of that discrepancy.

METHODOLOGY/PRINCIPAL FINDINGS

The hypothesis is tested that failure to distinguish the serine residues encoded by two disjunct clusters of codons (TCN, AGY) in amino acid analyses leads to this discrepancy. In one test, the two clusters of serine codons (Ser1, Ser2) are conceptually translated as separate amino acids. Analysis of the resulting 21-amino-acid data matrix shows striking increases in bootstrap support, in some cases matching that in nucleotide analyses. In a second approach, nucleotide and 20-amino-acid data sets are artificially altered through targeted deletions, modifications, and replacements, revealing the pivotal contributions of distinct Ser1 and Ser2 codons. We confirm that previous methods of coding nonsynonymous nucleotide change are robust and computationally efficient by introducing two new degeneracy coding methods. We demonstrate for degeneracy coding that neither compositional heterogeneity at the level of nucleotides nor codon usage bias between Ser1 and Ser2 clusters of codons (or their separately coded amino acids) is a major source of non-phylogenetic signal.

CONCLUSIONS

The incongruity in support between amino-acid and nucleotide analyses of the forementioned arthropod data set is resolved by showing that "standard" 20-amino-acid analyses yield lower node support specifically when serine provides crucial signal. Separate coding of Ser1 and Ser2 residues yields support commensurate with that found by degenerated nucleotides, without introducing phylogenetic artifacts. While exclusion of all serine data leads to reduced support for serine-sensitive nodes, these nodes are still recovered in the ML topology, indicating that the enhanced signal from Ser1 and Ser2 is not qualitatively different from that of the other amino acids.

Cap'n'collar differentiates the mandible from the maxilla in the beetle Tribolium castaneum

BACKGROUND

The biting mandible of the arthropods is thought to have evolved in the ancestor of the insects, crustaceans and myriapods: the Mandibulata. A unique origin suggests a common set of developmental genes will be required to pattern the mandible in different arthropods. To date we have functional studies on patterning of the mandibular segment of Drosophila melanogaster showing in particular the effects of the gene cap'n'collar (cnc), however, the dipteran head is far from representative of insects or of more distantly related mandibulates; Drosophila does not even possess a mandibular appendage. To study the development of a more representative insect mandible, we chose the red flour beetle Tribolium castaneum and investigated the function of the Tribolium orthologs of cap'n'collar (Tc-cnc) and the Hox gene Deformed (Tc-Dfd). In order to determine the function of Tc-cnc and Tc-Dfd, transcripts were knocked down by maternal RNA interference (RNAi). The effects of gene knockdown were examined in the developing embryos and larvae. The effect of Tc-cnc and Tc-Dfd knockdown on the expression of other genes was determined by using in situ hybridization on Tribolium embryos.

RESULTS

Our analyses show that Tc-cnc is required for specification of the identity of the mandibular segment of Tribolium and differentiates the mandible from maxillary identity. Loss of Tc-cnc function results in a transformation of the mandible to maxillary identity as well as deletion of the labrum. Tc-Dfd and the Tribolium homolog of proboscipedia (Tc-mxp = maxillopedia), Hox genes that are required to pattern the maxillary appendage, are expressed in a maxilla-like manner in the transformed mandible. Tribolium homologs of paired (Tc-prd) and Distal-less (Tc-Dll) that are expressed in the endites and telopodites of embryonic appendages are also expressed in a maxilla-like manner in the transformed mandible.We also show that Tc-Dfd is required to activate the collar of Tc-cnc expression in the mandibular segment but not the cap expression in the labrum. Tc-Dfd is also required for the activation of Tc-prd in the endites of the mandible and maxillary appendages.

CONCLUSIONS

Tc-cnc is necessary for patterning the mandibular segment of Tribolium. Together, Tc-cnc and Tc-Dfd cooperate to specify mandibular identity, as in Drosophila. Expression patterns of the homologs of cnc and Dfd are conserved in mandibulate arthropods suggesting that the mandible specifying function of cnc is likely to be conserved across the mandibulate arthropods.

Serine codon-usage bias in deep phylogenomics: pancrustacean relationships as a case study

Phylogenomic analyses of ancient relationships are usually performed using amino acid data, but it is unclear whether amino acids or nucleotides should be preferred. With the 2-fold aim of addressing this problem and clarifying pancrustacean relationships, we explored the signals in the 62 protein-coding genes carefully assembled by Regier et al. in 2010. With reference to the pancrustaceans, this data set infers a highly supported nucleotide tree that is substantially different to the corresponding, but poorly supported, amino acid one. We show that the discrepancy between the nucleotide-based and the amino acids-based trees is caused by substitutions within synonymous codon families (especially those of serine-TCN and AGY). We show that different arthropod lineages are differentially biased in their usage of serine, arginine, and leucine synonymous codons, and that the serine bias is correlated with the topology derived from the nucleotides, but not the amino acids. We suggest that a parallel, partially compositionally driven, synonymous codon-usage bias affects the nucleotide topology. As substitutions between serine codon families can proceed through threonine or cysteine intermediates, amino acid data sets might also be affected by the serine codon-usage bias. We suggest that a Dayhoff recoding strategy would partially ameliorate the effects of such bias. Although amino acids provide an alternative hypothesis of pancrustacean relationships, neither the nucleotides nor the amino acids version of this data set seems to bring enough genuine phylogenetic information to robustly resolve the relationships within group, which should still be considered unresolved.

Serotonin immunoreactive interneurons in the brain of the Remipedia: new insights into the phylogenetic affinities of an enigmatic crustacean taxon

BACKGROUND

Remipedia, a group of homonomously segmented, cave-dwelling, eyeless arthropods have been regarded as basal crustaceans in most early morphological and taxonomic studies. However, molecular sequence information together with the discovery of a highly differentiated brain led to a reconsideration of their phylogenetic position. Various conflicting hypotheses have been proposed including the claim for a basal position of Remipedia up to a close relationship with Malacostraca or Hexapoda. To provide new morphological characters that may allow phylogenetic insights, we have analyzed the architecture of the remipede brain in more detail using immunocytochemistry (serotonin, acetylated α-tubulin, synapsin) combined with confocal laser-scanning microscopy and image reconstruction techniques. This approach allows for a comprehensive neuroanatomical comparison with other crustacean and hexapod taxa.

RESULTS

The dominant structures of the brain are the deutocerebral olfactory neuropils, which are linked by the olfactory globular tracts to the protocerebral hemiellipsoid bodies. The olfactory globular tracts form a characteristic chiasm in the center of the brain. In Speleonectes tulumensis, each brain hemisphere contains about 120 serotonin immunoreactive neurons, which are distributed in distinct cell groups supplying fine, profusely branching neurites to 16 neuropilar domains. The olfactory neuropil comprises more than 300 spherical olfactory glomeruli arranged in sublobes. Eight serotonin immunoreactive neurons homogeneously innervate the olfactory glomeruli. In the protocerebrum, serotonin immunoreactivity revealed several structures, which, based on their position and connectivity resemble a central complex comprising a central body, a protocerebral bridge, W-, X-, Y-, Z-tracts, and lateral accessory lobes.

CONCLUSIONS

The brain of Remipedia shows several plesiomorphic features shared with other Mandibulata, such as deutocerebral olfactory neuropils with a glomerular organization, innervations by serotonin immunoreactive interneurons, and connections to protocerebral neuropils. Also, we provided tentative evidence for W-, X-, Y-, Z-tracts in the remipedian central complex like in the brain of Malacostraca, and Hexapoda. Furthermore, Remipedia display several synapomorphies with Malacostraca supporting a sister group relationship between both taxa. These homologies include a chiasm of the olfactory globular tract, which connects the olfactory neuropils with the lateral protocerebrum and the presence of hemiellipsoid bodies. Even though a growing number of molecular investigations unites Remipedia and Cephalocarida, our neuroanatomical comparison does not provide support for such a sister group relationship.

The complete mitochondrial genome of Pauropus longiramus (Myriapoda: Pauropoda): implications on early diversification of the myriapods revealed from comparative analysis

Myriapods are among the earliest arthropods and may have evolved to become part of the terrestrial biota more than 400 million years ago. A noticeable lack of mitochondrial genome data from Pauropoda hampers phylogenetic and evolutionary studies within the subphylum Myriapoda. We sequenced the first complete mitochondrial genome of a microscopic pauropod, Pauropus longiramus (Arthropoda: Myriapoda), and conducted comprehensive mitogenomic analyses across the Myriapoda. The pauropod mitochondrial genome is a circular molecule of 14,487 bp long and contains the entire set of thirty-seven genes. Frequent intergenic overlaps occurred between adjacent tRNAs, and between tRNA and protein-coding genes. This is the first example of a mitochondrial genome with multiple intergenic overlaps and reveals a strategy for arthropods to effectively compact the mitochondrial genome by overlapping and truncating tRNA genes with neighbor genes, instead of only truncating tRNAs. Phylogenetic analyses based on protein-coding genes provide strong evidence that the sister group of Pauropoda is Symphyla. Additionally, approximately unbiased (AU) tests strongly support the Progoneata and confirm the basal position of Chilopoda in Myriapoda. This study provides an estimation of myriapod origins around 555 Ma (95% CI: 444-704 Ma) and this date is comparable with that of the Cambrian explosion and candidate myriapod-like fossils. A new time-scale suggests that deep radiations during early myriapod diversification occurred at least three times, not once as previously proposed. A Carboniferous origin of pauropods is congruent with the idea that these taxa are derived, rather than basal, progoneatans.

A molecular phylogeny for the leaf-roller moths (Lepidoptera: Tortricidae) and its implications for classification and life history evolution

BACKGROUND

Tortricidae, one of the largest families of microlepidopterans, comprise about 10,000 described species worldwide, including important pests, biological control agents and experimental models. Understanding of tortricid phylogeny, the basis for a predictive classification, is currently provisional. We present the first detailed molecular estimate of relationships across the tribes and subfamilies of Tortricidae, assess its concordance with previous morphological evidence, and re-examine postulated evolutionary trends in host plant use and biogeography.

METHODOLOGY/PRINCIPAL FINDINGS

We sequenced up to five nuclear genes (6,633 bp) in each of 52 tortricids spanning all three subfamilies and 19 of the 22 tribes, plus up to 14 additional genes, for a total of 14,826 bp, in 29 of those taxa plus all 14 outgroup taxa. Maximum likelihood analyses yield trees that, within Tortricidae, differ little among data sets and character treatments and are nearly always strongly supported at all levels of divergence. Support for several nodes was greatly increased by the additional 14 genes sequenced in just 29 of 52 tortricids, with no evidence of phylogenetic artifacts from deliberately incomplete gene sampling. There is strong support for the monophyly of Tortricinae and of Olethreutinae, and for grouping of these to the exclusion of Chlidanotinae. Relationships among tribes are robustly resolved in Tortricinae and mostly so in Olethreutinae. Feeding habit (internal versus external) is strongly conserved on the phylogeny. Within Tortricinae, a clade characterized by eggs being deposited in large clusters, in contrast to singly or in small batches, has markedly elevated incidence of polyphagous species. The five earliest-branching tortricid lineages are all species-poor tribes with mainly southern/tropical distributions, consistent with a hypothesized Gondwanan origin for the family.

CONCLUSIONS/SIGNIFICANCE

We present the first robustly supported phylogeny for Tortricidae, and a revised classification in which all of the sampled tribes are now monophyletic.

Phylogenetics of the millipede genus Brachycybe Wood, 1864 (Diplopoda: Platydesmida: Andrognathidae): patterns of deep evolutionary history and recent speciation

The genus Brachycybe Wood is a little known group of millipedes comprising eight nominal species distributed throughout North America, Japan, South Korea, Taiwan, and China. The group's species are relatively morphologically homogenous and have been described primarily on the basis of differences in somatic morphology largely ignoring the often-used characters in millipede taxonomy and systematics - male genitalia (the gonopods). The objectives of this study were to survey male gonopods with the aim of evaluating inter-specific variation, assess existing species boundaries and phylogeny using molecular characters, examine the historical biogeography of the genus, and estimate the timing of lineage divergence using a molecular clock. We surveyed two mitochondrial genes (cytochrome c oxidase I and cytochrome b) and one nuclear protein-coding gene (glutamyl and prolyl-tRNA synthetase) each comprising 952, 746, and 555 aligned base pairs respectively. Phylogenetic inference coupled with an analysis of species delimitation using a generalized mixed Yule coalescent model recovered eight species, two of which were considered to be cryptic. Molecular dating analyses coupled with ancestral range reconstructions indicate that the group is quite ancient (age>50 million years) and its origins are likely traced back geographically to the mountains of California's Sierra Nevada and Coastal Ranges. The genus Brachycybe appears to have expanded its range at least twice out of present day California into eastern North America and at least once into Asia. This study highlights the need for integrative approaches to describe biodiversity and furthers the evidence for cryptic diversity even in groups where genitalia are generally thought to be rapidly evolving, diagnostic features.

Advances in insect phylogeny at the dawn of the postgenomic era

Most species on Earth are insects and thus, understanding their evolutionary relationships is key to understanding the evolution of life. Insect relationships are increasingly well supported, due largely to technological advances in molecular sequencing and phylogenetic computational analysis. In this postgenomic era, insect systematics will be furthered best by integrative methods aimed at hypothesis corroboration from molecular, morphological, and paleontological evidence. This review of the current consensus of insect relationships provides a foundation for comparative study and offers a framework to evaluate incoming genomic evidence. Notable recent phylogenetic successes include the resolution of Holometabola, including the identification of the enigmatic Strepsiptera as a beetle relative and the early divergence of Hymenoptera; the recognition of hexapods as a crustacean lineage within Pancrustacea; and the elucidation of Dictyoptera orders, with termites placed as social cockroaches. Regions of the tree that require further investigation include the earliest winged insects (Palaeoptera) and Polyneoptera (orthopteroid lineages).

Insect phylogenomics: exploring the source of incongruence using new transcriptomic data

The evolution of the diverse insect lineages is one of the most fascinating issues in evolutionary biology. Despite extensive research in this area, the resolution of insect phylogeny especially of interordinal relationships has turned out to be still a great challenge. One of the challenges for insect systematics is the radiation of the polyneopteran lineages with several contradictory and/or unresolved relationships. Here, we provide the first transcriptomic data for three enigmatic polyneopteran orders (Dermaptera, Plecoptera, and Zoraptera) to clarify one of the most debated issues among higher insect systematics. We applied different approaches to generate 3 data sets comprising 78 species and 1,579 clusters of orthologous genes. Using these three matrices, we explored several key mechanistic problems of phylogenetic reconstruction including missing data, matrix selection, gene and taxa number/choice, and the biological function of the genes. Based on the first phylogenomic approach including these three ambiguous polyneopteran orders, we provide here conclusive support for monophyletic Polyneoptera, contesting the hypothesis of Zoraptera + Paraneoptera and Plecoptera + remaining Neoptera. In addition, we employ various approaches to evaluate data quality and highlight problematic nodes within the Insect Tree that still exist despite our phylogenomic approach. We further show how the support for these nodes or alternative hypotheses might depend on the taxon- and/or gene-sampling.

Can deliberately incomplete gene sample augmentation improve a phylogeny estimate for the advanced moths and butterflies (Hexapoda: Lepidoptera)?

This paper addresses the question of whether one can economically improve the robustness of a molecular phylogeny estimate by increasing gene sampling in only a subset of taxa, without having the analysis invalidated by artifacts arising from large blocks of missing data. Our case study stems from an ongoing effort to resolve poorly understood deeper relationships in the large clade Ditrysia ( > 150,000 species) of the insect order Lepidoptera (butterflies and moths). Seeking to remedy the overall weak support for deeper divergences in an initial study based on five nuclear genes (6.6 kb) in 123 exemplars, we nearly tripled the total gene sample (to 26 genes, 18.4 kb) but only in a third (41) of the taxa. The resulting partially augmented data matrix (45% intentionally missing data) consistently increased bootstrap support for groupings previously identified in the five-gene (nearly) complete matrix, while introducing no contradictory groupings of the kind that missing data have been predicted to produce. Our results add to growing evidence that data sets differing substantially in gene and taxon sampling can often be safely and profitably combined. The strongest overall support for nodes above the family level came from including all nucleotide changes, while partitioning sites into sets undergoing mostly nonsynonymous versus mostly synonymous change. In contrast, support for the deepest node for which any persuasive molecular evidence has yet emerged (78-85% bootstrap) was weak or nonexistent unless synonymous change was entirely excluded, a result plausibly attributed to compositional heterogeneity. This node (Gelechioidea + Apoditrysia), tentatively proposed by previous authors on the basis of four morphological synapomorphies, is the first major subset of ditrysian superfamilies to receive strong statistical support in any phylogenetic study. A "more-genes-only" data set (41 taxa×26 genes) also gave strong signal for a second deep grouping (Macrolepidoptera) that was obscured, but not strongly contradicted, in more taxon-rich analyses.

Targeted amplicon sequencing (TAS): a scalable next-gen approach to multilocus, multitaxa phylogenetics

Next-gen sequencing technologies have revolutionized data collection in genetic studies and advanced genome biology to novel frontiers. However, to date, next-gen technologies have been used principally for whole genome sequencing and transcriptome sequencing. Yet many questions in population genetics and systematics rely on sequencing specific genes of known function or diversity levels. Here, we describe a targeted amplicon sequencing (TAS) approach capitalizing on next-gen capacity to sequence large numbers of targeted gene regions from a large number of samples. Our TAS approach is easily scalable, simple in execution, neither time-nor labor-intensive, relatively inexpensive, and can be applied to a broad diversity of organisms and/or genes. Our TAS approach includes a bioinformatic application, BarcodeCrucher, to take raw next-gen sequence reads and perform quality control checks and convert the data into FASTA format organized by gene and sample, ready for phylogenetic analyses. We demonstrate our approach by sequencing targeted genes of known phylogenetic utility to estimate a phylogeny for the Pancrustacea. We generated data from 44 taxa using 68 different 10-bp multiplexing identifiers. The overall quality of data produced was robust and was informative for phylogeny estimation. The potential for this method to produce copious amounts of data from a single 454 plate (e.g., 325 taxa for 24 loci) significantly reduces sequencing expenses incurred from traditional Sanger sequencing. We further discuss the advantages and disadvantages of this method, while offering suggestions to enhance the approach.