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.

Identification of transcriptome SNPs between Xiphophorus lines and species for assessing allele specific gene expression within F₁ interspecies hybrids

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for ≈70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged ≈5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.