Methods to study xenografted human cancer in genetically diverse mice.

Xenografting human cancer tissues into mice to test new cures against cancers is critical for understanding and treating the disease. However, only a few inbred strains of mice are used to study cancers, and derivatives of mainly one strain, mostly NOD/ShiLtJ, are used for therapy efficacy studies. As it has been demonstrated when human cancer cell lines or patient-derived tissues (PDX) are xenografted into mice, the neoplastic cells are human but the supporting cells that comprise the tumor (the stroma) are from the mouse. Therefore, results of studies of xenografted tissues are influenced by the host strain. We previously published that when the same neoplastic cells are xenografted into different mouse strains, the pattern of tumor growth, histology of the tumor, number of immune cells infiltrating the tumor, and types of circulating cytokines differ depending on the strain. Therefore, to better comprehend the behavior of cancer in vivo, one must xenograft multiple mouse strains. Here we describe and report a series of methods that we used to reveal the genes and proteins expressed when the same cancer cell line, MDA-MB-231, is xenografted in different hosts. First, using proteomic analysis, we show how to use the same cell line in vivo to reveal the protein changes in the neoplastic cell that help it adapt to its host. Then, we show how different hosts respond molecularly to the same cell line. We also find that using multiple strains can reveal a more suitable host than those traditionally used for a "difficult to xenograft" PDX. In addition, using complex trait genetics, we illustrate a feasible method for uncovering the alleles of the host that support tumor growth. Finally, we demonstrate that Diversity Outbred mice, the epitome of a model of mouse-strain genetic diversity, can be xenografted with human cell lines or PDX using 2-deoxy-D-glucose treatment.

Structural bioinformatics enhances the interpretation of somatic mutations in KDM6A found in human cancers.

The histone demethylase KDM6A has recently elicited significant attention because its mutations are associated with a rare congenital disorder (Kabuki syndrome) and various types of human cancers. However, distinguishing KDM6A mutations that are deleterious to the enzyme and their underlying mechanisms of dysfunction remain to be fully understood. Here, we report the results from a multi-tiered approach evaluating the impact of 197 KDM6A somatic mutations using information derived from combining conventional genomics data with computational biophysics. This comprehensive approach incorporates multiple scores derived from alterations in protein sequence, structure, and molecular dynamics. Using this method, we classify the KDM6A mutations into 136 damaging variants (69.0%), 32 tolerated variants (16.2%), and 29 variants of uncertain significance (VUS, 14.7%), which is a significant improvement from the previous classification based on the conventional tools (over 40% VUS). We further classify the damaging variants into 15 structural variants (SV), 88 dynamic variants (DV), and 33 structural and dynamic variants (SDV). Comparison with variant scoring methods used in current clinical diagnosis guidelines demonstrates that our approach provides a more comprehensive evaluation of damaging potential and reveals mechanisms of dysfunction. Thus, these results should be taken into consideration for clinical assessment of the damaging potential of each mutation, as they provide hypotheses for experimental validation and critical information for the development of mutant-specific drugs to fight diseases caused by KDM6A dysfunctions.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis.

Disruptor of telomeric silencing 1-like (DOT1L) is the only non-SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time-dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1-AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L-AF9 complex are less extensive and highly dynamic, resembling a swivel-like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

KrasG12D induces changes in chromatin territories that differentially impact early nuclear reprogramming in pancreatic cells.

BACKGROUND: Pancreatic ductal adenocarcinoma initiation is most frequently caused by Kras mutations. RESULTS: Here, we apply biological, biochemical, and network biology methods to validate GEMM-derived cell models using inducible KrasG12D expression. We describe the time-dependent, chromatin remodeling program that impacts function during early oncogenic signaling. We find that the KrasG12D-induced transcriptional response is dominated by downregulated expression concordant with layers of epigenetic events. More open chromatin characterizes the ATAC-seq profile associated with a smaller group of upregulated genes and epigenetic marks. RRBS demonstrates that promoter hypermethylation does not account for the silencing of the extensive gene promoter network. Moreover, ChIP-Seq reveals that heterochromatin reorganization plays little role in this early transcriptional program. Notably, both gene activation and silencing primarily depend on the marking of genes with a combination of H3K27ac, H3K4me3, and H3K36me3. Indeed, integrated modeling of all these datasets shows that KrasG12D regulates its transcriptional program primarily through unique super-enhancers and enhancers, and marking specific gene promoters and bodies. We also report chromatin remodeling across genomic areas that, although not contributing directly to cis-gene transcription, are likely important for KrasG12D functions. CONCLUSIONS: In summary, we report a comprehensive, time-dependent, and coordinated early epigenomic program for KrasG12D in pancreatic cells, which is mechanistically relevant to understanding chromatin remodeling events underlying transcriptional outcomes needed for the function of this oncogene.

