dipSPAdes: Assembler for Highly Polymorphic Diploid Genomes

Uncovering the viral aetiology of undiagnosed acute febrile illness in Uganda using metagenomic sequencing

Viruses associated with acute febrile illness in Africa cause a spectrum of clinical disease from mild to life-threatening. Routine diagnostic methods are insufficient to identify all viral pathogens in this region. In this study, 1281 febrile Ugandan patients were prospectively recruited as part of the CDC-UVRI Acute Febrile Illness Study and pre-screened for common pathogens. 210/1281 undiagnosed samples, and 20 additional samples from viral outbreaks were subjected to metagenomic sequencing. Viral pathogens were identified in 44/230 (19%), including respiratory, hepatitis, blood-borne, gastrointestinal and vector-borne viruses. Importantly, one case of Crimean-Congo haemorrhagic fever and two cases each of Rift Valley fever, dengue and yellow fever were detected in 7/230 (3%) of cases. Le Dantec virus, last reported in 1969, was also identified in one patient. The presence of high-consequence and (re-)emerging viruses of public health concern highlights the need for enhanced population-based diagnostic surveillance in the African region.

Transposable Element Diversity and Activity Patterns in Neotropical Salamanders

Transposable elements (TEs) compose a substantial proportion of the largest eukaryotic genomes. TE diversity has been hypothesized to be negatively correlated with genome size, yet empirical demonstrations of such a relationship in a phylogenetic context are largely lacking. Furthermore, the most abundant type of TEs in genomes varies across groups, and it is not clear if there are patterns of TE activity consistent with genome size among different taxa with large genome sizes. We use low-coverage sequencing of 16 species of Neotropical salamanders, which vary ∼7-fold in genome size, to estimate TE relative abundance and diversity for each species. We also compare the divergence of copies of each TE superfamily to estimate patterns of TE activity in each species. We find a negative relationship between TE diversity and genome size, which is consistent with the hypothesis that either competition among TEs or reduced selection against ectopic recombination may result in lower diversity in the largest genomes. We also find divergent activity patterns in the largest versus the smallest genomes, suggesting that the history of TE activity may explain differences in genome size. Our results suggest that both TE diversity and relative abundance may be predictable, at least within taxonomic groups.

Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Complete chloroplast genomes of three wild perennial Hordeum species from Central Asia: genome structure, mutation hotspot, phylogenetic relationships, and comparative analysis

Hordeum L. is widely distributed in mountain or plateau of subtropical and warm temperate regions around the world. Three wild perennial Hordeum species, including H. bogdanii, H. brevisubulatum, and H. violaceum, have been used as forage and for grassland ecological restoration in high-altitude areas in recent years. To date, the degree of interspecies sequence variation in the three Hordeum species within existing gene pools is still not well-defined. Herein, we sequenced and assembled chloroplast (cp) genomes of the three species. The results revealed that the cp genome of H. bogdanii showed certain sequence variations compared with the cp genomes of the other two species (H. brevisubulatum and H. violaceum), and the latter two were characterized by a higher relative affinity. Parity rule 2 plot (PR2) analysis illuminated that most genes of all ten Hordeum species were concentrated in nucleotide T and G. Numerous single nucleotide polymorphism (SNP) and insertion/deletion (In/Del) events were detected in the three Hordeum species. A series of hotspots regions (tRNA-GGU ~ tRNA-GCA, tRNA-UGU ~ ndhJ, psbE ~ rps18, ndhF ~ tRNA-UAG, etc.) were identified by mVISTA procedures, and the five highly polymorphic genes (tRNA-UGC, tRNA-UAA, tRNA-UUU, tRNA-UAC, and ndhA) were proved by the nucleotide diversity (Pi). Although the distribution and existence of cp simple sequence repeats (cpSSRs) were predicted in the three Hordeum cp genomes, no rearrangement was found between them. A similar phenomenon has been found in the cp genome of the other seven Hordeum species, which has been published so far. In addition, evolutionary relationships were reappraised based on the currently reported cp genome of Hordeum L. This study offers a framework for gaining a better understanding of the evolutionary history of Hordeum species through the re-examination of their cp genomes, and by identifying highly polymorphic genes and hotspot regions that could provide important insights into the genetic diversity and differentiation of these species.

A target capture approach for phylogenomic analyses at multiple evolutionary timescales in rosewoods (Dalbergia spp.) and the legume family (Fabaceae)

Understanding the genetic changes associated with the evolution of biological diversity is of fundamental interest to molecular ecologists. The assessment of genetic variation at hundreds or thousands of unlinked genetic loci forms a sound basis to address questions ranging from micro- to macroevolutionary timescales, and is now possible thanks to advances in sequencing technology. Major difficulties are associated with (i) the lack of genomic resources for many taxa, especially from tropical biodiversity hotspots; (ii) scaling the numbers of individuals analysed and loci sequenced; and (iii) building tools for reproducible bioinformatic analyses of such data sets. To address these challenges, we developed target capture probes for genomic studies of the highly diverse, pantropically distributed and economically significant rosewoods (Dalbergia spp.), explored the performance of an overlapping probe set for target capture across the legume family (Fabaceae), and built the general purpose bioinformatic pipeline CaptureAl. Phylogenomic analyses of Malagasy Dalbergia species yielded highly resolved and well supported hypotheses of evolutionary relationships. Population genomic analyses identified differences between closely related species and revealed the existence of a potentially new species, suggesting that the diversity of Malagasy Dalbergia species has been underestimated. Analyses at the family level corroborated previous findings by the recovery of monophyletic subfamilies and many well-known clades, as well as high levels of gene tree discordance, especially near the root of the family. The new genomic and bioinformatic resources, including the Fabaceae1005 and Dalbergia2396 probe sets, will hopefully advance systematics and ecological genetics research in legumes, and promote conservation of the highly diverse and endangered Dalbergia rosewoods.

Karyon: a computational framework for the diagnosis of hybrids, aneuploids, and other nonstandard architectures in genome assemblies

BACKGROUND

