Giles Goetz

A genetic basis for the phenotypic differentiation between siscowet and lean lake trout (Salvelinus namaycush)

In Lake Superior there are three principal forms of lake trout (Salvelinus namaycush): lean, siscowet and humper. Wild lean and siscowet differ in the shape and relative size of the head, size of the fins, location and size of the eyes, caudal peduncle shape and lipid content of the musculature. To investigate the basis for these phenotypic differences, lean and siscowet lake trout, derived from gametes of wild populations in Lake Superior, were reared communally under identical environmental conditions for 2.5 years. Fish were analysed for growth, morphometry and lipid content, and differences in liver transcriptomics were investigated using Roche 454 GS‐FLX pyrosequencing. The results demonstrate that key phenotypic differences between wild lean and siscowet lake trout such as condition factor, morphometry and lipid levels, persist in these two forms when reared in the laboratory under identical environmental conditions. This strongly suggests that these differences are genetic and not a result of environmental plasticity. Transcriptomic analysis involving the comparison of hepatic gene frequencies (RNA‐seq) and expression (quantitative reverse transcription–polymerase chain reaction (qPCR)) between the two lake trout forms, indicated two primary gene groups that were differentially expressed; those involving lipid synthesis, metabolism and transport (acyl‐CoA desaturase, acyl‐CoA binding protein, peroxisome proliferator‐activated receptor gamma, and apolipoproteins), and those involved with immunity (complement component C3, proteasome, FK506 binding protein 5 and C1q proteins). The results demonstrate that RNA‐seq can be used to identify differentially expressed genes; however, some discrepancies between RNA‐seq analysis and qPCR indicate that methods for deep sequencing may need to be refined and/or different RNA‐seq platforms utilized.

Temporal and spatial variability in nearshore bacterioplankton communities of Lake Michigan.

The spatial and temporal variability of bacterial communities were determined for the nearshore waters of Lake Michigan, an oligotrophic freshwater inland sea. A freshwater estuary and nearshore sites were compared six times during 2006 using denaturing gradient gel electrophoresis (DGGE). Bacterial composition clustered by individual site and date rather than by depth. Seven 16S rRNA gene clone libraries were constructed, yielding 2717 bacterial sequences. Spatial variability was detected among the DGGE banding patterns and supported by clone library composition. The clone libraries from deep waters and the estuary environment revealed highest overall bacterial diversity. Betaproteobacteria sequence types were the most dominant taxa, comprising 40.2-67.7% of the clone libraries. BAL 47 was the most abundant freshwater cluster of Betaproteobacteria, indicating widespread distribution of this cluster in the nearshore waters of Lake Michigan. Incertae sedis 5 and Oxalobacteraceae sequence types were prevalent in each clone library, displaying more diversity than previously described in other freshwater environments. Among the Oxalobacteraceae sequences, a globally distributed freshwater cluster was determined. The nearshore waters of Lake Michigan are a dynamic environment that experience forces similar to the coastal ocean environment and share common bacterial diversity with other freshwater habitats.

RNA-Seq Reveals an Integrated Immune Response in Nucleated Erythrocytes

Background Throughout the primary literature and within textbooks, the erythrocyte has been tacitly accepted to have maintained a unique physiological role; namely gas transport and exchange. In non-mammalian vertebrates, nucleated erythrocytes are present in circulation throughout the life cycle and a fragmented series of observations in mammals support a potential role in non-respiratory biological processes. We hypothesised that nucleated erythrocytes could actively participate via ligand-induced transcriptional re-programming in the immune response. Methodology/Principal Findings Nucleated erythrocytes from both fish and birds express and regulate specific pattern recognition receptor (PRR) mRNAs and, thus, are capable of specific pathogen associated molecular pattern (PAMP) detection that is central to the innate immune response. In vitro challenge with diverse PAMPs led to de novo specific mRNA synthesis of both receptors and response factors including interferon-alpha (IFNα) that exhibit a stimulus-specific polysomal shift supporting active translation. RNA-Seq analysis of the PAMP (Poly (I∶C), polyinosinic∶polycytidylic acid)-erythrocyte response uncovered diverse cohorts of differentially expressed mRNA transcripts related to multiple physiological systems including the endocrine, reproductive and immune. Moreover, erythrocyte-derived conditioned mediums induced a type-1 interferon response in macrophages thus supporting an integrative role for the erythrocytes in the immune response. Conclusions/Significance We demonstrate that nucleated erythrocytes in non-mammalian vertebrates spanning significant phylogenetic distance participate in the immune response. RNA-Seq studies highlight a mRNA repertoire that suggests a previously unrecognized integrative role for the erythrocytes in other physiological systems.

