Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish.

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

Strong population differentiation in lingcod (Ophiodon elongatus) is driven by a small portion of the genome.

Delimiting intraspecific genetic variation in harvested species is crucial to the assessment of population status for natural resource management and conservation purposes. Here, we evaluated genetic population structure in lingcod (Ophiodon elongatus), a commercially and recreationally important fishery species along the west coast of North America. We used 16,749 restriction site-associated DNA sequencing (RADseq) markers, in 611 individuals collected from across the bulk of the species range from Southeast Alaska to Baja California, Mexico. In contrast to previous population genetic work on this species, we found strong evidence for two distinct genetic clusters. These groups separated latitudinally with a break near Point Reyes off Northern California, and there was a high frequency of admixed individuals in close proximity to the break. F-statistics corroborate this genetic break between northern and southern sampling sites, although most loci are characterized by low FST values, suggesting high gene flow throughout most of the genome. Outlier analyses identified 182 loci putatively under divergent selection, most of which mapped to a single genomic region. When individuals were grouped by cluster assignment (northern, southern, and admixed), 71 loci were fixed between the northern and southern cluster, all of which were identified in the outlier scans. All individuals identified as admixed exhibited near 50:50 assignment to northern and southern clusters and were heterozygous for most fixed loci. Alignments of RADseq loci to a draft lingcod genome assembly and three other teleost genomes with chromosome-level assemblies suggest that outlier and fixed loci are concentrated on a single chromosome. Similar genomic patterns have been attributed to chromosomal inversions in diverse taxonomic groups. Regardless of the evolutionary mechanism, these results represent novel observations of genetic structure in lingcod and designate clear evolutionary units that could be used to inform fisheries management.

Cardiac remodeling in response to embryonic crude oil exposure involves unconventional NKX family members and innate immunity genes.

Cardiac remodeling results from both physiological and pathological stimuli. Compared with mammalian hearts, fish hearts show a broader array of remodeling changes in response to environmental influences, providing exceptional models for dissecting the molecular and cellular bases of cardiac remodeling. We recently characterized a form of pathological remodeling in juvenile pink salmon (Oncorhynchus gorbuscha) in response to crude oil exposure during embryonic cardiogenesis. In the absence of overt pathology (cardiomyocyte death or inflammatory infiltrate), cardiac ventricles in exposed fish showed altered shape, reduced thickness of compact myocardium and hypertrophic changes in spongy, trabeculated myocardium. Here, we used RNA sequencing to characterize molecular pathways underlying these defects. In juvenile ventricular cardiomyocytes, antecedent embryonic oil exposure led to dose-dependent upregulation of genes involved in innate immunity and two NKX homeobox transcription factors not previously associated with cardiomyocytes, nkx2.3 and nkx3.3 Absent from mammalian genomes, the latter is largely uncharacterized. In zebrafish embryos, nkx3.3 demonstrated a potent effect on cardiac morphogenesis, equivalent to that of nkx2.5, the primary transcription factor associated with ventricular cardiomyocyte identity. The role of nkx3.3 in heart growth is potentially linked to the unique regenerative capacity of fish and amphibians. Moreover, these findings support a cardiomyocyte-intrinsic role for innate immune response genes in pathological hypertrophy. This study demonstrates how an expanding mechanistic understanding of environmental pollution impacts - i.e. the chemical perturbation of biological systems - can ultimately yield new insights into fundamental biological processes.

Temporal Dynamics of DNA Methylation Patterns in Response to Rearing Juvenile Steelhead (Oncorhynchus mykiss) in a Hatchery versus Simulated Stream Environment.

Genetic selection is often implicated as the underlying cause of heritable phenotypic differences between hatchery and wild populations of steelhead trout (Oncorhynchus mykiss) that also differ in lifetime fitness. Developmental plasticity, which can also affect fitness, may be mediated by epigenetic mechanisms such as DNA methylation. Our previous study identified significant differences in DNA methylation between adult hatchery- and natural-origin steelhead from the same population that could not be distinguished by DNA sequence variation. In the current study, we tested whether hatchery-rearing conditions can influence patterns of DNA methylation in steelhead with known genetic backgrounds, and assessed the stability of these changes over time. Eyed-embryos from 22 families of Methow River steelhead were split across traditional hatchery tanks or a simulated stream-rearing environment for 8 months, followed by a second year in a common hatchery tank environment. Family assignments were made using a genetic parentage analysis to account for relatedness among individuals. DNA methylation patterns were examined in the liver, a relatively homogeneous organ that regulates metabolic processes and somatic growth, of juveniles at two time points: after eight months of rearing in either a tank or stream environment and after a subsequent year of rearing in a common tank environment. Further, we analyzed DNA methylation in the sperm of mature 2-year-old males from the earlier described treatments to assess the potential of environmentally-induced changes to be passed to offspring. Hepatic DNA methylation changes in response to hatchery versus stream-rearing in yearling fish were substantial, but few persisted after a second year in the tank environment. However, the early rearing environment appeared to affect how fish responded to developmental and environmental signals during the second year since novel DNA methylation differences were identified in the livers of hatchery versus stream-reared fish after a year of common tank rearing. Furthermore, we found profound differences in DNA methylation due to age, irrespective of rearing treatment. This could be due to smoltification associated changes in liver physiology after the second year of rearing. Although few rearing-treatment effects were observed in the sperm methylome, strong family effects were observed. These data suggest limited potential for intergenerational changes, but highlight the importance of understanding the effects of kinship among studied individuals in order to properly analyze and interpret DNA methylation data in natural populations. Our work is the first to study family effects and temporal dynamics of DNA methylation patterns in response to hatchery-rearing.

