Iterative Correction of Reference Nucleotides (iCORN) using second generation sequencing technology

Minimizing Reference Bias with an Impute-First Approach

Pangenome indexes reduce reference bias in sequencing data analysis. However, bias can be reduced further by using a personalized reference, e.g. a diploid human reference constructed to match a donor individual's alleles. We present a novel impute-first alignment framework that combines elements of genotype imputation and pangenome alignment. It begins by genotyping the individual using only a subsample of the input reads. It next uses a reference panel and efficient imputation algorithm to impute a personalized diploid reference. Finally, it indexes the personalized reference and applies a read aligner, which could be a linear or graph aligner, to align the full read set to the personalized reference. This framework achieves higher variant-calling recall (99.54% vs. 99.37%), precision (99.36% vs. 99.18%), and F1 (99.45% vs. 99.28%) compared to a graph pangenome aligner. The personalized reference is also smaller and faster to query compared to a pangenome index, making it an overall advantageous choice for whole-genome DNA sequencing experiments.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

Localization of Epigenetic Markers in Leishmania Chromatin

Eukaryotes use histone variants and post-translation modifications (PTMs), as well as DNA base modifications, to regulate DNA replication/repair, chromosome condensation, and gene expression. Despite the unusual organization of their protein-coding genes into large polycistronic transcription units (PTUs), trypanosomatid parasites also employ a “histone code” to control these processes, but the details of this epigenetic code are poorly understood. Here, we present the results of experiments designed to elucidate the distribution of histone variants and PTMs over the chromatin landscape of Leishmania tarentolae. These experiments show that two histone variants (H2A.Z and H2B.V) and three histone H3 PTMs (H3K4me3, H3K16ac, and H3K76me3) are enriched at transcription start sites (TSSs); while a histone variant (H3.V) and the trypanosomatid-specific hyper-modified DNA base J are located at transcription termination sites (TTSs). Reduced nucleosome density was observed at all TTSs and TSSs for RNA genes transcribed by RNA polymerases I (RNAPI) or RNAPIII; as well as (to a lesser extent) at TSSs for the PTUs transcribed by RNAPII. Several PTMs (H3K4me3, H3K16ac H3K20me2 and H3K36me3) and base J were enriched at centromeres, while H3K50ac was specifically associated with the periphery of these centromeric sequences. These findings significantly expand our knowledge of the epigenetic markers associated with transcription, DNA replication and/or chromosome segregation in these early diverging eukaryotes and will hopefully lay the groundwork for future studies to elucidate how they control these fundamental processes.

Possible stochastic sex determination in Bursaphelenchus nematodes

Sex determination mechanisms evolve surprisingly rapidly, yet little is known in the large nematode phylum other than for Caenorhabditis elegans, which relies on chromosomal XX-XO sex determination and a dosage compensation mechanism. Here we analyze by sex-specific genome sequencing and genetic analysis sex determination in two fungal feeding/plant-parasitic Bursaphelenchus nematodes and find that their sex differentiation is more likely triggered by random, epigenetic regulation than by more well-known mechanisms of chromosomal or environmental sex determination. There is no detectable difference in male and female chromosomes, nor any linkage to sexual phenotype. Moreover, the protein sets of these nematodes lack genes involved in X chromosome dosage counting or compensation. By contrast, our genetic screen for sex differentiation mutants identifies a Bursaphelenchus ortholog of tra-1, the major output of the C. elegans sex determination cascade. Nematode sex determination pathways might have evolved by “bottom-up” accretion from the most downstream regulator, tra-1.

In most species, sex is determined by genetic or environmental factors. Here, the authors present evidence that sex determination in Bursaphelenchus nematodes is instead likely to be regulated by a random, epigenetic mechanism.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveals expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri. We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most “missing” orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Paenarthrobacter sp. GOM3 Is a Novel Marine Species With Monoaromatic Degradation Relevance

Paenarthrobacter sp. GOM3, which is a strain that represents a new species-specific context within the genus Paenarthrobacter, is clearly a branched member independent of any group described thus far. This strain was recovered from marine sediments in the Gulf of Mexico, and despite being isolated from a consortium capable of growing with phenanthrene as a sole carbon source, this strain could not grow successfully in the presence of this substrate alone. We hypothesized that the GOM3 strain could participate in the assimilation of intermediate metabolites for the degradation of aromatic compounds. To date, there are no experimental reports of Paenarthrobacter species that degrade polycyclic aromatic hydrocarbons (PAHs) or their intermediate metabolites. In this work, we report genomic and experimental evidence of metabolic benzoate, gentisate, and protocatechuate degradation by Paenarthrobacter sp. GOM3. Gentisate was the preferred substrate with the highest volumetric consumption rate, and genomic analysis revealed that this strain possesses multiple gene copies for the specific transport of gentisate. Furthermore, upon analyzing the GOM3 genome, we found five different dioxygenases involved in the activation of aromatic compounds, suggesting its potential for complete remediation of PAH-contaminated sites in combination with strains capable of assimilating the upper PAH degradation pathway. Additionally, this strain was characterized experimentally for its pathogenic potential and in silico for its antimicrobial resistance. An overview of the potential ecological role of this strain in the context of other members of this taxonomic clade is also reported.

Plasmodium vinckei genomes provide insights into the pan-genome and evolution of rodent malaria parasites

Background

Rodent malaria parasites (RMPs) serve as tractable tools to study malaria parasite biology and host-parasite-vector interactions. Among the four RMPs originally collected from wild thicket rats in sub-Saharan Central Africa and adapted to laboratory mice, Plasmodium vinckei is the most geographically widespread with isolates collected from five separate locations. However, there is a lack of extensive phenotype and genotype data associated with this species, thus hindering its use in experimental studies.

Results

We have generated a comprehensive genetic resource for P. vinckei comprising of five reference-quality genomes, one for each of its subspecies, blood-stage RNA sequencing data for five P. vinckei isolates, and genotypes and growth phenotypes for ten isolates. Additionally, we sequenced seven isolates of the RMP species Plasmodium chabaudi and Plasmodium yoelii, thus extending genotypic information for four additional subspecies enabling a re-evaluation of the genotypic diversity and evolutionary history of RMPs.

The five subspecies of P. vinckei have diverged widely from their common ancestor and have undergone large-scale genome rearrangements. Comparing P. vinckei genotypes reveals region-specific selection pressures particularly on genes involved in mosquito transmission. Using phylogenetic analyses, we show that RMP multigene families have evolved differently across the vinckei and berghei groups of RMPs and that family-specific expansions in P. chabaudi and P. vinckei occurred in the common vinckei group ancestor prior to speciation. The erythrocyte membrane antigen 1 and fam-c families in particular show considerable expansions among the lowland forest-dwelling P. vinckei parasites. The subspecies from the highland forests of Katanga, P. v. vinckei, has a uniquely smaller genome, a reduced multigene family repertoire and is also amenable to transfection making it an ideal parasite for reverse genetics. We also show that P. vinckei parasites are amenable to genetic crosses.

Conclusions

Plasmodium vinckei isolates display a large degree of phenotypic and genotypic diversity and could serve as a resource to study parasite virulence and immunogenicity. Inclusion of P. vinckei genomes provide new insights into the evolution of RMPs and their multigene families. Amenability to genetic crossing and transfection make them also suitable for classical and functional genetics to study Plasmodium biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00995-5.

Additional description and genome analyses of Caenorhabditis auriculariae representing the basal lineage of genus Caenorhabditis

