Skip to main content

Advertisement

Expanded complement of Niemann-Pick type C2-like protein genes in Clonorchis sinensis suggests functions beyond sterol binding and transport

Abstract

Background

The parasitic flatworm Clonorchis sinensis inhabits the biliary tree of humans and other piscivorous mammals. This parasite can survive and thrive in the bile duct, despite exposure to bile constituents and host immune attack. Although the precise biological mechanisms underlying this adaptation are unknown, previous work indicated that Niemann-pick type C2 (NPC2)-like sterol-binding proteins might be integral in the host-parasite interplay. Expansions of this family in some invertebrates, such as arthropods, have shown functional diversification, including novel forms of chemoreception. Thus, here we curated the NPC2-like protein gene complement in C. sinensis, and predicted their conserved and/or divergent functional roles.

Methods

We used an established comparative genomic-bioinformatic approach to curate NPC2-like proteins encoded in published genomes of Korean and Chinese isolates of C. sinensis. Protein sequence and structural homology, presence of conserved domains and phylogeny were used to group and functionally classify NPC2-like proteins. Furthermore, transcription levels of NPC2-like protein-encoding genes were explored in different developmental stages and tissues.

Results

Totals of 35 and 32 C. sinensis NPC2-like proteins were predicted to be encoded in the genomes of the Korean and Chinese isolates, respectively. Overall, these proteins had low sequence homology and high variability of sequence alignment coverage when compared with curated NPC2s. Most C. sinensis proteins were predicted to retain a conserved ML domain and a conserved fold conformation, with a large cavity within the protein. Only one protein sequence retained the conserved amino acid residues required in bovine NPC2 to bind cholesterol. Non-canonical C. sinensis NPC2-like protein-coding domains clustered into four distinct phylogenetic groups with members of a group frequently encoded on the same genome scaffolds. Interestingly, NPC2-like protein-encoding genes were predicted to be variably transcribed in different developmental stages and adult tissues, with most being transcribed in the metacercarial stage.

Conclusions

The results of the present investigation confirms an expansion of NPC2-like proteins in C. sinensis, suggesting a diverse array of functions beyond sterol binding and transport. Functional explorations of this protein family should elucidate the mechanisms enabling the establishment and survival of C. sinensis and related flukes in the biliary systems of mammalian hosts.

Background

Parasitic flatworms (phylum Platyhelminthes; class Trematoda) are responsible for neglected tropical diseases (NTDs) that affect more than 750 million humans and other mammals throughout the world, particularly in Southeast Asia and the Western Pacific regions [1,2,3]. Important trematode species include Clonorchis sinensis and Opisthorchis spp. [3, 4]. The consumption of fish containing the infective developmental stage (metacercaria) of these liver flukes leads to an infection that, if chronic, can cause serious hepatobiliary diseases in humans, including clonorchiasis (C. sinensis) or opisthorchiasis (Opisthorchis spp.), particularly in Asia [3, 5,6,7]. For example, C. sinensis impacts more than 35 million people across China, Japan, Korea and Vietnam [8, 9], and has been classified as a Class I carcinogen by the International Agency for Research on Cancer (IARC) [10]. Despite the importance of this parasite, the molecular mechanisms that govern or modulate the interactions between C. sinensis and its host animals remain largely unknown.

When piscivorous mammals (e.g. humans, dogs and cats) consume fish infected by C. sinensis, metacercariae excyst in the duodenum, and juveniles migrate and develop to adult flukes in the biliary system [11, 12]. The adult flukes can alter biliary duct pathology via mechanical irritation and by releasing molecules to facilitate feeding [13]. Furthermore, as the adult flukes grow and migrate, they can obstruct the bile duct and elevate bile duct pressure [6]. Clonorchiasis often leads to chronic hepatobiliary illness and can induce cholangiocarcinoma (CCA), a malignant cancer of the biliary system [1, 5, 14]. To better understand the pathogenesis of clonorchiasis and CCA, and to assist in efforts to control the parasite causing these diseases, researchers have explored the function of C. sinensis proteins, including proteins likely to be important for establishment and survival in the human biliary system.

In bile, cholesterol and phospholipids are abundant, and lipid-binding proteins (LBPs) have been proposed to play an important role in maintaining the chemical homeostasis of liver flukes in the bile duct [15]. Interestingly, gene duplication events have led to more than 20 copies of genes encoding homologues of a lipid-binding protein, Niemann-Pick type C2 (NPC2), in the genomes of C. sinensis and O. viverrini [16,17,18]. In most eukaryotes, NPC2 is encoded by a single copy gene. In some mammals, this protein binds cholesterol and other lipids and transports them out of the lysosome to other parts of the cell [19]. In arthropods, duplication and genetic differentiation of NPC2 genes [20,21,22,23] have led to new molecular functions, including chemoreception via binding to semiochemical and other volatile compounds [23,24,25]. Related proteins that contain a conserved myeloid differentiation factor-2 (MD-2)-related lipid-binding (ML) domain are also reported to play diverse roles in lipid metabolism, innate immunity and/or chemoreception in arthropods [20, 26]. Despite the substantial expansion of the NPC2-like protein family in liver flukes, almost nothing is known, at the molecular level, about their biological functions in such parasites and/or the regulation of host-pathogen interactions. Here, we employed a bioinformatic workflow to predict, curate and annotate NPC2-like proteins encoded in the draft genomes of a Korean and a Chinese isolate of C. sinensis.

