High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis.

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Fixation of allelic gene expression landscapes and expression bias pattern shape the transcriptome of the clonal Amazon molly.

The Amazon molly is a unique clonal fish species that originated from an interspecies hybrid between Poecilia species P. mexicana and P. latipinna It reproduces by gynogenesis, which eliminates paternal genomic contribution to offspring. An earlier study showed that Amazon molly shows biallelic expression for a large portion of the genome, leading to two main questions: (1) Are the allelic expression patterns from the initial hybridization event stabilized or changed during establishment of the asexual species and its further evolution? (2) Is allelic expression biased toward one parental allele a stochastic or adaptive process? To answer these questions, the allelic expression of P. formosa siblings was assessed to investigate intra- and inter-cohort allelic expression variability. For comparison, interspecies hybrids between P. mexicana and P. latipinna were produced in the laboratory to represent the P. formosa ancestor. We have identified inter-cohort and intra-cohort variation in parental allelic expression. The existence of inter-cohort divergence suggests functional P. formosa allelic expression patterns do not simply reflect the atavistic situation of the first interspecies hybrid but potentially result from long-term selection of transcriptional fitness. In addition, clonal fish show a transcriptional trend representing minimal intra-clonal variability in allelic expression patterns compared to the corresponding hybrids. The intra-clonal similarity in gene expression translates to sophisticated genetic functional regulation at the individuum level. These findings suggest the parental alleles inherited by P. formosa form tightly regulated genetic networks that lead to a stable transcriptomic landscape within clonal individuals.

Particulate hexavalent chromium alters microRNAs in human lung cells that target key carcinogenic pathways.

Hexavalent chromium [Cr(VI)] is a global environmental pollutant and human lung carcinogen. However, the mechanisms of Cr(VI) carcinogenesis are not well defined. Cr(VI)-altered gene expression has been reported in the literature and is implicated in numerous mechanisms of Cr(VI) carcinogenesis. MicroRNAs (miRNAs) play a key role in controlling gene expression and are associated with carcinogenic mechanisms. To date no studies have evaluated global changes in miRNA expression in human cells after Cr(VI) exposure. We used RNA sequencing to evaluate how a particulate Cr(VI) compound (zinc chromate), the most potent form of Cr(VI), alters global miRNA expression after acute (24 h) or prolonged (72 and 120 h) exposure to 0.1, 0.2 and 0.3 µg/cm2 zinc chromate in an immortalized, non-cancerous human lung cell line (WTHBF-6). Particulate Cr(VI) significantly affected expression of miRNAs at all time points and concentrations tested. We also found the number of significantly downregulated miRNAs increased in a time- and concentration-dependent manner and many miRNAs were upregulated after 24 h exposure at the intermediate concentration tested. Pathway analyses of the differentially expressed miRNAs predicted miRNAs target pathways of Cr(VI) carcinogenesis in a time- and concentration-dependent manner. These data are the first to evaluate global changes in miRNA expression in human lung cells after Cr(VI) exposure and indicate miRNAs may play a key role in pathways of Cr(VI) carcinogenesis.

Global assessment of organ specific basal gene expression over a diurnal cycle with analyses of gene copies exhibiting cyclic expression patterns.

BACKGROUND: Studying functional divergences between paralogs that originated from genome duplication is a significant topic in investigating molecular evolution. Genes that exhibit basal level cyclic expression patterns including circadian and light responsive genes are important physiological regulators. Temporal shifts in basal gene expression patterns are important factors to be considered when studying genetic functions. However, adequate efforts have not been applied to studying basal gene expression variation on a global scale to establish transcriptional activity baselines for each organ. Furthermore, the investigation of cyclic expression pattern comparisons between genome duplication created paralogs, and potential functional divergence between them has been neglected. To address these questions, we utilized a teleost fish species, Xiphophorus maculatus, and profiled gene expression within 9 organs at 3-h intervals throughout a 24-h diurnal period. RESULTS: Our results showed 1.3-21.9% of genes in different organs exhibited cyclic expression patterns, with eye showing the highest fraction of cycling genes while gonads yielded the lowest. A majority of the duplicated gene pairs exhibited divergences in their basal level expression patterns wherein only one paralog exhibited an oscillating expression pattern, or both paralogs exhibit oscillating expression patterns, but each gene duplicate showed a different peak expression time, and/or in different organs. CONCLUSIONS: These observations suggest cyclic genes experienced significant sub-, neo-, or non-functionalization following the teleost genome duplication event. In addition, we developed a customized, web-accessible, gene expression browser to facilitate data mining and data visualization for the scientific community.

