This Is Normal Text
Total Page:16
File Type:pdf, Size:1020Kb
A COMPARATIVE ANALYSIS OF GENE EXPRESSION AMONG CASTES OF THE TERMITE RETICULITERMES FLAVIPES USING EXPRESSED SEQUENCE TAGS (ESTS) AND A MICROARRAY by MATTHEW MICHAEL STELLER B.S., UNIVERSITY OF MINNESOTA-DULUTH, 2005 A THESIS submitted in partial fulfillment of the requirements for the degree MASTER OF SCIENCE Department of Entomology College of Agriculture KANSAS STATE UNIVERSITY Manhattan, Kansas 2009 Approved by: Major Professor Srini Kambhampati Copyright MATTHEW MICHEAL STELLER 2009 Abstract Termites (Isoptera) are separated into morphologically and behaviorally specialized castes of sterile workers and soldiers, and the reproductive alates. Previous research on eusocial insects has indicated that caste differentiation has a genetic basis. Although much has been studied about the genetic basis of caste differentiation and behavior in the honey bee, Apis mellifera, termites remain comparatively understudied. Therefore, my objective was to compare the gene expression patterns of different castes of the termite Reticulitermes flavipes based on EST analyses and a microarray. Soldier, worker, and alate caste and two larval life stage cDNA libraries were constructed, and ~15,000 randomly chosen clones were sequenced to compile an EST database. Putative gene functions were assigned based on a BLASTX Swissprot search. Categorical expression patterns for each library were compared using the in silico methods of BLAST2GO and r-statistics. I chose 2,240 unique-ESTs based on their putative function and sequence quality, which I used to fabricate a Combimatrix microarray. I used the microarray to compare expression levels between workers and soldiers from Kansas and Florida populations. Seventy to ninety percent of the sequences from the ESTs of each caste and life stages had no significant similarity to those in existing databases. All libraries contained sequences with putative reproductive functions, which was unexpected in the non-reproductive soldier and worker castes. Sequences of interest that showed a putative bias among castes include a viral protein in soldiers and a possible chemosensory protein in alates, which may be involved in termite reproductive functionality or communication. The microarray showed increased expression in the soldier caste of a sequence that matched tropomyosin and an increased expression level of a sequence that matched a PDZ-domain containing protein in some worker samples. This study leads to several candidate genes of potential caste specific function, which can be further tested using functional analysis and between the different castes and life stages of R. flavipes. These genes include the sequences similar to pebIII and RGS-GAIP. I have also expanded upon the available sequences for this termite and utilized the r-statistic in silico method for the first time to putatively compare gene expression in the different castes of a eusocial insect. The in silico analysis allowed us to identify several genes which may show biased expression patterns in the different cDNA libraries and which may reveal caste-specific expression controlling important functions after further analysis. These candidates include: an alate-biased gene, which had a predicted function of neurotransmitter secretion and cholesterol absorption as well as a late larval-biased gene which was predicted to be involved in protein biosynthesis and ligase activity. Table of Contents List of Figures ............................................................................................................................... vii List of Tables ................................................................................................................................. ix Acknowledgements ........................................................................................................................ xi CHAPTER 1 - Analysis of Caste-Specific Gene Expression Patterns ........................................... 1 INTRODUCTION ....................................................................................................................... 1 MATERIALS AND METHODS ................................................................................................. 4 Termites ................................................................................................................................... 4 cDNA Libraries ........................................................................................................................ 4 Sequence Analyses ................................................................................................................... 5 In silico transcript abundance .................................................................................................. 5 Microarray Construction and Gene Selection .......................................................................... 6 Total RNA isolation, Labeled-cDNA synthesis, and Microarray Hybridization ..................... 7 Microarray Imaging ................................................................................................................. 8 Quantitative PCR for Validation of Microarray Results .......................................................... 8 RESULTS .................................................................................................................................. 10 General Library Analysis ....................................................................................................... 10 Sequence Analysis ................................................................................................................. 10 In Silico Analysis- Significant BLAST Hits .......................................................................... 11 In Silico Analysis- Protein Prediction ................................................................................... 11 Microarrays Combined .......................................................................................................... 12 Microarrays of Florida Samples ............................................................................................. 13 Quantitative PCR for Validation of Microarray Results ........................................................ 13 DISCUSSION ............................................................................................................................ 14 General Library Analysis ....................................................................................................... 14 In Silico Analysis ................................................................................................................... 15 Microarrays of Combined Samples ....................................................................................... 18 Microarray Analysis of Florida Samples ............................................................................... 20 CONCLUSION .......................................................................................................................... 20 v Figures and Tables ..................................................................................................................... 22 References ..................................................................................................................................... 66 Appendix -A Total BLAST Results .............................................................................................. 72 Total Microarray BLAST .......................................................................................................... 72 vi List of Figures Figure 1.1 Proportions of GO cellular component terms (level two) in each EST library sequences annotated by Blast2GO A) Soldiers; B) Alates; C) Workers; D) Early Larval; E) Late Larval ........................................................................................................ 22 Figure 1.2 Proportions of GO biological process terms (level two) of each EST library annotated by Blast2GO A) Soldiers; B) Alates; C) Workers; D) Early Larval; E) Late Larval .................................................................................................................................... 27 Figure 1.3 Proportional analysis of GO molecular function terms (level two) of each individual EST library annotated by Blast2GO A) Soldiers; B) Alates; C) Workers; D) Early Larval; E) Late Larval .............................................................................................. 32 Figure 1.4 Frequency of putatively reproductive GO terms in each EST library as a fraction to total number of unigenes. ................................................................................ 37 Figure 1.5 Measure of statistical significance, resultant significant r-values were validated using a chi-squared analysis. R-statistic measure of relative entropy and significance was done as in Stekel et al. 2000. ........................................................................................ 38 Figure 1.6 R-statistic comparison of contig composistion and alate singleton bias of contig545 of predicted neurotransmitter secretion or cholesterol absorption functioning. ........................................................................................................................... 39 Figure 1.7 R-statistic comparison of contig composistion and late larval singleton bias of contig366 of predicted protein biosynthesis or ligase activity functionality. .................. 40 Figure 1.8 R-statistic comparison of contig composistion