Recent technological developments have made genome sequencing and assembly highly accessible and widely used. However, the presence in sequenced organisms of certain genomic features such as high heterozygosity, polyploidy, aneuploidy, heterokaryosis, or extreme compositional biases can challenge current standard assembly procedures and result in highly fragmented assemblies. Hence, we hypothesized that genome databases must contain a nonnegligible fraction of low-quality assemblies that result from such type of intrinsic genomic factors.

FINDINGS

Here we present Karyon, a Python-based toolkit that uses raw sequencing data and de novo genome assembly to assess several parameters and generate informative plots to assist in the identification of nonchanonical genomic traits. Karyon includes automated de novo genome assembly and variant calling pipelines. We tested Karyon by diagnosing 35 highly fragmented publicly available assemblies from 19 different Mucorales (Fungi) species.

CONCLUSIONS

Our results show that 10 (28.57%) of the assemblies presented signs of unusual genomic configurations, suggesting that these are common, at least for some lineages within the Fungi.

Paraphysoderma sedebokerense GlnS III Is Essential for the Infection of Its Host Haematococcus lacustris

Glutamine synthetase (GlnS) is a key enzyme in nitrogen metabolism. We investigated the effect of the GlnS inhibitor glufosinate on the infection of H. lacustris by the blastocladialean fungus P. sedebokerense, assuming that interfering with the host nitrogen metabolism will affect the success of the parasite. Complete inhibition of infection, which could be bypassed by the GlnS product glutamine, was observed at millimolar concentrations of glufosinate. However, this effect of glufosinate was attributed to its direct interaction with the blastoclad and not the host, which results in development and growth inhibition of the blastoclad. In our P. sedebokerense draft genome, we found that the sequence of GlnS is related to another fungal GlnS, type III, found in many poor known phyla of fungi, including Blastocladiomycota and Chytridiomycota, and absent in the main subkingdom of fungi, the Dikarya. We further tested the ability of the blastoclad to utilize nitrate and ammonia as inorganic nitrogen sources and glutamine for growth. We found that P. sedebokerense equally use ammonia and glutamine and use also nitrate, but with less efficiency. Altogether, our results show that GlnS type III is mandatory for the development and growth of P. sedebokerense and could be an efficient target to develop strategies for the control of the fungal parasite of H. lacustris.

FrogCap: A modular sequence capture probe-set for phylogenomics and population genetics for all frogs, assessed across multiple phylogenetic scales

Despite the prevalence of high-throughput sequencing in phylogenetics, many relationships remain difficult to resolve because of conflicting signal among genomic regions. Selection of different types of molecular markers from different genomic regions is required to overcome these challenges. For evolutionary studies in frogs, we introduce the publicly available FrogCap suite of genomic resources, which is a large collection of ~15,000 markers that unifies previous genetic sequencing efforts. FrogCap is designed to be modular, such that subsets of markers and SNPs can be selected based on the desired phylogenetic scale. FrogCap uses a variety of marker types that include exons and introns, ultraconserved elements, and previously sequenced Sanger markers, which span up to 10,000 bp in alignment lengths; in addition, we demonstrate potential for SNP-based analyses. We tested FrogCap using 121 samples distributed across five phylogenetic scales, comparing probes designed using a consensus- or exemplar genome-based approach. Using the consensus design is more resilient to issues with sensitivity, specificity, and missing data than picking an exemplar genome sequence. We also tested the impact of different bait kit sizes (20,020 vs. 40,040) on depth of coverage and found triple the depth for the 20,020 bait  kit. We observed sequence capture success (i.e., missing data, sequenced markers/bases, marker length, and informative sites) across phylogenetic scales. The incorporation of different marker types is effective for deep phylogenetic relationships and shallow population genetics studies. Having demonstrated FrogCap's utility and modularity, we conclude that these new resources are efficacious for high-throughput sequencing projects across variable timescales.

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.

Viral capture sequencing detects unexpected viruses in the cerebrospinal fluid of adults with meningitis

OBJECTIVES

Many patients with meningitis have no aetiology identified leading to unnecessary antimicrobials and prolonged hospitalisation. We used viral capture sequencing to identify possible pathogenic viruses in adults with community-acquired meningitis.

METHODS

Cerebrospinal fluid (CSF) from 73 patients was tested by VirCapSeq-VERT, a probe set designed to capture viral targets using high throughput sequencing. Patients were categorised as suspected viral meningitis - CSF pleocytosis, no pathogen identified (n = 38), proven viral meningitis - CSF pleocytosis with a pathogen identified (n = 15) or not meningitis - no CSF pleocytosis (n = 20).

RESULTS

VirCapSeq-VERT detected virus in the CSF of 16/38 (42%) of those with suspected viral meningitis, including twelve individual viruses. A potentially clinically relevant virus was detected in 9/16 (56%). Unexpectedly Toscana virus, rotavirus and Saffold virus were detected and assessed to be potential causative agents.

CONCLUSION

VirCapSeq-VERT increases the probability of detecting a virus. Using this agnostic approach we identified Toscana virus and, for the first time in adults, rotavirus and Saffold virus, as potential causative agents in adult meningitis. Further work is needed to determine the prevalence of atypical viral candidates as well as the clinical impact of using sequencing methods in real time. This knowledge can help to reduce antimicrobial use and hospitalisations leading to both patient and health system benefits.

Genome analysis of Candida subhashii reveals its hybrid nature and dual mitochondrial genome conformations

Candida subhashii belongs to the CUG-Ser clade, a group of phylogenetically closely related yeast species that includes some human opportunistic pathogens, such as Candida albicans. Despite being present in the environment, C. subhashii was initially described as the causative agent of a case of peritonitis. Considering the relevance of whole-genome sequencing and analysis for our understanding of genome evolution and pathogenicity, we sequenced, assembled and annotated the genome of C. subhashii type strain. Our results show that C. subhashii presents a highly heterozygous genome and other signatures that point to a hybrid ancestry. The presence of functional pathways for assimilation of hydroxyaromatic compounds goes in line with the affiliation of this yeast with soil microbial communities involved in lignin decomposition. Furthermore, we observed that different clones of this strain may present circular or linear mitochondrial DNA. Re-sequencing and comparison of strains with differential mitochondrial genome topology revealed five candidate genes potentially associated with this conformational change: MSK1, SSZ1, ALG5, MRPL9 and OYE32.

