Skip to main content

Identification, analysis, and linkage mapping of expressed sequence tags from the Australian sheep blowfly



The Australian sheep blowfly Lucilia cuprina (Wiedemann) (Diptera: Calliphoridae) is a destructive pest of the sheep, a model organism for insecticide resistance research, and a valuable tool for medical and forensic professionals. However, genomic information on L. cuprina is still sparse.


We report here the construction of an embryonic and 2 larval cDNA libraries for L. cuprina. A total of 29,816 expressed sequence tags (ESTs) were obtained and assembled into 7,464 unique clusters. The sequence collection captures a great diversity of genes, including those related to insecticide resistance (e.g., 12 cytochrome P450s, 2 glutathione S transferases, and 6 esterases). Compared to Drosophila melanogaster, codon preference is different in 13 of the 18 amino acids encoded by redundant codons, reflecting the lower overall GC content in L. cuprina. In addition, we demonstrated that the ESTs could be converted into informative gene markers by capitalizing on the known gene structures in the model organism D. melanogaster. We successfully assigned 41 genes to their respective chromosomes in L. cuprina. The relative locations of these loci revealed high but incomplete chromosomal synteny between L. cuprina and D. melanogaster.


Our results represent the first major transcriptomic undertaking in L. cuprina. These new genetic resources could be useful for the blowfly and insect research community.


The Australian sheep blowfly Lucilia cuprina (Wiedemann) (Diptera: Calliphoridae) is an important biological tool for medical treatment and forensic investigation. Disinfected blowfly larvae are routinely used in maggot debridement therapy to promote wound healing [13]. The necrophagous nature of L. cuprina also makes it invaluable for forensic analysis, particularly in estimating postmortem interval [4]. In contrast to these beneficial roles, L. cuprina is the primary cause of flystrike in Australia and New Zealand [5, 6]. The practice of surgical mulesing, as well as various chemical insecticides, has been used to control this formidable pest. However, fly populations often evolved resistance rapidly [710]. Research into the genetic and biochemical mechanisms of resistance has provided some of the best examples of genetic adaptation to selection [1113].

