X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Dual RNAseq analyses at soma and germline levels reveal evolutionary innovations in the elephantiasis-agent Brugia malayi, and adaptation of its Wolbachia endosymbionts.

Brugia malayi is a human filarial nematode responsible for elephantiasis, a debilitating condition that is part of a broader spectrum of diseases called filariasis, including lymphatic filariasis and river blindness. Almost all filarial nematode species infecting humans live in mutualism with Wolbachia endosymbionts, present in somatic hypodermal tissues but also in the female germline which ensures their vertical transmission to the nematode progeny. These α-proteobacteria potentially provision their host with essential metabolites and protect the parasite against the vertebrate immune response. In the absence of Wolbachia wBm, B. malayi females become sterile, and the filarial nematode lifespan is greatly reduced. In order to better comprehend this symbiosis, we investigated the adaptation of wBm to the host nematode soma and germline, and we characterized these cellular environments to highlight their specificities. Dual RNAseq experiments were performed at the tissue-specific and ovarian developmental stage levels, reaching the resolution of the germline mitotic proliferation and meiotic differentiation stages. We found that most wBm genes, including putative effectors, are not differentially regulated between infected tissues. However, two wBm genes involved in stress responses are upregulated in the hypodermal chords compared to the germline, indicating that this somatic tissue represents a harsh environment to which wBm have adapted. A comparison of the B. malayi and C. elegans germline transcriptomes reveals a poor conservation of genes involved in the production of oocytes, with the filarial germline proliferative zone relying on a majority of genes absent from C. elegans. The first orthology map of the B. malayi genome presented here, together with tissue-specific expression enrichment analyses, indicate that the early steps of oogenesis are a developmental process involving genes specific to filarial nematodes, that likely result from evolutionary innovations supporting the filarial parasitic lifestyle.

Nearly Complete Genome Sequence of Brugia pahangi FR3.

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi, including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled.

Complete Genome Sequence of wBp, the Wolbachia Endosymbiont of Brugia pahangi FR3.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of wBp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi, from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp.

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.

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.

Drug Repurposing of Bromodomain Inhibitors as Potential Novel Therapeutic Leads for Lymphatic Filariasis Guided by Multispecies Transcriptomics.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi, its Wolbachia endosymbiont wBm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In wBm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti, the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi, expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro, suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis.IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.

Multispecies Transcriptomics Data Set of Brugia malayi, Its Wolbachia Endosymbiont wBm, and Aedes aegypti across the B. malayi Life Cycle.

Here, we present a comprehensive transcriptomics data set of Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host. This study samples from 16 stages across the entire B. malayi life cycle, including stage 1 through 4 larvae, adult males and females, embryos, immature microfilariae, and mature microfilariae.

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r2 = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r2 = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies.

Visualizing Bacterial Colony Morphologies Using Time-Lapse Imaging Chamber MOCHA.

Capturing microbial growth on a macroscopic scale is of great importance to further our understanding of microbial life. However, methods for imaging microbial life on a scale of millimeters to centimeters are often limited by designs that have poor environmental control, resulting in dehydration of the agar plate within just a few days. Here, we created MOCHA (microbial chamber), a simple but effective chamber that allows users to study microbial growth for extended periods (weeks) in a stable environment. Agar hydration is maintained with a double-decker design, in which two glass petri dishes are connected by a wick, allowing the lower plate to keep the upper plate hydrated. This flexible chamber allows the observation of a variety of microbiological phenomena, such as the growth and development of single bacterial and fungal colonies, interspecies interactions, swarming motility, and pellicle formation.IMPORTANCE Detailed study of microbial life on the colony scale of millimeters to centimeters has been lagging considerably behind microscopic inspection of microbes. One major reason for this is the lack of inexpensive instrumentation that can reproducibly capture images in a controlled environment. In this study, we present the design and use of a unique chamber that was used to produce several time-lapse movies that aimed to capture the diversity of microbial colony phenotypes over long periods.

Characterization of innate immunity genes in the parasitic nematode Brugia malayi.

The filarial nematode Brugia malayi is one of the causative agents of lymphatic filariasis, a neglected tropical disease that affects 120 million people worldwide. The limited effectiveness of available anthelmintics and the absence of a vaccine have prompted extensive research on the interaction between Brugia and its obligate bacterial endosymbiont, Wolbachia. Recent studies suggest that Wolbachia is able to manipulate its nematode host immunity but relatively little is known about the immune system of filarial nematodes. Therefore, elucidation of the mechanisms underlying the immune system of B. malayi may be useful for understanding how the symbiotic relationship is maintained and help in the identification of new drug targets. In order to characterize the main genetic pathways involved in B. malayi immunity, we exposed adult female worms to two bacterial lysates (Escherichia coli and Bacillus amyloliquefaciens), dsRNA and dsDNA. We performed transcriptome sequencing of worms exposed to each immune elicitor at two different timepoints. Gene expression analysis of untreated and immune-challenged worms was performed to characterize gene expression patterns associated with each type of immune stimulation. Our results indicate that different immune elicitors produced distinct expression patterns in B. malayi, with changes in the expression of orthologs of well-characterized C. elegans immune pathways such as insulin, TGF-ß, and p38 MAPK pathways, as well as C-type lectins and several stress-response genes.