Shared evolutionary footprints suggest mitochondrial oxidative damage underlies multiple complex I losses in fungi

Oxidative phosphorylation is among the most conserved mitochondrial pathways. However, one of the cornerstones of this pathway, the multi-protein complex NADH : ubiquinone oxidoreductase (complex I) has been lost multiple independent times in diverse eukaryotic lineages. The causes and consequences of these convergent losses remain poorly understood. Here, we used a comparative genomics approach to reconstruct evolutionary paths leading to complex I loss and infer possible evolutionary scenarios. By mining available mitochondrial and nuclear genomes, we identified eight independent events of mitochondrial complex I loss across eukaryotes, of which six occurred in fungal lineages. We focused on three recent loss events that affect closely related fungal species, and inferred genomic changes convergently associated with complex I loss. Based on these results, we predict novel complex I functional partners and relate the loss of complex I with the presence of increased mitochondrial antioxidants, higher fermentative capabilities, duplications of alternative dehydrogenases, loss of alternative oxidases and adaptation to antifungal compounds. To explain these findings, we hypothesize that a combination of previously acquired compensatory mechanisms and exposure to environmental triggers of oxidative stress (such as hypoxia and/or toxic chemicals) induced complex I loss in fungi.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Gigantic Genomes Provide Empirical Tests of Transposable Element Dynamics Models

Transposable elements (TEs) are a major determinant of eukaryotic genome size. The collective properties of a genomic TE community reveal the history of TE/host evolutionary dynamics and impact present-day host structure and function, from genome to organism levels. In rare cases, TE community/genome size has greatly expanded in animals, associated with increased cell size and changes to anatomy and physiology. Here, we characterize the TE landscape of the genome and transcriptome in an amphibian with a giant genome — the caecilianIchthyophis bannanicus, which we show has a genome size of 12.2 Gb. Amphibians are an important model system because the clade includes independent cases of genomic gigantism. The I. bannanicus genome differs compositionally from other giant amphibian genomes, but shares a low rate of ectopic recombination-mediated deletion. We examine TE activity using expression and divergence plots; TEs account for 15% of somatic transcription, and most superfamilies appear active. We quantify TE diversity in the caecilian, as well as other vertebrates with a range of genome sizes, using diversity indices commonly applied in community ecology. We synthesize previous models that integrate TE abundance, diversity, and activity, and test whether the caecilian meets model predictions for genomes with high TE abundance. We propose thorough, consistent characterization of TEs to strengthen future comparative analyses. Such analyses will ultimately be required to reveal whether the divergent TE assemblages found across convergent gigantic genomes reflect fundamental shared features of TE/host genome evolutionary dynamics.

HaploTypo: a variant-calling pipeline for phased genomes

Summary

An increasing number of phased (i.e. with resolved haplotypes) reference genomes are available. However, the most genetic variant calling tools do not explicitly account for haplotype structure. Here, we present HaploTypo, a pipeline tailored to resolve haplotypes in genetic variation analyses. HaploTypo infers the haplotype correspondence for each heterozygous variant called on a phased reference genome.

Availability and implementation

HaploTypo is implemented in Python 2.7 and Python 3.5, and is freely available at https://github.com/gabaldonlab/haplotypo, and as a Docker image.

Supplementary information

Supplementary data are available at Bioinformatics online.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Pathogenic Allodiploid Hybrids of Aspergillus Fungi

Interspecific hybridization substantially alters genotypes and phenotypes and can give rise to new lineages. Hybrid isolates that differ from their parental species in infection-relevant traits have been observed in several human-pathogenic yeasts and plant-pathogenic filamentous fungi but have yet to be found in human-pathogenic filamentous fungi. We discovered 6 clinical isolates from patients with aspergillosis originally identified as Aspergillus nidulans (section Nidulantes) that are actually allodiploid hybrids formed by the fusion of Aspergillus spinulosporus with an unknown close relative of Aspergillus quadrilineatus, both in section Nidulantes. Evolutionary genomic analyses revealed that these isolates belong to Aspergillus latus, an allodiploid hybrid species. Characterization of diverse infection-relevant traits further showed that A. latus hybrid isolates are genomically and phenotypically heterogeneous but also differ from A. nidulans, A. spinulosporus, and A. quadrilineatus. These results suggest that allodiploid hybridization contributes to the genomic and phenotypic diversity of filamentous fungal pathogens of humans.

Overlap graph-based generation of haplotigs for diploids and polyploids

MOTIVATION

Haplotype-aware genome assembly plays an important role in genetics, medicine and various other disciplines, yet generation of haplotype-resolved de novo assemblies remains a major challenge. Beyond distinguishing between errors and true sequential variants, one needs to assign the true variants to the different genome copies. Recent work has pointed out that the enormous quantities of traditional NGS read data have been greatly underexploited in terms of haplotig computation so far, which reflects that methodology for reference independent haplotig computation has not yet reached maturity.

RESULTS

We present POLYploid genome fitTEr (POLYTE) as a new approach to de novo generation of haplotigs for diploid and polyploid genomes of known ploidy. Our method follows an iterative scheme where in each iteration reads or contigs are joined, based on their interplay in terms of an underlying haplotype-aware overlap graph. Along the iterations, contigs grow while preserving their haplotype identity. Benchmarking experiments on both real and simulated data demonstrate that POLYTE establishes new standards in terms of error-free reconstruction of haplotype-specific sequence. As a consequence, POLYTE outperforms state-of-the-art approaches in various relevant aspects, where advantages become particularly distinct in polyploid settings.

AVAILABILITY AND IMPLEMENTATION

POLYTE is freely available as part of the HaploConduct package at https://github.com/HaploConduct/HaploConduct, implemented in Python and C++.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Complete Chloroplast Genome Sequencing and Comparative Analysis of Reed Canary Grass (Phalaris arundinacea) and Hardinggrass (P. aquatica)