Initiation of sex change and gonadal gene expression in black sea bass (Centropristis striata) exposed to exemestane, an aromatase inhibitor.

Many teleost fishes exhibit sequential hermaphroditism, where male or female gonads develop first and later undergo sex change. Model sex change species are characterized by social hierarchies and coloration changes, which enable experimental manipulations to better understand these processes. However, other species such as the protogynous black sea bass (Centropristis striata) do not exhibit these characteristics and instead receive research attention due to their importance in fisheries or aquaculture. Black sea bass social structure is unknown, which makes sex change sampling difficult, and few molecular resources are available. The purpose of the present study was to induce sex change using exemestane, an aromatase inhibitor, and assess gonadal gene expression using sex markers (amh, zpc2) and genes involved in steroidogenesis (cyp19a1a, cyp11b), estrogen signaling (esr1, esr2b), and apoptosis or atresia (aen, casp9, fabp11, parg, pdcd4, rif1). Overall, dietary exemestane treatment was effective, and most exposed females exhibited early histological signs of sex change and significantly higher rates of ovarian atresia relative to control females. Genes associated with atresia did not reflect this, however, as expression patterns in sex changing gonads were overall similar to those of ovaries, likely due to a whole ovary dilution effect of the RNA. Still, small but insignificant expression decreases during early sex change were detected for ovary-related genes (aen, casp9, fabp11, zpc2) and anti-apoptotic factors (parg, rif1). Exemestane treatment did not impact spermatogenesis or testicular gene expression, but testes were generally characterized by elevated steroidogenic enzyme and estrogen receptor mRNAs. Further research will be needed to understand these processes in black sea bass, using isolated ovarian follicles and multiple stages of sex change.

Characterization of Genetic and Epigenetic Variation in Sperm and Red Blood Cells from Adult Hatchery and Natural-Origin Steelhead, Oncorhynchus mykiss.

While the goal of most conservation hatchery programs is to produce fish that are genetically and phenotypically indistinguishable from the wild stocks they aim to restore, there is considerable evidence that salmon and steelhead reared in hatcheries differ from wild fish in phenotypic traits related to fitness. Some evidence suggests that these phenotypic differences have a genetic basis (e.g., domestication selection) but another likely mechanism that remains largely unexplored is that differences between hatchery and wild populations arise as a result of environmentally-induced heritable epigenetic change. As a first step toward understanding the potential contribution of these two possible mechanisms, we describe genetic and epigenetic variation in hatchery and natural-origin adult steelhead, Oncorhynchus mykiss, from the Methow River, WA. Our main objectives were to determine if hatchery and natural-origin fish could be distinguished genetically and whether differences in epigenetic programming (DNA methylation) in somatic and germ cells could be detected between the two groups. Genetic analysis of 72 fish using 936 SNPs generated by Restriction Site Associated DNA Sequencing (RAD-Seq) did not reveal differentiation between hatchery and natural-origin fish at a population level. We performed Reduced Representation Bisulfite Sequencing (RRBS) on a subset of 10 hatchery and 10 natural-origin fish and report the first genome-wide characterization of somatic (red blood cells (RBCs)) and germ line (sperm) derived DNA methylomes in a salmonid, from which we identified considerable tissue-specific methylation. We identified 85 differentially methylated regions (DMRs) in RBCs and 108 DMRs in sperm of steelhead reared for their first year in a hatchery environment compared to those reared in the wild. This work provides support that epigenetic mechanisms may serve as a link between hatchery rearing and adult phenotype in steelhead; furthermore, DMRs identified in germ cells (sperm) highlight the potential for these changes to be passed on to future generations.

A teleost androgen promotes development of primary ovarian follicles in coho salmon and rapidly alters the ovarian transcriptome.

