13th International Research Conference General Sir John Kotelawala Defence University Basic and Applied Sciences Sessions Paper ID: 204 Identification of Wolbachia like endosymbiont DNA in Setaria digitata genome and phylogenetic analysis of filarial

MSA Kothalawala1#, N Rashanthy1, TS Mugunamalwaththa1, WAS Darshika1, GLY Lakmali1, RS Dassanayake1, YINS Gunawardene2, NV Chandrasekharan1, Prashanth Suravajhala3, Kasun de Zoysa4

1Department of Chemistry, Faculty of Science, University of Colombo, Sri Lanka. 2Molecular Medicine Unit, Faculty of Medicine, University of Kelaniya, Ragama, Sri Lanka. 3Department of Biotechnology and Bioinformatics, Birla Institute of Scientific Research, Statue Circle, Jaipur 302001, RJ, India 4Department of Communication and Media Technologies, University of Colombo School of Computing, Sri Lanka.

#For correspondence;

Abstract - Setaria digitata is a Wolbachia-free BLAST2GO analysis was able to identify 6055 filarial parasite that causes cerebrospinal annotations and 95 metabolic pathways nematodiasis in non-permissive hosts such as within the S. digitata genome. Based on goats and sheep, leading to substantial FASTA36 and phylogenetic analyses, it can be economic losses in husbandry. concluded that ancestors of S. digitata were Therefore, there arises a considerable need for colonized with Wolbachia in the distant past, the development of new interventions to and suspected gene transfer may have brought disease control and eradication of this filarial Wolbachia DNA into the S. digitata nuclear parasite. Owing to the limited knowledge on S. genome prior to endosymbiont loss. digitata genome, it's host-parasite Keywords - Setaria digitata, Wolbachia, NGS relationship and the potential impact of the Wolbachia endosymbiont in filarial Introduction nematodes, this research was focused on the S. digitata is an ivory-color slimy filarial generation of the draft genome of S. digitata; parasitic-worm with a coiled tapering tail. It is identification of Wolbachia like endosymbiont classified under class , order DNA in S. digitata genome; phylogenetic and family . S. digitata analysis as well as functional annotation and resides in the peritoneal cavity of grazing metabolic pathway analysis of the genome. A hoofed (Shiny et al., 2011; Shin et al., draft genome of 78.8 Mbp size with a GC% of 2002) and cause cerebrospinal nematodiasis, 31.45% was generated for the S. digitata a neuropathological disorder that causes worms collected from the peritoneal cavity of dysfunction of the central nervous system. It slaughtered cattle, using NGS Illumina leads to lumbar paralysis and eventual death platform. FASTA36 sequence similarity of non-permissive domesticated hosts such as analysis was able to identify homologous goats, sheep and horses, which cause sequences of coxA and gatB Wolbachia MLST substantial economic loss in animal genes within the S. digitata genome, while husbandry in Asia and the Far East. However, phylogenetic analysis using Geneious Prime infection of natural hosts such as cattle and Software revealed that S. digitata is more buffalos by the S. digitata cause mild disease closely related to the filarial nematodes with conditions such as mild fibrinous peritonitis Wolbachia endosymbiont than Wolbachia- and not considered as parasitic in their natural free filarial nematodes. Furthermore, hosts. Human infections are also reported