Bacteroidales Diversity in Ring-Billed Gulls (Laurus delawarensis) Residing at Lake Michigan Beaches

ABSTRACT This study investigated the occurrence and diversity of Bacteroidales fecal bacteria in gulls residing in the Great Lakes region. Members of this bacterial order have been widely employed as human and bovine host-specific markers of fecal pollution; however, few studies have focused on gulls, which can be a major source of fecal indicator bacteria and pathogens at beaches. We found a low but consistent occurrence of Bacteroidales in gulls at five beaches in three different counties spanning the Wisconsin shoreline of Lake Michigan. The percentages of gulls positive for Bacteroidales were 4 to 8% at beaches in the southern part of the state and 8 to 50% at beaches in the north. Sequencing of 931 clones from seven gull Bacteroidales 16S rRNA gene libraries revealed a large amount of diversity in both individual and pooled gull fecal samples. Two libraries constructed from pooled gull fecal samples (n = 5 and n = 6) did not have a greater richness of sequences than individual samples, suggesting that even within a single gull diversity is high and an extensive sequencing effort is needed to characterize the populations. Estimates of the numbers of operational taxonomic units (OTUs) for the libraries obtained using different similarity levels revealed a large amount of microdiveristy with a limited number of OTUs at the 95% similarity level. Gull sequences were clustered by the beach from which they were collected, suggesting that there were geographic effects on the distribution of Bacteriodales. More than 53% of the 16S rRNA gene sequences from gulls at the southern beaches were associated with the family Porphyromonadaceae, primarily the genus Parabacteroides, whereas sequences from gulls at the northern beaches were comprised of Bacteroidaceae and Prevotellaceae sequences. Comparison of gull sequences with sequences from goose, canine, raccoon, and sewage sources revealed distinct clusters of closely related gull sequences; however, these sequences were widely dispersed across a dendrogram that included all other sources, including previously characterized gull Bacteroidales from other studies, suggesting that geographic influence or simply sample representation plays a greater role in the observed population structure than strictly the host gut environment.

High-throughput sequencing and pathway analysis reveal alteration of the pituitary transcriptome by 17α-ethynylestradiol (EE2) in female coho salmon, Oncorhynchus kisutch.

Considerable research has been done on the effects of endocrine disrupting chemicals (EDCs) on reproduction and gene expression in the brain, liver and gonads of teleost fish, but information on impacts to the pituitary gland are still limited despite its central role in regulating reproduction. The aim of this study was to further our understanding of the potential effects of natural and synthetic estrogens on the brain-pituitary-gonad axis in fish by determining the effects of 17α-ethynylestradiol (EE2) on the pituitary transcriptome. We exposed sub-adult coho salmon (Oncorhynchus kisutch) to 0 or 12 ng EE2/L for up to 6 weeks and effects on the pituitary transcriptome of females were assessed using high-throughput Illumina(®) sequencing, RNA-Seq and pathway analysis. After 1 or 6 weeks, 218 and 670 contiguous sequences (contigs) respectively, were differentially expressed in pituitaries of EE2-exposed fish relative to control. Two of the most highly up- and down-regulated contigs were luteinizing hormone β subunit (241-fold and 395-fold at 1 and 6 weeks, respectively) and follicle-stimulating hormone β subunit (-3.4-fold at 6 weeks). Additional contigs related to gonadotropin synthesis and release were differentially expressed in EE2-exposed fish relative to controls. These included contigs involved in gonadotropin releasing hormone (GNRH) and transforming growth factor-β signaling. There was an over-representation of significantly affected contigs in 33 and 18 canonical pathways at 1 and 6 weeks, respectively, including circadian rhythm signaling, calcium signaling, peroxisome proliferator-activated receptor (PPAR) signaling, PPARα/retinoid x receptor α activation, and netrin signaling. Network analysis identified potential interactions between genes involved in circadian rhythm and GNRH signaling, suggesting possible effects of EE2 on timing of reproductive events.

