EDP Sciences logo
Submit your paper
Search
Menu
All issues Volume 29 (2022) Parasite, 29 (2022) 22 Full HTML
Open Access
Issue
Parasite
Volume 29, 2022
Article Number 22
Number of page(s) 17
DOI https://doi.org/10.1051/parasite/2022022
Published online 27 April 2022

© P. Nooroong et al., published by EDP Sciences, 2022

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Introduction

Leucocytozoon sabrazesi Mathis & Léger, 1911 is an important blood parasite belonging to the phylum Apicomplexa, which commonly infects a wide range of avian species. In addition, L. sabrazesi has frequently been reported in both fighting cocks (Gallus gallus) and domestic chickens (Gallus gallus domesticus) [2, 32, 41]. Both black fly (Simuliidae) and culicoides midges (Ceratopogonidae) act as potential vectors for Leucocytozoon transmission [20, 21, 36, 46–48]. Leucocytozoon sabrazesi infections of domestic stock cause symptoms including lethargy, green feces, loss of appetite, anemia, and death. Further, infections are known to cause economic losses through increased chicken mortality and reduced egg production [3, 28, 35]. Notably, leucocytozoonosis or Leucocytozoon infection are reported in many kinds of birds around the world, including in Asia (Thailand), Africa, Europe, and North America [6, 10, 39, 40, 45].

A conventional diagnosis of Leucocytozoon infection is based on microscopic examination of the gametocytes in Giemsa-stained blood smears of the infected chickens. Currently, polymerase chain reaction (PCR) may be more reliable and widely used to diagnose the infection and be supplemented by the standard parasitological method, especially in the laboratory for high sensitivity and specificity even when blood smears are negative with low parasitemia. Although the diversity of hemosporidian parasites has been demonstrated based on mitochondrial genes, such as cytochrome b (cytb), in ecological and evolutionary studies [4, 17], there is little information about the genetic diversity of L. sabrazesi isolates in Thailand with Plasmodium spp. co-infection when using mitochondrial genes (cytb, coxI and coxIII). Therefore, this study aimed to investigate the mitochondrial genetic diversity of L. sabrazesi and Plasmodium spp. coinfections in chickens in Thailand at these three loci, including phylogenetic and biogeographic relationships. In addition, the phylogenetic relationship, haplotype diversity, entropy, and geographic and evolutionary distribution among the isolates identified in this work and those from other countries are presented.

Materials and methods

Ethics statement

Experimentation on animals was carried out under the following approval and permit from the Animal Care and Use Committee (IMBMU-ACUC), Institute of Molecular Biosciences, Mahidol University, Thailand. All suitable international, national and/or institutional guidelines for animal care and use were followed. Also, we received consent to collect chicken blood samples at the animal farm.

Blood sample collection

Thirty chickens (Gallus gallus domesticus) from the Bongti (14°04′20.8″N 98°59′50.1″E) and Tha Sao (14°10′27.2″N 99°07′15.6″E) districts in Kanchanaburi province, Thailand were collected via the brachial wing vein. The blood samples were kept in sterile 1.5-mL tubes containing lithium heparin to prevent coagulation, and stored at −80 °C until use.

Ficoll density gradient centrifugation

For Giemsa-stained blood smears, elongated gametocytes of Leucocytozoon sabrazesi were detected in 30 chicken blood samples. The blood samples were diluted with 0.1 M phosphate-buffered saline (PBS), pH 7.4 and overlayed with Ficoll-Paque (Sigma-Aldrich, Burlington, MA, USA). They were centrifuged at 400 ×g for 30 min at 25 °C. The gametocytes were gently harvested by inserting the pipette directly through the upper layer and later washed twice in the PBS solution.

Leucocytozoon sabrazesi DNA extraction

The genomic DNA of L. sabrazesi in blood samples was extracted by using an E.Z.N.A.® Tissue DNA Kit (OMEGA Bio-Tek, Norcross, GA, USA) following the protocol of Watthanadirek et al. [42, 43] and Junsiri [22] with some modifications. Briefly, 250 μL of blood samples were mixed thoroughly with 25 μL of proteinase K solution, and incubated at 70 °C for 10 min. Then, 250 μL of absolute ethanol were added and all lysates transferred to the HiBind® DNA Mini Column. The lysates were centrifuged at maximum speed for 1 min before adding 500 μL of HBC buffer. After adding 700 μL of DNA washing buffer, the genomic material was eluted with 50 μL of elution buffer. Finally, the extracted DNA solutions were stored at −20 °C until further use.

Molecular amplification of L. sabrazesi DNA