Methods

Inference of NPC2-like protein sequences

Two C. sinensis genomes, one assembled from a Korean isolate (gene accession numbers beginning with “Cs-k2”; BioProject ID: PRJNA386618) [17] and one from a Chinese isolate (gene accession numbers beginning with “csin”; BioProject ID: PRJNA72781) [18], and their gene annotations were downloaded from the WormBase ParaSite database (v.13; accessed 20 May 2019) [27]. NPC2-like homologues in C. sinensis were initially identified (Fig. 1, steps 1-3). First, homologues of nine curated NPC2 proteins available in the SWISS-PROT database (Table 1; accessed 20 May 2019) [28] were identified in the C. sinensis proteomes using BLASTp v.2.2.29 (E-value cut-off: 10) [29]. Second, a reciprocal BLASTp (E-value cut-off: 10) search of homologues of NPC2 proteins from SWISS-PROT against the NCBI non-redundant protein database (NCBI-nr; accessed 20 May 2019) [30] was performed. Clonorchis sinensis NPC2-like protein homologues matching proteins submitted to NCBI-nr and annotated as “Niemann-pick C2 protein”, “NPC intracellular cholesterol transport 2”, “Epididymal secretory protein E1” or “Epididymal secretory protein E1-like”, or that lacked a protein description (unnamed or hypothetical proteins), were retained. Third, we searched for conserved domains in the predicted C. sinensis proteome using the program InterProScan v.5.15.54 [31], utilising the Pfam database v.27.0 [32]. Proteins with a conserved ML domain (Pfam identifier PF02221) were retained. Numbers of C. sinensis proteins with NPC2 homologues and/or a conserved ML domain were displayed in a Venn diagram employing the Intervene tool [33]. Heatmaps of BLASTp and InterProScan bit scores were created using the R package ggtree v.1.16.6 [34].

Fig. 1
figure1

Summary of workflow and results for the prediction, curation, and annotation of NPC2-like proteins in the genomes of a Korean and Chinese isolate of Clonorchis sinensis. Initially, putative C. sinensis NPC2-like proteins were identified using BLASTp sequence homology searches against curated NPC2 proteins submitted to the SWISS-PROT database (step 1). Selected proteins with a reciprocal best match to NPC2-like proteins submitted to NCBI-nr database were then identified (step 2). Next, C. sinensis predicted proteins with a conserved myeloid differentiation factor-2 (MD-2)-related lipid-binding (ML) domain (PF02221) were identified using InterProScan (step 3). The intersect between NPC2 protein sequence homology and retention of a conserved ML domain was then used to select putative NPC2-like proteins. Reciprocal BLASTn sequence homology across gene regions was used to infer orthology and incomplete homologous sequences were complemented in both isolates of C. sinensis using BLAT and Exonerate (step 5). Phylogenetic relationships among identified C. sinensis NPC2-like protein-encoding genes was then determined by employing MACSE2, TrimAl and MrBayes (step 6). Last, all NPC2-like proteins identified were functionally annotated using SignalP-5.0, SWISS-PROT and InterProScan databases, as well as I-TASSER, to predict tertiary structures. Transcription in different development stages and adult tissues was inferred using HISAT2 and FeatureCounts (step 7). The number of Korean (blue) and Chinese (red) C. sinensis sequences retained in each step are indicated

Table 1 Nine curated Niemann-Pick C2 (NPC2) proteins in the SWISS-PROT database used for identification of C. sinensis NPC2-like protein homologues. Protein length, signal peptide position and the most similar solved crystal structure are given for each sequence

Curation of sequences

All inferred C. sinensis NPC2-like proteins were individually curated in three steps (Fig. 1). First, gene regions (including exonic and intronic nucleotide sequences) encoding C. sinensis NPC2-like proteins were extracted from each of the two draft genomes (representing Chinese and Korean isolates) via the corresponding “general feature format” (GFF) file employing the gffread tool v.0.11.4 [35]. A reciprocal BLASTn v.2.2.29 search against the C. sinensis NPC2-like gene sets (for both isolates) was then conducted to infer gene orthology. Second, C. sinensis NPC2-like proteins from each isolate (Chinese or Korean) were reciprocally mapped to the genome of the alternative isolate using BLAT v.34x12 [36], and a new gene model was inferred based on this mapping employing the program Exonerate v.2.2.0 [37].

Analysis of transcription

Available RNA-Seq data (testis, sucker, muscle, ovary, adult, 8-week adult, metacercaria, 2-week juvenile; ENA/SRA accession numbers: ERR604978–ERR604981, SRR189060, SRR6188894–SRR6188896) for C. sinensis were mapped to each reference genome using HISAT2 [38]. From these mapped data (stored in the BAM format), read counts were inferred using the corresponding GFF files and employing FeatureCounts v.1.6.4 [39]. Read counts from each library were then normalised to counts per million (CPM) using the edgeR package v.3.26.8 [40]. A heatmap matrix of CPM per gene per library was created using the ggtree R package v.1.16.6 [34].

