Detection of ribonucleotides embedded in DNA by Nanopore sequencing.

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Amplicon-Based Microbiome Profiling: From Second- to Third-Generation Sequencing for Higher Taxonomic Resolution.

The 16S rRNA amplicon-based sequencing approach represents the most common and cost-effective strategy with great potential for microbiome profiling. The use of second-generation sequencing (NGS) technologies has led to protocols based on the amplification of one or a few hypervariable regions, impacting the outcome of the analysis. Nowadays, comparative studies are necessary to assess different amplicon-based approaches, including the full-locus sequencing currently feasible thanks to third-generation sequencing (TGS) technologies. This study compared three different methods to achieve the deepest microbiome taxonomic characterization: (a) the single-region approach, (b) the multiplex approach, covering several regions of the target gene/region, both based on NGS short reads, and (c) the full-length approach, which analyzes the whole length of the target gene thanks to TGS long reads. Analyses carried out on benchmark microbiome samples, with a known taxonomic composition, highlighted a different classification performance, strongly associated with the type of hypervariable regions and the coverage of the target gene. Indeed, the full-length approach showed the greatest discriminating power, up to species level, also on complex real samples. This study supports the transition from NGS to TGS for the study of the microbiome, even if experimental and bioinformatic improvements are still necessary.

HaploCoV: unsupervised classification and rapid detection of novel emerging variants of SARS-CoV-2.

Accurate and timely monitoring of the evolution of SARS-CoV-2 is crucial for identifying and tracking potentially more transmissible/virulent viral variants, and implement mitigation strategies to limit their spread. Here we introduce HaploCoV, a novel software framework that enables the exploration of SARS-CoV-2 genomic diversity through space and time, to identify novel emerging viral variants and prioritize variants of potential epidemiological interest in a rapid and unsupervised manner. HaploCoV can integrate with any classification/nomenclature and incorporates an effective scoring system for the prioritization of SARS-CoV-2 variants. By performing retrospective analyses of more than 11.5 M genome sequences we show that HaploCoV demonstrates high levels of accuracy and reproducibility and identifies the large majority of epidemiologically relevant viral variants - as flagged by international health authorities - automatically and with rapid turn-around times.Our results highlight the importance of the application of strategies based on the systematic analysis and integration of regional data for rapid identification of novel, emerging variants of SARS-CoV-2. We believe that the approach outlined in this study will contribute to relevant advances to current and future genomic surveillance methods.

Yamanaka Factors in the Budding Tunicate Botryllus schlosseri Show a Shared Spatio-Temporal Expression Pattern in Chordates.

In vertebrates, the four transcription factors Sox2, c-Myc, Pou5f1 and Klf4 are involved in the differentiation of several tissues during vertebrate embryogenesis; moreover, they are normally co-expressed in embryonic stem cells and play roles in pluripotency, self-renewal, and maintenance of the undifferentiated state in adult cells. The in vitro forced co-expression of these factors, named Yamanaka factors (YFs), induces pluripotency in human or mouse fibroblasts. Botryllus schlosseri is a colonial tunicate undergoing continuous stem cell-mediated asexual development, providing a valuable model system for the study of pluripotency in the closest living relatives of vertebrates. In this study, we identified B. schlosseri orthologs of human Sox2 and c-Myc genes, as well as the closest homologs of the vertebrate-specific Pou5f1 gene, through an in-depth evolutionary analysis of the YF gene families in tunicates and other deuterostomes. Then, we studied the expression of these genes during the asexual cycle of B. schlosseri using in situ hybridization in order to investigate their possible involvement in tissue differentiation and in pluripotency maintenance. Our results show a shared spatio-temporal expression pattern consistent with the reported functions of these genes in invertebrate and vertebrate embryogenesis. Moreover, Myc, SoxB1 and Pou3 were expressed in candidate stem cells residing in their niches, while Pou2 was found expressed exclusively in the immature previtellogenic oocytes, both in gonads and circulating in the colonial vascular system. Our data suggest that Myc, SoxB1 and Pou3 may be individually involved in the differentiation of the same territories seen in other chordates, and that, together, they may play a role in stemness even in this colonial ascidian.

Finding Aquaporins in Annelids: An Evolutionary Analysis and a Case Study.