Recent studies using several teleost models have revealed that androgens increase the size of previtellogenic (primary and/or early secondary) ovarian follicles. To explore our hypothesis that androgens drive the development of primary follicles into early secondary follicles, and to determine the mechanisms underlying these androgenic effects, we exposed juvenile coho salmon to near-physiological and relatively sustained levels of the nonaromatizable androgen 11-ketotestosterone (11-KT). This resulted in significant growth of primary ovarian follicles after 10 and 20 days, with follicles after 20 days displaying a morphological phenotype characteristic of early secondary follicles (presence of cortical alveoli). Utilizing the same experimental approach, we then analyzed how 11-KT rapidly altered the ovarian transcriptome after 1 and 3 days of treatment. RNA-Seq analysis revealed that 69 (day 1) and 1,022 (day 3) contiguous sequences (contigs) were differentially expressed relative to controls. The differentially expressed contigs mapped to genes including those encoding proteins involved in gonadotropin, steroid hormone, and growth factor signaling, and in cell and ovarian development, including genes with putative androgen-response elements. Biological functions and canonical pathways identified as potentially altered by 11-KT include those involved in ovarian development, tissue differentiation and remodeling, and lipid metabolism. We conclude that androgens play a major role in stimulating primary ovarian follicle development and the transition into secondary growth.

Extending Immunological Profiling in the Gilthead Sea Bream, Sparus aurata, by Enriched cDNA Library Analysis, Microarray Design and Initial Studies upon the Inflammatory Response to PAMPs.

This study describes the development and validation of an enriched oligonucleotide-microarray platform for Sparus aurata (SAQ) to provide a platform for transcriptomic studies in this species. A transcriptome database was constructed by assembly of gilthead sea bream sequences derived from public repositories of mRNA together with reads from a large collection of expressed sequence tags (EST) from two extensive targeted cDNA libraries characterizing mRNA transcripts regulated by both bacterial and viral challenge. The developed microarray was further validated by analysing monocyte/macrophage activation profiles after challenge with two Gram-negative bacterial pathogen-associated molecular patterns (PAMPs; lipopolysaccharide (LPS) and peptidoglycan (PGN)). Of the approximately 10,000 EST sequenced, we obtained a total of 6837 EST longer than 100 nt, with 3778 and 3059 EST obtained from the bacterial-primed and from the viral-primed cDNA libraries, respectively. Functional classification of contigs from the bacterial- and viral-primed cDNA libraries by Gene Ontology (GO) showed that the top five represented categories were equally represented in the two libraries: metabolism (approximately 24% of the total number of contigs), carrier proteins/membrane transport (approximately 15%), effectors/modulators and cell communication (approximately 11%), nucleoside, nucleotide and nucleic acid metabolism (approximately 7.5%) and intracellular transducers/signal transduction (approximately 5%). Transcriptome analyses using this enriched oligonucleotide platform identified differential shifts in the response to PGN and LPS in macrophage-like cells, highlighting responsive gene-cassettes tightly related to PAMP host recognition. As observed in other fish species, PGN is a powerful activator of the inflammatory response in S. aurata macrophage-like cells. We have developed and validated an oligonucleotide microarray (SAQ) that provides a platform enriched for the study of gene expression in S. aurata with an emphasis upon immunity and the immune response.

Novel adverse outcome pathways revealed by chemical genetics in a developing marine fish.

Crude oil spills are a worldwide ocean conservation threat. Fish are particularly vulnerable to the oiling of spawning habitats, and crude oil causes severe abnormalities in embryos and larvae. However, the underlying mechanisms for these developmental defects are not well understood. Here, we explore the transcriptional basis for four discrete crude oil injury phenotypes in the early life stages of the commercially important Atlantic haddock (Melanogrammus aeglefinus). These include defects in (1) cardiac form and function, (2) craniofacial development, (3) ionoregulation and fluid balance, and (4) cholesterol synthesis and homeostasis. Our findings suggest a key role for intracellular calcium cycling and excitation-transcription coupling in the dysregulation of heart and jaw morphogenesis. Moreover, the disruption of ionoregulatory pathways sheds new light on buoyancy control in marine fish embryos. Overall, our chemical-genetic approach identifies initiating events for distinct adverse outcome pathways and novel roles for individual genes in fundamental developmental processes.

Sea lampreys elicit strong transcriptomic responses in the lake trout liver during parasitism.

