Retinoic acid receptor α activity in proximal tubules prevents kidney injury and fibrosis.

Chronic kidney disease (CKD) is characterized by a gradual loss of kidney function and affects ~13.4% of the global population. Progressive tubulointerstitial fibrosis, driven in part by proximal tubule (PT) damage, is a hallmark of late stages of CKD and contributes to the development of kidney failure, for which there are limited treatment options. Normal kidney development requires signaling by vitamin A (retinol), which is metabolized to retinoic acid (RA), an endogenous agonist for the RA receptors (RARα, ß, γ). RARα levels are decreased in a mouse model of diabetic nephropathy and restored with RA administration; additionally, RA treatment reduced fibrosis. We developed a mouse model in which a spatiotemporal (tamoxifen-inducible) deletion of RARα in kidney PT cells of adult mice causes mitochondrial dysfunction, massive PT injury, and apoptosis without the use of additional nephrotoxic substances. Long-term effects (3 to 4.5 mo) of RARα deletion include increased PT secretion of transforming growth factor ß1, inflammation, interstitial fibrosis, and decreased kidney function, all of which are major features of human CKD. Therefore, RARα's actions in PTs are crucial for PT homeostasis, and loss of RARα causes injury and a key CKD phenotype.

Novel genetically engineered mouse models for clear cell renal cell carcinoma.

Genetically engineered mouse models (GEMMs) are important immunocompetent models for research into the roles of individual genes in cancer and the development of novel therapies. Here we use inducible CRISPR-Cas9 systems to develop two GEMMs which aim to model the extensive chromosome p3 deletion frequently observed in clear cell renal cell carcinoma (ccRCC). We cloned paired guide RNAs targeting early exons of Bap1, Pbrm1, and Setd2 in a construct containing a Cas9D10A (nickase, hSpCsn1n) driven by tetracycline (tet)-responsive elements (TRE3G) to develop our first GEMM. The founder mouse was crossed with two previously established transgenic lines, one carrying the tet-transactivator (tTA, Tet-Off) and one with a triple-mutant stabilized HIF1A-M3 (TRAnsgenic Cancer of the Kidney, TRACK), both driven by a truncated, proximal tubule-specific γ-glutamyltransferase 1 (ggt or γGT) promoter, to create triple-transgenic animals. Our results indicate that this model (BPS-TA) induces low numbers of somatic mutations in Bap1 and Pbrm1 (but not in Setd2), known tumor suppressor genes in human ccRCC. These mutations, largely restricted to kidneys and testis, induced no detectable tissue transformation in a cohort of 13 month old mice (N = 10). To gain insights into the low frequencies of insertions and deletions (indels) in BPS-TA mice we analyzed wild type (WT, N = 7) and BPS-TA (N = 4) kidneys by RNAseq. This showed activation of both DNA damage and immune response, suggesting activation of tumor suppressive mechanisms in response to genome editing. We then modified our approach by generating a second model in which a ggt-driven, cre-regulated Cas9WT(hSpCsn1) was employed to introduce Bap1, Pbrm1, and Setd2 genome edits in the TRACK line (BPS-Cre). The BPS-TA and BPS-Cre lines are both tightly controlled in a spatiotemporal manner with doxycycline (dox) and tamoxifen (tam), respectively. In addition, whereas the BPS-TA line relies on paired guide RNAs (gRNAs), the BPS-Cre line requires only single gRNAs for gene perturbation. In the BPS-Cre we identified increased Pbrm1 gene-editing frequencies compared to the BPS-TA model. Whereas we did not detect Setd2 edits in the BPS-TA kidneys, we found extensive editing of Setd2 in the BPS-Cre model. Bap1 editing efficiencies were comparable between the two models. Although no gross malignancies were observed in our study, this is the first reported GEMM which models the extensive chromosome 3p deletion frequently observed in kidney cancer patients. Further studies are required (1) to model more extensive 3p deletions, e.g. impacting additional genes, and (2) to increase the cellular resolution, e.g. by employing single-cell RNAseq to ascertain the effects of specific combinatorial gene inactivation.

