Transcriptome of the Eastern Oyster Crassostrea Virginica in Response to Bacterial Challenge
Total Page:16
File Type:pdf, Size:1020Kb
University of Rhode Island DigitalCommons@URI Open Access Master's Theses 2013 Transcriptome of the Eastern Oyster Crassostrea Virginica in Response to Bacterial Challenge Ian McDowell University of Rhode Island, [email protected] Follow this and additional works at: https://digitalcommons.uri.edu/theses Recommended Citation McDowell, Ian, "Transcriptome of the Eastern Oyster Crassostrea Virginica in Response to Bacterial Challenge" (2013). Open Access Master's Theses. Paper 19. https://digitalcommons.uri.edu/theses/19 This Thesis is brought to you for free and open access by DigitalCommons@URI. It has been accepted for inclusion in Open Access Master's Theses by an authorized administrator of DigitalCommons@URI. For more information, please contact [email protected]. TRANSCRIPTOME OF THE EASTERN OYSTER CRASSOSTREA VIRGINICA IN RESPONSE TO BACTERIAL CHALLENGE BY IAN MCDOWELL A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN BIOLOGICAL AND ENVIRONMENTAL SCIENCES UNIVERSITY OF RHODE ISLAND 2013 ! MASTER OF SCIENCE THESIS OF IAN MCDOWELL APPROVED: Thesis Committee: Major Professor Marta Gomez-Chiarri Chris E. Lane Matt P. Hare Jason J. Kolbe Nasser H. Zawia DEAN OF THE GRADUATE SCHOOL ! UNIVERSITY OF RHODE ISLAND 2013 ABSTRACT The Eastern oyster Crassostrea virginica, an ecologically and economically important estuarine organism, suffers mortalities as high as 90-100% in affected areas due to Roseovarius Oyster Disease (ROD), caused by the bacterial pathogen, Roseovarius crassostreae. Advanced genotypic breeding techniques necessitate information regarding markers and genes associated with disease resistance. As yet, the host-pathogen interaction between C. virginica and R. crassostreae is poorly understood at the molecular level. The identification of potential genes and pathways responsible for an effective host defense response in the Eastern oyster to R. crassostreae is important not only to provide a basis for enhanced breeding techniques, but also to enhance understanding of innate immunity in a broader, evolutionary sense. The present study proposed to uncover not only genes and general processes potentially involved in disease resistance to ROD in the Eastern oyster, but also diversified gene families. To that end the present study entailed a disease challenge exposing ROD- resistant and ROD-susceptible families of oysters to R. crassostreae, high- throughput cDNA sequencing of samples from several timepoints throughout the disease challenge, assembly of sequence data into a reference transcriptome, analysis of the transcriptome through differential gene expression and gene family similarity clustering, and single nucleotide polymorphism (SNP) detection to identify candidate gene markers. Oyster resistance to R. crassostreae was found to involve extracellular matrix remodeling, cell adhesion, inflammation, metabolism, and other processes. Several gene families identified as putatively diversified and important in the oyster host defense response were enumerated and described, including serine proteases, serine protease inhibitors, c-type lectins, C1q domain-containing proteins, fibrinogen domain-containing proteins, scavenger receptors with class B SRCR domains, interferon-induced protein 44 (IFI44) family proteins, and GTPase of the immunity associated protein (GIMAP) family proteins. Further, similarity clustering of proteins and translated transcripts from diverse invertebrates suggested that GIMAP proteins are expanded in molluscs and IFI44 proteins are expanded in bivalves. ACKNOWLEDGEMENTS I thank my major professor, Marta Gomez-Chiarri, for her guidance, patience and expertise. I thank all of my labmates who have helped with this project and have helped to guide me in the lab including Chamilani Nikapitiya, Jess Piesz, Murni Karim, and Sae Bohm Sohn. I thank Ian Misner and Chris Lane who helped me greatly on the computational aspects of this project. I also thank the staff at Brown University’s Center for Computing and Visualization, especially Lingsheng Dong and Mark Howison, for their timely and immensely helpful support. I thank my family, friends, and very supportive girlfriend, Stacey McMillin. "#! TABLE OF CONTENTS ! ABSTRACT!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!""! ACKNOWLEDGEMENTS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#! TABLE OF CONTENTS!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!#! LIST OF TABLES !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!#"! LIST OF FIGURES !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!#"""! INTRODUCTION !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!$ METHODS!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!%! RESULTS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! &'! DISCUSSION !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ()! TABLES!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! '*! FIGURES !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! +'! APPENDIX !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *%, BIBLIOGRAPHY !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!$)-! , #! LIST OF TABLES Table 1. Assembly metrics for transcriptome assembly$$$$$$$$$$%59 Table 2. Top 50 most highly differentially expressed genes (DEGs) unique to GX_early up-regulated relative to control pool (CGX_late)$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$60 Table 3. Differentially expressed genes (DEGs) unique to GX_early down- regulated relative to control pool$$$$$$$$$$$%%%$$$$$$$$$$$$&' Table 4. Top 50 most highly differentially expressed genes (DEGs) among GX_early shared with F3L_early and/or F3L_late up-regulated relative to control pool $$$$$$$$$$$$$$$$$$$$$$$$%%$$$$$$$$$$$66 Table 5. Top 50 most highly differentially expressed genes (DEGs) among GX_early DEGs shared with F3L_early and/or F3L_late down-regulated relative to control pool $$$$$$$$$$$$$$%$$$$$%$$$$$$$$$%%%69 Table 6. Pfam families enriched among the differentially expressed genes (DEGs) from resistant (GX) and susceptible (F3L) oysters at early (1 and 5 d) and late (15 and 30 d) timepoints compared to control (CGX) oysters, with rank order of significance of enrichment and significance of enrichment (p-value)%%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$72 Table 7. Pfam domains enriched among the differentially expressed genes #"! (DEGs) from resistant (GX) and susceptible (F3L) oysters at early (1 and 5 d) and late (15 and 30 d) timepoints compared to control (CGX) oysters, with rank order of significance of enrichment and significance of enrichment (p-value)$$$$$$$$$$$$$$$$$$$$$$$%%$$$$$$..74 Table 8. TribeMCL clusters enriched for differentially expressed genes (DEGs) from resistant (GX) and susceptible (F3L) oysters at early (1 and 5 d) and late (15 and 30 d) timepoints, with rank order of significance of enrichment and significance of enrichment (p-value)$$$$$%%$$$$%%$$76 Table 9. Comparison of select “diversified” groups of genes in the Eastern oyster to the Pacific oyster C.gigas genome and differential expression$$$$$$$$$$$$$$$$$$$$..$$$$$$$%$$$$$$$%$() #""! LIST OF FIGURES Figure 1. Cumulative percent mortality in two families of oysters challenged with the bacterial pathogen Roseovarius crassostreae (F3L and GX) and in unchallenged controls (CF3L and CGX).$$$$$$$$$%85 Figure 2. Heatmap of all differentially expressed genes (DEGs) in two oyster families experimentally challenged with the bacterial pathogen Roseovarius crassostreae (GX and F3L) and unchallenged controls$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$%$$%%$..86 Figure 3. Numbers of differentially expressed contigs (DEGs) shared and unique between GX_early (resistant oysters – days 1 and 5), F3L_early (susceptible oysters – days 1 and 5), and F3L_late (susceptible oysters – days 15 and 30)%%$$$$$$$$$$$$$$%$$$$$$$$$$$%%$$$$$$.87 Figure 4. Functionally enriched Gene Ontology terms among DEGs from resistant (GX) and susceptible (F3L) oysters at early (1 and 5 d) and late (15 and 30 d) time points compared to control (CGX) oysters mapped by semantic similarity using SimRel method and REViGO$$$$$$$$%$$88 Figure 5. Numbers of non-redundant transcripts/proteins and gene models in selected putative diversified gene families in multiple organisms from diverse taxa!$$$$%$$$$$$$$$$$$%%%$$$$$$$$$$$$$$$$$..90 #"""! Figure 6. C. virginica GIMAP translated transcripts clustered by similarity with color and size reflective of differential expression$$$$$$$$%%$$.91 Figure 7. C. virginica IFI44 translated transcripts clustered by similarity with color and size reflective of differential