Despite its medical and agronomical importance and its historical status as one of the model organisms in insecticide resistance research, genomic information on L. cuprina is still relatively sparse. The haploid genome is approximately 810 mega bases [14], which is about 5 times the genome size of Drosophila melanogaster. Polytene in situ hybridization and genetic mapping studies have determined the basic organization of its 6 chromosomes [1517]. In 1993, Weller and Foster published a recombination map based on 72 morphological and enzyme markers, and this remains the most comprehensive linkage map of L. cuprina to date [18]. The chromosomal location of these markers indicates that the major linkage elements (i.e., Muller's Elements) remain relatively conserved in higher Diptera [18].

To combat this insect pest more intelligently, it is desirable to improve our knowledge of its genetic makeup. Molecular tools have become increasingly accessible to generate large amount of information in a cost-effective manner. As the per-base cost of DNA sequencing continues to fall, large-scale expressed sequence tag (EST) projects have been accomplished in many insect species, and thousands of ESTs have been deposited in public databases, including higher dipteran species such as the screwworm Cochliomyia hominivorax[19] and the tsetse fly Glossina morsitans[20].

This paper describes a similar gene discovery effort to identify transcripts expressed in preadult stages (embryonic and larval). An assembly of 7,464 unique gene clusters was produced from a total of 29,816 ESTs. The protein-coding contents of this non-redundant dataset were evaluated via a series of homology analyses. We short-listed a subset of these L. cuprina genes, which showed high sequence conservation, favorable gene structure (suitable exon/intron positions and sizes), and single correspondence in the D. melanogaster and Anopheles gambiae genomes. To demonstrate their usefulness in comparative mapping, we carried out chromosomal assignment of 41 genes to infer inter-chromosomal rearrangements. Comparison between L. cuprina and D. melanogaster revealed a high but incomplete chromosomal synteny. This newly generated EST dataset is a significant step in the systematic buildup of genomic resources for this important insect in agricultural and medical entomology.


Characteristics of EST assembly

We obtained 13,666 embryonic (JG422424-JG436089), 14,640 first-instar (JG407784-JG422423), and 880 third-instar (JG406904-JG407783) ESTs. The combined dataset (29,186 ESTs) was assembled into 7,464 unique clusters (Additional file 1), comprising 2,797 contigs and 4,667 singletons (Figure 1). The basic features of the sequence assembly are summarized in Figures 2 and 3. The majority (65.5%) of the 7,464 sequences were 500-900 bases in length (Figure 2), and the number of EST reads in a contig ranged between 2 and 3,700 (Figure 3), with an average of 8.8 reads (median = 3) per contig.

Figure 1

An overview of the acquisition, assembly, analysis, and application of L. cuprina -expressed sequence tags. A total of 29,816 ESTs from embryonic and larval libraries was assembled into 7,464 unique sequence clusters using the TGICL procedures. E-values from BLAST searches were arranged in ascending order from left to right, indicated by the darkness of shade.

Figure 2

Size distribution of the 7,464 non-redundant sequence clusters.

Figure 3

Frequency distribution of contig sizes.

Abundant transcripts in embryonic and larval stages

Since the cDNA libraries were not normalized, the number of reads in contigs can be used as an indication of levels of gene expression. Table 1 summarizes information about the most abundantly expressed genes, defined as contigs containing more than 100 ESTs. The mitochondria-derived transcripts dominate this list, accounting for a total of 5,275 ESTs, or 17% of the entire EST collection. In comparison to the published Lucilia sericata mitochondrial genome [21], our ESTs captured 12 protein coding genes and the 16S rRNA gene (Additional file 2). In addition to mitochondrial transcripts, various ribosomal protein genes (RpL6, RpL7A, RpL7, RpL4, and RpS3A), 18S rRNA, and elongation factors (ef1-α, ef2, and ef1-γ) were also amongst the most highly expressed genes (Table 1), reflecting the robust translation and protein synthesis processes in the embryonic and larval stages.

Table 1 The most highly expressed genes in the EST dataset, indicated by the number of ESTs in a contig

GC content and codon usage bias

Based on results from a set of 200 conserved genes (Additional file 3), the average GC content (mean ± standard deviation) per coding sequence (CDS) is 0.4344 ± 0.0433 in L. cuprina and 0.5654 ± 0.0418 in D. melanogaster. The effective number of codons (Nc) is 43.81 in L. cuprina and 40.89 in D. melanogaster. Compared to D. melanogaster, L. cuprina shows a different codon preference for 13 of the 18 amino acids encoded by redundant codons (Table 2). The most noticeable changes occur in the preferred codons for glutamine (Q), glutamic acid (E), and leucine (L).

Table 2 Codon usage comparison between L. cuprina and D. melanogaster based on 200 conserved genes

Protein-coding contents of EST assembly

To evaluate the protein coding contents of our ESTs, the 7,464 non-redundant sequences were subject to various homology searches against existing sequences (Figure 1). Our homology analyses showed that 5,257 (70%) of the non-redundant sequences had significant (E-value ≤ 1e-10) matches in the public domains (Additional file 4). The sequences that had recognizable homologs constituted 937 InterProScan and 494 Gene Ontology terms (Additional files 5 and 6), indicating that a great diversity of protein motifs and biological processes was represented in our dataset. We also estimated that ~78% (205 of 262) of the existing L. cuprina nucleotide sequences in Genbank were represented in our EST collection. Our Lucilia ESTs matched 3,280 unique D. melanogaster genes (from 3,409 unique polypeptides) at E-value ≤ 1e-10; this is equivalent to ~24.1% of the total gene count in D. melanogaster (assuming the total number of genes is ~13,600).

Identification of potential detoxification and insecticide target genes

We identified 12 cytochrome P450 (Cyp12a5, Cyp12d1, Cyp302a1 or disembodied, Cyp307a1 or spook, Cyp317a1, Cyp4d2, Cyp4g15, Cyp6a13, Cyp6d2, Cyp6d4, Cyp6d5, and Cyp9f2), 2 glutathione S-transferase (GstD1 and GstS1), and 6 esterase (Glt, Nrt, CG9289, alpha-Est5, CG9287, and alpha-Est7) homologs in the L. cuprina (Table 3). In addition to these detoxification gene families, we also identified ESTs encoding target proteins that have previously been implicated in insecticide resistance (see [2225]). These included homologs of the gamma-aminobutyric acid receptor-associated protein (GABA(A) receptor-associated protein; EST = GI: 333428695), glutamate receptor ionotropic kainate 2 (glutamate receptor 6; EST = GI: 333421827), and a probable sodium channel protein type 9 subunit alpha (EST = GI:333416352).

Table 3 Identification of potential detoxification genes in L. cuprina

Blast negatives in EST assembly

As of July 2010, 2,207 of the 7,464 unique gene clusters did not match any sequences in the public databases. Three hundred and sixty-five (or 16.5%) of these blast negatives had an ORF (minimum 20 codons), and the average length of their hypothetical polypeptide products was 126 amino acids (median = 126 amino acids; range = 20-584 amino acids) (details, see Additional file 7).

Anchor loci development and chromosomal synteny in higher Diptera

One aim of the current EST project was to identify single-copy genes that are highly conserved between species for synteny comparison. Reciprocal homology searches among L. cuprina, D. melanogaster, and A. gambiae yielded a set of reciprocal best-hit trios, and 298 of such trios had favorable intron position and size range in D. melanogaster (Additional file 8). This list of 298 orthologous groups facilitated our ongoing linkage map construction in L. cuprina. We successfully assigned 41 gene markers to 5 linkage groups using a male informative pedigree by scoring intron length polymorphisms (Figure 4). Comparison between L. cuprina and D. melanogaster based on 41 gene markers revealed a high level of synteny, although several deviations were also evident (Figure 5; Additional file 9). Deviations included inx3 RpL30, CG3564, RpS13, and RpL15., A potential translocation or fusion/dissociation event was identified between the smallest chromosome (Muller F in Drosophila) and an autosome (Muller D in Lucilia), as suggested by the location of RpS3A (Figure 5).

Figure 4

Mapping anchor loci by scoring intron length polymorphisms. Top right: Mapping pedigrees were generated using a backcross design initiated using the laboratory strain MI5 and the field strain Tara. Top left: Primers were first tested in the 4 backcross parents to identify intron length polymorphism. Bottom: Informative primer pairs were used to screen the backcross pedigree (TMM01).

Figure 5

Chromosomal synteny between L. cuprina and D. Melanogaster. Genomic locations of the anchor loci are shown on the 6 D. melanogaster chromosome arms (X, 2L, 2R, 3L, 3R, and 4). The corresponding chromosomal origins (II, III, IV, V, and VI) of the Lucilia homologs are indicated in brackets. Six instances of synteny violation are shown in red.


The main outcome of this project is the significant improvement of the gene inventory for the Australian sheep blowfly Lucilia cuprina. Amongst other applications, this new resource presents promising benefits to such areas as medical, forensic, pest control, and the understanding of genetic adaptation to insecticides.

Barring major gene expansion or contraction, and assuming that L. cuprina has the same number of genes as in D. melanogaster (~13,600) [26], the 7,464 unique gene clusters we found in our EST libraries would account for up to 55% of the genes present in the species. The actual percentage is much lower due to (but not limited to) the TGICL assembly parameters and the presence of 4.7% short (≤ 100 bases) sequences (Figure 2). An estimate of 24.1% gene coverage was obtained by limiting homology comparison to L. cuprina and D. melanogaster. However, fast-evolving genes and gene families that have been expanded in the blowfly lineage are under-represented in this analysis. Hence, the estimate of 24.1% could be considered the lower bound of total gene coverage. Nonetheless, this is a conservative yet reasonable estimation given that our cDNA libraries were not experimentally normalized and that only preadult developmental stages contributed to the transcript pool. The EST sequences contain a large number of recognizable protein motifs, as suggested by InterProScan results (Additional file 5), whose protein products are likely to participate in a myriad of biological and cellular processes, as also suggested by Gene Ontology analysis (Additional file 6).

Compared to D. melanogaster, L. cuprina appears to have low GC content and a different codon preference for many amino acids. Despite the fact that the comparison was based on 200 conserved gene homologs, the codon preferences for D. melanogaster are consistent with those reported by Vicario et al. [27]. The higher effective Nc in L. cuprina (43.81) than D. melanogaster (40.89) suggests a weaker selection constraint on codon usage in L. cuprina, at least for these highly conserved genes. It is noted that the 200 sequence pairs analyzed represent only a small fraction (1.5%) of the coding sequences in the 2 species; perhaps a different pattern might emerge when less-conserved gene homologs are included. Nevertheless, these results could be useful for training gene-finding algorithms and the analysis of the full genome sequence when it becomes available.

The acquisition of > 3,280 blowfly genes allows more sophisticated experimental systems to be developed in the future. Aside from the improvement in the knowledge about the genetic composition of the species, the dataset provides a foundation for designing gene-based microarrays for expression profiling. Furthermore, the plasmid collections can also serve as a permanent source of cDNA clones for protein expression, in situ hybridization, and even for transgenic manipulation such as those described in [2830]. The sequence knowledge of the housekeeping genes such as the ribosomal protein genes, tubulin, and actin could serve as internal controls for quantitative real-time PCR. In fact, the need for such reference genes was recently discussed in [31]. The availability of the L. cuprina cDNA sequences would also facilitate quantification of expression profiles of many genes of interest, bypassing the time-consuming gene discovery steps. It is expected that our EST collection will be invaluable for annotating the genic regions of the L. cuprina genome, when it is eventually sequenced. Conversely, the cDNA information could itself serve as a gene database, such that short peptides generated by the high-throughput proteome sequencing, similar to those reported in the brain tissues of another blowfly, Protophormia terraenovae[32], could be compared, forming a transcriptomic-proteomic feed-forward loop.

We identified genes that are related to insecticide resistance in L. cuprina (Table 3). Isolation of these homologs in L. cuprina would allow their expression patterns to be accurately measured (e.g., by real-time PCR), and their roles in insecticide resistance to be evaluated. PCR assays to screen for naturally occurring DNA polymorphisms (e.g., exon-primed intron-crossing (EPIC) markers) could also be developed to monitor the temporal and spatial distribution of different alleles. While many of their D. melanogaster homologs have been implicated in insecticide detoxification [3336], some of the genes identified are involved in other developmental processes such as ecdysone biosynthesis (disembodied and spook) [37, 38] and brain function/development (Cyp4g15) [39]. The proportions of the new L. cuprina homologs represent only a small fraction of these 3 detoxification gene families (see [4042]). With the advent of next-generation sequencing (NGS) technologies, large-scale genome or transcriptome sequencing has become increasingly popular. For example, transcriptomic analyses using NGS have now been reported in many non-model insect species [4348]. Similar approaches could be extended to L. cuprina and other related blowfly species, to enable a more comprehensive assessment of novel insecticide targets.

Another important application of our newly identified ESTs was to improve the genetic map of L. cuprina. ESTs can be converted to a set of anchor loci for linkage mapping, as has been repeatedly shown in other insects [49, 50]. We adopted a conservative "reciprocal best hit with strong homology" strategy in the selection of homologous markers, in which D. melanogaster served as the primary reference. A. gambiae, which diverged from the Lucilia and Drosophila lineages about 250 MYA, acted as an outgroup to improve the confidence in orthology calling, as sequence homology needed to reach the E-value of ≤ 1e-50 threshold to be considered further. In other words, we opted for confidence rather than sensitivity in our search for orthologous markers.

The EST-derived markers constituted a substantial proportion of anchor loci in the present study and were useful for inferring chromosomal synteny (Figure 5). Linkage assignment of 41 markers allowed us to conclude that chromosomal synteny is high between the Lucilia and Drosophila lineages. Our results are typical for higher dipteran species, as suggested by previous studies [17, 51, 52]. Several chromosomal fusion/dissociation events have nonetheless been documented within the Drosophila genus. For example, the fusion of Chr 4 (Muller's element F) to an autosome was found in Drosophila willistoni[53]. Moreover, comparison between mosquito and Drosophila reveals that patches of syntenic regions are scattered across many chromosome regions [54]. Our mapping results suggested that gene content on each of the Muller's elements in L. cuprina can, to a large extent, be predicted from the D. melanogaster map. However, the obvious cases of synteny violation (Figure 5; Additional file 9) would mean that direct extrapolation of linkage information from D. melanogaster would require extra caution. The interspecies differences should justify future de novo construction of linkage maps for L. cuprina, with denser markers.

The 298 putative orthologs effectively form a pipeline for future comparative mapping efforts (Additional file 8). Their chromosome addresses in D. melanogaster span virtually all regions of the genome, allowing flexible control over marker density for genomic regions of interest. Several chromosomal areas in Lucilia are of significant historical importance: the Scallop/Notch[55, 56] on chromosome II and the Rop-1[12, 57] regions on chromosome IV. In fact, several gene markers generated in the present study have already been utilized to understand the patterns of selective sweeps around the Rop-1 locus [58]. The marker pipeline also offers a starting point for fine scale mapping of the fitness modifier locus (M), which is believed to counter the fitness disadvantage of the diazinon-resistant flies in the absence of insecticide [5961]. Together with an appropriate genomic library, these newly acquired ESTs provide an ample supply of markers for positional cloning of the M locus.

The evolutionary origin and phylogeny relationship among blowfly species has been of great interest to many researchers, owing to its medical and forensic implications [62, 63]. With the much expanded gene repertoire, some of the L. cuprina genes identified here, especially those that show least similarity to other known sequences could be utilized to develop species diagnostic assays. The current EST sequences would greatly complement such an exploration.

While the assemblage of 29,816 ESTs into 7,464 was straight forward, the interpretation of the information contents requires regular re-adjustment, in light of the constantly expanding sequence databases in other species. In order to evaluate the coding components of the newly acquired sequences, they were sorted according to the level of homology to their counterparts in the Genbank reference protein database, producing a typical BLAST significance spectrum (Figure 1). It is anticipated that such a spectrum would change over time. As new sequences from other organisms become publicly accessible, it would simultaneously alter the structure of the existing sequence databases and hence the BLAST results. The recently released EST collections (116,737 reads) from 3 closely related taxa (Glossina, Cochliomyia and, Muscinae) clearly illustrate this notion (Figure 1).

Given that the number of non-redundant sequence clusters depends largely on the assembly settings, the "90% identity over 50 bases" requirement could be viewed as a balanced option, but might not be an optimized condition for all genes. One indication is the presence of residual sequence redundancy in the dataset, presumably due to the natural existence of splice variants, transcript isoforms, natural polymorphisms, or genuine gene families. Hence, it is worthwhile to disassemble relevant contigs that belong to the gene of interest and find the most appropriate parameters to reassemble these reads. Furthermore, we did not impose any restriction on the length of the sequences, i.e., removal of assembled contigs or reads less than a certain length (e.g., 200 bases), because such sequences could be part of the untranslated regions of many legitimate mRNA transcripts. As more similar EST sequences from closely related taxa become available, these short reads might ultimately be informative in the future. In summary, the TGICL assembly described in this paper only represents a generic, non-discriminatory clustering approach for the entire dataset, and re-assembling for the original ESTs might be necessary to produce the most accurate assembly for a given gene or a set of related genes.


We report the generation of 29,816 ESTs (7,464 unique clusters) from the Australian sheep blowfly Lucilia cuprina. Homology analyses revealed that the dataset captured a wide diversity of genes, including those related to insecticide resistance targets and detoxification gene families. Our data also indicate that L. cuprina coding sequences are AT rich and that codon usage patterns are distinct from that of D. melanogaster. In addition, a subset of putative orthologous genes was identified and mapped to the Lucilia linkage groups, which revealed a high but incomplete chromosomal synteny with D. melanogaster.


Construction and sequencing of cDNA libraries

Three cDNA libraries were constructed using RNA extracted from embryos, first-instar, and third-instar larvae. Construction of the embryonic cDNA library was previously described in Chen et al. (1998) [55]. Embryonic RNA was extracted from the non-modifier "seeking" strain using the Gibco-BRL mRNA Isolation System (Gaithersburg, MD). A unidirectional cDNA library was constructed in the Eco RI/Xho I sites of the λZAP II vector (Stratagene, La Jolla, CA), and sequenced from the 5' end using the T3 primer and the 3' end using the T7 primer. Two larval cDNA libraries (first and third instar) were made using total RNA from the inbred laboratory M15 strain. cDNA was prepared using the SMART® cDNA library construction kit (Clontech Laboratories, Inc.), directionally cloned into the αTripIEx2 vector via the Sfi-I A/Sfi-I B restriction sites, and transformed into BM25.8 competent cells. Plasmid cDNA clones were sequenced from the 5' end using the Sp6 primer. All sequencing was completed using the Sanger dideoxy sequencing method at the Australian Genome Research Facility (AGRF) in Brisbane, Australia. Original EST reads were subject to standard quality-trimming, vector-removal, and poly-A-clipping procedures. The output sequences were then assembled using the TGICL (TIGR Gene Indices clustering tools) algorithm [64] with the minimum threshold level set at 90% identity over a stretch of 50 bases.

Homology comparison

The GenBank non-redundant reference protein database (8,328,903 sequences) was downloaded from in July 2010. The complete L. sericata mitochondrial genome (GI:154623433) [21] was also retrieved from Genbank. The 18,648 EST sequences from C. hominivorax (primary screwworm), 18, 797 EST sequences (generated by 454 GS FLX; SRA: SRA012250) from Stomoxys calcitrans (the stable fly), and 79,292 EST sequences from G. morsitans (tsetse fly) were batch downloaded from GenBank via the species taxonomy page in the National Center for Biotechnology Information (NCBI) The D. melanogaster and Anopheles gambiae proteome sequences were retrieved from FlyBase and VectorBase, respectively. All sequences were converted into separate local databases using the NCBI standalone BLAST executables. Homology searches (BLASTX) were performed with E-value cut-off at 1e-10. To estimate gene coverage of our Lucilia dataset, we performed a BLASTX (E-value ≤ 1e-10) search against the D. melanogaster peptide database (r5.37). Due to the existence of isoforms in the BLASTX hit list, we extracted their corresponding gene identifiers (i.e., CG numbers and gene symbols) from FlyBase to estimate the number of unique genes. To identify putative orthologs among Lucilia, Drosophila and Anopheles, a more stringent BLAST E-value threshold (1e-50) was used to retain only the most conserved homolog pairs. EST sequences that failed to find a significant match in the reference protein database were then used to search (BLASTN) against the L. sericata mitochondrial genome, the C. hominivorax and the G. morsitans sequences with E-value cutoff at 1e-10. Finally, a TBLASTX search (at E-value < 1e-10) was performed for EST sequences that did not find sequence homology in both the BLASTX and BLASTN searches.

GC content and codon bias analyses

To compare the GC content and codon usage properties between L. cuprina and D. melanogaster, we confined our analyses to 200 homologous sequence pairs. These 200 homologous pairs are highly conserved (BLASTX E-value < 1e-50) at the amino acid level and have identifiable complete open reading frames (ORFs) (see Additional file 3). Putative ORFs in Lucilia were extracted using the GENSCAN program developed by Burge and Karlin [65]. The homologous ORFs (or CDS) in D. melanogaster were retrieved from FlyBase GC content and codon usage statistics were calculated using the GEECEE program The effective Nc was estimated using the CHIP program All 3 programs were available at BioManager maintained by Peter Reeves at the University of Sydney, Australia.

Pedigree construction

The M15 and the Tara strains were used to generate the male and female informative mapping families. The M15 strain is highly inbred and carries 1 visible phenotypic marker on each of its 5 autosomes. The Tara strain (provided by Garry Levot) is a more recent field-derived strain that originated from Tara, Queensland, Australia. It is morphologically wild type and displays high levels of resistance to diflubenzuron and tolerance to cyromazine. To generate the male informative family TMM1, an F1 male from a single pair mating between an M15 male and a Tara female was backcrossed to a virgin M15 female (Figure 4).

Choice of markers

cDNA sequences (previously characterized genes and ESTs) were converted into gene markers for linkage analysis and synteny comparison. Marker selection was based on (1) their physical locations in D. melanogaster, to ensure an even coverage of all Muller's elements; (2) that these genes contain intron(s) of suitable size (100-400 bases); and (3) that the intron positions are conserved in both D. melanogaster and A. gambiae. EPIC primers were designed using the Primer3 program (see Additional file 9 for primer information).

DNA isolation, polymerase chain reactions, and electrophoresis

DNA from all individuals in the pedigrees was extracted using DNAzol® reagent (Invitrogen; Cat. No. 10503027). In subsequent genotyping assays, 0.1% of the whole body DNA in 1 μL was used per PCR reaction. PCR was done in 25 μL reactions, which contained 1 μL of genomic DNA, 2.5 μL of 10 × reaction buffer, 3.0 μL of MgCl2 at 25 mM, 2.5 μL of dNTPs at 2 mM, 1 μL of each of the forward and reverse primer at 10 μL, 0.3 μL of Taq DNA polymerase (Fermentas; Cat. No. EP0402), and 13.7 μL of nuclease-free water. We used a touchdown thermo-cycling strategy for all PCR amplification, which involved an initial denaturation step at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 65°C for 30 s (reduce 0.5°C per cycle), and 72°C for 2 min, followed by another 30 cycles of 95°C for 30 s, 50°C for 30 s, and 72°C for 2 min. PCR amplicons were separated by electrophoresis.

The parents of the mapping crosses were first screened using EPIC primers for detectable size polymorphisms on a 1.2% agarose gel, which contained 1% agarose (Bioline; Cat. No. BIO-41025) and 0.2% Ultra-High Resolution Agarose (Scientifix; Cat. No. 9030A), at 250 V for 25 min. However, if the agarose electrophoresis did not reveal intron size polymorphism, the PCR amplicons were heat denatured (95°C for 3 min) and run on a 6% polyacrylamide gel (SequaGel®-6 system, National Diagnostics; Cat. No. EC-836 and EC-841) at 500 V for 1.25 h using the Gel-Scan 2000 system (Corbett Research). The polyacrylamide gels were stained with 1 × SYBRGold® I nucleic acid gel stain (Invitrogen; Cat. No. S-11494) to look for heteroduplex or single-strand conformation polymorphisms. If polymorphism was found in the parents of the mapping family, then identical procedures were applied to genotype the progeny.

Linkage analysis

Due to the general lack of meiotic crossing overs in dipteran males, markers on the same chromosome are transmitted together from the male parent to its progeny. Under our backcrossing schemes (Figure 4), all polymorphisms should have come from the Tara strain. Hence, the presence or absence of the Tara allele in the male informative family indicates the presence or absence of a specific Tara chromosome. Markers were assigned to the same linkage group if they shared identical segregation patterns in the male informative cross TMM1. Twenty-two backcross individuals from TMM1 were used in genotyping assays.


  1. 1.

    Paul AG, Ahmad NW, Lee HL, Ariff AM, Saranum M, Naicker AS, Osman Z: Maggot debridement therapy with Lucilia cuprina: a comparison with conventional debridement in diabetic foot ulcers. International wound journal. 2009, 6 (1): 39-46. 10.1111/j.1742-481X.2008.00564.x.

    PubMed  Article  Google Scholar 

  2. 2.

    Tantawi TI, Williams KA, Villet MH: An accidental but safe and effective use of Lucilia cuprina (Diptera: Calliphoridae) in maggot debridement therapy in Alexandria, Egypt. J Med Entomol. 2010, 47 (3): 491-494. 10.1603/ME09183.

    PubMed  Article  Google Scholar 

  3. 3.

    Sherman RA: Maggot therapy for foot and leg wounds. The international journal of lower extremity wounds. 2002, 1 (2): 135-142. 10.1177/1534734602001002009.

    PubMed  Article  Google Scholar 

  4. 4.

    Williams H: A model for the aging of fly larvae in forensic entomology. Forensic science international. 1984, 25 (3): 191-199. 10.1016/0379-0738(84)90193-2.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Heath AC, Bishop DM: Flystrike in New Zealand: An overview based on a 16-year study, following the introduction and dispersal of the Australian sheep blowfly, Lucilia cuprina Wiedemann (Diptera: Calliphoridae). Vet Parasitol. 2006, 137 (3-4): 333-344. 10.1016/j.vetpar.2006.01.006.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Watts JE, Murray MD, Graham NP: The blowfly strike problem of sheep in New South Wales. Aust Vet J. 1979, 55 (7): 325-334. 10.1111/j.1751-0813.1979.tb00419.x.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Blackman GG, Bakker JA: Resistance of the sheep blowfly Lucilia cuprina to insecticides in the Republic of South Africa. J S Afr Vet Assoc. 1975, 46 (4): 337-339.

    CAS  PubMed  Google Scholar 

  8. 8.

    Levot GW: Resistance and the control of sheep ectoparasites. International journal for parasitology. 1995, 25 (11): 1355-1362. 10.1016/0020-7519(95)00070-I.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    McKenzie JA: Dieldrin and diazinon resistance in populations of the Australian sheep blowfly, Lucilia cuprina, from sheep-grazing areas and rubbish tips. Aust J Biol Sci. 1984, 37 (5-6): 367-374.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Shanahan GJ: Resistance to dieldrin in Lucilia cuprina Wied., the Australian sheep blowfly. Nature. 1958, 181 (4612): 860-861. 10.1038/181860a0.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    McKenzie JA: Ecological and Evolutionary Aspects of Insecticide Resistance. 1996, Austin, Texas, U.S.A. R.G. Landes Company and Academic Press, Inc

    Google Scholar 

  12. 12.

    Newcomb RD, Campbell PM, Ollis DL, Cheah E, Russell RJ, Oakeshott JG: A single amino acid substitution converts a carboxylesterase to an organophosphorus hydrolase and confers insecticide resistance on a blowfly. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (14): 7464-7468. 10.1073/pnas.94.14.7464.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Newcomb RD, Gleeson DM, Yong CG, Russell RJ, Oakeshott JG: Multiple mutations and gene duplications conferring organophosphorus insecticide resistance have been selected at the Rop-1 locus of the sheep blowfly, Lucilia cuprina. J Mol Evol. 2005, 60 (2): 207-220. 10.1007/s00239-004-0104-x.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Ullerich FH, Schottke M: Karyotypes, constitutive heterochromatin, and genomic DNA values in the blowfly genera Chrysomya, Lucilia, and Protophormia (Diptera: Calliphoridae). Genome. 2006, 49 (6): 584-597. 10.1139/G06-013.

    PubMed  Article  Google Scholar 

  15. 15.

    Bedo DG: Polytene chromosomes of the Old World screwworm fly (Chrysomya bezziana) and its evolutionary relationships with Lucilia cuprina and Cochliomyia hominivorax (Diptera: Calliphoridae). Genome. 1992, 35 (2): 294-303. 10.1139/g92-045.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Perkins HD, Bedo DG, Howells AJ: Characterization and chromosomal distribution of a tandemly repeated DNA sequence from the Australian sheep blowfly, Lucilia cuprina. Chromosoma. 1992, 101 (5-4): 358-364.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Foster GG, Whitten MJ, Konovalov C, Arnold JTA, Maffi G: Autosomal genetic maps of the Australian Sheep Blowfly, Lucilia cuprina dorsalis R.-D. (Diptera: Calliphoridae), and possible correlations with the linkage maps of Musca domestica L. and Drosophila melanogaster (Mg.). Genetical Research. 1981, 37: 55-69. 10.1017/S0016672300020012.

    Article  Google Scholar 

  18. 18.

    Weller GL, Foster GG: Genetic maps of the sheep blowfly Lucilia cuprina: linkage-group correlations with other dipteran genera. Genome. 1993, 36 (3): 495-506. 10.1139/g93-068.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Guerrero FD, Dowd SE, Djikeng A, Wiley G, Macmil S, Saldivar L, Najar F, Roe BA: A database of expressed genes from Cochliomyia hominivorax (Diptera: Calliphoridae). Journal of medical entomology. 2009, 46 (5): 1109-1116. 10.1603/033.046.0518.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Attardo GM, Strickler-Dinglasan P, Perkin SA, Caler E, Bonaldo MF, Soares MB, El-Sayeed N, Aksoy S: Analysis of fat body transcriptome from the adult tsetse fly, Glossina morsitans morsitans. Insect molecular biology. 2006, 15 (4): 411-424. 10.1111/j.1365-2583.2006.00649.x.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Stevens JR, West H, Wall R: Mitochondrial genomes of the sheep blowfly, Lucilia sericata, and the secondary blowfly, Chrysomya megacephala. Med Vet Entomol. 2008, 22 (1): 89-91. 10.1111/j.1365-2915.2008.00710.x.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    ffrench-Constant RH, Steichen JC, Rocheleau TA, Aronstein K, Roush RT: A single-amino acid substitution in a gamma-aminobutyric acid subtype A receptor locus is associated with cyclodiene insecticide resistance in Drosophila populations. Proceedings of the National Academy of Sciences of the United States of America. 1993, 90 (5): 1957-1961. 10.1073/pnas.90.5.1957.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Hope M, Menzies M, Kemp D: Identification of a dieldrin resistance-associated mutation in Rhipicephalus (Boophilus) microplus (Acari: Ixodidae). J Econ Entomol. 2010, 103 (4): 1355-1359. 10.1603/EC09267.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Narahashi T: Neuronal Ion Channels as the Target Sites of Insecticides. Pharmacology & Toxicology. 1996, 79 (1): 1-14. 10.1111/j.1600-0773.1996.tb00234.x.

    CAS  Article  Google Scholar 

  25. 25.

    Raymond-Delpech V, Matsuda K, Sattelle BM, Rauh JJ, Sattelle DB: Ion channels: molecular targets of neuroactive insecticides. Invert Neurosci. 2005, 5 (3-4): 119-133. 10.1007/s10158-005-0004-9.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF, George RA, Lewis SE, Richards S, Ashburner M, Henderson SN, Sutton GG, Wortman JR, Yandell MD, Zhang Q, Chen LX, Brandon RC, Rogers YH, Blazej RG, Champe M, Pfeiffer BD, Wan KH, Doyle C, Baxter EG, Helt G, Nelson CR, et al: The genome sequence of Drosophila melanogaster. Science. 2000, 287 (5461): 2185-2195. 10.1126/science.287.5461.2185.

    PubMed  Article  Google Scholar 

  27. 27.

    Vicario S, Moriyama EN, Powell JR: Codon usage in twelve species of Drosophila. BMC evolutionary biology. 2007, 7: 226-10.1186/1471-2148-7-226.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Coates CJ, Howells AJ, O'Brochta DA, Atkinson PW: The 5' regulatory region from the Drosophila pseudoobscura hsp82 gene results in a high level of reporter gene expression in Lucilia cuprina embryos. Gene. 1996, 175 (1-2): 199-201. 10.1016/0378-1119(96)00149-7.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Heinrich JC, Li X, Henry RA, Haack N, Stringfellow L, Heath AC, Scott MJ: Germ-line transformation of the Australian sheep blowfly Lucilia cuprina. Insect molecular biology. 2002, 11 (1): 1-10. 10.1046/j.0962-1075.2001.00301.x.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Scott MJ, Heinrich JC, Li X: Progress towards the development of a transgenic strain of the Australian sheep blowfly (Lucilia cuprina) suitable for a male-only sterile release program. Insect biochemistry and molecular biology. 2004, 34 (2): 185-192. 10.1016/j.ibmb.2003.11.001.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Bagnall NH, Kotze AC: Evaluation of reference genes for real-time PCR quantification of gene expression in the Australian sheep blowfly, Lucilia cuprina. Med Vet Entomol. 2010, 24 (2): 176-181. 10.1111/j.1365-2915.2010.00866.x.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Inosaki A, Yasuda A, Shinada T, Ohfune Y, Numata H, Shiga S: Mass spectrometric analysis of peptides in brain neurosecretory cells and neurohemal organs in the adult blowfly, Protophormia terraenovae. Comparative biochemistry and physiology. 2010, 155 (2): 190-199. 10.1016/j.cbpa.2009.10.036.

    PubMed  Article  Google Scholar 

  33. 33.

    Le Goff G, Boundy S, Daborn PJ, Yen JL, Sofer L, Lind R, Sabourault C, Madi-Ravazzi L, ffrench-Constant RH: Microarray analysis of cytochrome P450 mediated insecticide resistance in Drosophila. Insect biochemistry and molecular biology. 2003, 33 (7): 701-708. 10.1016/S0965-1748(03)00064-X.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Low WY, Feil SC, Ng HL, Gorman MA, Morton CJ, Pyke J, McConville MJ, Bieri M, Mok YF, Robin C, Gooley PR, Parker MW, Batterham P: Recognition and detoxification of the insecticide DDT by Drosophila melanogaster glutathione S-transferase D1. Journal of molecular biology. 2010, 399 (3): 358-366. 10.1016/j.jmb.2010.04.020.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Willoughby L, Chung H, Lumb C, Robin C, Batterham P, Daborn PJ: A comparison of Drosophila melanogaster detoxification gene induction responses for six insecticides, caffeine and phenobarbital. Insect biochemistry and molecular biology. 2006, 36 (12): 934-942. 10.1016/j.ibmb.2006.09.004.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Chung H, Sztal T, Pasricha S, Sridhar M, Batterham P, Daborn PJ: Characterization of Drosophila melanogaster cytochrome P450 genes. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (14): 5731-5736. 10.1073/pnas.0812141106.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Chavez VM, Marques G, Delbecque JP, Kobayashi K, Hollingsworth M, Burr J, Natzle JE, O'Connor MB: The Drosophila disembodied gene controls late embryonic morphogenesis and codes for a cytochrome P450 enzyme that regulates embryonic ecdysone levels. Development. 2000, 127 (19): 4115-4126.

    CAS  PubMed  Google Scholar 

  38. 38.

    Warren JT, Petryk A, Marques G, Jarcho M, Parvy JP, Dauphin-Villemant C, O'Connor MB, Gilbert LI: Molecular and biochemical characterization of two P450 enzymes in the ecdysteroidogenic pathway of Drosophila melanogaster. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (17): 11043-11048. 10.1073/pnas.162375799.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Maibeche-Coisne M, Monti-Dedieu L, Aragon S, Dauphin-Villemant C: A new cytochrome P450 from Drosophila melanogaster, CYP4G15, expressed in the nervous system. Biochemical and biophysical research communications. 2000, 273 (3): 1132-1137. 10.1006/bbrc.2000.3058.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Low WY, Ng HL, Morton CJ, Parker MW, Batterham P, Robin C: Molecular evolution of glutathione S-transferases in the genus Drosophila. Genetics. 2007, 177 (3): 1363-1375. 10.1534/genetics.107.075838.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Ranson H, Claudianos C, Ortelli F, Abgrall C, Hemingway J, Sharakhova MV, Unger MF, Collins FH, Feyereisen R: Evolution of supergene families associated with insecticide resistance. Science. 2002, 298 (5591): 179-181. 10.1126/science.1076781.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Tijet N, Helvig C, Feyereisen R: The cytochrome P450 gene superfamily in Drosophila melanogaster: annotation, intron-exon organization and phylogeny. Gene. 2001, 262 (1-2): 189-198. 10.1016/S0378-1119(00)00533-3.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Vogel H, Heidel AJ, Heckel DG, Groot AT: Transcriptome analysis of the sex pheromone gland of the noctuid moth Heliothis virescens. BMC genomics. 2010, 11: 29-10.1186/1471-2164-11-29.

    PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Khajuria C, Zhu YC, Chen MS, Buschman LL, Higgins RA, Yao J, Crespo AL, Siegfried BD, Muthukrishnan S, Zhu KY: Expressed sequence tags from larval gut of the European corn borer (Ostrinia nubilalis): exploring candidate genes potentially involved in Bacillus thuringiensis toxicity and resistance. BMC genomics. 2009, 10: 286-10.1186/1471-2164-10-286.

    PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Walters JR, Harrison RG: EST analysis of male accessory glands from Heliconius butterflies with divergent mating systems. BMC genomics. 2008, 9: 592-10.1186/1471-2164-9-592.

    PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    O'Neil ST, Dzurisin JD, Carmichael RD, Lobo NF, Emrich SJ, Hellmann JJ: Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon. BMC genomics. 2010, 11: 310-10.1186/1471-2164-11-310.

    PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Molecular ecology. 2008, 17 (7): 1636-1647. 10.1111/j.1365-294X.2008.03666.x.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Ferguson L, Lee SF, Chamberlain N, Nadeau N, Joron M, Baxter S, Wilkinson P, Papanicolaou A, Kumar S, Kee TJ, Clark R, Davidson C, Glithero R, Beasley H, Vogel H, ffrench-Constant R, Jiggins C: Characterization of a hotspot for mimicry: assembly of a butterfly wing transcriptome to genomic sequence at the HmYbSb locus. Molecular ecology. 2010, 19 (Suppl 1): 240-254.

    PubMed  Article  Google Scholar 

  49. 49.

    Beldade P, Saenko SV, Pul N, Long AD: A gene-based linkage map for Bicyclus anynana butterflies allows for a comprehensive analysis of synteny with the lepidopteran reference genome. PLoS genetics. 2009, 5 (2): e1000366-10.1371/journal.pgen.1000366.

    PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Yasukochi Y, Ashakumary LA, Baba K, Yoshido A, Sahara K: A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects. Genetics. 2006, 173 (3): 1319-1328. 10.1534/genetics.106.055541.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Roethele JB, Feder JL, Berlocher SH, Kreitman ME, Lashkari DA: Toward a molecular genetic linkage map for the apple maggot fly (Diptera: Tephritidae): Comparison of alternative strategies. Ann Entomol Soc Am. 1997, 90 (4): 470-479.

    CAS  Article  Google Scholar 

  52. 52.

    Severson DW, Mori A, Kassner VA, Christensen BM: Comparative linkage maps for the mosquitoes, Aedes albopictus and Ae. aegypti, based on common RFLP loci. Insect molecular biology. 1995, 4 (1): 41-45. 10.1111/j.1365-2583.1995.tb00006.x.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Bhutkar A, Schaeffer SW, Russo SM, Xu M, Smith TF, Gelbart WM: Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes. Genetics. 2008, 179 (3): 1657-1680. 10.1534/genetics.107.086108.

    PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Zdobnov EM, von Mering C, Letunic I, Torrents D, Suyama M, Copley RR, Christophides GK, Thomasova D, Holt RA, Subramanian GM, Mueller HM, Dimopoulos G, Law JH, Wells MA, Birney E, Charlab R, Halpern AL, Kokoza E, Kraft CL, Lai Z, Lewis S, Louis C, Barillas-Mury C, Nusskern D, Rubin GM, Salzberg SL, Sutton GG, Topalis P, Wides R, Wincker P, et al: Comparative genome and proteome analysis of Anopheles gambiae and Drosophila melanogaster. Science. 2002, 298 (5591): 149-159. 10.1126/science.1077061.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Chen Z, Newsome T, McKenzie JA, Batterham P: Molecular characterization of the Notch homologue from the Australian sheep blowfly, Lucilia cuprina. Insect biochemistry and molecular biology. 1998, 28 (8): 601-612. 10.1016/S0965-1748(98)00032-0.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Davies AG, Game AY, Chen Z, Williams TJ, Goodall S, Yen JL, McKenzie JA, Batterham P: Scalloped wings is the Lucilia cuprina Notch homologue and a candidate for the modifier of fitness and asymmetry of diazinon resistance. Genetics. 1996, 143 (3): 1321-1337.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Russell RJ, Robin GC, Kostakos P, Newcomb RD, Boyce TM, Medveczky KM, Oakeshott JG: Molecular cloning of an alpha-esterase gene cluster on chromosome 3R of Drosophila melanogaster. Insect biochemistry and molecular biology. 1996, 26 (3): 235-247. 10.1016/0965-1748(95)00081-X.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Rose CJ, Chapman J, Marshall SD, Lee SF, Batterham P, Ross HA, Newcomb RD: Selective sweeps at the organophosphorus insecticide resistance locus, Rop-1, have affected variation across and beyond the {alpha}-esterase gene cluster in the Australian sheep blowfly, Lucilia cuprina. Molecular biology and evolution. 2011, 28 (6): 1835-1846. 10.1093/molbev/msr006.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    McKenzie JA, Game AY: Diazinon resistance in Lucilia cuprina; mapping of a fitness modifier. Heredity. 1987, 59: 371-381. 10.1038/hdy.1987.145.

    Article  Google Scholar 

  60. 60.

    McKenzie JA, Purvis A: Chromosomal localisation of fitness modifiers of diazinon resistance genotypes of Lucilia cuprina. Heredity. 1984, 53: 625-634. 10.1038/hdy.1984.120.

    Article  Google Scholar 

  61. 61.

    McKenzie JA, Whitten MJ, Adena MA: The effect of genetic background on the fitness of diazinon resistance genotypes of the Australian sheep blowfly, Lucilia cuprina. Heredity. 1982, 49: 1-9. 10.1038/hdy.1982.60.

    Article  Google Scholar 

  62. 62.

    Stevens JR: The evolution of myiasis in blowflies (Calliphoridae). International journal for parasitology. 2003, 33 (10): 1105-1113. 10.1016/S0020-7519(03)00136-X.

    PubMed  Article  Google Scholar 

  63. 63.

    Harvey ML, Gaudieri S, Villet MH, Dadour IR: A global study of forensically significant calliphorids: implications for identification. Forensic science international. 2008, 177 (1): 66-76. 10.1016/j.forsciint.2007.10.009.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19 (5): 651-652. 10.1093/bioinformatics/btg034.

    CAS  PubMed  Article  Google Scholar 

  65. 65.

    Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. Journal of molecular biology. 1997, 268 (1): 78-94. 10.1006/jmbi.1997.0951.

    CAS  PubMed  Article  Google Scholar 

Download references


We thank Aysha Hill-Williams and Alex Blasetti for fly maintenance, Yoganand Sundaravadanam for bioinformatics support, and Matthew Johnson and Peter Wilson for sequencing support. The project was funded by the Australian Wool Innovation grant to P. Batterham and M. Scott and the Australia Research Council grants to P. Batterham.

Author information



Corresponding author

Correspondence to Siu F Lee.

Additional information

Authors' contributions

SFL drafted the manuscript. ZC prepared the cDNA libraries and obtained the ESTs. SFL and ZC performed linkage mapping and quality control of the sequence assembly. AM performed the sequence assembly, gene ontology, and InterProScan. RTG performed homology searches. PB designed and supervised the project. All authors have read and approved the final manuscript.

Siu F Lee, Zhenzhong Chen contributed equally to this work.

Electronic supplementary material

non-redundant EST clusters

Additional file 1:L. cuprina non-redundant EST clusters. A sequence file containing 7,464 L. cuprina non-redundant EST clusters in FASTA format. (TXT 5 MB)

mitochondrial genes

Additional file 2:L. cuprina mitochondrial genes. A sequence file containing consensus sequences of L. cuprina mitochondrial genes in FASTA format. (TXT 13 KB)

Input file for GC and codon usage analyses

Additional file 3:. A sequence file containing 200 L. cuprina open reading frame sequences in FASTA format for GC and codon usage analyses. (TXT 149 KB)

Homology search results

Additional file 4:. A table containing BLAST hits of the non-redundant 7,464 L. cuprina EST sequences. (XLS 996 KB)

InterProScan results

Additional file 5:. Two spreadsheets containing InterProScan terms captured by the L. cuprina EST sequences. (XLS 350 KB)

Gene ontology results

Additional file 6:. Two spreadsheets containing Gene ontology terms captured by the L. cuprina EST sequences. (XLS 664 KB)

Additional file 7:BLAST-negative EST clusters with protein coding potential. A table containing a list of BLAST-negative EST clusters that have a hypothetical ORF of minimum 20 amino acids. (XLS 41 KB)

Potential orthologous genes among

Additional file 8: L. cuprina, D. melanogaster and A. gambiae. A table containing accession numbers corresponding to orthologous genes among L. cuprina, D. melanogaster, and A. gambiae. (XLS 68 KB)

Synteny between

Additional file 9: L. cuprina and D. melanogaster and primer information. A table containing chromosomal locations of L. cuprina genes and their corresponding primer sequences. (XLS 30 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Lee, S.F., Chen, Z., McGrath, A. et al. Identification, analysis, and linkage mapping of expressed sequence tags from the Australian sheep blowfly. BMC Genomics 12, 406 (2011).

Download citation


  • Codon Usage
  • Insecticide Resistance
  • Codon Preference
  • Chromosomal Synteny
  • Australian Sheep