Molecular genetic response of Xiphophorus maculatus-X. couchianus interspecies hybrid skin to UVB exposure

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.

Validity of Xiphophorus fish as models for human disease

Platyfish and swordtails of the genus Xiphophorus provide a well-established model for melanoma research and have become well known for this feature. Recently, modelling approaches for other human diseases in Xiphophorus have been developed or are emerging. This Review provides a comprehensive summary of these models and discusses how findings from basic biological and molecular studies and their translation to medical research demonstrate that Xiphophorus models have face, construct and predictive validity for studying a broad array of human diseases. These models can thus improve our understanding of disease mechanisms to benefit patients.

m6A mRNA Methylation Analysis Provides Novel Insights into Pigmentation in Sheep Skin

N6-methyladenosine (m⁶A) is the most universal post-transcriptional modification of mRNA which may play important roles in verious species. However, the potential roles of m⁶A in the pigmentation of skin are not completely understood. To explore the role of m6A modification in pigmentation of sheep skin, we used MeRIP-seq and RNA-seq to profile the skin transcriptome in black and white coat color (n=3). Our results showed that an average of 7701 m⁶A peaks were obtained for all samples and the average length was 305.89 bp. The GGACUU sequence was the most enrichment motif and shared in black skin and white skin. The m⁶A peaks were mainly enriched in the CDS, 3'UTR and 5'UTR, especially in CDS region near the stop codon of the transcript. 235 significantly differential peaks were found in black skin vs. white skin. The KEGG signaling pathways of downregulated and upregulated m⁶A peaks were mainly enriched in AGE-RAGE signaling pathway in diabetic complications, Viral carcinogenesis, Transcriptional misregulation in cancer, ABC transporters, Basal transcription factors and Thyroid hormone synthesis (P value <0.05). For RNA-seq, 71 differently expressed genes (DEGs) were scanned in black skin vs. white skin. DEGs were significantly enriched in tyrosine metabolism, melanogenesis, neuroactive ligand-receptor interaction pathway (P value <0.05). Combined m⁶A-seq and RNA-seq analysis showed that the hyper-up genes and hypo-up genes were both enriched in ErbB signaling pathway (P value <0.05). In conclusion, it provide a basis for further research into the functions of m⁶A methylation modifications in pigmentation.

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.

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.

Integration of eQTL and a Single-Cell Atlas in the Human Eye Identifies Causal Genes for Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of vision loss. To better understand disease pathogenesis and identify causal genes in GWAS loci for AMD risk, we present a comprehensive database of human retina and retinal pigment epithelium (RPE). Our database comprises macular and non-macular RNA sequencing (RNA-seq) profiles from 129 donors, a genome-wide expression quantitative trait loci (eQTL) dataset that includes macula-specific retina and RPE/choroid, and single-nucleus RNA-seq (NucSeq) from human retina and RPE with subtype resolution from more than 100,000 cells. Using NucSeq, we find enriched expression of AMD candidate genes in RPE cells. We identify 15 putative causal genes for AMD on the basis of co-localization of genetic association signals for AMD risk and eye eQTL, including the genes TSPAN10 and TRPM1. These results demonstrate the value of our human eye database for elucidating genetic pathways and potential therapeutic targets for ocular diseases.

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.

Gene expression variation and parental allele inheritance in a Xiphophorus interspecies hybridization model

Understanding the genetic mechanisms underlying segregation of phenotypic variation through successive generations is important for understanding physiological changes and disease risk. Tracing the etiology of variation in gene expression enables identification of genetic interactions, and may uncover molecular mechanisms leading to the phenotypic expression of a trait, especially when utilizing model organisms that have well-defined genetic lineages. There are a plethora of studies that describe relationships between gene expression and genotype, however, the idea that global variations in gene expression are also controlled by genotype remains novel. Despite the identification of loci that control gene expression variation, the global understanding of how genome constitution affects trait variability is unknown. To study this question, we utilized Xiphophorus fish of different, but tractable genetic backgrounds (inbred, F₁ interspecies hybrids, and backcross hybrid progeny), and measured each individual’s gene expression concurrent with the degrees of inter-individual expression variation. We found, (a) F₁ interspecies hybrids exhibited less variability than inbred animals, indicting gene expression variation is not affected by the fraction of heterozygous loci within an individual genome, and (b), that mixing genotypes in backcross populations led to higher levels of gene expression variability, supporting the idea that expression variability is caused by heterogeneity of genotypes of cis or trans loci. In conclusion, heterogeneity of genotype, introduced by inheritance of different alleles, accounts for the largest effects on global phenotypical variability.

