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

1,3,4-oxadiazoles as inhibitors of the atypical member of the BET family bromodomain factor 3 from Trypanosoma cruzi (TcBDF3)

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, affects millions globally, with increasing urban cases outside of Latin America. Treatment is based on two compounds, namely, benznidazole (BZ) and nifurtimox, but chronic cases pose several challenges. Targeting lysine acetylation, particularly bromodomain-containing proteins, shows promise as a novel antiparasitic target. Our research focuses on TcBDF3, a cytoplasmic protein, which is crucial for parasite differentiation that recognizes acetylated alpha-tubulin. In our previous study, A1B4 was identified as a high-affinity binder of TcBDF3, showing significant trypanocidal activity with low host toxicity in vitro. In this report, the binding of TcBDF3 to A1B4 was validated using differential scanning fluorescence, fluorescence polarization, and molecular modeling, confirming its specific interaction. Additionally, two new 1,3,4-oxadiazoles derived from A1B4 were identified, which exhibited improved trypanocide activity and cytotoxicity profiles. Furthermore, TcBDF3 was classified for the first time as an atypical divergent member of the bromodomain extraterminal family found in protists and plants. These results make TcBDF3 a unique target due to its localization and known functions not shared with higher eukaryotes, which holds promise for Chagas disease treatment.

Updated annotation and meta-analysis of Brugia malayi transcriptomics data reveals consistent transcriptional profiles across time and space with some study-specific differences in adult female worm transcriptional profiles

Genomics, transcriptomics, and proteomics have significantly advanced our understanding of obligately host-associated microbes, where interrogation of the biology is often limited by the complexity of the biological system and limited tools. This includes the causative agents of many neglected tropical diseases, including filarial nematodes. Therefore, numerous transcriptomics studies have been undertaken on filarial nematodes. Most of these transcriptomics studies focus on Brugia malayi, which causes lymphatic filariasis and is a laboratory model for human filarial disease. Here, we undertook a meta-analysis of the publicly available B. malayi transcriptomics data enabling the direct cross comparison of samples from almost a dozen studies. This reanalysis highlights the consistency of transcriptomics results across many different studies and experimental designs from across the globe for over a decade of research, across many different generations of a sequencing technology, library preparation protocols, and differential expression tools. Males and microfilariae across samples had similar expression profiles. However, female samples were clustered into two differential expression patterns that were significantly different from one another. Largely, we confirm previous results for all studies reanalyzed including tissue-specific gene expression and anti-Wolbachia doxycycline treatment of microfilaria. However, we did not detect previously reported differential expression upon in vitro or in vivo treatment with ivermectin, albendazole, and DEC, instead identifying a consistent lack of transcriptomic change upon exposure to these anthelminthic drugs. Updated annotation has been provided that denotes poorly supported genes including those overlapping rRNAs.

Identification of novel bromodomain inhibitors of Trypanosoma cruzi bromodomain factor 2 (TcBDF2) using a fluorescence polarization-based high-throughput assay

Bromodomains are structural folds present in all eukaryotic cells that bind to other proteins recognizing acetylated lysines. Most proteins with bromodomains are part of nuclear complexes that interact with acetylated histone residues and regulate DNA replication, transcription, and repair through chromatin structure remodeling. Bromodomain inhibitors are small molecules that bind to the hydrophobic pocket of bromodomains, interfering with the interaction with acetylated histones. Using a fluorescent probe, we have developed an assay to select inhibitors of the bromodomain factor 2 of Trypanosoma cruzi (TcBDF2) using fluorescence polarization. Initially, a library of 28,251 compounds was screened in an endpoint assay. The top 350-ranked compounds were further analyzed in a dose-response assay. From this analysis, seven compounds were obtained that had not been previously characterized as bromodomain inhibitors. Although these compounds did not exhibit significant trypanocidal activity, all showed bona fide interaction with TcBDF2 with dissociation constants between 1 and 3 µM validating these assays to search for bromodomain inhibitors.

Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

Filarial DAF-12 sense the host serum to resume iL3 development during infection