Workshop report: Cryopreservation of aquatic biomedical models.

The genetic resources of aquatic biomedical model organisms are the products of millions of years of evolution, decades of scientific development, and hundreds of millions of dollars of research funding investment. Genetic resources (e.g., specific alleles, transgenes, or combinations) of each model organism can be considered a form of scientific wealth that can be accumulated and exchanged, typically in the form of live animals or germplasm. Large-scale maintenance of live aquatic organisms that carry these genetic resources is inefficient, costly, and risky. In situ maintenance may be substantially enhanced and backed up by combining cryopreserved germplasm repositories and genetic information systems with live animal culture. Unfortunately, cryopreservation has not advanced much beyond the status of an exploratory research for most aquatic species, lacks widespread application, and methods for successful cryopreservation remain poorly defined. For most aquatic species biological materials other than sperm or somatic cells are not comprehensively banked to represent and preserve a broad range of genetic diversity for each species. Therefore, new approaches and standardization are needed for repository-level application to ensure reproducible recovery of cryopreserved materials. Additionally, development of new technologies is needed to address preservation of novel biological materials, such as eggs and embryos of aquatic species. To address these goals, the Office of Research Infrastructure Programs (ORIP) of the National Institutes of Health (NIH) hosted the Cryopreservation of Aquatic Biomedical Models Workshop on January 7 to 8, 2017, in conjunction with the 8th Aquatic Animal Models of Human Disease Conference in Birmingham, Alabama. The goals of the workshop were to assess the status of germplasm cryopreservation in various biomedical aquatic models and allow representatives of the scientific community to develop and prioritize a consensus of specific actionable recommendations that will move the field of cryopreservation of aquatic resources forward. This workshop included sessions devoted to new approaches for cryopreservation of aquatic species, discussion of current efforts and approaches in preservation of aquatic model germplasm, consideration of needs for standardization of methods to support reproducibility, and enhancement of repository development by establishment of scalable high-throughput technologies. The following three broad recommendations were forwarded from workshop attendees: 1: Establish a comprehensive, centralized unit ("hub") to programmatically develop training for and documentation of cryopreservation methods for aquatic model systems. This would include development of species-specific protocols and approaches, outreach programs, community development and standardization, freezing services and training of the next generation of experts in aquatic cryopreservation. 2: Provide mechanisms to support innovative technical advancements that will increase the reliability, reproducibility, simplicity, throughput, and efficiency of the cryopreservation process, including vitrification and pipelines for sperm, oocytes, eggs, embryos, larvae, stem cells, and somatic cells of all aquatic species. This recommendation encompasses basic cryopreservation knowledge and engineering technology, such as microfluidics and automated processing technologies. 3: Implement mechanisms that allow the various aquatic model stock centers to increase their planning, personnel, ability to secure genetic resources and to promote interaction within an integrated, comprehensive repository network for aquatic model species repositories.

On-Site Capabilities of a Mobile Laboratory for Aquatic Germplasm Cryopreservation.