Aquaporins (AQPs) are a family of membrane channels facilitating diffusion of water and small solutes into and out of cells. Despite their biological relevance in osmoregulation and ubiquitous distribution throughout metazoans, the presence of AQPs in annelids has been poorly investigated. Here, we searched and annotated Aqp sequences in public genomes and transcriptomes of annelids, inferred their evolutionary relationships through phylogenetic analyses and discussed their putative physiological relevance. We identified a total of 401 Aqp sequences in 27 annelid species, including 367 sequences previously unrecognized as Aqps. Similar to vertebrates, phylogenetic tree reconstructions clustered these annelid Aqps in four clades: AQP1-like, AQP3-like, AQP8-like and AQP11-like. We found no clear indication of the existence of paralogs exclusive to annelids; however, several gene duplications seem to have occurred in the ancestors of some Sedentaria annelid families, mainly in the AQP1-like clade. Three of the six Aqps annotated in Alitta succinea, an estuarine annelid showing high salinity tolerance, were validated by RT-PCR sequencing, and their similarity to human AQPs was investigated at the level of "key" conserved residues and predicted three-dimensional structure. Our results suggest a diversification of the structures and functions of AQPs in Annelida comparable to that observed in other taxa.

Evaluating the Efficiency of DNA Metabarcoding to Analyze the Diet of Hippocampus guttulatus (Teleostea: Syngnathidae).

Seahorses are considered a flagship species for conservation efforts and due to their conservation status, improving knowledge on their dietary composition while applying a non-invasive approach, could be useful. Using Hippocampus guttulatus as a case study, the present study represents pioneering research into investigating the diet of seahorses by NGS-based DNA metabarcoding of fecal samples. The study developed and tested the protocol for fecal DNA metabarcoding during the feeding trials where captive seahorses were fed on a diet of known composition; the process was subsequently applied on fecal samples collected from wild individuals. The analysis of samples collected during the feeding trials indicated the reliability of the applied molecular approach by allowing the characterization of the effectively ingested prey. In the field study, among detected prey species, results revealed that the majority of the seahorse samples contained taxa such as Amphipoda, Decapoda, Isopoda, and Calanoida, while less common prey taxa were Gastropoda and Polyplacophora. As only a small amount of starting fecal material is needed and the sampling procedure is neither invasive nor lethal. The present study indicates DNA metabarcoding as useful for investigating seahorse diet and could help define management and conservation actions.

An elongated COI fragment to discriminate botryllid species and as an improved ascidian DNA barcode.

Botryllids are colonial ascidians widely studied for their potential invasiveness and as model organisms, however the morphological description and discrimination of these species is very problematic, leading to frequent specimen misidentifications. To facilitate species discrimination and detection of cryptic/new species, we developed new barcoding primers for the amplification of a COI fragment of about 860 bp (860-COI), which is an extension of the common Folmer's barcode region. Our 860-COI was successfully amplified in 177 worldwide-sampled botryllid colonies. Combined with morphological analyses, 860-COI allowed not only discriminating known species, but also identifying undescribed and cryptic species, resurrecting old species currently in synonymy, and proposing the assignment of clade D of the model organism Botryllus schlosseri to Botryllus renierii. Importantly, within clade A of B. schlosseri, 860-COI recognized at least two candidate species against only one recognized by the Folmer's fragment, underlining the need of further genetic investigations on this clade. This result also suggests that the 860-COI could have a greater ability to diagnose cryptic/new species than the Folmer's fragment at very short evolutionary distances, such as those observed within clade A. Finally, our new primers simplify the amplification of 860-COI even in non-botryllid ascidians, suggesting their wider usefulness in ascidians.

Comparative Genomics Reveals Early Emergence and Biased Spatiotemporal Distribution of SARS-CoV-2.

