Search

Background: Long-read sequencing technologies were launched a few years ago, and in contrast with short-read sequencing technologies, they offered a promise of solving assembly problems for large and complex genomes. Moreover by providing long-range information, it could also solve haplotype phasing. However, existing long-read technologies still have several limitations that complicate their use for most research laboratories, as well as in large and/or complex genome projects. In 2014, Oxford Nanopore released the MinION® device, a small and low-cost single-molecule nanopore sequencer, which offers the possibility of sequencing long DNA fragments. Results: The assembly of long reads generated using the Oxford Nanopore MinION® instrument is challenging as existing assemblers were not implemented to deal with long reads exhibiting close to 30% of errors. Here, we presented a hybrid approach developed to take advantage of data generated using MinION® device. We sequenced a well-known bacterium, Acinetobacter baylyi ADP1 and applied our method to obtain a highly contiguous (one single contig) and accurate genome assembly even in repetitive regions, in contrast to an Illumina-only assembly. Our hybrid strategy was able to generate NaS (Nanopore Synthetic-long) reads up to 60 kb that aligned entirely and with no error to the reference genome and that spanned highly conserved repetitive regions. The average accuracy of NaS reads reached 99.99% without losing the initial size of the input MinION® reads. Conclusions: We described NaS tool, a hybrid approach allowing the sequencing of microbial genomes using the MinION® device. Our method, based ideally on 20x and 50x of NaS and Illumina reads respectively, provides an efficient and cost-effective way of sequencing microbial or small eukaryotic genomes in a very short time even in small facilities. Moreover, we demonstrated that although the Oxford Nanopore technology is a relatively new sequencing technology, currently with a high error rate, it is already useful in the generation of high-quality genome assemblies. [ABSTRACT FROM AUTHOR]

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. The invertebrate iridoviruses 22 (IIV22) and 25 (IIV25) were originally isolated from a single sample of blackfly larva ( Simulium spp., order Diptera) collected from the Ystwyth river near Aberystwyth, Wales. Recently, the genomes of IIV22 (197.7 kbp) and IIV25 (204.8 kbp) were sequenced and reported. Here, we describe the complete genome sequence of IIV22A, a variant that was isolated from the same pool of virions collected from the blackfly larva from which the IIV22 virion genome originated. The IIV22A genome, 196.5 kbp, is smaller than IIV22. Nevertheless, it contains 7 supplementary putative ORFs. Its analysis enables evaluation of the degree of genomic polymorphisms within an IIV isolate. Despite the occurrence of this IIV variant with IIV22 and IIV25 in a single blackfly larva and the features of their DNA polymerase, we found no evidence of lateral genetic transfers between the genomes of these two IIV species. [ABSTRACT FROM AUTHOR]

Background: Many plant-pathogenic fungi have a tendency towards genome size expansion, mostly driven by increasing content of transposable elements (TEs). Through comparative and evolutionary genomics, five members of the Leptosphaeria maculans-Leptosphaeria biglobosa species complex (class Dothideomycetes, order Pleosporales), having different host ranges and pathogenic abilities towards cruciferous plants, were studied to infer the role of TEs on genome shaping, speciation, and on the rise of better adapted pathogens. Results: L. maculans ‘brassicae’, the most damaging species on oilseed rape, is the only member of the species complex to have a TE-invaded genome (32.5%) compared to the other members genomes (<4%). These TEs had an impact at the structural level by creating large TE-rich regions and are suspected to have been instrumental in chromosomal rearrangements possibly leading to speciation. TEs, associated with species-specific genes involved in disease process, also possibly had an incidence on evolution of pathogenicity by promoting translocations of effector genes to highly dynamic regions and thus tuning the regulation of effector gene expression in planta. Conclusions: Invasion of L. maculans ‘brassicae’ genome by TEs followed by bursts of TE activity allowed this species to evolve and to better adapt to its host, making this genome species a peculiarity within its own species complex as well as in the Pleosporales lineage. [ABSTRACT FROM AUTHOR]

