Comparison of Xiphophorus and human melanoma transcriptomes reveals conserved pathway interactions

Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

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.

Xmrks the Spot: Fish Models for Investigating Epidermal Growth Factor Receptor Signaling in Cancer Research

Studies conducted in several fish species, e.g., Xiphophorus hellerii (green swordtail) and Xiphophorus maculatus (southern platyfish) crosses, Oryzias latipes (medaka), and Danio rerio (zebrafish), have identified an oncogenic role for the receptor tyrosine kinase, Xmrk, a gene product closely related to the human epidermal growth factor receptor (EGFR), which is associated with a wide variety of pathological conditions, including cancer. Comparative analyses of Xmrk and EGFR signal transduction in melanoma have shown that both utilize STAT5 signaling to regulate apoptosis and cell proliferation, PI3K to modulate apoptosis, FAK to control migration, and the Ras/Raf/MEK/MAPK pathway to regulate cell survival, proliferation, and differentiation. Further, Xmrk and EGFR may also modulate similar chemokine, extracellular matrix, oxidative stress, and microRNA signaling pathways in melanoma. In hepatocellular carcinoma (HCC), Xmrk and EGFR signaling utilize STAT5 to regulate cell proliferation, and Xmrk may signal through PI3K and FasR to modulate apoptosis. At the same time, both activate the Ras/Raf/MEK/MAPK pathway to regulate cell proliferation and E-cadherin signaling. Xmrk models of melanoma have shown that inhibitors of PI3K and MEK have an anti-cancer effect, and in HCC, that the steroidal drug, adrenosterone, can prevent metastasis and recover E-cadherin expression, suggesting that fish Xmrk models can exploit similarities with EGFR signal transduction to identify and study new chemotherapeutic drugs.

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.

Neural crest lineage analysis: from past to future trajectory

Since its discovery 150 years ago, the neural crest has intrigued investigators owing to its remarkable developmental potential and extensive migratory ability. Cell lineage analysis has been an essential tool for exploring neural crest cell fate and migration routes. By marking progenitor cells, one can observe their subsequent locations and the cell types into which they differentiate. Here, we review major discoveries in neural crest lineage tracing from a historical perspective. We discuss how advancing technologies have refined lineage-tracing studies, and how clonal analysis can be applied to questions regarding multipotency. We also highlight how effective progenitor cell tracing, when combined with recently developed molecular and imaging tools, such as single-cell transcriptomics, single-molecule fluorescence in situ hybridization and high-resolution imaging, can extend the scope of neural crest lineage studies beyond development to regeneration and cancer initiation.

Spontaneously occurring melanoma in animals and their relevance to human melanoma

In contrast to other cancer types, melanoma incidence has been increasing over the last 50 years, and while it still represents less than 5% of all cutaneous malignancies, melanoma accounts for the majority of skin cancer deaths, due to its propensity to metastasise. Whilst melanoma most commonly affects the skin, it can also arise in mucosal surfaces, the eye, and the brain. For new therapies to be developed, a better understanding of the genetic landscape, signalling pathways, and tumour–microenvironmental interactions is needed. This is where animal models are of critical importance. The mouse is the foremost used model of human melanoma. Arguably this is due to its plethora of benefits as a laboratory animal; however, it is important to note that unlike humans, melanocytes are not present at the dermal–epidermal junction in mice and mice do not develop melanoma without genetic manipulation. In contrast, there are numerous reports of animals that spontaneously develop melanoma, ranging from sharks and parrots to hippos and monkeys. In addition, several domesticated and laboratory‐bred animals spontaneously develop melanoma or UV‐induced melanoma, specifically, fish, opossums, pigs, horses, cats, and dogs. In this review, we look at spontaneously occurring animal ‘models’ of melanoma and discuss their relevance to the different types of melanoma found in humans. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland..

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.

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.