Caenorhabditis auriculariae, which was morphologically described in 1999, was re-isolated from a Platydema mushroom-associated beetle. Based on the re-isolated materials, some morphological characteristics were re-examined and ascribed to the species. In addition, to clarify phylogenetic relationships with other Caenorhabditis species and biological features of the nematode, the whole genome was sequenced and assembled into 109.5 Mb with 16,279 predicted protein-coding genes. Molecular phylogenetic analyses based on ribosomal RNA and 269 single-copy genes revealed the species is closely related to C. sonorae and C. monodelphis placing them at the most basal clade of the genus. C. auriculariae has morphological characteristics clearly differed from those two species and harbours a number of species-specific gene families, indicating its usefulness as a new outgroup species for Caenorhabditis evolutionary studies. A comparison of carbohydrate-active enzyme (CAZy) repertoires in genomes, which we found useful to speculate about the lifestyle of Caenorhabditis nematodes, suggested that C. auriculariae likely has a life-cycle with tight-association with insects.

Integration of Genomic and Transcriptomic Data to Elucidate Molecular Processes in Babesia divergens

Emerging pathogens have developed ingenious life cycles to facilitate their growth and survival in the host organism. Detailed knowledge of the life cycle of these pathogens is increasingly necessary if we are to design new strategies to prevent infection and transmission. Multi-omics platforms provide useful data at different biological levels, and integration of these data into current approaches can facilitate holistic assessment of emerging pathogens. In this chapter, we bring together various methods and apply an integrative approach for analysis of genomic and transcriptomic data in Babesia divergens, an Apicomplexa emerging parasite that invades red blood cells and causes redwater fever in cattle and the most severe form of babesiosis in humans in Europe. The integrative methodology described herein can be helpful to identify genes active at specific points during life cycle of Apicomplexa parasites.

Strain-specific genome evolution in Trypanosoma cruzi, the agent of Chagas disease

The protozoan Trypanosoma cruzi almost invariably establishes life-long infections in humans and other mammals, despite the development of potent host immune responses that constrain parasite numbers. The consistent, decades-long persistence of T. cruzi in human hosts arises at least in part from the remarkable level of genetic diversity in multiple families of genes encoding the primary target antigens of anti-parasite immune responses. However, the highly repetitive nature of the genome-largely a result of these same extensive families of genes-have prevented a full understanding of the extent of gene diversity and its maintenance in T. cruzi. In this study, we have combined long-read sequencing and proximity ligation mapping to generate very high-quality assemblies of two T. cruzi strains representing the apparent ancestral lineages of the species. These assemblies reveal not only the full repertoire of the members of large gene families in the two strains, demonstrating extreme diversity within and between isolates, but also provide evidence of the processes that generate and maintain that diversity, including extensive gene amplification, dispersion of copies throughout the genome and diversification via recombination and in situ mutations. Gene amplification events also yield significant copy number variations in a substantial number of genes presumably not required for or involved in immune evasion, thus forming a second level of strain-dependent variation in this species. The extreme genome flexibility evident in T. cruzi also appears to create unique challenges with respect to preserving core genome functions and gene expression that sets this species apart from related kinetoplastids.

Population genomic evidence that human and animal infections in Africa come from the same populations of Dracunculus medinensis

BACKGROUND

Guinea worm-Dracunculus medinensis-was historically one of the major parasites of humans and has been known since antiquity. Now, Guinea worm is on the brink of eradication, as efforts to interrupt transmission have reduced the annual burden of disease from millions of infections per year in the 1980s to only 54 human cases reported globally in 2019. Despite the enormous success of eradication efforts to date, one complication has arisen. Over the last few years, hundreds of dogs have been found infected with this previously apparently anthroponotic parasite, almost all in Chad. Moreover, the relative numbers of infections in humans and dogs suggests that dogs are currently the principal reservoir on infection and key to maintaining transmission in that country.

PRINCIPAL FINDINGS

In an effort to shed light on this peculiar epidemiology of Guinea worm in Chad, we have sequenced and compared the genomes of worms from dog, human and other animal infections. Confirming previous work with other molecular markers, we show that all of these worms are D. medinensis, and that the same population of worms are causing both infections, can confirm the suspected transmission between host species and detect signs of a population bottleneck due to the eradication efforts. The diversity of worms in Chad appears to exclude the possibility that there were no, or very few, worms present in the country during a 10-year absence of reported cases.

CONCLUSIONS

This work reinforces the importance of adequate surveillance of both human and dog populations in the Guinea worm eradication campaign and suggests that control programs aiming to interrupt disease transmission should stay aware of the possible emergence of unusual epidemiology as pathogens approach elimination.

Comparative genomic analysis identifies X-factor (haemin)-independent Haemophilus haemolyticus: a formal re-classification of 'Haemophilus intermedius'

The heterogeneous and highly recombinogenic genus Haemophilus comprises several species, some of which are pathogenic to humans. All share an absolute requirement for blood-derived factors during growth. Certain species, such as the pathogen Haemophilus influenzae and the commensal Haemophilus haemolyticus, are thought to require both haemin (X-factor) and nicotinamide adenine dinucleotide (NAD, V-factor), whereas others, such as the informally classified 'Haemophilus intermedius subsp. intermedius', and Haemophilus parainfluenzae, only require V-factor. These differing growth requirements are commonly used for species differentiation, although a number of studies are now revealing issues with this approach. Here, we perform large-scale phylogenomics of 240 Haemophilus spp. genomes, including five 'H. intermedius' genomes generated in the current study, to reveal that strains of the 'H. intermedius' group are in fact haemin-independent H. haemolyticus (hiHh). Closer examination of these hiHh strains revealed that they encode an intact haemin biosynthesis pathway, unlike haemin-dependent H. haemolyticus and H. influenzae, which lack most haemin biosynthesis genes. Our results suggest that the common ancestor of modern-day H. haemolyticus and H. influenzae lost key haemin biosynthesis loci, likely as a consequence of specialized adaptation to otorhinolaryngeal and respiratory niches during their divergence from H. parainfluenzae. Genetic similarity analysis demonstrated that the haemin biosynthesis loci acquired in the hiHh lineage were likely laterally transferred from a H. parainfluenzae ancestor, and that this event probably occurred only once in hiHh. This study further challenges the validity of phenotypic methods for differentiating among Haemophilus species, and highlights the need for whole-genome sequencing for accurate characterization of species within this taxonomically challenging genus.

Genomic and transcriptomic evidence for descent from Plasmodium and loss of blood schizogony in Hepatocystis parasites from naturally infected red colobus monkeys

Hepatocystis is a genus of single-celled parasites infecting, amongst other hosts, monkeys, bats and squirrels. Although thought to have descended from malaria parasites (Plasmodium spp.), Hepatocystis spp. are thought not to undergo replication in the blood–the part of the Plasmodium life cycle which causes the symptoms of malaria. Furthermore, Hepatocystis is transmitted by biting midges, not mosquitoes. Comparative genomics of Hepatocystis and Plasmodium species therefore presents an opportunity to better understand some of the most important aspects of malaria parasite biology. We were able to generate a draft genome for Hepatocystis sp. using DNA sequencing reads from the blood of a naturally infected red colobus monkey. We provide robust phylogenetic support for Hepatocystis sp. as a sister group to Plasmodium parasites infecting rodents. We show transcriptomic support for a lack of replication in the blood and genomic support for a complete loss of a family of genes involved in red blood cell invasion. Our analyses highlight the rapid evolution of genes involved in parasite vector stages, revealing genes that may be critical for interactions between malaria parasites and mosquitoes.

