Transcriptional background effects on a tumor driver gene in different pigment cell types of medaka.

The Xiphophorus melanoma receptor kinase gene, xmrk, is a bona fide oncogene driving melanocyte tumorigenesis of Xiphophorus fish. When ectopically expressed in medaka, it not only induces development of several pigment cell tumor types in different strains of medaka but also induces different tumor types within the same animal, suggesting its oncogenic activity has a transcriptomic background effect. Although the central pathways that xmrk utilizes to lead to melanomagenesis are well documented, genes and genetic pathways that modulate the oncogenic effect and alter the course of disease have not been studied so far. To understand how the genetic networks between different histocytes of xmrk-driven tumors are composed, we isolated two types of tumors, melanoma and xanthoerythrophoroma, from the same xmrk transgenic medaka individuals, established the transcriptional profiles of both xmrk-driven tumors, and compared (1) genes that are co-expressed with xmrk in both tumor types, and (2) differentially expressed genes and their associated molecular functions, between the two tumor types. Transcriptomic comparisons between the two tumor types show melanoma and xanthoerythrophoroma are characterized by transcriptional features representing varied functions, indicating distinct molecular interactions between the driving oncogene and the cell-type-specific transcriptomes. Melanoma tumors exhibit gene signatures that are relevant to proliferation and invasion, while xanthoerythrophoroma tumors are characterized by expression profiles related to metabolism and DNA repair. We conclude the transcriptomic backgrounds, exemplified by cell-type-specific genes that are downstream of xmrk effected signaling pathways, contribute the potential to change the course of tumor development and may affect overall tumor outcomes.

Oncogenic allelic interaction in Xiphophorus highlights hybrid incompatibility.

Mixing genomes of different species by hybridization can disrupt species-specific genetic interactions that were adapted and fixed within each species population. Such disruption can predispose the hybrids to abnormalities and disease that decrease the overall fitness of the hybrids and is therefore named as hybrid incompatibility. Interspecies hybridization between southern platyfish and green swordtails leads to lethal melanocyte tumorigenesis. This occurs in hybrids with tumor incidence following progeny ratio that is consistent with two-locus interaction, suggesting melanoma development is a result of negative epistasis. Such observations make Xiphophorus one of the only two vertebrate hybrid incompatibility examples in which interacting genes have been identified. One of the two interacting loci has been characterized as a mutant epidermal growth factor receptor. However, the other locus has not been identified despite over five decades of active research. Here we report the localization of the melanoma regulatory locus to a single gene, rab3d, which shows all expected features of the long-sought oncogene interacting locus. Our findings provide insights into the role of egfr regulation in regard to cancer etiology. Finally, they provide a molecular explainable example of hybrid incompatibility.

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.

Application of the Transcriptional Disease Signature (TDSs) to Screen Melanoma-Effective Compounds in a Small Fish Model.

Cell culture and protein target-based compound screening strategies, though broadly utilized in selecting candidate compounds, often fail to eliminate candidate compounds with non-target effects and/or safety concerns until late in the drug developmental process. Phenotype screening using intact research animals is attractive because it can help identify small molecule candidate compounds that have a high probability of proceeding to clinical use. Most FDA approved, first-in-class small molecules were identified from phenotypic screening. However, phenotypic screening using rodent models is labor intensive, low-throughput, and very expensive. As a novel alternative for small molecule screening, we have been developing gene expression disease profiles, termed the Transcriptional Disease Signature (TDS), as readout of small molecule screens for therapeutic molecules. In this concept, compounds that can reverse, or otherwise affect known disease-associated gene expression patterns in whole animals may be rapidly identified for more detailed downstream direct testing of their efficacy and mode of action. To establish proof of concept for this screening strategy, we employed a transgenic strain of a small aquarium fish, medaka (Oryzias latipes), that overexpresses the malignant melanoma driver gene xmrk, a mutant egfr gene, that is driven by a pigment cell-specific mitf promoter. In this model, melanoma develops with 100% penetrance. Using the transgenic medaka malignant melanoma model, we established a screening system that employs the NanoString nCounter platform to quantify gene expression within custom sets of TDS gene targets that we had previously shown to exhibit differential transcription among xmrk-transgenic and wild-type medaka. Compound-modulated gene expression was identified using an internet-accessible custom-built data processing pipeline. The effect of a given drug on the entire TDS profile was estimated by comparing compound-modulated genes in the TDS using an activation Z-score and Kolmogorov-Smirnov statistics. TDS gene probes were designed that target common signaling pathways that include proliferation, development, toxicity, immune function, metabolism and detoxification. These pathways may be utilized to evaluate candidate compounds for potential favorable, or unfavorable, effects on melanoma-associated gene expression. Here we present the logistics of using medaka to screen compounds, as well as, the development of a user-friendly NanoString data analysis pipeline to support feasibility of this novel TDS drug-screening strategy.