Cryopreservation of genetic material can become an important tool for user groups in imperiled fishes, wild fisheries, aquaculture, and biomedical research. Persistent challenges within aquatic species cryopreservation are standardization and reliable collection of diverse, high quality samples. The overall goal of this study was to work with different user groups and cryopreserve sperm on-site at their facilities to evaluate the uses and challenges of a mobile laboratory with high-throughput and quality control capabilities comparable to those of a specialized centralized facility. The objectives were to demonstrate collection and cryopreservation of sperm of: 1) large-bodied freshwater Blue Catfish (Ictalurus furcatus) for aquaculture; 2) small-bodied freshwater Xiphophorus for biomedical and imperiled repository development, and 3) saltwater Red Snapper (Lutjanus campechanus) for wild fisheries research. Over the course of this project, the mobile laboratory traveled more than 4,000 km collecting germplasm from more than 650 male fishes. A total of 136 Blue Catfish were processed in 2015 and 2016 resulting in a total of 6,146 0.5-mL French straws. A total of 521 males from 11 different species in the genus Xiphophorus were processed over 4 d in 2015 resulting in a total of 488 0.25-mL French straws. And, a total of 17 Red Snapper males were processed during 2015 resulting in a total of 316 0.5-mL French straws. This is the first development of a mobile laboratory with high-throughput capability for aquatic species. User groups would no longer be limited to germplasm resources that can only be shipped as samples or transported as live animals to a central cryopreservation facility. Mobile laboratories create opportunities to collect higher quality germplasm, provide access to new species, and enable direct cooperation, including training, with a wide variety of user groups and applications.

Waveband specific transcriptional control of select genetic pathways in vertebrate skin (Xiphophorus maculatus).

BACKGROUND: Evolution occurred exclusively under the full spectrum of sunlight. Conscription of narrow regions of the solar spectrum by specific photoreceptors suggests a common strategy for regulation of genetic pathways. Fluorescent light (FL) does not possess the complexity of the solar spectrum and has only been in service for about 60 years. If vertebrates evolved specific genetic responses regulated by light wavelengths representing the entire solar spectrum, there may be genetic consequences to reducing the spectral complexity of light. RESULTS: We utilized RNA-Seq to assess changes in the transcriptional profiles of Xiphophorus maculatus skin after exposure to FL ("cool white"), or narrow wavelength regions of light between 350 and 600 nm (i.e., 50 nm or 10 nm regions, herein termed "wavebands"). Exposure to each 50 nm waveband identified sets of genes representing discrete pathways that showed waveband specific transcriptional modulation. For example, 350-400 or 450-500 nm waveband exposures resulted in opposite regulation of gene sets marking necrosis and apoptosis (i.e., 350-400 nm; necrosis suppression, apoptosis activation, while 450-500 nm; apoptosis suppression, necrosis activation). Further investigation of specific transcriptional modulation employing successive 10 nm waveband exposures between 500 and 550 nm showed; (a) greater numbers of genes may be transcriptionally modulated after 10 nm exposures, than observed for 50 nm or FL exposures, (b) the 10 nm wavebands induced gene sets showing greater functional specificity than 50 nm or FL exposures, and (c) the genetic effects of FL are primarily due to 30 nm between 500 and 530 nm. Interestingly, many genetic pathways exhibited completely opposite transcriptional effects after different waveband exposures. For example, the epidermal growth factor (EGF) pathway exhibits transcriptional suppression after FL exposure, becomes highly active after 450-500 nm waveband exposure, and again, exhibits strong transcriptional suppression after exposure to the 520-530 nm waveband. CONCLUSIONS: Collectively, these results suggest one may manipulate transcription of specific genetic pathways in skin by exposure of the intact animal to specific wavebands of light. In addition, we identify genes transcriptionally modulated in a predictable manner by specific waveband exposures. Such genes, and their regulatory elements, may represent valuable tools for genetic engineering and gene therapy protocols.

Germ cell and tumor associated piRNAs in the medaka and Xiphophorus melanoma models.