Background: Metatranscriptomics is rapidly expanding our knowledge of gene expression patterns and pathway dynamics in natural microbial communities. However, to cope with the challenges of environmental sampling, various rRNA removal and cDNA synthesis methods have been applied in published microbial metatranscriptomic studies, making comparisons arduous. Whereas efficiency and biases introduced by rRNA removal methods have been relatively well explored, the impact of cDNA synthesis and library preparation on transcript abundance remains poorly characterized. The evaluation of potential biases introduced at this step is challenging for metatranscriptomic samples, where data analyses are complex, for example because of the lack of reference genomes. Results: Herein, we tested four cDNA synthesis and Illumina library preparation protocols on a simplified mixture of total RNA extracted from four bacterial species. In parallel, RNA from each microbe was tested individually. cDNA synthesis was performed on rRNA depleted samples using the TruSeq Stranded Total RNA Library Preparation, the SMARTer Stranded RNA-Seq, or the Ovation RNA-Seq V2 System. A fourth experiment was made directly from total RNA using the Encore Complete Prokaryotic RNA-Seq. The obtained sequencing data were analyzed for: library complexity and reproducibility; rRNA removal efficiency and bias; the number of genes detected; coverage uniformity; and the impact of protocols on expression biases. Significant variations, especially in organism representation and gene expression patterns, were observed among the four methods. TruSeq generally performed best, but is limited by its requirement of hundreds of nanograms of total RNA. The SMARTer method appears the best solution for smaller amounts of input RNA. For very low amounts of RNA, the Ovation System provides the only option; however, the observed biases emphasized its limitations for quantitative analyses. Conclusions: cDNA and library preparation methods may affect the outcome and interpretation of metatranscriptomic data. The most appropriate method should be chosen based on the available quantity of input RNA and the quantitative or non-quantitative objectives of the study. When low amounts of RNA are available, as in most metatranscriptomic studies, the SMARTer method seems to be the best compromise to obtain reliable results. This study emphasized the difficulty in comparing metatranscriptomic studies performed using different methods. [ABSTRACT FROM AUTHOR]

Background: Microbes drive the biogeochemistry that fuels the planet. Microbial viruses modulate their hosts directly through mortality and horizontal gene transfer, and indirectly by re-programming host metabolisms during infection. However, our ability to study these virus-host interactions is limited by methods that are low-throughput and heavily reliant upon the subset of organisms that are in culture. One way forward are culture-independent metagenomic approaches, but these novel methods are rarely rigorously tested, especially for studies of environmental viruses, air microbiomes, extreme environment microbiology and other areas with constrained sample amounts. Here we perform replicated experiments to evaluate Roche 454, Illumina HiSeq, and Ion Torrent PGM sequencing and library preparation protocols on virus metagenomes generated from as little as 10pg of DNA. Results: Using %G + C content to compare metagenomes, we find that (i) metagenomes are highly replicable, (ii) some treatment effects are minimal, e.g., sequencing technology choice has 6-fold less impact than varying input DNA amount, and (iii) when restricted to a limited DNA concentration (<1μg), changing the amount of amplification produces little variation. These trends were also observed when examining the metagenomes for gene function and assembly performance, although the latter more closely aligned to sequencing effort and read length than preparation steps tested. Among Illumina library preparation options, transposon-based libraries diverged from all others and adaptor ligation was a critical step for optimizing sequencing yields. Conclusions: These data guide researchers in generating systematic, comparative datasets to understand complex ecosystems, and suggest that neither varied amplification nor sequencing platforms will deter such efforts. [ABSTRACT FROM AUTHOR]

Background: Improving knowledge of biodiversity will benefit conservation biology, enhance bioremediation studies, and could lead to new medical treatments. However there is no standard approach to estimate and to compare the diversity of different environments, or to study its past, and possibly, future evolution. Presentation of the hypothesis: We argue that there are two conditions for significant progress in the identification and quantification of biodiversity. First, integrative metagenomic studies - aiming at the simultaneous examination (or even better at the integration) of observations about the elements, functions and evolutionary processes captured by the massive sequencing of multiple markers - should be preferred over DNA barcoding projects and over metagenomic projects based on a single marker. Second, such metagenomic data should be studied with novel inclusive network-based approaches, designed to draw inferences both on the many units and on the many processes present in the environments. Testing the hypothesis: We reached these conclusions through a comparison of the theoretical foundations of two molecular approaches seeking to assess biodiversity: metagenomics (mostly used on prokaryotes and protists) and DNA barcoding (mostly used on multicellular eukaryotes), and by pragmatic considerations of the issues caused by the 'species problem' in biodiversity studies. Implications of the hypothesis: Evolutionary gene networks reduce the risk of producing biodiversity estimates with limited explanatory power, biased either by unequal rates of LGT, or difficult to interpret due to (practical) problems caused by type I and type II grey zones. Moreover, these networks would easily accommodate additional (meta)transcriptomic and (meta)proteomic data. [ABSTRACT FROM AUTHOR]