Hepatocystis parasites are single-celled organisms, closely related to the Plasmodium species which cause malaria. But Hepatocystis are distinct–unlike Plasmodium they are thought not to replicate in the blood and cause little or no disease in their mammalian hosts. They are transmitted from one host to the next, not by mosquitoes, but by biting midges. In this study we generated a genome sequence for Hepatocystis–the first time this data has ever been produced and analysed for this species. We compared genome sequences of Hepatocystis and Plasmodium, confirming that Hepatocystis is descended from Plasmodium. We strengthened support for the absence of replication in the blood and, in line with this finding, discovered that genes involved in interaction with red blood cells have been lost in Hepatocystis. Our analyses revealed rapid evolution of genes which are active when the parasite is in the insect vector, highlighting those which might be important for understanding interaction between malaria parasites and mosquitoes. Hepatocystis has a fascinating evolutionary story and is a powerful comparator for understanding malaria parasite biology.

Population genomics of Vibrionaceae isolated from an endangered oasis reveals local adaptation after an environmental perturbation

Background

In bacteria, pan-genomes are the result of an evolutionary “tug of war” between selection and horizontal gene transfer (HGT). High rates of HGT increase the genetic pool and the effective population size (N_e), resulting in open pan-genomes. In contrast, selective pressures can lead to local adaptation by purging the variation introduced by HGT and mutation, resulting in closed pan-genomes and clonal lineages. In this study, we explored both hypotheses, elucidating the pan-genome of Vibrionaceae isolates after a perturbation event in the endangered oasis of Cuatro Ciénegas Basin (CCB), Mexico, and looking for signals of adaptation to the environments in their genomes.

Results

We obtained 42 genomes of Vibrionaceae distributed in six lineages, two of them did not showed any close reference strain in databases. Five of the lineages showed closed pan-genomes and were associated to either water or sediment environment; their high N_e estimates suggest that these lineages are not from a recent origin. The only clade with an open pan-genome was found in both environments and was formed by ten genetic groups with low N_e, suggesting a recent origin. The recombination and mutation estimators (r/m) ranged from 0.005 to 2.725, which are similar to oceanic Vibrionaceae estimations. However, we identified 367 gene families with signals of positive selection, most of them found in the core genome; suggesting that despite recombination, natural selection moves the Vibrionaceae CCB lineages to local adaptation, purging the genomes and keeping closed pan-genome patterns. Moreover, we identify 598 SNPs associated with an unstructured environment; some of the genes associated with these SNPs were related to sodium transport.

Conclusions

Different lines of evidence suggest that the sampled Vibrionaceae, are part of the rare biosphere usually living under famine conditions. Two of these lineages were reported for the first time. Most Vibrionaceae lineages of CCB are adapted to their micro-habitats rather than to the sampled environments. This pattern of adaptation is concordant with the association of closed pan-genomes and local adaptation.

Nearly Complete Genome Sequence of Brugia malayi Strain FR3

Lymphatic filariasis affects ∼120 million people and can result in elephantiasis and hydrocele. Here, we report the nearly complete genome sequence of the best-studied causative agent of lymphatic filariasis, Brugia malayi. The assembly contains four autosomes, an X chromosome, and only eight gaps but lacks a contiguous sequence for the known Y chromosome.

Lymphatic filariasis affects ∼120 million people and can result in elephantiasis and hydrocele. Here, we report the nearly complete genome sequence of the best-studied causative agent of lymphatic filariasis, Brugia malayi. The assembly contains four autosomes, an X chromosome, and only eight gaps but lacks a contiguous sequence for the known Y chromosome.

Characterization of a Novel Functional Trimeric Catechol 1,2-Dioxygenase From a Pseudomonas stutzeri Isolated From the Gulf of Mexico

Catechol 1,2 dioxygenases (C12DOs) have been studied for its ability to cleavage the benzene ring of catechol, the main intermediate in the degradation of aromatic compounds derived from aerobic degradation of hydrocarbons. Here we report the genome sequence of the marine bacterium Pseudomonas stutzeri GOM2, isolated from the southwestern Gulf of Mexico, and the biochemical characterization of its C12DO (PsC12DO). The catA gene, encoding PsC12DO of 312 amino acid residues, was cloned and expressed in Escherichia coli. Many C12DOs have been described as dimeric enzymes including those present in Pseudomonas species. The purified PsC12DO enzyme was found as an active trimer, with a molecular mass of 107 kDa. Increasing NaCl concentration in the enzyme reaction gradually reduced activity; in high salt concentrations (0.7 M NaCl) quaternary structural analysis determined that the enzyme changes to a dimeric arrangement and causes a 51% decrease in specific activity on catechol substrate. In comparison with other C12DOs, our enzyme showed a broad range of action for PsC12DO in solutions with pH values ranging from neutral to alkaline (70%). The enzyme is still active after incubation at 50°C for 30 min and in low temperatures to long term storage after 6 weeks at 4°C (61%). EDTA or Ca2+ inhibitors cause no drastic changes on residual activity; nevertheless, the activity of the enzyme was affected by metal ions Fe3+, Zn2+ and was completely inhibited by Hg2+. Under optimal conditions the k cat and K m values were 16.13 s-1 and 13.2 μM, respectively. To our knowledge, this is the first report describing the characterization of a marine C12DOs from P. stutzeri isolated from the Gulf of Mexico that is active in a trimeric state. We consider that our enzyme has important features to be used in environments in presence of EDTA, metals and salinity conditions.

Reduced Genome of the Gut Symbiotic Bacterium "Candidatus Benitsuchiphilus tojoi" Provides Insight Into Its Possible Roles in Ecology and Adaptation of the Host Insect

Diverse animals, including insects, harbor microbial symbionts within their gut, body cavity, or cells. The subsocial parastrachiid stinkbug Parastrachia japonensis is well-known for its peculiar ecological and behavioral traits, including its prolonged non-feeding diapause period and maternal care of eggs/nymphs in an underground nest. P. japonensis harbors a specific bacterial symbiont within the gut cavity extracellularly, which is vertically inherited through maternal excretion of symbiont-containing white mucus. Thus far, biological roles of the symbiont in the host lifecycle has been little understood. Here we sequenced the genome of the uncultivable gut symbiont “Candidatus Benitsuchiphilus tojoi.” The symbiont has an 804 kb circular chromosome encoding 606 proteins and a 14.5 kb plasmid encoding 13 proteins. Phylogenetic analysis indicated that the bacterium is closely related to other obligate insect symbionts belonging to the Gammaproteobacteria, including Buchnera of aphids and Blochmannia of ants, and the most closely related to Ishikawaella, an extracellular gut symbiont of plataspid stinkbugs. These data suggested that the symbiont genome has evolved like highly reduced gamma-proteobacterial symbiont genomes reported from a variety of insects. The presence of genes involved in biosynthesis pathways for amino acids, vitamins, and cofactors in the genome implicated the symbiont as a nutritional mutualist, supplementing essential nutrients to the host. Interestingly, the symbiont’s plasmid encoded genes for thiamine and carotenoid synthesis pathways, suggesting the possibility of additional functions of the symbiont for protecting the host against oxidative stress and DNA damage. Finally, possible involvement of the symbiont in uric acid metabolism during diapause is discussed.

Genetic Markers in S. Paratyphi C Reveal Primary Adaptation to Pigs

