Transcriptional patterns of sexual dimorphism and in host developmental programs in the model parasitic nematode Heligmosomoides bakeri.

BACKGROUND: Heligmosomoides bakeri (often mistaken for Heligmosomoides polygyrus) is a promising model for parasitic nematodes with the key advantage of being amenable to study and manipulation within a controlled laboratory environment. While draft genome sequences are available for this worm, which allow for comparative genomic analyses between nematodes, there is a notable lack of information on its gene expression. METHODS: We generated biologically replicated RNA-seq datasets from samples taken throughout the parasitic life of H. bakeri. RNA from tissue-dwelling and lumen-dwelling worms, collected under a dissection microscope, was sequenced on an Illumina platform. RESULTS: We find extensive transcriptional sexual dimorphism throughout the fourth larval and adult stages of this parasite and identify alternative splicing, glycosylation, and ubiquitination as particularly important processes for establishing and/or maintaining sex-specific gene expression in this species. We find sex-linked differences in transcription related to aging and oxidative and osmotic stress responses. We observe a starvation-like signature among transcripts whose expression is consistently upregulated in males, which may reflect a higher energy expenditure by male worms. We detect evidence of increased importance for anaerobic respiration among the adult worms, which coincides with the parasite's migration into the physiologically hypoxic environment of the intestinal lumen. Furthermore, we hypothesize that oxygen concentration may be an important driver of the worms encysting in the intestinal mucosa as larvae, which not only fully exposes the worms to their host's immune system but also shapes many of the interactions between the host and parasite. We find stage- and sex-specific variation in the expression of immunomodulatory genes and in anthelmintic targets. CONCLUSIONS: We examine how different the male and female worms are at the molecular level and describe major developmental events that occur in the worm, which extend our understanding of the interactions between this parasite and its host. In addition to generating new hypotheses for follow-up experiments into the worm's behavior, physiology, and metabolism, our datasets enable future more in-depth comparisons between nematodes to better define the utility of H. bakeri as a model for parasitic nematodes in general.

Trickle infection with Heligmosomoides polygyrus results in decreased worm burdens but increased intestinal inflammation and scarring.

Introduction: Intestinal roundworms cause chronic debilitating disease in animals, including humans. Traditional experimental models of these types of infection use a large single-dose infection. However, in natural settings, hosts are exposed to parasites on a regular basis and when mice are exposed to frequent, smaller doses of Heligmosomoides polygyrus, the parasites are cleared more quickly. Whether this more effective host response has any negative consequences for the host is not known. Results: Using a trickle model of infection, we found that worm clearance was associated with known resistance-related host responses: increased granuloma and tuft cell numbers, increased levels of granuloma IgG and decreased intestinal transit time, as well as higher serum IgE levels. However, we found that the improved worm clearance was also associated with an inflammatory phenotype in and around the granuloma, increased smooth muscle hypertrophy/hyperplasia, and elevated levels of Adamts gene expression. Discussion: To our knowledge, we are the first to identify the involvement of this protein family of matrix metalloproteinases (MMPs) in host responses to helminth infections. Our results highlight the delicate balance between parasite clearance and host tissue damage, which both contribute to host pathology. When continually exposed to parasitic worms, improved clearance comes at a cost.

Draft Hybrid Genome Assembly of a Canadian Cyclospora cayetanensis Isolate.

The apicomplexan parasite Cyclospora cayetanensis causes foodborne gastrointestinal disease in humans. Here, we report the first hybrid assembly for C. cayetanensis, which uses both Illumina MiSeq and Oxford Nanopore Technologies MinION platforms to generate genomic sequence data. The final genome assembly consists of 44,586,677 bases represented in 313 contigs.

Characterization of the ß-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance.

BACKGROUND: The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding ß-tubulin proteins, the principal targets of BZ drugs. METHODOLOGY AND PRINCIPAL FINDINGS: First, the ß-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different ß-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations. SIGNIFICANCE: This study demonstrated the presence of at least seven ß-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the ß-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of ß-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.

Benchmarking hybrid assemblies of Giardia and prediction of widespread intra-isolate structural variation.