The cytb, coxI, and coxIII genes of L. sabrazesi were amplified by nested PCR using the specific primers: Hemo_cytbF (5′–CATATATTAAGAGAATTATGGAG–3′) and Hemo_cytbR (5′–ATAAAATGYTAAGAAATACCATTC–3′) (GenBank accession number AB299369) for the first step of amplification. At the second step of amplification, Ls_cytbF (5′–CACC TAATCACATGGGTTTGTGGA–3′) and Ls_cytbR (5′–GCTTTGGGCTAAGAATAATACC–3′) for the cytb gene, PgCoxIF (5′–CACCGCGTACTTTGGACCGAAAAA–3′) and PgCoxIR (5′–CATCCAGTACCACCACCAAA–3′) for the coxI gene, as well as CoxIII F (5′–CACCTAA CAT TCT ACA TGA TGT AGT–3′) and CoxIII R (5′–GTAAAAGCACACTTATCTAG–3′) for the coxIII gene were used in this study. The 4 base sequences (CACC) were added to the 5′ end of forward primers with overhang sequence (GTGG) in a pET100/D-TOPO® vector (Invitrogen, Waltham, MA, USA) to certify a cloning direction. The PCR reaction mixture contained 10× Standard Taq Reaction Buffer, 10 mM of each deoxynucleotide triphosphate (dNTPs), 10 μM of forward and reverse primers, 0.625 U of Taq DNA polymerase (NEB, UK), RNase-free water, and 1 μg of DNA template. The thermal cycling was performed in Mastercycler® nexus Thermal Cycler (Germany) with 95 °C for 2 min, followed by 35 cycles of 95 °C for 30 s, 55 °C for 30 s, 68 °C for 1 min, and then 68 °C for 5 min. The RNase-free water and confirmed L. sabrazesi DNA samples were used as negative and positive controls, respectively. The PCR products were separated by 1% agarose gel electrophoresis and stained with SYBR green fluorescence dye, then visualized under ultraviolet light. The positive samples were purified using a PureDireX PCR Clean-Up & Gel Extraction Kit (Bio-Helix Co., Taiwan).

Cloning and sequencing of the L. sabrazesi cytb, coxI and coxIII genes

The PCR products were purified using a PureDireX PCR Clean-Up & Gel Extraction Kit. The 5′ blunt end of purified PCR products was ligated into a pET100/D-TOPO® vector (Invitrogen Life Technologies, Carlsbad, CA, USA). The ligation products were heat-shocked and transformed into chemically competent Escherichia coli host strain TOP10 cells (Invitrogen Life Technologies). The transformed E. coli cultures were spread on Luria-Bertani (LB) agar plates containing 100 μg of ampicillin and incubated at 37 °C for 16 h. The positive bacterial colonies were picked and cultured in LB media containing ampicillin with shaking at 37 °C for 16 h. The plasmids were extracted from bacterial cultures using an AxyPrep Plasmid Miniprep Kit (Axygen Bioscience, Union City, CA, USA) before sequencing.

Sequence and in silico analysis

The presence of cytb, coxI and coxIII inserts was confirmed by Sanger sequencing. All sequences were submitted and deposited in the National Center for Biotechnology Information (NCBI) GenBank database. The sequences were also analyzed by BLAST (https://blast.ncbi.nlm.nih.gov). All nucleotide and amino acid sequences were analyzed by the computer programs MEGA 7.0.26 [24] and Jalview [10]. For nucleic acid substitution analysis, nucleotide diversity was determined using DnaSP software, V.6.0 [26]. All base substitutions were determined as synonymous and nonsynonymous substitutions in nucleotides and amino acid sequences were assessed using PROVEAN analysis [9] as compensation of physicochemical properties of amino acid replacement. In addition, the haplotype analysis was determined through DnaSP software, V.6.0 [26] before visualization of the mutational occurrence of haplotypes from different geographic distribution, and the relationships among haplotypes were visualized with a TCS network in the popART program [25].

Multiple sequence alignment and phylogenetic analysis

The cytb, coxI and coxIII sequences were employed for sequence alignment and phylogenetic analysis. Multiple sequence alignments were conducted with the MUSCLE algorithm [12]. All aligned DNA sequences were used to construct the molecular phylogenetic trees using neighbour-joining (NJ), maximum likelihood (ML), maximum parsimony (MP) and Bayesian analysis (BA) [19]. The reliability of the internal branching pattern of the phylogenetic tree was determined in each clade by statistical calculation of 1000 replicates using the bootstrapping method [13] and MrBayes program for posterior probability. The evolutionary distances were evaluated by the Kimura 2-parameter method [23]. Similarity (as a percentage) was also analyzed by using a sequence identity matrix in BioEdit software V.7.0.5.3 [16].

Entropy analysis

The entropy values for nucleotide and amino acid variation were assessed with Shannon’s entropy (H(x)) plot method in BioEdit software, V.7.0.5.3 [1, 16].

Results

Determination of L. sabrazesi mitochondrial gene sequences

The DNA sequences of L. sabrazesi cytb, coxI and coxIII were partially amplified by nested PCR. The quality of PCR products was evaluated by the ratio of optical density (OD260/280) of 1.8–2.0, which showed no contamination of the products. The lengths of cytb, coxI and coxIII sequences Thailand strain were 248, 588 and 294 bp, respectively. All DNA sequences of L. sabrazesi investigated in this study were submitted and deposited in the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/genbank/) under accession numbers MZ634375 to MZ634390 for the cytb gene, MZ634391 to MZ634403 for the coxI gene, and MZ634404 to MZ634417 for the coxIII gene (Table 1).

Table 1

The L. sabrazesi and Plasmodium spp. mitochondrial nucleotide sequences amplified in Thailand deposited in GenBank.

Phylogenetic analysis

The L. sabrazesi cytb sequences obtained in this work were aligned with other sequences retrieved from GenBank including sequences from Thailand, Malaysia, Myanmar, China, USA, Uganda, Congo, Sri Lanka, Brazil, Philippines, UK and Japan. Our sequences detected in this study were positioned in the same clade as L. sabrazesi (Fig. 1). The Thailand coxI sequences were determined in the different clades in the phylogenetic tree together with other sequences of P. gallinaceum and P. juxtanucleare (Fig. 2), while the phylogenetic tree constructed from coxIII sequences was positioned in the same clade as L. sabrazesi (Fig. 3). Not only the phylogenetic tree constructed from each mitochondrial sequence, but also the concatenated genes from all mitochondrial sequences were used to construct the phylogenetic tree, which showed that our sequences were grouped and positioned in the same clade as L. sabrazesi (Fig. 4). Moreover, the reliability of bootstrap frequencies and Bayesian posterior probabilities of all phylogenies are displayed with the highest values on each branch.