Comparison of Xiphophorus and human melanoma transcriptomes reveals conserved pathway interactions.

Comparative analysis of human and animal model melanomas can uncover conserved pathways and genetic changes that are relevant for the biology of cancer cells. Spontaneous melanoma in Xiphophorus interspecies backcross hybrid progeny may be informative in identifying genes and functional pathways that are similarly related to melanoma development in all vertebrates, including humans. To assess functional pathways involved in the Xiphophorus melanoma, we performed gene expression profiling of the melanomas produced in interspecies BC1 and successive backcross generations (i.e., BC5 ) of the cross: X. hellerii × [X. maculatus Jp 163 A × X. hellerii]. Using RNA-Seq, we identified genes that are transcriptionally co-expressed with the driver oncogene, xmrk. We determined functional pathways in the fish melanoma that are also present in human melanoma cohorts that may be related to dedifferentiation based on the expression levels of pigmentation genes. Shared pathways between human and Xiphophorus melanomas are related to inflammation, cell migration, cell proliferation, pigmentation, cancer development, and metastasis. Our results suggest xmrk co-expressed genes are associated with dedifferentiation and highlight these signaling pathways as playing important roles in melanomagenesis.

A three year study of metal levels in skin biopsies of whales in the Gulf of Mexico after the Deepwater Horizon oil crisis.

In response to the explosion of the Deepwater Horizon and the massive release of oil that followed, we conducted three annual research voyages to investigate how the oil spill would impact the marine offshore environment. Most investigations into the ecological and toxicological impacts of the Deepwater Horizon Oil crisis have mainly focused on the fate of the oil and dispersants, but few have considered the release of metals into the environment. From studies of previous oil spills, other marine oil industries, and analyses of oil compositions, it is evident that metals are frequently encountered. Several metals have been reported in the MC252 oil from the Deepwater Horizon oil spill, including the nonessential metals aluminum, arsenic, chromium, nickel, and lead; genotoxic metals, such as these are able to damage DNA and can bioaccumulate in organisms resulting in persistent exposure. In the Gulf of Mexico, whales are the apex species; hence we collected skin biopsies from sperm whales (Physeter macrocephalus), short-finned pilot whales (Globicephala macrorhynchus), and Bryde's whales (Balaenoptera edeni). The results from our three-year study of monitoring metal levels in whale skin show (1) genotoxic metals at concentrations higher than global averages previously reported and (2) patterns for MC252-relevant metal concentrations decreasing with time from the oil spill.

Molecular genetic analysis of the melanoma regulatory locus in Xiphophorus interspecies hybrids.