15 13th International Research Conference General Sir John Kotelawala Defence University Basic and Applied Sciences Sessions causing allergic reactions, eye lesions, potential impact of the endosymbiont and to abscesses, enlarged lymph nodes and lung gain insight into the host-parasite inflammation (Taylor and Hoerauf, 1999; relationship. Gunawardene and Dassanayake, 2015; Sundar S. digitata also considered as a model and Souza, 2015) organism for human lymphatic filariasis Most of the filarial nematodes mutually (HLF), due to their close resemblance to associate with Wolbachia and without Wuchereria bancrofti, the primary causative Wolbachia, such filarial nematodes cannot agent of HLF, in morphology, histology and survive (Werren et al., 2008). The reason is antigenic properties (Gunawardene and that the genome of Wolbachia carries the Dassanayake, 2015). Hence, generation of a genes required for the metabolism of heam, draft genome of S. digitata and complete riboflavin, FAD, glutathione and nucleotides functional analysis and metabolic pathway whereas its filarial host does not (Fenn et al., reconstructions of S. digitata will have a huge 2006). However, S. digitata does not harbor impact on the development of novel drugs Wolbachia endosymbiont, but undergo and/or vaccines for human filariasis and other biosynthetic pathways of heam, riboflavin and filarial diseases. nucleotides. The genes of those pathways Methodology should be present in its genome itself (Voronin et al., 2015). There are two possible Adult worm samples of S. digitata were explanations for the independent survival of S. collected from the peritoneal cavity of cattle digitata. It can be hypothesized that either the (Bos taurus) from the western province. S. ancestors of S. digitata were colonized with digitata worms were washed thoroughly in Wolbachia in the distant past and horizontal PBS and preserved in 80% ethanol at −20˚C gene transfer (HGT) may have brought before analysis. Genomic DNA of adult S. Wolbachia DNA into the nuclear genome of S. digitata worms was extracted by DNA Micro digitata prior to endosymbiont loss, or the kit, QIAGEN. Confirmation of extracted S. endosymbiotic relationship between digitata DNA was done by carrying out a series Wolbachia and its filarial host is dispensable. of PCRs using S. digitata ARV1 specific Studies showed that Wolbachia free filarial primers. Loa loa evolved own DNA After performing quality control (QC), the sequences to code heam and riboflavin passed sample was proceeded with the Next biosynthesis pathways. For some pathways, Generation Sequencing (NGS) library they have gained partial gene sequences from construction. Sequencing libraries were Wolbachia, and this indicates the horizontal constructed from the extracted DNA using the gene transfer within Loa loa and Wolbachia at TruSeq™ DNA PCR-Free Kit. Purified libraries some point during the evolution (Desjardins et were loaded into an Illumina HiSeq4000 for al., 2013). Therefore, independent survival of paired-end sequencing. Sequenced data (base S. digitata also can be explained by horizontal call files) were converted to raw FASTQ files gene transfer of Wolbachia gene fragments to using the sequencer soon after the initial the S. digitata genome and/or evolution of sequencing. The reads were filtered before their own DNA sequences to code those assembly such that for a pair of PE reads each pathways (Mcnulty et al., 2010). Therefore, read should have more than 90 % of bases identification of Wolbachia like endosymbiont with base quality greater than or equal to Q20. DNA in S. digitata genome and phylogenetic Generated 150bp reads were analyzed for K- analysis of S. digitata and other filarial mers using JELLYFISH. Once the optimum k- nematodes, are important to understand the

16 13th International Research Conference General Sir John Kotelawala Defence University Basic and Applied Sciences Sessions mer size is identified, a de novo draft assembly 31.77% and a Q20 value of 99.18%. A draft was built using SOAPdenovo2. genome of 78,774,594 bases belongs to a total of 2,075 contigs was generated after the In order to identify Wolbachia like assembly. The percentage of guanine-cytosine endosymbiont DNA in S. digitata genome, base pairs (GC%) in the assembled draft FASTA file containing 2075 S. digitata contigs genome is 31.45%. In a typical random library, was blasted against several Wolbachia specific it is expected to see a roughly normal reference sequences such as Wolbachia distribution of GC content where the central surface protein (WSP), Wolbachia-specific 16S peak corresponds to the overall GC content. An rRNA and Wolbachia MLST genes (coxA, gatB, unusually shaped distribution could indicate a fbpA, ftsZ, and hcpA) using FASTA36 program. contaminated library or some other kinds of a Later, using Wolbachia marker sequences biased subset. Here, the obtained histogram (gatB, coxA, ftsZ, fbpA, hcpA, WSP and 16S showed a normal distribution and therefore, it rRNA), BLASTx was performed against S. can conclude that there is no contamination. digitata, Loa loa, Brugia malayi, W. bancrofti, Five Wolbachia marker gene sequences were Onchocerca volvulus and Wolbachia selected for identification of Wolbachia like endosymbiont of B. malayi (wBm) protein endosymbiont DNA in S. digitata genome sequences available on the NCBI non- because of their presence in the Rickettsiales redundant (NR) protein sequence database. order. These marker genes have a single copy, Resulted best hit for each organism for each a wide spatial distribution and a strong marker gene was recorded, and the aligned stabilizing selection within the Wolbachia region was obtained. Those aligned sequences genome. Homologous sequences were found of each organism for each marker gene were only for coxA and gatB. For coxA, homologous used to generate a phylogenetic tree using sequence was found in the S. digitata Geneious Prime software. Firstly, sequences contig_633 and for gatB, homologous were loaded into the program and alignment sequences were found in the S. digitata was done using MUSCLE algorithm and the contig_915 and contig_724, among which phylogenetic tree was generated using the contig_915 has the higher similarity. Sequence Neighbor-joining (NJ) method. Wolbachia identity was higher than 50% and e-value was endosymbiont of B. malayi (wBm) was used as smaller than 10-4, and therefore it can be the out-group. considered that sequence similarity is Lastly, BLAST2GO BASIC software was used accountable. for downstream annotations and metabolic After performing BLASTx against protein pathway reconstruction. sequences of S. digitata, Loa loa, Brugia malayi, Results and Discussion Wuchereria bancrofti, Onchocerca volvulus After the initial sequencing using the Illumina and wBm available on the NCBI NR protein platform, a total of 97,586,942 reads were database using the Wolbachia marker obtained with a GC percentage of 31.67% and sequences (gatB, coxA, ftsZ, fbpA, hcpA, WSP a Q20 value of 96.09%. Q20 value gives the and 16S rRNA), only coxA gave hits for all six percentage of bases called that have a quality organisms. Therefore, best hit of each score of 20 or above. Phred quality score is a organism for coxA was used and aligned using numerical value that expresses the accuracy of MUSCLE algorithm to create the phylogenetic each nucleotide. Higher the Q number will be tree (Figure 1). Neighbor-joining (NJ) method higher the accuracy. After the quality control was used here because it is ideal for phylogeny and pre-processing, a total of 79,292,174 construction from the sequence data, it is reads were obtained with a GC percentage of rapid, and it does not assume an equal rate of