thumbnail Fig. 1

Phylogenetic tree of the cytb gene sequences in this study (bold face) and those taken from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.

thumbnail Fig. 2

Phylogenetic tree of the coxI gene sequences in this study (bold face) and those obtained from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.

thumbnail Fig. 3

Phylogenetic tree of the coxIII gene sequences in this study (bold face) and those taken from GenBank.The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.

thumbnail Fig. 4

Phylogenetic tree of the concatenated gene sequences in this study (bold face) and those obtained from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.

Similarity analysis

All DNA samples from chickens were positive for all three mitochondrial genes. The similarity among the Thailand cytb, coxI and coxIII sequences taken in this study was 99.9–100%, 97–100% and 98–100%, respectively (Tables S1–S3), while the similarity of those compared only between Plasmodium sequences obtained from GenBank by BLAST was 85–100%, 86–100% and 84–100%, respectively (Tables S1–S3). Interestingly, one sequence in this study showed 100% similarity co-infection of L. sabrazesi and P. gallinaceum (Tables 1 and S2). As well, our two sequences exhibited 100% similarity of co-infection of L. sabrazesi and P. juxtanucleare (Tables 1 and S2). For amino acid sequencing of L. sabrazesi, the similarity among the Thailand cytb, coxI and coxIII sequences taken in this study was 98–100%, 97–100% and 96–100%, respectively, whereas the similarity of these compared with other sequences obtained from GenBank by BLAST was 74–100%, 60–99% and 60–100%, respectively (Tables S4–S6).

Entropy analysis

The similarity analysis from Simplot showed higher nucleotide variation in Plasmodium spp. than in L. sabrazesi. The entropy analysis of the cytb, coxI and coxIII genes showed more variation of nucleic acid sequences than amino acid sequences. To analyze nucleic acid entropy, cytb, coxI and coxIII sequences showed 81 peaks with entropy values ranging from 0.11691 to 0.93764, 174 peaks with entropy values ranging from 0.13579 to 1.06709, and 125 peaks with entropy values ranging from 0.14614 to 0.94469, respectively. Entropy analysis of amino acid sequences exhibited that the charts showed 24 peaks with entropy values ranging from 0.11691 to 1.05331 for cytb, 62 peaks with entropy values ranging from 0.13579 to 1.61397 for coxI, and 46 peaks with entropy values ranging from 0.14614 to 1.18722 for coxIII (Fig. 5). The coxI gene was found to be more diverse than cytb and coxIII and this is consistent with multiple sequence alignment which showed more similarity among amino acid sequences than nucleic acid sequences (Supplementary Figs. 1–3). The nucleic acid variation from multiple sequence alignment correlated to high nucleic acid diversity in the coxI gene caused by nucleic acid sequences of L. sabrazesi. Besides the coxI gene, both the cytb and coxIII genes exhibited higher genetic diversity in Plasmodium spp. than in L. sabrazesi (Tables 2 and 3).

thumbnail Fig. 5

Entropy analysis of L. sabrazesi cytb, coxI and coxIII gene sequences. Entropy plot of multiple nucleic acid sequence alignment of the cytb (A), coxI (B) and coxIII (C) genes. The red peaks indicate the high variation at each position of the nucleic acid sequences. Entropy plot of multiple amino acid sequence alignment of CYTb (D), COXI (E) and COXIII (F). The red peaks indicate the high variation at each position of amino acid sequences.

Table 2

Comparison of nucleotide sequence analyses of three mitochondrial and concatenated genes of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand and other countries.

Table 3

Polymorphism and genetic diversity of the three mitochondrial and concatenated genes of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand and other countries.

Nucleic acid substitution analysis

Each nucleic acid substitution of cytb, coxI and coxIII was validated as transition from purine to purine and from pyrimidine to pyrimidine. In addition, the percentage of base composition of these genes indicated the number of A and T bases greater than G and C contents. However, most base substitutions were indicated as the synonymous substitutions (Fig. 6). Moreover, the synonymous frequency (Ks) of these genes was higher than non-synonymous frequency values (Ka). The Ka/Ks ratios of cytb, coxI, coxIII and concatenated genes were 0.13, 0.168, 0.227 and 0.181, respectively (Table 2). While all results of the evolutionary estimation of Tajima’s D values exhibited minus values, only coxI showed statistical significance, which determined an excess of low frequency polymorphisms relative to expectations under the neutral model of evolution (p < 0.10) (Table 2). In addition to Tajima D values, the Fu’s Fs statistic based on the distribution of haplotypes displayed minus values, indicating an excess of rare haplotypes over what would be expected under neutrality; especially coxI exhibited significant negative values of both Tajima D and Fu’s Fs statistic (p < 0.10) (Table 2). Each base non-synonymous substitution was analyzed in regards to the compensation of physicochemical properties of amino acid replacement. The cytb gene was found to have two positions of hydrophobic amino acid replacement from I15V and L66V. In the case of coxI, there were five amino acid replacements in the L. sabrazesi population, including R6I, Y32C, K56N, S99T and A113V, while Plasmodium spp. were found to have 23 amino acid substitutions, including R6K, Y32K, R50 K, N58T, N58K, N61K, K66I, L67H, I71M, S73F, L74F, F81L, C93W, P95S, K97E, P99A, K102R, I103L, Q111H, G116E, L117F, F119I, P123A, S123A, F126C and F126Y. The coxIII gene was found to have six amino acid replacements, including T4H, T4P, L5I, L32I, S55F, I69T and I79S. However, all amino acid replacements exhibiting the most conservative replacements occurred by non-synonymous substitution.