Protein annotation

Signal peptides and their cleavage sites were inferred for C. sinensis NPC2-like proteins using the SignalP-5.0 [41]. The structures of the mature C. sinensis NPC2-like proteins (i.e. without their predicted signal peptide domain) were modelled and annotated using the program I-TASSER v.4.4 [42] and compared with the crystal structures of NPC2s from Bos taurus (cow; Protein Data Bank (PDB) accession 2HKA chains A and C) [19] and Camponotus japonicus (Japanese carpenter ant; PDB accession 3WEA chain A) [43] employing the program UCSF Chimera v.1.9 [44].

Phylogenetic analysis

The mRNA sequences encoding C. sinensis NPC2-like proteins were extracted from each assembled genome sequence via their corresponding GFF files employing gffread. Codons of C. sinensis NPC2-like sequences were aligned using a translated protein sequence alignment employing MACSE v.2.03 [45]. Gaps were removed from the nucleotide alignment using trimAl v.1.4.rev15 [46] using the -gappyout option. The Akaike Information Criteria (AIC) test in ModelFinder [47] selected the general time reversible model of evolution for subsequent phylogenetic analyses. Bayesian phylogenetic inference (BI) was determined using Markov chain Monte Carlo (MCMC) analysis in MrBayes [48]. Two million generations of MCMC analysis were performed, and trees were recorded every 200th generation. At this point, the standard deviation of split frequencies was < 0.01, and the potential scale reduction factor (PSRF) approached 1. Consensus trees (50% majority rule) were generated using the final 75% of trees. Trees were annotated and enhanced using the ggtree R package v.1.16.6 [34], and nodal support values on trees were indicated as posterior probabilities (pp).

Results

NPC2-like proteins of Clonorchis sinensis

Based on predicted protein sequence homology, 72 of 14,538 Korean C. sinensis proteins and 71 of 13,634 Chinese C. sinensis proteins were homologous to one or more NPC2 proteins in SWISS-PROT (Fig. 1, Additional file 1: Table S1). Sequence homology between NPC2 proteins and their best matched C. sinensis proteins was mostly low (20.0–47.6 % amino acid identity), and sequence alignment coverage was highly variable (14–100%). In total, 55 Korean and 56 Chinese NPC2 homologues were reciprocal BLASTp matches (E-value cut-off: 10) to NCBI-nr proteins annotated as NPC2-like or those without a protein description and were thus retained (Fig. 1, Additional file 1: Table S1). Of those, 35 Korean and 27 Chinese C. sinensis proteins contained at least one ML domain, with an average conserved domain length of 107 amino acids (Fig. 1, Additional file 1: Table S1). Two conserved ML domains were predicted in three Korean C. sinensis protein sequences (Cs-k2.gene14549, Cs-k2.gene14290 and Cs-k2.gene14112). Based on protein sequence homology to NPC2 and the presence of a conserved ML domain (Fig. 1), 35 and 27 NPC2-like proteins were predicted to be encoded in the Korean and Chinese C. sinensis gene sets, respectively. One Chinese protein sequence (csin112467) that shared no significant sequence homology to NPC2 proteins, but contained a conserved ML domain, was also retained for subsequent curation (Fig. 1).

Gene models

Reciprocal nucleotide alignments of Korean and Chinese C. sinensis NPC2-like gene regions (including introns) and mapping to the alternative genome assemblies identified paired orthologues between the two isolates. Nine pairs of NPC2-like gene models were consistent in gene model structure and overall sequence length. Based on reciprocal nucleotide matches across exonic and intronic regions, we identified 11 additional orthologous pairs with unresolved variation in the first exon positions and lengths. Three of the Korean C. sinensis genes identified by reciprocal nucleotide alignments (i.e. Cs-k2.gene992, Cs-k2.gene8673 and Cs-k2.gene14547) were removed as they did not encode a conserved ML domain or share amino acid sequence homology with NPC2 proteins from SWISS-PROT. In contrast, four Chinese C. sinensis genes (csin101111, csin103126, csin111538 and csin111895) were added based on mapping of Korean NPC2-like gene sequences to the genomic region encoding these genes. Six Korean gene models, although aligning to the Chinese genome, only aligned partially to the corresponding Chinese gene model in this genomic region, and available sequence data were insufficient to resolve a single gene model for both isolates. Eight Korean C. sinensis genes were not identified in the Chinese genome, and four Chinese C. sinensis genes were not identified in the Korean genome. Taken together, 35 Korean and 32 Chinese NPC2-like protein genes were retained.

Tertiary structures and functional annotation