Inactivation of the Euchromatic Histone-Lysine N-Methyltransferase 2 Pathway in Pancreatic Epithelial Cells Antagonizes Cancer Initiation and Pancreatitis-Associated Promotion by Altering Growth and Immune Gene Expression Networks.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive, painful disease with a 5-year survival rate of only 9%. Recent evidence indicates that distinct epigenomic landscapes underlie PDAC progression, identifying the H3K9me pathway as important to its pathobiology. Here, we delineate the role of Euchromatic Histone-lysine N-Methyltransferase 2 (EHMT2), the enzyme that generates H3K9me, as a downstream effector of oncogenic KRAS during PDAC initiation and pancreatitis-associated promotion. EHMT2 inactivation in pancreatic cells reduces H3K9me2 and antagonizes Kras G12D -mediated acinar-to-ductal metaplasia (ADM) and Pancreatic Intraepithelial Neoplasia (PanIN) formation in both the Pdx1-Cre and P48 Cre/+ Kras G12D mouse models. Ex vivo acinar explants also show impaired EGFR-KRAS-MAPK pathway-mediated ADM upon EHMT2 deletion. Notably, Kras G12D increases EHMT2 protein levels and EHMT2-EHMT1-WIZ complex formation. Transcriptome analysis reveals that EHMT2 inactivation upregulates a cell cycle inhibitory gene expression network that converges on the Cdkn1a/p21-Chek2 pathway. Congruently, pancreas tissue from Kras G12D animals with EHMT2 inactivation have increased P21 protein levels and enhanced senescence. Furthermore, loss of EHMT2 reduces inflammatory cell infiltration typically induced during Kras G12D -mediated initiation. The inhibitory effect on Kras G12D -induced growth is maintained in the pancreatitis-accelerated model, while simultaneously modifying immunoregulatory gene networks that also contribute to carcinogenesis. This study outlines the existence of a novel KRAS-EHMT2 pathway that is critical for mediating the growth-promoting and immunoregulatory effects of this oncogene in vivo, extending human observations to support a pathophysiological role for the H3K9me pathway in PDAC.

Molecular mechanics and dynamic simulations of well-known Kabuki syndrome-associated KDM6A variants reveal putative mechanisms of dysfunction.

BACKGROUND: Kabuki syndrome is a genetic disorder that affects several body systems and presents with variations in symptoms and severity. The syndrome is named for a common phenotype of faces resembling stage makeup used in a Japanese traditional theatrical art named kabuki. The most frequent cause of this syndrome is mutations in the H3K4 family of histone methyltransferases while a smaller percentage results from genetic alterations affecting the histone demethylase, KDM6A. Because of the rare presentation of the latter form of the disease, little is known about how missense changes in the KDM6A protein sequence impact protein function. RESULTS: In this study, we use molecular mechanic and molecular dynamic simulations to enhance the annotation and mechanistic interpretation of the potential impact of eleven KDM6A missense variants found in Kabuki syndrome patients. These variants (N910S, D980V, S1025G, C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W, and R1351Q) are predicted to be pathogenic, likely pathogenic or of uncertain significance by sequence-based analysis. Here, we demonstrate, for the first time, that although Kabuki syndrome missense variants are found outside the functionally critical regions, they could affect overall function by significantly disrupting global and local conformation (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), chemical environment (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), and/or molecular dynamics of the catalytic domain (all variants). In addition, our approaches predict that many mutations, in particular C1153R, could allosterically disrupt the key enzymatic interactions of KDM6A. CONCLUSIONS: Our study demonstrates that the KDM6A Kabuki syndrome variants may impair histone demethylase function through various mechanisms that include altered protein integrity, local environment, molecular interactions and protein dynamics. Molecular dynamics simulations of the wild type and the variants are critical to gain a better understanding of molecular dysfunction. This type of comprehensive structure- and MD-based analyses should help develop improved impact scoring systems to interpret the damaging effects of variants in this protein and other related proteins as well as provide detailed mechanistic insight that is not currently predictable from sequence alone.

Genome-wide map of proximity linkage to renin proximal promoter in rat.