Stimulation of growth and changes in the hepatic transcriptome by 17β-estradiol in the yellow perch (Perca flavescens)

The effects of dietary 17beta-estradiol (E(2)) on growth and liver transcriptomics were investigated in the yellow perch (Perca flavescens). After a 3-mo treatment, E(2) significantly stimulated an increase in length and weight of juvenile male and female perch relative to control animals. The increase was significantly greater in females compared with males. Separate, unnormalized cDNA libraries were constructed from equal quantities of RNA from 6 male and 6 female livers of E(2)-treated and control perch, and 3,546 and 3,719 expressed sequence tags (ESTs) were obtained, respectively. To characterize E(2)-regulated transcripts, EST frequencies between libraries were calculated within contiguous sequences that were assembled from the combined ESTs of both libraries. Frequencies were also determined in EST transcript groupings produced by aligning all of the ESTs from both libraries at the nucleotide level. From these analyses, there were 28 annotated transcripts that were regulated by 75% between libraries and for which there were at least 5 ESTs of the same transcript between libraries. Regulation of a subset (14) of these transcripts was confirmed by quantitative reverse transcription-polymerase chain reaction (QPCR). Transcripts that were upregulated by E(2) included reproduction-related proteins, binding proteins, and proteases and protease inhibitors. While not part of the transcript frequency analysis, QPCR showed significant upregulation of estrogen receptor esr1 and of insulin-like growth factor I (IGF-I) in E(2) livers. E(2)-downregulated transcripts represented a variety of functional categories including components of the respiratory chain, lipid transport and metabolism, glycolysis, amino acid and nitrogen metabolism, binding proteins, a hydrolytic enzyme, and a transcriptional regulator. In perch it appears that exogenous estrogen drastically shifts liver metabolism toward the production of lipoproteins and carbohydrate binding proteins, and that the growth-promoting action may involve an increase in hepatic IGF-I production.

Pathogen-associated gene expression profiles in rainbow trout macrophages.

Pathogens can be distinctively recognized by the cells of the immune system through interactions between the Pathogen-Associated Molecular Patterns (PAMPs) that they produce and the innate immune receptors of leukocytes. The present paper reports on the PAMP-modulated expression of a group of genes expressed in trout macrophages. The genes were identified in subtracted libraries from lipopolysaccharide (LPS)-stimulated macrophages and their expression was analyzed using quantitative real time PCR following stimulation of the cells with E. coli LPS, poly (I:C) and zymosan; representing Gram-negative bacteria, viruses and fungi, respectively. Genes (SPINT1L, DDIT4L, STEAP4, and TNFAIP3), the expression of which was induced by LPS and zymosan, were not significantly up-regulated by poly(I:C) and the opposite was found for transcripts (HMGB1 and PSMB9) up-regulated by poly(I:C). Overall, the differences in gene expression were greater at a later stage of macrophage activation (24 h) at a time when stimulation with poly(I:C) resulted in substantially different responses as compared to LPS and, to a lesser extent, zymosan.

Gene expression and pathologic alterations in juvenile rainbow trout due to chronic dietary TCDD exposure.