The tertiary structures of 35 Korean and 32 Chinese C. sinensis NPC2-like proteins were modelled using I-TASSER (Additional file 1: Table S2). The I-TASSER model confidence (C-) scores ranged from -5 (lowest confidence) to 1.29 (highest confidence). For comparison, I-TASSER models were also inferred for nine curated NPC2 proteins from SWISS-PROT (Table 1). For these proteins, the model C-scores ranged from 1.30 to 1.46. Eleven Korean and nine Chinese C. sinensis NPC2-like proteins had predicted structures with C-scores of ≥ 1. These high-confidence models were retained for further analysis. Two Korean and three Chinese NPC2-like proteins had C-scores of < -4 and were thus low-confidence predictions. Proteins Cs-k2.gene6404 (Korean isolate) and csin102672 (Chinese isolate) had the highest paired ortholog C-scores (1.26). Chinese NPC2-like proteins csin107773 and csin111438 had the highest (1.29) and lowest C-score (-5), respectively. Of note, six Korean and Chinese NCP2-like proteins with a C-score of < -4 were not predicted to retain a signal peptide region (Additional file 1: Table S2). Based on protein structurally close to the target in PDB, 2HKA chain A (NPC2 from cow) and 3WEA chain A (NPC2 from the Japanese carpenter ant) had the highest and second-highest structural similarity to C. sinensis NPC2-like models (representing 44 and 5 of all structures predicted, respectively; Additional file 1: Table S2). The most commonly predicted ligand (73%) was cholesterol (C3S).

Phylogenetic relationships

The phylogenetic relationship among C. sinensis NPC2-like proteins was determined using aligned coding domains; the resulting tree (Fig. 2a) was annotated with experimental data (Fig. 2b-h). NPC2-like proteins of C. sinensis clustered in four well-supported groups (pp = 0.81–1.0) that contained 27 paired orthologues, with eight and four proteins being unique to the Korean and Chinese isolates, respectively (Fig. 2b). The Korean and Chinese C. sinensis proteins most similar to NPC2 proteins from SWISS-PROT clustered within group 2 (E-value 1.23E−20–2.91E−24), whereas C. sinensis sequences that were least similar to NPC2 proteins from SWISS-PROT were within group 1 (E-value 2.64E−10–4.65) or group 4 (E-value 6.12E−23–5.28). Interestingly, NPC2-like proteins were frequently encoded on the same genome scaffold (Fig. 2c). For example, four groups of two to six Korean NPC2-like proteins within group 4 were encoded on the same scaffold, whereas one and two groups of Korean proteins (n = 2–5) in group 3 and group 1, respectively, were encoded on the same scaffold.

Fig. 2
figure2

Phylogenetic relationships among curated NPC2 proteins identified in two genomes of Clonorchis sinensis (a). Curated SWISS-PROT NPC2 proteins and known PDB structures for NPC2-like proteins were used as an outgroup. For each C. sinensis gene, their predicted Chinese (red) and Korean (blue) orthologues are presented as pairs in the phylogenetic tree. For each gene, the phylogenetic grouping (b) and the encoding scaffold (c) are listed. For all encoded protein sequences, the presence of a signal peptide (d), homology to canonical SWISS-PROT NPC2 proteins (e) and the presence of a conserved Pfam ML domain (PF02221) (f) are indicated. For representative pairs of orthologous proteins, the predicted quality of the I-TASSER model (C-score; “confidence score”) and the best matched PDB structure is shown (g). Transcription levels for each NPC2-like protein-encoding gene in available developmental and tissue-specific RNA-Seq libraries are included as a heatmap (h)

Based on available C. sinensis transcriptomic sequence data, NPC2-like protein-encoding genes within group 2 showed evidence of transcription in all developmental stages (metacercaria, juvenile and adult) and adult tissues (testis, ovary, muscle and sucker) (Fig. 2h). Most NPC2-like protein-encoding genes in C. sinensis showed moderate to high transcription in the metacercarial stage. Genes in group 1 showed the lowest (overall) transcription levels across all stages and tissues. The gene Cs-k2.gene7762 (Korean isolate) showed the highest transcription overall, and was highly transcribed in all stages and tissues studied.

Employing curated data sets (Fig. 2 and Additional file 1: Table S2), the paired orthologues Cs-k2.gene6404 and csin102672 were inferred to encode canonical NPC2 proteins. Transcriptomic evidence supported their constitutive transcription in all stages and in adult tissues. In addition, the proteins encoded by these genes were the most similar to curated NPC2 proteins from SWISS-PROT, and their predicted tertiary structures had the highest C-scores.

Evidence for structural conservation

Predicted high-confidence (C-score: ≥ 1) structures of 21 C. sinensis NPC2 and NPC2-like proteins (11 and 10 for Korean and Chinese isolates, respectively) were aligned with the two most similar PDB reference structures (2HKA and 3WEA) to assess conservation (Fig. 3). Most proteins were predicted to retain a conserved Ig-like β-sandwich fold conformation with seven-stranded β-sandwich folds fixed by three disulfide bonds (Cys-8-Cys-121, Cys-23-Cys-28, and Cys-74-Cys-80) and a large cavity in the interior of a protein barrel (Fig. 3a–c). In mammals, NPC2 binds cholesterol in the deep hydrophobic tunnel created by the βa and βb-βc loops (Fig. 3c) [19]. Predicted C. sinensis NPC2 and NPC2-like models were more similar to the resolved structure in the absence of bound cholesterol sulphate (2HKA chain A; Fig. 3b, c) than the resolved structure with an open pocket in the presence of bound cholesterol sulphate (2HKA chain C; Fig. 3d). Importantly, only Cs-k2.gene6404 and csin102672 proteins retained the three amino acid residues (Val-105, Tyr-109 and Phe-73) that are required for cholesterol binding [19] (Fig. 3a, c), further supporting their annotation as canonical NPC2 proteins.