Salmonella enterica with the identical antigenic formula 6,7:c:1,5 can be differentiated biochemically and by disease syndrome. One grouping, Salmonella Paratyphi C, is currently considered a typhoidal serovar, responsible for enteric fever in humans. The human-restricted typhoidal serovars (S. Typhi and Paratyphi A, B and C) typically display high levels of genome degradation and are cited as an example of convergent evolution for host adaptation in humans. However, S. Paratyphi C presents a different clinical picture to S. Typhi/Paratyphi A, in a patient group with predisposition, raising the possibility that its natural history is different, and that infection is invasive salmonellosis rather than enteric fever. Using whole genome sequencing and metabolic pathway analysis, we compared the genomes of 17 S. Paratyphi C strains to other members of the 6,7:c:1,5 group and to two typhoidal serovars: S. Typhi and Paratyphi A. The genome degradation observed in S. Paratyphi C was much lower than S. Typhi/Paratyphi A, but similar to the other 6,7:c:1,5 strains. Genomic and metabolic comparisons revealed little to no overlap between S. Paratyphi C and the other typhoidal serovars, arguing against convergent evolution and instead providing evidence of a primary adaptation to pigs in accordance with the 6,7:c:1.5 strains.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage

BACKGROUND

In short-read DNA sequencing experiments, the read coverage is a key parameter to successfully assemble the reads and reconstruct the sequence of the input DNA. When coverage is very low, the original sequence reconstruction from the reads can be difficult because of the occurrence of uncovered gaps. Reference guided assembly can then improve these assemblies. However, when the available reference is phylogenetically distant from the sequencing reads, the mapping rate of the reads can be extremely low. Some recent improvements in read mapping approaches aim at modifying the reference according to the reads dynamically. Such approaches can significantly improve the alignment rate of the reads onto distant references but the processing of insertions and deletions remains challenging.

RESULTS

Here, we introduce a new algorithm to update the reference sequence according to previously aligned reads. Substitutions, insertions and deletions are performed in the reference sequence dynamically. We evaluate this approach to assemble a western-grey kangaroo mitochondrial amplicon. Our results show that more reads can be aligned and that this method produces assemblies of length comparable to the truth while limiting error rate when classic approaches fail to recover the correct length. Finally, we discuss how the core algorithm of this method could be improved and combined with other approaches to analyse larger genomic sequences.

CONCLUSIONS

We introduced an algorithm to perform dynamic alignment of reads on a distant reference. We showed that such approach can improve the reconstruction of an amplicon compared to classically used bioinformatic pipelines. Although not portable to genomic scale in the current form, we suggested several improvements to be investigated to make this method more flexible and allow dynamic alignment to be used for large genome assemblies.

Evolutionary analysis of the most polymorphic gene family in falciparum malaria

The var gene family of the human malaria parasite Plasmodium falciparum encode proteins that are crucial determinants of both pathogenesis and immune evasion and are highly polymorphic. Here we have assembled nearly complete var gene repertoires from 2398 field isolates and analysed a normalised set of 714 from across 12 countries. This therefore represents the first large scale attempt to catalogue the worldwide distribution of var gene sequences We confirm the extreme polymorphism of this gene family but also demonstrate an unexpected level of sequence sharing both within and between continents. We show that this is likely due to both the remnants of selective sweeps as well as a worrying degree of recent gene flow across continents with implications for the spread of drug resistance. We also address the evolution of the var repertoire with respect to the ancestral genes within the Laverania and show that diversity generated by recombination is concentrated in a number of hotspots. An analysis of the subdomain structure indicates that some existing definitions may need to be revised From the analysis of this data, we can now understand the way in which the family has evolved and how the diversity is continuously being generated. Finally, we demonstrate that because the genes are distributed across the genome, sequence sharing between genotypes acts as a useful population genetic marker.

Transcriptional and genomic parallels between the monoxenous parasite Herpetomonas muscarum and Leishmania

Trypanosomatid parasites are causative agents of important human and animal diseases such as sleeping sickness and leishmaniasis. Most trypanosomatids are transmitted to their mammalian hosts by insects, often belonging to Diptera (or true flies). These are called dixenous trypanosomatids since they infect two different hosts, in contrast to those that infect just insects (monoxenous). However, it is still unclear whether dixenous and monoxenous trypanosomatids interact similarly with their insect host, as fly-monoxenous trypanosomatid interaction systems are rarely reported and under-studied-despite being common in nature. Here we present the genome of monoxenous trypanosomatid Herpetomonas muscarum and discuss its transcriptome during in vitro culture and during infection of its natural insect host Drosophila melanogaster. The H. muscarum genome is broadly syntenic with that of human parasite Leishmania major. We also found strong similarities between the H. muscarum transcriptome during fruit fly infection, and those of Leishmania during sand fly infections. Overall this suggests Drosophila-Herpetomonas is a suitable model for less accessible insect-trypanosomatid host-parasite systems such as sand fly-Leishmania.

Genetics of Adaptation of the Ascomycetous Fungus Podospora anserina to Submerged Cultivation

Podospora anserina is a model ascomycetous fungus which shows pronounced phenotypic senescence when grown on solid medium but possesses unlimited lifespan under submerged cultivation. In order to study the genetic aspects of adaptation of P. anserina to submerged cultivation, we initiated a long-term evolution experiment. In the course of the first 4 years of the experiment, 125 single-nucleotide substitutions and 23 short indels were fixed in eight independently evolving populations. Six proteins that affect fungal growth and development evolved in more than one population; in particular, in the G-protein alpha subunit FadA, new alleles fixed in seven out of eight experimental populations, and these fixations affected just four amino acid sites, which is an unprecedented level of parallelism in experimental evolution. Parallel evolution at the level of genes and pathways, an excess of nonsense and missense substitutions, and an elevated conservation of proteins and their sites where the changes occurred suggest that many of the observed fixations were adaptive and driven by positive selection.

Is reliance on an inaccurate genome sequence sabotaging your experiments?

Advances in genomics have made whole genome studies increasingly feasible across the life sciences. However, new technologies and algorithmic advances do not guarantee flawless genomic sequences or annotation. Bias, errors, and artifacts can enter at any stage of the process from library preparation to annotation. When planning an experiment that utilizes a genome sequence as the basis for the design, there are a few basic checks that, if performed, may better inform the experimental design and ideally help avoid a failed experiment or inconclusive result.

Comparative and functional genomics of the protozoan parasite Babesia divergens highlighting the invasion and egress processes

Babesiosis is considered an emerging disease because its incidence has significantly increased in the last 30 years, providing evidence of the expanding range of this rare but potentially life-threatening zoonotic disease. Babesia divergens is a causative agent of babesiosis in humans and cattle in Europe. The recently sequenced genome of B. divergens revealed over 3,741 protein coding-genes and the 10.7-Mb high-quality draft become the first reference tool to study the genome structure of B. divergens. Now, by exploiting this sequence data and using new computational tools and assembly strategies, we have significantly improved the quality of the B. divergens genome. The new assembly shows better continuity and has a higher correspondence to B. bovis chromosomes. Moreover, we present a differential expression analysis using RNA sequencing of the two different stages of the asexual lifecycle of B. divergens: the free merozoite capable of invading erythrocytes and the intraerythrocytic parasite stage that remains within the erythrocyte until egress. Comparison of mRNA levels of both stages identified 1,441 differentially expressed genes. From these, around half were upregulated and the other half downregulated in the intraerythrocytic stage. Orthogonal validation by real-time quantitative reverse transcription PCR confirmed the differential expression. A moderately increased expression level of genes, putatively involved in the invasion and egress processes, were revealed in the intraerythrocytic stage compared with the free merozoite. On the basis of these results and in the absence of molecular models of invasion and egress for B. divergens, we have proposed the identified genes as putative molecular players in the invasion and egress processes. Our results contribute to an understanding of key parasitic strategies and pathogenesis and could be a valuable genomic resource to exploit for the design of diagnostic methods, drugs and vaccines to improve the control of babesiosis.

