Skip to main content

Major SCP/TAPS protein expansion in Lucilia cuprina is associated with novel tandem array organisation and domain architecture

Abstract

Background

Larvae of the Australian sheep blowfly, Lucilia cuprina, parasitise sheep by feeding on skin excretions, dermal tissue and blood, causing severe damage known as flystrike or myiasis. Recent advances in -omic technologies and bioinformatic data analyses have led to a greater understanding of blowfly biology and should allow the identification of protein families involved in host-parasite interactions and disease. Current literature suggests that proteins of the SCP (Sperm-Coating Protein)/TAPS (Tpx-1/Ag5/PR-1/Sc7) (SCP/TAPS) superfamily play key roles in immune modulation, cross-talk between parasite and host as well as developmental and reproductive processes in parasites.

Methods

Here, we employed a bioinformatics workflow to curate the SCP/TAPS protein gene family in L. cuprina. Protein sequence, the presence and number of conserved CAP-domains and phylogeny were used to group identified SCP/TAPS proteins; these were compared to those found in Drosophila melanogaster to make functional predictions. In addition, transcription levels of SCP/TAPS protein-encoding genes were explored in different developmental stages.

Results

A total of 27 genes were identified as belonging to the SCP/TAPS gene family: encoding 26 single-domain proteins each with a single CAP domain and a solitary double-domain protein containing two conserved cysteine-rich secretory protein/antigen 5/pathogenesis related-1 (CAP) domains. Surprisingly, 16 SCP/TAPS predicted proteins formed an extended tandem array spanning a 53 kb region of one genomic region, which was confirmed by MinION long-read sequencing. RNA-seq data indicated that these 16 genes are highly transcribed in all developmental stages (excluding the embryo).

Conclusions

Future work should assess the potential of selected SCP/TAPS proteins as novel targets for the control of L. cuprina and related parasitic flies of major socioeconomic importance.

Graphical Abstract

Background

Flystrike (cutaneous myiasis), caused by the Australian sheep blowfly, Lucilia cuprina, is a disease of major economic significance in livestock animals worldwide [1,2,3,4,5,6]. Traditionally, the control of this disease has relied heavily on insecticide use, management procedures (e.g., surgical mulesing) or a combination of these measures [7]. However, due to resistance against most insecticides in current use [8,9,10] and pressure from animal welfare groups to eliminate surgical mulesing [11], there is a need for an alternative approach to treatment or control. A possible pathway to new intervention strategies could be to interfere with biological processes or pathways in L. cuprina and/or disrupting the parasitic relationship between the fly and its host.

Available genomic resources for L. cuprina [12] provide an excellent opportunity to explore biological processes and pathways in the fly as well as host-fly interactions and the pathogenesis of flystrike at the molecular level. Indeed, previous work [12] eluded to a number of groups of molecules, such as GTPases, gustatory and odorant receptors and serine proteases, some of which could be involved in parasitism (the fly-host interplay) and might be potential intervention targets. One such group is the cysteine-rich secretory protein/antigen 5/pathogenesis related-1 (CAP) or sperm-coating protein/TAPS (Tpx-1/Ag5/PR-1/Sc7) (SCP/TAP) protein superfamily, typified by a conserved CAP domain (Pfam: PF00188, InterProScan: IPR014044) [13].

Most SCP/TAPS proteins characterised to date are excretory/secretory and appear to function extracellularly [14, 15]. Many are implicated in male reproductive processes [16, 17], host-parasite interactions [18,19,20,21] and host immune modulation [22,23,24,25]. SCP/TAPS proteins are also components of venoms produced by wasps, ants and some reptiles and are present in the saliva of some blood-feeding insects [26,27,28,29]. Although the functions of many such proteins are enigmatic, some are proposed to be proteases, protease inhibitors, signaling molecules, ion channel regulators or membrane disrupters [14]. The best characterised insect, Drosophila melanogaster (vinegar fly), has 26 SCP/TAPS proteins (representing two groups), 15 of which are preferentially expressed in adult males and proposed to be involved in male reproduction [16]. Despite this information, almost nothing is known about SCP/TAPS proteins in L. cuprina. In the present study, we identified the full complement of SCP/TAPS genes and inferred proteins in L. cuprina, studied the structure and organisation of these genes and explored their transcription profiles in different stages/sexes of the fly as a foundation for future investigations of their functional roles and possible potential as intervention targets.

Methods

Annotation and 3D modelling of SCP/TAPS proteins

