Dominik R. Laetsch

No evidence for extensive horizontal gene transfer in the genome of the tardigrade Hypsibius dujardini

Significance Tardigrades, also known as moss piglets or water bears, are renowned for their ability to withstand extreme environmental challenges. A recently published analysis of the genome of the tardigrade Hypsibius dujardini by Boothby et al. concluded that horizontal acquisition of genes from bacterial and other sources might be key to cryptobiosis in tardigrades. We independently sequenced the genome of H. dujardini and detected a low level of horizontal gene transfer. We show that the extensive horizontal transfer proposed by Boothby et al. was an artifact of a failure to eliminate contaminants from sequence data before assembly. Tardigrades are meiofaunal ecdysozoans that are key to understanding the origins of Arthropoda. Many species of Tardigrada can survive extreme conditions through cryptobiosis. In a recent paper [Boothby TC, et al. (2015) Proc Natl Acad Sci USA 112(52):15976–15981], the authors concluded that the tardigrade Hypsibius dujardini had an unprecedented proportion (17%) of genes originating through functional horizontal gene transfer (fHGT) and speculated that fHGT was likely formative in the evolution of cryptobiosis. We independently sequenced the genome of H. dujardini. As expected from whole-organism DNA sampling, our raw data contained reads from nontarget genomes. Filtering using metagenomics approaches generated a draft H. dujardini genome assembly of 135 Mb with superior assembly metrics to the previously published assembly. Additional microbial contamination likely remains. We found no support for extensive fHGT. Among 23,021 gene predictions we identified 0.2% strong candidates for fHGT from bacteria and 0.2% strong candidates for fHGT from nonmetazoan eukaryotes. Cross-comparison of assemblies showed that the overwhelming majority of HGT candidates in the Boothby et al. genome derived from contaminants. We conclude that fHGT into H. dujardini accounts for at most 1–2% of genes and that the proposal that one-sixth of tardigrade genes originate from functional HGT events is an artifact of undetected contamination.

The submarine volcano eruption at the island of El Hierro: physical-chemical perturbation and biological response

On October 10 2011 an underwater eruption gave rise to a novel shallow submarine volcano south of the island of El Hierro, Canary Islands, Spain. During the eruption large quantities of mantle-derived gases, solutes and heat were released into the surrounding waters. In order to monitor the impact of the eruption on the marine ecosystem, periodic multidisciplinary cruises were carried out. Here, we present an initial report of the extreme physical-chemical perturbations caused by this event, comprising thermal changes, water acidification, deoxygenation and metal-enrichment, which resulted in significant alterations to the activity and composition of local plankton communities. Our findings highlight the potential role of this eruptive process as a natural ecosystem-scale experiment for the study of extreme effects of global change stressors on marine environments.

Comparative genomics of the tardigrades Hypsibius dujardini and Ramazzottius varieornatus

Tardigrada, a phylum of meiofaunal organisms, have been at the center of discussions of the evolution of Metazoa, the biology of survival in extreme environments, and the role of horizontal gene transfer in animal evolution. Tardigrada are placed as sisters to Arthropoda and Onychophora (velvet worms) in the superphylum Panarthropoda by morphological analyses, but many molecular phylogenies fail to recover this relationship. This tension between molecular and morphological understanding may be very revealing of the mode and patterns of evolution of major groups. Limnoterrestrial tardigrades display extreme cryptobiotic abilities, including anhydrobiosis and cryobiosis, as do bdelloid rotifers, nematodes, and other animals of the water film. These extremophile behaviors challenge understanding of normal, aqueous physiology: how does a multicellular organism avoid lethal cellular collapse in the absence of liquid water? Meiofaunal species have been reported to have elevated levels of horizontal gene transfer (HGT) events, but how important this is in evolution, and particularly in the evolution of extremophile physiology, is unclear. To address these questions, we resequenced and reassembled the genome of H. dujardini, a limnoterrestrial tardigrade that can undergo anhydrobiosis only after extensive pre-exposure to drying conditions, and compared it to the genome of R. varieornatus, a related species with tolerance to rapid desiccation. The 2 species had contrasting gene expression responses to anhydrobiosis, with major transcriptional change in H. dujardini but limited regulation in R. varieornatus. We identified few horizontally transferred genes, but some of these were shown to be involved in entry into anhydrobiosis. Whole-genome molecular phylogenies supported a Tardigrada+Nematoda relationship over Tardigrada+Arthropoda, but rare genomic changes tended to support Tardigrada+Arthropoda.