A challenge to understanding enhancer-gene relationships is that enhancers are not always sequentially close to the gene they regulate. Physical proximity mapping through sequencing can provide an unbiased view of the chromatin close to the proximal promoter of the renin gene ( Ren). Our objective was to determine genomic regions that physically interact with the renin proximal promoter, using two different genetic backgrounds, the Dahl salt sensitive and normotensive SS-13BN, which have been shown to have different regulation of plasma renin in vivo. The chromatin conformation capture method with sequencing focused at the Ren proximal promoter in rat-derived cardiac endothelial cells was used. Cells were fixed, chromatin close to the Ren promoter was captured, and fragments were sequenced. The clustering of mapped reads produced a genome-wide map of chromatin in contact with the Ren promoter. The largest number of contacts was found on chromosome 13, the chromosome with Ren, and contacts were found on all other chromosomes except chromosome X. These contacts were significantly enriched with genes positively correlated with Ren expression and with mapped quantitative trait loci associated with blood pressure, cardiovascular, and renal phenotypes. The results were reproducible in an independent biological replicate. The findings reported here represent the first map between a critical cardiovascular gene and physical interacting loci throughout the genome and will provide the basis for several new directions of research.

Mechanisms of Mas1 Receptor-Mediated Signaling in the Vascular Endothelium.

OBJECTIVE: Angiotensin II (AngII) has been shown to regulate angiogenesis and at high pathophysiological doses to cause vasoconstriction through the AngII receptor type 1. Angiotensin 1 to 7 (Ang-(1-7)) acting through the Mas1 receptor can act antagonistically to high pathophysiological levels of AngII by inducing vasodilation, whereas the effects of Ang-(1-7) signaling on angiogenesis are less defined. To complicate the matter, there is growing evidence that a subpressor dose of AngII produces phenotypes similar to Ang-(1-7). APPROACH AND RESULTS: This study shows that low-dose Ang-(1-7), acting through the Mas1 receptor, promotes angiogenesis and vasodilation similar to a low, subpressor dose of AngII acting through AngII receptor type 1. In addition, we show through in vitro tube formation that Ang-(1-7) augments the angiogenic response in rat microvascular endothelial cells. Using proteomic and genomic analyses, downstream components of Mas1 receptor signaling were identified, including Rho family of GTPases, phosphatidylinositol 3-kinase, protein kinase D1, mitogen-activated protein kinase, and extracellular signal-related kinase signaling. Further experimental antagonism of extracellular signal-related kinases 1/2 and p38 mitogen-activated protein kinase signaling inhibited endothelial tube formation and vasodilation when stimulated with equimolar, low doses of either AngII or Ang-(1-7). CONCLUSIONS: These results significantly expand the known Ang-(1-7)/Mas1 receptor signaling pathway and demonstrate an important distinction between the pathological effects of elevated and suppressed AngII compared with the beneficial effects of AngII normalization and Ang-(1-7) administration. The observed convergence of Ang-(1-7)/Mas1 and AngII/AngII receptor type 1 signaling at low ligand concentrations suggests a nuanced regulation in vasculature. These data also reinforce the importance of mitogen-activated protein kinase/extracellular signal-related kinase signaling in maintaining vascular function.

Breaking the Cycle: Estrous Variation Does Not Require Increased Sample Size in the Study of Female Rats.

Despite the striking differences between male and female physiology, female physiology is understudied. In response, the National Institutes of Health is promulgating new policies to increase the use of female organisms in preclinical research. Females are commonly believed to have greater variability than males because of the estrous cycle, but recent studies call this belief into question. Effects of estrous cycling on mean arterial pressure were assessed in female Dahl S rats using telemetry and vaginal cytometry and found that estrous cycling did not affect mean arterial pressure magnitude or variance. Data from the PhysGen arm of the Program for Genomic Applications was used to compare male and female variance and coefficient of variation in 142 heart, lung, vascular, kidney, and blood phenotypes, each measured in hundreds to thousands of individual rats from over 50 inbred strains. Seventy-four of 142 phenotypes from this data set demonstrated a sex difference in variance; however, 59% of these phenotypes exhibited greater variance in male rats rather than female. Remarkably, a retrospective power analysis demonstrated that only 16 of 74 differentially variable phenotypes would be detected when using an experimental cohort large enough to detect a difference in magnitude. No overall difference in coefficient of variation between male and female rats was detected when analyzing these 142 phenotypes. We conclude that variability of 142 traits in male and female rats is similar, suggesting that differential treatment of males and females for the purposes of experimental design is unnecessary until proven otherwise, rather than the other way around.

Tumor Necrosis Factor α Regulates Endothelial Progenitor Cell Migration via CADM1 and NF-kB.