Total RNA was isolated separately from eggs, mixed larval stages (1st, 2nd and 3rd instars), an adult male and an adult female of L. cuprina and sequenced using established methods [12]. RNA-sequencing (RNA-seq) data were mapped to an enhanced version of L. cuprina genomic scaffolds (GenBank accession no. JRES00000000.1) using the program TopHat version 2.0.14 [30] and assembled into transcripts using Cufflinks version 2.2.1 [31], and SCP/TAPS protein-coding genes were identified using TransDecoder [32] (https://github.com/TransDecoder/TransDecoder/wiki). Transcripts with a conserved CAP domain (Pfam ID "PF00188" and/or InterProScan ID "IPR014044") were identified using InterProScan 5 [33] and aligned to genomic scaffolds using BLAT version 35.1 [34] employing the “-fine” parameter. Additionally, GeneMark version 4.2.9 [35] employing the “-ES” parameter was used to infer gene models for the eight scaffolds that contained BLAT hits. To infer gene models using Exonerate version 2.2.0 [36], the positions of BLAT hits were first extracted to restrict gene prediction to candidate gene regions and then extended up- and down-stream to predict complete, full-length gene models. These regions were analysed employing the “Coding2genome:bestfit” model. Finally, MAKER3 predictions (Anstead et al., unpublished), BLAT hits, gene predictions inferred using GeneMark and Exonerate, genomic scaffolds and RNA-seq transcripts were displayed and manually inspected in the Integrative Genomics Viewer (IGV; version 2.4.3) [37]. Manual curation of SCP/TAPS coding regions was performed by analysing each gene prediction in IGV and determining the optimal gene model based on the following criteria: (1) each gene model had to have an open reading frame (ORF) with unambiguous intron-exon boundaries supported by RNA-seq expression data; (2) coding regions could not be smaller than 400 bp, in accordance with the minimal length of the CAP domain (IPRO14044; 150 amino acids); (3) a minimum of five RNA-seq reads had to be mapped to the coding region. Alternative transcripts were not included in this analysis. Three dimensional (3D) structures of individual curated amino acid sequences were modelled using the program I-TASSER [38] and displayed using PyMOL (The PyMOL Molecular Graphics System, version 1.2r3pre, Schrödinger, LLC.)

Confirmation of scaffold assembly

The assembly of selected scaffolds was confirmed by mapping data from 5 kb and 3 kb mate pair libraries and 100 bp paired-end libraries using Bowtie2 v.2.1.0 [39], and read-mapping coverage was calculated using BEDtools [40] and deepTools [41]. Next, high-molecular-weight genomic DNA (gDNA) was extracted from one adult male L. cuprina, which was digested in CTAB buffer containing proteinase K (20 mg/ml, Clontech, Hilden, Germany) and incubated for 24 h at 37 °C, followed by standard phenol/chloroform extraction [42]. Samples were not vortexed, and wide-bore pipette tips were used to reduce the chance of shearing DNA. Genomic DNA was sequenced using the MinION Mk1B SQK-LSK108 1D gDNA kit for long reads (> 100 kb) without size selection. Base calling was performed using Albacore version 2.2.7 (Oxford Nanopore) (https://github.com/zhaoc1/nanoflow) within the MinKNOW software. Reads that passed base calling quality control were aligned to genomic scaffolds using GraphMap [43] and Minimap2 [44], and the consensus sequences were analysed using NCBI motif finders [45] to identify CAP domains. Intronic regions of all 16 CAP genes located on the array were analysed for duplicated intronic sequences using BLAST [46] and MUSCLE [47] to confirm the integrity of the scaffold assembly. Finally, ten arthropod genomes (Aedes aegypti, Ae. albopictus, Anopheles albimanus, Glossina austeni, G. brevipalpis, G. fuscipes, G. pallidipes, G. palpalis, Musca domestica and Stomoxys calcitrans) were obtained from VectorBase (https://vectorbase.org/) and the Drosophila melanogaster genome from FlyBase [48] and were interrogated to identify SCP/TAPS tandem array genes. All listed genes containing the CAP domain (IPR014044) were downloaded (7 March 2019) and gene position data were analysed to determine tandem array SCP/TAPS proteins.

Transcription profiling

The transcription profiles for SCP/TAPS genes were determined using RSEM [49] employing stringent criteria. Data from paired-end RNA-seq libraries (i.e., female, male and mixed larvae) were mapped using TopHat, employing stringent settings of 100% identity and not allowing multi-mapping and insertion/deletion events, but permitting gaps. These conditions were used to ensure accurate measurements of gene expression as mapping using standard parameters can make it difficult to distinguish between highly similar genes. RSEM transcription levels were normalised and represented as TPM (transcripts per million) per gene across the three RNA-seq libraries.

Phylogenetic analysis of SCP/TAPS proteins

Phylogenetic analysis of D. melanogaster [48] and L. cuprina SCP/TAPS amino acid sequence data was performed using Geneious v.7.1.9 [50]. SCP/TAPS amino acid sequences were aligned using MUSCLE [47]. Alignment positions representing a gap in more than 50% of all aligned sequenced were manually removed. Bayesian inference trees were built using the program MrBayes v.3.2.2 [51] and rooted using a Caenorhabditis elegans SCP/TAPS gene (NP_001256323.1) as an outgroup. MrBayes parameters were set, as described previously [52]. Trees were displayed and modified using the programs FigTree v.1.4.1 (https://tree.bio.ed.ac.uk/software/figtree/) and Inkscape (https://www.inkscape.org/en/).

Molecular amplification of selected SCP/TAPS genes

Total RNA was extracted from a laboratory strain (LS) of L. cuprina (provided by Dr Peter James, University of Queensland, Australia) from all life stages (i.e., egg, all three larval stages and both male and female adults). Total RNA was extracted using the TRIsure RNA extraction protocol (Bioline Pty Ltd, New South Wales, Australia) according to the manufacturer’s instructions and reverse transcription was performed using 2 μg of total RNA per sample using GoScript Reverse Transcription System (Promega). Primers (Additional file 1: Table S1) were designed manually or using the program Geneious v.7.1.9 and 3’ and 5’ RACE protocols (Clontech Laboratories Inc., Palo Alto, CA, USA) were employed in accordance with the manufacturer’s directions to amplify selected proteins.

Results

Lucilia cuprina has an expanded SCP/TAPS gene set representing three distinct groups

A total of 75 transcripts containing a conserved CAP domain were identified and mapped to six genomic scaffolds (Additional file 2: Data S1). Based on these mapped regions, 27 genes belonging to the SCP/TAPS family were predicted (Additional file 3: Data S2). Curation of the inferred gene models identified one gene coding for a solitary double-domain (DD) protein (denoted as LCSc0-DD1) with two conserved CAP domains and 26 genes coding for single-domain (SD) proteins with a solitary CAP domain (Additional file 4: Data S3). The majority (81%) of these SD proteins displayed a moderate to high degree of amino acid similarity within the CAP domain (65–100%), whereas the DD protein displayed moderate ( = 47%) sequence similarity to the SD proteins (Fig. 1). The proximal CAP domain aligned poorly with SD proteins (below the 30% threshold) and is therefore not displayed in Fig. 1. All 26 SD proteins contained conserved cysteine residues critical for the formation of disulphide bonds, a feature common to those proteins belonging to the SCP/TAPS family [15]. The majority (approximately 65%) had cysteine residues located in conserved positions (C1–C2 between 20 and 30, C3–C4 between 106–117 and C5–C6 between 224 and 233), while the remaining proteins displayed a greater variety in cysteine residue position.

Fig. 1
figure 1

Alignment of 27 L. cuprina SCP/TAPS. The 27 SCP/TAPS identified in L. cuprina were aligned using MUSCLE (47) and manually filtered to remove alignment positions with > 50% gap frequency. Conserved cysteine residues are highlighted in yellow and indicated by an arrow. Black dots represent conserved residues common to the CAP domain. Red bars below the sequence alignment indicate CAP domain location. The sequence logo indicates the prevalence of an amino acid at each position based on the size of the letter, and numbers indicate the amino acid position in the alignment. Residues that are highly conserved (≥ 90%) are highlighted based on biochemical properties (see Geneious version 7.1.9). Regions with low (< 30%) similarity (such as the proximal CAP domain of LCSc0-DD1) are not displayed here

Phylogenetic analyses of the 27 SCP/TAPS proteins suggested three main clades: those encoded by genes on scaffold xfSc0000049, those on scaffold Sc0000031 and those on scaffold Sc0000002 (Fig. 4). SD protein LCSc8-SD4 did not group with any other SD proteins, while DD protein LCSc0-DD1 did not cluster with any of the 26 SD proteins (Fig. 4). Phylogenetic analyses comparing SCP/TAPS proteins found in D. melanogaster and L. cuprina demonstrated that the 16 tandem array genes formed a clade with 2 D. melanogaster proteins Agr and Agr2 (Fig. 4), and syntenic analysis of the associated genes revealed the presence of the same gene (rutabaga; FlyBase ID: FBgn0003301) upstream of both the 16 tandem array genes in L. cuprina and the Agr and Agr2 genes in D. melanogaster. Furthermore, three L. cuprina genes had single-copy orthologs in D. melanogaster (i.e., LCSc8-SD4 to CG43775, LCSc19-SD5 to CG8483 and LCSc31-SD10 to CG9400).

Most SCP/TAPS protein-encoding genes are organised in a large tandem array

Sixteen SCP/TAPS genes of L. cuprina formed a tandem array of 53 kb on scaffold xfSc0000049 (Fig. 2), with all the encoded proteins (n = 16) displaying a moderate to high degree (60–95%) of amino acid similarity upon pairwise comparison. Additionally, the tandem array proteins possessed a unique genomic architecture in that each of the 16 transcripts encoded ~250 amino acid residues spanning three similarly sized exons and intronic domains of similar length, except for LCSc49-SD13 and LCSc49-SD22 (Additional file 3: Data S2). The CAP domain was conserved among proteins encoded by these 16 genes relative to other L. cuprina SCP/TAPS proteins exhibiting 76%–100% amino acid similarity upon pairwise comparison (Fig. 1). All SCP/TAPS in this array contained the CRISP motif [14] with up to ten conserved cysteine residues required for disulphide bond formation [15]. These residues were positioned in similar locations along the transcript, which is consistent with the location of these residues found in other species [53]. Alignment of L. cuprina RNA-seq libraries (i.e., female, male and mixed larvae) indicated that each gene was transcribed. Stringent mapping of transcriptomic reads to the entire scaffold containing the putative tandem array demonstrated that each gene was supported by RNA-seq data. All genes had at least ten RNA-seq reads mapped across the entire coding domain. Predicted protein models verified the core structural elements of the SCP/TAPS proteins (i.e., the α-β-α sandwich and stabilizing disulphide bonds). These structural elements are largely comprised of conserved CAP residues and appear to be internal while variable sequence regions more commonly coincide with amino acids located on the protein surface (Additional file 5: Figure S1, Additional file 6: Data S4).

Fig. 2
figure 2

Representation of the SCP/TAPS tandem array in L. cuprina. A 53 kb region of scaffold xfSc0000049 containing a tandem array of 16 SCP/TAPS genes. Genes (black arrows with transcription orientation indicated by arrow), Oxford Nanopore MinION mapped long reads (blue) with associated coverage plot (grey) are pictured. A mapped long read is shown depicted with a red outline highlighting that a single read spans the entire tandem array

When the data from the 5 kb and 3 kb mate-pair libraries and 100 bp paired-end reads (SUB3252906) were mapped to scaffold xfSc0000049, consistent coverage was observed across predicted coding regions. Long-read sequencing data using MinION sequencing (Oxford Nanopore Technologies, Oxford, UK) and subsequent read-mapping supported the integrity of scaffold xfSc0000049. One 61,501 bp read (ID: 102324e2-a088-4652-a4d2-7164e8fb4f52) mapped across the entire scaffold, providing additional support for this region (Fig. 2). Coverage analysis using the software package plotCoverage revealed a mean coverage of 12.4 reads and 80% of all positions had a coverage of greater than 9 reads across the scaffold xfSc0000049 region (xfSc0000049: 608647-664041). In addition, CIGAR string extraction showed 84.7% of all mapped nucleotide positions were matches. The consensus sequence shared 96.1% sequence identity with the reference scaffold and all 16 CAP domains were present in the consensuses sequences when analysed using NCBI motif finders. Finally, the analysis of 11 related arthropod genomes did not identify SCP/TAPS genes arranged in a tandem array to the extent identified in this study (Additional file 7: Table S2).

A novel double-domain CAP protein identified in L. cuprina

One of the 27 SCP/TAPS protein genes inferred here encoded a unique double-domain protein (LCSc0-DD1). Two independent gene prediction tools (i.e., GeneMark and Exonerate) inferred slightly different models for this unique gene. Of these models, only five exons were supported by RNA-seq data, suggesting the predicted gene model inferred by GeneMark was accurate (Additional file 8: Data S5). This gene model was curated manually to ensure that the appropriate splice junctions and open reading frame across the five predicted exons was evident. This model was fully supported by the MAKER3 prediction (Anstead 2020, unpublished) and BLAT-mapped transcripts. Exonerate predictions also supported the curated model; however, a putative downstream exon was also predicted in both cases (Fig. 3, Additional file 8: Data S5).

Fig. 3
figure 3

Representation of LCSc0-DD1 protein in L. cuprina genomic scaffolds. A predicted gene model (grey) consisting of five exons encoding a double-domain protein (LCSc0-DD1); supported by complementary data from the L. cuprina genome assembly, MAKER 3 predictions (blue), BLAT transcript mapping (red), gene prediction tools GeneMark and Exonerate (green and pink, respectively) and RNA-seq data from all life stages with associated coverage plot. Raw files supplied in Additional file 8: Data S5

Predicted gene models were verified employing conventional PCR strategies using primers flanking the conserved CAP domains. Transcripts containing both CAP domains were successfully amplified from a male cDNA library; however, a complete transcript was not amplified using conventional and/or RACE amplification strategies. The presence of two CAP domains in a single cDNA amplicon provides strong support for the presence of a double-domain protein in L. cuprina (Additional file 9: Data S6).

Diverse transcription profiles for SCP/TAPS protein-encoding genes

Although transcription profiles for the individual SCP/TAPS genes varied among libraries (Additional file 10: Table S3), all genes were transcribed in at least one of the RNA-seq libraries (TPM: 0.4 to 456.8). Eleven of the 16 genes (i.e., LCSc49-SD14, LCSc49-SD15, LCSc49-SD16, LCSc49-SD17, LCSc49-SD18, LCSc49-SD19, LCSc49-SD20, LCSc49-SD23, LCSc49-SD24, LCSc49-SD25 and LCSc49-SD26) in the tandem array were predominantly transcribed in mixed-larval samples, whereas two tandem-array genes (LCSc49-SD13 and LCSc49-SD22) were expressed in male and female RNA-seq libraries. Two genes (LCSc49-SD11 and LCSc49-SD12) in the tandem array were transcribed in all three RNA-seq libraries, and one gene (LCSc49-SD21) was specifically transcribed in the library representing adult males (Fig. 4). The TPM (transcripts per million) values for genes in the tandem array (Additional file 10: Table S3) were significantly higher (P < 0.05) (mean: 46.6) compared with all other SCP/TAPS genes (mean: 4.3), indicating higher transcription levels for genes within the tandem array relative to other SCP/TAPS protein genes in L. cuprina. Additionally, two genes (i.e., LCSc2-SD1 and LCSc31-SD9) that were not encoded in the array had TPM values > 10 (Additional file 10: Table S3). Furthermore, LCSc0-DD1 also displayed low levels of transcription across the three libraries (female: 0.2, male: 1.88, mixed larvae: 1.34) indicating that the double-domain protein was predominantly transcribed in adult male and mixed larval samples, albeit at low levels (Fig. 4).

Fig. 4
figure 4

Phylogenetic relationships among curated SCP/TAPS proteins. a Phylogenetic tree of 27 SCP/TAPS predicted protein sequences derived from L. cuprina. Colours indicate the scaffold from which the SCP/TAPS gene originated. b Heatmap displaying log10(TPM+1)-scaled transcription levels for each SCP/TAPS gene across the adult female, adult male and mixed larval libraries. Unadjusted TPM values can be obtained from Additional file 10: Table S3. c Phylogenetic tree of SCP/TAPS in both L. cuprina (red) and D. melanogaster (blue). d Indicates a Double Domain or Single Domain protein and e gene name of SCP/TAPS proteins from both species in the order in which they cluster. The curly bracket indicates genes that form the predicted tandem array and the square brackets show single-copy orthologs and f the encoding scaffold or chromosome. Trees were rooted using the C. elegans SCP/TAPS protein NP_001256323.1 as an outgroup and constructed using Bayesian inference (MrBayes v3.3.3) with branch labels indicating posterior probability values

Discussion

Here, we annotated and curated all SCP/TAPS protein-encoding genes in L. cuprina, classified them using a phylogenetic approach, established their genomic locations and organisation, assessed the transcription profiles of individual genes and explored the domain architectures of encoded proteins.

SCP/TAPS proteins share characteristic structural components that typify this protein superfamily [13]. One such feature is the α-β-α sandwich structure, which usually comprises ~50% of the protein structure [21] and is conserved among species. The variability in the remaining 50% of the protein structure suggests that the CAP domain forms a versatile molecular framework, enabling the diversification of SCP/TAPS protein function [21]. This proposal of functional diversification is supported, to some extent, by marked variation in transcription profiles between adult and larval stages and sexes of the fly; in adult flies, some of these molecules might associate with basic reproductive and developmental processes, as shown in Drosophila, whereas in larvae, others might assume roles in regulating protease activity, protease inhibition and/or a modulation of the host immune responses during larvae establishment and parasitism [19, 20]. Clearly, these are exciting areas worthy of pursuit using a combination of technological approaches, including immunoproteomics, transcriptomics and in vivo experiments in sheep.

Lucilia cuprina tandem array genes (TAGs) were variable, with predicted models suggesting that variable sequence regions (non-CAP residues) comprise much of the protein surface, while conserved residues form the internal protein structure (Additional file 5: Figure S1). It would be tempting to suggest that these proteins are under selective pressure, but future work is required to assess dN/dS ratios linked to tandem arrays of SCP/TAPS protein and other (e.g., yolk protein) genes [54] in L. cuprina and other flies to explore the evolutionary mechanisms leading to and maintaining these arrays. The solitary double-domain SCP/TAPS protein (LCSc0-DD1) characterised here is of particular interest; like solitary CAP domain proteins, double-domain proteins have also been associated with many crucial parasitic functions [15, 24, 55,56,57,58,59,60]. Given this association in other parasitic organisms, the presence of a double-domain protein in the L. cuprina genome warrants further investigation of the suspected role of LCSc0-DD1 in parasitism and whether LCSc0-DD1 could be a suitable vaccine target.

Analyses of the SCP/TAPS gene family in L. cuprina identified an extended tandem array of 16 SCP/TAPS genes located in a 53 kb domain of scaffold xfSc0000049. Although SCP/TAPS protein genes are represented in triplets in the genome of D. melanogaster [16] or aggregate as a series of four genes in the genome of the blood fluke Schistosoma mansoni [15], the tandem array identified here spans 16 genes and is significantly larger than aggregates seen in known SCP/TAPS tandem arrays. This appears to be unique to L. cuprina, as our analysis of 11 related arthropod genomes did not identify similar TAGs of the extent seen here (Additional file 7: Table S2). It is important to note that this analysis is preliminary, and future work should investigate a wider array of genomes in greater detail to establish whether or not this array is unique to L. cuprina. Nevertheless, our data suggest a large-scale gene duplication or replication event, resulting in a number of very similar and structurally conserved proteins forming this tandem array. An analysis of this region employing long-read sequencing data and RNA-seq data provided convincing evidence for the presence and transcription of 16 separate genes containing the SCP/TAPS domain. Tandem duplication events have been shown to drive the development and maintenance of gene superfamilies [61] and are a common mechanism used to adapt to harsh environments [62,63,64,65]. Alternatively, tandem array proteins may be a mechanism to promote production of large amounts of protein product, representing an alternative to other means of gene upregulation [66]. Given that L. cuprina TAGs are predominantly transcribed in the parasitic life stages (i.e., larval stages) and are transcribed at higher levels relative to other SCP/TAPS protein genes of L. cuprina, it is plausible that large-scale tandem duplications of these genes are a means by which L. cuprina has upregulated the production of these proteins. As SCP/TAPS proteins are commonly detected in excretory/secretory products [67] and proposed to modulate host immune responses [21], it is plausible that L. cuprina transcribes and expresses some that play a key role in host-larval interactions, reflected in their upregulation during parasitism. Investigating the cross-talk between fly larvae and the host animal could provide clues as to how to disrupt this relationship and might identify novel intervention targets for flystrike.

TAG function may also be linked to digestion as D. melanogaster proteins (Agr and Agr2) cluster phylogenetically with the 16 SCP/TAPS gene tandem array in L. cuprina and are known to be highly upregulated in digestive tissues during larval stages [48, 68, 69]. Synteny analysis of the associated genes revealed the presence of the same gene (rutabaga) upstream of both the L. cuprina tandem array and the Agr and Agr2 genes in D. melanogaster, suggesting a close evolutionary relationship between the encoded proteins. Agr is highly expressed in larval salivary glands, while Agr2 is predominately expressed in larval hindgut and midgut, signifying a role in larval digestion [48, 69]. The need for L. cuprina larvae to rapidly and effectively digest tissue could explain the necessity for a large-scale duplication/replication event resulting in a tandem array. Interestingly, transcription profiles for SCP/TAPS protein genes in L. cuprina were distinct concordant with those of D. melanogaster, suggesting that at least some of those transcribed and expressed in larval stages are required for a parasitic existence.

Despite advances in knowledge of the structural biology of CAP domains [20], their functions remain unclear, yet cysteine-rich secretory proteins have previously been the subject of successful vaccine candidate identification in several species of hookworm [70,71,72]. Hookworm L3 larvae are known to secrete two main CAP-like proteins, namely Ancylostoma secreted proteins 1 and 2 (i.e., ASP-1 and ASP-2) with two and one CAP domains, respectively [73]. Previous work showed that hamsters vaccinated with a recombinant form of Ac-ASP-2 exhibited a reduction in Ancylostoma ceylanicum burdens compared with control hamsters, while the full-length Ac-ASP-1 did not [71]. Intriguingly, immunisation with CAP domain-truncated form of Ac-ASP-1 was effective at reducing burdens compared with the full-length Ac-ASP-1, which did not [74,75,76], suggesting that the CAP domain and oligomerisation of the double-domain unit play a part in immunogenicity in the host animal. It is important to note that the inability to amplify a complete LCSc0-DD1 transcript using conventional and RACE PCR protocols could be in part due to the low transcription levels of this isoform or due to suboptimal primer design given the low GC content in the L. cuprina genome. Assessing the presence of other putative isoforms derived from LCSc0-DD1 should be the subject of further investigations as this might aid in elucidating the variability of double-domain structural isoforms and enhance our understanding of single- and double-CAP domain differences.

Conclusions

In conclusion, future investigations of the molecular processes underlying the biological functions of these intriguing proteins should be facilitated by advances in proteomic and transcriptomic data analyses as well as continuous improvement in genome sequencing and analyses. In this context, elucidating the biological function of SCP/TAPS through the use of functional genomics tools such as the CRISPR-cas9 genome-editing system [77] would be a logical next step. Initial success in full knockout of blowfly genes has been achieved [78], supporting the validity of such an approach. Establishing such a knockout system for L. cuprina could enable functional assessments of SCP/TAPS proteins as vaccine or insecticide targets providing prospects toward the development of new interventions against flystrike.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its Additional files. Any additional data are available from the corresponding author upon request. Accession codes: This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession JRES00000000.1 (BioProject: PRJNA419080, BioSample: SAMN08042461). The genome version described in this paper is version PPXR01000000. The RNAseq and Nanopore fastq files have been deposited in the NCBI short read archive (SRA) under the accession number PRJNA419080.

Abbreviations

DD:

Double-domain protein

LS:

Lab strain

SCP/TAPS:

Sperm-Coating Proteins/Tpx1/Ag5/PR-1/Sc7

SD:

Single-domain protein

TPM:

Transcripts per million

References

  1. Anstead CA, Batterham P, Korhonen PK, Young ND, Hall RS, Bowles VM, et al. A blow to the fly—Lucilia cuprina draft genome and transcriptome to support advances in biology and biotechnology. Biotechnol Adv. 2016;34:605–20.

    Article  CAS  PubMed  Google Scholar 

  2. Bisdorff B, Milnes A, Wall R. Prevalence and regional distribution of scab, lice and blowfly strike in Great Britain. Vet Rec. 2006;158:749–52.

    Article  CAS  PubMed  Google Scholar 

  3. Hall M, Wall R. Myiasis of humans and domestic animals. Adv Parasit. 1995;35:257–334.

    Article  CAS  Google Scholar 

  4. Azeredo-Espin AM, Lessinger AC. Genetic approaches for studying myiasis-causing flies: molecular markers and mitochondrial genomics. Genetica. 2006;126:111–31.

    Article  CAS  PubMed  Google Scholar 

  5. Williams KA, Richards CS, Villet MH. Predicting the geographic distribution of Lucilia sericata and Lucilia cuprina (Diptera: Calliphoridae) in South Africa. Afr Invertebr. 2014;55:157–70.

    Article  Google Scholar 

  6. Bisdorff B, Wall R. Sheep blowfly strike risk and management in Great Britain: a survey of current practice. Med Vet Entomol. 2008;22:303–8.

    Article  CAS  PubMed  Google Scholar 

  7. Anstead CA, Perry T, Richards S, Korhonen PK, Young ND, Bowles VM, et al. The battle against flystrike—past research and new prospects through genomics. Adv Parasit. 2017;98:227–81.

    Article  Google Scholar 

  8. Sandeman RM, Levot GW, Heath ACG, James PJ, Greeff JC, Scott MJ, et al. Control of the sheep blowfly in Australia and New Zealand – are we there yet? Int J Parasitol. 2014;44:879–91.

    Article  CAS  PubMed  Google Scholar 

  9. Levot GW. Cyromazine resistance detected in Australian sheep blowfly. Aust Vet J. 2012;90:433–7.

    Article  CAS  PubMed  Google Scholar 

  10. Levot GW. Resistance and the control of sheep ectoparasites. Int J Parasitol. 1995;25:1355–62.

    Article  CAS  PubMed  Google Scholar 

  11. Sneddon J, Rollin B. Mulesing and animal ethics. J Agric Environ Ethics. 2010;23:371–86.

    Article  Google Scholar 

  12. Anstead CA, Korhonen PK, Young ND, Hall RS, Jex AR, Murali SC, et al. Lucilia cuprina genome unlocks parasitic fly biology to underpin future interventions. Nat Commun. 2015;6:7344.

    Article  CAS  PubMed  Google Scholar 

  13. Cantacessi C, Campbell BE, Visser A, Geldhof P, Nolan MJ, Nisbet AJ, et al. A portrait of the “SCP/TAPS” proteins of eukaryotes—developing a framework for fundamental research and biotechnological outcomes. Biotechnol Adv. 2009;27:376–88.

    Article  CAS  PubMed  Google Scholar 

  14. Gibbs GM, Roelants K, O’Bryan MK. The CAP superfamily: cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins—roles in reproduction, cancer, and immune defense. Endocr Rev. 2008;29:865–97.

    Article  CAS  PubMed  Google Scholar 

  15. Chalmers IW, McArdle AJ, Coulson RMR, Wagner MA, Schmid R, Hirai H, et al. Developmentally regulated expression, alternative splicing and distinct sub-groupings in members of the Schistosoma mansoni venom allergen-like (SmVAL) gene family. BMC Genomics. 2008;9:89.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Kovalick GE, Griffin DL. Characterization of the SCP/TAPS gene family in Drosophila melanogaster. Insect Biochem Mol Biol. 2005;35:825–35.

    Article  CAS  PubMed  Google Scholar 

  17. Da Ros VG, Maldera JA, Willis WD, Cohen DJ, Goulding EH, Gelman DM, et al. Impaired sperm fertilizing ability in mice lacking Cysteine-RIch Secretory Protein 1 (CRISP1). Dev Biol. 2008;320:12–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Teixeira PJPL, Thomazella DPT, Vidal RO, Prado PFV, Reis O, Baroni RM, et al. The fungal pathogen Moniliophthora perniciosa has genes similar to plant PR-1 that are highly expressed during its interaction with cacao. PLoS ONE. 2012;7:e45929.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cantacessi C, Hofmann A, Young ND, Broder U, Hall RS, Loukas A, et al. Insights into SCP/TAPS proteins of liver flukes based on large-scale bioinformatic analyses of sequence datasets. PLoS ONE. 2012;7:e31164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Osman A, Wang CK, Winter A, Loukas A, Tribolet L, Gasser RB, et al. Hookworm SCP/TAPS protein structure—a key to understanding host-parasite interactions and developing new interventions. Biotechnol Adv. 2012;30:652–7.

    Article  CAS  PubMed  Google Scholar 

  21. Wilbers RHP, Schneiter R, Holterman MHM, Drurey C, Smant G, Asojo OA, et al. Secreted venom allergen-like proteins of helminths: Conserved modulators of host responses in animals and plants. PLoS Pathog. 2018;14:e1007300.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Moyle M, Foster DL, McGrath DE, Brown SM, Laroche Y, Meutter J, et al. A hookworm glycoprotein that inhibits neutrophil function is a ligand of the integrin CD11b/CD18. J Biol Chem. 1994;269:10008–15.

    Article  CAS  PubMed  Google Scholar 

  23. Winkler B, Bolwig C, Seppälä U, Spangfort MD, Ebner C, Wiedermann U. Allergen-specific immunosuppression by mucosal treatment with recombinant Ves v 5, a major allergen of Vespula vulgaris venom, in a murine model of wasp venom allergy. Immunology. 2003;110:376–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Bower MA, Constant SL, Mendez S. Necator americanus: The Na-ASP-2 protein secreted by the infective larvae induces neutrophil recruitment in vivo and in vitro. Exp Parasitol. 2008;118:569–75.

    Article  CAS  PubMed  Google Scholar 

  25. Asojo OA, Loukas A, Inan M, Barent R, Huang J, Plantz B, et al. Crystallization and preliminary X-ray analysis of Na-ASP-1, a multi-domain pathogenesis-related-1 protein from the human hookworm parasite Necator americanus. Acta Crystallogr F. 2005;61:391–4.

    Article  CAS  Google Scholar 

  26. Lu G, Villalba M, Coscia MR, Hoffman DR, King TP. Sequence analysis and antigenic cross-reactivity of a venom allergen, antigen 5, from hornets, wasps, and yellow jackets. J Immunol Res. 1993;150:2823.

    CAS  Google Scholar 

  27. King TP, Spangfort MD. Structure and biology of stinging insect venom allergens. Int Arch Allergy Immunol. 2000;123:99–106.

    Article  CAS  PubMed  Google Scholar 

  28. Charlab R, Valenzuela JG, Rowton ED, Ribeiro JMC. Toward an understanding of the biochemical and pharmacological complexity of the saliva of a hematophagous sand fly Lutzomyia longipalpis. Proc Natl Acad Sci. 1999;96:15155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Yamazaki Y, Morita T. Structure and function of snake venom cysteine-rich secretory proteins. Toxicon. 2004;44:227–31.

    Article  CAS  PubMed  Google Scholar 

  30. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tang S, Lomsadze A, Borodovsky M. Identification of protein coding regions in RNA transcripts. Nucleic Acids Res. 2015;43:e78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Kent WJ. BLAT - the BLAST-like alignment tool. Genome Res. 2002;12.

  35. Borodovsky M, Lomsadze A. Eukaryotic gene prediction using GeneMark.hmm-E and GeneMark-ES. Curr Protoc Bioinformatics. 2011;4:4.6.1-10.

    Google Scholar 

  36. Slater GS, Birney E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 2005;6:31.

    Article  CAS  Google Scholar 

  37. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat methods. 2012;9:357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Quinlan AR. BEDTools: the swiss-army tool for genome feature analysis. Curr Protoc Bioinform. 2014;47:11.2.1-34.

    Article  Google Scholar 

  41. Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:W187–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytolog Bull. 1987;19:11–5.

    Google Scholar 

  43. Sović I, Šikić M, Wilm A, Fenlon SN, Chen S, Nagarajan N. Fast and sensitive mapping of nanopore sequencing reads with GraphMap. Nat Commun. 2016;7:11307.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34:3094–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45:D200–3.

    Article  CAS  PubMed  Google Scholar 

  46. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Thurmond J, Goodman JL, Strelets VB, Attrill H, Gramates LS, Marygold SJ, et al. FlyBase 2.0: the next generation. Nucleic Acids Res. 2018;47:D759-65.

    Article  PubMed Central  CAS  Google Scholar 

  49. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12:323.

    Article  CAS  Google Scholar 

  50. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Stroehlein AJ, Young ND, Korhonen PK, Chang BCH, Sternberg PW, La Rosa G, et al. Analyses of compact Trichinella Kinomes reveal a MOS-like protein kinase with a unique N-terminal domain. G3 (Bethesda). 2016;6:2847–56.

    Article  CAS  Google Scholar 

  53. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40:D290-301.

    Article  CAS  PubMed  Google Scholar 

  54. Scott MJ, Atapattu A, Schiemann AH, Concha C, Henry R, Carey BL, et al. Organisation and expression of a cluster of yolk protein genes in the Australian sheep blowfly, Lucilia cuprina. Genetica. 2011;139:63–70.

    Article  CAS  PubMed  Google Scholar 

  55. Murray J, Gregory WF, Gomez-Escobar N, Atmadja AK, Maizels RM. Expression and immune recognition of Brugia malayi VAL-1, a homologue of vespid venom allergens and Ancylostoma secreted proteins. Mol Biochem Parasitol. 2001;118:89–96.

    Article  CAS  PubMed  Google Scholar 

  56. Hawdon JM, Narasimhan S, Hotez PJ. Ancylostoma secreted protein 2: cloning and characterization of a second member of a family of nematode secreted proteins from Ancylostoma caninum. Mol Biochem Parasitol. 1999;99:149–65.

    Article  CAS  PubMed  Google Scholar 

  57. Hawdon JM, Jones BF, Hoffman DR, Hotez PJ. Cloning and characterization of Ancylostoma-secreted protein. A novel protein associated with the transition to parasitism by infective hookworm larvae. J Biol Chem. 1996;271:6672–8.

    Article  CAS  PubMed  Google Scholar 

  58. Cantacessi C, Gasser RB. SCP/TAPS proteins in helminths - where to from now? Mol Cell Probes. 2012;26:54–9.

    Article  CAS  PubMed  Google Scholar 

  59. Ding X, Shields J, Allen R, Hussey RS. Molecular cloning and characterisation of a venom allergen AG5-like cDNA from Meloidogyne incognita. Int J Parasitol. 2000;30:77–81.

    Article  CAS  PubMed  Google Scholar 

  60. Stroehlein AJ, Young ND, Hall RS, Korhonen PK, Hofmann A, Sternberg PW, et al. CAP protein superfamily members in Toxocara canis. Parasit Vectors. 2016;9:360.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Reams AB, Neidle EL. Selection for gene clustering by tandem duplication. Annu Rev Microbiol. 2004;58:119–42.

    Article  CAS  PubMed  Google Scholar 

  62. Harms CT, Armour SL, DiMaio JJ, Middlesteadt LA, Murray D, Negrotto DV, et al. Herbicide resistance due to amplification of a mutant acetohydroxyacid synthase gene. Mol Gen Genet. 1992;233:427–35.

    Article  CAS  PubMed  Google Scholar 

  63. Lenormand T, Guillemaud T, Bourguet D, Raymond M. Evaluating gene flow using selected markers: a case study. Genetics. 1998;149:1383–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Anderson RP, Roth JR. Tandem genetic duplications in phage and bacteria. Annu Rev Microbiol. 1977;31:473–505.

    Article  CAS  PubMed  Google Scholar 

  65. Stark GR. Regulation and mechanisms of mammalian gene amplification. Adv Cancer Res. 1993;61:87–113.

    Article  CAS  PubMed  Google Scholar 

  66. Pan D, Zhang L. Tandemly arrayed genes in vertebrate genomes. Comp Funct Genom. 2008;2008:545269.

    Article  CAS  Google Scholar 

  67. Maizels RM, Gomez-Escobar N, Gregory WF, Murray J, Zang X. Immune evasion genes from filarial nematodes. Int J Parasitol. 2001;31:889–98.

    Article  CAS  PubMed  Google Scholar 

  68. Kovalick GE, Schreiber MC, Dickason AK, Cunningham RA. Structure and expression of the antigen 5-related gene of Drosophila melanogaster. Insect Biochem Mol Biol. 1998;28:491–500.

    Article  CAS  PubMed  Google Scholar 

  69. Megraw T, Kaufman TC, Kovalick GE. Sequence and expression of Drosophila antigen 5-related 2, a new member of the CAP gene family. Gene. 1998;222:297–304.

    Article  CAS  PubMed  Google Scholar 

  70. Goud GN, Bottazzi ME, Zhan B, Mendez S, Deumic V, Plieskatt J, et al. Expression of the Necator americanus hookworm larval antigen Na-ASP-2 in Pichia pastoris and purification of the recombinant protein for use in human clinical trials. Vaccine. 2005;23:4754–64.

    Article  CAS  PubMed  Google Scholar 

  71. Goud GN, Zhan B, Ghosh K, Loukas A, Hawdon J, Dobardzic A, et al. Cloning, yeast expression, isolation, and vaccine testing of recombinant Ancylostoma-secreted protein (ASP)-1 and ASP-2 from Ancylostoma ceylanicum. J Infect Dis. 2004;189:919–29.

    Article  CAS  PubMed  Google Scholar 

  72. Hotez PJ, Diemert D, Bacon KM, Beaumier C, Bethony JM, Bottazzi ME, et al. The human hookworm vaccine. Vaccine. 2013;31:B227–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Asojo O. Structure of a two-CAP-domain protein from the human hookworm parasite Necator americanus. Acta Crystallogr D. 2011;67:455–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Sen L, Ghosh K, Bin Z, Qiang S, Thompson MG, Hawdon JM, et al. Hookworm burden reductions in BALB/c mice vaccinated with recombinant Ancylostoma secreted proteins (ASPs) from Ancylostoma duodenale, Ancylostoma caninum and Necator americanus.. Vaccine. 2000;18:1096–102.

    Article  CAS  PubMed  Google Scholar 

  75. Ghosh K, Hawdon J, Hotez P. Vaccination with alum-precipitated recombinant Ancylostoma-secreted protein 1 protects mice against challenge infections with infective hookworm (Ancylostoma caninum) larvae. J Infect Dis. 1996;174:1380–3.

    Article  CAS  PubMed  Google Scholar 

  76. Ghosh K, Hotez PJ. Antibody-dependent reductions in mouse hookworm burden after vaccination with Ancylostoma caninum secreted protein 1. J Infect Dis. 1999;180:1674–81.

    Article  CAS  PubMed  Google Scholar 

  77. Port F, Chen H-M, Lee T, Bullock SL. Optimized CRISPR/Cas tools for efficient germline and somatic genome engineering in Drosophila. PNAS. 2014;111:2967–76.

    Article  CAS  Google Scholar 

  78. Paulo DF, Williamson ME, Arp AP, Li F, Sagel A, Skoda SR, et al. Specific gene disruption in the major livestock pests Cochliomyia hominivorax and Lucilia cuprina using CRISPR/Cas9. G3 (Bethesda). 2019;9:3045–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Funding from Australian Wool Innovation (AWI; to C.A.A. and T.P., and P.B. and R.B.G.) is gratefully acknowledged. AWI is grateful for its funding, which is primarily provided by Australian woolgrowers through a wool levy, and by the Australian Government which provides a matching contribution for eligible R&D activities. Funding also provided from the National Health and Medical Research Council (NHMRC) of Australia (to R.B.G. and N.D.Y.), the Australian Research Council (ARC) of Australia (to R.B.G.), and Melbourne Water Corporation (to R.B.G.). N.D.Y. holds an NHMRC Career Development Fellowship.

Author information

Authors and Affiliations

Authors

Contributions

YDJP, CAA, and TP planned the experimental design. YDJP performed the experiments, analysed the data and interpreted the results with major contributions in computational analysis from AJS, RSH, CAA, SK and NDY. Drafting of the manuscript was performed by YDJP with editing assistance from CAA, AJS, TP, NDY, SK, RBG, PB, and RSH, RG assisted with MinION sequencing. CAA, TP, PB, and RBG secured funding. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yair D. J. Prawer or Clare A. Anstead.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Table S1.

Primers used to verify the LCSc0-DD1 protein.

Additional file 2: Data S1.

Subset of all L. cuprina genomic scaffolds that contain mapped CAP transcripts.

Additional file 3: Data S2.

GFF file for all 27 SCP/TAPS protein genes.

Additional file 4: Data S3.

List of SCP/TAPS proteins predicted for L. cuprina with gene identifiers, scaffold location and predicted amino acid sequences. Gene identifiers are composed of the species of origin (LC: L. cuprina), scaffold number (Sc0: Scaffold Sc0000000) and whether the amino acid sequence produces a single or double-domain protein (SD or DD)

Additional file 5: Figure S1.

A protein model displaying the surface of LCSc49-SD12 and indicating the surface of the protein for the conserved CAP domain (red), for non-CAP domain residues (white), and cysteine residues that form disulphide bonds (yellow). Cysteine residue C3-C4 are not visible on the protein surface.

Additional file 6: Data S4. A:

L. cuprina SCP/TAPS best-predicted protein models in .pdb file format, produced with program I-TASSER and B: accompanying summary statists indicating accuracy of the predicted models. C-score is a confidence score ranging [-5,2]. TM-score and RMSD measure the similarity to the native structure. TM-score ranging [0,1] with a value >0.5 implying the model has correct topology.

Additional file 7: Table S2.

List of arthropod genomes with their respective number of predicted SCP/TAPS proteins as well as the largest number of SCP/TAPS genes clustered together in a tandem array.

Additional file 8: Data S5.

Output files loaded into IGV to manually determine SCP/TAPS gene models A: Mapped CAP transcripts using BLAT B: Gene prediction tool GeneMark C: Gene prediction tool Exonerate.

Additional file 9: Data S6.

LCSc0-DD1 amplicon sequencing data. A: Representation of the largest mapped amplicon to LCSc0-DD1, indicating two distinct CAP domains (red square brackets) and primer locations (red arrows). B: and C: LCSc0-DD1 amplicon sequence data in gb file format; sequenced in forward and reverse direction, respectively.

Additional file 10: Table S3.

Transcription levels (in transcripts per million, TPM) and read counts of SCP/TAPS genes for adult female, adult male and mixed larval libraries. Tandem array SCP/TAPS are indicated with an asterisk.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Prawer, Y.D.J., Stroehlein, A.J., Young, N.D. et al. Major SCP/TAPS protein expansion in Lucilia cuprina is associated with novel tandem array organisation and domain architecture. Parasites Vectors 13, 598 (2020). https://doi.org/10.1186/s13071-020-04476-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-020-04476-6

Keywords