BACKGROUND: The sea lamprey (Petromyzon marinus) is a jawless vertebrate that parasitizes fish as an adult and, with overfishing, was responsible for the decline in lake trout (Salvelinus namaycush) populations in the Great Lakes. While laboratory studies have looked at the rates of wounding on various fish hosts, there have been few investigations on the physiological effects of lamprey wounding on the host. In the current study, two morphotypes of lake trout, leans and siscowets, were parasitized in the laboratory by sea lampreys and the liver transcriptomes of parasitized and nonparasitized fish were analyzed by RNA-seq (DESeq2 and edgeR) to determine which genes and gene pathways (Ingenuity Pathway Analysis) were altered by lamprey parasitism. RESULTS: Overall, genes encoding molecules involved in catalytic (e.g., enzymatic) and binding activities (factors and regulators) predominated the regulated gene lists. In siscowets, the top upregulated gene was growth arrest and DNA-damage-inducible protein and for leans it was interleukin-18-binding protein. In leans, the most significantly downregulated gene was UDP-glucuronosyltransferase 2A2 - DESeq2 or phosphotriesterase related - edgeR. For siscowets, the top downregulated gene was C-C motif chemokine 19 - DESeq2 or GTP-binding protein Rhes - edgeR. Gene pathways associated with inflammatory-related responses or factors (cytokines, chemokines, oxidative stress, apoptosis) were regulated following parasitism in both morphotypes. However, pathways related to energy metabolism (glycolysis, gluconeogenesis, lipolysis, lipogenesis) were also regulated. These pathways or the intensity or direction (up/downregulation) of regulation were different between leans and siscowets. Finally, one of the most significantly downregulated pathways in both leans and siscowets was the kynurenine (tryptophan degradation) pathway. CONCLUSIONS: The results indicate a strong transcriptional response in the lake trout to lamprey parasitism that entails genes involved in the regulation of inflammation and cellular damage. Responses to energy utilization as well as hydromineral balance also occurred indicating an adjustment in the host to energy demands and osmotic imbalances during parasitism. Given the role of the kynurenine pathway in promoting immunotolerance in mammals, the downregulation observed in this pathway during parasitism may signify an attempt by the host to inhibit any feedback suppression of the immune response to the lamprey.

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 2 x 10(6) 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.

Ontogenetic changes in embryonic and brain gene expression in progeny produced from migratory and resident Oncorhynchus mykiss.

Little information has been gathered regarding the ontogenetic changes that contribute to differentiation between resident and migrant individuals, particularly before the onset of gross morphological and physiological changes in migratory individuals. The aim of this study was to evaluate gene expression during early development in Oncorhynchus mykiss populations with different life histories, in a tissue known to integrate environmental cues to regulate complex developmental processes and behaviours. We sampled offspring produced from migrant and resident parents, collecting whole embryos prior to the beginning of first feeding, and brain tissue at three additional time points over the first year of development. RNA sequencing for 32 individuals generated a reference transcriptome of 30 177 genes that passed count thresholds. Differential gene expression between migrant and resident offspring was observed for 1982 genes. The greatest number of differentially expressed genes occurred at 8 months of age, in the spring a full year before the obvious physiological transformation from stream-dwelling parr to sea water-adaptable smolts begins for migrant individuals. Sex and age exhibited considerable effects on differential gene expression between migrants and resident offspring. Differential gene expression was observed in genes previously associated with migration, but also in genes previously unassociated with early life history divergence. Pathway analysis revealed coordinated differential expression in genes related to phototransduction, which could modulate photoperiod responsiveness and variation in circadian rhythms. The role for early differentiation in light sensitivity and biological rhythms is particularly intriguing in understanding early brain processes involved in differentiation of migratory and resident life history types.

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.

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.

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.

Differential gene expression associated with dietary methylmercury (MeHg) exposure in rainbow trout (Oncorhynchus mykiss) and zebrafish (Danio rerio).

The objective of this study was to identify and evaluate conserved biomarkers that could be used in most species of teleost fish at most life-stages. We investigated the effects of sublethal methylmercury (MeHg) exposure on developing rainbow trout and zebrafish. Juvenile rainbow trout and young adult zebrafish were fed food with MeHg added at 0, 0.5, 5, and 50 ppm. Atomic absorption spectrometry was applied to measure whole body total Hg levels, and pathologic analysis was performed to identify MeHg-induced toxicity. Fish at 6 weeks were sampled from each group for microarray analysis using RNA from whole fish. MeHg-exposed trout and zebrafish did not show overt signs of toxicity or pathology, nor were significant differences seen in mortality, length, mass, or condition factor. The accumulation of MeHg in trout and zebrafish exhibited dose- and time-dependent patterns during 6 weeks, and zebrafish exhibited greater assimilation of total Hg than rainbow trout. The dysregulated genes in MeHg-treated fish have multiple functional annotations, such as iron ion homeostasis, glutathione transferase activity, regulation of muscle contraction, troponin I binding and calcium-dependent protein binding. Genes were selected as biomarker candidates based on their microarray data and their expression was evaluated by QPCR. Unfortunately, these genes are not good consistent biomarkers for both rainbow trout and zebrafish from QPCR evaluation using individual fish. Our conclusion is that biomarker analysis for aquatic toxicant assessment using fish needs to be based on tissue-, sex- and species-specific consideration.

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.

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.

Stimulation of growth and changes in the hepatic transcriptome by 17beta-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.

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.