Massive encapsulation of larval Anguillicoloides crassus in the intestinal wall of Japanese eels

BackgroundWithin the last 25 years, after the introduction of the swimbladder nematode Anguillicoloides crassus from East-Asia to Europe, a body of work has aggregated on the host parasite interactions in the acquired host Anguilla anguilla. Despite the emerging evolutionary interest there is still a lack of knowledge about host parasite relations of A. crassus in its natural host Anguilla japonica. We examined the Anguillicoloides infections of wild-caught Japanese eels as well as from aquacultured specimens in Taiwan with respect to the fate of migratory L3 larvae and performed infection experiments with Japanese eels.ResultsInside the intestinal wall of cultured eels, where the infective pressure was higher than among wild eels, we found large numbers of granuloma-like cysts. In a few eels these cysts contained nematodes still recognizable as L3 larvae of A. crassus, while in most cases the content of these capsules was degraded to amorphous matter. Occurrence of these objects was correlated with the number of encapsulated larvae in the swimbladder wall. We were able to show, that the cysts contained disintegrated L3 larvae by amplification and subsequent sequencing of large subunit ribosomal rRNA. Furthermore we identified repeated infections with high doses of larvae as prerequisites for the processes of encapsulation in infection experiments.ConclusionUnder high infective pressure a large percentage of L3 larvae of A. crassus coming from the gut lumen are eliminated by the natural host within its intestinal tissue. It is possible to reproduce this condition in infection experiments. We provide a fast, easy and reliable PCR-based method for identification of encapsulated swimbladder parasites.

Comparative genomics of apomictic root-knot nematodes : hybridization, ploidy, and dynamic genome change

Abstract The root-knot nematodes (genus Meloidogyne) are important plant parasites causing substantial agricultural losses. The Meloidogyne incognita group (MIG) of species, most of which are obligatory apomicts (mitotic parthenogens), are extremely polyphagous and important problems for global agriculture. While understanding the genomic basis for their variable success on different crops could benefit future agriculture, analyses of their genomes are challenging due to complex evolutionary histories that may incorporate hybridization, ploidy changes, and chromosomal fragmentation. Here, we sequence 19 genomes, representing five species of key root-knot nematodes collected from different geographic origins. We show that a hybrid origin that predated speciation within the MIG has resulted in each species possessing two divergent genomic copies. Additionally, the apomictic MIG species are hypotriploids, with a proportion of one genome present in a second copy. The hypotriploid proportion varies among species. The evolutionary history of the MIG genomes is revealed to be very dynamic, with noncrossover recombination both homogenizing the genomic copies, and acting as a mechanism for generating divergence between species. Interestingly, the automictic MIG species M. floridensis differs from the apomict species in that it has become homozygous throughout much of its genome.

Macrophage origin limits functional plasticity in helminth-bacterial co-infection