Combined Metabolomics and Genome-Wide Transcriptomics Analyses Show Multiple HIF1α-Induced Changes in Lipid Metabolism in Early Stage Clear Cell Renal Cell Carcinoma.

The accumulation of lipids is a hallmark of human clear cell renal cell carcinoma (ccRCC). Advanced ccRCC tumors frequently show increased lipid biosynthesis, but the regulation of lipid metabolism in early stage ccRCC tumors has not been studied. Here, we performed combined transcriptomics and metabolomics on a previously characterized transgenic mouse model (TRAnsgenic Cancer of the Kidney, TRACK) of early stage ccRCC. We found that in TRACK kidneys, HIF1α activation increases transcripts of lipid receptors (Cd36, ACVRL1), lipid storage genes (Hilpda and Fabp7), and intracellular levels of essential fatty acids, including linoleic acid and linolenic acid. Feeding the TRACK mice a high-fat diet enhances lipid accumulation in the kidneys. These results show that HIF1α increases the uptake and storage of dietary lipids in this early stage ccRCC model. By then analyzing early stage human ccRCC specimens, we found similar increases in CD36 transcripts and increases in linoleic and linolenic acid relative to normal kidney samples. CD36 mRNA levels decreased, while FASN transcript levels increased with increasing ccRCC tumor stage. These results suggest that an increase in the lipid biosynthesis pathway in advanced ccRCC tumors may compensate for a decreased capacity of these advanced ccRCCs to scavenge extracellular lipids.

Role of NADH Dehydrogenase (Ubiquinone) 1 Alpha Subcomplex 4-Like 2 in Clear Cell Renal Cell Carcinoma.

PURPOSE: We delineated the functions of the hypoxia-inducible factor-1α (HIF1α) target NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4-like 2 (NDUFA4L2) in clear cell renal cell carcinoma (ccRCC) and characterized NDUFA4L2 as a novel molecular target for ccRCC treatment. EXPERIMENTAL DESIGN: We evaluated normal kidney and ccRCC patient microarray and RNAseq data from Oncomine and The Cancer Genome Atlas for NDUFA4L2 mRNA levels and the clinical implications of high NDUFA4L2 expression. In addition, we examined normal kidney and ccRCC patient tissue samples, human ccRCC cell lines, and murine models of ccRCC for NDUFA4L2 mRNA and protein expression. Utilizing short hairpin RNA, we performed NDUFA4L2 knockdown experiments and analyzed the proliferation, clonogenicity, metabolite levels, cell structure, and autophagy in ccRCC cell lines in culture. RESULTS: We found that NDUFA4L2 mRNA and protein are highly expressed in ccRCC samples but undetectable in normal kidney tissue samples, and that NDUFA4L2 mRNA expression correlates with tumor stage and lower overall survival. In addition, we demonstrated that NDUFA4L2 is an HIF1α target in ccRCC and that NDUFA4L2 knockdown has a profound antiproliferative effect, alters metabolic pathways, and causes major stress in cultured RCC cells. CONCLUSIONS: Collectively, our data show that NDUFA4L2 is a novel molecular target for ccRCC treatment. Clin Cancer Res; 22(11); 2791-801. ©2016 AACR.

Analyses of the transcriptome and metabolome demonstrate that HIF1α mediates altered tumor metabolism in clear cell renal cell carcinoma.

Hypoxia inducible factor 1 alpha (HIF1α) is a transcription factor that is frequently stabilized and active in human clear cell renal cell carcinoma (ccRCC). We have found that constitutively active HIF1α is sufficient to cause neoplastic transformation in a murine model of ccRCC termed the TRACK model. RNA sequencing (RNAseq) and untargeted metabolomics analyses of samples from TRACK kidneys demonstrate that HIF1α activates the transcription of genes that cause increased glucose uptake, glycolysis, and lactate production, as well as a decrease in the flux of pyruvate entering the tricarboxylic acid (TCA) cycle and a decrease in oxidative phosphorylation; these changes are identical to those observed in human ccRCC samples. These studies show that a constitutively active HIF1α promotes tumorigenesis in TRACK mice by mediating a metabolic switch to aerobic glycolysis, i.e., the Warburg effect, and suggest that TRACK mice are a valid model to test novel therapies targeting metabolic changes to inhibit human ccRCC.

