No signal of deleterious mutation accumulation in conserved gene sequences of extant asexual hexapods.

Loss of sex and recombination is generally assumed to impede the effectiveness of purifying selection and to result in the accumulation of slightly deleterious mutations. Empirical evidence for this has come from several studies investigating mutational load in a small number of individual genes. However, recent whole transcriptome based studies have yielded inconsistent results, hence questioning the validity of the assumption of mutational meltdown in asexual populations. Here, we study the effectiveness of purifying selection in eight asexual hexapod lineages and their sexual relatives, as present in the 1 K Insect Transcriptome Evolution (1KITE) project, covering eight hexapod groups. We analyse the accumulation of slightly deleterious nonsynonymous and synonymous point mutations in 99 single copy orthologue protein-coding loci shared among the investigated taxa. While accumulation rates of nonsynonymous mutations differed between genes and hexapod groups, we found no effect of reproductive mode on the effectiveness of purifying selection acting at nonsynonymous and synonymous sites. Although the setup of this study does not fully rule out nondetection of subtle effects, our data does not support the established consensus of asexual lineages undergoing 'mutational meltdown'

No signal of deleterious mutation accumulation in conserved gene sequences of extant asexual hexapods.

Loss of sex and recombination is generally assumed to impede the effectiveness of purifying selection and to result in the accumulation of slightly deleterious mutations. Empirical evidence for this has come from several studies investigating mutational load in a small number of individual genes. However, recent whole transcriptome based studies have yielded inconsistent results, hence questioning the validity of the assumption of mutational meltdown in asexual populations. Here, we study the effectiveness of purifying selection in eight asexual hexapod lineages and their sexual relatives, as present in the 1 K Insect Transcriptome Evolution (1KITE) project, covering eight hexapod groups. We analyse the accumulation of slightly deleterious nonsynonymous and synonymous point mutations in 99 single copy orthologue protein-coding loci shared among the investigated taxa. While accumulation rates of nonsynonymous mutations differed between genes and hexapod groups, we found no effect of reproductive mode on the effectiveness of purifying selection acting at nonsynonymous and synonymous sites. Although the setup of this study does not fully rule out nondetection of subtle effects, our data does not support the established consensus of asexual lineages undergoing 'mutational meltdown'

A minimum reporting standard for multiple sequence alignments

Multiple sequence alignments (MSAs) play a pivotal role in studies of molecular sequence data, but nobody has developed a minimum reporting standard (MRS) to quantify the completeness of MSAs. We present an MRS that relies on four simple completeness metrics. The metrics are implemented in AliStat, a program developed to support the MRS. A survey of published MSAs illustrates the benefits and unprecedented transparency offered by the MRS.CSIRO

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution

Background: The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. Results: Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. Conclusions: The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus)

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.Danielle N. Stringer, Terry Bertozzi, Karen Meusemann, Steven Delean, Michelle T. Guzik, Simon M. Tierney ... et al

Potential contribution of ancient introgression to the evolution of a derived reproductive strategy in ricefishes

Transitions from no parental care to extensive care are costly and involve major changes in life history, behaviour and morphology. Nevertheless, in Sulawesi ricefishes, pelvic brooding evolved from transfer brooding in two distantly related lineages within the genera Adrianichthys and Oryzias, respectively. Females of pelvic brooding species carry their eggs attached to their belly until the fry hatches. Despite their phylogenetic distance, both pelvic brooding lineages share a set of external morphological traits. A recent study found no direct gene flow between pelvic brooding lineages, suggesting independent evolution of the derived reproductive strategy. Convergent evolution can, however, also rely on repeated sorting of pre-existing variation of an admixed ancestral population, especially when subjected to similar external selection pressures. We thus used a multi-species coalescent (MSC) model and D-statistics to identify gene-tree - species-tree incongruencies, to evaluate the evolution of pelvic brooding with respect to inter-specific gene flow not only between pelvic brooding lineages, but between pelvic brooding lineages and other Sulawesi ricefish lineages. We found a general network-like evolution in Sulawesi ricefishes and as previously reported, no gene flow between the pelvic brooding lineages. Instead, we found hybridization between the ancestor of pelvic brooding Oryzias and the common ancestor of the Oryzias species from the Lake Poso area. We further detected signs of introgression within the confidence interval of a quantitative trait locus (QTL) associated with pelvic brooding in O. eversi. Our results hint towards a contribution of ancient standing genetic variation to the evolution of pelvic brooding in Oryzias

Improved Phylogenetic Analyses Corroborate a Plausible Position of Martialis heureka in the Ant Tree of Life