Development of spontaneous melanoma in Xiphophorus interspecies backcross hybrid progeny, (X. hellerii × [X. maculatus Jp 163 A × X. hellerii]) is due to Mendelian segregation of a oncogene (xmrk) and a molecularly uncharacterized locus, called R(Diff), on LG5. R(Diff) is thought to suppresses the activity of xmrk in healthy X. maculatus Jp 163 A parental species that rarely develop melanoma. To better understand the molecular genetics of R(Diff), we utilized RNA-Seq to study allele-specific gene expression of spontaneous melanoma tumors and corresponding normal skin samples derived from 15 first generation backcross (BC1 ) hybrids and 13 fifth generation (BC5 ) hybrids. Allele-specific expression was determined for all genes and assigned to parental allele inheritance for each backcross hybrid individual. Results showed that genes residing in a 5.81 Mbp region on LG5 were exclusively expressed from the X. hellerii alleles in tumor-bearing BC1 hybrids. This observation indicates this region is consistently homozygous for X. hellerii alleles in tumor bearing animals, and therefore defines this region to be the R(Diff) locus. The R(Diff) locus harbors 164 gene models and includes the previously characterized R(Diff) candidate, cdkn2x. Twenty-one genes in the R(Diff) region show differential expression in the tumor samples compared to normal skin tissue. These results further characterize the R(Diff) locus and suggest tumor suppression may require a multigenic region rather than a single gene variant. Differences in gene expression between tumor and normal skin tissue in this region may indicate interactions among several genes are required for backcross hybrid melanoma development.

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.

Molecular genetic response to varied wavelengths of light in Xiphophorus maculatus skin.

Xiphophorus fishes represent a model often utilized to study UVB induced tumorigenesis. Recently, varied genetic responses to UVB exposure have been documented in the skin of female and male Xiphophorus, as have differences in UVB response in the skin of different parental species and for interspecies hybrids produced from crossing them. Additionally, it has been shown that exposure to "cool white" fluorescent light induces a shift in the genetic profiles of Xiphophorus skin that is nearly as robust as the UVB response, but involves a fundamentally different set of genes. Given these results and the use of Xiphophorus interspecies hybrids as an experimental model for UVB inducible melanoma, it is of interest to characterize genes that may be transcriptionally modulated in a wavelength specific manner. The global molecular genetic response of skin upon exposure of the intact animal to specific wavelengths of light has not been investigated. Herein, we report results of RNA-Seq experiments from the skin of male Xiphophorus maculatus Jp 163 B following exposure to varied 50nm wavelengths of light ranging from 300-600nm. We identify two specific wavelength regions, 350-400nm (88 genes) and 500-550nm (276 genes), that exhibit transcriptional modulation of a significantly greater number of transcripts than any of the other 50nm regions in the 300-600nm range. Observed functional sets of genes modulated within these two transcriptionally active light regions suggest different mechanisms of gene modulation.

Sex-specific molecular genetic response to UVB exposure in Xiphophorus maculatus skin.

In both Xiphophorus fishes and humans, males are reported to have a higher incidence of melanoma than females. To better understand sex-specific differences in the molecular genetic response to UVB, we performed RNA-Seq experiments in skin of female and male Xiphophorus maculatus Jp 163 B following UVB doses of 8 or 16kJ/m(2) exposure. Male X. maculatus differentially express a significantly larger number of transcripts following exposure to 16kJ/m(2) UVB (1293 genes) compared to 8kJ/m(2) UVB (324 genes). Female skin showed differential gene expression in a larger number of transcripts following 8kJ/m(2) UVB (765) than did males; however, both females and males showed similar numbers of differentially expressed genes at 16kJ/m(2) UVB (1167 and1293, respectively). Although most modulated transcripts after UVB exposure represented the same dominant pathways in both females and males (e.g., DNA repair, circadian rhythm, and fatty acid biosynthesis), we identified genes in several pathways that exhibited opposite modulation in female vs. male skin (e.g., synaptic development, cell differentiation, wound healing, and glucose metabolism). The oppositely modulated genes appear related through uncoupling protein 3 (UCP3) that is involved with the regulation of fatty acid oxidation and serves to balance glucose and lipid metabolism. Overall, these results identify gender-specific differences in UVB-induced genetic profiles in the skin of females and males and show female and male X. maculatus respond to UVB differently through pathways involved in reactive oxygen species, wound healing, and energy homeostasis.

Cortisol release in response to UVB exposure in Xiphophorus fish.