Effective systems for the analysis of molecular data are fundamental for monitoring the spread of infectious diseases and studying pathogen evolution. The rapid identification of emerging viral strains, and/or genetic variants potentially associated with novel phenotypic features is one of the most important objectives of genomic surveillance of human pathogens and represents one of the first lines of defense for the control of their spread. During the COVID 19 pandemic, several taxonomic frameworks have been proposed for the classification of SARS-Cov-2 isolates. These systems, which are typically based on phylogenetic approaches, represent essential tools for epidemiological studies as well as contributing to the study of the origin of the outbreak. Here, we propose an alternative, reproducible, and transparent phenetic method to study changes in SARS-CoV-2 genomic diversity over time. We suggest that our approach can complement other systems and facilitate the identification of biologically relevant variants in the viral genome. To demonstrate the validity of our approach, we present comparative genomic analyses of more than 175,000 genomes. Our method delineates 22 distinct SARS-CoV-2 haplogroups, which, based on the distribution of high-frequency genetic variants, fall into four major macrohaplogroups. We highlight biased spatiotemporal distributions of SARS-CoV-2 genetic profiles and show that seven of the 22 haplogroups (and of all of the four haplogroup clusters) showed a broad geographic distribution within China by the time the outbreak was widely recognized-suggesting early emergence and widespread cryptic circulation of the virus well before its isolation in January 2020. General patterns of genomic variability are remarkably similar within all major SARS-CoV-2 haplogroups, with UTRs consistently exhibiting the greatest variability, with s2m, a conserved secondary structure element of unknown function in the 3'-UTR of the viral genome showing evidence of a functional shift. Although several polymorphic sites that are specific to one or more haplogroups were predicted to be under positive or negative selection, overall our analyses suggest that the emergence of novel types is unlikely to be driven by convergent evolution and independent fixation of advantageous substitutions, or by selection of recombined strains. In the absence of extensive clinical metadata for most available genome sequences, and in the context of extensive geographic and temporal biases in the sampling, many questions regarding the evolution and clinical characteristics of SARS-CoV-2 isolates remain open. However, our data indicate that the approach outlined here can be usefully employed in the identification of candidate SARS-CoV-2 genetic variants of clinical and epidemiological importance.

Next generation sequencing of SARS-CoV-2 genomes: challenges, applications and opportunities.

Various next generation sequencing (NGS) based strategies have been successfully used in the recent past for tracing origins and understanding the evolution of infectious agents, investigating the spread and transmission chains of outbreaks, as well as facilitating the development of effective and rapid molecular diagnostic tests and contributing to the hunt for treatments and vaccines. The ongoing COVID-19 pandemic poses one of the greatest global threats in modern history and has already caused severe social and economic costs. The development of efficient and rapid sequencing methods to reconstruct the genomic sequence of SARS-CoV-2, the etiological agent of COVID-19, has been fundamental for the design of diagnostic molecular tests and to devise effective measures and strategies to mitigate the diffusion of the pandemic. Diverse approaches and sequencing methods can, as testified by the number of available sequences, be applied to SARS-CoV-2 genomes. However, each technology and sequencing approach has its own advantages and limitations. In the current review, we will provide a brief, but hopefully comprehensive, account of currently available platforms and methodological approaches for the sequencing of SARS-CoV-2 genomes. We also present an outline of current repositories and databases that provide access to SARS-CoV-2 genomic data and associated metadata. Finally, we offer general advice and guidelines for the appropriate sharing and deposition of SARS-CoV-2 data and metadata, and suggest that more efficient and standardized integration of current and future SARS-CoV-2-related data would greatly facilitate the struggle against this new pathogen. We hope that our 'vademecum' for the production and handling of SARS-CoV-2-related sequencing data, will contribute to this objective.

LC-HRMS-Database Screening Metrics for Rapid Prioritization of Samples to Accelerate the Discovery of Structurally New Natural Products.

In order to accelerate the isolation and characterization of structurally new or novel secondary metabolites, it is crucial to develop efficient strategies that prioritize samples with greatest promise early in the workflow so that resources can be utilized in a more efficient and cost-effective manner. We have developed a metrics-based prioritization approach using exact LC-HRMS, which uses data for 24 618 marine natural products held in the PharmaSea database. Each sample was evaluated and allocated a metric score by a software algorithm based on the ratio of new masses over the total (sample novelty), ratio of known masses over the total (chemical novelty), number of peaks above a defined peak area threshold (sample complexity), and peak area (sample diversity). Samples were then ranked and prioritized based on these metric scores. To validate the approach, eight marine sponges and six tunicate samples collected from the Fiji Islands were analyzed, metric scores calculated, and samples targeted for isolation and characterization of new compounds. Structures of new compounds were elucidated by spectroscopic techniques, including 1D and 2D NMR, MS, and MS/MS. Structures were confirmed by computer-assisted structure elucidation methods (CASE) using the ACD/Structure Elucidator Suite.

Extensive mitochondrial gene rearrangements in Ctenophora: insights from benthic Platyctenida.