Babesiosis has long been recognized as an economically important disease of cattle, but only in the last 40 years has Babesia been recognized as an important pathogen in humans. Babesiosis in humans is caused by one of several species (B. microti, B. divergens, B. duncani and B. venatorum). The complete Babesia lifecycle requires two hosts, the ixodid ticks and a vertebrate host. It is the parasite's ability to first recognize and then invade host erythrocytes that is central to the pathogenesis of babesiosis. Once inside the cell, the parasite begins a cycle of maturation and growth, resulting in merozoites that egress from the red blood cells (RBCs) and seek new, uninfected RBCs to invade, perpetuating the infection. To better understand this asexual lifecycle, the authors focused on the parasite genome and transcriptome of the asexual erythrocytic forms of B. divergens. Through this functional and comparative genomic approach, the authors have identified genes putatively involved in invasion, gliding motility, moving junction formation and egress, providing new insights into the molecular mechanisms of these processes necessary for B. divergens to survive and propagate during its life cycle.

Dramatic changes in gene expression in different forms of Crithidia fasciculata reveal potential mechanisms for insect-specific adhesion in kinetoplastid parasites

Kinetoplastids are a group of parasites that includes several medically-important species. These human-infective species are transmitted by insect vectors in which the parasites undergo specific developmental transformations. For each species, this includes a stage in which parasites adhere to insect tissue via a hemidesmosome-like structure. Although this structure has been described morphologically, it has never been molecularly characterized. We are using Crithidia fasciculata, an insect parasite that produces large numbers of adherent parasites inside its mosquito host, as a model kinetoplastid to investigate both the mechanism of adherence and the signals required for differentiation to an adherent form. An advantage of C. fasciculata is that adherent parasites can be generated both in vitro, allowing a direct comparison to cultured swimming forms, as well as in vivo within the mosquito. Using RNAseq, we identify genes associated with adherence in C. fasciculata. As almost all of these genes have orthologs in other kinetoplastid species, our findings may reveal shared mechanisms of adherence, allowing investigation of a crucial step in parasite development and disease transmission. In addition, dual-RNAseq allowed us to explore the interaction between the parasites and the mosquito. Although the infection is well-tolerated, anti-microbial peptides and other components of the mosquito innate immune system are upregulated. Our findings indicate that C. fasciculata is a powerful model system for probing kinetoplastid-insect interactions.

Kinetoplastids are single-celled parasites that cause devastating human diseases worldwide. Although this group includes many species that infect a variety of hosts, they have a great deal of shared biology. One relatively unexplored aspect of the kinetoplastid life cycle is their ability to adhere to insect tissue. For pathogenic species, adherence is critical for transmission by insect vectors. We have used an insect parasite called Crithidia fasciculata as a model kinetoplastid to reveal shared mechanisms of insect adherence. We have compared gene expression profiles of motile, non-adherent C. fasciculata to those of C. fasciculata adhered to non-living substrates and those attached to the hindgut of mosquitoes. Through this analysis, we have identified a large number of candidate proteins that may mediate adhesion in these and related parasites. In addition, our findings suggest that the mosquito immune system is responding to the presence of parasites in the gut. These results establish a new, robust system to explore the interaction between kinetoplastids and their insect hosts.

BAUM: improving genome assembly by adaptive unique mapping and local overlap-layout-consensus approach

Motivation

It is highly desirable to assemble genomes of high continuity and consistency at low cost. The current bottleneck of draft genome continuity using the second generation sequencing (SGS) reads is primarily caused by uncertainty among repetitive sequences. Even though the single-molecule real-time sequencing technology is very promising to overcome the uncertainty issue, its relatively high cost and error rate add burden on budget or computation. Many long-read assemblers take the overlap-layout-consensus (OLC) paradigm, which is less sensitive to sequencing errors, heterozygosity and variability of coverage. However, current assemblers of SGS data do not sufficiently take advantage of the OLC approach.

Results

Aiming at minimizing uncertainty, the proposed method BAUM, breaks the whole genome into regions by adaptive unique mapping; then the local OLC is used to assemble each region in parallel. BAUM can (i) perform reference-assisted assembly based on the genome of a close species (ii) or improve the results of existing assemblies that are obtained based on short or long sequencing reads. The tests on two eukaryote genomes, a wild rice Oryza longistaminata and a parrot Melopsittacus undulatus, show that BAUM achieved substantial improvement on genome size and continuity. Besides, BAUM reconstructed a considerable amount of repetitive regions that failed to be assembled by existing short read assemblers. We also propose statistical approaches to control the uncertainty in different steps of BAUM.

Availability and implementation

http://www.zhanyuwang.xin/wordpress/index.php/2017/07/21/baum.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mitochondrial genomics of human pathogenic parasite Leishmania (Viannia) panamensis

Background

The human parasite Leishmania (V.) panamensis is one of the pathogenic species responsible for cutaneous leishmaniasis in Central and South America. Despite its importance in molecular parasitology, its mitochondrial genome, divided into minicircles and maxicircles, haven’t been described so far.

Methods

Using NGS-based sequencing (454 and ILLUMINA), and combining de novo genome assembly and mapping strategies, we report the maxicircle kDNA annotated genome of L. (V.) panamensis, the first reference of this molecule for the subgenus Viannia. A comparative genomics approach is performed against other Leishmania and Trypanosoma species.

Results

The results show synteny of mitochondrial genes of L. (V.) panamensis with other kinetoplastids. It was also possible to identify nucleotide variants within the coding regions of the maxicircle, shared among some of them and others specific to each strain. Furthermore, we compared the minicircles kDNA sequences of two strains and the results show that the conserved and divergent regions of the minicircles exhibit strain-specific associations.

A Genetically Tractable, Natural Mouse Model of Cryptosporidiosis Offers Insights into Host Protective Immunity

Cryptosporidium is a leading cause of diarrheal disease and an important contributor to early childhood mortality, malnutrition, and growth faltering. Older children in high endemicity regions appear resistant to infection, while previously unexposed adults remain susceptible. Experimental studies in humans and animals support the development of disease resistance, but we do not understand the mechanisms that underlie protective immunity to Cryptosporidium. Here, we derive an in vivo model of Cryptosporidium infection in immunocompetent C57BL/6 mice by isolating parasites from naturally infected wild mice. Similar to human cryptosporidiosis, this infection causes intestinal pathology, and interferon-γ controls early infection while T cells are critical for clearance. Importantly, mice that controlled a live infection were resistant to secondary challenge and vaccination with attenuated parasites provided protection equal to live infection. Both parasite and host are genetically tractable and this in vivo model will facilitate mechanistic investigation and rational vaccine design.

•

We isolated and sequenced Cryptosporidium tyzzeri, a natural mouse pathogen

•

C. tyzzeri can be genetically manipulated using CRISPR-driven homologous repair

•

C. tyzzeri models human cryptosporidiosis with T cell- and IFN-γ-dependent resolution

•

Mice develop protective immunity following both live infection and vaccination

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

Genomic and Functional Analysis of Emerging Virulent and Multidrug-Resistant Escherichia coli Lineage Sequence Type 648