Phenotypical variability is a multi-factorial phenomenon. Although it has been shown that inheriting certain gene is associated with lower phenotypical variability, how genome complexity affect phenotypical variability is still unclear. To study this question, we used inbred Xiphophorus fish, backcross interspecies hybrids, and F₁ interspecies hybrids between select Xiphophorus species to model genetic composition with minimum, medium, and maximum heterozygosity respectively, and measured their global gene expression variability. We found gene expression variation is not affected by the percentage of heterozygous loci in individual genome, but instead related to heterogeneity of genotype at local or remote loci.

Functional implications of Rab27 GTPases in Cancer

Background

The Rab27 family of small GTPases promotes the progression of breast cancer, melanoma, and other human cancers. In this review, we discuss the role of Rab27 GTPases in cancer progression and the potential applications of these targets in cancer treatment.

Main body

Elevated expression of Rab27 GTPases is associated with poor prognosis and cancer metastasis. Moreover, these GTPases govern a variety of oncogenic functions, including cell proliferation, cell motility, and chemosensitivity. In addition, small GTPases promote tumor growth and metastasis by enhancing exosome secretion, which alters intracellular microRNA levels, signaling molecule expression, and the tumor microenvironment.

Conclusion

Rab27 GTPases may have applications as prognostic markers and therapeutic targets in cancer treatment.

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.

Characterization of basal gene expression trends over a diurnal cycle in Xiphophorus maculatus skin, brain and liver

Evolutionarily conserved diurnal circadian mechanisms maintain oscillating patterns of gene expression based on the day-night cycle. Xiphophorus fish have been used to evaluate transcriptional responses after exposure to various light sources and it was determined that each source incites distinct genetic responses in skin tissue. However, basal expression levels of genes that show oscillating expression patterns in day-night cycle, may affect the outcomes of such experiments, since basal gene expression levels at each point in the circadian path may influence the profile of identified light responsive genes. Lack of knowledge regarding diurnal fluctuations in basal gene expression patterns may confound the understanding of genetic responses to external stimuli (e.g., light) since the dynamic nature of gene expression implies animals subjected to stimuli at different times may be at very different stages within the continuum of genetic homeostasis. We assessed basal gene expression changes over a 24-hour period in 200 select Xiphophorus gene targets known to transcriptionally respond to various types of light exposure. We identified 22 genes in skin, 36 genes in brain and 28 genes in liver that exhibit basal oscillation of expression patterns. These genes, including known circadian regulators, produced the expected expression patterns over a 24-hour cycle when compared to circadian regulatory genes identified in other species, especially human and other vertebrate animal models. Our results suggest the regulatory network governing diurnal oscillating gene expression is similar between Xiphophorus and other vertebrates for the three Xiphophorus organs tested. In addition, we were able to categorize light responsive gene sets in Xiphophorus that do, and do not, exhibit circadian based oscillating expression patterns.

Fluorescent light exposure incites acute and prolonged immune responses in zebrafish (Danio rerio) skin

Artificial light produces an emission spectrum that is considerably different than the solar spectrum. Artificial light has been shown to affect various behavior and physiological processes in vertebrates. However, there exists a paucity of data regarding the molecular genetic effects of artificial light exposure. Previous studies showed that one of the commonly used fluorescent light source (FL; 4100K or "cool white") can affect signaling pathways related to maintenance of circadian rhythm, cell cycle progression, chromosome segregation, and DNA repair/recombination in the skin of male Xiphophorus maculatus. These observations raise questions concerning the kinetics of the FL induced gene expression response, and which biological functions become modulated at various times after light exposure. To address these questions, we exposed zebrafish to 4100K FL and utilized RNA-Seq to assess gene expression changes in skin at various times (1 to 12h) after FL exposure. We found 4100K FL incites a robust early (1-2h) transcriptional response, followed by a more protracted late response (i.e., 4-12h). The early transcriptional response involves genes associated with cell migration/infiltration and cell proliferation as part of an overall increase in immune function and inflammation. The protracted late transcriptional response occurs within gene sets predicted to maintain and perpetuate the inflammatory response, as well as suppression of lipid, xenobiotic, and melatonin metabolism.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4735-5) contains supplementary material, which is available to authorized users.