thumbnail Fig. 6

Nucleic acid substitution rate and base composition of cytb, coxI, coxIII and concatenated sequences among Leucocytozoon spp. and Plasmodium spp. Tables showing the transition and transversion from nucleotide substitution in cytb (A), coxI (B.), coxIII (C) and concatenated (D) genes. Graph incidating the synonymous and non-synonymous substitutions of cytb (E), coxI (F), coxIII (G) and concatenated (H) genes of Leucocytozoon spp. and Plasmodium spp.

Haplotype diversity

The TCS Network tool was used to construct the haplotype network of the cytb, coxI and coxIII gene sequences of Leucocytozoon spp. and Plasmodium spp. The haplotype of each gene was estimated together with geographic distribution, consistently displaying high variation from multiple sequence alignment. The coxI gene showed a greater number of nucleotide variations and higher diversity than coxIII and cytb. However, L. sabrazesi harbored 4, 8 and 10 haplotypes of cytb, coxI and coxIII, respectively. For L. sabrazesi cytb gene Thailand strain, our findings showed that most sequences are found in haplotype #1 and some sequences are found in haplotypes #3 and #4 obtained from Myanmar and Malaysia (Fig. 7, Tables 2 and 3). In the case of coxI, L. sabrazesi Thailand strain contained seven haplotypes, including haplotypes #1 to #5 and #10 to #11 formed the nearest clade with haplotype #14 of L. sabrazesi Malaysia strain. In addition to the coxI gene of L. sabrazesi, haplotype #9 of P. gallinaceum from Thailand formed the nearest branch to haplotype #30 of P. gallinaceum from the Philippines. Five haplotypes of P. juxtanucleare from Thailand, including haplotypes #6 to #8, #12 and #13 also formed the nearest branch to haplotype #28 of P. juxtanucleare from Japan (Fig. 8, Tables 2 and 3). Additionally, nine haplotypes of L. sabrazesi coxIII gene Thailand strain exhibited the nearest branch to haplotype #10 of L. sabrazesi Malaysia strain (Fig. 9, Tables 2 and 3). The concatenated gene comprising eight haplotypes in L. sabrazesi Thailand strain also grouped together with haplotype #9 in L. sabrazesi Malaysia strain (Fig. 10, Tables 2 and 3).

thumbnail Fig. 7

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. cytb gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).

thumbnail Fig. 8

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. coxI gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).

thumbnail Fig. 9

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. coxIII gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).

thumbnail Fig. 10

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. concatenated gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).

Discussion

Leucocytozoonosis caused by the hemoprotozoan L. sabrazesi is an important insect-borne disease of chickens and causes high economic losses to chicken industries worldwide, including in Thailand. In general, genetic diversity is a survival strategy which is employed by parasites to evade the immune responses of avian hosts (chickens, ducks and birds) [5, 37]. There have been studies of genetic diversity of Leucocytozoon sp. based on the mitochondrial gene sequences in several countries, and almost all of these studies focused on the cytb gene [8, 18, 29, 44]. However, there has been no information available regarding the genetic diversity and phylogeny of L. sabrazesi mitochondrial genes in Thailand until now. In the present study, we used the cytb, coxI, coxIII and concatenated genes in the chicken population sampled in Thailand to ascertain the genetic diversity of L. sabrazesi and their co-infections in these regions.

The molecular detection and DNA sequencing displayed the highest similarity of both cytb and coxIII genes of L. sabrazesi. Interestingly, this is the first report of co-infection between L. sabrazesi and P. gallinaceum and that of L. sabrazesi and P. juxtanucleare in the leucocytes of chickens in Thailand. Notably, the coxI gene has the ability to cross-react and could be used to detect infection of L. sabrazesi and Plasmodium spp. Our findings are consistent with the report obtained by Pacheco et al. [32]. A phylogenetic analysis was carried out to display the relationship between individual and multi-locus genes of mitochondria determining the detection of L. sabrazesi. Moreover, the coxI gene has been employed to detect te infection of P. gallinaceum and P. juxtanucleare in chickens from Bongti and Tha Sao districts in Kanchanaburi province located near the Chacheongsao province of Thailand which are reported about P. gallinaceum [34] and near at the border of Myanmar which are reported regarding P. juxtanucleare in chickens [44]. Regarding three mitochondrial nucleotide sequences, our results indicated the highest sequence similarity to L. sabrazesi and some co-infected with P. gallinaceum and P. juxtanucleare.