The pathogenic extended-spectrum-beta-lactamase (ESBL)-producing Escherichia coli lineage ST648 is increasingly reported from multiple origins. Our study of a large and global ST648 collection from various hosts (87 whole-genome sequences) combining core and accessory genomics with functional analyses and in vivo experiments suggests that ST648 is a nascent and generalist lineage, lacking clear phylogeographic and host association signals. By including large numbers of ST131 (n = 107) and ST10 (n = 96) strains for comparative genomics and phenotypic analysis, we demonstrate that the combination of multidrug resistance and high-level virulence are the hallmarks of ST648, similar to international high-risk clonal lineage ST131. Specifically, our in silico, in vitro, and in vivo results demonstrate that ST648 is well equipped with biofilm-associated features, while ST131 shows sophisticated signatures indicative of adaption to urinary tract infection, potentially conveying individual ecological niche adaptation. In addition, we used a recently developed NFDS (negative frequency-dependent selection) population model suggesting that ST648 will increase significantly in frequency as a cause of bacteremia within the next few years. Also, ESBL plasmids impacting biofilm formation aided in shaping and maintaining ST648 strains to successfully emerge worldwide across different ecologies. Our study contributes to understanding what factors drive the evolution and spread of emerging international high-risk clonal lineages.

Stage-specific transcriptome of Bursaphelenchus xylophilus reveals temporal regulation of effector genes and roles of the dauer-like stages in the lifecycle

The pine wood nematode Bursaphelenchus xylophilus is the causal agent of pine wilt disease, one of the most devastating forest diseases in East Asian and West European countries. The lifecycle of B. xylophilus includes four propagative larval stages and gonochoristic adults which are involved in the pathogenicity, and two stages of dispersal larvae involved in the spread of the disease. To elucidate the ecological roles of each developmental stage in the pathogenic life cycle, we performed a comprehensive transcriptome analysis using RNA-seq generated from all developmental stages of B. xylophilus and compared transcriptomes between stages. We found more than 9000 genes are differentially expressed in at least one stage of the life cycle including genes involved in general nematode biology such as reproduction and moulting but also effector genes likely to be involved in parasitism. The dispersal-stage transcriptome revealed its analogy to C. elegans dauer and the distinct roles of the two larval stages from each other regarding survival and transmission. This study provides important insights and resources to understand B. xylophilus parasitic biology.

lptG contributes to changes in membrane permeability and the emergence of multidrug hypersusceptibility in a cystic fibrosis isolate of Pseudomonas aeruginosa

PURPOSE

In the lungs of cystic fibrosis patients, Pseudomonas aeruginosa is exposed to a myriad of antibiotics leading to alterations in antibiotic susceptibility. This study identifies mutations resulting in hypersusceptibility in isogenic mutants of a P. aeruginosa clinical isolate, PA34.

METHODS

PA34 was exposed to subinhibitory concentrations of doripenem or meropenem during growth to mid-log phase. Antibiotic susceptibility of surviving colonies was determined by agar dilution. Two carbapenem-resistant colonies hypersusceptible to non-carbapenem antibiotics were selected for further analysis. Antibiotic resistance gene expression was evaluated by RT-rtPCR and OprD production by SDS-PAGE. PA34 and isogenic mutants were evaluated with whole genome sequencing. Sequence variants were confirmed by Sanger sequencing, and cognate genes in eight carbapenem-resistant clinical isolates hypersusceptible to non-carbapenem antibiotics were sequenced. Lipopolysaccharide preparations of PA34 and hypersusceptible mutants were evaluated with ProQ-Emerald stain.

RESULTS

Isogenic mutants showed 4- to 8-fold MIC increase for imipenem, meropenem, and doripenem. However, they were hypersusceptible (≥4-fold MIC decrease) to aminoglycosides, fluoroquinolones, and non-carbapenem β-lactams. Expression of ampC or mex-opr efflux pumps was unchanged, but OprD production was decreased. Mutations causing Q86H AlgU and G77C LptG amino acid substitutions and nonsense mutations within OprD were observed in both mutants. Lipopolysaccharide modifications were observed between isogenic mutants and PA34. Non-synonymous mutations in LptF or LptG were observed in 6/8 hypersusceptible clinical isolates resistant to carbapenem antibiotics.

CONCLUSION

Evaluation of hypersusceptible mutants identified the association between lptG and a hypersusceptible phenotype. Modifications in lipopolysaccharide profiles suggests LptG modification interferes with lipopolysaccharide transport and contributes to hypersusceptibility.

Draft Genome Assemblies of Clinical Isolates of Klebsiella pneumoniae V9011662 and Enterobacter hormaechei Entb306

Enterobacter hormaechei and Klebsiella pneumoniae are pathogenic Enterobacteriaceae that have been associated with the spread of antibiotic resistance. Here, we report draft genome assemblies of an Enterobacter hormaechei clinical isolate and a multidrug-resistant clinical isolate of Klebsiella pneumoniae.

Enterobacter hormaechei and Klebsiella pneumoniae are pathogenic Enterobacteriaceae that have been associated with the spread of antibiotic resistance. Here, we report draft genome assemblies of an Enterobacter hormaechei clinical isolate and a multidrug-resistant clinical isolate of Klebsiella pneumoniae.

Integrated transcriptomic and proteomic analysis of pathogenic mycobacteria and their esx-1 mutants reveal secretion-dependent regulation of ESX-1 substrates and WhiB6 as a transcriptional regulator

The mycobacterial type VII secretion system ESX-1 is responsible for the secretion of a number of proteins that play important roles during host infection. The regulation of the expression of secreted proteins is often essential to establish successful infection. Using transcriptome sequencing, we found that the abrogation of ESX-1 function in Mycobacterium marinum leads to a pronounced increase in gene expression levels of the espA operon during the infection of macrophages. In addition, the disruption of ESX-1-mediated protein secretion also leads to a specific down-regulation of the ESX-1 substrates, but not of the structural components of this system, during growth in culture medium. This effect is observed in both M. marinum and M. tuberculosis. We established that down-regulation of ESX-1 substrates is the result of a regulatory process that is influenced by the putative transcriptional regulator whib6, which is located adjacent to the esx-1 locus. In addition, the overexpression of the ESX-1-associated PE35/PPE68 protein pair resulted in a significantly increased secretion of the ESX-1 substrate EsxA, demonstrating a functional link between these proteins. Taken together, these data show that WhiB6 is required for the secretion-dependent regulation of ESX-1 substrates and that ESX-1 substrates are regulated independently from the structural components, both during infection and as a result of active secretion.

VirGenA: a reference-based assembler for variable viral genomes

Characterization of the within-host genetic diversity of viral pathogens is required for selection of effective treatment of some important viral infections, e.g. HIV, HBV and HCV. Despite the technical ability of detection, there are conflicting data regarding the clinical significance of low-frequency variants, partially because of the difficulty of their distinguishing from experimental artifacts. The issue of cross-contamination is relevant for all highly sensitive techniques, including deep sequencing: even trace contamination leads to a significant increase of false positives in identified SNVs. Determination of infections by multiple genotypes of some viruses, the incidence of which can be considerable, especially in risk groups, is also clinically significant in some cases. We developed a new viral reference-guided assembler, VirGenA, that can separate mixtures of strains of different intraspecies genetic groups (genotypes, subtypes, clades, etc.) and assemble a separate consensus sequence for each group in a mixture. It produced long assemblies for mixture components of extremely low frequencies (<1%) allowing detection of cross-contamination of samples by divergent genotypes. We tested VirGenA on both clinical and simulated data. On both types of data, VirGenA shows better or similar results than the existing de novo assemblers. Cross-platform implementation (including source code) is freely available at https://github.com/gFedonin/VirGenA/releases.

Comparative genomics of the major parasitic worms

Parasitic nematodes (roundworms) and platyhelminths (flatworms) cause debilitating chronic infections of humans and animals, decimate crop production and are a major impediment to socioeconomic development. Here we report a broad comparative study of 81 genomes of parasitic and non-parasitic worms. We have identified gene family births and hundreds of expanded gene families at key nodes in the phylogeny that are relevant to parasitism. Examples include gene families that modulate host immune responses, enable parasite migration though host tissues or allow the parasite to feed. We reveal extensive lineage-specific differences in core metabolism and protein families historically targeted for drug development. From an in silico screen, we have identified and prioritized new potential drug targets and compounds for testing. This comparative genomics resource provides a much-needed boost for the research community to understand and combat parasitic worms.