Rapid reprogramming of the macrophage activation phenotype is considered important in the defense against consecutive infection with diverse infectious agents. However, in the setting of persistent, chronic infection the functional importance of macrophage-intrinsic adaptation to changing environments vs. recruitment of new macrophages remains unclear. Here we show that resident peritoneal macrophages expanded by infection with the nematode Heligmosomoides polygyrus bakeri altered their activation phenotype in response to infection with Salmonella enterica ser. Typhimurium in vitro and in vivo. The nematode-expanded resident F4/80high macrophages efficiently upregulated bacterial induced effector molecules (e.g. MHC-II, NOS2) similarly to newly recruited monocyte-derived macrophages. Nonetheless, recruitment of blood monocyte-derived macrophages to Salmonella infection occurred with equal magnitude in co-infected animals and caused displacement of the nematode-expanded, tissue resident-derived macrophages from the peritoneal cavity. Global gene expression analysis revealed that although nematode-expanded resident F4/80high macrophages made an anti-bacterial response, this was muted as compared to newly recruited F4/80low macrophages. However, the F4/80high macrophages adopted unique functional characteristics that included enhanced neutrophil-stimulating chemokine production. Thus, our data provide important evidence that plastic adaptation of MΦ activation does occur in vivo, but that cellular plasticity is outweighed by functional capabilities specific to the tissue origin of the cell.

KinFin: Software for Taxon-Aware Analysis of Clustered Protein Sequences

The field of comparative genomics is concerned with the study of similarities and differences between the information encoded in the genomes of organisms. A common approach is to define gene families by clustering protein sequences based on sequence similarity, and analyze protein cluster presence and absence in different species groups as a guide to biology. Due to the high dimensionality of these data, downstream analysis of protein clusters inferred from large numbers of species, or species with many genes, is nontrivial, and few solutions exist for transparent, reproducible, and customizable analyses. We present KinFin, a streamlined software solution capable of integrating data from common file formats and delivering aggregative annotation of protein clusters. KinFin delivers analyses based on systematic taxonomy of the species analyzed, or on user-defined, groupings of taxa, for example, sets based on attributes such as life history traits, organismal phenotypes, or competing phylogenetic hypotheses. Results are reported through graphical and detailed text output files. We illustrate the utility of the KinFin pipeline by addressing questions regarding the biology of filarial nematodes, which include parasites of veterinary and medical importance. We resolve the phylogenetic relationships between the species and explore functional annotation of proteins in clusters in key lineages and between custom taxon sets, identifying gene families of interest. KinFin can easily be integrated into existing comparative genomic workflows, and promotes transparent and reproducible analysis of clustered protein data.

Characterisation of the UK honey bee (Apis mellifera) metagenome.

The European honey bee (Apis mellifera) plays a major role in pollination and food production, but is under threat from emerging pathogens and agro-environmental insults. As with other organisms, honey bee health is a complex product of environment, host genetics and associated microbes (commensal, opportunistic and pathogenic). Improved understanding of bee genetics and their molecular ecology can help manage modern challenges to bee health and production. Sampling bee and cobiont genomes, we characterised the metagenome of 19 honey bee colonies across Britain. Low heterozygosity was observed in bees from many Scottish colonies, sharing high similarity to the native dark bee, A. mellifera mellifera. Apiaries exhibited high diversity in the composition and relative abundance of individual microbiome taxa. Most non-bee sequences derived from known honey bee commensal bacteria or known pathogens, e.g. Lotmaria passim (Trypanosomatidae), and Nosema spp. (Microsporidia). However, DNA was also detected from numerous additional bacterial, plant (food source), protozoan and metazoan organisms. To classify sequences from cobionts lacking genomic information, we developed a novel network analysis approach clustering orphan contigs, allowing the identification of a pathogenic gregarine. Our analyses demonstrate the power of high-throughput, directed metagenomics in agroecosystems identifying potential threats to honey bees present in their microbiota.

An extracellular Argonaute protein mediates export of repeat-associated small RNAs into vesicles in parasitic nematodes