Genetic variation of three mitochondrial genes commonly occurred in Plasmodium spp., while coxI showed high genetic variation in Leucocytozoon spp. However, these genes were found to have higher transition than transversion rates, and caused mutational bias to high A-T content and were proned to express the evolutionary saturation for divergence of parasites, which are consistent with the analysis of hemosporidian mitochondrial genomes [33]. Moreover, the lack of mitochondrial sequences from Leucocytozoon spp. and Plasmodium spp. directly affected the evolutionary analysis. These genes displayed Ka/Ks ratios less than one and minus values, indicating purifying selection [30]. Tajima’s D results indicated minus values, but only coxI indicated selective sweep, which was consistent with the negative value of Fu’s Fs statistic which determined the population expansion under statistical significance [15]. In addition, the cytb and coxIII genes indicated minus values of Ka/Ks ratios that determined purifying selection, but both Tajima’s D and Fu’s Fs were negative and not significant, indicating neutrality or perhaps these values can result in indirect selection from balancing selection on a nearby locus (linked genes) [38]. All evolutionary analyses reflected that hemosporidian organisms passed through the important obstacle of evolution like genetic drift before performing population expansion later [7, 14]. In addition, some variations affected haplotype distribution, which occurred from the polyphyletic relationship of genus Leucocytozoon spp. and likely displayed as an ancestor of avian parasites [27]. In addition, only partial nucleotide sequences exhibited the number of synonymous greater than non-synonymous substitution, and amino acid replacement caused by non-synonymous substitution did not show lethal effects to L. sabrazesi and mitochondrial genome variation caused by the host switching during their life cycle [33]. However, the number of non-synonymous substitutions affecting amino acid replacements exhibited a higher number of conservative than radical amino acid replacements, reflecting the purifying selection of mitochondrial genes [11]. In addition to nucleotide substitution, the non-synonymous substitutions which caused the amino acid substitutions were estimated concerning the compensation of amino acids by physicochemical properties through PROVEAN program. We found that all amino acid substitutions did not affect their function based on comparisons from the NCBI database. Moreover, the multiple amino acid sequence alignment also consistently displayed compensation of physicochemical properties of amino acid replacement through BLOSUM 62 score and exhibited a higher number of conservative than radical amino acid replacements with the dark blue color (Supplementary Figs. 1B, 2B, 3B).

In this study, the entropy and multiple sequence analysis showed the genetic variation in Plasmodium spp. to be greater than in Leucocytozoon spp. This indicated greater diversification of malaria parasites and the paraphyletic relationship among avian hemosporidians [31]. However, Leucocytozoon spp. displayed a higher number of haplotypes than Plasmodium spp., these values were affected in populations of Leucocytozoon sp. but not L. sabrazesi [37]. Similarly, haplotype diversity indicated the close genetic relationship among L. sabrazesi. detected in Thailand, Malaysia and Myanmar [44].

Conclusions

This study is the first report on the genetic diversity of L. sabrazesi based on the mitochondrial genes including cytb, coxI, coxIII and concatenated sequences in Thailand. The co-infection between L. sabrazesi either P. gallinaceum or P. juxtanucleare in chickens in Thailand was investigated. The advantage of cross-PCR amplification of the coxI gene is that it can discriminate co-infection, which is not verified by microscopic examination. Even though the phylogenetic relationship and evolutionary distribution showed high genetic variation and haplotype diversity in the coxI, coxIII and cytb genes, they still indicated purifying selection, which occurred together with population expansion after genetic drift events in switching-host hemosporidian populations. These findings could help to improve the understanding of molecular phylogenetics and diversity among these mitochondrial sequences of L. sabrazesi Thailand strain. Our findings could therefore be beneficial for the development of immunodiagnostic tools and vaccine strategies for chicken leucocytozoonosis.

Supplementary materials

thumbnail Supplementary Figure 1.

Alignment of nucleic acid sequences of the cytb gene among L. sabrazesi and Plasmodium spp. The highest similarity of nucleotide positions is represented with dark blue color, while white color represents the least similarity of each nucleic acid position (A). Multiple amino acid sequence alignment of CYTb protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).

thumbnail Supplementary Figure 2:

Alignment of nucleic acid sequences of coxI gene among L. sabrazesi and Plasmodium sp.p The highest similarity of nucleotide positions is represented with dark blue color, while white color represents the least similarity of each nucleic acid position (A). Multiple amino acid sequence alignment of COXI protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).

thumbnail Supplementary Figure 3:

Alignment of nucleic acid sequences of coxIII gene among L. sabrazesi and Plasmodium sp. The highest similarity of nucleotide positions is representd with dark blue color, while white color represents the least similarity of each nucleic acid position(A). Multiple amino acid sequence alignment of COXIII protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).

Table S1: Similarity of the cytb gene sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences taken from GenBank.

Table S2: Similarity of the cox I gene sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences obtained from GenBank.

Table S3: Similarity of the cox III gene sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences taken from GenBank.

Table S4: Similarity of the cytb amino acid sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences taken from GenBank.

Table S5: Similarity of the cox I amino acid sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences taken from GenBank.

Table S6: Similarity of the cox III amino acid sequences of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand compared with other sequences taken from GenBank.

Access here

Acknowledgments

This research project was supported financially by Mahidol University (Basic Research Fund: fiscal year 2021) and a Research Grant from the National Research Council of Thailand (NRCT) and the Mid-Career Research Grant co-founded by the National Research Council of Thailand (NRCT) and Mahidol University [grant number NRCT5-RSA63015-23] to Panat Anuracpreeda.