Genome-Wide Profiling of TRACK Kidneys Shows Similarity to the Human ccRCC Transcriptome.

UNLABELLED: Renal cell carcinoma (RCC) is the most common cancer arising from the kidney in adults, with clear cell RCC (ccRCC) representing the majority of all RCCs. Expression of a human HIF1α triple-mutant (P402A, P564A, and N803A) construct in the proximal tubule cells of C57BL/6 mice [TRAnsgenic model of Cancer of the Kidney (TRACK); ref. 1] mimics the histologic changes found in early stage human ccRCC. To better understand the genomic landscape, a high-throughput sequence analysis was performed with cDNA libraries (RNAseq) derived from TRACK transgenic positive (TG(+)) kidney cortex along with human ccRCC transcripts from the Oncomine and The Cancer Genome Atlas databases. Importantly, the expression profiles of TRACK TG(+) kidneys show significant similarities with those observed in human ccRCC, including increased expression of genes involved in glycolysis and the tricarboxylic acid cycle. Some of the transcripts overexpressed in both the TRACK mouse model and human ccRCC include ANKRD37, CA9, EGLN3, HK2, NDUFA4L2, and SLC16A3. These data suggest that constitutive activation of HIF1α in kidney proximal tubule cells transcriptionally reprograms the regulation of metabolic pathways in the kidney and that HIF1α is a major contributor to the altered metabolism observed in human ccRCC. IMPLICATIONS: TRACK (GGT-HIF1αM3) kidney mRNA profiles show similarities to human ccRCC transcriptome and phenotypes associated with the Warburg effect.

The role of HIF1α in renal cell carcinoma tumorigenesis.

UNLABELLED: The transcription factor HIF1α is implicated in the development of clear cell renal cell carcinoma (ccRCC). Although HIF1α was initially believed to be essential for ccRCC development, recent studies hypothesize an oncogenic role for HIF2α in ccRCC, but a tumor suppressor role for HIF1α, leading to uncertainty as to the precise roles of the different HIF transcription factors in this disease. Using evidence available from studies with human ccRCC cell lines, mouse xenografts, murine models of ccRCC, and human ccRCC specimens, we evaluate the roles of HIF1α and HIF2α in the pathogenesis of ccRCC. We present a convergence of clinical and mechanistic data supporting an important role for HIF1α in promoting tumorigenesis in a clinically important and large subset of ccRCC. This indicates that current understanding of the exact roles of HIF1α and HIF2α is incomplete and that further research is required to determine the diverse roles of HIF1α and HIF2α in ccRCC. KEY MESSAGES: The TRACK mouse ccRCC model with constitutively active HIF1α but not HIF2α expressed in proximal tubules develops RCC. HIF1α protein is expressed in the majority of human ccRCC specimens. Elevated HIF1α in ccRCC correlates with a worse prognosis. Many publications do not support a tumor suppressor role for HIF1α in ccRCC. HIF1α, but not HIF2α, is expressed in some types of cancer stem cells.

Activation of HIF2α in kidney proximal tubule cells causes abnormal glycogen deposition but not tumorigenesis.

Renal cell carcinoma (RCC) is the most common primary cancer arising from the kidney in adults, with clear cell renal cell carcinoma (ccRCC) representing approximately 75% of all RCCs. Increased expression of the hypoxia-induced factors-1α (HIF1α) and HIF2α has been suggested as a pivotal step in ccRCC carcinogenesis, but this has not been thoroughly tested. Here, we report that expression of a constitutively activated form of HIF2α (P405A, P530A, and N851A, named as HIF2αM3) in the proximal tubules of mice is not sufficient to promote ccRCC by itself, nor does it enhance HIF1αM3 oncogenesis when coexpressed with constitutively active HIF1αM3. Neoplastic transformation in kidneys was not detected at up to 33 months of age, nor was increased expression of Ki67 (MKI67), γH2AX (H2AFX), or CD70 observed. Furthermore, the genome-wide transcriptome of the transgenic kidneys does not resemble human ccRCC. We conclude that a constitutively active HIF2α is not sufficient to cause neoplastic transformation of proximal tubules, arguing against the idea that HIF2α activation is critical for ccRCC tumorigenesis.