Mobile small RNAs are an integral component of the arms race between plants and fungal parasites, and several studies suggest microRNAs could similarly operate between parasitic nematodes and their animal hosts. However, whether and how specific sequences are selected for export by parasites is unknown. Here we use density gradient purification and proteinase K sensitivity analysis to demonstrate that a specific Argonaute protein (exWAGO) is secreted in extracellular vesicles (EVs) released by the gastrointestinal nematode Heligmosomodies bakeri, at multiple copies per EV. Phylogenetic and gene expression analyses demonstrate exWAGO is highly conserved and abundantly expressed in related parasites, including the human hookworm and proteomic analyses confirm this is the only Argonaute secreted by rodent parasites. In contrast, exWAGO orthologues in species from the free-living genus Caenorhabditis are highly diverged. By re-sequencing and re-annotating the H. bakeri genome, and sequencing multiple small RNA libraries, we determined that the most abundant small RNAs released from the nematode parasite are not microRNAs but rather secondary small interfering RNAs (siRNAs) that are produced by RNA-dependent RNA Polymerases. We further identify distinct evolutionary properties of the siRNAs resident in free-living or parasitic nematodes versus those exported in EVs by the parasite and show that the latter are specifically associated with exWAGO. Together this work identifies an Argonaute protein as a mediator of RNA export and suggests rhabditomorph nematode parasites may have co-opted a novel nematode-unique pathway to communicate with their hosts.Extracellular RNA has been proposed to mediate communication between cells and organisms however relatively little is understood regarding how specific sequences are selected for export. Here we describe a specific Argonaute protein (exWAGO) that is secreted in extracellular vesicles (EVs) released by the gastrointestinal nematode Heligmosomoides bakeri, at multiple copies per EV. Phylogenetic and gene expression analyses demonstrate exWAGO orthologues are highly conserved and abundantly expressed in related parasites but highly diverged in free-living genus Caenorhabditis. We show that the most abundant small RNAs released from the nematode parasite are not microRNAs as previously thought, but rather secondary small interfering RNAs (siRNAs) that are produced by RNA-dependent RNA Polymerases. The siRNAs that are released in EVs have distinct evolutionary properties compared to those resident in free-living or parasitic nematodes. Immunoprecipitation of exWAGO demonstrates that it specifically associates with siRNAs from transposons and newly evolved repetitive elements that are packaged in EVs and released into the host environment. Together this work demonstrates molecular and evolutionary selectivity in the small RNA sequences that are released in EVs into the host environment and identifies a novel Argonaute protein as the mediator of this.

Shared transcriptional control and disparate gain and loss of aphid parasitism genes and loci acquired via horizontal gene transfer

Background Aphids are a diverse group of taxa that contain hundreds of agronomically important species, which vary in their host range and pathogenicity. However, the genome evolution underlying agriculturally important aphid traits is not well understood. Results We generated highly-contiguous draft genome assemblies for two aphid species: the narrow host range Myzus cerasi, and the cereal specialist Rhopalosiphum padi. Using a de novo gene prediction pipeline on both these genome assemblies, and those of three related species (Acyrthosiphon pisum, D. noxia and M. persicae), we show that aphid genomes consistently encode similar gene numbers, and in the case of A. pisum, fewer and larger genes than previously reported. We compare gene content, gene duplication, synteny, horizontal gene transfer events, and putative effector repertoires between these five species to understand the genome evolution of globally important plant parasites. Aphid genomes show signs of relatively distant gene duplication, and substantial, relatively recent, gene birth, and are characterized by disparate gain and loss of genes acquired by horizontal gene transfer (HGT). Such HGT events account for approximately 1% of loci, and contribute to the protein-coding content of aphid species analysed. Putative effector repertoires, originating from duplicated loci, putative HGT events and other loci, have an unusual genomic organisation and evolutionary history. We identify a highly conserved effector-pair that is tightly genetically-linked in all aphid species. In R. padi, this effector pair is tightly transcriptionally-linked, and shares a transcriptional control mechanism with a subset of approximately 50 other putative effectors distributed across the genome. Conclusions This study extends our current knowledge on the evolution of aphid genomes and reveals evidence for a shared control mechanism, which underlies effector expression, and ultimately plant parasitism.