Fig. 3
figure3

Assessment of structural conservation in Clonorchis sinensis NPC2 and NPC2-like proteins. Predicted structures of 21 C. sinensis NPC2 and NPC2-like proteins were aligned to assess conservation relative to two reference NPC2 protein structures (2HKA and 3WEA). a Alignment of NPC2 and NPC2-like sequences and conserved barrel with seven-stranded β-sandwich folds (shown in green) fixed by three disulfide bonds (Cys-8-Cys-121, Cys-23-Cys-28, and Cys-74-Cys-80; indicated by black lines below the alignment). b Positioning of loops in C. sinensis NPC2 models was modelled on a closed sterol binding pocket. Conserved residues are shown in red, variable residues are shown in blue and fewer variable residues are shown in white. c Hydrophobic tunnel from an opening created by three β-sandwich loops and highlighting three amino acid residues important for binding cholesterol (Val-105, i; Tyr-109, ii; and Phe-73, iii). d A fully open pocket in the presence of bound cholesterol sulphate (orange)

Discussion

This study conducted comparative genomic analyses to explore the substantial expansion of a NPC2-like protein family in C. sinensis. Using a bioinformatic workflow, we curated 35 and 32 C. sinensis NPC2-like proteins representing the Korean and Chinese isolates, respectively, clustered them into four distinct phylogenetic groups, predicted their tertiary structures and recorded transcriptional levels for the genes encoding these proteins in distinct developmental stages and tissues. These data provide evidence for the presence of a structurally conserved canonical form of NPC2 in C. sinensis and were used to explore the functional implications of genetic variation among members of the C. sinensis NPC2-like protein family.

Earlier characterisations of human NPC2 homologues [20] in liver flukes revealed an expansion of this protein family [16,17,18]. The present study predicted a higher number of NPC2-like proteins in C. sinensis compared with earlier studies and provides evidence of relatively high levels of genetic conservation between the Korean and Chinese isolates. To date, investigations of other platyhelminths or most other eukaryotes have reported only one or two copies of NPC2-like proteins [20], which likely retain a conserved role in sterol transport from the late endosome and/or lysosomes, together with conserved NPC1 proteins (C. sinensis gene identifiers Cs-k2.gene5262 and csin107525) [49]. The genetic mechanisms in C. sinensis that have led to a relatively rapid expansion and genetic diversification beyond a single, canonical NPC2-like gene remain to be investigated in detail. Clustering of C. sinensis NCP2-like proteins on the same draft genomic scaffolds (for each isolate) suggests gene duplication mechanisms are likely to be involved, including unequal crossing-over [50], retrotransposition [51], duplicated DNA transposition [52] and/or polyploidisation [53]. In arthropods, gene duplication events appear to have led to a similar expansion of ML domain-containing proteins (which includes NPC2-like proteins). For example, MD-2 or NPC2-like protein family expansions in insects have been recorded in Anopheles gambiae (13 copies), Aedes aegypti (15 copies), Tribolium castaneum (8 copies) and Drosophila melanogaster (8 copies) [20, 21]. The extent of NPC2 gene duplication events in other flatworms remains to be determined. Further characterisation of NPC2-like proteins encoded in all available genomes of all members of the phylum Platyhelminthes should provide useful insights into the evolution of this family of proteins as well as their functional roles in free-living and parasitic taxa.

Molecular characterisations of several arthropod ML-proteins support neofunctionalisation arising from gene duplication, with diversified ML proteins playing crucial roles in steroid biosynthesis [54], immunity [55] and chemoreception [23,24,25,26]. For example, arthropod ML proteins can act as receptors (e.g. “pattern recognition receptors”) or co-receptors for various ligands to modulate innate immune signalling pathways [21]. In addition, several NPC2-like proteins are highly expressed in the chemosensory organs of ants [43], ticks [24] and spiders [56], where they are reported to play a key role in chemoreception, by acting as carriers of semiochemicals [23,24,25]. As C. sinensis is taxonomically and evolutionarily distinct from arthropods, the gene expansion events in arthropods and liver flukes appear to be independent as they do not share common ancestry. Therefore, it is unlikely that the function of C. sinensis NPC2-like proteins can be inferred from amino acid sequence homology. However, the predominant transcription of most NPC2-like C. sinensis proteins in the metacercarial stage does provide support for a role in chemoreception; the infective metacercarial stage is exposed to a hostile environment which requires molecular mechanisms facilitating survival, adaptation, migration and development [57], and there is increasing evidence that bile stimulates the expression of particular genes involved in these processes [58]. Importantly, chemoreception has been shown to play a criticial role in the ability of newly excysted juveniles to locate the ampulla of Vater and migrate into the biliary system [15, 59]. The abundance of transcripts of most NPC2-like protein genes in the metacercarial stage might be indicative of an important role for these proteins during excystation, initial growth and development and/or migration. Whether these NPC2-like proteins are expressed and/or linked to the chemotactic behaviour of C. sinensis warrants detailed investigation. In addition to roles in chemotaxis, a high level of transcription of some of the NPC2-encoding genes in adult tissues (including reproductive tissues and sucker) might suggest broader roles in lipid metabolism, feeding and/or reproduction [16, 20]. In other parasites, including intracellular protists [60, 61] and schistosomes [62], which cannot synthesise cholesterol, these molecules may have functions in nutrient uptake, immune evasion and/or energy storage. If C. sinensis is unable to synthesise cholesterol, as is the case for the related liver fluke O. viverrini [16], canonical NPC2 in C. sinensis might assume similar functions.