Nematode parasites enter their definitive host at the developmentally arrested infectious larval stage (iL3), and the ligand-dependent nuclear receptor DAF-12 contributes to trigger their development to adulthood. Here, we characterized DAF-12 from the filarial nematodes Brugia malayi and Dirofilaria immitis and compared them with DAF-12 from the non-filarial nematodes Haemonchus contortus and Caenorhabditis elegans. Interestingly, Dim and BmaDAF-12 exhibit high sequence identity and share a striking higher sensitivity than Hco and CelDAF-12 to the natural ligands Δ4- and Δ7-dafachronic acids (DA). Moreover, sera from different mammalian species activated specifically Dim and BmaDAF-12 while the hormone-depleted sera failed to activate the filarial DAF-12. Accordingly, hormone-depleted serum delayed the commencement of development of D. immitis iL3 in vitro. Consistent with these observations, we show that spiking mouse charcoal stripped-serum with Δ4-DA at the concentration measured in normal mouse serum restores its capacity to activate DimDAF-12. This indicates that DA present in mammalian serum participate in filarial DAF-12 activation. Finally, analysis of publicly available RNA sequencing data from B. malayi showed that, at the time of infection, putative gene homologs of the DA synthesis pathways are coincidently downregulated. Altogether, our data suggest that filarial DAF-12 have evolved to specifically sense and survive in a host environment, which provides favorable conditions to quickly resume larval development. This work sheds new light on the regulation of filarial nematodes development while entering their definitive mammalian host and may open the route to novel therapies to treat filarial infections.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Cost effective, experimentally robust differential-expression analysis for human/mammalian, pathogen and dual-species transcriptomics

As sequencing read length has increased, researchers have quickly adopted longer reads for their experiments. Here, we examine 14 pathogen or host–pathogen differential gene expression data sets to assess whether using longer reads is warranted. A variety of data sets was used to assess what genomic attributes might affect the outcome of differential gene expression analysis including: gene density, operons, gene length, number of introns/exons and intron length. No genome attribute was found to influence the data in principal components analysis, hierarchical clustering with bootstrap support, or regression analyses of pairwise comparisons that were undertaken on the same reads, looking at all combinations of paired and unpaired reads trimmed to 36, 54, 72 and 101 bp. Read pairing had the greatest effect when there was little variation in the samples from different conditions or in their replicates (e.g. little differential gene expression). But overall, 54 and 72 bp reads were typically most similar. Given differences in costs and mapping percentages, we recommend 54 bp reads for organisms with no or few introns and 72 bp reads for all others. In a third of the data sets, read pairing had absolutely no effect, despite paired reads having twice as much data. Therefore, single-end reads seem robust for differential-expression analyses, but in eukaryotes paired-end reads are likely desired to analyse splice variants and should be preferred for data sets that are acquired with the intent to be community resources that might be used in secondary data analyses.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

The Immunoglobulin Superfamily Receptome Defines Cancer-Relevant Networks Associated with Clinical Outcome

Cell surface receptors and their interactions play a central role in physiological and pathological signaling. Despite its clinical relevance, the immunoglobulin superfamily (IgSF) remains uncharacterized and underrepresented in databases. Here, we present a systematic extracellular protein map, the IgSF interactome. Using a high-throughput technology to interrogate most single transmembrane receptors for binding to 445 IgSF proteins, we identify over 500 interactions, 82% previously undocumented, and confirm more than 60 receptor-ligand pairs using orthogonal assays. Our study reveals a map of cell-type-specific interactions and the landscape of dysregulated receptor-ligand crosstalk in cancer, including selective loss of function for tumor-associated mutations. Furthermore, investigation of the IgSF interactome in a large cohort of cancer patients identifies interacting protein signatures associated with clinical outcome. The IgSF interactome represents an important resource to fuel biological discoveries and a framework for understanding the functional organization of the surfaceome during homeostasis and disease, ultimately informing therapeutic development.

A systematic strategy for the investigation of vaccines and drugs targeting bacteria

Infectious and epidemic diseases induced by bacteria have historically caused great distress to people, and have even resulted in a large number of deaths worldwide. At present, many researchers are working on the discovery of viable drug and vaccine targets for bacteria through multiple methods, including the analyses of comparative subtractive genome, core genome, replication-related proteins, transcriptomics and riboswitches, which plays a significant part in the treatment of infectious and pandemic diseases. The 3D structures of the desired target proteins, drugs and epitopes can be predicted and modeled through target analysis. Meanwhile, molecular dynamics (MD) analysis of the constructed drug/epitope-protein complexes is an important standard for testing the suitability of these screened drugs and vaccines. Currently, target discovery, target analysis and MD analysis are integrated into a systematic set of drug and vaccine analysis strategy for bacteria. We hope that this comprehensive strategy will help in the design of high-performance vaccines and drugs.