BACKGROUND: A growing number of studies report an abnormal expression of Piwi-interacting RNAs (piRNAs) and the piRNA processing enzyme Piwi in many cancers. Whether this finding is an epiphenomenon of the chaotic molecular biology of the fast dividing, neoplastically transformed cells or is functionally relevant to tumorigenesisis is difficult to discern at present. To better understand the role of piRNAs in cancer development small laboratory fish models can make a valuable contribution. However, little is known about piRNAs in somatic and neoplastic tissues of fish. RESULTS: To identify piRNA clusters that might be involved in melanoma pathogenesis, we use several transgenic lines of medaka, and platyfish/swordtail hybrids, which develop various types of melanoma. In these tumors Piwi, is expressed at different levels, depending on tumor type. To quantify piRNA levels, whole piRNA populations of testes and melanomas of different histotypes were sequenced. Because no reference piRNA cluster set for medaka or Xiphophorus was yet available we developed a software pipeline to detect piRNA clusters in our samples and clusters were selected that were enriched in one or more samples. We found several loci to be overexpressed or down-regulated in different melanoma subtypes as compared to hyperpigmented skin. Furthermore, cluster analysis revealed a clear distinction between testes, low-grade and high-grade malignant melanoma in medaka. CONCLUSIONS: Our data imply that dysregulation of piRNA expression may be associated with development of melanoma. Our results also reinforce the importance of fish as a suitable model system to study the role of piRNAs in tumorigenesis.

X. couchianus and X. hellerii genome models provide genomic variation insight among Xiphophorus species.

BACKGROUND: Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species. RESULTS: We present sequencing, assembly and annotation of two new genomes representing Xiphophorus couchianus and Xiphophorus hellerii. The final X. couchianus and X. hellerii assemblies have total sizes of 708 Mb and 734 Mb and correspond to 98 % and 102 % of the X. maculatus Jp 163 A genome size, respectively. The rates of single nucleotide change range from 1 per 52 bp to 1 per 69 bp among the three genomes and the impact of putatively damaging variants are presented. In addition, a survey of transposable elements allowed us to deduce an ancestral TE landscape, uncovered potential active TEs and document a recent burst of TEs during evolution of this genus. CONCLUSIONS: Two new Xiphophorus genomes and their corresponding transcriptomes were efficiently assembled, the former using a novel guided assembly approach. Three assembled genome sequences within this single vertebrate order of new world live-bearing fishes will accelerate our understanding of relationship between environmental adaptation and genome evolution. In addition, these genome resources provide capability to determine allele specific gene regulation among interspecies hybrids produced by crossing any of the three species that are known to produce progeny predisposed to tumor development.

Xiphophorus and Medaka Cancer Models.

Besides recently developed zebrafish cancer models, other fish species have been employed for many years as cancer models in laboratory studies. Two models, namely in Xiphophorus and medaka have proven useful in providing important clues to cancer etiology. Medaka is a complementary model to zebrafish in many areas of research since it offers similar resources and experimental tools. Xiphophorus provides the advantages of a natural ("evolutionary mutant") model with established genetics. Xiphophorus hybrids can develop spontaneous and radiation or carcinogen induced cancers. This chapter describes the tumor models in both species, which mainly focus on melanoma, and summarizes the main findings and future research directions.

Genomic and physiological footprint of the Deepwater Horizon oil spill on resident marsh fishes.

The biological consequences of the Deepwater Horizon oil spill are unknown, especially for resident organisms. Here, we report results from a field study tracking the effects of contaminating oil across space and time in resident killifish during the first 4 mo of the spill event. Remote sensing and analytical chemistry identified exposures, which were linked to effects in fish characterized by genome expression and associated gill immunohistochemistry, despite very low concentrations of hydrocarbons remaining in water and tissues. Divergence in genome expression coincides with contaminating oil and is consistent with genome responses that are predictive of exposure to hydrocarbon-like chemicals and indicative of physiological and reproductive impairment. Oil-contaminated waters are also associated with aberrant protein expression in gill tissues of larval and adult fish. These data suggest that heavily weathered crude oil from the spill imparts significant biological impacts in sensitive Louisiana marshes, some of which remain for over 2 mo following initial exposures.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation.

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

Comprehensive phylogenetic analysis of all species of swordtails and platies (Pisces: Genus Xiphophorus) uncovers a hybrid origin of a swordtail fish, Xiphophorus monticolus, and demonstrates that the sexually selected sword originated in the ancestral lineage of the genus, but was lost again secondarily.