17 13th International Research Conference General Sir John Kotelawala Defence University Basic and Applied Sciences Sessions evolution amongst all lineages. The resulted Remaining 1435 BLAST hit generated contigs tree is simply an estimate and is unlikely to were mapped to retrieve GO terms. Out of represent the true evolutionary tree of these 1435 mapped contigs, only 1280 contigs organisms (Figure 2). Resulted tree was a generated GO annotations while 155 contigs rooted tree with 6 taxa and 11 nodes. The were only mapped. nodes are where lineages diverge, representing a speciation event from a common ancestor. The root node represents the most recent common ancestor of all of the taxa represented on the tree. Wolbachia was used as the outgroup because it has Figure 1. Six coxA BLAST hits after alignment using contrasting characteristics relative to the Geneious Software. Source: Geneious Software other included taxa. Wolbachia is a bacterium and all other taxa are nematodes. Several clades can be seen in the resulted phylogenetic tree. W. bancrofti and B. malayi forms a clade while O. volvulus and S. digitata form a separate clade. These four taxa all together form a clade that is related to the L. loa, a Figure 2. Phylogenetic tree constructed using Neighbour Wolbachia-free filarial nematode. joining method for coxA sequence data. Source: Geneious Software According to the evolutionary distance between taxa based on the sum of the branch Majority of the BLAST hits were belong to length (patristic distance), S. digitata is closely filarial nematodes with Wolbachia, like Brugia related to the O. volvulus, a filarial nematode malayi (4998 hits), Wuchereria bancrofti with Wolbachia endosymbiont while most (3311 hits), Brugia pahangi (1940 hits) and distance taxon is Wolbachia. It was expected to Onchocerca ochengi (1673 hits). In here 6055 be closely related to the L. loa since it also a number of total GO annotations have been Wolbachia free nematode just like S. digitata. done. Overall, 1280 of 1435 (89.1%) mapped Based on the generated phylogenetic tree, it BLAST hits were annotated by at least one of can be suggested that Wolbachia genes have the three categories of the GO function been laterally transferred into S. digitata classification. The 1280 mapped contigs could genome since S. digitata is more closely be annotated to 6055 GO terms, among which related to the filarial nematodes with 2297 (grouped in 9 subcategories), 2206 Wolbachia endosymbiont than Wolbachia- (grouped in 8 subcategories) and 1552 free filarial nematode, and since S. digitata and (grouped in 5 subcategories) GO terms could all other Wolbachia containing nematodes be grouped to the biological process category, form a clade separately from Wolbachia-free molecular function category and the cellular L.loa. component category, respectively. Out of 2075 total contigs put through the As an alternative method of categorizing BLAST2GO pipeline, 530 contigs did not contigs/sequences by biochemical function, generate a BLAST hit. Here BLASTx was sequences were assigned to biological performed against the NCBI non-redundant pathways using the KEGG database. Total of database with a cut-off of 1e-3. Out of 1545 246 contigs had been clustered into 95 BLAST hit generated contigs, 110 contigs pathways, in which the most over-represented generated a BLAST hit without further pathways are Biosynthesis of antibiotics downstream Gene Ontology (GO) annotation. (22genes), Phosphatidylinositol signalling

