Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury
- Rafaela M Maia†1,
- Valeria Valente†1,
- Marco AV Cunha1,
- Josane F Sousa1,
- Daniela D Araujo1, 5,
- Wilson A SilvaJr2,
- Marco A Zago3,
- Emmanuel Dias-Neto4, 6, 9,
- Sandro J Souza4,
- Andrew JG Simpson4, 7,
- Nadia Monesi1, 8,
- Ricardo GP Ramos1,
- Enilza M Espreafico1 and
- Maria L Paçó-Larson1Email author
© Maia et al; licensee BioMed Central Ltd. 2007
Received: 13 November 2006
Accepted: 24 July 2007
Published: 24 July 2007
The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome.
Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury.
Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.
Genome sequence determination of the model organism Drosophila melanogaster was a landmark that launched a new era for functional genomic studies in complex organisms. The almost complete version of the euchromatic DNA sequence was first released in March 2000 due to a collaborative effort of the Drosophila Genome Projects and Celera Genomics . Using gene prediction softwares in combination with searches of protein and EST databases, initial in silico analyses indicated the existence of 13,601 protein-coding genes (PCG), an extraordinarily small number of genes when compared to the approximately 19.000 PCG encoded in the C.elegans genome .
After the release 1, an intensive collective work took place in order to improve sequence quality and annotation, fill in the gaps, and correct the assembly. With the aim of generating the information necessary to define the transcripts encoded in the genome, the Berkeley Drosophila Genome Project (BDGP) initiated a high throughput production of both EST and full length cDNA sequences based on conventional and normalized cDNA libraries from different tissues and developmental stages . This effort was followed by non-BDGP projects with a major contribution from the Exelixis Drosophila melanogaster EST project, which has adopted sequencing of random primed libraries of mixed stage embryos, imaginal disks, and adult heads to increase the transcription units coverage . Currently, there are about 39,346 full length mRNA and 532,557 EST sequences available in the NCBI database, totalizing approximately 16,681 clusters according to UniGene . Since the year 2000, several subsequent genome versions have been released, each one improved by BDGP and annotated by FlyBase . Release 3.2, considered the first finished version, was published in March 2004 and provided a complete revision of all gene models and other genome features , estimating a total number of 13,792 PCGs plus 527 non-protein coding genes (tRNAs, rRNAs, microRNAs, sn/snoRNAs). Release 4.3, the last annotated genome version published in March 2006, includes a total of 14,816 genes and is available for searches by gene annotation, BLAST or sequence ID at the FlyBase website .
During the last few years, an enthusiastic debate about the number of PCGs in the organisms with sequenced genomes has arisen. For D. melanogaster, estimates varied from the initial ~ 13,600 coding gene predictions  to about 16,000 gene predictions, based on microarray expression data . A careful computational and experimental analysis carried to validate the Drosophila genome annotation has recently concluded that the D.melanogaster genome in fact contains approximately 14,000 protein-coding genes, although some genes presenting unusual features that make them refractory to prediction methods may remain to be discovered . However, the truthful notion about the complexity of the D. melanogaster transcriptome is still under construction. In this respect, it has been inferred from DNA oligonucleotide microarrays, with unique sequences tiled throughout the genome and across predicted splice junctions, that over 40% of the Drosophila genes contain one or more alternative exons . Additionally, a transcription map with a 35 bp resolution of the initial 24 hours of development indicates that 30% of the transcribed regions are still unannotated. Approximately 23% of these are intronic and 7% correspond to intergenic regions. Based on manual and computational surveys designed to identify coregulated expression patterns between unannotated and annotated genome regions, it was estimated that 29% of the unannotated regions are part of transcripts incompletely annotated or potential alternative exons of known genes . Therefore, correcting and refining the genome annotation is a reiterative task, which is continuously being done and depends on experimental data for final validation, especially for the identification of rare transcripts and alternative splice variants. With the aim of covering the diversity of transcripts expressed in Drosophila the generation of EST information from different sources is currently under way [2, 3].
Here we use the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is based on low stringency RT-PCR, to generate D. melanogaster expressed sequence information. ORESTES are preferentially derived from the central coding portions of the transcript and frequently identify less abundant messages [12, 13]. Such approach was previously applied for human transcriptome characterization, validating a large percentage of genes and identifying 219 unannotated transcribed sequences on chromosome 22 . More recently, a large-scale analysis of ORESTES derived from head, neck and thyroid tumors pointed to 788 putative new alternative splicing isoforms. A subset of 34 was submitted to experimental validation resulting in the confirmation of 23 (68%) new alternate exons .
Analysis of 1,303 Drosophila ORESTES clusters revealed 68 potential transcribed regions unannotated in the current version of the genome (release 4.3). Experimental validation of 38 (~ 50%) of this unannotated ORESTES revealed 17 new exons that most likely belong to low abundance transcripts. Using the ORESTES information together with a PCR based approach we obtained the complete coding sequence of a new serine protease which mRNA expression is induced upon infection. Our data reinforce the importance of PCR based methodologies for refining the Drosophila transcriptome, particularly for the identification of previously unannotated low copy transcripts.
Results and discussion
ORESTES in silico analysis
50% of the unannotated exons detected by ORESTES belong to low abundance transcripts
The majority of the validated ORESTES (14) detected a single band in northen blots, which could indicate the presence of unique transcripts. Three ORESTES detected more than one mRNA species (Figure 2, ORE-3, ORE-15, ORE-16), which could either constitute isoforms of the same gene or mRNAs encoded by different genes sharing common exons. Ten of these 17 validated ORESTES detected transcripts in all analysed developmental stages, namely: embryo, larvae and adult (Figure 2; ORE-1, -2, -3, -5, -6, -9, -11, -12, -14, -15). The other seven (Figure 2; ORE -4, -7, -8, -10, -13, -16, -17) detected stage specific transcripts. ORE-4 and ORE-10 detect mRNAs of 4.5 and 11.5 kb, respectively, which are mainly expressed in embryos. The 8.5 kb mRNA detected by ORE-17 is abundant in embryos and present in lower amounts in adults. ORE-12 detects an 11.0 kb transcript present mainly in embryo and adults. The transcripts of about 0.9 kb and 3.5 kb detected by ORE-8 and ORE-13, respectively, were only detected in larvae and adults. ORE-7 detects a 0.6 kb mRNA present at higher levels in larvae that is also expressed in embryos. ORE-3 detected one mRNA of about 4.5 kb exclusively expressed in embryos, and a 2.1 kb mRNA present at all stages. ORE-15 detected three mRNAs of 9.8, 5.0 and 1.8 kb. The 9.8 kb mRNA is mainly detected in larvae; the 5.0 kb RNA is much more abundant in embryos, but also detected in adults; and the 1.8 kb RNA is exclusively expressed in adults. ORE-16 detected 3 different transcripts of about 4.0, 3.1 and 2.2 kb, all of them observed in adult flies.
Genomic mapping of the validated ORESTES and sizes of the respective transcripts
Closest annotated gene
Predicted transcript sizesβ of the closest gene
1.8; 5.0; 9.8
2.2; 3.1; 4.0
Cloning a new D. melanogaster serine-protease – SP212
SP proteins are proteolytic enzymes that require serine for their catalytic activity. They are ubiquitous peptidases, which perform a wide array of important physiological functions, including digestion, blood coagulation, fibrinolysis, cellular and humoral immunity, fertilization and embryonic development . In a previous work that intended to map all SPs and SP related proteins in the genome of D. melanogaster, Ross and colleagues  performed a series of similarity searches (PSI-BLAST) and found a total of 211 GenBank entries encoding SPs and SPHs (SP homologs – proteins in which one or more residues of the catalytic triad are missing). Following the nomenclature criteria used by these authors, we named the additional SP gene characterized here as SP212. PROSITE analysis  revealed that D. melanogaster SP212 presents the conserved catalytic triad ordered His, Asp, Ser (HDS), a characteristic of SP chymotrypsin family. These residues form two diads, Ser-His and His-Asp, that operate in concert for the acyl mechanism of catalysis . These catalytic residues are, as in most SPs, embedded into highly conserved motifs: SAAHC, DIAL and GDSGGG. SP212 also presents the signal peptide sequence in its amino terminal end with the cleavage between residues 17 and 18, an indication that this peptidase might be secreted (figure 3B).
SP212 gene is upregulated in response to infection
The analysis of a relatively small set of ORESTES allowed the validation of 17 novel exons present in low abundance transcripts, and led to the cloning of a new serine peptidase that is induced during the defense response. These results illustrate the importance of PCR-based approaches as complementary tools for the identification of transcribed regions in sequenced genomes. The final determination of any transcriptome is not a trivial task and one might envisage the occurrence of rare transcripts that will be missed by conventional cloning.
Biological samples, ORESTES preparation and sequencing
Dechorionated embryos, larvae plus prepupae and pupae as well as adult flies were collected from an isogenic y, w1118 stock of Drosophila melanogaster, immediately frozen in liquid nitrogen and stored at -80 C until use. Total RNA was extracted with TRIZOL and poly(A)+ RNA was isolated (MiniMacs; Miltenyi Biotec). The RNA quality was assessed by northern blot hybridization using a Drosophila α1-tubulin probe . High quality total RNA preparations were further treated with DNase I (Promega). The absence of DNA contaminants was assessed by Southern blot hybridization of PCR amplification products using Drosophila mitochondrial DNA primers. Template preparations were performed as described by Dias-Neto and colleagues  with some minor modifications as follows: 15 ng of purified mRNAs were used for cDNA synthesis and amplification, using RT-PCR beads (Amershan-Pharmacia Biotech, USA) and a set of randomly selected oligonucleotide primers (15 to 20 mers). ORESTES profiles were generated after a cDNA synthesis step at 42°C for 60 mins immediately followed by cDNA denaturation at 75°C and amplification by PCR using a multiple annealing step. The annealing was performed for 10 secs at each temperature and the temperatures varied from 66°C to 44°C (with progressive reductions of 2°C), within each cycle. Primer extension was performed at 72°C for 1 min and denaturing at 95°C for 45 secs in 40 cycles. A final extension step at 72°C for 7 mins was undertaken. Profiles composed of a DNA smear were size selected in order to separate amplification products of distinct size ranges, varying from 0.3 to 1.5 kb. The fragments were ligated into pUC18 using Sureclone (Amersham-Pharmacia) and the recombinant plasmids used for bacterial (DH5α E.coli) transformation. The resulting colonies were grown overnight in liquid media and used as templates for PCR using vector primers. One microliter of the resulting PCR product was used for DNA sequencing using standard protocols of the ThermoSequenase II dye terminator cycle sequencing kit (Amershan-Pharmacia Biotech) and the reactions run on a MegaBACE 1000 automated sequencer.
A set of 10,092 ORESTES from different developmental stages of Drosophila melanogaster were generated. These sequences were submitted to an automated protocol for checking sequence quality, trimming to exclude vector and primer sequences, removing mtDNA, rRNA, bacterial and yeast and masking repetitive elements resulting in a total of 9,081 ORESTES (GeneBank_Accn EG974084 to EG974151 and ES688489 to ES697501): 360 from embryos at various stages (DE), 2,207 from larvae plus prepupae and pupae (DL), 4,490 from adults (DA) and 2,024 derived from a RNA mixture from embryos, larvae, pre-pupae, pupa and adults (DP). To assess the quality values of the sequences we used phred [36, 37]. "N" nucleotides were trimmed by using a PERL script (cleanN). Sequences corresponding to pUC18 and primers were identified by crossmatch (cleanup_vector). rRNA and mtDNA sequences were identified using a program based on FASTA3  and bacterial and yeast DNA sequences were identified by BLAST . The databases for rRNA, mtDNA, bacterial DNA (E.coli) and yeast were compiled from GenBank. These ORESTES were processed with the assembly tool Cap3 . Clustering of the 9,081 ORESTES resulted in 1,303 non-redundant clusters: 575 contigs plus 728 singletons. The 1,303 obtained clusters were submitted to automated BLASTn search against the release 1 annotated genes  resulting in the identification of 176 unannotated ORESTES clusters. Clustering of the ORESTES generated from each library resulted in 113 non-redundant clusters from DE, 430 from DL, 887 from DA and 729 from DP libraries.
RNA extraction and northern blot analysis
Drosophila total RNA was extracted from embryo, larvae plus prepupae and pupae or adults. After homogenization in lysis buffer (10 mM Tris-HCl, 2% SDS, 50 mM EDTA, 5% ethanol, pH 9.0), total RNA was extracted by adding 10 vol of Trizol, following the manufacter's instructions (Invitrogen). RNA PolyA+ was obtained using the Oligotex kit (Qiagen). The RNA was fractionated in 1% agarose formaldehyde-denaturing agarose gels and blotted to nylon membranes (Hybond N, Amershan, UK). RT-PCR, cloning, northern blotting, probe labelling, hybridization and post hybridization washes were performed essentially as described in Sambrook et al. . The final washes were performed at 65°C in the presence of 0.1× SSC and 0.2% SDS. Primers used for SP212 cloning: Foward-exon1-1: 5' TCA GTC TTA TTT GCC CAC CG 3'; Foward-exon1-2: 5' GTG TCG CTA ATC GCC TTG G 3'; Foward-exon2-1: 5' GCC TCC GCC GTG GGT TCC 3'; Reverse-exon2-1: 5' CCG ATC CTG TAA GCT GTC G 3'; Reverse-exon3-ORESTES: 5' TGC CAC GGG ATG AGG TAG G 3'; Reverse-exon4-1: 5' ACA GGG TCC AGA TCC ATC G 3'; Reverse-exon4-2: 5' CTG ATT AAG CTG GCA GGT GC 3'.
The insects were challenged by pricking with sharpened needles, which had been previously dipped into concentrated cultures of microorganisms. In the asceptic pricking the needles were first disinfected by ethanol. Microorganisms. E.coli (Gram-) was cultured in LB medium and Staphylococcus aureus ( Gram+) was grown in BHI medium (a gift from P.S.R. Coelho). Aspergillus fumigatus was grown on Sabouraud-agar medium. Spores and hyphae were harvested in saline (a gift from Dr. C. Maffei).
We thank Fernanda S. Zanola, Adriana A. Marques, Cristiane Ayres Ferreira, Alexandre C. de Oliveira, Benedita O. de Souza and Cirlei A.V. Saraiva for their dedicated technical assistance, Valdir Mazzucatto for the maintenance of the Drosophila room and Juçara Parra for administrative assistance. The FAPESP/Ludwig Institute for Cancer Research Consortium, as well as a FAPESP grant (MLPL) and stipends from FAEPA-FMRP supported the work. Sequencing and bio-informatics analysis were carried out at the Centro de Terapia Celular supoported by FAPESP. RMM received a fellowship from CAPES. VV and JFS were supported by FAPESP fellowships.
- Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF, et al: The genome sequence of Drosophila melanogaster. Science. 2000, 287: 2185-2195. 10.1126/science.287.5461.2185.PubMedView Article
- BDGP EST Submitted Collections. [http://www.fruitfly.org/EST/faq.html#cdna-1]
- Non-BDGP Fly EST Projects. [http://www.fruitfly.org/EST/otherEST.html]
- UniGene: An Organized View of the Transcriptome. [http://www.ncbi.nlm.nih.gov/UniGene/]
- Flybase: A database of the Drosophila Genome. [http://flybase.bio.indiana.edu/]
- Celniker SE, Wheeler DA, Kronmiller B, Carlson JW, Halpern A, Patel S, Adams M, Champe M, Dugan SP, Frise E, Hodgson A, et al: Finishing a whole-genome shotgun: release 3 of the Drosophila melanogaster euchromatic genome sequence. Genome Biology. 2002, 3: 1-14. 10.1186/gb-2002-3-12-research0079.View Article
- Misra S, Crosby MA, Mungall CJ, Matthews BB, Campbell KS, Hradecky P, Huang Y, Kaminker JS, Millburn GH, Prochnik SE, et al: Annotation of the Drosophila melanogaster euchromatic genome: a systematic review. Genome Biol. 2002, 3: research0083.1-0083.22. 10.1186/gb-2002-3-12-research0083.View Article
- Hild M, Beckmann B, Haas SA, Koch B, Solovyev V, Busold C, Fellenberg K, Boutros M, Vingron M, Sauer F, et al: An integrated gene annotation and transcriptional profiling approach towards the full gene content of the Drosophila genome. Genome Biology. 2003, 5: R3: 1-R3.
- Yandell M, Bailey AM, Misra S, Shu SQ, Wiel C, Evans-Holm M, Celniker SE, Rubin GM: A computational and experimental approach to validating annotations and gene predictions in the Drosophila melanogaster genome. Proc Natl Acad Sci USA. 2005, 102: 1566-1571. 10.1073/pnas.0409421102.PubMed CentralPubMedView Article
- Stolc V, Gauhar Z, Mason C, Halasz G, van Batenburg MF, Rifkin SA, Hua S, Herreman T, Tongprasit W, Barbano PE, et al: A gene expression map for the euchromatic genome of Drosophila melanogaster. Science. 2004, 306: 655-660. 10.1126/science.1101312.PubMedView Article
- Manak JR, Dike S, Sementchenko V, Kapranov P, Biemar F, Long J, Cheng J, Bell I, Ghosh S, Piccolboni A, Gingeras TR: Biological function of unannotated transcription during the early development of Drosophila melanogaster. Nature Genetics. 2004, 10: 1151-1158. 10.1038/nm1104-1151.
- Dias Neto E, Correa RG, Verjovski-Almeida S, Briones MR, Nagai MA, da Silva W, Zago MA, Bordin S, Costa FF, Goldman GH, Carvalho AF, et al: Shotgun sequencing of the human transcriptome with OFR expressed sequence tags. Proc Natl Acad Sci USA. 2000, 97: 3491-3496. 10.1073/pnas.97.7.3491.PubMed CentralPubMedView Article
- Camargo AA, Samaia HP, Dias-Neto E, Simao DF, Migotto IA, Briones MR, Costa FF, Nagai MA, Verjovski-Almeida S, Zago MA, et al: The contribution of 700,000 ORF sequence tags to the definition of the human. Proc Natl Acad Sci USA. 2001, 98: 12103-12108. 10.1073/pnas.201182798.PubMed CentralPubMedView Article
- de Souza SJ, Camargo AA, Briones MR, Costa FF, Nagai MA, Verjovski-Almeida S, Zago MA, Andrade LE, Carrer H, El-Dorry HF: Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags. Proc Natl Acad Sci USA. 2000, 97: 12690-3. 10.1073/pnas.97.23.12690.PubMed CentralPubMedView Article
- Reis EM, Ojopi EPB, Alberto FL, Rahal P, Tsukumo F, Mancini UM, Guimarães GS, Thompson GMA, Camacho C, Miracca E, et al: Large-scale Transcriptome Analyses Reveal New Genetic Marker Candidates of Head, Neck, and Thyroid Cancer. Cancer Res. 2005, 65: 1693-9. 10.1158/0008-5472.CAN-04-3506.PubMedView Article
- D. melanogaster BLAT Search. [http://genome.ucsc.edu/cgi-bin/hgBlat]
- Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T: The Small RNA Profile during Drosophila melanogaster Development. Developmental Cell. 2003, 5: 337-350. 10.1016/S1534-5807(03)00228-4.PubMedView Article
- Stapleton M, Carlson J, Brokstein P, Yu C, Champe M, George R, Guarin H, Kronmiller B, Pacleb J, Park S, Wan K, et al: Drosophila full-length cDNA resource. Drosophila. 2002, 3 (12): research0080.1-0080.8.
- Stapleton M, Liao G, Brokstein P, Hong L, Carninci P, Shiraki T, Hayashizaki Y, Champe M, Pacleb J, Wan K, et al: The Drosophila Gene Collection: Identification of Putative Full-Length cDNAs for 70% of D. melanogaster Genes. Genome Research. 2002, 12: 1294-1300. 10.1101/gr.269102.PubMed CentralPubMedView Article
- The Ludwig-FAPESP Transcript Finishing Initiative, Sogayar MC, Camargo AA: Transcript Finishing Initiative for Closing Gaps in the Human Transcriptome. Genome Research. 2004, 14: 1413-1423. 10.1101/gr.2111304.PubMed CentralView Article
- Rawlings ND, Barrett AJ: Evolutionary families of peptidases. Biochem J. 1993, 290: 205-218.PubMed CentralPubMedView Article
- Ross J, Jiang H, Kanost MR, Wang Y: Serine proteases and their homologs in the Drosophila melanogaster genome: an initial analysis of sequence conservation and phylogenetic relationships. Gene. 2003, 304: 117-131. 10.1016/S0378-1119(02)01187-3.PubMedView Article
- PROSITE: Database of protein families and domains. [http://www.expasy.org/prosite]
- Perona JJ, Craik CS: Structural basis of substrate specificity in the serine proteases. Protein Sci. 1995, 4: 337-360.PubMed CentralPubMedView Article
- Toolbox at the EBI European Bioinformatics. [http://www.ebi.ac.uk/clustalw/]
- Jiang H, Kanost MR: The clip-domain family of serine proteinases in arthropods. Insect Biochem Mol Biol. 2000, 30: 95-105. 10.1016/S0965-1748(99)00113-7.PubMedView Article
- De Gregorio E, Spellman PT, Rubin GM, Lemaitre B: Genome-wide analysis of the Drosophila immune response by using oligonucleotide microarrays. Proc Natl Acad Sci USA. 2001, 98: 12590-12595. 10.1073/pnas.221458698.PubMed CentralPubMedView Article
- De Gregorio E, Spellman PT, Tzou P, Rubin GM, Lemaitre B: The Toll and Imd pathways are the major regulators of the immune response in Drosophila. EMBO Journal. 2002, 21: 2568-2579. 10.1093/emboj/21.11.2568.PubMed CentralPubMedView Article
- Lemaitre B, Reichhart J-M, Hoffmann JA: Drosophila host defense: Differential induction of antimicrobial peptide genes after infection by various classes of microorganisms. Proc Natl Acad Sci USA. 1997, 94: 14614-14619. 10.1073/pnas.94.26.14614.PubMed CentralPubMedView Article
- Jang I-H, Chosa N, Kim S-H, Nam H-J, Lemaitre B, Ochiai M, Kambris Z, Brun S, Hashimoto C, Ashida M, et al: Spätzle-Processing Enzyme Required for Toll Signaling Activation in Drosophila Innate Immunity. Developmental Cell. 2006, 10: 45-55. 10.1016/j.devcel.2005.11.013.PubMedView Article
- Tang H, Kambris Z, Lemaitre B, Hashimoto C: Two Proteases Defining a Melanization Cascade in the Immune System of Drosophila. J Biol Chem. 2006, 281: 28097-28104. 10.1074/jbc.M601642200.PubMedView Article
- Castillejo-López C, Häcker U: The serine protease Sp7 is expressed in blood cells and regulates the melanization reaction in Drosophila. Biochemical and Biophysical Research Communications. 2005, 338: 1075-1082. 10.1016/j.bbrc.2005.10.042.PubMedView Article
- Scherfer C, Qazi MR, Takahashi K, Ueda R, Dushay MS, Theopold U, Lemaitre B: The Toll immune-regulated Drosophila protein Fondue is involved in hemolymph clotting and puparium formation. Developmental Biology. 2006, 295: 156-163. 10.1016/j.ydbio.2006.03.019.PubMedView Article
- Kambris Z, Brun S, Jang I-H, Nam I-J, Romeo Y, Takahashi K, Lee W-J, Ueda R, Lemaitre B: Drosophila Immunity: A Large-Scale In Vivo RNAi Screen Identifies Five Serine Proteases Required for Toll Activation. Current Biology. 2006, 16: 808-813. 10.1016/j.cub.2006.03.020.PubMedView Article
- Kalfayan L, Wensink PC: Developmental regulation of Drosophila alpha-tubulin genes. Cell. 1982, 29: 91-98. 10.1016/0092-8674(82)90093-9.PubMedView Article
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.PubMedView Article
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView Article
- Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci USA. 1988, 85: 2444-2448. 10.1073/pnas.85.8.2444.PubMed CentralPubMedView Article
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView Article
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 376-382. 10.1101/gr.9.9.868.View Article
- Sambrook J, Fitsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Spring Harbor, Cold Spring Harbor Press
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.