References

  1. Alzohairy A. 2011. BioEdit: An important software for molecular biology. GERF Bulletin of Biosciences, 2, 60–61. [Google Scholar]
  2. Anderson JR, Trainer DO, Defoliart GR. 1962. Natural and experimental transmission of the waterfowl parasite, Leucocytozoon simondi M. & L, in Wisconsin. Zoonoses Research, 1, 155–164. [PubMed] [Google Scholar]
  3. Anisuzzaman A. 2018. Prevalence and pathology of haemoprotozoan infection in chicken. Bangladesh Journal of Agricultural Research, 6, 79–85. [Google Scholar]
  4. Bensch SJ, Pérez-Tris J, Waldenstro M, Hellgren O. 2004. Linkage between nuclear and mitochondrial DNA sequences in avian malaria parasites: Multiple cases of cryptic speciation. Evolution, 58, 1617–1621. [CrossRef] [PubMed] [Google Scholar]
  5. Bensch S, Hellgren O, Pérez-Tris J. 2009. MalAvi: a public database of malaria parasites and related haemosporidians in avian hosts based on mitochondrial cytochrome b lineages. Molecular Ecology Resources, 9, 1353–1358. [CrossRef] [PubMed] [Google Scholar]
  6. Buranapim N, Chaiwisit P, Wangkawan A, Tiwananthagorn S. 2019. A survey on blood parasites of birds in Chiang Mai province. Veterinary Integrative Sciences, 17, 65–73. [Google Scholar]
  7. Chang HH, Moss EL, Park DJ, Ndiaye D, Mboup S, Volkman SK, Sabeti PC, Wirth DF, Neafsey DE, Hartl DL. 2013. Malaria life cycle intensifies both natural selection and random genetic drift. Proceedings of the National Academy of Sciences of the United States of America, 110, 20129–20134. [CrossRef] [PubMed] [Google Scholar]
  8. Chawengkirttikul R, Junsiri W, Watthanadirek A, Poolsawat N, Minsakorn S, Srionrod N, Anuracpreeda P. 2021. Molecular detection and genetic diversity of Leucocytozoon sabrazesi in chickens in Thailand. Scientific Reports, 11, 16686. [CrossRef] [PubMed] [Google Scholar]
  9. Choi Y, Chan AP. 2015. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics (Oxford, England), 31, 2745–2747. [CrossRef] [PubMed] [Google Scholar]
  10. Clamp M, Cuff J, Searle SM, Barton GJ. 2004. The Jalview Java alignment editor. Bioinformatics, 20, 426–427. [CrossRef] [PubMed] [Google Scholar]
  11. Dagan T, Talmor Y, Graur D. 2002. Ratios of radical to conservative amino acid replacement are affected by mutational and compositional factors and may not be indicative of positive Darwinian selection. Molecular Biology and Evolution, 19, 1022–1025. [CrossRef] [PubMed] [Google Scholar]
  12. Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research, 32, 1792–1797. [CrossRef] [PubMed] [Google Scholar]
  13. Felsenstein J. 1985. Confidence limits on phylogenies: An approach using the bootstrap. Evolution, 39, 783–791. [CrossRef] [PubMed] [Google Scholar]
  14. Ferretti L, Raineri E, Ramos-Onsins S. 2012. Neutrality tests for sequences with missing data. Genetics, 191, 1397–1401. [CrossRef] [PubMed] [Google Scholar]
  15. Fu FX. 1997. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics, 147, 915–925. [CrossRef] [PubMed] [Google Scholar]
  16. Hall TA. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series, 41, 95–98. [Google Scholar]
  17. Hellgren O, Waldenström J, Bensch S. 2004. A new PCR assay for simultaneous studies of Leucocytozoon, Plasmodium, and Haemoproteus from avian blood. Journal of Parasitology, 90, 797–802. [CrossRef] [PubMed] [Google Scholar]
  18. Hellgren O, Križanauskienė A, Valkiūnas G, Bensch S. 2007. Diversity and phylogeny of mitochondrial cytochrome b lineages from six morphospecies of avian Haemoproteus (Haemosporida: Haemoproteidae). Journal of Parasitology, 93, 889–896. [CrossRef] [PubMed] [Google Scholar]
  19. Huelsenbeck JP, Ronquist F. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics, 17, 754–755. [CrossRef] [PubMed] [Google Scholar]
  20. Ishtiaq F, Gering E, Rappole JH, Rahmani AR, Jhala YV, Dove CJ, Milensky C, Olson SL, Peirce MA, Fleischer RC. 2007. Prevalence and diversity of avian hematozoan parasites in Asia: a regional survey. Journal of Wildlife Diseases, 43, 382–398. [CrossRef] [PubMed] [Google Scholar]
  21. Jumpato W, Tangkawanit U, Wongpakam K, Pramual P. 2019. Molecular detection of Leucocytozoon (Apicomplexa: Haemosporida) in black flies (Diptera: Simuliidae) from Thailand. Acta Tropica, 190, 228–234. [CrossRef] [PubMed] [Google Scholar]
  22. Junsiri W, Watthanadirek A, Poolsawat N, Kaewmongkol S, Jittapalapong S, Chawengkirttikul R, Anuracpreeda P. 2020. Molecular detection and genetic diversity of Anaplasma marginale based on the major surface protein genes in Thailand. Acta Tropica, 105338. [CrossRef] [PubMed] [Google Scholar]
  23. Kimura M. 1980. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution, 16, 111–120. [CrossRef] [PubMed] [Google Scholar]
  24. Kumar S, Stecher G, Tamura K. 2016. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Molecular Biology and Evolution, 33, 1870–1874. [CrossRef] [PubMed] [Google Scholar]
  25. Leigh J, Bryant D. 2015. PopART: Full-feature software for haplotype network construction. Methods in Ecology and Evolution, 6, 1110–1116. [CrossRef] [Google Scholar]
  26. Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25, 1451–1452. [CrossRef] [PubMed] [Google Scholar]
  27. Martinsen ES, Perkins SL, Schall JJ. 2008. A three-genome phylogeny of malaria parasites (Plasmodium and closely related genera): Evolution of life-history traits and host switches. Molecular Phylogenetics and Evolution, 47, 261–273. [CrossRef] [PubMed] [Google Scholar]
  28. Mirzaei F, Siyadatpanah A, Norouzi R, Pournasir S, Nissapatorn V, Pereira MD. 2020. Blood parasites in domestic birds in Central Iran. Veterinary Sciences, 7, 126. [CrossRef] [Google Scholar]
  29. Morii T, Shiihara T, Lee YC, Manuel MF, Nakamura K, Iijima T, Horii K. 1981. Seroimmunological and parasitological surveys of Leucocytozoon caulleryi infection in chickens in several Asian countries. International Journal for Parasitology, 11, 187–190. [CrossRef] [PubMed] [Google Scholar]
  30. Nekrutenko A, Makova KD, Li WH. 2002. The K(A)/K(S) ratio test for assessing the protein-coding potential of genomic regions: an empirical and simulation study. Genome Research, 12, 198–202. [CrossRef] [PubMed] [Google Scholar]
  31. Outlaw DC, Ricklefs RE. 2011. Rerooting the evolutionary tree of malaria parasites. Proceedings of the National Academy of Sciences of the United States of America, 108, 13183–13187. [CrossRef] [PubMed] [Google Scholar]
  32. Pacheco MA, Cepeda AS, Bernotienė R, Lotta IA, Matta NE, Valkiūnas G, Escalante AA. 2018. Primers targeting mitochondrial genes of avian haemosporidians: PCR detection and differential DNA amplification of parasites belonging to different genera. International Journal for Parasitology, 48, 657–670. [CrossRef] [PubMed] [Google Scholar]
  33. Pacheco MA, Matta NE, Valkiūnas G, Parker PG, Mello B, Stanley CE Jr, Lentino M, Garcia A, Maria A, Cranfield M, Kosakovsky P, Sergei L, Escalante AA. 2018. Mode and rate of evolution of haemosporidian mitochondrial genomes: Timing the radiation of avian parasites. Molecular Biology and Evolution, 35, 383–403. [CrossRef] [PubMed] [Google Scholar]
  34. Pattaradilokrat S, Tiyamanee W, Simpalipan P, Kaewthamasorn M, Saiwichai T, Li J, Harnyuttanakorn P. 2015. Molecular detection of the avian malaria parasite Plasmodium gallinaceum in Thailand. Veterinary Parasitology, 210, 1–9. [CrossRef] [PubMed] [Google Scholar]
  35. Piratae S, Vaisusuk K, Chatan W. 2021. Prevalence and molecular identification of Leucocytozoon spp. in fighting cocks (Gallus gallus) in Thailand. Parasitology Research, 120, 2149–2155. [CrossRef] [PubMed] [Google Scholar]
  36. Pramual P, Tangkawanit U, Kunprom C, Vaisusuk K, Chatan W, Wongpakam K, Thongboonma S. 2020. Seasonal population dynamics and a role as natural vector of Leucocytozoon of black fly, Simulium chumpornense Takaoka & Kuvangkadilok. Acta Tropica, 211, 105617. [CrossRef] [PubMed] [Google Scholar]
  37. Reeves AB, Smith MM, Meixell BW, Fleskes JP, Ramey AM. 2015. Genetic diversity and host specificity varies across three genera of blood parasites in ducks of the Pacific Americas Flyway. PLoS One, 10, e0116661. [CrossRef] [PubMed] [Google Scholar]
  38. Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio V, Wallace DC. 2004. Effects of purifying and adaptive selection on regional variation in human mtDNA. Science, 303, 223–226. [CrossRef] [PubMed] [Google Scholar]
  39. Sato Y, Hagihara M, Yamaguchi T, Yukawa M, Murata K. 2007. Phylogenetic comparison of Leucocytozoon spp. from wild birds of Japan. Journal of Veterinary Medical Science, 69, 55–59. [CrossRef] [PubMed] [Google Scholar]
  40. Singjam S, Ruksachat N. 2011. Case Report: Outbreak of leucocytozoonosis in captive wild birds. Khao Kor Wildlife Captive Breeding Center, Veterinary Research and Development Center. DLD Thailand, Report number 8, 1–7. [Google Scholar]
  41. Suprihati E, Yuniarti W. 2017. The phylogenetics of Leucocytozoon caulleryi infecting broiler chickens in endemic areas in Indonesia. Veterinary World, 10, 1324–1328. [CrossRef] [PubMed] [Google Scholar]
  42. Watthanadirek A, Chawengkirttikul R, Poolsawat N, Junsiri W, Boonmekam D, Reamtong O, Anuracpreeda P. 2019. Recombinant expression and characterization of major surface protein 4 from Anaplasma marginale. Acta Tropica, 197, 105047. [CrossRef] [PubMed] [Google Scholar]
  43. Watthanadirek A, Junsiri W, Minsakorn S, Poolsawat N, Srionrod N, Khumpim P, Chawengkirttikul R, Anuracpreeda P. 2021. Molecular and recombinant characterization of major surface protein 5 from Anaplasma marginale. Acta Tropica, 220, 105933. [CrossRef] [PubMed] [Google Scholar]
  44. Win SY, Chel HM, Hmoon MM, Htun LL, Bawm S, Win MM, Murata S, Nonaka N, Nakao R, Katakura K. 2020. Detection and molecular identification of Leucocytozoon and Plasmodium species from village chickens in different areas of Myanmar. Acta Tropica, 212, 105719. [CrossRef] [PubMed] [Google Scholar]
  45. Xuan MNT, Kaewlamun W, Saiwichai T, Thanee S, Poofery J, Tiawsirisup S, Poofery J, Tiawsirisup S, Channumsin M, Kaewthamasorn M. 2021. Development and application of a novel multiplex PCR assay for the differentiation of four haemosporidian parasites in the chicken Gallus gallus domesticus. Veterinary Parasitology, 293, 109431. [CrossRef] [PubMed] [Google Scholar]
  46. Zhao W, Cai B, Qi Y, Liu S, Hong L, Lu M, Chen X, Qiu C, Peng W, Li J, Su XZ. 2014. Multi-strain infections and “relapse” of Leucocytozoon sabrazesi gametocytes in domestic chickens in southern China. PLoS One, 9, e94877. [CrossRef] [PubMed] [Google Scholar]
  47. Zhao W, Liu J, Xu R, Zhang C, Pang Q, Chen X, Liu S, Hong L, Yuan J, Li X, Chen Y, Li J, Su XZ. 2015. The gametocytes of Leucocytozoon sabrazesi Infect chicken thrombocytes, not other blood cells. PLoS One, 10, e0133478. [CrossRef] [PubMed] [Google Scholar]
  48. Zhao W, Pang Q, Xu R, Liu J, Liu S, Li J, Su XZ. 2016. Monitoring the prevalence of Leucocytozoon sabrazesi in southern China and testing tricyclic compounds against gametocytes. PLoS One, 11, e0161869. [CrossRef] [PubMed] [Google Scholar]