Background: Whole-genome physical maps facilitate genome sequencing, sequence assembly, mapping of candidate genes, and the design of targeted genetic markers. An automated protocol was used to construct a Vitis vinifera 'Cabernet Sauvignon' physical map. The quality of the result was addressed with regard to the effect of high heterozygosity on the accuracy of contig assembly. Its usefulness for the genome-wide mapping of genes for disease resistance, which is an important trait for grapevine, was then assessed. Results: The physical map included 29,727 BAC clones assembled into 1,770 contigs, spanning 715,684 kbp, and corresponding to 1.5-fold the genome size. Map inflation was due to high heterozygosity, which caused either the separation of allelic BACs in two different contigs, or local mis-assembly in contigs containing BACs from the two haplotypes. Genetic markers anchored 395 contigs or 255,476 kbp to chromosomes. The fully automated assembly and anchorage procedures were validated by BAC-by-BAC blast of the end sequences against the grape genome sequence, unveiling 7.3% of chimerical contigs. The distribution across the physical map of candidate genes for non-host and host resistance, and for defence signalling pathways was then studied. NBS-LRR and RLK genes for host resistance were found in 424 contigs, 133 of them (32%) were assigned to chromosomes, on which they are mostly organised in clusters. Non-host and defence signalling genes were found in 99 contigs dispersed without a discernable pattern across the genome. Conclusion: Despite some limitations that interfere with the correct assembly of heterozygous clones into contigs, the 'Cabernet Sauvignon' physical map is a useful and reliable intermediary step between a genetic map and the genome sequence. This tool was successfully exploited for a quick mapping of complex families of genes, and it strengthened previous clues of co-localisation of major NBS-LRR clusters and disease resistance loci in grapevine. [ABSTRACT FROM AUTHOR]

Background: The island of Madagascar and surrounding volcanic and coralline islands are considered to form a biodiversity hotspot with large numbers of unique taxa. The origin of this endemic fauna can be explained by two different factors: vicariance or over-water-dispersal. Deciphering which factor explains the current distributional pattern of a given taxonomic group requires robust phylogenies as well as estimates of divergence times. The lineage of Indian Ocean scops-owls (Otus: Strigidae) includes six or seven species that are endemic to Madagascar and portions of the Comoros and Seychelles archipelagos; little is known about the species limits, biogeographic affinities and relationships to each other. In the present study, using DNA sequence data gathered from six loci, we examine the biogeographic history of the Indian Ocean scops-owls. We also compare the pattern and timing of colonization of the Indian Ocean islands by scops-owls with divergence times already proposed for other bird taxa. Results: Our analyses revealed that Indian Ocean islands scops-owls do not form a monophyletic assemblage: the Seychelles Otus insularis is genetically closer to the South-East Asian endemic O. sunia than to species from the Comoros and Madagascar. The Pemba Scops-owls O. pembaensis, often considered closely related to, if not conspecific with O. rutilus of Madagascar, is instead closely related to the African mainland O. senegalensis. Relationships among the Indian Ocean taxa from the Comoros and Madagascar are unresolved, despite the analysis of over 4000 bp, suggesting a diversification burst after the initial colonization event. We also highlight one case of putative backcolonization to the Asian mainland from an island ancestor (O. sunia). Our divergence date estimates, using a Bayesian relaxed clock method, suggest that all these events occurred during the last 3.6 myr; albeit colonization of the Indian Ocean islands were not synchronous, O. pembaensis diverged from O. senegalensis about 1.7 mya while species from Madagascar and the Comoro diverged from their continental sister-group about 3.6 mya. We highlight that our estimates coincide with estimates of diversification from other bird lineages. Conclusion: Our analyses revealed the occurrence of multiple synchronous colonization events of the Indian Ocean islands by scops-owls, at a time when faunistic exchanges involving Madagascar was common as a result of lowered sea-level that would have allowed the formation of steppingstone islands. Patterns of diversification that emerged from the scops-owls data are: 1) a star-like pattern concerning the order of colonization of the Indian Ocean islands and 2) the high genetic distinctiveness among all Indian Ocean taxa, reinforcing their recognition as distinct species. [ABSTRACT FROM AUTHOR]