Shortly after the discovery of endothelial progenitor cells (EPCs) in 1997, many clinical trials were conducted using EPCs as a cellular based therapy with the goal of restoring damaged organ function by inducing growth of new blood vessels (angiogenesis). Results were disappointing, largely because the cellular and molecular mechanisms of EPC-induced angiogenesis were not clearly understood. Following injection, EPCs must migrate to the target tissue and engraft prior to induction of angiogenesis. In this study EPC migration was investigated in response to tumor necrosis factor α (TNFα), a pro-inflammatory cytokine, to test the hypothesis that organ damage observed in ischemic diseases induces an inflammatory signal that is important for EPC homing. In this study, EPC migration and incorporation were modeled in vitro using a coculture assay where TNFα treated EPCs were tracked while migrating toward vessel-like structures. It was found that TNFα treatment of EPCs increased migration and incorporation into vessel-like structures. Using a combination of genomic and proteomic approaches, NF-kB mediated upregulation of CADM1 was identified as a mechanism of TNFα induced migration. Inhibition of NF-kB or CADM1 significantly decreased migration of EPCs in vitro suggesting a role for TNFα signaling in EPC homing during tissue repair. Stem Cells 2016;34:1922-1933.

Targeting the endothelial progenitor cell surface proteome to identify novel mechanisms that mediate angiogenic efficacy in a rodent model of vascular disease.

Endothelial progenitor cells (EPCs) promote angiogenesis, and clinical trials suggest autologous EPC-based therapy may be effective in treatment of vascular diseases. Albeit promising, variability in the efficacy of EPCs associated with underlying disease states has hindered the realization of EPC-based therapy. Here we first identify and characterize EPC dysfunction in a rodent model of vascular disease (SS/Mcwi rat) that exhibits impaired angiogenesis. To identify molecular candidates that mediate the angiogenic potential of these cells, we performed a broad analysis of cell surface protein expression using chemical labeling combined with mass spectrometry. Analysis revealed EPCs derived from SS/Mcwi rats express significantly more type 2 low-affinity immunoglobulin Fc-gamma (FCGR2) and natural killer 2B4 (CD244) receptors compared with controls. Genome-wide sequencing (RNA-seq) and qt-PCR confirmed isoforms of CD244 and FCGR2a transcripts were increased in SS/Mcwi EPCs. EPCs with elevated expression of FCGR2a and CD244 receptors are predicted to increase the probability of SS/Mcwi EPCs being targeted for death, providing a mechanistic explanation for their reduced angiogenic efficacy in vivo. Pathway analysis supported this contention, as "key" molecules annotated to cell death paths were differentially expressed in the SS/Mcwi EPCs. We speculate that screening and neutralization of cell surface proteins that "tag" and impair EPC function may provide an alternative approach to utilizing incompetent EPCs in greater numbers, as circulating EPCs are depleted in patients with vascular disease. Overall, novel methods to identify putative targets for repair of EPCs using discovery-based technologies will likely provide a major advance in the field of regenerative medicine.

Bone marrow mononuclear cells induce beneficial remodeling and reduce diastolic dysfunction in the left ventricle of hypertensive SS/MCWi rats.

Bone marrow mononuclear cells (BMMNCs) increase capillary density and reduce fibrosis in rodents after myocardial infarction, resulting in an overall improvement in left ventricular function. Little is known about the effectiveness of BMMNC therapy in hypertensive heart disease. In the current study, we show that delivery of BMMNCs from hypertension protected SS-13(BN)/MCWi donor rats, but not BMMNC from hypertension susceptible SS/MCWi donor rats, resulted in 57.2 and 83.4% reductions in perivascular and interstitial fibrosis, respectively, as well as a 60% increase in capillary-to-myocyte count in the left ventricles (LV) of hypertensive SS/MCWi recipients. These histological changes were associated with improvements in LV compliance and relaxation (103 and 46.4% improvements, respectively). Furthermore, improved diastolic function in hypertensive SS/MCWi rats receiving SS-13(BN)/MCWi derived BMMNCs was associated with lower clinical indicators of heart failure, including reductions in end diastolic pressure (65%) and serum brain natriuretic peptide levels (49.9%) with no improvements observed in rats receiving SS/MCWi BMMNCs. SS/MCWi rats had a lower percentage of endothelial progenitor cells in their bone marrow relative to SS-13(BN)/MCWi rats. These results suggest that administration of BMMNCs can prevent or reverse pathological remodeling in hypertensive heart disease, which contributes to ameliorating diastolic dysfunction and associated symptomology. Furthermore, the health and hypertension susceptibility of the BMMNC donor are important factors influencing therapeutic efficacy, possibly via differences in the cellular composition of bone marrow.

Characterization of the genomic structure and function of regions influencing renin and angiogenesis in the SS rat.