BACKGROUND: Currently available short read genome assemblies of the tetraploid protozoan parasite Giardia intestinalis are highly fragmented, highlighting the need for improved genome assemblies at a reasonable cost. Long nanopore reads are well suited to resolve repetitive genomic regions resulting in better quality assemblies of eukaryotic genomes. Subsequent addition of highly accurate short reads to long-read assemblies further improves assembly quality. Using this hybrid approach, we assembled genomes for three Giardia isolates, two with published assemblies and one novel, to evaluate the improvement in genome quality gained from long reads. We then used the long reads to predict structural variants to examine this previously unexplored source of genetic variation in Giardia. METHODS: With MinION reads for each isolate, we assembled genomes using several assemblers specializing in long reads. Assembly metrics, gene finding, and whole genome alignments to the reference genomes enabled direct comparison to evaluate the performance of the nanopore reads. Further improvements from adding Illumina reads to the long-read assemblies were evaluated using gene finding. Structural variants were predicted from alignments of the long reads to the best hybrid genome for each isolate and enrichment of key genes was analyzed using random genome sampling and calculation of percentiles to find thresholds of significance. RESULTS: Our hybrid assembly method generated reference quality genomes for each isolate. Consistent with previous findings based on SNPs, examination of heterozygosity using the structural variants found that Giardia BGS was considerably more heterozygous than the other isolates that are from Assemblage A. Further, each isolate was shown to contain structural variant regions enriched for variant-specific surface proteins, a key class of virulence factor in Giardia. CONCLUSIONS: The ability to generate reference quality genomes from a single MinION run and a multiplexed MiSeq run enables future large-scale comparative genomic studies within the genus Giardia. Further, prediction of structural variants from long reads allows for more in-depth analyses of major sources of genetic variation within and between Giardia isolates that could have effects on both pathogenicity and host range.

Genomic analysis of the mesophilic Thermotogae genus Mesotoga reveals phylogeographic structure and genomic determinants of its distinct metabolism.

The genus Mesotoga, the only described mesophilic Thermotogae lineage, is common in mesothermic anaerobic hydrocarbon-rich environments. Besides mesophily, Mesotoga displays lineage-specific phenotypes, such as no or little H2 production and dependence on sulfur-compound reduction, which may influence its ecological role. We used comparative genomics of 18 Mesotoga strains (pairwise 16S rRNA identity >99%) and a transcriptome of M. prima to investigate how life at moderate temperatures affects phylogeography and to interrogate the genomic features of its lineage-specific metabolism. We propose that Mesotoga accomplish H2 oxidation and thiosulfate reduction using a sulfide dehydrogenase and a hydrogenase-complex and that a pyruvate:ferredoxin oxidoreductase acquired from Clostridia is responsible for oxidizing acetate. Phylogenetic analysis revealed three distinct Mesotoga lineages (89.6%-99.9% average nucleotide identity [ANI] within lineages, 79.3%-87.6% ANI between lineages) having different geographic distribution patterns and high levels of intra-lineage recombination but little geneflow between lineages. Including data from metagenomes, phylogeographic patterns suggest that geographical separation historically has been more important for Mesotoga than hyperthermophilic Thermotoga and we hypothesize that distribution of Mesotoga is constrained by their anaerobic lifestyle. Our data also suggest that recent anthropogenic activities and environments (e.g., wastewater treatment, oil exploration) have expanded Mesotoga habitats and dispersal capabilities.

Genomic insights into temperature-dependent transcriptional responses of Kosmotoga olearia, a deep-biosphere bacterium that can grow from 20 to 79 °C.

Temperature is one of the defining parameters of an ecological niche. Most organisms thrive within a temperature range that rarely exceeds ~30 °C, but the deep subsurface bacterium Kosmotoga olearia can grow over a temperature range of 59 °C (20-79 °C). To identify genes correlated with this flexible phenotype, we compared transcriptomes of K. olearia cultures grown at its optimal 65 °C to those at 30, 40, and 77 °C. The temperature treatments affected expression of 573 of 2224 K. olearia genes. Notably, this transcriptional response elicits re-modeling of the cellular membrane and changes in metabolism, with increased expression of genes involved in energy and carbohydrate metabolism at high temperatures and up-regulation of amino acid metabolism at lower temperatures. At sub-optimal temperatures, many transcriptional changes were similar to those observed in mesophilic bacteria at physiologically low temperatures, including up-regulation of typical cold stress genes and ribosomal proteins. Comparative genomic analysis of additional Thermotogae genomes indicates that one of K. olearia's strategies for low-temperature growth is increased copy number of some typical cold response genes through duplication and/or lateral acquisition. At 77 °C one-third of the up-regulated genes are of hypothetical function, indicating that many features of high-temperature growth are unknown.

Draft Genome Sequences of Kosmotoga sp. Strain DU53 and Kosmotoga arenicorallina S304.

Here, we announce the draft genome sequences of two thermophilic Thermotogae bacteria: Kosmotoga sp. strain DU53, isolated from a continental oil reservoir, and Kosmotoga arenicorallina, isolated from hydrothermal sediments. The sequences will provide further insight into evolution of the Kosmotogales.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae.

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.