In this study, our established bioinformatics workflow [63] assisted in the prediction of a conserved set (family) of C. sinensis NPC2-like proteins. Evidence that two pairs of NPC2-like protein orthologues (Cs-k2.gene11598/csin11538 and Cs-k2.gene9726/csin101111) undergo positive selection [17] lends support for a recent expansion of this protein family. Interestingly, a small number of gene encoding these proteins (n = 4–8) were present exclusively in the Korean or Chinese isolate; whether these genes/proteins have evolved recently and are indeed isolate-specific remains to be established. Resolving the final copy numbers of NPC2-like protein genes in C. sinensis will necessitate the sequencing and assembly of complete, chromosome-contiguous genomes using ‘third-generation’ technologies [64]. This would also pave the way to detailed comparative genome analyses.

Conclusions

In the present study, we used comparative genomics and functional annotation tools to explore a novel family of NPC2-like proteins in C. sinensis. Based on our findings, we propose that these proteins are involved in binding and transportation of sterols and/or other lipids throughout the life-cycle of this parasite. The curated set of these proteins presented herein provides a solid foundation for future investigations of the molecular functions of NPC2-like proteins in C. sinensis and other liver flukes.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its additional file.

Abbreviations

AIC:

Akaike information criteria

BI:

Bayesian inference

BLAST:

basic local alignment tool

BLASTn:

nucleotide BLAST

BLASTp:

protein BLAST

BLAT:

BLAST-like alignment tool

C-score:

I-TASSER model confidence score

CCA:

cholangiocarcinoma

CPM:

counts per million

ENA:

European Nucleotide Archive

GFF:

general feature format

I-TASSER:

iterative threading assembly refinement

IARC:

International Agency for Research on Cancer

LBPs:

lipid-binding proteins

MCMC:

Markov chain Monte Carlo

MD-2:

myeloid differentiation factor-2

ML:

MD-2-related lipid-binding domain

NCBI-nr:

National Center for Biotechnology Information nonredundant database

NPC2:

Niemann-pick type C2

NTD:

neglected tropical disease

PDB:

protein data bank

pp:

posterior probability

PSRF:

potential scale reduction factor

RMSD:

root-mean-square deviation of atomic positions

RNA-Seq:

RNA sequencing

SRA:

sequence read archive

TM-score:

I-TASSER topology modelling score

