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, F1 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) F1 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.

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.

Analysis of the putative tumor suppressor gene cdkn2ab in pigment cells and melanoma of Xiphophorus and medaka.

In humans, the CDKN2A locus encodes two transcripts, INK4A and ARF. Inactivation of either one by mutations or epigenetic changes is a frequent signature of malignant melanoma and one of the most relevant entry points for melanomagenesis. To analyze whether cdkn2ab, the fish ortholog of CDKN2A, has a similar function as its human counterpart, we studied its action in fish models for human melanoma. Overexpression of cdkn2ab in a Xiphophorus melanoma cell line led to decreased proliferation and induction of a senescence-like phenotype, indicating a melanoma-suppressive function analogous to mammals. Coexpression of Xiphophorus cdkn2ab in medaka transgenic for the mitfa:xmrk melanoma-inducing gene resulted in full suppression of melanoma development, whereas CRISPR/Cas9 knockout of cdkn2ab resulted in strongly enhanced tumor growth. In summary, this provides the first functional evidence that cdkn2ab acts as a potent tumor suppressor gene in fish melanoma models.

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.

Transcriptional control analyses of the Xiphophorus melanoma oncogene.

Melanoma development in interspecific hybrids of Xiphophorus is induced by the overexpression of the mutationally activated receptor tyrosine kinase Xmrk in pigment cells. Based on the melanocyte specificity of the transcriptional upregulation, a pigment cell-specific promoter region was postulated for xmrk, the activity of which is controlled in healthy purebred fish by the molecularly still unidentified regulator locus R. However, as yet the xmrk promoter region is still poorly characterized. In order to contribute to a better understanding of xmrk expression regulation, we performed a functional analysis of the entire putative gene regulatory region of the oncogene using conventional plasmid-based reporter systems as well as a newly established method employing BAC-derived luciferase reporter constructs in melanoma and non-melanoma cell lines. Using the melanocyte-specific mitfa promoter as control, we could demonstrate that our in vitro system is able to reliably monitor regulation of transcription through cell type-specific regulatory sequences. We found that sequences within 200kb flanking the xmrk oncogene do not lead to any specific transcriptional activation in melanoma compared to control cells. Hence, xmrk reporter constructs fail to faithfully reproduce the endogenous transcriptional regulation of the oncogene. Our data therefore strongly indicate that the melanocyte-specific transcription of xmrk is not the consequence of pigment cell-specific cis-regulatory elements in the promoter region. This hints at additional regulatory mechanisms involved in transcriptional control of the oncogene, thereby suggesting a key role for epigenetic mechanisms in oncogenic xmrk overexpression and thereby in tumor development in Xiphophorus.

Burst BMP triggered receptor kinase activity drives Smad1 mediated long-term target gene oscillation in C2C12 cells.

Bone Morphogenetic Proteins (BMPs) are important growth factors that regulate many cellular processes. During embryogenesis they act as morphogens and play a critical role during organ development. They influence cell fates via concentration-gradients in the embryos where cells transduce this extracellular information into gene expression profiles and cell fate decisions. How receiving cells decode and quantify BMP2/4 signals is hardly understood. There is little data on the quantitative relationships between signal input, transducing molecules, their states and location, and ultimately their ability to integrate graded systemic inputs and generate qualitative responses. Understanding this signaling network on a quantitative level should be considered a prerequisite for efficient pathway modulation, as the BMP pathway is a prime target for therapeutic invention. Hence, we quantified the spatial distribution of the main signal transducer of the BMP2/4 pathway in response to different types and levels of stimuli in c2c12 cells. We found that the subcellular localization of Smad1 is independent of ligand concentration. In contrast, Smad1 phosphorylation levels relate proportionally to BMP2 ligand concentrations and they are entirely located in the nucleus. Interestingly, we found that BMP2 stimulates target gene expression in non-linear, wave-like forms. Amplitudes showed a clear concentration-dependency, for sustained and transient stimulation. We found that even burst-stimulation triggers gene-expression wave-like modulations that are detectable for at least 30 h. Finally, we show here that target gene expression oscillations depend on receptor kinase activity, as the kinase drives further expression pulses without receptor reactivation and the target gene expression breaks off after inhibitor treatment in c2c12 cells.

Expression regulation triggers oncogenicity of xmrk alleles in the Xiphophorus melanoma system.

The Xiphophorus melanoma model has gained attention in biomedical research as a genetic model for tumor formation. Melanoma development in interspecific hybrids of Xiphophorus is connected to pigment cell specific overexpression of the mutationally activated receptor tyrosine kinase Xmrk. In purebred fish the oncogenic function of xmrk is suppressed by a so far unknown regulator locus R. To test the hypothesis that R is involved in transcriptional regulation of xmrk and consequently acts upstream of the xmrk signal, we performed a quantitative analysis of xmrk transcript levels in normal and melanoma tissues of different Xiphophorus genotypes carrying either a highly tumorigenic or a non-tumorigenic xmrk allele. Our results demonstrate that expression of the tumorigenic xmrk allele is highly increased in malignant melanomas compared to benign lesions, macromelanophore spots, and healthy skin. Transcription of the non-tumorigenic xmrk allele in pigment cells, in contrast, is not influenced by the presence or absence of R. These findings strongly indicate that differential transcriptional regulation of the xmrk promoter determines the tumorigenic potential of xmrk alleles in the Xiphophorus melanoma system, thereby supporting the hypothesis that R suppresses the oncogenic function of xmrk on the level of transcriptional control.