18 13th International Research Conference General Sir John Kotelawala Defence University Basic and Applied Sciences Sessions system (11 genes) and Purine metabolism Gunawardene, Y. I. N. S. and Dassanayake, R. (2015) (10genes). ‘Setaria digitata in advancing our Knowledge of human lymphatic filariasis Setaria digitata in Conclusion advancing our knowledge of human lymphatic According to the BLAST results of the filariasis’, (May). BLAST2GO functional analysis, the majority of Mcnulty, S. N. et al. (2010) ‘Endosymbiont DNA in the S. digitata sequences have a higher Endobacteria-Free Filarial Nematodes Indicates sequence similarity to Wolbachia containing Ancient Horizontal Genetic Transfer’, 5(6). filarial nematodes like Brugia malayi, Shiny, C. et al. (2011) ‘Existence of endosymbiont Wuchereria bancrofti, Brugia pahangi and Wolbachia in Setaria digitata’, 3(2), pp. 108–121. Onchocerca ochengi than Wolbachia free Shin, S. et al. (2002) ‘Ocular Infection of Cattle with filarial nematodes like Loa loa and Onchocerca Setaria digitata’ pp 3- 6 flexuosa. Based on this result it can be concluded that S. digitata is more closely Sundar, S. T. B. and Souza, P. E. D. (2015) ‘Morphological characterization of Setaria worms related to Wolbachia containing nematodes collected from cattle’, Journal of Parasitic Diseases. than Wolbachia free nematodes. Phylogenetic Springer India, 39(3), pp. 572–576. analysis also revealed that S. digitata is more closely related to the Wolbachia containing Taylor, M. J. and Hoerauf, A. (1999) ‘Wolbachia Bacteria of Filarial Nematodes’, 4(1987), pp. 437– nematodes. 442. According to the FASTA36 and BLASTx Voronin, D. et al. (2015) ‘Veterinary Parasitology sequence similarity analysis, partial Absence of Wolbachia endobacteria in Sri Lankan sequences of Wolbachia marker genes (coxA isolates of the nematode parasite of animals Setaria and gatB) were found within the S. digitata digitata’, Veterinary Parasitology. Elsevier B.V., genome. This result provided the 207(3–4), pp. 350–354. bioinformatics evidence for the presence of Werren, J. H. et al. (2008) ‘Wolbachia: master Wolbachia like DNA sequences in the S. manipulators of invertebrate biology’, 6, pp 741- digitata genome. 751. Therefore, based on the BALST results, Acknowledgment phylogenetic analysis and the sequence Professor R.S. Dassanayake and Ms. N similarity analysis, it can be concluded that as Rashanthy have been acknowledged for the 90% of the filarial nematodes studied up to guidance and support. The study was funded date contain the Wolbachia, ancestors of S. by the grant, AP/3/2017/CG/34, awarded to digitata may have colonized with Wolbachia in Professor R.S. Dassanayake by the University the distant past, and HGT may have brought of Colombo, Sri Lanka, to promote Wolbachia DNA into the nuclear genome prior collaborative research. to endosymbiont loss. Author References

Desjardins, C. A. et al. (2013) ‘Genomics of Loa loa, M.S.A. Kothalawala, BSc (Honours) in a Wolbachia -free filarial parasite of humans’, Biochemistry and Molecular Biology, Nature Publishing Group. Nature Publishing Group, University of Colombo, Sri Lanka 45(5), pp. 495–500. Fenn, K. et al. (2006) ‘Phylogenetic Relationships of the Wolbachia of Nematodes and Arthropods’, 2(10).

19