Xiphophorus fishes are comprised of 26 known species. Interspecies hybridization between select species has been utilized to produce experimental models to study melanoma development. Xiphophorus melanoma induction protocols utilize ultraviolet light (UVB) to induce DNA damage and associated downstream tumorigenesis. However, the impact of induced stress caused by the UVB treatment of the experimental animals undergoing tumor induction protocols has not been assessed. Stress is an adaptive physiological response to excessive or unpredictable environmental stimuli. The stress response in fishes may be measured by an assay of cortisol released into the water. Here, we present results from investigations of stress response during an experimental treatment and UVB exposure in Xiphophorus maculatus Jp 163 B, Xiphophorus couchianus, and F1 interspecies hybrids produced from the mating X. maculatus Jp 163 B×X. couchianus. Overall, cortisol release rates for males and females after UVB exposure showed no statistical differences. At lower UVB doses (8 and 16kJ/m(2)), X. couchianus exhibited 2 fold higher levels of DNA damage then either X. maculatus or the F1 hybrid. However, based on the cortisol release rates, none of the fish types tested induced a primary stress response at the UVB lower doses (8 and 16kJ/m(2)). In contrast, at a very high UVB dose (32kJ/m(2)) both X. maculatus and the F1 hybrid showed a 5 fold increase in the cortisol release rate. To determine the effect of pigmentation on UVB induced stress, wild type and albino Xiphophorus hellerii were exposed to UVB (32kJ/m(2)). Albino X. hellerii exhibited 3.7 fold increase in the cortisol release while wild type X. hellerii did not exhibit a significant cortisol response to UVB. Overall, the data suggest the rather low UVB doses often employed in tumor induction protocols do not induce a primary stress response in Xiphophorus fishes.

UVB-induced gene expression in the skin of Xiphophorus maculatus Jp 163 B.

Xiphophorus fish and interspecies hybrids represent long-standing models to study the genetics underlying spontaneous and induced tumorigenesis. The recent release of the Xiphophorus maculatus genome sequence will allow global genetic regulation studies of genes involved in the inherited susceptibility to UVB-induced melanoma within select backcross hybrids. As a first step toward this goal, we report results of an RNA-Seq approach to identify genes and pathways showing modulated transcription within the skin of X. maculatus Jp 163 B upon UVB exposure. X. maculatus Jp 163 B were exposed to various doses of UVB followed by RNA-Seq analysis at each dose to investigate overall gene expression in each sample. A total of 357 genes with a minimum expression change of 4-fold (p-adj<0.05) were identified as responsive to UVB. The molecular genetic response of Xiphophorus skin to UVB exposure permitted assessment of; (1) the basal expression level of each transcript for each skin sample, (2) the changes in expression levels for each gene in the transcriptome upon exposure to increasing doses of UVB, and (3) clusters of genes that exhibit similar patterns of change in expression upon UVB exposure. These data provide a foundation for understanding the molecular genetic response of fish skin to UVB exposure.

Characterization and differential expression of CPD and 6-4 DNA photolyases in Xiphophorus species and interspecies hybrids.

Among the many Xiphophorus interspecies hybrid tumor models are those that exhibit ultraviolet light (UVB) induced melanoma. In previous studies, assessment of UVB induced DNA damage and nucleotide excision DNA repair has been performed in parental lines and interspecies hybrids. Species and hybrid specific differences in the levels of DNA damage induced and the dark repair rates for cyclobutane pyrimidine dimers (CPDs) and 6-4 pyrimidine pyrimidine photoproducts (6-4PPs) have been reported. However, UVB induced DNA lesions in Xiphophorus fishes are thought to primarily be repaired via light dependent CPD and 6-4PP specific photolyases. Photolyases are of evolutionary interest since they are ancient and presumably function solely to ameliorate the deleterious effects of UVB exposure. Herein, we report results from detailed studies of CPD and 6-4PP photolyase gene expression within several Xiphophorus tissues. We determined photolyase gene expression patterns before and after exposure to fluorescent light in X. maculatus, X. couchianus, and for F1 interspecies hybrids produced from crossing these two parental lines (X. maculatus Jp 163 B×X. couchianus). We present novel results showing these two photolyase genes exhibit species, tissue, and hybrid-specific differences in basal and light induced gene expression.