Background: Circum-Antarctic waters harbour a rare example of a marine species flock -- the Notothenioid fish, most species of which are restricted to the continental shelf. It remains an open question as to how they survived Pleistocene climatic fluctuations characterised by repeated advances of continental glaciers as far as the shelf break that probably resulted in a loss of habitat for benthic organisms. Pelagic ecosystems, on the other hand, might have flourished during glacial maxima due to the northward expansion of Antarctic polar waters. In order to better understand the role of ecological traits in Quaternary climatic fluctuations, we performed demographic analyses of populations of four fish species from the tribe Trematominae, including both fully benthic and pelagic species using the mitochondrial cytochrome b gene and an intron from the nuclear S7 gene. Results: Nuclear and cytoplasmic markers showed differences in the rate and time of population expansions as well as the likely population structure. Neutrality tests suggest that such discordance comes from different coalescence dynamics of each marker, rather than from selective pressure. Demographic analyses based on intraspecific DNA diversity suggest a recent population expansion in both benthic species, dated by the cyt b locus to the last glacial cycle, whereas the population structure of pelagic feeders either did not deviate from a constant-size model or indicated that the onset of the major population expansion of these species by far predated those of the benthic species. Similar patterns were apparent even when comparing previously published data on other Southern Ocean organisms, but we observed considerable heterogeneity within both groups with regard to the onset of major demographic events and rates. Conclusion: Our data suggest benthic and pelagic species reacted differently to the Pleistocene ice-sheet expansions that probably significantly reduced the suitable habitat for benthic species. However, the asynchronous timing of major demographic events observed in different species within both "ecological guilds", imply that the species examined here may have different population and evolutionary histories, and that more species should be analysed in order to more precisely assess the role of life history in the response of organisms to climatic changes. [ABSTRACT FROM AUTHOR]

Background: Morphological data suggest that, unlike most other groups of marine organisms, scleractinian corals of the genus Stylophora are more diverse in the western Indian Ocean and in the Red Sea than in the central Indo-Pacific. However, the morphology of corals is often a poor predictor of their actual biodiversity: hence, we conducted a genetic survey of Stylophora corals collected in Madagascar, Okinawa, the Philippines and New Caledonia in an attempt to find out the true number of species in these various locations., Results: A molecular phylogenetic analysis of the mitochondrial ORF and putative control region concurs with a haploweb analysis of nuclear ITS2 sequences in delimiting three species among our dataset: species A and B are found in Madagascar whereas species C occurs in Okinawa, the Philippines and New Caledonia. Comparison of ITS1 sequences from these three species with data available online suggests that species C is also found on the Great Barrier Reef, in Malaysia, in the South China Sea and in Taiwan, and that a distinct species D occurs in the Red Sea. Shallow-water morphs of species A correspond to the morphological description of Stylophora madagascarensis, species B presents the morphology of Stylophora mordax, whereas species C comprises various morphotypes including Stylophora pistillata and Stylophora mordax., Conclusions: Genetic analysis of the coral genus Stylophora reveals species boundaries that are not congruent with morphological traits. Of the four hypotheses that may explain such discrepancy (phenotypic plasticity, morphological stasis, morphological convergence, and interspecific hybridization), the first two appear likely to play a role but the fourth one is rejected since mitochondrial and nuclear markers yield congruent species delimitations. The position of the root in our molecular phylogenies suggests that the center of origin of Stylophora is located in the western Indian Ocean, which probably explains why this genus presents a higher biodiversity in the westernmost part of its area of distribution than in the "Coral Triangle".

Background: Massively parallel DNA sequencing instruments are enabling the decoding of whole genomes at significantly lower cost and higher throughput than classical Sanger technology. Each of these technologies have been estimated to yield assemblies with more problematic features than the standard method. These problems are of a different nature depending on the techniques used. So, an appropriate mix of technologies may help resolve most difficulties, and eventually provide assemblies of high quality without requiring any Sanger-based input., Results: We compared assemblies obtained using Sanger data with those from different inputs from New Sequencing Technologies. The assemblies were systematically compared with a reference finished sequence. We found that the 454 GSFLX can efficiently produce high continuity when used at high coverage. The potential to enhance continuity by scaffolding was tested using 454 sequences from circularized genomic fragments. Finally, we explore the use of Solexa-Illumina short reads to polish the genome draft by implementing a technique to correct 454 consensus errors., Conclusion: High quality drafts can be produced for small genomes without any Sanger data input. We found that 454 GSFLX and Solexa/Illumina show great complementarity in producing large contigs and supercontigs with a low error rate.