Trypanosoma cruzi Genome Assemblies: Challenges and Milestones of Assembling a Highly Repetitive and Complex Genome

Trypanosoma cruzi present one of the most complex parasite genomes sequenced to date. Among its features are 600-kb-long repetitive multigene families' clusters, hybrid strains, and aneuploidies, which hampered genome assembly completeness and contiguity. Several approaches, such as Sanger sequencing in 2005, next-generation sequencing in 2011 and third-generation sequencing in 2018, were used to improve draft assemblies of different strains of this parasite. Hence, the study of T. cruzi genome assemblies' history is an excellent way to describe the evolution of genome sequencing methodologies and compare their efficiency and limitations to assembly complex genomes. In this book chapter, we summarize the principal findings and methodologies of T. cruzi genome assembly projects to date, highlighting the improvements and limitations of each approach.

A Guide to Next Generation Sequence Analysis of Leishmania Genomes

Next generation sequencing (NGS) technology transformed Leishmania genome studies and became an indispensable tool for Leishmania researchers. Recent Leishmania genomics analyses facilitated the discovery of various genetic diversities including single nucleotide polymorphisms (SNPs), copy number variations (CNVs), somy variations, and structural variations in detail and provided valuable insights into the complexity of the genome and gene regulation. Many aspects of Leishmania NGS analyses are similar to those of related pathogens like trypanosomes. However, the analyses of Leishmania genomes face a unique challenge because of the presence of frequent aneuploidy. This makes characterization and interpretation of read depth and somy a key part of Leishmania NGS analyses because read depth affects the accuracy of detection of all genetic variations. However, there are no general guidelines on how to explore and interpret the impact of aneuploidy, and this has made it difficult for biologists and bioinformaticians, especially for beginners, to perform their own analyses and interpret results across different analyses. In this guide we discuss a wide range of topics essential for Leishmania NGS analyses, ranging from how to set up a computational environment for genome analyses, to how to characterize genetic variations among Leishmania samples, and we will particularly focus on chromosomal copy number variation and its impact on genome analyses.

RGAAT: A Reference-based Genome Assembly and Annotation Tool for New Genomes and Upgrade of Known Genomes

The rapid development of high-throughput sequencing technologies has led to a dramatic decrease in the money and time required for de novo genome sequencing or genome resequencing projects, with new genome sequences constantly released every week. Among such projects, the plethora of updated genome assemblies induces the requirement of version-dependent annotation files and other compatible public dataset for downstream analysis. To handle these tasks in an efficient manner, we developed the reference-based genome assembly and annotation tool (RGAAT), a flexible toolkit for resequencing-based consensus building and annotation update. RGAAT can detect sequence variants with comparable precision, specificity, and sensitivity to GATK and with higher precision and specificity than Freebayes and SAMtools on four DNA-seq datasets tested in this study. RGAAT can also identify sequence variants based on cross-cultivar or cross-version genomic alignments. Unlike GATK and SAMtools/BCFtools, RGAAT builds the consensus sequence by taking into account the true allele frequency. Finally, RGAAT generates a coordinate conversion file between the reference and query genomes using sequence variants and supports annotation file transfer. Compared to the rapid annotation transfer tool (RATT), RGAAT displays better performance characteristics for annotation transfer between different genome assemblies, strains, and species. In addition, RGAAT can be used for genome modification, genome comparison, and coordinate conversion. RGAAT is available at https://sourceforge.net/projects/rgaat/ and https://github.com/wushyer/RGAAT_v2 at no cost.

FORGe: prioritizing variants for graph genomes

There is growing interest in using genetic variants to augment the reference genome into a graph genome, with alternative sequences, to improve read alignment accuracy and reduce allelic bias. While adding a variant has the positive effect of removing an undesirable alignment score penalty, it also increases both the ambiguity of the reference genome and the cost of storing and querying the genome index. We introduce methods and a software tool called FORGe for modeling these effects and prioritizing variants accordingly. We show that FORGe enables a range of advantageous and measurable trade-offs between accuracy and computational overhead.

Genomic and transcriptomic comparisons of closely related malaria parasites differing in virulence and sequestration pattern

Background: Malaria parasite species differ greatly in the harm they do to humans. While P. falciparum kills hundreds of thousands per year, P. vivax kills much less often and P. malariae is relatively benign. Strains of the rodent malaria parasite Plasmodium chabaudi show phenotypic variation in virulence during infections of laboratory mice. This make it an excellent species to study genes which may be responsible for this trait. By understanding the mechanisms which underlie differences in virulence we can learn how parasites adapt to their hosts and how we might prevent disease. Methods: Here we present a complete reference genome sequence for a more virulent P. chabaudi strain, PcCB, and perform a detailed comparison with the genome of the less virulent PcAS strain. Results: We found the greatest variation in the subtelomeric regions, in particular amongst the sequences of the pir gene family, which has been associated with virulence and establishment of chronic infection. Despite substantial variation at the sequence level, the repertoire of these genes has been largely maintained, highlighting the requirement for functional conservation as well as diversification in host-parasite interactions. However, a subset of pir genes, previously associated with increased virulence, were more highly expressed in PcCB, suggesting a role for this gene family in virulence differences between strains. We found that core genes involved in red blood cell invasion have been under positive selection and that the more virulent strain has a greater preference for reticulocytes, which has elsewhere been associated with increased virulence. Conclusions: These results provide the basis for a mechanistic understanding of the phenotypic differences between Plasmodium chabaudi strains, which might ultimately be translated into a better understanding of malaria parasites affecting humans.

Genomic and transcriptomic comparisons of closely related malaria parasites differing in virulence and sequestration pattern

Background: Malaria parasite species differ greatly in the harm they do to humans. While P. falciparum kills hundreds of thousands per year, P. vivax kills much less often and P. malariae is relatively benign. Strains of the rodent malaria parasite Plasmodium chabaudi show phenotypic variation in virulence during infections of laboratory mice. This make it an excellent species to study genes which may be responsible for this trait. By understanding the mechanisms which underlie differences in virulence we can learn how parasites adapt to their hosts and how we might prevent disease.

Methods: Here we present a complete reference genome sequence for a more virulent P. chabaudi strain, PcCB, and perform a detailed comparison with the genome of the less virulent PcAS strain.

Results: We found the greatest variation in the subtelomeric regions, in particular amongst the sequences of the pir gene family, which has been associated with virulence and establishment of chronic infection. Despite substantial variation at the sequence level, the repertoire of these genes has been largely maintained, highlighting the requirement for functional conservation as well as diversification in host-parasite interactions. However, a subset of pir genes, previously associated with increased virulence, were more highly expressed in PcCB, suggesting a role for this gene family in virulence differences between strains. We found that core genes involved in red blood cell invasion have been under positive selection and that the more virulent strain has a greater preference for reticulocytes, which has elsewhere been associated with increased virulence.

Conclusions: These results provide the basis for a mechanistic understanding of the phenotypic differences between Plasmodium chabaudi strains, which might ultimately be translated into a better understanding of malaria parasites affecting humans.

Genomic Analysis of Colombian Leishmania panamensis strains with different level of virulence