Impaired regulation of renin in Dahl salt-sensitive rats (SS/JRHsdMcwi, SS) contributes to attenuated angiogenesis in this strain. This study examined angiogenic function and genomic structure of regions surrounding the renin gene using subcongenic strains of the SS and BN/NHsdMcwi (BN) rat to identify important genomic variations between SS and BN involved in angiogenesis. Three candidate regions on Chr 13 were studied: two congenic strains containing 0.89 and 2.62 Mb portions of BN Chr 13 that excluded the BN renin allele and a third strain that contained a 2.02 Mb overlapping region that included the BN renin allele. Angiogenesis induced by electrical stimulation of the tibialis anterior muscle was attenuated in the SS compared with the BN. Congenics carrying the SS renin allele had impaired angiogenesis, while strains carrying the BN renin allele had angiogenesis restored. The exception was a congenic including a region of BN genome 0.4 Mb distal to renin that restored both renin regulation and angiogenesis. This suggests that there is a distant regulatory element in the BN capable of restoring normal regulation of the SS renin allele. The importance of ANG II in the restored angiogenic response was demonstrated by blocking with losartan. Sequencing of the 4.05 Mb candidate region in SS and BN revealed a total of 8,850 SNPs and other sequence variants. An analysis of the genes and their variants in the region suggested a number of pathways that may explain the impaired regulation of renin and angiogenesis in the SS rat.

Creation and characterization of a renin knockout rat.

The renin-angiotensin system plays an important role in the control of blood pressure (BP) and renal function. To illuminate the importance of renin in the context of a disease background in vivo, we used zinc-finger nucleases (ZFNs) designed to target the renin gene and create a renin knockout in the SS/JrHsdMcwi (SS) rat. ZFN against renin caused a 10-bp deletion in exon 5, resulting in a frameshift mutation. Plasma renin activity was undetectable in the Ren-/- rat, and renin protein was absent from the juxtaglomerular cells in the kidney. Body weight was lower in the Ren-/- rats (than in the Ren+/- or wild-type littermates), and conscious BP on low-salt diet (0.4% NaCl) was 58 ± 2 mm Hg in the Ren-/- male rats versus 117 mm Hg in the Ren+/- littermates, a reduction of almost 50 mm Hg. Blood urea nitrogen (BUN) and plasma creatinine levels were elevated in the Ren-/- strain (BUN 112 ± 7 versus 23 ± 2 mg/dL and creatinine 0.53 ± 0.02 versus 0.26 ± 0.02 mg/dL), and kidney morphology was abnormal with a rudimentary inner renal medulla, cortical interstitial fibrosis, thickening of arterial walls, and abnormally shaped glomeruli. The development of the first rat knockout in the renin-angiotensin system demonstrates the efficacy of the ZFN technology for creating knockout rats for cardiovascular disease on any genetic background and emphasizes the role of renin in BP regulation and kidney function even in the low-renin SS rat.

Role of the renin angiotensin system on bone marrow-derived stem cell function and its impact on skeletal muscle angiogenesis.

Autologous bone marrow cell (BMC) transplantation has been shown as a potential approach to treat various ischemic diseases. However, under many conditions BMC dysfunction has been reported, leading to poor cell engraftment and a failure of tissue revascularization. We have previously shown that skeletal muscle angiogenesis induced by electrical stimulation (ES) is impaired in the SS/Mcwi rats and that this effect is related to a dysregulation of the renin angiotensin system (RAS) that is normalized by the replacement of chromosome 13 derived from the Brown Norway rat (SS-13(BN)/Mcwi consomic rats). The present study explored bone marrow-derived endothelial cell (BM-EC) function in the SS/Mcwi rat and its impact on skeletal muscle angiogenesis induced by ES. SS/Mcwi rats were randomized to receive BMC from: SS/Mcwi; SS-13(BN)/Mcwi; SS/Mcwi rats infused with saline or ANG II (3 ng kg(-1) min(-1)). BMC were injected in the stimulated tibialis anterior muscle of SS/Mcwi rats. Vessel density was evaluated in unstimulated and stimulated muscles after 7 days of ES. BMC isolated from SS/Mcwi or SS/Mcwi rats infused with saline failed to restore angiogenesis induced by ES. However, BMC isolated from SS-13(BN)/Mcwi and SS/Mcwi rats infused with ANG II effectively restored the angiogenesis response in the SS/Mcwi recipient. Furthermore, ANG II infusion increased the capacity of BM-EC to induce endothelial cell tube formation in vitro and slightly increased VEGF protein expression. This study suggests that dysregulation of the RAS in the SS/Mcwi rat contributes to impaired BM-EC function and could impact the angiogenic therapeutic potential of BMC.