There are 22 species in the Phalaris genera that distribute almost all over the temperate regions of the world. Among them, reed canary grass (Phalaris arundinacea, tetraploid and hexaploid) and hardinggrass (P. aquatica, tetraploid) have been long cultivated as forage grass and have received attention as bio-energy materials in recent years. We aimed to facilitate inter-species/ploidies comparisons, and to illuminate the degree of sequence variation within existing gene pools, chloroplast (cp) genomes of three Phalaris cytotypes (P. aquatica/4x, P. arundinacea/4x and P. arundinacea/6x) were sequenced and assembled. The result indicated that certain sequence variations existed between the cp genomes of P. arundinacea and P. aquatica. Several hotspot regions (atpI~atpH, trnT-UGU~ndhJ, rbcL~psaI, and ndhF~rpl32) were found, and variable genes (infA, psaI, psbK, etc.) were detected. SNPs (single nucleotide polymorphisms) and/or indels (insertions and deletions) were confirmed by the high Ka/Ks and Pi value. Furthermore, distribution and presence of cp simple sequence repeats (cpSSRs) were identified in the three Phalaris cp genomes, although little difference was found between hexaploid and tetraploid P. arundinacea, and no rearrangement was detected among the three Phalaris cp genomes. The evolutionary relationship and divergent time among these species were discussed. The RNA-seq revealed several differentially expressed genes (DEGs), among which psaA, psaB, and psbB related to leaf color were further verified by leaf color differences.

Whole Genome Sequences of 23 Species from the Drosophila montium Species Group (Diptera: Drosophilidae): A Resource for Testing Evolutionary Hypotheses

Large groups of species with well-defined phylogenies are excellent systems for testing evolutionary hypotheses. In this paper, we describe the creation of a comparative genomic resource consisting of 23 genomes from the species-rich Drosophila montium species group, 22 of which are presented here for the first time. The montium group is well-positioned for clade genomics. Within the montium clade, evolutionary distances are such that large numbers of sequences can be accurately aligned while also recovering strong signals of divergence; and the distance between the montium group and D. melanogaster is short enough so that orthologous sequence can be readily identified. All genomes were assembled from a single, small-insert library using MaSuRCA, before going through an extensive post-assembly pipeline. Estimated genome sizes within the montium group range from 155 Mb to 223 Mb (mean = 196 Mb). The absence of long-distance information during the assembly process resulted in fragmented assemblies, with the scaffold NG50s varying widely based on repeat content and sample heterozygosity (min = 18 kb, max = 390 kb, mean = 74 kb). The total scaffold length for most assemblies is also shorter than the estimated genome size, typically by 5-15%. However, subsequent analysis showed that our assemblies are highly complete. Despite large differences in contiguity, all assemblies contain at least 96% of known single-copy Dipteran genes (BUSCOs, n = 2,799). Similarly, by aligning our assemblies to the D. melanogaster genome and remapping coordinates for a large set of transcriptional enhancers (n = 3,457), we showed that each montium assembly contains orthologs for at least 91% of D. melanogaster enhancers. Importantly, the genic and enhancer contents of our assemblies are comparable to that of far more contiguous Drosophila assemblies. The alignment of our own D. serrata assembly to a previously published PacBio D. serrata assembly also showed that our longest scaffolds (up to 1 Mb) are free of large-scale misassemblies. Our genome assemblies are a valuable resource that can be used to further resolve the montium group phylogeny; study the evolution of protein-coding genes and cis-regulatory sequences; and determine the genetic basis of ecological and behavioral adaptations.

Parachute geckos free fall into synonymy: Gekko phylogeny, and a new subgeneric classification, inferred from thousands of ultraconserved elements

Recent phylogenetic studies of gekkonid lizards have revealed unexpected, widespread paraphyly and polyphyly among genera, unclear generic boundaries, and a tendency towards the nesting of taxa exhibiting specialized, apomorphic morphologies within geographically widespread "generalist" clades. This is especially true in Australasia, where monophyly of Gekko proper has been questioned with respect to phenotypically ornate flap-legged geckos of the genus Luperosaurus, the Philippine false geckos of the genus Pseudogekko, and even the elaborately "derived" parachute geckos of the genus Ptychozoon. Here we employ sequence capture targeting 5060 ultraconserved elements (UCEs) to infer phylogenomic relationships among 42 representative ingroup gekkonine lizard taxa. We analyze multiple datasets of varying degrees of completeness (10, 50, 75, 95, and 100 percent complete with 4715, 4051, 3376, 2366, and 772 UCEs, respectively) using concatenated maximum likelihood and multispecies coalescent methods. Our sampling scheme addresses four persistent systematic questions in this group: (1) Are Luperosaurus and Ptychozoon monophyletic, and are any of these named species truly nested within Gekko? (2) Are prior phylogenetic estimates of Sulawesi's L. iskandari as the sister taxon to Melanesian G. vittatus supported by our genome-scale dataset? (3) Is the high-elevation L. gulat of Palawan Island correctly placed within Gekko? (4) And, finally, where do the enigmatic taxa P. rhacophorus and L. browni fall in a higher-level gekkonid phylogeny? We resolve these issues; confirm with strong support some previously inferred findings (placement of Ptychozoon taxa within Gekko; the sister taxon relationship between L. iskandari and G. vittatus); resolve the systematic position of unplaced taxa (L. gulat, and L. browni); and transfer L. iskandari, L. gulat, L. browni, and all members of the genus Ptychozoon to the genus Gekko. Our unexpected and novel systematic inference of the placement of Ptychozoon rhacophorus suggests that this species is not grouped with Ptychozoon or even Luperosaurus (as previously expected) but may, in fact, be most closely related to several Indochinese species of Gekko. With our resolved and strongly supported phylogeny, we present a new classification emphasizing the most inclusive, original generic name (Gekko) for these ~60 taxa, arranged into seven subgenera.

Metagenomic next-generation sequencing aids the diagnosis of viral infections in febrile returning travellers

OBJECTIVES