Clonal polymorphism and high heterozygosity in the celibate genome of the Amazon molly

The extreme rarity of asexual vertebrates in nature is generally explained by genomic decay due to absence of meiotic recombination, thus leading to extinction of such lineages. We explore features of a vertebrate asexual genome, the Amazon molly, Poecilia formosa, and find few signs of genetic degeneration but unique genetic variability and ongoing evolution. We uncovered a substantial clonal polymorphism and, as a conserved feature from its interspecific hybrid origin, a 10-fold higher heterozygosity than in the sexual parental species. These characteristics seem to be a principal reason for the unpredicted fitness of this asexual vertebrate. Our data suggest that asexual vertebrate lineages are scarce not because they are at a disadvantage, but because the genomic combinations required to bypass meiosis and to make up a functioning hybrid genome are rarely met in nature.

Allele-specific expression at the androgen receptor alpha gene in a hybrid unisexual fish, the Amazon molly (Poecilia formosa)

The all-female Amazon molly (Poecilia formosa) is the result of a hybridization of the Atlantic molly (P. mexicana) and the sailfin molly (P. latipinna) approximately 120,000 years ago. As a gynogenetic species, P. formosa needs to copulate with heterospecific males including males from one of its bisexual ancestral species. However, the sperm only triggers embryogenesis of the diploid eggs. The genetic information of the sperm donor typically will not contribute to the next generation of P. formosa. Hence, P. formosa possesses generally one allele from each of its ancestral species at any genetic locus. This raises the question whether both ancestral alleles are equally expressed in P. formosa. Allele-specific expression (ASE) has been previously assessed in various organisms, e.g., human and fish, and ASE was found to be important in the context of phenotypic variability and disease. In this study, we utilized Real-Time PCR techniques to estimate ASE of the androgen receptor alpha (arα) gene in several distinct tissues of Amazon mollies. We found an allelic bias favoring the maternal ancestor (P. mexicana) allele in ovarian tissue. This allelic bias was not observed in the gill or the brain tissue. Sequencing of the promoter regions of both alleles revealed an association between an Indel in a known CpG island and differential expression. Future studies may reveal whether our observed cis-regulatory divergence is caused by an ovary-specific trans-regulatory element, preferentially activating the allele of the maternal ancestor.

Transcriptome assembly and candidate genes involved in nutritional programming in the swordtail fish Xiphophorus multilineatus

BACKGROUND

Nutritional programming takes place in early development. Variation in the quality and/or quantity of nutrients in early development can influence long-term health and viability. However, little is known about the mechanisms of nutritional programming. The live-bearing fish Xiphophorus multilineatus has the potential to be a new model for understanding these mechanisms, given prior evidence of nutritional programming influencing behavior and juvenile growth rate. We tested the hypotheses that nutritional programming would influence behaviors involved in energy homeostasis as well gene expression in X. multilineatus.

METHODS

We first examined the influence of both juvenile environment (varied in nutrition and density) and adult environment (varied in nutrition) on behaviors involved in energy acquisition and energy expenditure in adult male X. multilineatus. We also compared the behavioral responses across the genetically influenced size classes of males. Males stop growing at sexual maturity, and the size classes of can be identified based on phenotypes (adult size and pigment patterns). To study the molecular signatures of nutritional programming, we assembled a de novo transcriptome for X. multilineatus using RNA from brain, liver, skin, testis and gonad tissues, and used RNA-Seq to profile gene expression in the brains of males reared in low quality (reduced food, increased density) and high quality (increased food, decreased density) juvenile environments.

RESULTS

We found that both the juvenile and adult environments influenced the energy intake behavior, while only the adult environment influenced energy expenditure. In addition, there were significant interactions between the genetically influenced size classes and the environments that influenced energy intake and energy expenditure, with males from one of the four size classes (Y-II) responding in the opposite direction as compared to the other males examined. When we compared the brains of males of the Y-II size class reared in a low quality juvenile environment to males from the same size class reared in high quality juvenile environment, 131 genes were differentially expressed, including metabolism and appetite master regulator agrp gene.

DISCUSSION

Our study provides evidence for nutritional programming in X. multilineatus, with variation across size classes of males in how juvenile environment and adult diet influences behaviors involved in energy homeostasis. In addition, we provide the first transcriptome of X. multilineatus, and identify a group of candidate genes involved in nutritional programming.

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 (BC₁ ) hybrids and 13 fifth generation (BC₅ ) 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 BC₁ 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.

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.