References

  1. 1.

    WHO. International statistical classification of diseases and related health problems. WHO International Classification of Diseases. Geneva: World Health Organization; 2004. https://www.who.int/classifications/icd/icdonlineversions/en/. Accessed 28 Nov 2019.

  2. 2.

    Khan S, Toledano M, Taylor-Robinson S. Epidemiology, risk factors, and pathogenesis of cholangiocarcinoma. HBP. 2008;10:77–82.

  3. 3.

    Fürst T, Keiser J, Utzinger J. Global burden of human food-borne trematodiasis: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12:210–21.

  4. 4.

    Keiser J, Utzinger J. Emerging foodborne trematodiasis. Emerg Infect Dis. 2005;11:1507–14.

  5. 5.

    Sripa B, Kaewkes S, Sithithaworn P, Mairiang E, Laha T, Smout M, et al. Liver fluke induces cholangiocarcinoma. PLoS Med. 2007;4:e201.

  6. 6.

    Qian MB, Utzinger J, Keiser J, Zhou XN. Clonorchiasis. The Lancet. 2016;387:800–10.

  7. 7.

    Xia J, Jiang SC, Peng HJ. Association between liver fluke infection and hepatobiliary pathological changes: a systematic review and meta-analysis. PLoS One. 2015;10:e0132673.

  8. 8.

    Lun ZR, Gasser RB, Lai DH, Li AX, Zhu XQ, Yu XB, et al. Clonorchiasis: a key foodborne zoonosis in China. Lancet Infect Dis. 2005;5:31–41.

  9. 9.

    Young ND, Campbell BE, Hall RS, Jex AR, Cantacessi C, Laha T, et al. Unlocking the transcriptomes of two carcinogenic parasites, Clonorchis sinensis and Opisthorchis viverrini. PLoS Negl Trop Dis. 2010;4:e719.

  10. 10.

    Zhang J, Chiodini R, Badr A, Zhang G. The impact of next-generation sequencing on genomics. J Genet Genomics. 2011;38:95–109.

  11. 11.

    Won J, Ju JW, Kim SM, Shin Y, Chung S, Pak JH. Clonorchis sinensis infestation promotes three-dimensional aggregation and invasion of cholangiocarcinoma cells. PloS One. 2014;9:e110705.

  12. 12.

    Kaewkes S. Taxonomy and biology of liver flukes. Acta Trop. 2003;88:177–86.

  13. 13.

    Kim HG, Han J, Kim MH, Cho KH, Shin IH, Kim GH, et al. Prevalence of clonorchiasis in patients with gastrointestinal disease: a Korean nationwide multicenter survey. World J Gastroenterol. 2009;15:86–94.

  14. 14.

    Watanapa P, Watanapa W. Liver fluke-associated cholangiocarcinoma. Br J Surg. 2002;89:962–70.

  15. 15.

    Li S, Im Kim T, Yoo WG, Cho PY, Kim TS, Hong SJ. Bile components and amino acids affect survival of the newly excysted juvenile Clonorchis sinensis in maintaining media. Parasitol Res. 2008;103:1019–24.

  16. 16.

    Young ND, Nagarajan N, Lin SJ, Korhonen PK, Jex AR, Hall RS, et al. The Opisthorchis viverrini genome provides insights into life in the bile duct. Nat Commun. 2014;5:4378.

  17. 17.

    Wang D, Korhonen PK, Gasser RB, Young ND. Improved genomic resources and new bioinformatic workflow for the carcinogenic parasite Clonorchis sinensis: biotechnological implications. Biotechnol Adv. 2018;36:894–904.

  18. 18.

    Wang X, Chen W, Huang Y, Sun J, Men J, Liu H, et al. The draft genome of the carcinogenic human liver fluke Clonorchis sinensis. Genome Biol. 2011;12:R107.

  19. 19.

    Xu S, Benoff B, Liou HL, Lobel P, Stock AM. Structural basis of sterol binding by NPC2, a lysosomal protein deficient in Niemann-Pick type C2 disease. J Biol Chem. 2007;282:23525–31.

  20. 20.

    Inohara N, Nuñez G. ML - a conserved domain involved in innate immunity and lipid metabolism. Trends Biochem Sci. 2002;27:219–21.

  21. 21.

    Shi XZ, Zhong X, Yu XQ. Drosophila melanogaster NPC2 proteins bind bacterial cell wall components and may function in immune signal pathways. Insect Biochem Mol Biol. 2012;42:545–56.

  22. 22.

    Storch J, Xu Z. Niemann-Pick C2 (NPC2) and intracellular cholesterol trafficking. Biochim Biophys Acta. 2009;1791:671–8.

  23. 23.

    Pelosi P, Iovinella I, Felicioli A, Dani FR. Soluble proteins of chemical communication: an overview across arthropods. Front Physiol. 2014;5:320.

  24. 24.

    Iovinella I, Ban L, Song L, Pelosi P, Dani FR. Proteomic analysis of castor bean tick Ixodes ricinus: a focus on chemosensory organs. Insect Biochem Mol Biol. 2016;78:58–68.

  25. 25.

    Renthal R, Manghnani L, Bernal S, Qu Y, Griffith WP, Lohmeyer K, et al. The chemosensory appendage proteome of Amblyomma americanum (Acari: Ixodidae) reveals putative odorant-binding and other chemoreception-related proteins. Insect Sci. 2017;24:730–42.

  26. 26.

    Bryant CE, Spring DR, Gangloff M, Gay NJ. The molecular basis of the host response to lipopolysaccharide. Nat Rev Microbiol. 2010;8:8–14.

  27. 27.

    Howe KL, Bolt BJ, Shafie M, Kersey P, Berriman M. WormBase ParaSite - a comprehensive resource for helminth genomics. Mol Biochem Parasitol. 2017;215:2–10.

  28. 28.

    Bairoch A, Apweiler R. The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999. Nucleic Acids Res. 1999;27:49–54.

  29. 29.

    Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.

  30. 30.

    Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5.

  31. 31.

    Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.

  32. 32.

    Sonnhammer EL, Eddy SR, Durbin R. Pfam: a comprehensive database of protein domain families based on seed alignments. Proteins. 1997;28:405–20.

  33. 33.

    Khan A, Mathelier A. Intervene: a tool for intersection and visualization of multiple gene or genomic region sets. BMC Bioinform. 2017;18:287.

  34. 34.

    Yu G, Smith DK, Zhu H, Guan Y, Lam TTY. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.

  35. 35.

    GFF/GTF utility providing format conversions, region filtering, FASTA sequence extraction and more. https://github.com/gpertea/gffread. Accessed 28 Nov 2019.

  36. 36.

    Kent WJ. BLAT - the BLAST-like alignment tool. Genome Res. 2002;12:656–64.

  37. 37.

    Slater GSC, Birney E. Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005;6:31.

  38. 38.

    Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.

  39. 39.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2013;30:923–30.

  40. 40.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

  41. 41.

    Armenteros JJA, Tsirigos KD, Sonderby CK, Petersen TN, Winther O, Brunak S, et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37:420–3.

  42. 42.

    Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.

  43. 43.

    Ishida Y, Tsuchiya W, Fujii T, Fujimoto Z, Miyazawa M, Ishibashi J, et al. Niemann-Pick type C2 protein mediating chemical communication in the worker ant. Proc Natl Acad Sci USA. 2014;111:3847–52.

  44. 44.

    Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera - a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.

  45. 45.

    Ranwez V, Douzery EJ, Cambon C, Chantret N, Delsuc F. MACSE v2: toolkit for the alignment of coding sequences accounting for frameshifts and stop codons. Mol Biol Evol. 2018;35:2582–4.

  46. 46.

    Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.

  47. 47.

    Kalyaanamoorthy S, Minh BQ, Wong TK, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587.

  48. 48.

    Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.

  49. 49.

    Subramanian K, Balch WE. NPC1/NPC2 function as a tag team duo to mobilize cholesterol. Proc Natl Acad Sci USA. 2008;105:15223–4.

  50. 50.

    Zhang J. Evolution by gene duplication: an update. Trends Ecol Evol. 2003;18:292–8.

  51. 51.

    Brosius J. Retroposons - seeds of evolution. Science. 1991;251:753.

  52. 52.

    Bailey JA, Liu G, Eichler EE. An Alu transposition model for the origin and expansion of human segmental duplications. Am J Hum Genet. 2003;73:823–34.

  53. 53.

    Ohno S. Gene duplication and the uniqueness of vertebrate genomes circa 1970–1999. Semin Cell Dev Biol. 1999;10:517–22.

  54. 54.

    Huang X, Warren JT, Buchanan J, Gilbert LI, Scott MP. Drosophila Niemann-Pick type C-2 genes control sterol homeostasis and steroid biosynthesis: a model of human neurodegenerative disease. Development. 2007;134:3733–42.

  55. 55.

    Venugopal PG, Nutman TB, Semnani RT. Activation and regulation of toll-like receptors (TLRs) by helminth parasites. Immunol Res. 2009;43:252–63.

  56. 56.

    Vizueta J, Frías-López C, Macías-Hernández N, Arnedo MA, Sánchez-Gracia A, Rozas J. Evolution of chemosensory gene families in arthropods: insight from the first inclusive comparative transcriptome analysis across spider appendages. Genome Biol Evol. 2017;9:178–96.

  57. 57.

    Lu Y, Yoo WG, Dai F, Lee JY, Pak JH, Sohn WM, et al. Characterization of a novel organic solute transporter homologue from Clonorchis sinensis. PLoS Negl Trop Dis. 2018;12:e0006459.

  58. 58.

    Kim TI, Cho PY, Yoo WG, Li S, Hong SJ. Bile-induced genes in Clonorchis sinensis metacercariae. Parasitol Res. 2008;103:1377–82.

  59. 59.

    Kim TI, Yoo WG, Kwak BK, Seok JW, Hong SJ. Tracing of the bile-chemotactic migration of juvenile Clonorchis sinensis in rabbits by PET-CT. PLoS Negl Trop Dis. 2011;5:e1414.

  60. 60.

    Bolaños J, Betanzos A, Javier-Reyna R, García-Rivera G, Huerta M, Pais-Morales J, et al. EhNPC1 and EhNPC2 proteins participate in trafficking of exogenous cholesterol in Entamoeba histolytica trophozoites: relevance for phagocytosis. PLoS Pathog. 2016;12:e1006089.

  61. 61.

    Semini G, Paape D, Paterou A, Schroeder J, Barrios-Llerena M, Aebischer T. Changes to cholesterol trafficking in macrophages by Leishmania parasites infection. MicrobiologyOpen. 2017;6:e469.

  62. 62.

    Nawaratna SS, Gobert GN, Willis C, Mulvenna J, Hofmann A, McManus DP, Jones MK. Lysosome-associated membrane glycoprotein (LAMP) - preliminary study on a hidden antigen target for vaccination against schistosomiasis. Sci Rep. 2015;5:15069.

  63. 63.

    Stroehlein AJ, Young ND, Gasser RB. Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups. Sci Rep. 2018;8:6808.

  64. 64.

    Korhonen PK, Young ND, Gasser RB. Making sense of genomes of parasitic worms: tackling bioinformatic challenges. Biotechnol Adv. 2016;34:663–86.

Download references

Acknowledgements

Not applicable.

Funding

This research was funded by grants from the National Health and Medical Research Council (NHMRC) and the Australian Research Council (ARC) to RBG and NDY. Other support from Yourgene Bioscience is gratefully acknowledged (RBG). NDY is supported by a NHMRC Career Development Fellowship.

Author information

Conceived and designed the study and supervised the project: NDY, AJS, BCHC and RBG. Undertook the study and data analysis: MA, AJS, RSH, NDY and RBG. Contributed to the interpretation of findings and writing of the manuscript: MA, NDY, AJS, BCHC and RBG. All authors read and approved the final manuscript.

Correspondence to Robin B. Gasser or Neil D. Young.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Anari, M., Stroehlein, A.J., Hall, R.S. et al. Expanded complement of Niemann-Pick type C2-like protein genes in Clonorchis sinensis suggests functions beyond sterol binding and transport. Parasites Vectors 13, 38 (2020). https://doi.org/10.1186/s13071-020-3910-0

Download citation

Keywords

  • Clonorchis sinensis
  • Comparative genomics
  • Niemann-pick type C2
  • NPC2
  • Functional protein annotation
  • Adaptation

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.