Travel-associated infections are challenging to diagnose because of the broad spectrum of potential aetiologies. As a proof-of-principle study, we used MNGS to identify viral pathogens in clinical samples from returning travellers in a single center to explore its suitability as a diagnostic tool.

METHODS

Plasma samples from 40 returning travellers presenting with a fever of ≥38°C were sequenced using MNGS on the Illumina MiSeq platform and compared with standard-of-care diagnostic assays.

RESULTS

In total, 11/40 patients were diagnosed with a viral infection. Standard of care diagnostics revealed 5 viral infections using plasma samples; dengue virus 1 (n = 2), hepatitis E (n = 1), Ebola virus (n = 1) and hepatitis A (n = 1), all of which were detected by MNGS. Three additional patients with Chikungunya virus (n = 2) and mumps virus were diagnosed by MNGS only. Respiratory infections detected by nasal/throat swabs only were not detected by MNGS of plasma. One patient had infection with malaria and mumps virus during the same admission.

CONCLUSIONS

MNGS analysis of plasma samples improves the sensitivity of diagnosis of viral infections and has potential as an all-in-one diagnostic test. It can be used to identify infections that have not been considered by the treating physician, co-infections and new or emerging pathogens.

SUMMARY

Next generation sequencing (NGS) has potential as an all-in-one diagnostic test. In this study we used NGS to diagnose returning travellers with acute febrile illness in the UK, highlighting cases where the diagnosis was missed using standard methods.

Revisiting the Phylogenetic History of Helminths Through Genomics, the Case of the New Echinococcus oligarthrus Genome

The first parasitic helminth genome sequence was published in 2007; since then, only ∼200 genomes have become available, most of them being draft assemblies. Nevertheless, despite the medical and economical global impact of helminthic infections, parasite genomes in public databases are underrepresented. Recently, through an integrative approach involving morphological, genetic, and ecological aspects, we have demonstrated that the complete life cycle of Echinococcus oligarthrus (Cestoda: Taeniidae) is present in South America. The neotropical E. oligarthrus parasite is capable of developing in any felid species and producing human infections. Neotropical echinococcosis is poorly understood yet and requires a complex medical examination to provide the appropriate intervention. Only a few cases of echinococcosis have been unequivocally identified and reported as a consequence of E. oligarthrus infections. Regarding phylogenetics, the analyses of mitogenomes and nuclear datasets have resulted in discordant topologies, and there is no unequivocal taxonomic classification of Echinococcus species so far. In this work, we sequenced and assembled the genome of E. oligarthrus that was isolated from agoutis (Dasyprocta azarae) naturally infected and performed the first comparative genomic study of a neotropical Echinococcus species. The E. oligarthrus genome assembly consisted of 86.22 Mb which showed ∼90% identity and 76.3% coverage with Echinococcus multilocularis and contained the 85.0% of the total expected genes. Genetic variants analysis of whole genome revealed a higher rate of intraspecific genetic variability (23,301 SNPs; 0.22 SNPs/kb) rather than for the genomes of E. multilocularis and Echinococcus canadensis G7 but lower with respect to Echinococcus granulosus G1. Comparative genomics against E. multilocularis, E. granulosus G1, and E. canadensis G7 revealed 38,762, 125,147, and 170,049 homozygous polymorphic sites, respectively, indicating a higher genetic distance between E. oligarthrus and E. granulosus sensu lato species. The SNP distribution in chromosomes revealed a higher SNP density in the longest chromosomes. Phylogenetic analysis using whole-genome SNPs demonstrated that E. oligarthrus is one of the basal species of the genus Echinococcus and is phylogenetically closer to E. multilocularis. This work sheds light on the Echinococcus phylogeny and settles the basis to study sylvatic Echinococcus species and their developmental evolutionary features.

Coordinated host-pathogen transcriptional dynamics revealed using sorted subpopulations and single macrophages infected with Candida albicans

The outcome of fungal infections depends on interactions with innate immune cells. Within a population of macrophages encountering Candida albicans, there are distinct host-pathogen trajectories; however, little is known about the molecular heterogeneity that governs these fates. Here we developed an experimental system to separate interaction stages and single macrophage cells infected with C. albicans from uninfected cells and assessed transcriptional variability in the host and fungus. Macrophages displayed an initial up-regulation of pathways involved in phagocytosis and proinflammatory response after C. albicans exposure that declined during later time points. Phagocytosed C. albicans shifted expression programs to survive the nutrient poor phagosome and remodeled the cell wall. The transcriptomes of single infected macrophages and phagocytosed C. albicans displayed a tightly coordinated shift in gene expression co-stages and revealed expression bimodality and differential splicing that may drive infection outcome. This work establishes an approach for studying host-pathogen trajectories to resolve heterogeneity in dynamic populations.

The outcomes of the interactions between individual host cells and pathogens are heterogeneous. Here, the authors assess transcriptional variability in both host and pathogen during infection of macrophages with the fungus Candida albicans, using sorted subpopulations and single macrophages.

High-Quality Draft Genome Sequence of the Microcolonial Black Fungus Aeminium ludgeri DSM 106916

Aeminium ludgeri is an extremotolerant microcolonial black fungus isolated from a biodeteriorated limestone art piece in the Old Cathedral of Coimbra, Portugal (a UNESCO World Heritage Site). The high-quality draft genome sequence of Aeminium ludgeri presented here represents the first sequenced genome for both the recently described fungal family Aeminiaceae and the genus Aeminium.

Aeminium ludgeri is an extremotolerant microcolonial black fungus isolated from a biodeteriorated limestone art piece in the Old Cathedral of Coimbra, Portugal (a UNESCO World Heritage Site). The high-quality draft genome sequence of Aeminium ludgeri presented here represents the first sequenced genome for both the recently described fungal family Aeminiaceae and the genus Aeminium.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies

Background

Recent developments in third-gen long read sequencing and diploid-aware assemblers have resulted in the rapid release of numerous reference-quality assemblies for diploid genomes. However, assembly of highly heterozygous genomes is still problematic when regional heterogeneity is so high that haplotype homology is not recognised during assembly. This results in regional duplication rather than consolidation into allelic variants and can cause issues with downstream analysis, for example variant discovery, or haplotype reconstruction using the diploid assembly with unpaired allelic contigs.