Martialinae are pale, eyeless and probably hypogaeic predatory ants. Morphological character sets suggest a close relationship to the ant subfamily Leptanillinae. Recent analyses based on molecular sequence data suggest that Martialinae are the sister group to all extant ants. However, by comparing molecular studies and different reconstruction methods, the position of Martialinae remains ambiguous. While this sister group relationship was well supported by Bayesian partitioned analyses, Maximum Likelihood approaches could not unequivocally resolve the position of Martialinae. By re-analysing a previous published molecular data set, we show that the Maximum Likelihood approach is highly appropriate to resolve deep ant relationships, especially between Leptanillinae, Martialinae and the remaining ant subfamilies. Based on improved alignments, alignment masking, and tree reconstructions with a sufficient number of bootstrap replicates, our results strongly reject a placement of Martialinae at the first split within the ant tree of life. Instead, we suggest that Leptanillinae are a sister group to all other extant ant subfamilies, whereas Martialinae branch off as a second lineage. This assumption is backed by approximately unbiased (AU) tests, additional Bayesian analyses and split networks. Our results demonstrate clear effects of improved alignment approaches, alignment masking and data partitioning. We hope that our study illustrates the importance of thorough, comprehensible phylogenetic analyses using the example of ant relationships

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths.

Butterflies and moths (Lepidoptera) are one of the major superradiations of insects, comprising nearly 160,000 described extant species. As herbivores, pollinators, and prey, Lepidoptera play a fundamental role in almost every terrestrial ecosystem. Lepidoptera are also indicators of environmental change and serve as models for research on mimicry and genetics. They have been central to the development of coevolutionary hypotheses, such as butterflies with flowering plants and moths' evolutionary arms race with echolocating bats. However, these hypotheses have not been rigorously tested, because a robust lepidopteran phylogeny and timing of evolutionary novelties are lacking. To address these issues, we inferred a comprehensive phylogeny of Lepidoptera, using the largest dataset assembled for the order (2,098 orthologous protein-coding genes from transcriptomes of 186 species, representing nearly all superfamilies), and dated it with carefully evaluated synapomorphy-based fossils. The oldest members of the Lepidoptera crown group appeared in the Late Carboniferous (∼300 Ma) and fed on nonvascular land plants. Lepidoptera evolved the tube-like proboscis in the Middle Triassic (∼241 Ma), which allowed them to acquire nectar from flowering plants. This morphological innovation, along with other traits, likely promoted the extraordinary diversification of superfamily-level lepidopteran crown groups. The ancestor of butterflies was likely nocturnal, and our results indicate that butterflies became day-flying in the Late Cretaceous (∼98 Ma). Moth hearing organs arose multiple times before the evolutionary arms race between moths and bats, perhaps initially detecting a wide range of sound frequencies before being co-opted to specifically detect bat sonar. Our study provides an essential framework for future comparative studies on butterfly and moth evolution

Phylogenomics of the superfamily Dytiscoidea (Coleoptera: Adephaga) with an evaluation of phylogenetic conflict and systematic error.

The beetle superfamily Dytiscoidea, placed within the suborder Adephaga, comprises six families. The phylogenetic relationships of these families, whose species are aquatic, remain highly contentious. In particular the monophyly of the geographically disjunct Aspidytidae (China and South Africa) remains unclear. Here we use a phylogenomic approach to demonstrate that Aspidytidae are indeed monophyletic, as we inferred this phylogenetic relationship from analyzing nucleotide sequence data filtered for compositional heterogeneity and from analyzing amino-acid sequence data. Our analyses suggest that Aspidytidae are the sister group of Amphizoidae, although the support for this relationship is not unequivocal. A sister group relationship of Hygrobiidae to a clade comprising Amphizoidae, Aspidytidae, and Dytiscidae is supported by analyses in which model assumptions are violated the least. In general, we find that both concatenation and the applied coalescent method are sensitive to the effect of among-species compositional heterogeneity. Four-cluster likelihood-mapping suggests that despite the substantial size of the dataset and the use of advanced analytical methods, statistical support is weak for the inferred phylogenetic placement of Hygrobiidae. These results indicate that other kinds of data (e.g. genomic meta-characters) are possibly required to resolve the above-specified persisting phylogenetic uncertainties. Our study illustrates various data-driven confounding effects in phylogenetic reconstructions and highlights the need for careful monitoring of model violations prior to phylogenomic analysis

iPhy: an integrated phylogenetic workbench for supermatrix analyses

<p>Abstract</p> <p>Background</p> <p>The increasing availability of molecular sequence data means that the accuracy of future phylogenetic studies is likely to by limited by systematic bias and taxon choice rather than by data. In order to take advantage of increasing datasets, user-friendly tools are required to facilitate phylogenetic analyses and to reduce duplication of dataset assembly efforts. Current phylogenetic pipelines are dependency-heavy and have significant technical barriers to use.</p> <p>Results</p> <p>Here we present iPhy, a web application that lets non-technical users assemble, share and analyse DNA sequence datasets for multigene phylogenetic investigations. Built on a simple client-server architecture, iPhy eases the collection of gene sets for analysis, facilitates alignment and reliably generates phylogenetic analysis-ready data files. Phylogenetic trees generated in external programs can be imported and stored, and iPhy integrates with iTol to allow trees to be displayed with rich data annotation. The datasets collated in iPhy can be shared through the client interface. We show how systematic biases can be addressed by using explicit criteria when selecting sequences for analysis from a large dataset. A representative instance of iPhy can be accessed at iphy.bio.ed.ac.uk, but the toolkit can also be deployed on a local server for advanced users.</p> <p>Conclusions</p> <p>iPhy provides an easy-to-use environment for the assembly, analysis and sharing of large phylogenetic datasets, while encouraging best practices in terms of phylogenetic analysis and taxon selection.</p