The goal of this project was to use functional genomic methods to identify molecular biomarkers as indicators of the impact of TCDD exposure in rainbow trout. Specifically, we investigated the effects of chronic dietary TCDD exposure on whole juvenile rainbow trout global gene expression associated with histopathological analysis. Juvenile rainbow trout were fed Biodiet starter with TCDD added at 0, 0.1, 1, 10 and 100 ppb (ngTCDD/g food), and fish were sampled from each group at 7, 14, 28 and 42 days after initiation of feeding. 100 ppb TCDD caused 100% mortality at 39 days. Fish fed with 100 ppb TCDD food had TCDD accumulation of 47.37 ppb (ngTCDD/g fish) in whole fish at 28 days. Histological analysis from TCDD-treated trout sampled from 28 and 42 days revealed that obvious lesions were found in skin, oropharynx, liver, gas bladder, intestine, pancreas, nose and kidney. In addition, TCDD caused anemia in peripheral blood, decreases in abdominal fat, increases of remodeling of fin rays, edema in pericardium and retrobulbar hemorrhage in the 100 ppb TCDD-treated rainbow trout compared to the control group at 28 days. Dose- and time-dependent global gene expression analyses were performed using the cGRASP 16,000 (16K) cDNA microarray. TCDD-responsive whole body transcripts identified in the microarray experiments have putative functions involved in various biological processes including growth, cell proliferation, metabolic process, and immune system processes. Nine microarray-identified genes were selected for QPCR validation. CYP1A3 and CYP1A1 were common up-regulated genes and HBB1 was a common down-regulated gene among each group based on microarray data, and their QPCR validations are consistent with microarray data for the 10 and 100 ppb TCDD treatment groups after 28 days exposure (p<0.05). In addition, in the 100 ppb group at 28 days, expression of complement component C3-1 and trypsin-1 precursor have a more than 10-fold induction from the microarray experiments, and their QPCR validations are consistent and showed significant induction in the 100 ppb group at 28 days (p<0.05). Overall, lesion in nasal epithelium is a novel and significant result in this study, and TCDD-responsive rainbow trout transcripts identified in the present study may lead to the development of new molecular biomarkers for assessing the potential impacts of environmental TCDD on rainbow trout populations.

Histopathologic Alterations Associated with Global Gene Expression Due to Chronic Dietary TCDD Exposure in Juvenile Zebrafish

The goal of this project was to investigate the effects and possible developmental disease implication of chronic dietary TCDD exposure on global gene expression anchored to histopathologic analysis in juvenile zebrafish by functional genomic, histopathologic and analytic chemistry methods. Specifically, juvenile zebrafish were fed Biodiet starter with TCDD added at 0, 0.1, 1, 10 and 100 ppb, and fish were sampled following 0, 7, 14, 28 and 42 d after initiation of the exposure. TCDD accumulated in a dose- and time-dependent manner and 100 ppb TCDD caused TCDD accumulation in female (15.49 ppb) and male (18.04 ppb) fish at 28 d post exposure. Dietary TCDD caused multiple lesions in liver, kidney, intestine and ovary of zebrafish and functional dysregulation such as depletion of glycogen in liver, retrobulbar edema, degeneration of nasal neurosensory epithelium, underdevelopment of intestine, and diminution in the fraction of ovarian follicles containing vitellogenic oocytes. Importantly, lesions in nasal epithelium and evidence of endocrine disruption based on alternatively spliced vasa transcripts are two novel and significant results of this study. Microarray gene expression analysis comparing vehicle control to dietary TCDD revealed dysregulated genes involved in pathways associated with cardiac necrosis/cell death, cardiac fibrosis, renal necrosis/cell death and liver necrosis/cell death. These baseline toxicological effects provide evidence for the potential mechanisms of developmental dysfunctions induced by TCDD and vasa as a biomarker for ovarian developmental disruption.

An Enriched European Eel Transcriptome Sheds Light upon Host-Pathogen Interactions with Vibrio vulnificus

Infectious diseases are one of the principal bottlenecks for the European eel recovery. The aim of this study was to develop a new molecular tool to be used in host-pathogen interaction experiments in the eel. To this end, we first stimulated adult eels with different pathogen-associated molecular patterns (PAMPs), extracted RNA from the immune-related tissues and sequenced the transcriptome. We obtained more than 2x106 reads that were assembled and annotated into 45,067 new descriptions with a notable representation of novel transcripts related with pathogen recognition, signal transduction and the immune response. Then, we designed a DNA-microarray that was used to analyze the early immune response against Vibrio vulnificus, a septicemic pathogen that uses the gills as the portal of entry into the blood, as well as the role of the main toxin of this species (RtxA13) on this early interaction. The gill transcriptomic profiles obtained after bath infecting eels with the wild type strain or with a mutant deficient in rtxA13 were analyzed and compared. Results demonstrate that eels react rapidly and locally against the pathogen and that this immune-response is rtxA13-dependent as transcripts related with cell destruction were highly up-regulated only in the gills from eels infected with the wild-type strain. Furthermore, significant differences in the immune response against the wild type and the mutant strain also suggest that host survival after V. vulnificus infection could depend on an efficient local phagocytic activity. Finally, we also found evidence of the presence of an interbranchial lymphoid tissue in European eel gills although further experiments will be necessary to identify such tissue.