BACKGROUND: Males in some species of the genus Xiphophorus, small freshwater fishes from Meso-America, have an extended caudal fin, or sword - hence their common name "swordtails". Longer swords are preferred by females from both sworded and - surprisingly also, non-sworded (platyfish) species that belong to the same genus. Swordtails have been studied widely as models in research on sexual selection. Specifically, the pre-existing bias hypothesis was interpreted to best explain the observed bias of females in presumed ancestral lineages of swordless species that show a preference for assumed derived males with swords over their conspecific swordless males. However, many of the phylogenetic relationships within this genus still remained unresolved. Here we construct a comprehensive molecular phylogeny of all 26 known Xiphophorus species, including the four recently described species (X. kallmani, X. mayae, X. mixei and X. monticolus). We use two mitochondrial and six new nuclear markers in an effort to increase the understanding of the evolutionary relationships among the species in this genus. Based on the phylogeny, the evolutionary history and character state evolution of the sword was reconstructed and found to have originated in the common ancestral lineage of the genus Xiphophorus and that it was lost again secondarily. RESULTS: We estimated the evolutionary relationships among all known species of the genus Xiphophorus based on the largest set of DNA markers so far. The phylogeny indicates that one of the newly described swordtail species, Xiphophorus monticolus, is likely to have arisen through hybridization since it is placed with the southern platyfish in the mitochondrial phylogeny, but with the southern swordtails in the nuclear phylogeny. Such discordance between these two types of markers is a strong indication for a hybrid origin. Additionally, by using a maximum likelihood approach the possession of the sexually selected sword trait is shown to be the most likely ancestral state for the genus Xiphophorus. Further, we provide a well supported estimation of the phylogenetic relationships between the previously unresolved northern swordtail groups. CONCLUSIONS: This comprehensive molecular phylogeny of the entire genus Xiphophorus provides evidence that a second swordtail species, X. monticolus, arose through hybridization. Previously, we demonstrated that X. clemenciae, another southern swordtail species, arose via hybridization. These findings highlight the potential key role of hybridization in the evolution of this genus and suggest the need for further investigations into how hybridization contributes to speciation more generally.

RNA-Seq reveals complex genetic response to Deepwater Horizon oil release in Fundulus grandis.

BACKGROUND: The release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future. RESULTS: Our de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set. CONCLUSION: RNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.

Lineage-specific co-evolution of the Egf receptor/ligand signaling system.

BACKGROUND: The epidermal growth factor receptor (Egfr) with its numerous ligands has fundamental roles in development, cell differentiation and physiology. Dysfunction of the receptor-ligand system contributes to many human malignancies. Consistent with such various tasks, the Egfr gene family has expanded during vertebrate evolution as a consequence of several rounds of whole genome duplication. Of particular interest is the effect of the fish-specific whole genome duplication (FSGD) on the ligand-receptor system, as it has supplied this largest group of vertebrates with additional opportunities for sub- and/or neofunctionalization in this signaling system. RESULTS: We identified the predicted components of the Egf receptor-ligand signaling system in teleost fishes (medaka, platyfish, stickleback, pufferfishes and zebrafish). We found two duplicated egfr genes, egfra and egfrb, in all available teleost genomes. Surprisingly only one copy for each of the seven Egfr ligands could be identified in most fishes, with zebrafish hbegf being the only exception. Special focus was put on medaka, for which we more closely investigated all Egf receptors and Egfr ligands. The different expression patterns of egfra, egfrb and their ligands in medaka tissues and embryo stages suggest differences in role and function. Preferential co-expression of different subsets of Egfr ligands corroborates the possible subfunctionalization and specialization of the two receptors in adult tissues. Bioinformatic analyses of the ligand-receptor interface between Egfr and its ligands show a very weak evolutionary conservation within this region. Using in vitro analyses of medaka Egfra, we could show that this receptor is only activated by medaka ligands, but not by human EGF. Altogether, our data suggest a lineage-specific Egfr/Egfr ligand co-evolution. CONCLUSIONS: Our data indicate that medaka Egfr signaling occurs via its two copies, Egfra and Egfrb, each of them being preferentially coexpressed with different subsets of Egfr ligands. This fish-specific occurrence of Egf receptor specialization offers unique opportunities to study the functions of different Egf receptor-ligand combinations and their biological outputs in vertebrates. Furthermore, our results strongly support the use of homologous ligands in future studies, as sufficient cross-specificity is very unlikely for this ligand/receptor system.