Results

A new pipeline—Purge Haplotigs—was developed specifically for third-gen sequencing-based assemblies to automate the reassignment of allelic contigs, and to assist in the manual curation of genome assemblies. The pipeline uses a draft haplotype-fused assembly or a diploid assembly, read alignments, and repeat annotations to identify allelic variants in the primary assembly. The pipeline was tested on a simulated dataset and on four recent diploid (phased) de novo assemblies from third-generation long-read sequencing, and compared with a similar tool. After processing with Purge Haplotigs, haploid assemblies were less duplicated with minimal impact on genome completeness, and diploid assemblies had more pairings of allelic contigs.

Conclusions

Purge Haplotigs improves the haploid and diploid representations of third-gen sequencing based genome assemblies by identifying and reassigning allelic contigs. The implementation is fast and scales well with large genomes, and it is less likely to over-purge repetitive or paralogous elements compared to alignment-only based methods. The software is available at https://bitbucket.org/mroachawri/purge_haplotigs under a permissive MIT licence.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2485-7) contains supplementary material, which is available to authorized users.

SKESA: strategic k-mer extension for scrupulous assemblies

SKESA is a DeBruijn graph-based de-novo assembler designed for assembling reads of microbial genomes sequenced using Illumina. Comparison with SPAdes and MegaHit shows that SKESA produces assemblies that have high sequence quality and contiguity, handles low-level contamination in reads, is fast, and produces an identical assembly for the same input when assembled multiple times with the same or different compute resources. SKESA has been used for assembling over 272,000 read sets in the Sequence Read Archive at NCBI and for real-time pathogen detection. Source code for SKESA is freely available at https://github.com/ncbi/SKESA/releases .

Trichoplax genomes reveal profound admixture and suggest stable wild populations without bisexual reproduction

The phylum Placozoa officially consists of only a single described species, Trichoplax adhaerens, although several lineages can be separated by molecular markers, geographical distributions and environmental demands. The placozoan 16S haplotype H2 (Trichoplax sp. H2) is the most robust and cosmopolitan lineage of placozoans found to date. In this study, its genome was found to be distinct but highly related to the Trichoplax adhaerens reference genome, for remarkably unique reasons. The pattern of variation and allele distribution between the two lineages suggests that both originate from a single interbreeding event in the wild, dating back at least several decades ago, and both seem not to have engaged in sexual reproduction since. We conclude that populations of certain placozoan haplotypes remain stable for long periods without bisexual reproduction. Furthermore, allelic variation within and between the two Trichoplax lineages indicates that successful bisexual reproduction between related placozoan lineages might serve to either counter accumulated negative somatic mutations or to cope with changing environmental conditions. On the other hand, enrichment of neutral or beneficial somatic mutations by vegetative reproduction, combined with rare sexual reproduction, could instantaneously boost genetic variation, generating novel ecotypes and eventually species.

Comparative genomics and the nature of placozoan species

Placozoans are a phylum of nonbilaterian marine animals currently represented by a single described species, Trichoplax adhaerens, Schulze 1883. Placozoans arguably show the simplest animal morphology, which is identical among isolates collected worldwide, despite an apparently sizeable genetic diversity within the phylum. Here, we use a comparative genomics approach for a deeper appreciation of the structure and causes of the deeply diverging lineages in the Placozoa. We generated a high-quality draft genome of the genetic lineage H13 isolated from Hong Kong and compared it to the distantly related T. adhaerens. We uncovered substantial structural differences between the two genomes that point to a deep genomic separation and provide support that adaptation by gene duplication is likely a crucial mechanism in placozoan speciation. We further provide genetic evidence for reproductively isolated species and suggest a genus-level difference of H13 to T. adhaerens, justifying the designation of H13 as a new species, Hoilungia hongkongensis nov. gen., nov. spec., now the second described placozoan species and the first in a new genus. Our multilevel comparative genomics approach is, therefore, likely to prove valuable for species distinctions in other cryptic microscopic animal groups that lack diagnostic morphological characters, such as some nematodes, copepods, rotifers, or mites.

Draft Genome Sequence of a Highly Heterozygous Yeast Strain from the Metschnikowia pulcherrima Subclade, UCD127

Metschnikowia strain UCD127 was isolated from soil in Ireland and sequenced. It is a highly heterozygous diploid strain with 385,000 single nucleotide polymorphisms (SNPs). Its ribosomal DNA has the highest similarity to that of M. chrysoperlae, but its ACT1 and TEF1 loci and mitochondrial genome show affinity to those of M. fructicola, whose genome is significantly larger.

Opisthorchis viverrini Draft Genome - Biomedical Implications and Future Avenues

Opisthorchiasis is a neglected tropical disease of major proportion, caused by the carcinogenic, Asian liver fluke, Opisthorchis viverrini. This hepatobiliary disease is known to be associated with malignant cancer (cholangiocarcinoma, CCA) and affects millions of people in Southeast Asia. No vaccine is available, and only one drug (praziquantel) is routinely employed against the parasite. Despite technological advances, little is known about the molecular biology of the fluke itself and the disease complex that it causes in humans. The advent of high-throughput nucleic acid sequencing and bioinformatic technologies is enabling researchers to gain global insights into the molecular pathways and processes in parasites. The principal aims of this chapter are to (1) review molecular research of O. viverrini and opisthorchiasis; (2) provide an account of recent advances in the sequencing and characterization of the genome and transcriptomes of O. viverrini; (3) describe the complex life of this worm in the biliary system of the definitive (human) host and how the fluke interacts with this host and causes disease at the molecular level; (4) discuss the implications of systems biological research and (5) consider how progress in genomics and informatics might enable explorations of O. viverrini and related worms and the discovery of new interventions against opisthorchiasis and CCA.

Genome Sequence of the Brown Rot Fungal Pathogen Monilinia laxa