The establishment of Leishmania infection in mammalian hosts and the subsequent manifestation of clinical symptoms require internalization into macrophages, immune evasion and parasite survival and replication. Although many of the genes involved in these processes have been described, the genetic and genomic variability associated to differences in virulence is largely unknown. Here we present the genomic variation of four Leishmania (Viannia) panamensis strains exhibiting different levels of virulence in BALB/c mice and its application to predict novel genes related to virulence. De novo DNA sequencing and assembly of the most virulent strain allowed comparative genomics analysis with sequenced L. (Viannia) panamensis and L. (Viannia) braziliensis strains, and showed important variations at intra and interspecific levels. Moreover, the mutation detection and a CNV search revealed both base and structural genomic variation within the species. Interestingly, we found differences in the copy number and protein diversity of some genes previously related to virulence. Several machine-learning approaches were applied to combine previous knowledge with features derived from genomic variation and predict a curated set of 66 novel genes related to virulence. These genes can be prioritized for validation experiments and could potentially become promising drug and immune targets for the development of novel prophylactic and therapeutic interventions.

Genome organization and DNA accessibility control antigenic variation in trypanosomes

Many evolutionarily distant pathogenic organisms have evolved similar survival strategies to evade the immune responses of their hosts. These include antigenic variation, through which an infecting organism prevents clearance by periodically altering the identity of proteins that are visible to the immune system of the host1. Antigenic variation requires large reservoirs of immunologically diverse antigen genes, which are often generated through homologous recombination, as well as mechanisms to ensure the expression of one or very few antigens at any given time. Both homologous recombination and gene expression are affected by three-dimensional genome architecture and local DNA accessibility2,3. Factors that link three-dimensional genome architecture, local chromatin conformation and antigenic variation have, to our knowledge, not yet been identified in any organism. One of the major obstacles to studying the role of genome architecture in antigenic variation has been the highly repetitive nature and heterozygosity of antigen-gene arrays, which has precluded complete genome assembly in many pathogens. Here we report the de novo haplotype-specific assembly and scaffolding of the long antigen-gene arrays of the model protozoan parasite Trypanosoma brucei, using long-read sequencing technology and conserved features of chromosome folding4. Genome-wide chromosome conformation capture (Hi-C) reveals a distinct partitioning of the genome, with antigen-encoding subtelomeric regions that are folded into distinct, highly compact compartments. In addition, we performed a range of analyses-Hi-C, fluorescence in situ hybridization, assays for transposase-accessible chromatin using sequencing and single-cell RNA sequencing-that showed that deletion of the histone variants H3.V and H4.V increases antigen-gene clustering, DNA accessibility across sites of antigen expression and switching of the expressed antigen isoform, via homologous recombination. Our analyses identify histone variants as a molecular link between global genome architecture, local chromatin conformation and antigenic variation.

Whole genome sequencing and microsatellite analysis of the Plasmodium falciparum E5 NF54 strain show that the var, rifin and stevor gene families follow Mendelian inheritance

Background

Plasmodium falciparum exhibits a high degree of inter-isolate genetic diversity in its variant surface antigen (VSA) families: P. falciparum erythrocyte membrane protein 1, repetitive interspersed family (RIFIN) and subtelomeric variable open reading frame (STEVOR). The role of recombination for the generation of this diversity is a subject of ongoing research. Here the genome of E5, a sibling of the 3D7 genome strain is presented. Short and long read whole genome sequencing (WGS) techniques (Ilumina, Pacific Bioscience) and a set of 84 microsatellites (MS) were employed to characterize the 3D7 and non-3D7 parts of the E5 genome. This is the first time that VSA genes in sibling parasites were analysed with long read sequencing technology.

Results

Of the 5733 E5 genes only 278 genes, mostly var and rifin/stevor genes, had no orthologues in the 3D7 genome. WGS and MS analysis revealed that chromosomal crossovers occurred at a rate of 0–3 per chromosome. var, stevor and rifin genes were inherited within the respective non-3D7 or 3D7 chromosomal context. 54 of the 84 MS PCR fragments correctly identified the respective MS as 3D7- or non-3D7 and this correlated with var and rifin/stevor gene inheritance in the adjacent chromosomal regions. E5 had 61 var and 189 rifin/stevor genes. One large non-chromosomal recombination event resulted in a new var gene on chromosome 14. The remainder of the E5 3D7-type subtelomeric and central regions were identical to 3D7.

Conclusions

The data show that the rifin/stevor and var gene families represent the most diverse compartments of the P. falciparum genome but that the majority of var genes are inherited without alterations within their respective parental chromosomal context. Furthermore, MS genotyping with 54 MS can successfully distinguish between two sibling progeny of a natural P. falciparum cross and thus can be used to investigate identity by descent in field isolates.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2503-2) contains supplementary material, which is available to authorized users.

RNAEditor: easy detection of RNA editing events and the introduction of editing islands

RNA editing of adenosine residues to inosine ('A-to-I editing') is the most common RNA modification event detectible with RNA sequencing (RNA-seq). While not directly detectable, inosine is read by next-generation sequencers as guanine. Therefore, mapping RNA-seq reads to their corresponding reference genome can detect potential editing events by identifying 'A-to-G' conversions. However, one must exercise caution when searching for editing sites, as A-to-G conversions also arise from sequencing errors as well as mutations. To address these complexities, several algorithms and software products have been developed to accurately identify editing events. Here, we survey currently available methods to analyze RNA editing events and introduce a new easy-to-use bioinformatics tool 'RNAEditor' for the detection of RNA editing events. During the development of RNAEditor, we noticed editing often happened in clusters, which we named 'editing islands'. We developed a clustering algorithm to find editing islands and included it in RNAEditor. RNAEditor is freely available at http://rnaeditor.uni-frankfurt.de. We anticipate that RNAEditor will provide biologists with an easy-to-use tool for studying RNA editing events and the newly defined editing islands.

Biology and genome of a newly discovered sibling species of Caenorhabditis elegans

A 'sibling' species of the model organism Caenorhabditis elegans has long been sought for use in comparative analyses that would enable deep evolutionary interpretations of biological phenomena. Here, we describe the first sibling species of C. elegans, C. inopinata n. sp., isolated from fig syconia in Okinawa, Japan. We investigate the morphology, developmental processes and behaviour of C. inopinata, which differ significantly from those of C. elegans. The 123-Mb C. inopinata genome was sequenced and assembled into six nuclear chromosomes, allowing delineation of Caenorhabditis genome evolution and revealing unique characteristics, such as highly expanded transposable elements that might have contributed to the genome evolution of C. inopinata. In addition, C. inopinata exhibits massive gene losses in chemoreceptor gene families, which could be correlated with its limited habitat area. We have developed genetic and molecular techniques for C. inopinata; thus C. inopinata provides an exciting new platform for comparative evolutionary studies.

HECIL: A Hybrid Error Correction Algorithm for Long Reads with Iterative Learning

Second-generation DNA sequencing techniques generate short reads that can result in fragmented genome assemblies. Third-generation sequencing platforms mitigate this limitation by producing longer reads that span across complex and repetitive regions. However, the usefulness of such long reads is limited because of high sequencing error rates. To exploit the full potential of these longer reads, it is imperative to correct the underlying errors. We propose HECIL-Hybrid Error Correction with Iterative Learning-a hybrid error correction framework that determines a correction policy for erroneous long reads, based on optimal combinations of decision weights obtained from short read alignments. We demonstrate that HECIL outperforms state-of-the-art error correction algorithms for an overwhelming majority of evaluation metrics on diverse, real-world data sets including E. coli, S. cerevisiae, and the malaria vector mosquito A. funestus. Additionally, we provide an optional avenue of improving the performance of HECIL's core algorithm by introducing an iterative learning paradigm that enhances the correction policy at each iteration by incorporating knowledge gathered from previous iterations via data-driven confidence metrics assigned to prior corrections.