Deconvoluting Wavelengths Leading to Fluorescent Light Induced Inflammation and Cellular Stress in Zebrafish (Danio rerio).

Fluorescent light (FL) has been shown to induce a cellular immune and inflammatory response that is conserved over 450 MY of evolutionary divergence and among vertebrates having drastically different lifestyles such as Mus musculus, Danio rerio, Oryzias latipes and Xiphophorus maculatus. This surprising finding of an inflammation and immune response to FL not only holds for direct light receiving organs (skin) but is also observed within internal organs (brain and liver). Light responsive genetic circuitry initiated by the IL1B regulator induces a highly conserved acute phase response in each organ assessed for all of biological models surveyed to date; however, the specific light wavelengths triggering this response have yet to be determined so investigation of mechanisms and/or light specific molecule(s) leading to this response are difficult to assess. To understand how specific light wavelengths are received in both external and internal organs, zebrafish were exposed to specific 50 nm light wavebands spanning the visible spectrum from 300-600 nm and the genetic responses to each waveband exposure were assessed. Surprisingly, the induced cellular stress response previously observed following FL exposure is not triggered by the lower "damaging" wavelengths of light (UVB and UVA from 300-400 nm) but instead is maximally induced by higher wavelengths ranging from 450-500 nm in skin to 500-600 nm in both brain and liver).

Intra-Strain Genetic Variation of Platyfish (Xiphophorus maculatus) Strains Determines Tumorigenic Trajectory.

Xiphophorus interspecies hybrids represent a valuable model system to study heritable tumorigenesis, and the only model system that exhibits both spontaneous and inducible tumors. Types of tumorigenesis depend on the specific pedigree of the parental species, X. maculatus, utilized to produce interspecies hybrids. Although the ancestors of the two currently used X. maculatus parental lines, Jp163 A and Jp163 B, were originally siblings produced by the same mother, backcross interspecies hybrid progeny between X. hellerii and X. maculatus Jp163 A develop spontaneous melanoma initiating at the dorsal fin due to segregation of an oncogene and a regulator encoded by the X. maculatus genome, while the backcross hybrid progeny with X. hellerii or X. couchianus and Jp163 B exhibit melanoma on the flanks of their bodies, especially after treatment with ultraviolet light. Therefore, dissecting the genetic differences between these two closely related lines may lead to better understanding of functional molecular differences associated with tumorigenic mechanisms. For this purpose, comparative genomic analyses were undertaken to establish genetic variants between these two X. maculatus lines. Surprisingly, given the heritage of these two fish lines, we found genetic variants are clustered together in select chromosomal regions. Among these variants are non-synonymous mutations located in 381 genes. The non-random distribution of genetic variants between these two may highlight ancestral chromosomal recombination patterns that became fixed during subsequent inbreeding. Employing comparative transcriptomics, we also determined differences in the skin transcriptional landscape between the two lines. The genetic differences observed are associated with pathways highlighting fundamental cellular functions including inter-cellular and microenvironment-cellular interactions, and DNA repair. These results collectively lead to the conclusion that diverged functional genetic baselines are present between Jp163 A and B strains. Further, disruption of these fixed genetic baselines in the hybrids may give rise to spontaneous or inducible mechanisms of tumorigenesis.

Genomic approaches in the identification of hypoxia biomarkers in model fish species.

