<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2024.1418032</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Dual RNA-seq in filarial nematodes and <italic>Wolbachia</italic> endosymbionts using RNase H based ribosomal RNA depletion</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Cantin</surname> <given-names>Lindsey J.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2425356/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Gregory</surname> <given-names>Vanessa</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2733896/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Blum</surname> <given-names>Laura N.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Foster</surname> <given-names>Jeremy M.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1018378/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Biochemistry and Microbiology Division, New England BioLabs</institution>, <addr-line>Ipswich, MA</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Applications and Product Development, New England BioLabs</institution>, <addr-line>Ipswich, MA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Takema Fukatsu, National Institute of Advanced Industrial Science and Technology (AIST), Japan</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Kenneth Pfarr, University Hospital Bonn, Germany</p>
<p>Yunsheng Wang, Hunan Agricultural University, China</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Lindsey J. Cantin <email>lcantin&#x00040;neb.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>20</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1418032</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>04</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>05</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2024 Cantin, Gregory, Blum and Foster.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Cantin, Gregory, Blum and Foster</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium <italic>Wolbachia</italic> 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&#x02013;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 <italic>Wolbachia</italic> 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 <italic>Brugia malayi</italic> containing the <italic>Wolbachia</italic> symbiont (<italic>w</italic>Bm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5&#x02032;-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 <italic>Wolbachia</italic> mRNA reads, resulting in a 16&#x02013;1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the <italic>Brugia malayi</italic> depletion probes can remove rRNA from the filarial worm <italic>Dirofilaria immitis</italic> and the majority of rRNA from the more distantly related free living nematode <italic>Caenorhabditis elegans</italic>. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.</p></abstract>
<kwd-group>
<kwd>filariasis</kwd>
<kwd><italic>Wolbachia</italic></kwd>
<kwd>dual RNA sequencing</kwd>
<kwd>symbiosis</kwd>
<kwd>RNA enrichment</kwd>
<kwd>ribosomal depletion</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="1"/>
<equation-count count="0"/>
<ref-count count="75"/>
<page-count count="14"/>
<word-count count="9964"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Microbial Symbioses</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>Filariasis is an infectious disease caused by parasitic nematodes. The disease is one of the leading causes of disability worldwide, with an estimated 2 billion people at risk of infection in the tropical and sub-tropical regions of Africa, Asia, and South America (Taylor et al., <xref ref-type="bibr" rid="B60">2010</xref>; Bhalla et al., <xref ref-type="bibr" rid="B9">2013</xref>). Filarial parasites are transmitted through blood feeding insect vectors and can infect both humans and animals (Chandy et al., <xref ref-type="bibr" rid="B13">2011</xref>). The pathology of the disease, from host preference to disease outcomes, is determined by the filarial species present. Lymphatic filariasis can lead to elephantiasis and disfigurement, which is caused by <italic>Brugia malayi</italic> (<italic>B. malayi</italic>) and <italic>Wuchereria bancrofti</italic> infections in humans (World Health Organization, <xref ref-type="bibr" rid="B70">2019</xref>; Medeiros et al., <xref ref-type="bibr" rid="B41">2021</xref>). Subcutaneous filariasis is caused by <italic>Loa loa, Onchocerca volvulus</italic>, and <italic>Mansonella streptocerca</italic> infections and can lead to skin issues and blindness (WHO Expert Committee on Onchocerciasis Control, <xref ref-type="bibr" rid="B69">1993</xref>; Boussinesq, <xref ref-type="bibr" rid="B10">2006</xref>; Mediannikov and Ranque, <xref ref-type="bibr" rid="B42">2018</xref>). <italic>Dirofilaria immitis</italic> (<italic>D. immitis</italic>), also known as heartworm, is a filarial nematode that primarily infects canid species, leading to lung and heart damage (Sim&#x000F3;n et al., <xref ref-type="bibr" rid="B53">2012</xref>; Anvari et al., <xref ref-type="bibr" rid="B5">2020</xref>). Combinations of albendazole, ivermectin, and diethylcarbamazine are used to control the spread of the disease by targeting microfilariae (Campbell, <xref ref-type="bibr" rid="B12">1991</xref>; Richards et al., <xref ref-type="bibr" rid="B50">2001</xref>; Mackenzie et al., <xref ref-type="bibr" rid="B37">2002</xref>; Molyneux et al., <xref ref-type="bibr" rid="B44">2003</xref>). However, drugs that kill adult stages are necessary to treat active infections, as the worms can survive long-term in the vertebrate hosts.</p>
<p>The majority of filarial nematodes contain <italic>Wolbachia</italic> as an obligate intracellular bacterial symbiont at all stages of their life cycle (Mclaren et al., <xref ref-type="bibr" rid="B40">1975</xref>; Sironi et al., <xref ref-type="bibr" rid="B54">1995</xref>; Bandi et al., <xref ref-type="bibr" rid="B6">1998</xref>, <xref ref-type="bibr" rid="B7">2001</xref>; Quek et al., <xref ref-type="bibr" rid="B47">2022</xref>). <italic>Wolbachia</italic> are critical to the development, reproduction, and long-term survival of adult nematodes, most likely by providing metabolites from biological pathways absent in the worms (Foster et al., <xref ref-type="bibr" rid="B21">2005</xref>; Wu et al., <xref ref-type="bibr" rid="B71">2009</xref>; Darby et al., <xref ref-type="bibr" rid="B18">2012</xref>; Li and Carlow, <xref ref-type="bibr" rid="B34">2012</xref>; Luck et al., <xref ref-type="bibr" rid="B36">2016</xref>). However, the molecular basis for this symbiosis is still poorly understood. Anti-<italic>Wolbachia</italic> treatments with antibiotics, such as doxycycline and rifampicin, are a promising strategy for filariasis eradication by sterilizing and killing adult worms (Taylor et al., <xref ref-type="bibr" rid="B59">2005</xref>, <xref ref-type="bibr" rid="B61">2019</xref>; Bazzocchi et al., <xref ref-type="bibr" rid="B8">2008</xref>; Hoerauf et al., <xref ref-type="bibr" rid="B25">2008</xref>; Specht et al., <xref ref-type="bibr" rid="B56">2008</xref>; Supali et al., <xref ref-type="bibr" rid="B58">2008</xref>; Coulibaly et al., <xref ref-type="bibr" rid="B17">2009</xref>; Mand et al., <xref ref-type="bibr" rid="B38">2009</xref>; Wanji et al., <xref ref-type="bibr" rid="B68">2009</xref>; Johnston et al., <xref ref-type="bibr" rid="B27">2014</xref>, <xref ref-type="bibr" rid="B28">2021</xref>; Aljayyoussi et al., <xref ref-type="bibr" rid="B2">2017</xref>; Clare et al., <xref ref-type="bibr" rid="B16">2019</xref>). Gene expression analysis for both the nematode and the bacteria during different life-cycle stages can help to identify the biological relationship between them and to determine putative drug targets.</p>
<p>RNA-sequencing is a powerful approach that can be used to identify and quantify genes expressed under certain conditions or throughout the life cycle (Wang et al., <xref ref-type="bibr" rid="B66">2009</xref>; Hrdlickova et al., <xref ref-type="bibr" rid="B26">2017</xref>). Ribosomal RNA (rRNA) makes up between 80 and 90% of total RNA in a cell (Zhao et al., <xref ref-type="bibr" rid="B74">2018</xref>; Deng et al., <xref ref-type="bibr" rid="B19">2022</xref>). Therefore, rRNA must be removed prior to sequencing to analyze the smaller fraction of messenger RNA (mRNA) and non-coding RNAs (ncRNA) involved in interesting cellular functions. Typically, rRNA is removed during library preparation using Poly(A) enrichment or ribosomal depletion methods (Sultan et al., <xref ref-type="bibr" rid="B57">2014</xref>; Zhao et al., <xref ref-type="bibr" rid="B75">2014</xref>; Schuierer et al., <xref ref-type="bibr" rid="B52">2017</xref>). Poly(A) enrichment is commonly used when studying eukaryotes by selecting for RNA species containing a poly(A) tail, usually mRNA and some ncRNAs, while removing all others, including rRNA. However, poly(A) tails are mostly absent from bacterial RNAs, resulting in a loss of all bacterial signal from libraries prepared with this method (Prezza et al., <xref ref-type="bibr" rid="B46">2020</xref>; Xiang et al., <xref ref-type="bibr" rid="B72">2022</xref>).</p>
<p>Ribosomal depletion can be used to remove only the rRNA, allowing for eukaryotic mRNA, bacterial mRNA, and ncRNA to all be sequenced within the same library. Ribosomal depletion typically requires sequence specific probes for either the pulldown or targeted degradation of rRNA (Adiconis et al., <xref ref-type="bibr" rid="B1">2013</xref>; Wahl et al., <xref ref-type="bibr" rid="B64">2022</xref>). Commonly studied model organisms, such as mouse, rat, human and bacteria, have commercially available kits for ribosomal depletion. There are limited depletion options for filarial nematodes and other <italic>Wolbachia</italic> containing hosts, preventing dual RNA-sequencing of the hosts and bacterial endosymbionts together (Kumar et al., <xref ref-type="bibr" rid="B31">2012</xref>, <xref ref-type="bibr" rid="B32">2016</xref>; Grote et al., <xref ref-type="bibr" rid="B23">2017</xref>; Chung et al., <xref ref-type="bibr" rid="B15">2019</xref>). 5&#x02032; phosphate dependent exonucleases, such as Terminator and Xrn1, have also been used to selectively degrade large rRNA molecules containing a 5&#x02032; monophosphate (He et al., <xref ref-type="bibr" rid="B24">2010</xref>; Kang et al., <xref ref-type="bibr" rid="B29">2011</xref>; Grote et al., <xref ref-type="bibr" rid="B23">2017</xref>; Wangsanuwat et al., <xref ref-type="bibr" rid="B67">2020</xref>). However, this method does not remove the 5S rRNA and is greatly affected by the secondary structure of RNAs, lowering its efficiency and increasing off-target effects compared to probe-based depletion.</p>
<p>Here, we present ribosomal depletion probes for <italic>B. malayi</italic> and its <italic>Wolbachia</italic> endosymbiont, <italic>w</italic>Bm, to perform dual RNA-seq. We also provide a proof of principal for the updated publicly available NEBNext Custom RNA Depletion design tool (<ext-link ext-link-type="uri" xlink:href="https://depletion-design.neb.com">https://depletion-design.neb.com</ext-link>), incorporating a clustering step to remove redundancy in the probe set. We compared libraries prepared with our custom filarial nematode depletion probes and RNase H digestion to RNA-seq libraries prepared with no treatment (total RNA), poly(A) enrichment, 5&#x02032; phosphate dependent exonuclease digestion, and commercially available NEBNext rRNA depletion with a combination of probes designed for Human/Mouse/Rat and Bacteria. We find almost complete rRNA removal and an increase in <italic>Wolbachia</italic> mRNA reads with the <italic>B. malayi</italic> specific probes, improving upon results seen with the other methods mentioned above. The custom depletion probes were also tested for their ability to remove rRNA from total RNA of the closely related filarial nematode <italic>D. immitis</italic> and the distantly related free-living nematode <italic>Caenorhabditis elegans</italic> (<italic>C. elegans</italic>). We found an almost 70% reduction in rRNA reads in the <italic>C. elegans</italic> library and an over 99% reduction in the <italic>D. immitis</italic> library, showing that this custom probe design can be used as a pan-filariae ribosomal depletion method for dual RNA-seq with their <italic>Wolbachia</italic> endosymbionts.</p></sec>
<sec sec-type="materials and methods" id="s2">
<title>2 Materials and methods</title>
<sec>
<title>2.1 <italic>B. malayi</italic> tissue collection</title>
<p>Thirty live adult <italic>B. malayi</italic> females were received from the NIH-NIAID Research Reagent Resource Center (FR3, Athens, GA; Michalski et al., <xref ref-type="bibr" rid="B43">2011</xref>). Immediately after arrival, the worms were placed at 37&#x000B0;C with 5% CO<sub>2</sub> in freshly prepared RPMI 1640 media containing 10% heat inactivated fetal bovine serum (Thermo Fisher Scientific, Waltham, MA), 2 mM L-glutamine, 5 g/L glucose, 100 ug/mL streptomycin, 100 U/mL penicillin and 250 ng/mL amphotericin (Millipore Sigma, Burlington, MA). The worms were incubated overnight for a minimum of 16 h. The worms were sorted into groups of 3 and washed twice with 1X PBS. All washes were removed, and the worms were snap frozen in liquid nitrogen and stored in 1.5 mL LoBind tubes at &#x02212;80&#x000B0;C.</p>
</sec>
<sec>
<title>2.2 <italic>D. immitis</italic> tissue collection</title>
<p>Five live adult female <italic>D. immitis</italic> were shipped overnight from the FR3. Upon arrival, the worms were rinsed two times in 1X PBS. All liquid was removed. The worms were placed individually in 1.5 mL DNA LoBind tubes and flash frozen in liquid nitrogen. The samples were stored at &#x02212;80&#x000B0;C until future use.</p>
</sec>
<sec>
<title>2.3 <italic>C. elegans</italic> tissue collection</title>
<p>N2 strain <italic>C. elegans</italic> nematodes were maintained on 10 cm Nematode Growth Medium (NGM) agar plates seeded with OP50 <italic>E. coli</italic> (Chaudhuri et al., <xref ref-type="bibr" rid="B14">2011</xref>). Once the plate was starved, the worms were recovered by rinsing the plates with 10 mL of M9 buffer (22 mM KH<sub>2</sub>PO<sub>4</sub>, 42 mM Na<sub>2</sub>HPO<sub>4</sub>, 86 mM NaCl) and transferring the liquid containing the worms to a 15 mL conical tube. The worms were pelleted by spinning at 400 &#x000D7; <italic>g</italic> for 5 min. To wash the worm pellet, the supernatant was removed and another 10 mL of M9 buffer was added to the tube to resuspend the pellet. The worms were washed a total of 5 times. After the final spin, 9 mL of the supernatant was discarded. The worm pellet was resuspended in the remaining buffer and transferred to a 1.5 mL DNA LoBind tube. The tube was spun at 400 &#x000D7; <italic>g</italic> for 5 min and the supernatant was removed. The <italic>C. elegans</italic> pellet was flash frozen in liquid nitrogen and stored at &#x02212;80&#x000B0;C until future use.</p>
</sec>
<sec>
<title>2.4 Total RNA extraction, sequencing, and analysis</title>
<p>Total RNA was extracted separately from individual aliquots for each nematode species, including <italic>B. malayi</italic> (three adult female worms), <italic>C. elegans</italic> (worm pellet from a starved 10 cm plate, mixed sex), and <italic>D. immitis</italic> (one adult female worm). For each sample, 0.5 mL of TRIzol (Thermo Fisher Scientific, Waltham, MA) was added to the 1.5 mL LoBind tube containing the frozen tissue. The tissue and TRIzol mix were then moved to a sterile 2 mL Dounce Homogenizer. The nematodes were homogenized using 10 strokes of the A pestle and 10&#x02013;20 strokes of the B pestle. The samples were incubated at room temperature for 5 min. We added 100 &#x003BC;L of chloroform to each sample and mixed by manually shaking for 10 s. After an additional 3-min incubation at room temperature, the aqueous and organic phases were separated by spinning at 12,000 &#x000D7; <italic>g</italic> for 15 min at 4&#x000B0;C. The upper clear aqueous phase containing the RNA was transferred to a clean 1.5 mL DNA LoBind tube. The volume in each tube was measured with a P1000 pipette and an equal volume of &#x0003E;95% ethanol was added. We used the Monarch Total RNA Miniprep Kit (New England BioLabs, Ipswich, MA) to complete the cleanup of the RNA, using the total RNA purification from TRIzol-extracted samples protocol as described, with the inclusion of the optional DNase I treatment. RNA quality was measured by running 1 &#x003BC;L of the purified RNA sample on a Bioanalyzer using the RNA 6000 pico kit (Agilent, Santa Clara, CA) with the Eukaryote Total RNA Pico Series II assay. All RNA used for further experiments had an RNA Integrity Number (RIN) of 8 or above.</p>
<p>Total RNA-seq libraries for each species were prepared using the NEBNext Ultra II RNA Library Prep Kit for Illumina (New England BioLabs, Ipswich, MA) using the manufacturer&#x00027;s instructions. We had three technical replicates for <italic>B. malayi</italic> and one replicate for the <italic>D. immitis</italic> and <italic>C. elegans</italic> samples. Briefly, RNA samples containing 10 ng of total RNA were fragmented for 15 min at 94&#x000B0;C prior to first and second strand synthesis. The cDNA was purified with 1.8X NEBNext Sample Purification Beads. The cDNA ends were repaired and ligated to Illumina adapters diluted 25-fold. After adapter ligation, the USER enzyme was used to selectively degrade the second cDNA strand, allowing for directional RNA-seq. The libraries were amplified by PCR using NEBNext Multiplex Oligos for Illumina (Dual Index Primers Set 1; New England BioLabs, Ipswich, MA), with a total of 10 cycles. The libraries were purified a final time using 0.9X NEBNext Sample Purification Beads and 1 uL was run on the Bioanalyzer using the High Sensitivity DNA Kit (Agilent, Santa Clara, CA) for validation and quantification. The libraries were pooled and sequenced on a single flow cell of the Illumina NextSeq 550 at a depth of at least 40 million paired-end reads per sample.</p>
<p>Each analysis tool used default options unless otherwise noted. The FASTQ reads were trimmed to remove adapter sequences with the &#x02013;paired option using Cutadapt (v1.16; Martin, <xref ref-type="bibr" rid="B39">2011</xref>). The quality of each library after trimming was calculated using FastQC (v0.11.9; Andrews, <xref ref-type="bibr" rid="B4">2023</xref>). The trimmed reads were aligned to the genome for each respective species using HISAT2 (v2.1.0; Kim et al., <xref ref-type="bibr" rid="B30">2019</xref>) with the &#x02013;fr &#x02013;rna-strandedness RF and -k 30 options. The trimmed reads were mapped to their reference genome, using the <italic>B. malayi</italic>-4.0 (GCA_000002995.5), ICBAS_JMDir_1.0 (GCA_024305405.1), and WBcel235 (GCA_000002985.3) genomes for <italic>B. malayi, D. immitis</italic>, and <italic>C. elegans</italic> reads, respectively. The SAM alignment files were converted to BAM files, keeping only the primary alignments and mapped reads using Samtools view -b -F 260 (v1.15.1; Li et al., <xref ref-type="bibr" rid="B33">2009</xref>). Samtools view was also used to down-sample to 40 million paired reads per library. BamCoverage from deepTools (v3.5.1; Ram&#x000ED;rez et al., <xref ref-type="bibr" rid="B49">2014</xref>) was used to create bedGraph files with 50 base pair (bp) bins, skipping all 0 values. Bedtools (v2.29.2; Quinlan and Hall, <xref ref-type="bibr" rid="B48">2010</xref>) merge with the options &#x0201C;-d 1 -c 4 -o sum&#x0201D; was used to combine the bedGraph bins directly beside one another. To identify rRNA regions with a high number of aligned sequencing reads, we kept all genomic regions with &#x0003E;1,000 reads. The regions were aligned to the nucleotide (nt) database using BLASTN (v2.13.0; Altschul et al., <xref ref-type="bibr" rid="B3">1990</xref>; Camacho et al., <xref ref-type="bibr" rid="B11">2009</xref>), to identify high coverage regions that match available rRNA sequences. A bed file was created for each species containing the genomic coordinates for regions that match rRNA sequences with BLAST and any genes that were assigned as rRNA in the reference annotation. The genomic sequences of these regions, representing the putative rRNA sequences, were obtained using the reference genomes and Bedtools getfasta. The <italic>B. malayi</italic> rRNA sequences were used to design custom rRNA depletion probes.</p>
</sec>
<sec>
<title>2.5 Ribosomal RNA probe design</title>
<p>Putative ribosomal RNA sequences from <italic>B. malayi</italic> and <italic>w</italic>Bm were fed into the NEBNext Custom RNA Depletion Design Tool (<ext-link ext-link-type="uri" xlink:href="https://depletion-design.neb.com">https://depletion-design.neb.com</ext-link>), a public tool which utilizes ChipD to select candidate probes at each position based on predicted melting temperature (Dufour et al., <xref ref-type="bibr" rid="B20">2010</xref>). The result is a set of probes which closely tile the reference sequences. The initial probe set (<italic>n</italic> = 2,020) was clustered to collapse highly similar probes using vsearch where pairwise similarity was defined by (<italic>number of matching columns</italic>)/(<italic>the shortest sequence length</italic>; &#x02013;cluster_fast &#x02013;iddef 0 &#x02013;qmask none; Rognes et al., <xref ref-type="bibr" rid="B51">2016</xref>).</p>
<p>Several similarity thresholds (50, 60, and 70%) were evaluated before selecting the probe set clustered at 60% similarity (<italic>n</italic> = 377). This similarity threshold was chosen because it had ample tiling over all rRNA sequences, with the additional advantage of being synthesized in a single oligo pool (max per pool = 384 oligos). Candidate probe sets were aligned back to the putative rRNA sequences with bbmap (v39.06; <ext-link ext-link-type="uri" xlink:href="https://sourceforge.net/projects/bbmap">https://sourceforge.net/projects/bbmap</ext-link>) to assess coverage. Probe efficacy was predicted <italic>in silico</italic> using bbduk to partition undepleted RNA reads by whether they had high k-mer similarity to the probe set or not. The selected probe set had high similarity to 88.69% of RNA reads, compared to 88.73% for the unclustered output. Following the development of this method, the NEBNext Custom RNA Depletion Design Tool was updated to include a clustering option to support the use of the tool with less well-characterized organisms. We ordered an oligo pool from IDT (Coralville, IA) containing the 377 <italic>B. malayi</italic> rRNA depletion probes at a concentration of 50 pmol/oligo. The oligo pool was resuspended in 25 &#x003BC;L of 10 mM Tris, 0.1 mM EDTA, at pH 7.5 (2 &#x003BC;M per probe final concentration) and stored at &#x02212;20&#x000B0;C until future use.</p>
</sec>
<sec>
<title>2.6 RNA enrichment and sequencing</title>
<p>The total RNA, described in Section 2.4, was used to make <italic>B. malayi</italic> RNA libraries prepared with various ribosomal depletion methods. We used 10 ng of total RNA as the starting material for each library. One sample was treated with Terminator 5&#x02032;-Phosphate-Dependant Exonuclease (Biosearch Technologies, Teddington, UK) using the manufacturer&#x00027;s standard procedure with Reaction Buffer A. The sample was incubated at 30&#x000B0;C for 60 min and the reaction was terminated by adding 1 &#x003BC;L of 100 mM EDTA. The sample was purified using the Monarch RNA Cleanup Kit (50 ug; New England BioLabs, Ipswich, MA), with a final elution of 6 &#x003BC;L.</p>
<p>We also tested different probe sets for RNase H based ribosomal depletion. The NEBNext rRNA Depletion Kit v2 for Human/Mouse/Rat in combination with the NEBNext rRNA depletion probes for bacteria (both from New England BioLabs, Ipswich, MA) were used to deplete the rRNA from <italic>B. malayi</italic> total RNA. We diluted 10 ng of RNA in nuclease-free water to final volume of 9 &#x003BC;L. We added 2 &#x003BC;L of the probe hybridization buffer, 2 &#x003BC;L of the human/mouse/rat v2 rRNA depletion solution and 2 &#x003BC;L of the bacteria rRNA depletion solution. The sample was then treated following the manufacturer&#x00027;s instructions. Briefly, the probes were annealed to the rRNA by incubating at 95&#x000B0;C for 2 min followed by a ramp down to 22&#x000B0;C at a rate of 0.1&#x000B0;C per second. The sample was then treated with RNase H for 30 min at 50&#x000B0;C to degrade the RNA in the RNA/DNA hybrids, followed by a DNase I incubation for 30 min at 37&#x000B0;C. The remaining RNA was purified using 1.8X of the NEBNext RNA Sample Purification Beads and eluted in 7 &#x003BC;L of nuclease-free water.</p>
<p>For the third RNA depletion method, we tested our custom <italic>B. malayi</italic> probe design with three technical replicates. We added 2 &#x003BC;L of our custom oligo pool, described in Section 2.5, to 2 &#x003BC;L of Hybridization Buffer from the NEBNext RNA Depletion Core Reagent Set (New England BioLabs, Ipswich, MA) per replicate. The hybridization mix was then added to 10 ng of total <italic>B. malayi</italic> RNA in 11 &#x003BC;L of nuclease-free water. The samples were then treated following manufacturer&#x00027;s instructions and as described above. We also used the custom <italic>B. malayi</italic> probes on <italic>C. elegans</italic> and <italic>D. immitis</italic> total RNA to test efficacy across nematode species. Both of these libraries started with 10 ng of total RNA, described in Section 2.4, and followed the same RNase H protocol described above using the NEBNext RNA Depletion Core Reagent Set.</p>
<p>Five &#x003BC;L from each elution for all three depletion methods, including the <italic>C. elegans</italic> and <italic>D. immitis</italic> samples, were taken into library preparation using the NEBNext Ultra II RNA Library Prep Kit for Illumina, following the same methods described in Section 2.4. All seven depletion libraries were pooled and sequenced using the Illumina NextSeq 550 to a depth of at least 40 million paired-end reads.</p>
</sec>
<sec>
<title>2.7 RNA-seq analysis</title>
<p>We used default settings for all analysis tools unless noted otherwise. All <italic>B. malayi</italic> RNA-seq libraries, including the 3 total (untreated), Terminator, Human/Mouse/Rat/Bacteria depleted, 3 custom depleted, and a publicly available Poly(A) enriched library (SRR3111490; Chung et al., <xref ref-type="bibr" rid="B15">2019</xref>), were trimmed using Cutadapt (v1.16; Martin, <xref ref-type="bibr" rid="B39">2011</xref>) with the &#x02013;paired option and aligned to the <italic>B. malayi</italic>-4.0 (GCA_000002995.5) reference genome, containing the <italic>Wolbachia</italic> endosymbiont (<italic>w</italic>Bm) chromosome, using HISAT2 (v2.1.0; Kim et al., <xref ref-type="bibr" rid="B30">2019</xref>) with the &#x0201C;&#x02013;fr &#x02013;rna-strandedness RF and -k 30&#x0201D; options. Samtools view (v1.15.1; Li et al., <xref ref-type="bibr" rid="B33">2009</xref>) -b -F 260 was used on each alignment file to convert to a BAM file and remove all non-primary alignments and unmapped reads. As each library has a different number of sequencing reads, the BAM files were subsampled to 40 million reads each by dividing 40,000,000 by the total number of reads and using this value as the downsampling factor with the -s option. To visualize the read alignments, we used bamCoverage from DeepTools (v3.5.1; Ram&#x000ED;rez et al., <xref ref-type="bibr" rid="B49">2014</xref>) to make bigwig files from the filtered BAMs with 50 bp bins and normalization using counts per million (CPM). The bigwig files were visualized using IGV (v2.11.9; Thorvaldsd&#x000F3;ttir et al., <xref ref-type="bibr" rid="B62">2013</xref>).</p>
<p>RSeQC (v2.6.4.3; Wang et al., <xref ref-type="bibr" rid="B65">2012</xref>) geneBody_coverage.py was used to visualize the distribution of reads across gene bodies to determine if there is a 5&#x02032; or 3&#x02032; bias with different depletion methods. FeatureCounts (v2.0.1; Liao et al., <xref ref-type="bibr" rid="B35">2014</xref>) was used to count the number of reads over annotated genomic regions, including protein coding, non-coding and ribosomal genes, with the &#x0201C;-s 2 &#x02013;primary -t gene -g &#x02018;gene_id&#x0201D;&#x00027; options set. An additional column was added to the count output files, where each gene id was matched to the assigned &#x0201C;gene_biotype&#x0201D; from the annotation file. The sum of the reads assigned to each gene type was calculated and then divided by the total number of reads to quantify the proportion of reads mapped to rRNA, protein coding, ncRNA, or unannotated regions. The feature count files were also used to plot the correlation between the gene counts for the total vs. the custom depletion libraries. Fragments per million (FPM) was calculated for all three replicates for each library type by first taking the sum of all reads in each sample divided by 1 million to obtain the scaling factor, unique to each individual replicate. The counts for each gene were divided by the scaling factor and the average FPM was calculated for the total and custom depletion replicate libraries, separately. The log10 of the average FPMs was taken and visualized using ggplot2 (Valero-Mora, <xref ref-type="bibr" rid="B63">2010</xref>). Correlations between the custom probe library and the Poly(A) or Terminator library were performed using the same methods, without averaging the Poly(A) and Terminator counts, as only one replicate was used.</p>
<p>Samtools view -L was used to identify which reads mapped to the <italic>w</italic>Bm chromosome using bed files for the whole chromosome and for just the <italic>w</italic>Bm rRNA regions. The reads mapping to the <italic>w</italic>Bm rRNA genes were subtracted from the total number of reads mapping to the whole chromosome, where the remaining reads correspond to <italic>Wolbachia</italic> reads mapping to genes of interest.</p>
<p>To analyze the rRNA depletion capacity of our custom probes on different nematode species, the <italic>C. elegans</italic> and <italic>D. immitis</italic> total and custom depletion libraries were compared. The same trimming, alignment, filtering and downsampling methods were used with their respective reference genomes: WBcel235 (GCA_000002985.3) for <italic>C. elegans</italic> and ICBAS_JMDir_1.0 (GCA_024305405.1) for <italic>D. immitis</italic>. BamCoverage was also used to make bigWig files with 50 bp bins and CPM normalization for visualization in the IGV browser. Samtools view -L was used to count reads mapping to the rRNA region bed files described in Section 2.4.</p></sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec>
<title>3.1 RNA enrichment methods</title>
<p>We tested various RNA-seq enrichment and depletion methods in three nematode species, two of which contain <italic>Wolbachia</italic> as a bacterial endosymbiont. The different library preparation methods should result in sequencing of distinct RNA species (<xref ref-type="fig" rid="F1">Figure 1</xref>, <xref ref-type="table" rid="T1">Table 1</xref>). Sequencing of <italic>B. malayi</italic> total untreated RNA should result in 80&#x02013;95% ribosomal sequences, with a small number of mRNA, ncRNA, and bacterial sequences (<xref ref-type="fig" rid="F1">Figure 1A</xref>). This would require large amounts of deep sequencing data to quantitatively study the protein coding gene expression. Poly(A) enrichment is one of the most commonly used methods for eukaryotic RNA-sequencing. This usually involves using Oligo(d)T beads to bind and pull-down the poly(A) tail of certain types of RNA molecules, including mRNA and some types of ncRNA, while removing all rRNA and bacterial RNAs (<xref ref-type="fig" rid="F1">Figure 1B</xref>). However, poly(A) enrichment will remove all bacterial sequences.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Schematic of RNA enrichment and depletion strategies tested in this study. <bold>(A)</bold> Total RNA sequencing. <bold>(B)</bold> Poly(A) enrichment with magnetic d(T) beads. <bold>(C)</bold> Terminator 5&#x02032; Phosphate Dependent Exonuclease treatment. <bold>(D)</bold> RNase H-based ribosomal depletion using sequence specific DNA probes. RNA species are color-coded with rRNA in dark blue, mRNA in teal, and ncRNA in yellow. The Terminator exonuclease is shown in purple, RNase H in green and DNA probes in red.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1418032-g0001.tif"/>
</fig>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Expected sequencing outcomes from RNA enrichment and depletion methods.</p></caption>
<table frame="box" rules="all">
<thead>
<tr style="background-color:#919498;color:#ffffff">
<th valign="top" align="left"><bold>Library preparation</bold></th>
<th valign="top" align="left"><bold>RNA species</bold></th>
<th valign="top" align="center"><bold>Bacterial mRNA</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Total (no depletion)</td>
<td valign="top" align="left">mRNA, ncRNA, and all rRNA</td>
<td valign="top" align="center">&#x0002B;</td>
</tr> <tr>
<td valign="top" align="left">Poly(A)</td>
<td valign="top" align="left">mRNA</td>
<td valign="top" align="center">-</td>
</tr> <tr>
<td valign="top" align="left">Terminator exonuclease</td>
<td valign="top" align="left">mRNA, ncRNA, and 5S rRNA</td>
<td valign="top" align="center">&#x0002B;</td>
</tr> <tr>
<td valign="top" align="left">Ribosomal depletion</td>
<td valign="top" align="left">mRNA and ncRNA</td>
<td valign="top" align="center">&#x0002B;</td>
</tr></tbody>
</table>
</table-wrap>
<p>The removal of rRNA should enable the user to sequence all ncRNAs and mRNAs, from both eukaryotic and bacteria together in one sequencing library. This can be done using species specific depletion probes or with a phosphate-dependent exonuclease. Here, we tested Terminator 5&#x02032;-Phosphate Dependent Exonuclease treatment prior to library preparation, which should remove all large rRNA molecules, while retaining mRNA, ncRNA and 5S rRNA, from both the nematodes and <italic>Wolbachia</italic> (<xref ref-type="fig" rid="F1">Figure 1C</xref>). Species specific depletion probes can be designed to remove certain RNA molecules from the total RNA pool based on their nucleotide sequence. Here, we used DNA probes and RNase H to selectively degrade RNA molecules within an RNA/DNA hybrid (<xref ref-type="fig" rid="F1">Figure 1D</xref>). When using DNA probes with sequences complementary to rRNA, we expect to sequence mRNA and ncRNA from both eukaryotes and bacteria. Comparing the sequencing outcomes of these library preparation methods allowed us to determine the most efficient technique for dual RNA-sequencing in samples containing both eukaryotic and prokaryotic RNAs.</p>
</sec>
<sec>
<title>3.2 Ribosomal depletion in <italic>B. malayi</italic></title>
<p>We extracted high quality total RNA from <italic>B. malayi</italic> with an RNA Integrity of 9.2 (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1A</xref>). Two distinct peaks can be seen in the Bioanalyzer trace representing the abundant 18S and 28S rRNAs. After depletion with the custom (Nema) probe set, we observed a significant reduction in the ribosomal peaks (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1B</xref>). We prepared and sequenced five <italic>B. malayi</italic> RNA-seq libraries using the methods mentioned above. After sequencing, we aligned our reads to the <italic>B. malayi</italic> reference genome (Foster et al., <xref ref-type="bibr" rid="B22">2020</xref>). The read coverage was mapped over annotated genes. No significant 5&#x02032; or 3&#x02032; bias was observed in any of our libraries, as the quality of the starting total RNA was high (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1C</xref>), although, the Terminator library has a slight 5&#x02032; bias.</p>
<p>There are a sizeable number of reads mapping to the ribosomal tandem repeat genes in the total RNA sample (<xref ref-type="fig" rid="F2">Figure 2A</xref>). The Terminator library also has a pile up of reads over this region. However, there is a slightly different pattern of read peaks, as some of these regions have been depleted. The NEBNext Human/Mouse/Rat/Bacterial depletion probes remove a large portion of the rRNA reads, leaving 400&#x02013;600 bp regions undepleted. The Poly(A) and Nema probes both remove almost all of the ribosomal reads. When we observe the read alignments at a smaller scale, we can see reads aligned to protein coding genes surrounding the tandem repeat in every sample besides the total RNA (<xref ref-type="fig" rid="F2">Figure 2B</xref>). However, the Nema probes and Poly(A) libraries have the highest level of mRNA signal. Additionally, we find virtually no reads aligned to the rRNA genes in the Nema probes library, while there appears to be some rRNA carryover in the Poly(A) library.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Ribosomal depletion across library preparation methods in <italic>B. malayi</italic>. <bold>(A)</bold> IGV trace showing the custom probes and RNA-seq libraries mapped to the ribosomal tandem repeat region on <italic>B. malayi</italic> Chr2. All BAM coverage files were normalized using CPM in 50 bp bins and shown at a scale of 0&#x02013;7,500 reads. Protein coding gene annotations are shown in teal and ribosomal gene annotations are shown in dark blue. <bold>(B)</bold> IGV trace showing the same genomic region and coverage files as <bold>(A)</bold> with a scale of 0&#x02013;50 reads. <bold>(C)</bold> Proportion of reads mapped to each gene annotation type, with unannotated regions in gray, non-coding genes in yellow, protein coding genes in teal. H/M/R/B probes = library made with NEBNext Human/Mouse/Rat rRNA depletion kit combined with the NEBNext Bacteria rRNA depletion probes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1418032-g0002.tif"/>
</fig>
<p>Each read was assigned to a genomic feature using the <italic>B. malayi</italic> reference genome annotation. We quantified the proportion of reads that map to each feature type for each library preparation method (<xref ref-type="fig" rid="F2">Figure 2C</xref>). All annotated genes, besides protein coding and rRNA genes, were considered non-coding genes. Unannotated reads refer to those mapped to regions that do not contain any known genomic features. RNA-seq experiments often aim to study protein coding (mRNA) or non-coding (ncRNA) gene expression. Around 81% of reads in the total library mapped to rRNA genes, which makes the large majority of the reads unusable for downstream analyses. The Poly(A) library had 1.2% of reads mapping to rRNA and almost no reads mapping to ncRNAs. Unsurprisingly, this results in a substantial increase in mRNA reads. Therefore, Poly(A) enrichment is very useful for studying protein coding genes but will not allow for the study of ncRNAs. The Terminator library still had 33% of reads mapping to rRNA and an increase in reads mapping to unannotated regions, with only 32% of reads mapping to protein coding genes. The Human/Mouse/Rat/Bacteria depletion library had 27% of reads mapping to rRNA genes. While this value is significantly lower than in the total RNA library, there are still a large number of rRNAs that need to be removed when using probes designed for distantly related species. The Nema depletion library had the lowest number of rRNA reads at 0.23%. Notably, this value was even lower than the number of rRNA reads found in the Poly(A) depletion library. In this case, the RNase H method is 5 times more efficient at reducing rRNA carryover than using d(T) magnetic beads. The Nema depletion library also had 6% of reads mapping to ncRNA genes, allowing for the examination of these less studied but functionally important RNA species.</p>
</sec>
<sec>
<title>3.3 Enrichment of <italic>w</italic>Bm RNA-seq reads after ribosomal depletion</title>
<p>The Nema depletion probes allow for a more comprehensive analysis of <italic>B. malayi</italic> gene expression. Next, we wanted to determine if this method can also be used to enhance studies of <italic>Wolbachia</italic> gene expression. The <italic>B. malayi</italic> reference genome contains the <italic>w</italic>Bm chromosome. When looking at the reads mapping to the bacteria, we find many reads mapping to rRNA genes in the total and Terminator libraries (<xref ref-type="fig" rid="F3">Figure 3A</xref>). The Terminator enzyme does not appear to have any effective degradation of the bacterial rRNA. When looking at a smaller scale, we observe some rRNA reads present in the Poly(A) library as well (<xref ref-type="fig" rid="F3">Figure 3B</xref>). The Nema depletion library has virtually no reads mapping to the rRNA regions but does contain expression signals in surrounding protein coding genes. Finally, the Human/Mouse/Rat/Bacteria depletion probes do remove the <italic>w</italic>Bm rRNAs, but this does not appear to increase the number of mRNA reads present.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Ribosomal depletion across library preparation methods in <italic>w</italic>Bm. <bold>(A)</bold> IGV trace showing the custom probes and RNA-seq libraries aligned to an rRNA annotated region in the <italic>w</italic>Bm chromosome. All BAM coverage files were normalized using CPM in 50 bp bins. The scale is set to 0&#x02013;2,000 reads. Protein coding gene annotations are shown in teal and ribosomal gene annotations are shown in dark blue. <bold>(B)</bold> IGV trace showing the same regions and tracks as <bold>(A)</bold>, with a scale of 0&#x02013;25 reads. <bold>(C)</bold> Percentage of reads mapped to the <italic>w</italic>Bm chromosome from all mapped reads in each library. Reads mapped to rRNA genes are shown in dark blue and reads mapped to all other regions of the <italic>w</italic>Bm chromosome are shown in teal. H/M/R/B probes = library made with NEBNext Human/Mouse/Rat rRNA depletion kit combined with the NEBNext Bacteria rRNA depletion probes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1418032-g0003.tif"/>
</fig>
<p>We extracted all reads that mapped to the <italic>w</italic>Bm chromosome from each library. The reads were classified as either rRNA or other <italic>Wolbachia</italic> genes, representing the protein coding and non-coding genes. The proportion of reads mapping to either gene type was compared across libraries (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 2</xref>). The total, Poly(A), and Terminator libraries had over 98% of <italic>Wolbachia</italic> reads mapping to ribosomal genes. Both probe depletion libraries had the opposite proportion, with over 99% of reads mapping outside of the rRNA genes. Next, we determined the percentage of reads mapping to the <italic>w</italic>Bm chromosome out of all the reads within each library (<xref ref-type="fig" rid="F3">Figure 3C</xref>). The total and Terminator libraries had the highest number of <italic>Wolbachia</italic> reads, with around 2 and 3%, respectively. However, these reads do not provide useful information about <italic>Wolbachia</italic> biology as they are primarily rRNA. The Human/Mouse/Rat/Bacteria depletion library had only 0.07% of total reads mapped to <italic>Wolbachia</italic>. The low number of overall <italic>Wolbachia</italic> reads using this method does not allow for accurate quantification of bacterial gene expression. Finally, the Nema depletion library had around 1% of all reads mapping to <italic>Wolbachia</italic>, with almost all of the reads mapping outside of rRNA regions. These data show that the custom filarial nematode depletion probes are the best option for performing rRNA depletion for dual RNA-seq between the worms and their bacterial endosymbiont.</p>
</sec>
<sec>
<title>3.4 Correlation of gene expression between nema depletion and total RNA-seq</title>
<p>We investigated the effect of library preparation method on individual gene counts. All gene counts were normalized by calculating the number of fragments per million mapped fragments (FPM). The genes were not normalized for length, as we directly compared the same genes across enrichment and depletion methods. We plotted the Log<sub>10</sub>(FPM) of each annotated gene for the Total and Nema depleted libraries (<xref ref-type="fig" rid="F4">Figure 4</xref>). All gene types, excluding rRNAs, generally had higher read counts in the Nema deplete library. As expected, most rRNA genes had between 100 and 300 times more gene counts in the Total library and no rRNA genes had higher read counts in the Nema depleted sample.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Correlation of gene counts between <italic>B. malayi</italic> Total and Nema depletion libraries. Plot showing the log10 FPM gene counts for the total (x-axis) and custom probe depleted (y-axis) libraries, with rRNA genes in dark blue, protein coding genes in teal, non-coding RNAs in yellow and <italic>Wolbachia</italic> non-rRNA genes in pink.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1418032-g0004.tif"/>
</fig>
<p>There are some protein coding and non-coding genes that appear to have slightly higher read counts in the Total library. We investigated these genes to ensure that they were not being depleted as a result of off target probe binding. We found no evidence of sequence similarity between these genes and the rRNA probes, as the probe sequences do not map to their genome regions <italic>in silico</italic>. Most genes, such as Bm15518, have the same distribution of reads across the gene body, but with slightly lower coverage in the Nema library (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 3A</xref>). This could be the result of technical artifacts during library preparation and sequencing, as many of these genes had variable gene counts across technical replicates. Bm9361 was the only gene where the distribution of reads differed across the gene body (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 3B</xref>). In the Nema deplete library, a 200 nucleotide sequence around the 3&#x02032; end of the gene appears to be depleted. This may represent an off-target effect of the probe set. However, the overwhelming majority of genes have higher gene counts in the Nema deplete library, without any off-target probe binding.</p>
<p>We also compared the Nema depletion library gene counts to those of the Poly(A) and Terminator libraries. The Nema depletion and Poly(A) libraries had similar levels of protein coding gene expression (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 4A</xref>). The Poly(A) library had higher counts in most rRNA genes, almost zero reads in <italic>Wolbachia</italic> genes and minimal reads in most ncRNA genes. The non-coding genes with similar counts in the Poly(A) and Nema depletion libraries probably represent the ncRNAs that contain a poly(A) tail. The Terminator library has higher counts in all rRNA genes, with less enrichment of protein coding and <italic>Wolbachia</italic> genes (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 4B</xref>). There is an enrichment of a small subset of non-coding and protein coding genes in the Terminator library. This could be the result of technical artifacts from sequencing or biological artifacts as a result of their cap or secondary structures.</p>
</sec>
<sec>
<title>3.5 Ribosomal depletion in other nematode species</title>
<p><italic>C. elegans</italic> is a commonly studied model organism in clade V of the nematode phylum (Parkinson et al., <xref ref-type="bibr" rid="B45">2004</xref>). <italic>D. immitis</italic> is a more closely related filarial parasite in the same group, clade III, as <italic>B. malayi</italic>. It is estimated that clades III and V separated over 300 million years ago (Xie et al., <xref ref-type="bibr" rid="B73">2022</xref>). We determined the efficiency of the depletion probes designed with <italic>B. malayi</italic> sequences on these other nematode species. Total and Nema depleted libraries were made with RNA from each species. The libraries were aligned to their species&#x00027; respective reference genomes, along with the <italic>B. malayi</italic> probe sequences. We observed read and probe alignments over the annotated rRNA genes. In <italic>C. elegans</italic>, the probes have a patchy alignment to the rRNA genes, with many gaps present (<xref ref-type="fig" rid="F5">Figure 5A</xref>). We still see rRNA depletion, however there are two distinct regions, of 500 and 250 bps, where the probes do not align and the rRNA is not depleted. The probes have better alignment to the <italic>D. immitis</italic> rRNA genes and almost all rRNA appears to be removed in the Nema depleted library (<xref ref-type="fig" rid="F5">Figure 5B</xref>). When observing read alignment over protein coding genes, we find an increase in reads for the Nema depleted libraries in both species (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figures 5A</xref>, <xref ref-type="supplementary-material" rid="SM1">B</xref>).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Ribosomal depletion with <italic>B. malayi</italic> custom depletion probes in other nematode species. <bold>(A)</bold> IGV track showing probe mapping in the <italic>C. elegans</italic> genome over ribosomal genes. Additional tracks include BAM coverage files normalized using CPM showing <italic>C. elegans</italic> total RNA (blue) and Nema depletion (pink) libraries at a scale of 0&#x02013;50,000 reads. <bold>(B)</bold> IGV track showing probe mapping in the <italic>D. immitis</italic> genome over ribosomal genes. Additional tracks show BAM coverage files normalized using CPM of the <italic>D. immitis</italic> total RNA (blue) and Nema depletion (pink) libraries at a scale of 0&#x02013;25,000. Protein coding gene annotations are shown in teal and rRNA annotations are shown in dark blue. <bold>(C)</bold> Proportion of reads mapping to ribosomal genes vs. other genomic regions in the <italic>C. elegans</italic> and <italic>D. immitis</italic> total and Nema depletion libraries. Reads mapping to ribosomal genes are represented in dark blue. Reads mapping to other genomic regions are shown in teal.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-15-1418032-g0005.tif"/>
</fig>
<p>We quantified the proportion of reads aligning to rRNA gene annotations from each library (<xref ref-type="fig" rid="F5">Figure 5C</xref>). All reads mapping to ribosomal genes were called rRNA and the remaining mapped reads were classified as other. The <italic>C. elegans</italic> total RNA library had 92.9% of reads mapped to rRNA, with only 7.1% of reads mapped to other genes. The Nema depletion library had 29.5% of reads mapping to rRNA and 70.5% mapped to other regions, indicating that the <italic>B. malayi</italic> probes are not as effective at removing these sequences in this distantly related nematode species. In the <italic>D. immitis</italic> total library, 83% of reads mapped to annotated rRNA genes, with 17% mapping to other genes. The Nema depletion library had only 0.3% of reads mapping to rRNA regions and 99.7% mapping to other regions. The <italic>B. malayi</italic> rRNA probes efficiently removed these sequences at a similar rate across both filarial species. This probe design can most likely work as a pan-filarial nematode ribosomal depletion method.</p></sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>Enrichment and depletion methods in RNA-sequencing are important for studying biologically interesting gene expression in species across life cycle stages and conditions. Here, we presented ribosomal depletion strategies to study filarial parasites and their endosymbiotic bacteria. We tested poly(A) enrichment, Terminator exonuclease depletion, and species-specific probe depletion, including commercially available and custom probes, on nematode total RNA samples. The custom probes were designed against the <italic>B. malayi</italic> and <italic>w</italic>Bm ribosomal sequences using the NEBNext Custom RNA Depletion design tool, resulting in 377 non-redundant oligo sequences. We obtained complete removal of rRNA using these probes in <italic>B. malayi</italic> samples, which was a significant improvement compared to all other methods tested. There was less rRNA contamination with this RNase H-based method when compared to poly(A) bead enrichment. We found that the probes have negligible off-target effects on other gene sequences. The Terminator enzyme is advertised as a method to remove large rRNA sequences. While we did observe a reduction in the rRNA reads, over 30% of reads still mapped to ribosomal genes.</p>
<p>The Terminator method also appears to enrich and deplete certain protein coding and non-coding genes. Commercially available 5&#x02032; Phosphate-dependent exonucleases have been shown to be affected by the secondary structure of RNA. This method also had a slight 5&#x02032; sequencing bias when compared to the other library preparation methods (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1C</xref>). These data show that probe based ribosomal depletion methods should be preferred over current 5&#x02032; Phosphate-dependent exonuclease methods.</p>
<p>The NEBNext Human/Mouse/Rat probes combined with the NEBNext Bacteria probes also removed a large proportion of ribosomal reads. However, there were specific regions of these genes that were not removed, leaving 30% of reads in rRNA genes. These regions likely correspond to sequences that have diverged significantly over time between vertebrates and nematodes. The addition of the Bacteria probes did efficiently remove the <italic>Wolbachia</italic> rRNAs.</p>
<p>The custom probes allowed for sequencing of non-coding and bacterial RNAs alongside the protein-coding RNAs. One limitation of this method is the necessity of an annotated reference genome or a total RNA-seq library to identify the abundant sequences. However, the nema depletion was the only method that resulted in abundant levels of <italic>Wolbachia</italic> sequences, not including rRNA reads. Therefore, this method can be used for dual RNA-seq studies in filarial nematodes containing <italic>Wolbachia</italic>. Sequencing both the eukaryotic and prokaryotic sequences in the same library will help to avoid technical artifacts and batch effects that could arise across <italic>B. malayi</italic> and <italic>Wolbachia</italic> libraries that must otherwise be prepared separately, using different methods. Additionally, these probes can be used to enrich for only <italic>Wolbachia</italic> sequences by combining rRNA depletion with poly(A) depletion. The probe-based rRNA depletion can also help to improve reference genome annotations of non-coding genes, as many of these sequences are removed when commonly used poly(A) enrichment is performed. This library preparation method is helpful in situations where the RNA is of low quality. Degraded RNAs are difficult to sequence using either poly(A) enrichment or Terminator depletion, as these methods rely on intact 3&#x02032; and 5&#x02032; ends, respectively. The probes are designed to bind over all regions of the rRNA, therefore those sequences can be removed even if the whole RNA molecule is fragmented.</p>
<p>The custom <italic>B. malayi</italic> probes did not completely remove rRNAs from the distantly related <italic>C. elegans</italic> samples. A small number of additional probes could be added to the pool to increase efficiency. However, the <italic>B. malayi</italic> probes were able to completely remove rRNAs from <italic>D. immitis</italic> samples. <italic>B. malayi</italic> and <italic>D. immitis</italic> are in distinct groups in the Filarioidea superfamily (Small et al., <xref ref-type="bibr" rid="B55">2014</xref>). Therefore, we expect that these probes can be used as a high-quality ribosomal depletion method across all filarial parasites and facilitate efficient dual RNA-seq of those species harboring the <italic>Wolbachia</italic> endosymbiont.</p></sec>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: <ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/">https://www.ncbi.nlm.nih.gov/</ext-link>, PRJNA1100530.</p></sec>
<sec sec-type="ethics-statement" id="s6">
<title>Ethics statement</title>
<p>The manuscript presents research on animals that do not require ethical approval for their study.</p></sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>LC: Writing &#x02013; review &#x00026; editing, Writing &#x02013; original draft, Visualization, Validation, Supervision, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. VG: Writing &#x02013; review &#x00026; editing, Writing &#x02013; original draft, Visualization, Methodology, Investigation, Formal analysis, Data curation. LB: Writing &#x02013; review &#x00026; editing, Writing &#x02013; original draft, Software, Methodology, Formal analysis, Conceptualization. JF: Writing &#x02013; review &#x00026; editing, Writing &#x02013; original draft, Supervision, Funding acquisition, Conceptualization.</p></sec>
</body>
<back>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<ack><p>We thank New England BioLabs for supporting parasitology research and the summer student internship program, particularly Thomas Evans, Salvatore Russello, and the late Donald Comb. We would also like to thank the NEB sequencing core for processing the samples. Finally, we acknowledge the Filariasis Research Reagent Resource Center for supplying the filarial nematodes used in this study.</p>
</ack>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>LC, VG, LB, and JF were employed by New England BioLabs.</p>
<p>The reviewer KP declared a past collaboration with the author JF to the handling editor.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s10">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2024.1418032/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmicb.2024.1418032/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Presentation_1.PPTX" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.presentationml.presentation" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_1.XLSX" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.XLSX" id="SM3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Adiconis</surname> <given-names>X.</given-names></name> <name><surname>Borges-Rivera</surname> <given-names>D.</given-names></name> <name><surname>Satija</surname> <given-names>R.</given-names></name> <name><surname>DeLuca</surname> <given-names>D. S.</given-names></name> <name><surname>Busby</surname> <given-names>M. A.</given-names></name> <name><surname>Berlin</surname> <given-names>A. M.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Comparative analysis of RNA sequencing methods for degraded or low-input samples</article-title>. <source>Nat. Methods</source> <volume>10</volume>, <fpage>623</fpage>&#x02013;<lpage>629</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.2483</pub-id><pub-id pub-id-type="pmid">23685885</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aljayyoussi</surname> <given-names>G.</given-names></name> <name><surname>Tyrer</surname> <given-names>H. E.</given-names></name> <name><surname>Ford</surname> <given-names>L.</given-names></name> <name><surname>Sjoberg</surname> <given-names>H.</given-names></name> <name><surname>Pionnier</surname> <given-names>N.</given-names></name> <name><surname>Waterhouse</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Short-course, high-dose rifampicin achieves <italic>Wolbachia</italic> depletion predictive of curative outcomes in preclinical models of lymphatic filariasis and onchocerciasis</article-title>. <source>Sci. Rep</source>. <volume>7</volume>:<fpage>210</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-00322-5</pub-id><pub-id pub-id-type="pmid">28303006</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname> <given-names>S. F.</given-names></name> <name><surname>Gish</surname> <given-names>W.</given-names></name> <name><surname>Miller</surname> <given-names>W.</given-names></name> <name><surname>Myers</surname> <given-names>E. W.</given-names></name> <name><surname>Lipman</surname> <given-names>D. J.</given-names></name></person-group> (<year>1990</year>). <article-title>Basic local alignment search tool</article-title>. <source>J. Mol. Biol</source>. <volume>215</volume>, <fpage>403</fpage>&#x02013;<lpage>410</lpage>. <pub-id pub-id-type="doi">10.1016/S0022-2836(05)80360-2</pub-id><pub-id pub-id-type="pmid">2231712</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Andrews</surname> <given-names>S.</given-names></name></person-group> (<year>2023</year>). <source>FastQC</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://github.com/s-andrews/FastQC">https://github.com/s-andrews/FastQC</ext-link> (accessed October 31, 2023).</citation>
</ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anvari</surname> <given-names>D.</given-names></name> <name><surname>Narouei</surname> <given-names>E.</given-names></name> <name><surname>Daryani</surname> <given-names>A.</given-names></name> <name><surname>Sarvi</surname> <given-names>S.</given-names></name> <name><surname>Moosazadeh</surname> <given-names>M.</given-names></name> <name><surname>Ziaei Hezarjaribi</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>The global status of <italic>Dirofilaria immitis</italic> in dogs: a systematic review and meta-analysis based on published articles</article-title>. <source>Res. Vet. Sci</source>. <volume>131</volume>, <fpage>104</fpage>&#x02013;<lpage>116</lpage>. <pub-id pub-id-type="doi">10.1016/j.rvsc.2020.04.002</pub-id><pub-id pub-id-type="pmid">32330696</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bandi</surname> <given-names>C.</given-names></name> <name><surname>Anderson</surname> <given-names>T. J.</given-names></name> <name><surname>Genchi</surname> <given-names>C.</given-names></name> <name><surname>Blaxter</surname> <given-names>M. L.</given-names></name></person-group> (<year>1998</year>). <article-title>Phylogeny of <italic>Wolbachia</italic> in filarial nematodes</article-title>. <source>Proc. R Soc. B Biol. Sci</source>. <volume>265</volume>, <fpage>2407</fpage>&#x02013;<lpage>2413</lpage>. <pub-id pub-id-type="doi">10.1098/rspb.1998.0591</pub-id><pub-id pub-id-type="pmid">9921679</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bandi</surname> <given-names>C.</given-names></name> <name><surname>Trees</surname> <given-names>A. J.</given-names></name> <name><surname>Brattig</surname> <given-names>N. W.</given-names></name></person-group> (<year>2001</year>). <article-title><italic>Wolbachia</italic> in filarial nematodes: evolutionary aspects and implications for the pathogenesis and treatment of filarial diseases</article-title>. <source>Vet. Parasitol</source>. <volume>98</volume>, <fpage>215</fpage>&#x02013;<lpage>238</lpage>. <pub-id pub-id-type="doi">10.1016/S0304-4017(01)00432-0</pub-id><pub-id pub-id-type="pmid">11516587</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bazzocchi</surname> <given-names>C.</given-names></name> <name><surname>Mortarino</surname> <given-names>M.</given-names></name> <name><surname>Grandi</surname> <given-names>G.</given-names></name> <name><surname>Kramer</surname> <given-names>L. H.</given-names></name> <name><surname>Genchi</surname> <given-names>C.</given-names></name> <name><surname>Bandi</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Combined ivermectin and doxycycline treatment has microfilaricidal and adulticidal activity against <italic>Dirofilaria immitis</italic> in experimentally infected dogs</article-title>. <source>Int. J. Parasitol</source>. <volume>38</volume>, <fpage>1401</fpage>&#x02013;<lpage>1410</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijpara.2008.03.002</pub-id><pub-id pub-id-type="pmid">18433753</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Bhalla</surname> <given-names>D.</given-names></name> <name><surname>Dumas</surname> <given-names>M.</given-names></name> <name><surname>Preux</surname> <given-names>P. M.</given-names></name></person-group> (<year>2013</year>). <article-title>&#x0201C;Chapter 18&#x02014;neurological manifestations of filarial infections,&#x0201D;</article-title> in <source>Handbook of Clinical Neurology</source>, eds. H. H. Garcia, H. B. Tanowitz, and O. H. Del Brutto (Elsevier), <fpage>235</fpage>&#x02013;<lpage>242</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.sciencedirect.com/science/article/pii/B9780444534903000182">https://www.sciencedirect.com/science/article/pii/B9780444534903000182</ext-link> (accessed March 25, 2024).</citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boussinesq</surname> <given-names>M.</given-names></name></person-group> (<year>2006</year>). <article-title>Loiasis</article-title>. <source>Ann. Trop. Med. Parasitol</source>. <volume>100</volume>, <fpage>715</fpage>&#x02013;<lpage>731</lpage>. <pub-id pub-id-type="doi">10.1179/136485906X112194</pub-id><pub-id pub-id-type="pmid">17227650</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Camacho</surname> <given-names>C.</given-names></name> <name><surname>Coulouris</surname> <given-names>G.</given-names></name> <name><surname>Avagyan</surname> <given-names>V.</given-names></name> <name><surname>Ma</surname> <given-names>N.</given-names></name> <name><surname>Papadopoulos</surname> <given-names>J.</given-names></name> <name><surname>Bealer</surname> <given-names>K.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>BLAST&#x0002B;: architecture and applications</article-title>. <source>BMC Bioinform</source>. <volume>10</volume>:<fpage>421</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2105-10-421</pub-id><pub-id pub-id-type="pmid">20003500</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Campbell</surname> <given-names>W. C.</given-names></name></person-group> (<year>1991</year>). <article-title>Ivermectin as an antiparasitic agent for use in humans</article-title>. <source>Annu. Rev. Microbiol</source>. <volume>45</volume>, <fpage>445</fpage>&#x02013;<lpage>474</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.mi.45.100191.002305</pub-id><pub-id pub-id-type="pmid">1741621</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chandy</surname> <given-names>A.</given-names></name> <name><surname>Thakur</surname> <given-names>A. S.</given-names></name> <name><surname>Singh</surname> <given-names>M. P.</given-names></name> <name><surname>Manigauha</surname> <given-names>A.</given-names></name></person-group> (<year>2011</year>). <article-title>A review of neglected tropical diseases: filariasis</article-title>. <source>Asian Pac. J. Trop. Med</source>. <volume>4</volume>, <fpage>581</fpage>&#x02013;<lpage>586</lpage>. <pub-id pub-id-type="doi">10.1016/S1995-7645(11)60150-8</pub-id><pub-id pub-id-type="pmid">21803313</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chaudhuri</surname> <given-names>J.</given-names></name> <name><surname>Parihar</surname> <given-names>M.</given-names></name> <name><surname>Pires-daSilva</surname> <given-names>A.</given-names></name></person-group> (<year>2011</year>). <article-title>An introduction to worm lab: from culturing worms to mutagenesis</article-title>. <source>J. Vis. Exp.</source> <volume>47</volume>:<fpage>2293</fpage>. <pub-id pub-id-type="doi">10.3791/2293</pub-id><pub-id pub-id-type="pmid">21248706</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chung</surname> <given-names>M.</given-names></name> <name><surname>Teigen</surname> <given-names>L. E.</given-names></name> <name><surname>Libro</surname> <given-names>S.</given-names></name> <name><surname>Bromley</surname> <given-names>R. E.</given-names></name> <name><surname>Olley</surname> <given-names>D.</given-names></name> <name><surname>Kumar</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Drug repurposing of bromodomain inhibitors as potential novel therapeutic leads for lymphatic filariasis guided by multispecies transcriptomics</article-title>. <source>mSystems</source> <volume>4</volume>:<fpage>19</fpage>. <pub-id pub-id-type="doi">10.1128/mSystems.00596-19</pub-id><pub-id pub-id-type="pmid">31796568</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clare</surname> <given-names>R. H.</given-names></name> <name><surname>Bardelle</surname> <given-names>C.</given-names></name> <name><surname>Harper</surname> <given-names>P.</given-names></name> <name><surname>Hong</surname> <given-names>W. D.</given-names></name> <name><surname>B&#x000F6;rjesson</surname> <given-names>U.</given-names></name> <name><surname>Johnston</surname> <given-names>K. L.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Industrial scale high-throughput screening delivers multiple fast acting macrofilaricides</article-title>. <source>Nat. Commun</source>. <volume>10</volume>:<fpage>11</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-018-07826-2</pub-id><pub-id pub-id-type="pmid">30602718</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Coulibaly</surname> <given-names>Y. I.</given-names></name> <name><surname>Dembele</surname> <given-names>B.</given-names></name> <name><surname>Diallo</surname> <given-names>A. A.</given-names></name> <name><surname>Lipner</surname> <given-names>E. M.</given-names></name> <name><surname>Doumbia</surname> <given-names>S. S.</given-names></name> <name><surname>Coulibaly</surname> <given-names>S. Y.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>A randomized trial of doxycycline for <italic>Mansonella perstans</italic> infection</article-title>. <source>N. Engl. J. Med</source>. <volume>361</volume>, <fpage>1448</fpage>&#x02013;<lpage>1458</lpage>. <pub-id pub-id-type="doi">10.1056/NEJMoa0900863</pub-id><pub-id pub-id-type="pmid">19812401</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Darby</surname> <given-names>A. C.</given-names></name> <name><surname>Armstrong</surname> <given-names>S. D.</given-names></name> <name><surname>Bah</surname> <given-names>G. S.</given-names></name> <name><surname>Kaur</surname> <given-names>G.</given-names></name> <name><surname>Hughes</surname> <given-names>M. A.</given-names></name> <name><surname>Kay</surname> <given-names>S. M.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Analysis of gene expression from the <italic>Wolbachia</italic> genome of a filarial nematode supports both metabolic and defensive roles within the symbiosis</article-title>. <source>Genome Res</source>. <volume>22</volume>, <fpage>2467</fpage>&#x02013;<lpage>2477</lpage>. <pub-id pub-id-type="doi">10.1101/gr.138420.112</pub-id><pub-id pub-id-type="pmid">22919073</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deng</surname> <given-names>Z. L.</given-names></name> <name><surname>M&#x000FC;nch</surname> <given-names>P. C.</given-names></name> <name><surname>Mreches</surname> <given-names>R.</given-names></name> <name><surname>McHardy</surname> <given-names>A. C.</given-names></name></person-group> (<year>2022</year>). <article-title>Rapid and accurate identification of ribosomal RNA sequences via deep learning</article-title>. <source>Nucl. Acids Res</source>. <volume>50</volume>:<fpage>e60</fpage>. <pub-id pub-id-type="doi">10.1093/nar/gkac112</pub-id><pub-id pub-id-type="pmid">35188571</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dufour</surname> <given-names>Y. S.</given-names></name> <name><surname>Wesenberg</surname> <given-names>G. E.</given-names></name> <name><surname>Tritt</surname> <given-names>A. J.</given-names></name> <name><surname>Glasner</surname> <given-names>J. D.</given-names></name> <name><surname>Perna</surname> <given-names>N. T.</given-names></name> <name><surname>Mitchell</surname> <given-names>J. C.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>chipD: a web tool to design oligonucleotide probes for high-density tiling arrays</article-title>. <source>Nucl. Acids Res</source>. <volume>38</volume>, <fpage>W321</fpage>&#x02013;<lpage>W325</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkq517</pub-id><pub-id pub-id-type="pmid">20529880</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Foster</surname> <given-names>J.</given-names></name> <name><surname>Ganatra</surname> <given-names>M.</given-names></name> <name><surname>Kamal</surname> <given-names>I.</given-names></name> <name><surname>Ware</surname> <given-names>J.</given-names></name> <name><surname>Makarova</surname> <given-names>K.</given-names></name> <name><surname>Ivanova</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2005</year>). <article-title>The <italic>Wolbachia</italic> genome of <italic>Brugia malayi</italic>: endosymbiont evolution within a human pathogenic nematode</article-title>. <source>PLoS Biol</source>. <volume>3</volume>:<fpage>e121</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pbio.0030121</pub-id><pub-id pub-id-type="pmid">15780005</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Foster</surname> <given-names>J. M.</given-names></name> <name><surname>Grote</surname> <given-names>A.</given-names></name> <name><surname>Mattick</surname> <given-names>J.</given-names></name> <name><surname>Tracey</surname> <given-names>A.</given-names></name> <name><surname>Tsai</surname> <given-names>Y. C.</given-names></name> <name><surname>Chung</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Sex chromosome evolution in parasitic nematodes of humans</article-title>. <source>Nat. Commun</source>. <volume>11</volume>:<fpage>1964</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-020-15654-6</pub-id><pub-id pub-id-type="pmid">32327641</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grote</surname> <given-names>A.</given-names></name> <name><surname>Voronin</surname> <given-names>D.</given-names></name> <name><surname>Ding</surname> <given-names>T.</given-names></name> <name><surname>Twaddle</surname> <given-names>A.</given-names></name> <name><surname>Unnasch</surname> <given-names>T. R.</given-names></name> <name><surname>Lustigman</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Defining <italic>Brugia malayi</italic> and <italic>Wolbachia</italic> symbiosis by stage-specific dual RNA-seq</article-title>. <source>PLoS Negl. Trop. Dis</source>. <volume>11</volume>:<fpage>e0005357</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pntd.0005357</pub-id><pub-id pub-id-type="pmid">28358880</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>He</surname> <given-names>S.</given-names></name> <name><surname>Wurtzel</surname> <given-names>O.</given-names></name> <name><surname>Singh</surname> <given-names>K.</given-names></name> <name><surname>Froula</surname> <given-names>J. L.</given-names></name> <name><surname>Yilmaz</surname> <given-names>S.</given-names></name> <name><surname>Tringe</surname> <given-names>S. G.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>Validation of two ribosomal RNA removal methods for microbial metatranscriptomics</article-title>. <source>Nat. Methods</source> <volume>7</volume>, <fpage>807</fpage>&#x02013;<lpage>812</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.1507</pub-id><pub-id pub-id-type="pmid">20852648</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hoerauf</surname> <given-names>A.</given-names></name> <name><surname>Specht</surname> <given-names>S.</given-names></name> <name><surname>B&#x000FC;ttner</surname> <given-names>M.</given-names></name> <name><surname>Pfarr</surname> <given-names>K.</given-names></name> <name><surname>Mand</surname> <given-names>S.</given-names></name> <name><surname>Fimmers</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title><italic>Wolbachia</italic> endobacteria depletion by doxycycline as antifilarial therapy has macrofilaricidal activity in onchocerciasis: a randomized placebo-controlled study</article-title>. <source>Med. Microbiol. Immunol</source>. <volume>197</volume>, <fpage>295</fpage>&#x02013;<lpage>311</lpage>. <pub-id pub-id-type="doi">10.1007/s00430-007-0062-1</pub-id><pub-id pub-id-type="pmid">17999080</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hrdlickova</surname> <given-names>R.</given-names></name> <name><surname>Toloue</surname> <given-names>M.</given-names></name> <name><surname>Tian</surname> <given-names>B.</given-names></name></person-group> (<year>2017</year>). <article-title>RNA-seq methods for transcriptome analysis</article-title>. <source>WIREs RNA</source> <volume>8</volume>:<fpage>e1364</fpage>. <pub-id pub-id-type="doi">10.1002/wrna.1364</pub-id><pub-id pub-id-type="pmid">27198714</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johnston</surname> <given-names>K. L.</given-names></name> <name><surname>Ford</surname> <given-names>L.</given-names></name> <name><surname>Umareddy</surname> <given-names>I.</given-names></name> <name><surname>Townson</surname> <given-names>S.</given-names></name> <name><surname>Specht</surname> <given-names>S.</given-names></name> <name><surname>Pfarr</surname> <given-names>K.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Repurposing of approved drugs from the human pharmacopoeia to target <italic>Wolbachia</italic> endosymbionts of onchocerciasis and lymphatic filariasis</article-title>. <source>Int. J. Parasitol. Drugs Drug Resist</source>. <volume>4</volume>, <fpage>278</fpage>&#x02013;<lpage>286</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijpddr.2014.09.001</pub-id><pub-id pub-id-type="pmid">25516838</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johnston</surname> <given-names>K. L.</given-names></name> <name><surname>Hong</surname> <given-names>W. D.</given-names></name> <name><surname>Turner</surname> <given-names>J. D.</given-names></name> <name><surname>O&#x00027;Neill</surname> <given-names>P. M.</given-names></name> <name><surname>Ward</surname> <given-names>S. A.</given-names></name> <name><surname>Taylor</surname> <given-names>M. J.</given-names></name></person-group> (<year>2021</year>). <article-title>Anti-<italic>Wolbachia</italic> drugs for filariasis</article-title>. <source>Trends Parasitol</source>. <volume>37</volume>, <fpage>1068</fpage>&#x02013;<lpage>1081</lpage>. <pub-id pub-id-type="doi">10.1016/j.pt.2021.06.004</pub-id><pub-id pub-id-type="pmid">34229954</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kang</surname> <given-names>Y.</given-names></name> <name><surname>Norris</surname> <given-names>M. H.</given-names></name> <name><surname>Zarzycki-Siek</surname> <given-names>J.</given-names></name> <name><surname>Nierman</surname> <given-names>W. C.</given-names></name> <name><surname>Donachie</surname> <given-names>S. P.</given-names></name> <name><surname>Hoang</surname> <given-names>T. T.</given-names></name></person-group> (<year>2011</year>). <article-title>Transcript amplification from single bacterium for transcriptome analysis</article-title>. <source>Genome Res</source>. <volume>21</volume>, <fpage>925</fpage>&#x02013;<lpage>935</lpage>. <pub-id pub-id-type="doi">10.1101/gr.116103.110</pub-id><pub-id pub-id-type="pmid">21536723</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>D.</given-names></name> <name><surname>Paggi</surname> <given-names>J. M.</given-names></name> <name><surname>Park</surname> <given-names>C.</given-names></name> <name><surname>Bennett</surname> <given-names>C.</given-names></name> <name><surname>Salzberg</surname> <given-names>S. L.</given-names></name></person-group> (<year>2019</year>). <article-title>Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype</article-title>. <source>Nat. Biotechnol</source>. <volume>37</volume>, <fpage>907</fpage>&#x02013;<lpage>915</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-019-0201-4</pub-id><pub-id pub-id-type="pmid">31375807</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>N.</given-names></name> <name><surname>Creasy</surname> <given-names>T.</given-names></name> <name><surname>Sun</surname> <given-names>Y.</given-names></name> <name><surname>Flowers</surname> <given-names>M.</given-names></name> <name><surname>Tallon</surname> <given-names>L. J.</given-names></name> <name><surname>Dunning Hotopp</surname> <given-names>J. C.</given-names></name></person-group> (<year>2012</year>). <article-title>Efficient subtraction of insect rRNA prior to transcriptome analysis of <italic>Wolbachia</italic>-<italic>Drosophila</italic> lateral gene transfer</article-title>. <source>BMC Res. Not</source>. <volume>5</volume>:<fpage>230</fpage>. <pub-id pub-id-type="doi">10.1186/1756-0500-5-230</pub-id><pub-id pub-id-type="pmid">22583543</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>N.</given-names></name> <name><surname>Lin</surname> <given-names>M.</given-names></name> <name><surname>Zhao</surname> <given-names>X.</given-names></name> <name><surname>Ott</surname> <given-names>S.</given-names></name> <name><surname>Santana-Cruz</surname> <given-names>I.</given-names></name> <name><surname>Daugherty</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Efficient enrichment of bacterial mRNA from host-bacteria total RNA samples</article-title>. <source>Sci. Rep</source>. <volume>6</volume>:<fpage>34850</fpage>. <pub-id pub-id-type="doi">10.1038/srep34850</pub-id><pub-id pub-id-type="pmid">27713560</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>H.</given-names></name> <name><surname>Handsaker</surname> <given-names>B.</given-names></name> <name><surname>Wysoker</surname> <given-names>A.</given-names></name> <name><surname>Fennell</surname> <given-names>T.</given-names></name> <name><surname>Ruan</surname> <given-names>J.</given-names></name> <name><surname>Homer</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>The sequence alignment/map format and SAMtools</article-title>. <source>Bioinformatics</source> <volume>25</volume>, <fpage>2078</fpage>&#x02013;<lpage>2079</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btp352</pub-id><pub-id pub-id-type="pmid">19505943</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Carlow</surname> <given-names>C. K. S.</given-names></name></person-group> (<year>2012</year>). <article-title>Characterization of transcription factors that regulate the type IV secretion system and riboflavin biosynthesis in <italic>Wolbachia</italic> of <italic>Brugia malayi</italic></article-title>. <source>PLoS ONE</source> <volume>7</volume>:<fpage>e51597</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0051597</pub-id><pub-id pub-id-type="pmid">23251587</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liao</surname> <given-names>Y.</given-names></name> <name><surname>Smyth</surname> <given-names>G. K.</given-names></name> <name><surname>Shi</surname> <given-names>W.</given-names></name></person-group> (<year>2014</year>). <article-title>featureCounts: an efficient general purpose program for assigning sequence reads to genomic features</article-title>. <source>Bioinformatics</source> <volume>30</volume>, <fpage>923</fpage>&#x02013;<lpage>930</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btt656</pub-id><pub-id pub-id-type="pmid">24227677</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luck</surname> <given-names>A. N.</given-names></name> <name><surname>Yuan</surname> <given-names>X.</given-names></name> <name><surname>Voronin</surname> <given-names>D.</given-names></name> <name><surname>Slatko</surname> <given-names>B. E.</given-names></name> <name><surname>Hamza</surname> <given-names>I.</given-names></name> <name><surname>Foster</surname> <given-names>J. M.</given-names></name></person-group> (<year>2016</year>). <article-title>Heme acquisition in the parasitic filarial nematode <italic>Brugia malayi</italic></article-title>. <source>FASEB J</source>. <volume>30</volume>, <fpage>3501</fpage>&#x02013;<lpage>3514</lpage>. <pub-id pub-id-type="doi">10.1096/fj.201600603R</pub-id><pub-id pub-id-type="pmid">27363426</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Mackenzie</surname> <given-names>C. D.</given-names></name> <name><surname>Malecela</surname> <given-names>M.</given-names></name> <name><surname>Mueller</surname> <given-names>I.</given-names></name> <name><surname>Homeida</surname> <given-names>M. A.</given-names></name></person-group> (<year>2002</year>). <article-title>&#x0201C;Approaches to the control and elimination of the clinically important filarial diseases,&#x0201D;</article-title> in <source>The Filaria</source> (<publisher-loc>Boston, MA</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>155</fpage>&#x02013;<lpage>165</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://link.springer.com/chapter/10.1007/0-306-47661-4_13">https://link.springer.com/chapter/10.1007/0-306-47661-4_13</ext-link> (accessed October 24, 2023).</citation>
</ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mand</surname> <given-names>S.</given-names></name> <name><surname>Pfarr</surname> <given-names>K.</given-names></name> <name><surname>Sahoo</surname> <given-names>P. K.</given-names></name> <name><surname>Satapathy</surname> <given-names>A. K.</given-names></name> <name><surname>Specht</surname> <given-names>S.</given-names></name> <name><surname>Klarmann</surname> <given-names>U.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>Macrofilaricidal activity and amelioration of lymphatic pathology in bancroftian filariasis after 3 weeks of doxycycline followed by single-dose diethylcarbamazine</article-title>. <source>Am. J. Trop. Med. Hyg</source>. <volume>81</volume>, <fpage>702</fpage>&#x02013;<lpage>711</lpage>. <pub-id pub-id-type="doi">10.4269/ajtmh.2009.09-0155</pub-id><pub-id pub-id-type="pmid">19815891</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>M.</given-names></name></person-group> (<year>2011</year>). <article-title>Cutadapt removes adapter sequences from high-throughput sequencing reads</article-title>. <source>EMBnet.j</source> <volume>17</volume>, <fpage>10</fpage>&#x02013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.14806/ej.17.1.200</pub-id><pub-id pub-id-type="pmid">28715235</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mclaren</surname> <given-names>D. J.</given-names></name> <name><surname>Worms</surname> <given-names>M. J.</given-names></name> <name><surname>Laurence</surname> <given-names>B. R.</given-names></name> <name><surname>Simpson</surname> <given-names>M. G.</given-names></name></person-group> (<year>1975</year>). <article-title>Micro-organisms in filarial larvae (Nematoda)</article-title>. <source>Trans. R Soc. Trop. Med. Hyg</source>. <volume>69</volume>, <fpage>509</fpage>&#x02013;<lpage>514</lpage>. <pub-id pub-id-type="doi">10.1016/0035-9203(75)90110-8</pub-id><pub-id pub-id-type="pmid">1228988</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Medeiros</surname> <given-names>Z. M.</given-names></name> <name><surname>Vieira</surname> <given-names>A. V. B.</given-names></name> <name><surname>Xavier</surname> <given-names>A. T.</given-names></name> <name><surname>Bezerra</surname> <given-names>G. S. N.</given-names></name> <name><surname>Lopes</surname> <given-names>M. d. F. C.</given-names></name> <name><surname>Bonfim</surname> <given-names>C. V.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Lymphatic filariasis: a systematic review on morbidity and its repercussions in countries in the Americas</article-title>. <source>Int. J. Environ. Res. Publ. Health</source> <volume>19</volume>:<fpage>316</fpage>. <pub-id pub-id-type="doi">10.3390/ijerph19010316</pub-id><pub-id pub-id-type="pmid">35010576</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mediannikov</surname> <given-names>O.</given-names></name> <name><surname>Ranque</surname> <given-names>S.</given-names></name></person-group> (<year>2018</year>). <article-title>Mansonellosis, the most neglected human filariasis</article-title>. <source>N. Microb. N. Infect</source>. <volume>26</volume>, <fpage>S19</fpage>&#x02013;<lpage>S22</lpage>. <pub-id pub-id-type="doi">10.1016/j.nmni.2018.08.016</pub-id><pub-id pub-id-type="pmid">30402239</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Michalski</surname> <given-names>M. L.</given-names></name> <name><surname>Griffiths</surname> <given-names>K. G.</given-names></name> <name><surname>Williams</surname> <given-names>S. A.</given-names></name> <name><surname>Kaplan</surname> <given-names>R. M.</given-names></name> <name><surname>Moorhead</surname> <given-names>A. R.</given-names></name></person-group> (<year>2011</year>). <article-title>The NIH-NIAID filariasis research reagent resource center</article-title>. <source>PLoS Negl. Trop. Dis</source>. <volume>5</volume>:<fpage>e1261</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pntd.0001261</pub-id><pub-id pub-id-type="pmid">22140585</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Molyneux</surname> <given-names>D. H.</given-names></name> <name><surname>Bradley</surname> <given-names>M.</given-names></name> <name><surname>Hoerauf</surname> <given-names>A.</given-names></name> <name><surname>Kyelem</surname> <given-names>D.</given-names></name> <name><surname>Taylor</surname> <given-names>M. J.</given-names></name></person-group> (<year>2003</year>). <article-title>Mass drug treatment for lymphatic filariasis and onchocerciasis</article-title>. <source>Trends Parasitol</source>. <volume>19</volume>, <fpage>516</fpage>&#x02013;<lpage>522</lpage>. <pub-id pub-id-type="doi">10.1016/j.pt.2003.09.004</pub-id><pub-id pub-id-type="pmid">14580963</pub-id></citation></ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Parkinson</surname> <given-names>J.</given-names></name> <name><surname>Mitreva</surname> <given-names>M.</given-names></name> <name><surname>Whitton</surname> <given-names>C.</given-names></name> <name><surname>Thomson</surname> <given-names>M.</given-names></name> <name><surname>Daub</surname> <given-names>J.</given-names></name> <name><surname>Martin</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2004</year>). <article-title>A transcriptomic analysis of the phylum Nematoda</article-title>. <source>Nat. Genet</source>. <volume>36</volume>, <fpage>1259</fpage>&#x02013;<lpage>1267</lpage>. <pub-id pub-id-type="doi">10.1038/ng1472</pub-id><pub-id pub-id-type="pmid">15543149</pub-id></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Prezza</surname> <given-names>G.</given-names></name> <name><surname>Heckel</surname> <given-names>T.</given-names></name> <name><surname>Dietrich</surname> <given-names>S.</given-names></name> <name><surname>Homberger</surname> <given-names>C.</given-names></name> <name><surname>Westermann</surname> <given-names>A. J.</given-names></name> <name><surname>Vogel</surname> <given-names>J.</given-names></name></person-group> (<year>2020</year>). <article-title>Improved bacterial RNA-seq by Cas9-based depletion of ribosomal RNA reads</article-title>. <source>RNA</source> <volume>26</volume>, <fpage>1069</fpage>&#x02013;<lpage>1078</lpage>. <pub-id pub-id-type="doi">10.1261/rna.075945.120</pub-id><pub-id pub-id-type="pmid">32345633</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quek</surname> <given-names>S.</given-names></name> <name><surname>Cook</surname> <given-names>D. A. N.</given-names></name> <name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Marriott</surname> <given-names>A. E.</given-names></name> <name><surname>Steven</surname> <given-names>A.</given-names></name> <name><surname>Johnston</surname> <given-names>K. L.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title><italic>Wolbachia</italic> depletion blocks transmission of lymphatic filariasis by preventing chitinase-dependent parasite exsheathment</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A</source>. <volume>119</volume>:<fpage>e2120003119</fpage>. <pub-id pub-id-type="doi">10.1073/pnas.2120003119</pub-id><pub-id pub-id-type="pmid">35377795</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quinlan</surname> <given-names>A. R.</given-names></name> <name><surname>Hall</surname> <given-names>I. M.</given-names></name></person-group> (<year>2010</year>). <article-title>BEDTools: a flexible suite of utilities for comparing genomic features</article-title>. <source>Bioinformatics</source> <volume>26</volume>, <fpage>841</fpage>&#x02013;<lpage>842</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btq033</pub-id><pub-id pub-id-type="pmid">20110278</pub-id></citation></ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ram&#x000ED;rez</surname> <given-names>F.</given-names></name> <name><surname>D&#x000FC;ndar</surname> <given-names>F.</given-names></name> <name><surname>Diehl</surname> <given-names>S.</given-names></name> <name><surname>Gr&#x000FC;ning</surname> <given-names>B. A.</given-names></name> <name><surname>Manke</surname> <given-names>T.</given-names></name></person-group> (<year>2014</year>). <article-title>deepTools: a flexible platform for exploring deep-sequencing data</article-title>. <source>Nucl. Acids Res</source>. <volume>42</volume>, <fpage>W187</fpage>&#x02013;<lpage>W191</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gku365</pub-id><pub-id pub-id-type="pmid">24799436</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Richards</surname> <given-names>F. O.</given-names></name> <name><surname>Boatin</surname> <given-names>B.</given-names></name> <name><surname>Sauerbrey</surname> <given-names>M.</given-names></name> <name><surname>S&#x000E9;k&#x000E9;t&#x000E9;li</surname> <given-names>A.</given-names></name></person-group> (<year>2001</year>). <article-title>Control of onchocerciasis today: status and challenges</article-title>. <source>Trends Parasitol</source>. <volume>17</volume>, <fpage>558</fpage>&#x02013;<lpage>563</lpage>. <pub-id pub-id-type="doi">10.1016/S1471-4922(01)02112-2</pub-id><pub-id pub-id-type="pmid">11756018</pub-id></citation></ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rognes</surname> <given-names>T.</given-names></name> <name><surname>Flouri</surname> <given-names>T.</given-names></name> <name><surname>Nichols</surname> <given-names>B.</given-names></name> <name><surname>Quince</surname> <given-names>C.</given-names></name> <name><surname>Mah&#x000E9;</surname> <given-names>F</given-names></name></person-group>. (<year>2016</year>). <article-title>VSEARCH: a versatile open source tool for metagenomics</article-title>. <source>PeerJ</source> <volume>4</volume>:<fpage>e2584</fpage>. <pub-id pub-id-type="doi">10.7717/peerj.2584</pub-id><pub-id pub-id-type="pmid">27781170</pub-id></citation></ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schuierer</surname> <given-names>S.</given-names></name> <name><surname>Carbone</surname> <given-names>W.</given-names></name> <name><surname>Knehr</surname> <given-names>J.</given-names></name> <name><surname>Petitjean</surname> <given-names>V.</given-names></name> <name><surname>Fernandez</surname> <given-names>A.</given-names></name> <name><surname>Sultan</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>A comprehensive assessment of RNA-seq protocols for degraded and low-quantity samples</article-title>. <source>BMC Genomics</source> <volume>18</volume>:<fpage>442</fpage>. <pub-id pub-id-type="doi">10.1186/s12864-017-3827-y</pub-id><pub-id pub-id-type="pmid">28583074</pub-id></citation></ref>
<ref id="B53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sim&#x000F3;n</surname> <given-names>F.</given-names></name> <name><surname>Siles-Lucas</surname> <given-names>M.</given-names></name> <name><surname>Morch&#x000F3;n</surname> <given-names>R.</given-names></name> <name><surname>Gonz&#x000E1;lez-Miguel</surname> <given-names>J.</given-names></name> <name><surname>Mellado</surname> <given-names>I.</given-names></name> <name><surname>Carret&#x000F3;n</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Human and animal dirofilariasis: the emergence of a zoonotic mosaic</article-title>. <source>Clin. Microbiol. Rev</source>. <volume>25</volume>, <fpage>507</fpage>&#x02013;<lpage>544</lpage>. <pub-id pub-id-type="doi">10.1128/CMR.00012-12</pub-id><pub-id pub-id-type="pmid">22763636</pub-id></citation></ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sironi</surname> <given-names>M.</given-names></name> <name><surname>Bandi</surname> <given-names>C.</given-names></name> <name><surname>Sacchi</surname> <given-names>L.</given-names></name> <name><surname>Sacco</surname> <given-names>B. D.</given-names></name> <name><surname>Damiani</surname> <given-names>G.</given-names></name> <name><surname>Genchi</surname> <given-names>C.</given-names></name></person-group> (<year>1995</year>). <article-title>Molecular evidence for a close relative of the arthropod endosymbiont <italic>Wolbachia</italic> in a filarial worm</article-title>. <source>Mol. Biochem. Parasitol</source>. <volume>74</volume>, <fpage>223</fpage>&#x02013;<lpage>227</lpage>. <pub-id pub-id-type="doi">10.1016/0166-6851(95)02494-8</pub-id><pub-id pub-id-type="pmid">8719164</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Small</surname> <given-names>S. T.</given-names></name> <name><surname>Tisch</surname> <given-names>D. J.</given-names></name> <name><surname>Zimmerman</surname> <given-names>P. A.</given-names></name></person-group> (<year>2014</year>). <article-title>Molecular epidemiology, phylogeny and evolution of the filarial nematode <italic>Wuchereria bancrofti</italic></article-title>. <source>Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis</source>. <volume>2014</volume>, <fpage>33</fpage>&#x02013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.meegid.2014.08.018</pub-id><pub-id pub-id-type="pmid">25176600</pub-id></citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Specht</surname> <given-names>S.</given-names></name> <name><surname>Mand</surname> <given-names>S.</given-names></name> <name><surname>Marfo-Debrekyei</surname> <given-names>Y.</given-names></name> <name><surname>Debrah</surname> <given-names>A. Y.</given-names></name> <name><surname>Konadu</surname> <given-names>P.</given-names></name> <name><surname>Adjei</surname> <given-names>O.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Efficacy of 2- and 4-week rifampicin treatment on the <italic>Wolbachia</italic> of <italic>Onchocerca volvulus</italic></article-title>. <source>Parasitol. Res</source>. <volume>103</volume>, <fpage>1303</fpage>&#x02013;<lpage>1309</lpage>. <pub-id pub-id-type="doi">10.1007/s00436-008-1133-y</pub-id><pub-id pub-id-type="pmid">18679718</pub-id></citation></ref>
<ref id="B57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sultan</surname> <given-names>M.</given-names></name> <name><surname>Amstislavskiy</surname> <given-names>V.</given-names></name> <name><surname>Risch</surname> <given-names>T.</given-names></name> <name><surname>Schuette</surname> <given-names>M.</given-names></name> <name><surname>D&#x000F6;kel</surname> <given-names>S.</given-names></name> <name><surname>Ralser</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Influence of RNA extraction methods and library selection schemes on RNA-seq data</article-title>. <source>BMC Genomics</source> <volume>15</volume>:<fpage>675</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-675</pub-id><pub-id pub-id-type="pmid">25113896</pub-id></citation></ref>
<ref id="B58">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Supali</surname> <given-names>T.</given-names></name> <name><surname>Djuardi</surname> <given-names>Y.</given-names></name> <name><surname>Pfarr</surname> <given-names>K. M.</given-names></name> <name><surname>Wibowo</surname> <given-names>H.</given-names></name> <name><surname>Taylor</surname> <given-names>M. J.</given-names></name> <name><surname>Hoerauf</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Doxycycline treatment of <italic>Brugia malayi</italic>-infected persons reduces microfilaremia and adverse reactions after diethylcarbamazine and albendazole treatment</article-title>. <source>Clin. Infect. Dis</source>. <volume>46</volume>, <fpage>1385</fpage>&#x02013;<lpage>1393</lpage>. <pub-id pub-id-type="doi">10.1086/586753</pub-id><pub-id pub-id-type="pmid">18419441</pub-id></citation></ref>
<ref id="B59">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taylor</surname> <given-names>M. J.</given-names></name> <name><surname>Bandi</surname> <given-names>C.</given-names></name> <name><surname>Hoerauf</surname> <given-names>A.</given-names></name></person-group> (<year>2005</year>). <article-title><italic>Wolbachia</italic> bacterial endosymbionts of filarial nematodes</article-title>. <source>Adv. Parasitol</source>. <volume>60</volume>, <fpage>245</fpage>&#x02013;<lpage>284</lpage>. <pub-id pub-id-type="doi">10.1016/S0065-308X(05)60004-8</pub-id><pub-id pub-id-type="pmid">16230105</pub-id></citation></ref>
<ref id="B60">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taylor</surname> <given-names>M. J.</given-names></name> <name><surname>Hoerauf</surname> <given-names>A.</given-names></name> <name><surname>Bockarie</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Lymphatic filariasis and onchocerciasis</article-title>. <source>Lancet</source> <volume>376</volume>, <fpage>1175</fpage>&#x02013;<lpage>1185</lpage>. <pub-id pub-id-type="doi">10.1016/S0140-6736(10)60586-7</pub-id><pub-id pub-id-type="pmid">20739055</pub-id></citation></ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taylor</surname> <given-names>M. J.</given-names></name> <name><surname>von Geldern</surname> <given-names>T. W.</given-names></name> <name><surname>Ford</surname> <given-names>L.</given-names></name> <name><surname>H&#x000FC;bner</surname> <given-names>M. P.</given-names></name> <name><surname>Marsh</surname> <given-names>K.</given-names></name> <name><surname>Johnston</surname> <given-names>K. L.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Preclinical development of an oral anti-<italic>Wolbachia</italic> macrolide drug for the treatment of lymphatic filariasis and onchocerciasis</article-title>. <source>Sci. Transl. Med</source>. <volume>11</volume>:<fpage>eaau2086</fpage>. <pub-id pub-id-type="doi">10.1126/scitranslmed.aau2086</pub-id><pub-id pub-id-type="pmid">30867321</pub-id></citation></ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thorvaldsd&#x000F3;ttir</surname> <given-names>H.</given-names></name> <name><surname>Robinson</surname> <given-names>J. T.</given-names></name> <name><surname>Mesirov</surname> <given-names>J. P.</given-names></name></person-group> (<year>2013</year>). <article-title>Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration</article-title>. <source>Brief Bioinform</source>. <volume>14</volume>, <fpage>178</fpage>&#x02013;<lpage>192</lpage>. <pub-id pub-id-type="doi">10.1093/bib/bbs017</pub-id><pub-id pub-id-type="pmid">22517427</pub-id></citation></ref>
<ref id="B63">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Valero-Mora</surname> <given-names>P. M.</given-names></name></person-group> (<year>2010</year>). <article-title>ggplot2: elegant graphics for data analysis</article-title>. <source>J. Stat. Softw</source>. <volume>35</volume>, <fpage>1</fpage>&#x02013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.18637/jss.v035.b01</pub-id></citation>
</ref>
<ref id="B64">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wahl</surname> <given-names>A.</given-names></name> <name><surname>Huptas</surname> <given-names>C.</given-names></name> <name><surname>Neuhaus</surname> <given-names>K.</given-names></name></person-group> (<year>2022</year>). <article-title>Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing</article-title>. <source>Sci. Rep</source>. <volume>12</volume>:<fpage>5765</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-022-09710-y</pub-id><pub-id pub-id-type="pmid">35388078</pub-id></citation></ref>
<ref id="B65">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>W.</given-names></name></person-group> (<year>2012</year>). <article-title>RSeQC: quality control of RNA-seq experiments</article-title>. <source>Bioinformatics</source> <volume>28</volume>, <fpage>2184</fpage>&#x02013;<lpage>2185</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/bts356</pub-id><pub-id pub-id-type="pmid">22743226</pub-id></citation></ref>
<ref id="B66">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Gerstein</surname> <given-names>M.</given-names></name> <name><surname>Snyder</surname> <given-names>M.</given-names></name></person-group> (<year>2009</year>). <article-title>RNA-Seq: a revolutionary tool for transcriptomics</article-title>. <source>Nat. Rev. Genet</source>. <volume>10</volume>, <fpage>57</fpage>&#x02013;<lpage>63</lpage>. <pub-id pub-id-type="doi">10.1038/nrg2484</pub-id><pub-id pub-id-type="pmid">19015660</pub-id></citation></ref>
<ref id="B67">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wangsanuwat</surname> <given-names>C.</given-names></name> <name><surname>Heom</surname> <given-names>K. A.</given-names></name> <name><surname>Liu</surname> <given-names>E.</given-names></name> <name><surname>O&#x00027;Malley</surname> <given-names>M. A.</given-names></name> <name><surname>Dey</surname> <given-names>S. S.</given-names></name></person-group> (<year>2020</year>). <article-title>Efficient and cost-effective bacterial mRNA sequencing from low input samples through ribosomal RNA depletion</article-title>. <source>BMC Genomics</source> <volume>21</volume>:<fpage>717</fpage>. <pub-id pub-id-type="doi">10.1186/s12864-020-07134-4</pub-id><pub-id pub-id-type="pmid">33066726</pub-id></citation></ref>
<ref id="B68">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wanji</surname> <given-names>S.</given-names></name> <name><surname>Tendongfor</surname> <given-names>N.</given-names></name> <name><surname>Nji</surname> <given-names>T.</given-names></name> <name><surname>Esum</surname> <given-names>M.</given-names></name> <name><surname>Che</surname> <given-names>J. N.</given-names></name> <name><surname>Nkwescheu</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>Community-directed delivery of doxycycline for the treatment of onchocerciasis in areas of co-endemicity with loiasis in Cameroon</article-title>. <source>Parasit. Vect</source>. <volume>2</volume>:<fpage>39</fpage>. <pub-id pub-id-type="doi">10.1186/1756-3305-2-39</pub-id><pub-id pub-id-type="pmid">19712455</pub-id></citation></ref>
<ref id="B69">
<citation citation-type="web"><person-group person-group-type="author"><collab>WHO Expert Committee on Onchocerciasis Control</collab></person-group> (<year>1993</year>). <source>Onchocerciasis and Its Control: Report of a WHO Expert Committee on Onchocerciasis Control</source>. Geneva: World Health Organization. Available online at: <ext-link ext-link-type="uri" xlink:href="https://iris.who.int/handle/10665/37346">https://iris.who.int/handle/10665/37346</ext-link> (accessed October 24, 2023).</citation>
</ref>
<ref id="B70">
<citation citation-type="journal"><person-group person-group-type="author"><collab>World Health Organization</collab></person-group> (<year>2019</year>). <article-title>Global programme to eliminate lymphatic filariasis: progress report, 2018</article-title>. <source>Wkly Epidemiol. Rec</source>. <volume>94</volume>, <fpage>457</fpage>&#x02013;<lpage>470</lpage>.</citation>
</ref>
<ref id="B71">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>B.</given-names></name> <name><surname>Novelli</surname> <given-names>J.</given-names></name> <name><surname>Foster</surname> <given-names>J.</given-names></name> <name><surname>Vaisvila</surname> <given-names>R.</given-names></name> <name><surname>Conway</surname> <given-names>L.</given-names></name> <name><surname>Ingram</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>The heme biosynthetic pathway of the obligate <italic>Wolbachia</italic> endosymbiont of <italic>Brugia malayi</italic> as a potential anti-filarial drug target</article-title>. <source>PLoS Negl. Trop. Dis</source>. <volume>3</volume>:<fpage>e475</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pntd.0000475</pub-id><pub-id pub-id-type="pmid">19597542</pub-id></citation></ref>
<ref id="B72">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiang</surname> <given-names>X.</given-names></name> <name><surname>Poli</surname> <given-names>D.</given-names></name> <name><surname>Degnan</surname> <given-names>B. M.</given-names></name> <name><surname>Degnan</surname> <given-names>S. M.</given-names></name></person-group> (<year>2022</year>). <article-title>Ribosomal RNA-depletion provides an efficient method for successful dual RNA-seq expression profiling of a marine sponge holobiont</article-title>. <source>Mar. Biotechnol</source>. <volume>24</volume>, <fpage>722</fpage>&#x02013;<lpage>732</lpage>. <pub-id pub-id-type="doi">10.1007/s10126-022-10138-8</pub-id><pub-id pub-id-type="pmid">35895230</pub-id></citation></ref>
<ref id="B73">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xie</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>S.</given-names></name> <name><surname>Wu</surname> <given-names>S.</given-names></name> <name><surname>Gao</surname> <given-names>S.</given-names></name> <name><surname>Meng</surname> <given-names>Q.</given-names></name> <name><surname>Wang</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Genome of the giant panda roundworm illuminates its host shift and parasitic adaptation</article-title>. <source>Genom. Proteom. Bioinformat</source>. <volume>20</volume>, <fpage>366</fpage>&#x02013;<lpage>381</lpage>. <pub-id pub-id-type="doi">10.1016/j.gpb.2021.08.002</pub-id><pub-id pub-id-type="pmid">34487863</pub-id></citation></ref>
<ref id="B74">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>S.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Gamini</surname> <given-names>R.</given-names></name> <name><surname>Zhang</surname> <given-names>B.</given-names></name> <name><surname>von Schack</surname> <given-names>D.</given-names></name></person-group> (<year>2018</year>). <article-title>Evaluation of two main RNA-seq approaches for gene quantification in clinical RNA sequencing: polyA&#x0002B; selection versus rRNA depletion</article-title>. <source>Sci. Rep</source>. <volume>8</volume>:<fpage>4781</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-23226-4</pub-id><pub-id pub-id-type="pmid">29556074</pub-id></citation></ref>
<ref id="B75">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>W.</given-names></name> <name><surname>He</surname> <given-names>X.</given-names></name> <name><surname>Hoadley</surname> <given-names>K. A.</given-names></name> <name><surname>Parker</surname> <given-names>J. S.</given-names></name> <name><surname>Hayes</surname> <given-names>D. N.</given-names></name> <name><surname>Perou</surname> <given-names>C. M.</given-names></name></person-group> (<year>2014</year>). <article-title>Comparison of RNA-Seq by poly (A) capture, ribosomal RNA depletion, and DNA microarray for expression profiling</article-title>. <source>BMC Genomics</source> <volume>15</volume>:<fpage>419</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-419</pub-id><pub-id pub-id-type="pmid">24888378</pub-id></citation></ref>
</ref-list>
</back>
</article>