BACKGROUND: Complete mitochondrial (mt) genomes have been sequenced for thousands of animals and represent a molecule of choice for many evolutionary studies. Nevertheless, some animal groups have remained under-sampled. Ctenophora (comb jellies) is one such example, with only two complete mt sequences determined hitherto for this phylum, which encompasses ca. 150-200 described species. This lack of data derives from the extremely fast mt evolutionary rate in this lineage, complicating primer design and DNA amplification. Indeed, in the two ctenophore mt genomes sequenced to date, i.e. those of Mnemiopsis leidyi (order Lobata) and Pleurobrachia bachei (order Cydippida), both rRNA and protein coding genes exhibit an extraordinary size reduction and have highly derived sequences. Additionally, all tRNAs, and the atp6 and atp8 genes are absent. In order to determine whether these characteristics are shared by other ctenophores, we obtained the complete mt genomes of three benthic ctenophores belonging to the so far unsampled order of Platyctenida: Coeloplana loyai, Coeloplana yulianicorum and Vallicula multiformis. RESULTS: The mt genomes of benthic ctenophores reveal the same peculiarities found in Mnemiopsis and Pleurobrachia, demonstrating that the fast evolutionary rate is a general trait of the ctenophore mt genomes. Our results also indicate that this high evolutionary rate not only affects the nucleotide substitution but also gene rearrangements. Indeed, gene order was highly rearranged among representatives of the different taxonomic orders in which it was close to random, but also quite variable within Platyctenida, in which the genera Coeloplana and Vallicula share only four conserved synteny blocks. However, the two congeneric Coeloplana species display exactly the same gene order. Because of the extreme evolutionary rate, our phylogenetic analyses were unable to resolve the phylogenetic position of ctenophores within metazoans or the relationships among the different Ctenophora orders. Comparative sequence-analyses allowed us to correct the annotation of the Pleurobrachia mt genome, confirming the absence of tRNAs, the presence of both rRNA genes, and the existence of a reassignment of codon TGA from tryptophan to serine for this species. CONCLUSIONS: Since Platyctenida is an early diverging lineage among Ctenophora, our findings suggest that the mt traits described above are ancestral characteristics of this phylum.

Fixation, description and DNA barcode of a neotype for Botryllus schlosseri (Pallas, 1766) (Tunicata, Ascidiacea).

Botryllus schlosseri is a widespread colonial ascidian commonly considered cosmopolitan and amply used as model for researches ranging from developmental biology to immunobiology. Recently, molecular data lead to hypothesize that the species named B. schlosseri may consist of more than a single taxon. Indeed, five highly divergent clades, named A-E, have been genetically identified and are referred as cryptic species. In this context, and lacking both a type and a detailed morphological description, we believe that it is necessary, as a taxonomic reference point, to designate a neotype and re-describe the species. Therefore, a sample from the Lagoon of Venice (Adriatic Sea, Italy) was deposited as neotype in the Natural History Museum of Venice (Italy), preserved both in formalin and in 90% ethanol. Here we provide a morphological description of the suggested neotype of B. schlosseri that takes into account several developmental stages (oozooid, zooid of first blastogenic generations, and mature zooid) and is carefully compared with the previous descriptions of samples coming from other European and non-European localities. Finally, we associate our morphological description to a "DNA barcode", consisting in a long fragment of the mitochondrial COI gene. Our description is associated to clade A, although at now we cannot guarantee that this association is univocal.

Back to solitude: Solving the phylogenetic position of the Diazonidae using molecular and developmental characters.

The order Aplousobranchia (Chordata, Ascidiacea) contains approximately 1500 species distributed worldwide. Their phylogeny, however, remains unclear, with unresolved family relationships. While most Aplousobranchia are colonial, debates exist concerning the phylogenetic position of families such as the Diazonidae and Cionidae, which exhibit a solitary lifestyle and share morphological characteristics with both Aplousobranchia and Phlebobranchia orders. To clarify the phylogenetic position of the Diazonidae and Cionidae, we determined the complete mitochondrial sequence of the solitary diazonid Rhopalaea idoneta. The phylogenetic reconstruction based on the 13 mitochondrial protein coding genes strongly supports a positioning of Diazonidae well-nested within the Aplousobranchia rather than a positioning as a sister clade of the Aplousobranchia. In addition, we examined the regenerative ability of R. idoneta. Similar to colonial Aplousobranchia, R. idoneta was found to be able to completely regenerate its thorax. Ciona, also known to possess high regenerative abilities, is the Aplousobranchia sister clade rather than a member of the Phlebobranchia. Our results thus indicate that the colonial lifestyle was acquired in the Aplousobranchia, starting from a Ciona-like solitary ancestor and secondarily lost in Diazonidae representatives such as Rhopalaea. The solitary lifestyle of Rhopalaea is thus a derived characteristic rather than an ancestral trait.

ANISEED 2015: a digital framework for the comparative developmental biology of ascidians.

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr.