Eutrophication leading to hypoxic water conditions has become a major problem in aquatic systems worldwide. Monitoring the levels and biological effects of lowered oxygen levels in aquatic systems may provide data useful in management of natural aquatic environments. Fishes represent an economically important resource that is subject to hypoxia exposure effects. Due to the extreme diversity of fish species and their habitats, fishes in general have evolved unique capabilities to modulate gene expression patterns in response to hypoxic stress. Recent studies have attempted to document quantitative changes in gene expression patterns induced in various fish species in response to reduced dissolved oxygen levels. From a management perspective, the goal of these studies is to provide a more complete characterization of hypoxia responsive genes in fish, as molecular indicators (biomarkers) of ecosystem hypoxic stress.The molecular genetic response to hypoxia is highly complex and overlaps with other stress responses making it difficult to identify hypoxia specific responses using traditional single gene or low throughput approaches. Therefore, recent approaches have been aimed at developing functional genomic (e.g. high density microarray and real-time PCR) and proteomic (two-dimensional fluorescence difference in gel electrophoresis coupled with mass spectrometry based peptide identification) technologies that employ fish species. Many of the fish species utilized in these studies do not have the advantages of underlying genome resources (i.e., genome or transcriptome sequences). Efforts have attempted to establish correlations between discreet molecular responses elicited by fish in response to hypoxia and changes in the genetic profiles of stressed organs or tissues. Notable progress in these areas has been made using several different versions of either cDNA or oligonucleotide based microarrays to profile changes in gene expression patterns in response to hypoxic stress.Due to these efforts, hundreds of hypoxia responsive genes have been identified both from laboratory reared aquaria fish and from feral fish derived from both fresh and saltwater habitats. Herein, we review these reports and the emergence of hypoxia biomarker development in aquatic species. We also include some of our own recent results using the medaka (Oryzias latipes) as a model to define genetic profiles of hypoxia exposure.

Expression Signatures of Cisplatin- and Trametinib-Treated Early-Stage Medaka Melanomas.

Small aquarium fish models provide useful systems not only for a better understanding of the molecular basis of many human diseases, but also for first-line screening to identify new drug candidates. For testing new chemical substances, current strategies mostly rely on easy to perform and efficient embryonic screens. Cancer, however, is a disease that develops mainly during juvenile and adult stage. Long-term treatment and the challenge to monitor changes in tumor phenotype make testing of large chemical libraries in juvenile and adult animals cost prohibitive. We hypothesized that changes in the gene expression profile should occur early during anti-tumor treatment, and the disease-associated transcriptional change should provide a reliable readout that can be utilized to evaluate drug-induced effects. For the current study, we used a previously established medaka melanoma model. As proof of principle, we showed that exposure of melanoma developing fish to the drugs cisplatin or trametinib, known cancer therapies, for a period of seven days is sufficient to detect treatment-induced changes in gene expression. By examining whole body transcriptome responses we provide a novel route toward gene panels that recapitulate anti-tumor outcomes thus allowing a screening of thousands of drugs using a whole-body vertebrate model. Our results suggest that using disease-associated transcriptional change to screen therapeutic molecules in small fish model is viable and may be applied to pre-clinical research and development stages in new drug discovery.

Fluorescent Light Incites a Conserved Immune and Inflammatory Genetic Response within Vertebrate Organs (Danio Rerio, Oryzias Latipes and Mus Musculus).

Fluorescent light (FL) has been utilized for ≈60 years and has become a common artificial light source under which animals, including humans, spend increasing amounts of time. Although the solar spectrum is quite dissimilar in both wavelengths and intensities, the genetic consequences of FL exposure have not been investigated. Herein, we present comparative RNA-Seq results that establish expression patterns within skin, brain, and liver for Danio rerio, Oryzias latipes, and the hairless mouse (Mus musculus) after exposure to FL. These animals represent diurnal and nocturnal lifestyles, and ≈450 million years of evolutionary divergence. In all three organisms, FL induced transcriptional changes of the acute phase response signaling pathway and modulated inflammation and innate immune responses. Our pathway and gene clustering analyses suggest cellular perception of oxidative stress is promoting induction of primary up-stream regulators IL1B and TNF. The skin and brain of the three animals as well as the liver of both fish models all exhibit increased inflammation and immune responses; however, the mouse liver suppressed the same pathways. Overall, the conserved nature of the genetic responses observed after FL exposure, among fishes and a mammal, suggest the presence of light responsive genetic circuitry deeply embedded in the vertebrate genome.