Monilinia laxa (phylum Ascomycota) is a plant pathogen responsible for the brown rot blossom blight disease in stone fruit trees of the Rosaceae family, such as apricots. We report here the genome sequence of strain 8L of this species, which was assembled into 618 scaffolds, having a total size of 40.799 Mb and encoding 9,567 unique protein-coding genes.

A 200K SNP chip reveals a novel Pacific salmon louse genotype linked to differential efficacy of emamectin benzoate

Antiparasitic drugs such as emamectin benzoate (EMB) are relied upon to reduce the parasite load, particularly of the sea louse Lepeophtheirus salmonis, on farmed salmon. The decline in EMB treatment efficacy for this purpose is an important issue for salmon producers around the world, and particularly for those in the Atlantic Ocean where widespread EMB tolerance in sea lice is recognized as a significant problem. Salmon farms in the Northeast Pacific Ocean have not historically experienced the same issues with treatment efficacy, possibly due to the relatively large population of endemic salmonid hosts that serve to both redistribute surviving lice and dilute populations potentially under selection by introducing naïve lice to farms. Frequent migration of lice among farmed and wild hosts should limit the effect of farm-specific selection pressures on changes to the overall allele frequencies of sea lice in the Pacific Ocean. A previous study using microsatellites examined L. salmonis oncorhynchi from 10 Pacific locations from wild and farmed hosts and found no population structure. Recently however, a farm population of sea lice was detected where EMB bioassay exposure tolerance was abnormally elevated. In response, we have developed a Pacific louse draft genome that complements the previously-released Atlantic louse sequence. These genomes were combined with whole-genome re-sequencing data to design a highly sensitive 201,279 marker SNP array applicable for both subspecies (90,827 validated Pacific loci; 153,569 validated Atlantic loci). Notably, kmer spectrum analysis of the re-sequenced samples indicated that Pacific lice exhibit a large within-individual heterozygosity rate (average of 1 in every 72 bases) that is markedly higher than that of Atlantic individuals (1 in every 173 bases). The SNP chip was used to produce a high-density map for Atlantic sea louse linkage group 5 that was previously shown to be associated with EMB tolerance in Atlantic lice. Additionally, 478 Pacific louse samples from farmed and wild hosts obtained between 2005 and 2014 were also genotyped on the array. Clustering analysis allowed us to detect the apparent emergence of an otherwise rare genotype at a high frequency among the lice collected from two farms in 2013 that had reported elevated EMB tolerance. This genotype was not observed in louse samples collected from the same farm in 2010, nor in any lice sampled from other locations prior to 2013. However, this genotype was detected at low frequencies in louse samples from farms in two locations reporting elevated EMB tolerance in 2014. These results suggest that a rare genotype present in Pacific lice may be locally expanded in farms after EMB treatment. Supporting this hypothesis, 437 SNPs associated with this genotype were found to be in a region of linkage group 5 that overlaps the region associated with EMB resistance in Atlantic lice. Finally, five of the top diagnostic SNPs within this region were used to screen lice that had been subjected to an EMB survival assay, revealing a significant association between these SNPs and EMB treatment outcome. To our knowledge this work is the first report to identify a genetic link to altered EMB efficacy in L. salmonis in the Pacific Ocean.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

Epichloë hybrida, sp. nov., an emerging model system for investigating fungal allopolyploidy

Endophytes of the genus Epichloë (Clavicipitaceae, Ascomycota) frequently occur within cool-season grasses and form interactions with their hosts that range from mutualistic to antagonistic. Many Epichloë species have arisen via interspecific hybridization, resulting in species with two or three subgenomes that retain all or nearly all of their original parental genomes, a process termed allopolyploidization. Here, we characterize Epichloë hybrida, sp. nov., a mutualistic species that has increasingly become a model system for investigating allopolyploidy in fungi. The Epichloë species so far identified as the closest known relatives of the two progenitors of E. hybrida are E. festucae var. lolii and E. typhina. We confirm that the nuclear genome of E. hybrida contains two homeologs of most protein-coding genes from E. festucae and E. typhina, with genome-wide gene expression analysis indicating a slight bias in overall gene expression from the E. typhina subgenome. Mitochondrial DNA is detectable only from E. festucae, whereas ribosomal DNA is detectable only from E. typhina. Inheriting ribosomal DNA from just one parent might be expected to preferentially favor interactions with ribosomal proteins from the same parent, but we find that ribosomal protein genes from both parental subgenomes are nearly all expressed equally in E. hybrida. Finally, we provide a comprehensive set of resources for this model system that are intended to facilitate further study of fungal hybridization by other researchers.

HaploMerger2: rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly

SUMMARY

De novo assembly is a difficult issue for heterozygous diploid genomes. The advent of high-throughput short-read and long-read sequencing technologies provides both new challenges and potential solutions to the issue. Here, we present HaploMerger2 (HM2), an automated pipeline for rebuilding both haploid sub-assemblies from the polymorphic diploid genome assembly. It is designed to work on pre-existing diploid assemblies, which are typically created by using de novo assemblers. HM2 can process any diploid assemblies, but it is especially suitable for diploid assemblies with high heterozygosity (≥3%), which can be difficult for other tools. This pipeline also implements flexible and sensitive assembly error detection, a hierarchical scaffolding procedure and a reliable gap-closing method for haploid sub-assemblies. Using HM2, we demonstrate that two haploid sub-assemblies reconstructed from a real, highly-polymorphic diploid assembly show greatly improved continuity.

AVAILABILITY AND IMPLEMENTATION

Source code, executables and the testing dataset are freely available at https://github.com/mapleforest/HaploMerger2/releases/.

CONTACT

hshengf2@mail.sysu.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Draft Genome Sequence of Dermatophagoides pteronyssinus, the European House Dust Mite

Dermatophagoides pteronyssinus is the European dust mite and a major source of human allergens. Here, we present the first draft genome sequence of the mite, as well as the ab initio gene prediction and functional analyses that will facilitate comparative genomic analyses with other mite species.

metaSPAdes: a new versatile metagenomic assembler

