Loss of sexual reproduction is considered an evolutionary dead end for metazoans, but bdelloid rotifers challenge this view as they appear to have persisted asexually for millions of years 1 . Neither male sex organs nor meiosis have ever been observed in these microscopic animals: oocytes are formed through mitotic divisions, with no reduction of chromosome number and no indication of chromosome pairing 2 . However, current evidence does not exclude that they may engage in sex on rare, cryptic occasions. Here we report the genome of a bdelloid rotifer, Adineta vaga (Davis, 1873) 3 , and show that its structure is incompatible with conventional meiosis. At gene scale, the genome of A. vaga is tetraploid and comprises both anciently duplicated segments and less divergent allelic regions. However, in contrast to sexual species, the allelic regions are rearranged and sometimes even found on the same chromosome. Such structure does not allow meiotic pairing; instead, we find abundant evidence of gene conversion, which may limit the accumulation of deleterious mutations in the absence of meiosis. Gene families involved in resistance to oxidation, carbohydrate metabolism and defence against transposons are significantly expanded, which may explain why transposable elements cover only 3% of the assembled sequence. Furthermore, 8% of the genes are likely to be of non-metazoan origin and were probably acquired horizontally. This apparent convergence between bdelloids and prokaryotes sheds new light on the evolutionary significance of sex.With more than 460 described species 4 , bdelloid rotifers ( Fig. 1) represent the highest metazoan taxonomic rank in which males, hermaphrodites and meiosis are unknown. Such persistence and diversification of an ameiotic clade of animals are in contradiction with the supposed long-term disadvantages of asexuality, making bdelloids an 'evolutionary scandal' 5 . Another unusual feature of bdelloid rotifers is their extreme resistance to desiccation at any stage of their life cycle 6 , enabling these microscopic animals to dwell in ephemeral freshwater habitats such as mosses, lichens and forest litter; this ability is presumably the source of their extreme resistance to ionizing radiation 7 .We assembled the genome of a clonal A. vaga lineage into separate haplotypes with a N 50 of 260 kilobases (kb) (that is, half of the assembly was composed of fragments longer than 260 kb). Assembly size was 218 megabases (Mb) but 26 Mb of the sequence had twice the average sequencing coverage, suggesting that some nearly identical regions were not resolved during assembly ( Supplementary Fig. 3); hence, the total genome size is likely to be 244 Mb, which corresponds to the estimate obtained independently using fluorometry (Supplementary Note C2). Annotation of the complete assembly (including all haplotypes) yielded 49,300 genes. Intragenomic sequence comparisons revealed numerous homologous blocks with conserved gene order (colinear regions). For each such block we computed the per-site synonymous d...

BackgroundThe duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.ResultsWe found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.ConclusionsOur results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12915-017-0399-x) contains supplementary material, which is available to authorized users.

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

The program phase is widely used for Bayesian inference of haplotypes from diploid genotypes; however, manually creating phase input files from sequence alignments is an error-prone and time-consuming process, especially when dealing with numerous variable sites and/or individuals. Here, a web tool called seqphase is presented that generates phase input files from fasta sequence alignments and converts phase output files back into fasta. During the production of the phase input file, several consistency checks are performed on the dataset and suitable command line options to be used for the actual phase data analysis are suggested. seqphase was written in perl and is freely accessible over the Internet at the address http://www.mnhn.fr/jfflot/seqphase.

Genomic analyses of microbial populations in their natural environment remain limited by the difficulty to assemble full genomes of individual species. Consequently, the chromosome organization of microorganisms has been investigated in a few model species, but the extent to which the features described can be generalized to other taxa remains unknown. Using controlled mixes of bacterial and yeast species, we developed meta3C, a metagenomic chromosome conformation capture approach that allows characterizing individual genomes and their average organization within a mix of organisms. Not only can meta3C be applied to species already sequenced, but a single meta3C library can be used for assembling, scaffolding and characterizing the tridimensional organization of unknown genomes. By applying meta3C to a semi-complex environmental sample, we confirmed its promising potential. Overall, this first meta3C study highlights the remarkable diversity of microorganisms chromosome organization, while providing an elegant and integrated approach to metagenomic analysis.DOI: http://dx.doi.org/10.7554/eLife.03318.001

Describing biological diversity is a challenging endeavour, especially for the small, cryptic animals that make up the meiofauna. The field of DNA taxonomy, i.e., the use of DNA to delineate species boundaries, is rapidly growing and changing; herein we review the recent advances in the acquisition of DNA sequence data and the analytical tools for DNA-based species delimitation, with a focus on applications to the meiofauna. After providing general guidelines on the data collection and analysis steps (sampling design, sequencing, phasing of nuclear markers, and sequence alignment), we explain the rationale and usage of several widely used or promising methods developed for delineating species from single-locus data sets (distance-based DNA barcoding, automated barcode gap discovery, K/θ, generalized mixed Yule-coalescent models, Poisson tree process model, and haplowebs). As it is increasingly recognised that several loci are required to delineate species accurately, we then briefly outline multilocus species delimitation approaches (Structure, Structurama, Bayesian phylogenetics & phylogeography, SpedeSTEM, O'Meara's heuristic search, and several newly published Bayesian approaches based on the multispecies coalescent).

Closing gaps in draft genome assemblies can be costly and time-consuming, and published genomes are therefore often left ‘unfinished.’ Here we show that genome-wide chromosome conformation capture (3C) data can be used to overcome these limitations, and present a computational approach rooted in polymer physics that determines the most likely genome structure using chromosomal contact data. This algorithm—named GRAAL—generates high-quality assemblies of genomes in which repeated and duplicated regions are accurately represented and offers a direct probabilistic interpretation of the computed structures. We first validated GRAAL on the reference genome of Saccharomyces cerevisiae, as well as other yeast isolates, where GRAAL recovered both known and unknown complex chromosomal structural variations. We then applied GRAAL to the finishing of the assembly of Trichoderma reesei and obtained a number of contigs congruent with the know karyotype of this species. Finally, we showed that GRAAL can accurately reconstruct human chromosomes from either fragments generated in silico or contigs obtained from de novo assembly. In all these applications, GRAAL compared favourably to recently published programmes implementing related approaches.

Six variable sequence markers are developed and analyzed to find out species boundaries in Hawaiian corals of the genus Pocillopora: the putative mitochondrial control region; a recently discovered, hypervariable mitochondrial open reading frame; the internal transcribed spacer 2 (ITS2), located in the nuclear ribosomal DNA; three nuclear introns of calmodulin, elongation factor-1alpha and the ATP synthase beta subunit. Using the first two markers, we identify five distinct mitochondrial lineages and these lineages are compatible with morphology. The situation is more complex with nuclear markers since more than two haplotypes are observed in some individuals. To detect clusters of individuals, haplotype networks are constructed with additional connections drawn between co-occurring haplotypes to delineate potential fields for recombination: few clusters of nuclear haplotypes are found to correspond to clusters of individuals, but those that are detected (mostly in the ITS2 dataset) are also compatible with morphology.