Generation of a mouse model of Von Hippel-Lindau kidney disease leading to renal cancers by expression of a constitutively active mutant of HIF1α.

Renal cancers are highly aggressive and clinically challenging, but a transgenic mouse model to promote pathologic studies and therapeutic advances has yet to be established. Here, we report the generation of a transgenic mouse model of von Hippel-Lindau (VHL) renal cancer termed the TRACK model (transgenic model of cancer of the kidney). TRACK mice specifically express a mutated, constitutively active HIF1α in kidney proximal tubule (PT) cells. Kidney histologies displayed by TRACK mice are highly similar to histologies seen in patients with VHL disease, including areas of distorted tubular structure, cells with clear cytoplasm and increased glycogen and lipid deposition, multiple renal cysts, and early onset of clear cell renal cell carcinoma (ccRCC). Distorted tubules in TRACK mice exhibit higher levels of CA-IX, Glut1, and VEGF than tubules in nontransgenic control mice. Furthermore, these tubules exhibit increased numbers of endothelial cells, increased cell proliferation, and increased expression of the human ccRCC marker CD70(TNFSF7). Moreover, PT cells in kidney tubules from TRACK mice exhibit increased genomic instability, as monitored by elevated levels of γH2AX. Our findings establish that activated HIF1α in murine kidney PT cells is sufficient to promote cell proliferation, angiogenesis, genomic instability, and other phenotypic alterations characteristic of human VHL kidney disease, establishing the TRACK mouse as a valid preclinical model of human renal cell carcinoma.

Degradation of p53 by human Alphapapillomavirus E6 proteins shows a stronger correlation with phylogeny than oncogenicity.

BACKGROUND: Human Papillomavirus (HPV) E6 induced p53 degradation is thought to be an essential activity by which high-risk human Alphapapillomaviruses (alpha-HPVs) contribute to cervical cancer development. However, most of our understanding is derived from the comparison of HPV16 and HPV11. These two viruses are relatively distinct viruses, making the extrapolation of these results difficult. In the present study, we expand the tested strains (types) to include members of all known HPV species groups within the Alphapapillomavirus genus. PRINCIPAL FINDINGS: We report the biochemical activity of E6 proteins from 27 HPV types representing all alpha-HPV species groups to degrade p53 in human cells. Expression of E6 from all HPV types epidemiologically classified as group 1 carcinogens significantly reduced p53 levels. However, several types not associated with cancer (e.g., HPV53, HPV70 and HPV71) were equally active in degrading p53. HPV types within species groups alpha 5, 6, 7, 9 and 11 share a most recent common ancestor (MRCA) and all contain E6 ORFs that degrade p53. A unique exception, HPV71 E6 ORF that degraded p53 was outside this clade and is one of the most prevalent HPV types infecting the cervix in a population-based study of 10,000 women. Alignment of E6 ORFs identified an amino acid site that was highly correlated with the biochemical ability to degrade p53. Alteration of this amino acid in HPV71 E6 abrogated its ability to degrade p53, while alteration of this site in HPV71-related HPV90 and HPV106 E6s enhanced their capacity to degrade p53. CONCLUSIONS: These data suggest that the alpha-HPV E6 proteins' ability to degrade p53 is an evolved phenotype inherited from a most recent common ancestor of the high-risk species that does not always segregate with carcinogenicity. In addition, we identified an amino-acid residue strongly correlated with viral p53 degrading potential.

Identification of a novel human papillomavirus (HPV97) related to HPV18 and HPV45.

Human papillomavirus (HPV) type 97 was identified and the genome was cloned from cervicovaginal cells of a Costa Rican woman with a normal Pap smear. The HPV97 L1 open reading frame (ORF) was most closely related to HPV45 (84% identity) and HPV18 (79% identity), placing it into the high-risk alpha7 species. Ectopic expression of the HPV97 E6 and E7 proteins significantly decreased steady state p53 and pRb levels using an in vitro cotransfection assay, respectively. These data suggest that HPV97 shares a most recent common ancestor with HPV18 and HPV45 and should be evaluated in cancer specimens from different geographic populations.

Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses.