While metagenomics has emerged as a technology of choice for analyzing bacterial populations, the assembly of metagenomic data remains challenging, thus stifling biological discoveries. Moreover, recent studies revealed that complex bacterial populations may be composed from dozens of related strains, thus further amplifying the challenge of metagenomic assembly. metaSPAdes addresses various challenges of metagenomic assembly by capitalizing on computational ideas that proved to be useful in assemblies of single cells and highly polymorphic diploid genomes. We benchmark metaSPAdes against other state-of-the-art metagenome assemblers and demonstrate that it results in high-quality assemblies across diverse data sets.

The Echinococcus canadensis (G7) genome: a key knowledge of parasitic platyhelminth human diseases

BACKGROUND

The parasite Echinococcus canadensis (G7) (phylum Platyhelminthes, class Cestoda) is one of the causative agents of echinococcosis. Echinococcosis is a worldwide chronic zoonosis affecting humans as well as domestic and wild mammals, which has been reported as a prioritized neglected disease by the World Health Organisation. No genomic data, comparative genomic analyses or efficient therapeutic and diagnostic tools are available for this severe disease. The information presented in this study will help to understand the peculiar biological characters and to design species-specific control tools.

RESULTS

We sequenced, assembled and annotated the 115-Mb genome of E. canadensis (G7). Comparative genomic analyses using whole genome data of three Echinococcus species not only confirmed the status of E. canadensis (G7) as a separate species but also demonstrated a high nucleotide sequences divergence in relation to E. granulosus (G1). The E. canadensis (G7) genome contains 11,449 genes with a core set of 881 orthologs shared among five cestode species. Comparative genomics revealed that there are more single nucleotide polymorphisms (SNPs) between E. canadensis (G7) and E. granulosus (G1) than between E. canadensis (G7) and E. multilocularis. This result was unexpected since E. canadensis (G7) and E. granulosus (G1) were considered to belong to the species complex E. granulosus sensu lato. We described SNPs in known drug targets and metabolism genes in the E. canadensis (G7) genome. Regarding gene regulation, we analysed three particular features: CpG island distribution along the three Echinococcus genomes, DNA methylation system and small RNA pathway. The results suggest the occurrence of yet unknown gene regulation mechanisms in Echinococcus.

CONCLUSIONS

This is the first work that addresses Echinococcus comparative genomics. The resources presented here will promote the study of mechanisms of parasite development as well as new tools for drug discovery. The availability of a high-quality genome assembly is critical for fully exploring the biology of a pathogenic organism. The E. canadensis (G7) genome presented in this study provides a unique opportunity to address the genetic diversity among the genus Echinococcus and its particular developmental features. At present, there is no unequivocal taxonomic classification of Echinococcus species; however, the genome-wide SNPs analysis performed here revealed the phylogenetic distance among these three Echinococcus species. Additional cestode genomes need to be sequenced to be able to resolve their phylogeny.

Comparing and Evaluating Metagenome Assembly Tools from a Microbiologist's Perspective - Not Only Size Matters!

With the constant improvement in cost-efficiency and quality of Next Generation Sequencing technologies, shotgun-sequencing approaches -such as metagenomics- have nowadays become the methods of choice for studying and classifying microorganisms from various habitats. The production of data has dramatically increased over the past years and processing and analysis steps are becoming more and more of a bottleneck. Limiting factors are partly the availability of computational resources, but mainly the bioinformatics expertise in establishing and applying appropriate processing and analysis pipelines. Fortunately, a large diversity of specialized software tools is nowadays available. Nevertheless, choosing the most appropriate methods for answering specific biological questions can be rather challenging, especially for non-bioinformaticians. In order to provide a comprehensive overview and guide for the microbiological scientific community, we assessed the most common and freely available metagenome assembly tools with respect to their output statistics, their sensitivity for low abundant community members and variability in resulting community profiles as well as their ease-of-use. In contrast to the highly anticipated "Critical Assessment of Metagenomic Interpretation" (CAMI) challenge, which uses general mock community-based assembler comparison we here tested assemblers on real Illumina metagenome sequencing data from natural communities of varying complexity sampled from forest soil and algal biofilms. Our observations clearly demonstrate that different assembly tools can prove optimal, depending on the sample type, available computational resources and, most importantly, the specific research goal. In addition, we present detailed descriptions of the underlying principles and pitfalls of publically available assembly tools from a microbiologist's perspective, and provide guidance regarding the user-friendliness, sensitivity and reliability of the resulting phylogenetic profiles.

Exploring viral infection using single-cell sequencing

Single-cell sequencing (SCS) has emerged as a valuable tool to study cellular heterogeneity in diverse fields, including virology. By studying the viral and cellular genome and/or transcriptome, the dynamics of viral infection can be investigated at single cell level. Most studies have explored the impact of cell-to-cell variation on the viral life cycle from the point of view of the virus, by analyzing viral sequences, and from the point of view of the cell, mainly by analyzing the cellular host transcriptome. In this review, we will focus on recent studies that use single-cell sequencing to explore viral diversity and cell variability in response to viral replication.

Redundans: an assembly pipeline for highly heterozygous genomes

Many genomes display high levels of heterozygosity (i.e. presence of different alleles at the same loci in homologous chromosomes), being those of hybrid organisms an extreme such case. The assembly of highly heterozygous genomes from short sequencing reads is a challenging task because it is difficult to accurately recover the different haplotypes. When confronted with highly heterozygous genomes, the standard assembly process tends to collapse homozygous regions and reports heterozygous regions in alternative contigs. The boundaries between homozygous and heterozygous regions result in multiple assembly paths that are hard to resolve, which leads to highly fragmented assemblies with a total size larger than expected. This, in turn, causes numerous problems in downstream analyses such as fragmented gene models, wrong gene copy number, or broken synteny. To circumvent these caveats we have developed a pipeline that specifically deals with the assembly of heterozygous genomes by introducing a step to recognise and selectively remove alternative heterozygous contigs. We tested our pipeline on simulated and naturally-occurring heterozygous genomes and compared its accuracy to other existing tools. Our method is freely available at https://github.com/Gabaldonlab/redundans.

Understanding Haemonchus contortus Better Through Genomics and Transcriptomics

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.