Genetic Signature of Resistance to White Band Disease in the Caribbean Staghorn Coral Acropora cervicornis.

Coral reefs are declining worldwide due to multiple factors including rising sea surface temperature, ocean acidification, and disease outbreaks. Over the last 30 years, White Band Disease (WBD) alone has killed up to 95% of the Caribbean`s dominant shallow-water corals--the staghorn coral Acropora cervicornis and the elkhorn coral A. palmata. Both corals are now listed on the US Endangered Species Act, and while their recovery has been slow, recent transmission surveys indicate that more than 5% of staghorn corals are disease resistant. Here we compared transcriptome-wide gene expression between resistant and susceptible staghorn corals exposed to WBD using in situ transmission assays. We identified constitutive gene expression differences underlying disease resistance that are independent from the immune response associated with disease exposure. Genes involved in RNA interference-mediated gene silencing, including Argonaute were up-regulated in resistant corals, whereas heat shock proteins (HSPs) were down-regulated. Up-regulation of Argonaute proteins indicates that post-transcriptional gene silencing plays a key, but previously unsuspected role in coral immunity and disease resistance. Constitutive expression of HSPs has been linked to thermal resilience in other Acropora corals, suggesting that the down-regulation of HSPs in disease resistant staghorn corals may confer a dual benefit of thermal resilience.

Evaluating topological conflict in centipede phylogeny using transcriptomic data sets.

Relationships between the five extant orders of centipedes have been considered solved based on morphology. Phylogenies based on samples of up to a few dozen genes have largely been congruent with the morphological tree apart from an alternative placement of one order, the relictual Craterostigmomorpha, consisting of two species in Tasmania and New Zealand. To address this incongruence, novel transcriptomic data were generated to sample all five orders of centipedes and also used as a test case for studying gene-tree incongruence. Maximum likelihood and Bayesian mixture model analyses of a data set composed of 1,934 orthologs with 45% missing data, as well as the 389 orthologs in the least saturated, stationary quartile, retrieve strong support for a sister-group relationship between Craterostigmomorpha and all other pleurostigmophoran centipedes, of which the latter group is newly named Amalpighiata. The Amalpighiata hypothesis, which shows little gene-tree incongruence and is robust to the influence of among-taxon compositional heterogeneity, implies convergent evolution in several morphological and behavioral characters traditionally used in centipede phylogenetics, such as maternal brood care, but accords with patterns of first appearances in the fossil record.

RNA-seq profiles of immune related genes in the staghorn coral Acropora cervicornis infected with white band disease.

Coral diseases are among the most serious threats to coral reefs worldwide, yet most coral diseases remain poorly understood. How the coral host responds to pathogen infection is an area where very little is known. Here we used next-generation RNA-sequencing (RNA-seq) to produce a transcriptome-wide profile of the immune response of the Staghorn coral Acropora cervicornis to White Band Disease (WBD) by comparing infected versus healthy (asymptomatic) coral tissues. The transcriptome of A. cervicornis was assembled de novo from A-tail selected Illumina mRNA-seq data from whole coral tissues, and parsed bioinformatically into coral and non-coral transcripts using existing Acropora genomes in order to identify putative coral transcripts. Differentially expressed transcripts were identified in the coral and non-coral datasets to identify genes that were up- and down-regulated due to disease infection. RNA-seq analyses indicate that infected corals exhibited significant changes in gene expression across 4% (1,805 out of 47,748 transcripts) of the coral transcriptome. The primary response to infection included transcripts involved in macrophage-mediated pathogen recognition and ROS production, two hallmarks of phagocytosis, as well as key mediators of apoptosis and calcium homeostasis. The strong up-regulation of the enzyme allene oxide synthase-lipoxygenase suggests a key role of the allene oxide pathway in coral immunity. Interestingly, none of the three primary innate immune pathways--Toll-like receptors (TLR), Complement, and prophenoloxydase pathways, were strongly associated with the response of A. cervicornis to infection. Five-hundred and fifty differentially expressed non-coral transcripts were classified as metazoan (n = 84), algal or plant (n = 52), fungi (n = 24) and protozoans (n = 13). None of the 52 putative Symbiodinium or algal transcript had any clear immune functions indicating that the immune response is driven by the coral host, and not its algal symbionts.