Cite this article as: Nooroong P, Watthanadirek A, Minsakorn S, Poolsawat N, Junsiri W, Srionrod N, Sangchuai S, Chawengkirttikul R & Anuracpreeda P. 2022. Molecular genetic diversity and bioinformatic analysis of Leucocytozoon sabrazesi based on the mitochondrial genes cytb, coxI and coxIII and co-infection of Plasmodium spp. Parasite 29, 22.

All Tables

Table 1

The L. sabrazesi and Plasmodium spp. mitochondrial nucleotide sequences amplified in Thailand deposited in GenBank.

In the text
Table 2

Comparison of nucleotide sequence analyses of three mitochondrial and concatenated genes of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand and other countries.

In the text
Table 3

Polymorphism and genetic diversity of the three mitochondrial and concatenated genes of Leucocytozoon spp. and Plasmodium spp. as detected in chicken samples in Thailand and other countries.

In the text

All Figures

thumbnail Fig. 1

Phylogenetic tree of the cytb gene sequences in this study (bold face) and those taken from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.
In the text
thumbnail Fig. 2

Phylogenetic tree of the coxI gene sequences in this study (bold face) and those obtained from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.
In the text
thumbnail Fig. 3

Phylogenetic tree of the coxIII gene sequences in this study (bold face) and those taken from GenBank.The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.
In the text
thumbnail Fig. 4

