KRAS-retroviral fusion transcripts and gene amplification in arsenic-transformed, human prostate CAsE-PE cancer cells.

CAsE-PE cells are an arsenic-transformed, human prostate epithelial line containing oncogenic mutations in KRAS compared to immortalized, normal KRAS parent cells, RWPE-1. We previously reported increased copy number of mutated KRAS in CAsE-PE cells, suggesting gene amplification. Here, KRAS flanking genomic and transcriptomic regions were sequenced in CAsE-PE cells for insight into KRAS amplification. Comparison of DNA-Seq and RNA-Seq showed increased reads from background aligning to all KRAS exons in CAsE-PE cells, while a uniform DNA-Seq read distribution occurred in RWPE-1 cells with normal transcript expression. We searched for KRAS fusions in DNA and RNA sequencing data finding a portion of reads aligning to KRAS and viral sequence. After generation of cDNA from total RNA, short and long KRAS probes were generated to hybridize cDNA and KRAS enriched fragments were PacBio sequenced. More KRAS reads were captured from CAsE-PE cDNA versus RWPE-1 by each probe set. Only CAsE-PE cDNA showed KRAS viral fusion transcripts, primarily mapping to LTR and endogenous retrovirus sequences on either 5'- or 3'-ends of KRAS. Most KRAS viral fusion transcripts contained 4 to 6 exons but some PacBio sequences were in unusual orientations, suggesting viral insertions within the gene body. Additionally, conditioned media was extracted for potential retroviral particles. RNA-Seq of culture media isolates identified KRAS retroviral fusion transcripts in CAsE-PE media only. Truncated KRAS transcripts suggested multiple retroviral integration sites occurred within the KRAS gene producing KRAS retroviral fusions of various lengths. Findings suggest activation of endogenous retroviruses in arsenic carcinogenesis should be explored.

Arsenite malignantly transforms human prostate epithelial cells in vitro by gene amplification of mutated KRAS.

Inorganic arsenic is an environmental human carcinogen of several organs including the urinary tract. RWPE-1 cells are immortalized, non-tumorigenic, human prostate epithelia that become malignantly transformed into the CAsE-PE line after continuous in vitro exposure to 5µM arsenite over a period of months. For insight into in vitro arsenite transformation, we performed RNA-seq for differential gene expression and targeted sequencing of KRAS. We report >7,000 differentially expressed transcripts in CAsE-PE cells compared to RWPE-1 cells at >2-fold change, q<0.05 by RNA-seq. Notably, KRAS expression was highly elevated in CAsE-PE cells, with pathway analysis supporting increased cell proliferation, cell motility, survival and cancer pathways. Targeted DNA sequencing of KRAS revealed a mutant specific allelic imbalance, 'MASI', frequently found in primary clinical tumors. We found high expression of a mutated KRAS transcript carrying oncogenic mutations at codons 12 and 59 and many silent mutations, accompanied by lower expression of a wild-type allele. Parallel cultures of RWPE-1 cells retained a wild-type KRAS genotype. Copy number analysis and sequencing showed amplification of the mutant KRAS allele. KRAS is expressed as two splice variants, KRAS4a and KRAS4b, where variant 4b is more prevalent in normal cells compared to greater levels of variant 4a seen in tumor cells. 454 Roche sequencing measured KRAS variants in each cell type. We found KRAS4a as the predominant transcript variant in CAsE-PE cells compared to KRAS4b, the variant expressed primarily in RWPE-1 cells and in normal prostate, early passage, primary epithelial cells. Overall, gene expression data were consistent with KRAS-driven proliferation pathways found in spontaneous tumors and malignantly transformed cell lines. Arsenite is recognized as an important environmental carcinogen, but it is not a direct mutagen. Further investigations into this in vitro transformation model will focus on genomic events that cause arsenite-mediated mutation and overexpression of KRAS in CAsE-PE cells.

HAfTs are novel lncRNA transcripts from aflatoxin exposure.

The transcriptome can reveal insights into precancer biology. We recently conducted RNA-Seq analysis on liver RNA from male rats exposed to the carcinogen, aflatoxin B1 (AFB1), for 90 days prior to liver tumor onset. Among >1,000 differentially expressed transcripts, several novel, unannotated Cufflinks-assembled transcripts, or HAfTs (Hepatic Aflatoxin Transcripts) were found. We hypothesized PCR-cloning and RACE (rapid amplification of cDNA ends) could further HAfT identification. Sanger data was obtained for 6 transcripts by PCR and 16 transcripts by 5'- and 3'-RACE. BLAST alignments showed, with two exceptions, HAfT transcripts were lncRNAs, >200nt without apparent long open reading frames. Six rat HAfT transcripts were classified as 'novel' without RefSeq annotation. Sequence alignment and genomic synteny showed each rat lncRNA had a homologous locus in the mouse genome and over half had homologous loci in the human genome, including at least two loci (and possibly three others) that were previously unannotated. While HAfT functions are not yet clear, coregulatory roles may be possible from their adjacent orientation to known coding genes with altered expression that include 8 HAfT-gene pairs. For example, a unique rat HAfT, homologous to Pvt1, was adjacent to known genes controlling cell proliferation. Additionally, PCR and RACE Sanger sequencing showed many alternative splice variants and refinements of exon sequences compared to Cufflinks assembled transcripts and gene prediction algorithms. Presence of multiple splice variants and short tandem repeats found in some HAfTs may be consequential for secondary structure, transcriptional regulation, and function. In summary, we report novel, differentially expressed lncRNAs after exposure to the genotoxicant, AFB1, prior to neoplastic lesions. Complete cloning and sequencing of such transcripts could pave the way for a new set of sensitive and early prediction markers for chemical hepatocarcinogens.