Human papillomavirus type 16 (HPV16) is the primary etiological agent of cervical cancer, the second most common cancer in women worldwide. Complete genomes of 12 isolates representing the major lineages of HPV16 were cloned and sequenced from cervicovaginal cells. The sequence variations within the open reading frames (ORFs) and noncoding regions were identified and compared with the HPV16R reference sequence. This whole-genome approach gives us unprecedented precision in detailing sequence-level changes that are under selection on a whole-viral-genome scale. Of 7,908 base pair nucleotide positions, 313 (4.0%) were variable. Within the 2,452 amino acids (aa) comprising 8 ORFs, 243 (9.9%) amino acid positions were variable. In order to investigate the molecular evolution of HPV16 variants, maximum likelihood models of codon substitution were used to identify lineages and amino acid sites under selective pressure. Five codon sites in the E5 (aa 48, 65) and E6 (aa 10, 14, 83) ORFs were demonstrated to be under diversifying selective pressure. The E5 ORF had the overall highest nonsynonymous/synonymous substitution rate (omega) ratio (M3 = 0.7965). The E2 gene had the next-highest omega ratio (M3 = 0.5611); however, no specific codons were under positive selection. These data indicate that the E6 and E5 ORFs are evolving under positive Darwinian selection and have done so in a relatively short time period. Whether response to selective pressure upon the E5 and E6 ORFs contributes to the biological success of HPV16, its specific biological niche, and/or its oncogenic potential remains to be established.

Codetection of a mixed population of candHPV62 containing wild-type and disrupted E1 open-reading frame in a 45-year-old woman with normal cytology.

We have cloned, sequenced, and characterized the complete genome of a novel human papillomavirus (HPV), candHPV62. During cloning, 2 candHPV62 viral isolates were recovered from a single cervical sample; 1 had all anticipated HPV open-reading frames (ORFs) intact, whereas the other exhibited an E1 frame-shift mutation. Further experiments indicated that the 2 strains were equivalent in abundance. It appears that an early mutation occurred within the E1 ORF, which was transcomplemented by an intact E1 protein. A search of the HPV database identified disruption of the E1 ORF in the cloned reference isolates of HPV16, HPV53, HPV56, and HPV72. These data suggest that disruption of the E1 ORF in genital HPVs is not uncommon.

Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix.

The distributions of human papillomavirus (HPV) types detected in cervical adenocarcinomas and squamous cell tumors differ. However, whether the distributions of intratypic HPV variants seen in these two histological forms of cervical disease differ is unknown. Our objective was to compare the distribution of HPV intratypic variants observed in squamous cell carcinomas (SCC) and cervical tumors of glandular origin (e.g., adenocarcinomas; AC) for two HPV types commonly observed in cervical tumors, HPV16 and HPV18. Participants in a multicenter case-control study of AC and SCC conducted in the eastern United States were studied. A total of 85 HPV16 and/or HPV18 positive individuals (31 diagnosed with AC, 43 diagnosed with SCC, and 11 population controls) were included. For HPV16-positive individuals, both the noncoding long control region and the E6 open reading frame were sequenced, and classified into phylogenetic-based lineage groups (European, Asian-American, African1, and African2). For HPV18-positive individuals, the long control region region only was sequenced and classified into known intratypic lineages (European, Asian-Amerindian, and African). The distribution of these different intratypic lineages among AC cases, SCC cases, and population controls was compared using standard methods. Non-European HPV16 and/or HPV18 intratypic variants were observed in 42% of ACs compared with 16% of SCCs and 18% of population controls (P = 0.04). Intratypic variants from the Asian-American lineage of HPV16 accounted for the differences seen between histological groups. The differences observed between AC and SCC cases were strongest for HPV16, and persisted in analysis restricted to Caucasian women, suggesting that the effect cannot be explained by differences in the ethnic make-up of AC versus SCC cases. Cervical AC and SCC differ not only with respect to the distribution of HPV types detected but also with respect to intratypic variants observed. Non-European HPV16 and/or HPV18 variants are commonly seen in AC. A possible hormonal mechanism is suggested to explain the observed findings.