Phylogenetic tree of the concatenated gene sequences in this study (bold face) and those obtained from GenBank. The boostrap values calculated from NJ, ML, MP and BA are labeled on each branch.
In the text
thumbnail Fig. 5

Entropy analysis of L. sabrazesi cytb, coxI and coxIII gene sequences. Entropy plot of multiple nucleic acid sequence alignment of the cytb (A), coxI (B) and coxIII (C) genes. The red peaks indicate the high variation at each position of the nucleic acid sequences. Entropy plot of multiple amino acid sequence alignment of CYTb (D), COXI (E) and COXIII (F). The red peaks indicate the high variation at each position of amino acid sequences.
In the text
thumbnail Fig. 6

Nucleic acid substitution rate and base composition of cytb, coxI, coxIII and concatenated sequences among Leucocytozoon spp. and Plasmodium spp. Tables showing the transition and transversion from nucleotide substitution in cytb (A), coxI (B.), coxIII (C) and concatenated (D) genes. Graph incidating the synonymous and non-synonymous substitutions of cytb (E), coxI (F), coxIII (G) and concatenated (H) genes of Leucocytozoon spp. and Plasmodium spp.
In the text
thumbnail Fig. 7

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. cytb gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).
In the text
thumbnail Fig. 8

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. coxI gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).
In the text
thumbnail Fig. 9

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. coxIII gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).
In the text
thumbnail Fig. 10

TCS network of haplotypes based on Leucocytozoon spp. and Plasmodium spp. concatenated gene sequences (A) detected in Thailand and other countries. The number of bars on lines between a haplotype and another represent the number of nucleotide mutation (B).
In the text
thumbnail Supplementary Figure 1.

Alignment of nucleic acid sequences of the cytb gene among L. sabrazesi and Plasmodium spp. The highest similarity of nucleotide positions is represented with dark blue color, while white color represents the least similarity of each nucleic acid position (A). Multiple amino acid sequence alignment of CYTb protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).
In the text
thumbnail Supplementary Figure 2:

Alignment of nucleic acid sequences of coxI gene among L. sabrazesi and Plasmodium sp.p The highest similarity of nucleotide positions is represented with dark blue color, while white color represents the least similarity of each nucleic acid position (A). Multiple amino acid sequence alignment of COXI protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).
In the text
thumbnail Supplementary Figure 3:

Alignment of nucleic acid sequences of coxIII gene among L. sabrazesi and Plasmodium sp. The highest similarity of nucleotide positions is representd with dark blue color, while white color represents the least similarity of each nucleic acid position(A). Multiple amino acid sequence alignment of COXIII protein among L. sabrazesi and Plasmodium spp. The highest similarity of physicochemical properties (BLOSUM score 62) of each amino acid position is represented with blue color, while white color represents the least similarity of each amino acid position (B).
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.