Morphological Differences between Larvae of the Ciona intestinalis Species Complex: Hints for a Valid Taxonomic Definition of Distinct Species.

The cosmopolitan ascidian Ciona intestinalis is the most common model species of Tunicata, the sister-group of Vertebrata, and widely used in developmental biology, genomics and evolutionary studies. Recently, molecular studies suggested the presence of cryptic species hidden within the C. intestinalis species, namely C. intestinalis type A and type B. So far, no substantial morphological differences have been identified between individuals belonging to the two types. Here we present morphometric, immunohistochemical, and histological analyses, as well as 3-D reconstructions, of late larvae obtained by cross-fertilization experiments of molecularly determined type A and type B adults, sampled in different seasons and in four different localities. Our data point to quantitative and qualitative differences in the trunk shape of larvae belonging to the two types. In particular, type B larvae exhibit a longer pre-oral lobe, longer and relatively narrower total body length, and a shorter ocellus-tail distance than type A larvae. All these differences were found to be statistically significant in a Discriminant Analysis. Depending on the number of analyzed parameters, the obtained discriminant function was able to correctly classify > 93% of the larvae, with the remaining misclassified larvae attributable to the existence of intra-type seasonal variability. No larval differences were observed at the level of histology and immunohistochemical localization of peripheral sensory neurons. We conclude that type A and type B are two distinct species that can be distinguished on the basis of larval morphology and molecular data. Since the identified larval differences appear to be valid diagnostic characters, we suggest to raise both types to the rank of species and to assign them distinct names.

Tissue-specific mtDNA abundance from exome data and its correlation with mitochondrial transcription, mass and respiratory activity.

Eukaryotic cells contain a population of mitochondria, variable in number and shape, which in turn contain multiple copies of a tiny compact genome (mtDNA) whose expression and function is strictly coordinated with the nuclear one. mtDNA copy number varies between different cell or tissues types, both in response to overall metabolic and bioenergetics demands and as a consequence or cause of specific pathological conditions. Here we present a novel and reliable methodology to assess the effective mtDNA copy number per diploid genome by investigating off-target reads obtained by whole-exome sequencing (WES) experiments. We also investigate whether and how mtDNA copy number correlates with mitochondrial mass, respiratory activity and expression levels. Analyzing six different tissues from three age- and sex-matched human individuals, we found a highly significant linear correlation between mtDNA copy number estimated by qPCR and the frequency of mtDNA off target WES reads. Furthermore, mtDNA copy number showed highly significant correlation with mitochondrial gene expression levels as measured by RNA-Seq as well as with mitochondrial mass and respiratory activity. Our methodology makes thus feasible, at a large scale, the investigation of mtDNA copy number in diverse cell-types, tissues and pathological conditions or in response to specific treatments.

Ascidian mitogenomics: comparison of evolutionary rates in closely related taxa provides evidence of ongoing speciation events.

Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (dN) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri, where most pairs show a dN value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1. Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Bioinformatics in Italy: BITS 2012, the ninth annual meeting of the Italian Society of Bioinformatics.

The BITS2012 meeting, held in Catania on May 2-4, 2012, brought together almost 100 Italian researchers working in the field of Bioinformatics, as well as students in the same or related disciplines. About 90 original research works were presented either as oral communication or as posters, representing a landscape of Italian current research in bioinformatics. This preface provides a brief overview of the meeting and introduces the manuscripts that were accepted for publication in this supplement, after a strict and careful peer-review by an International board of referees.

The genome sequence of the colonial chordate, Botryllus schlosseri.

Botryllus schlosseri is a colonial urochordate that follows the chordate plan of development following sexual reproduction, but invokes a stem cell-mediated budding program during subsequent rounds of asexual reproduction. As urochordates are considered to be the closest living invertebrate relatives of vertebrates, they are ideal subjects for whole genome sequence analyses. Using a novel method for high-throughput sequencing of eukaryotic genomes, we sequenced and assembled 580 Mbp of the B. schlosseri genome. The genome assembly is comprised of nearly 14,000 intron-containing predicted genes, and 13,500 intron-less predicted genes, 40% of which could be confidently parceled into 13 (of 16 haploid) chromosomes. A comparison of homologous genes between B. schlosseri and other diverse taxonomic groups revealed genomic events underlying the evolution of vertebrates and lymphoid-mediated immunity. The B. schlosseri genome is a community resource for studying alternative modes of reproduction, natural transplantation reactions, and stem cell-mediated regeneration. DOI:http://dx.doi.org/10.7554/eLife.00569.001.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes.

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.