Abstract 4165: Novel genetic rat models to identify factors that modulate cardiac and tumor radiation sensitivity

Purpose/Objectives: Over 50% of breast cancer patients receive radiation therapy, but radiation doses can be limited by normal tissue toxicity. Radiation therapy can improve breast cancer-specific survival, but cardiac morbidity can be increased in patients with left-sided tumors. We used rat genetic models to identify targets to improve the therapeutic ratio of radiation. We assessed the radiation responsiveness of mammary tumors and the heart in genetically similar consomic rats conducive to genetic mapping.

Materials/Methods: Female SS rats and SS.BN3 consomic rats, which are genetically identical to SS rats except that chromosome 3 is inherited from the BN strain, have previously been shown to exhibit different vascular dynamics and breast tumor growth. Human MDA-MD-231 cells or syngeneic mammary tumor cells developed from DMBA-induced mammary tumors were implanted orthotopically into immunodeficient or immunocompetent SS and SS.BN3 rats, respectively, and tumors were treated locally with mock or 5x4 Gy. To examine cardiac toxicity, adult female SS and SS.BN3 rats received image-guided localized whole-heart radiation to a dose of 24 Gy or 9 Gy x 5 (AP and 2 lateral fields, weighted 1:1:1). Echocardiograms with strain analysis were performed at baseline, 3 months and 5 months. The Student's t-test was used to compare values.

Results: The BN strain-derived genetic variant(s) on rat chromosome 3 is important for tumor radiation sensitivity. Tumors in SS.BN3 rats were significantly more radiosensitive than tumors in the parental SS strain. A supra-additive effect was seen with both tumor cell lines, with recurrence-free survival of 30% vs. 67% at 137 days in SS vs SS.BN3 rats (p=0.02) in xenografts, and recurrence-free survival of 9% vs. 100% at 141 days in SS vs. SS.BN3 rats (p<0.0001) in syngeneic tumors. The SS female rats that received 24 Gy exhibited enhanced cardiac toxicity compared to SS.BN3 rats, with larger pericardial effusions in SS vs SS.BN3 rats (p<0.05), significantly elevated end diastolic volume (EDV) and end systolic volume (ESV) at 5 months (EDV: 0.62 vs. 0.49 mL, p<0.01; ESV: 0.12 vs. 0.03 mL, p<0.01), and increased cardiac mortality (5/11 SS vs. 0/7 SS.BN3 rats). Fractionated heart radiation yielded similar results. Taken together, the SS.BN3 tumors are more sensitive to radiation, while the hearts of SS.BN3 rats are protected against radiation toxicity, when compared to the SS strain.

Conclusions: These results demonstrate that genetic variants on rat chromosome 3 alter the sensitivity to radiation therapy, enhancing tumor responses to radiation and protecting the heart, thus improving the therapeutic ratio. Gene expression analysis and genetic mapping will be performed to identify the causative target(s). This project has the potential to enhance the effectiveness and toxicity profile of radiation therapy in breast cancer.

Citation Format: Rachel A. Schlaak, Anne Frei, Aronne M. Schottstaedt, Brian L. Fish, Leanne Harmann, Tracy Gasperetti, Michael J. Flister, Meetha Medhora, Jennifer L. Strande, Carmen R. Bergom. Novel genetic rat models to identify factors that modulate cardiac and tumor radiation sensitivity [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4165.

Abstract 408: Evidence of DLL4, NNAT, and SLC35C2 in suppression of breast cancer initiation, growth, and metastasis

Breast cancer affects 1 in 8 women, resulting in 40,000 deaths annually. In most cases, a single cause of breast cancer cannot be found, but rather multiple environmental and genetic factors contribute to overall disease susceptibility. This, combined with complex gene interaction in both malignant tumor cells and nonmalignant tumor microenvironment (TME) cells, poses significant challenges in sifting through the many variants that contribute to the ~31% of breast cancer risk that is heritable. Here, we narrowed the regions associated with breast cancer risk on rat chromosome 3 (RNO3) by introgressing portions of RNO3 derived from the BN rat (protective strain) onto the genomic background of the SS rat (susceptible strain). These SS.BN3 congenics were then phenotyped for DMBA-induced mammary tumor incidence, latency, and multiplicity, which revealed two loci in close proximity that contribute to mammary tumor risk: chr3:95-130Mb (QTL1) and chr3: 154-177Mb (QTL2). By comparing these data with a previous study (Flister et al. Breast Cancer Res Treat. 2017 Aug;165(1):53-64.), we concluded that QTL1 is dependent on the host TME, whereas QTL2 directly modifies breast tumorigenesis and cancer cell proliferation. By combining the congenic mapping studies with genomic and transcriptomic sequencing, and functional analysis, we have now localized the top three candidate modifiers on RNO3: DLL4 (QTL1; TME modifier), NNAT (QTL2; cancer cell modifier), and SLC35C2 (QTL2; cancer cell modifier). Collectively, these data demonstrate the effects of several novel breast cancer modifiers, as well as highlight the potential interactions between modifiers of the malignant cancer cells and the nonmalignant host TME.

Citation Format: Cody Plasterer, Angela Lemke, Dana Murphy, Carmen Bergom, Amit Joshi, Hallgeir Rui, Michael J. Flister. Evidence of DLL4, NNAT, and SLC35C2 in suppression of breast cancer initiation, growth, and metastasis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 408.

Erratum: Methods for detecting host genetic modifiers of tumor vascular function using dynamic near-infrared fluorescence imaging: errata

[This corrects the article on p. 543 in vol. 9, PMID: 29552392.].

mTORC1 and mTORC2 differentially promote natural killer cell development

Natural killer (NK) cells are innate lymphoid cells that are essential for innate and adaptive immunity. Mechanistic target of rapamycin (mTOR) is critical for NK cell development; however, the independent roles of mTORC1 or mTORC2 in regulating this process remain unknown. Ncr1^iCre-mediated deletion of Rptor or Rictor in mice results in altered homeostatic NK cellularity and impaired development at distinct stages. The transition from the CD27⁺CD11b⁻ to the CD27⁺CD11b⁺ stage is impaired in Rptor cKO mice, while, the terminal maturation from the CD27⁺CD11b⁺ to the CD27⁻CD11b⁺ stage is compromised in Rictor cKO mice. Mechanistically, Raptor-deficiency renders substantial alteration of the gene expression profile including transcription factors governing early NK cell development. Comparatively, loss of Rictor causes more restricted transcriptome changes. The reduced expression of T-bet correlates with the terminal maturation defects and results from impaired mTORC2-Akt^S473-FoxO1 signaling. Collectively, our results reveal the divergent roles of mTORC1 and mTORC2 in NK cell development.

Genetic Modifiers of the Breast Tumor Microenvironment

Multiple nonmalignant cell types in the tumor microenvironment (TME) impact breast cancer risk, metastasis, and response to therapy, yet most heritable mechanisms that influence TME cell function and breast cancer outcomes are largely unknown. Breast cancer risk is ∼30% heritable and >170 genetic loci have been associated with breast cancer traits. However, the majority of candidate genes have poorly defined mechanistic roles in breast cancer biology. Research indicates that breast cancer risk modifiers directly impact cancer cells, yet it is equally plausible that some modifier alleles impact the nonmalignant TME. The objective of this review is to examine the list of current breast cancer candidate genes that may modify breast cancer risk and outcome through the TME.

Characterization of Coding/Noncoding Variants for SHROOM3 in Patients with CKD

Background Interpreting genetic variants is one of the greatest challenges impeding analysis of rapidly increasing volumes of genomic data from patients. For example, SHROOM3 is an associated risk gene for CKD, yet causative mechanism(s) of SHROOM3 allele(s) are unknown.Methods We used our analytic pipeline that integrates genetic, computational, biochemical, CRISPR/Cas9 editing, molecular, and physiologic data to characterize coding and noncoding variants to study the human SHROOM3 risk locus for CKD.Results We identified a novel SHROOM3 transcriptional start site, which results in a shorter isoform lacking the PDZ domain and is regulated by a common noncoding sequence variant associated with CKD (rs17319721, allele frequency: 0.35). This variant disrupted allele binding to the transcription factor TCF7L2 in podocyte cell nuclear extracts and altered transcription levels of SHROOM3 in cultured cells, potentially through the loss of repressive looping between rs17319721 and the novel start site. Although common variant mechanisms are of high utility, sequencing is beginning to identify rare variants involved in disease; therefore, we used our biophysical tools to analyze an average of 112,849 individual human genome sequences for rare SHROOM3 missense variants, revealing 35 high-effect variants. The high-effect alleles include a coding variant (P1244L) previously associated with CKD (P=0.01, odds ratio=7.95; 95% CI, 1.53 to 41.46) that we find to be present in East Asian individuals at an allele frequency of 0.0027. We determined that P1244L attenuates the interaction of SHROOM3 with 14-3-3, suggesting alterations to the Hippo pathway, a known mediator of CKD.Conclusions These data demonstrate multiple new SHROOM3-dependent genetic/molecular mechanisms that likely affect CKD.

Abstract P2-11-18: The use of consomic animal models to identify genetic factors that modulate radiation-induced cardiac toxicity

Purpose/Objectives: Radiation therapy is used by more than 50% of breast cancer patients, but radiation doses can be limited by normal tissue side effects. For example, breast cancer radiation therapy can improve breast cancer-specific survival, but increase cardiac deaths in those with left-sided cancers. Identifying genetic factors that can enhance tumor radiation sensitivity while decreasing normal tissue toxicities has the potential to improve the therapeutic ratio of radiation therapy – leading to more cures and less long-term toxicities. The use of animal models with differing genetic backgrounds to assess radiation toxicity, followed by genetic mapping of radiosensitivity phenotypes, has the potential to identify new targets that can predict cardiac toxicity from radiation therapy. This project examines how genetic host factors alter normal tissue toxicity risks from breast cancer radiation.

Materials/Methods: Inbred female SS rats and SS.BN3 consomic rats, that are genetically identical to SS rats except that chromosome 3 is inherited from the BN strain, have previously been shown to exhibit different vascular dynamics and breast tumor growth. For this study, adult female SS and SS.BN3 rats received image-guided whole heart radiation to a dose of 21 Gy (3 fields, AP and 2 laterals). Cardiac troponin was serially measured at 2, 6, and 12 weeks, and echocardiograms with strain analysis were performed at baseline and 3 months. The Student's t-test was used to compare values.

Results: The SS female rats exhibited enhanced cardiac toxicity compared to SS.BN3 rats, with cardiac troponin levels elevated at 12 weeks (0.32 ng/ml vs.0.08 ng/ml for SS vs. SS.BN3, p=0.01), and moderate to severe pericardial effusions seen in 6 of 9 SS rats vs. 2 of 7 SS.BN3 rats. At 3 months post-radiation, echocardiograms revealed increased left ventricular posterior wall thickness at end diastole (LVPWd) in SS vs. SS.BN3 rats (0.25 vs. 0.20 cm, p=0.002) and increased left ventricular mass (LVM) in SS vs. SS.BN3 rats (1.54 vs. 1.28 g, p<0.001). Taken together, the SS female rats are more sensitive to cardiac irradiation than SS.BN3.

Conclusions: These results demonstrate that genetic variant on rat chromosome 3 alter the radiosensitivity to single fraction cardiac radiation therapy. Gene expression analysis and genetic mapping will be performed to identify the causative target(s). These models will also be expanded to test whether similar results are seen with fractionated cardiac radiation therapy. This project has the potential to enhance the effectiveness and toxicity profile of radiation therapy in breast cancer.

Citation Format: Bergom C, Schlaak R, Frei A, Fish BL, Harmann L, Gasperetti T, Schottstaedt AM, Flister MJ, Medhora M, Strande JL. The use of consomic animal models to identify genetic factors that modulate radiation-induced cardiac toxicity [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-11-18.

Methods for detecting host genetic modifiers of tumor vascular function using dynamic near-infrared fluorescence imaging

Vascular supply is a critical component of the tumor microenvironment (TME) and is essential for tumor growth and metastasis, yet the endogenous genetic modifiers that impact vascular function in the TME are largely unknown. To identify the host TME modifiers of tumor vascular function, we combined a novel genetic mapping strategy [Consomic Xenograft Model] with near-infrared (NIR) fluorescence imaging and multiparametric analysis of pharmacokinetic modeling. To detect vascular flow, an intensified cooled camera based dynamic NIR imaging system with 785 nm laser diode based excitation was used to image the whole-body fluorescence emission of intravenously injected indocyanine green dye. Principal component analysis was used to extract the spatial segmentation information for the lungs, liver, and tumor regions-of-interest. Vascular function was then quantified by pK modeling of the imaging data, which revealed significantly altered tissue perfusion and vascular permeability that were caused by host genetic modifiers in the TME. Collectively, these data demonstrate that NIR fluorescent imaging can be used as a non-invasive means for characterizing host TME modifiers of vascular function that have been linked with tumor risk, progression, and response to therapy.

Genetic Fine-Mapping and Identification of Candidate Genes and Variants for Adiposity Traits in Outbred Rats

OBJECTIVE

Obesity is a major risk factor for multiple diseases and is in part heritable, yet the majority of causative genetic variants that drive excessive adiposity remain unknown. Here, outbred heterogeneous stock (HS) rats were used in controlled environmental conditions to fine-map novel genetic modifiers of adiposity.

METHODS

Body weight and visceral fat pad weights were measured in male HS rats that were also genotyped genome-wide. Quantitative trait loci (QTL) were identified by genome-wide association of imputed single-nucleotide polymorphism (SNP) genotypes using a linear mixed effect model that accounts for unequal relatedness between the HS rats. Candidate genes were assessed by protein modeling and mediation analysis of expression for coding and noncoding variants, respectively.

RESULTS

HS rats exhibited large variation in adiposity traits, which were highly heritable and correlated with metabolic health. Fine-mapping of fat pad weight and body weight revealed three QTL and prioritized five candidate genes. Fat pad weight was associated with missense SNPs in Adcy3 and Prlhr and altered expression of Krtcap3 and Slc30a3, whereas Grid2 was identified as a candidate within the body weight locus.

CONCLUSIONS

These data demonstrate the power of HS rats for identification of known and novel heritable mediators of obesity traits.

SmgGDS is a transient nucleolar protein that protects cells from nucleolar stress and promotes the cell cycle by regulating DREAM complex gene expression

The chaperone protein and guanine nucleotide exchange factor SmgGDS (RAP1GDS1) is a key promoter of cancer cell proliferation and tumorigenesis. SmgGDS undergoes nucleocytoplasmic shuttling, suggesting that it has both cytoplasmic and nuclear functions that promote cancer. Previous studies indicate that SmgGDS binds cytoplasmic small GTPases and promotes their trafficking to the plasma membrane. In contrast, little is known about the functions of SmgGDS in the nucleus, or how these nuclear functions might benefit cancer cells. Here we show unique nuclear localization and regulation of gene transcription pathways by SmgGDS. Strikingly, SmgGDS depletion significantly reduces expression of over 600 gene products that are targets of the DREAM complex, which is a transcription factor complex that regulates expression of proteins controlling the cell cycle. The cell cycle regulators E2F1, MYC, MYBL2 (B-Myb) and FOXM1 are among the DREAM targets that are diminished by SmgGDS depletion. E2F1 is well known to promote G1 cell cycle progression, and the loss of E2F1 in SmgGDS-depleted cells provides an explanation for previous reports that SmgGDS depletion characteristically causes a G1 cell cycle arrest. We show that SmgGDS localizes in nucleoli, and that RNAi-mediated depletion of SmgGDS in cancer cells disrupts nucleolar morphology, signifying nucleolar stress. We show that nucleolar SmgGDS interacts with the RNA polymerase I transcription factor upstream binding factor (UBF). The RNAi-mediated depletion of UBF diminishes nucleolar localization of SmgGDS and promotes proteasome-mediated degradation of SmgGDS, indicating that nucleolar sequestration of SmgGDS by UBF stabilizes SmgGDS protein. The ability of SmgGDS to interact with UBF and localize in the nucleolus is diminished by expressing DiRas1 or DiRas2, which are small GTPases that bind SmgGDS and act as tumor suppressors. Taken together, our results support a novel nuclear role for SmgGDS in protecting malignant cells from nucleolar stress, thus promoting cell cycle progression and tumorigenesis.

Host genetic modifiers of nonproductive angiogenesis inhibit breast cancer

PURPOSE

Multiple aspects of the tumor microenvironment (TME) impact breast cancer, yet the genetic modifiers of the TME are largely unknown, including those that modify tumor vascular formation and function.

METHODS

To discover host TME modifiers, we developed a system called the Consomic/Congenic Xenograft Model (CXM). In CXM, human breast cancer cells are orthotopically implanted into genetically engineered consomic xenograft host strains that are derived from two parental strains with different susceptibilities to breast cancer. Because the genetic backgrounds of the xenograft host strains differ, whereas the inoculated tumor cells are the same, any phenotypic variation is due to TME-specific modifier(s) on the substituted chromosome (consomic) or subchromosomal region (congenic). Here, we assessed TME modifiers of growth, angiogenesis, and vascular function of tumors implanted in the SSIL2Rγ and SS.BN3IL2Rγ CXM strains.

RESULTS

Breast cancer xenografts implanted in SS.BN3IL2Rγ (consomic) had significant tumor growth inhibition compared with SSIL2Rγ (parental control), despite a paradoxical increase in the density of blood vessels in the SS.BN3IL2Rγ tumors. We hypothesized that decreased growth of SS.BN3IL2Rγ tumors might be due to nonproductive angiogenesis. To test this possibility, SSIL2Rγ and SS.BN3IL2Rγ tumor vascular function was examined by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), micro-computed tomography (micro-CT), and ex vivo analysis of primary blood endothelial cells, all of which revealed altered vascular function in SS.BN3IL2Rγ tumors compared with SSIL2Rγ. Gene expression analysis also showed a dysregulated vascular signaling network in SS.BN3IL2Rγ tumors, among which DLL4 was differentially expressed and co-localized to a host TME modifier locus (Chr3: 95-131 Mb) that was identified by congenic mapping.

CONCLUSIONS

Collectively, these data suggest that host genetic modifier(s) on RNO3 induce nonproductive angiogenesis that inhibits tumor growth through the DLL4 pathway.

Abstract 348: Interleukin 13 is Required for Neonatal Heart Regeneration

Shortly after birth neonatal mice can fully regenerate their hearts, but this potential is lost in the first week of life. Cell cycle entry of existing cardiomyocytes is thought to be an essential mechanism enabling neonatal mouse heart regeneration. In previous studies we found that the cytokine interleukin 13 (IL13) was a an upstream regulator of differentially expressed gene networks during neonatal heart regeneration and stimulated cell cycle activity of cultured rat cardiomyocytes, suggesting that this factor might be important in neonatal heart regeneration in vivo . In the present study, we subjected wildtype and IL13 knockout mice to ventricular apical resection at one day of age and assessed heart regeneration 21 days post resection (dpr). Compared to wildtype controls, IL13 knockout mice failed to regenerate their hearts as determined by extensive scar formation at the ventricular apex. To gain insight into the mechanism of impaired regeneration, we quantified cardiomyocyte proliferation and expression of macrophage markers at 7 dpr. We found no difference in gene expression of macrophage markers in IL13 knockout mice compared to wildtype. Interestingly, IL13 knockout mice demonstrate a significant increase cardiomyocyte cell cycle activity as determined by phosphorylated Histone H3 (pH3) staining. This seemingly contradictory result appears to be due to an underlying developmental defect in IL13 knockout hearts. Cardiomyocytes in IL13 knockout mice appeared large and disorganized. Cardiomyocytes from IL13 knockout unoperated mice showed decreased pH3 staining and had increased expression marker of hypertrophic growth such as Nppb and Nppa. Histologically, hearts from IL13 knockout mice appeared to have a dilated cardiomyopathy phenotype. Collectively our data suggests that during heart development IL13 influences proliferative versus hypertrophic growth. We surmise that following neonatal apical resection in IL13 knockout mice the significant increase in cardiomyocyte proliferation is a compensatory attempt to repair the injury, but the underlying cardiomyocyte phenotype inhibits complete regeneration. These data are the first to report a role for IL13 in normal heart development and neonatal heart regeneration.

Abstract B07: Utilizing consomic xenograft models to identify genetic variants in the tumor microenvironment that determine breast cancer radiation responses

Progress in elucidating the molecular basis of breast cancer has allowed for treatment breakthroughs such as anti-estrogen and Her2-targeted therapy. It has also shaped the approaches to both surgical and systemic therapy. However, no similar use of molecular information has been utilized to better direct the use of radiation therapy. The development of predictive tools for the radiosensitivity of tumors could allow for personally tailored radiation doses, with treatment de-escalation for radiosensitive tumors, or dose escalation or the use of adjunct treatments in the case of radioresistant tumors. Communication between malignant tumor cells and the tumor microenvironment (TME) underlies most aspects of tumor biology, including chemotherapy and radiation resistance. We have developed a Consomic Xenograft Model (CXM), which maps germline variants that impact only the TME, as well as a species-specific RNA-seq (SSRS) protocol which allows detection of expression changes in the malignant and nonmalignant cellular compartments of tumor xenografts, in parallel and without cell-sorting. Here we utilize these unique techniques to identify genetic variants in the TME that can affect radiation sensitivity. In CXM, human triple negative breast cancer MDA-MD-231 cells are orthotopically implanted into immunodeficient (IL2Rγ-/-) consomic rat strains, which are rat strains in which an entire chromosome is introgressed into the isogenic background of another inbred strain by selective breeding. Because the strain backgrounds are different but the tumor cells are not varied, the observed changes in tumor progression are due to genetic differences in the non-malignant TME. We hypothesized that the tumors in SS.BN3 rats (identical to SS rats but with BN strain chromosome 3) would be more sensitive to radiation due to increased tumor vascularity via CD31 staining, and increased tumor blood volume capacity, as measured by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Our studies demonstrate differential responses to radiation in the CXM model comparing parental SS (IL2Rγ) rats to SS.BN3 (IL2Rγ) rats treated with fractionated radiation therapy (4 Gray x 3), with altered tumor growth kinetics and tumor recurrence rates. A difference was seen in time to 5-fold increase in tumor growth, with 44 vs. >130 days for SS versus SS.BN3 rats (supra-additive, p<0.05). There was a recurrence-free survival of 30% vs. 67% at 130 days, with a median time to recurrence of 57 days vs. time not reached (>130 days) in the SS versus SS.BN3 rats (p=0.02). These results suggest that genetic determinants in the TME affect the radiation sensitivity of genetically identical tumor cells. Using SSRS, we identified a number of candidates on rat chromosome 3 that may potentially influence radiation sensitivity by altering the tumor vasculature. Future studies will further dissect the pathways responsible for the changes in radiation sensitivity. Determining TME factors that affect the radiation sensitivity of tumors has the potential to allow for more tailored and effective radiation treatments in breast cancer.

Citation Format: Carmen Bergom, Michael Straza, Amy Rymaszewski, Anne Frei, Angela Lemke, Shirng-Wern Tsaih, Howard Jacob, Michael J. Flister. Utilizing consomic xenograft models to identify genetic variants in the tumor microenvironment that determine breast cancer radiation responses. [abstract]. In: Proceedings of the AACR Special Conference: Function of Tumor Microenvironment in Cancer Progression; 2016 Jan 7–10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2016;76(15 Suppl):Abstract nr B07.

Abstract PR16: New tools for mapping genetic modifiers of cancer risk in the tumor microenvironment

The majority of heritable breast cancer risk is unknown. One potential source of “missing heritability” is genetic modifiers in the tumor microenvironment (TME). Although genetic modifiers in the TME have long been suspected, they have rarely been studied and are largely unknown. Here, we used two new techniques: the Consomic Xenograft Model (CXM) and species-specific RNA-seq (SSRS) to map genetic modifiers in the TME. In CXM, human breast cancer xenografts are implanted in immunodeficient consomic rat strains and tracked for tumor progression. Because the rat strains vary by one chromosome (i.e., consomic), whereas the malignant tumor cells do not differ, any observed changes in tumor phenotypes are due to genetic modifiers in the TME and can be localized to one chromosome. The SSRS method uses probabilistic mapping of RNAseq reads to a joint human and rat transcriptome to assess differential expression (DE) in malignant (human) tumor cells and the nonmalignant (rat) TME. Validation of SSRS revealed >99.4% specificity in calling human or rat reads, which was significantly better than conventional RNA-seq. Using CXM, we found that BN-derived genetic variant(s) on rat chromosome 3 significantly reduced growth of MDA-MB-231-Luc (231Luc+) tumors by 49% (P<0.05) in the SS.BN3IL2Rγ CXM strain compared with parental SSIL2Rγ. This coincided with a 3.1-fold (P<0.001) decrease in blood vascular invasion by 231Luc+ tumor cells and 7.3-fold (P<0.05) lower metastatic burden in the lungs in SS.BN3IL2Rγ compared with SSIL2Rγ, despite a paradoxical 27% (P<0.05) increase in blood vessel density (BVD) in SS.BN3IL2Rγ rats. The tumor-associated blood vessels in SS.BN3IL2Rγ rats appeared collapsed and dysfunctional, possibly explaining the decreased tumor growth and metastasis, despite increased BVD. Lymphatic vasculature and lymphogenous metastasis were completely unaffected by the SS.BN3IL2Rγ background, suggesting that the causative variant(s) on BN rat chromosome 3 are vascular cell-type specific. We used SSRS to begin identifying the TME-specific mediators on rat chromosome 3 (RNO3) that inhibit growth and hematogenous metastasis of human 231Luc+ breast cancer xenografts implanted in the SS.BN3IL2Rγ. Compared with SSIL2Rγ tumors, we identified a network of 539 DE transcripts in the TME of SS.BN3IL2Rγ rats, of which 28% (150 genes) reside on RNO3, which was significantly higher (>4-fold; P<0.001) than any other rat chromosome. Moreover, a two-sample Kolmogorov-Smirnov test revealed that the difference in distributions of adjusted p-values for RN03 versus the rest of the genome was highly significantly higher for DE genes (P=3.152e-08) or DE transcripts (P=3.441e-16). Compared with other rat chromosomes, RNO3 also had by far the highest incidence of alternative isoform usage (91% of all instances). Pathway analysis of DE genes using DAVID revealed that the two most significant GO clusters were extracellular matrix (49 genes; P<10-20) and blood vessel development (43 genes; P<10-17), which recapitulated the vascular defects observed in the SS.BN3IL2Rγ tumors. Collectively, our data demonstrate that CXM and SSRS can be used to detect genetic modifiers in the TME.

Citation Format: Michael J. Flister, Alexander Stoddard, Shirng-Wern Tsaih, Angela Lemke, Jozef Lazar, Howard Jacob. New tools for mapping genetic modifiers of cancer risk in the tumor microenvironment. [abstract]. In: Proceedings of the AACR Special Conference on Translation of the Cancer Genome; Feb 7-9, 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 1):Abstract nr PR16.

Abstract PR08: New tools for mapping genetic modifiers of cancer risk in the tumor microenvironment

The majority of heritable breast cancer risk is unknown. One potential source of “missing heritability” is genetic modifiers in the tumor microenvironment (TME). Although genetic modifiers in the TME have long been suspected, they have rarely been studied and are largely unknown. Here, we used two new techniques: the Consomic Xenograft Model (CXM) and species-specific RNA-seq (SSRS) to map genetic modifiers in the TME. In CXM, human breast cancer xenografts are implanted in immunodeficient consomic rat strains and tracked for tumor progression. Because the rat strains vary by one chromosome (i.e., consomic), whereas the malignant tumor cells do not differ, any observed changes in tumor phenotypes are due to genetic modifiers in the TME and can be localized to one chromosome. The SSRS method uses probabilistic mapping of RNAseq reads to a joint human and rat transcriptome to assess differential expression (DE) in malignant (human) tumor cells and the nonmalignant (rat) TME. Validation of SSRS revealed >99.4% specificity in calling human or rat reads, which was significantly better than conventional RNA-seq. Using CXM, we found that BN-derived genetic variant(s) on rat chromosome 3 significantly reduced growth of MDA-MB-231-Luc (231Luc+) tumors by 49% (P<0.05) in the SS.BN3IL2Rγ CXM strain compared with parental SSIL2Rγ. This coincided with a 3.1-fold (P<0.001) decrease in blood vascular invasion by 231Luc+ tumor cells and 7.3-fold (P<0.05) lower metastatic burden in the lungs in SS.BN3IL2Rγ compared with SSIL2Rγ, despite a paradoxical 27% (P<0.05) increase in blood vessel density (BVD) in SS.BN3IL2Rγ rats. The tumor-associated blood vessels in SS.BN3IL2Rγ rats appeared collapsed and dysfunctional, possibly explaining the decreased tumor growth and metastasis, despite increased BVD. Lymphatic vasculature and lymphogenous metastasis were completely unaffected by the SS.BN3IL2Rγ background, suggesting that the causative variant(s) on BN rat chromosome 3 are vascular cell-type specific. We used SSRS to begin identifying the TME-specific mediators on rat chromosome 3 (RNO3) that inhibit growth and hematogenous metastasis of human 231Luc+ breast cancer xenografts implanted in the SS.BN3IL2Rγ. Compared with SSIL2Rγ tumors, we identified a network of 539 DE transcripts in the TME of SS.BN3IL2Rγ rats, of which 28% (150 genes) reside on RNO3, which was significantly higher (>4-fold; P<0.001) than any other rat chromosome. Moreover, a two-sample Kolmogorov-Smirnov test revealed that the difference in distributions of adjusted p-values for RN03 versus the rest of the genome was highly significantly higher for DE genes (P=3.152e-08) or DE transcripts (P=3.441e-16). Compared with other rat chromosomes, RNO3 also had by far the highest incidence of alternative isoform usage (91% of all instances). Pathway analysis of DE genes using DAVID revealed that the two most significant GO clusters were extracellular matrix (49 genes; P<10-20) and blood vessel development (43 genes; P<10-17), which recapitulated the vascular defects observed in the SS.BN3IL2Rγ tumors. Collectively, our data demonstrate that CXM and SSRS can be used to detect genetic modifiers in the TME.

Citation Format: Michael J. Flister, Alexander Stoddard, Shirng-Wern Tsaih, Angela Lemke, Jozef Lazar, Howard Jacob. New tools for mapping genetic modifiers of cancer risk in the tumor microenvironment. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr PR08.

2015 Guidelines for Establishing Genetically Modified Rat Models for Cardiovascular Research

The rat has long been a key physiological model for cardiovascular research, most of the inbred strains having been previously selected for susceptibility or resistance to various cardiovascular diseases (CVD). These CVD rat models offer a physiologically relevant background on which candidates of human CVD can be tested in a more clinically translatable experimental setting. However, a diverse toolbox for genetically modifying the rat genome to test molecular mechanisms has only recently become available. Here, we provide a high-level description of several strategies for developing genetically modified rat models of CVD.

MicroRNA-135b Regulates Leucine Zipper Tumor Suppressor 1 in Cutaneous Squamous Cell Carcinoma

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin malignancy and it presents a therapeutic challenge in organ transplant recipient patients. Despite the need, there are only a few targeted drug treatment options. Recent studies have revealed a pivotal role played by microRNAs (miRNAs) in multiple cancers, but only a few studies tested their function in cSCC. Here, we analyzed differential expression of 88 cancer related miRNAs in 43 study participants with cSCC; 32 immunocompetent, 11 OTR patients, and 15 non-lesional skin samples by microarray analysis. Of the examined miRNAs, miR-135b was the most upregulated (13.3-fold, 21.5-fold; p=0.0001) in both patient groups. Similarly, the miR-135b expression was also upregulated in three cSCC cell lines when evaluated by quantitative real-time PCR. In functional studies, inhibition of miR-135b by specific anti-miR oligonucleotides resulted in upregulation of its target gene LZTS1 mRNA and protein levels and led to decreased cell motility and invasion of both primary and metastatic cSCC cell lines. In contrast, miR-135b overexpression by synthetic miR-135b mimic induced further down-regulation of LZTS1 mRNA in vitro and increased cancer cell motility and invasiveness. Immunohistochemical evaluation of 67 cSCC tumor tissues demonstrated that miR-135b expression inversely correlated with LZTS1 staining intensity and the tumor grade. These results indicate that miR-135b functions as an oncogene in cSCC and provide new understanding into its pathological role in cSCC progression and invasiveness.

SH2B3 Is a Genetic Determinant of Cardiac Inflammation and Fibrosis

BACKGROUND

Genome-wide association studies are powerful tools for nominating pathogenic variants, but offer little insight as to how candidate genes affect disease outcome. Such is the case for SH2B adaptor protein 3 (SH2B3), which is a negative regulator of multiple cytokine signaling pathways and is associated with increased risk of myocardial infarction (MI), but its role in post-MI inflammation and fibrosis is completely unknown.

METHODS AND RESULTS

Using an experimental model of MI (left anterior descending artery occlusion/reperfusion injury) in wild-type and Sh2b3 knockout rats (Sh2b3(em2Mcwi)), we assessed the role of Sh2b3 in post-MI fibrosis, leukocyte infiltration, angiogenesis, left ventricle contractility, and inflammatory gene expression. Compared with wild-type, Sh2b3(em2Mcwi) rats had significantly increased fibrosis (2.2-fold; P<0.05) and elevated leukocyte infiltration (>2-fold; P<0.05), which coincided with decreased left ventricle fractional shortening (-Δ11%; P<0.05) at 7 days post left anterior descending artery occlusion/reperfusion injury. Despite an increased angiogenic potential in Sh2b3(em2Mcwi) rats (1.7-fold; P<0.05), we observed no significant differences in left ventricle capillary density between wild-type and Sh2b3(em2Mcwi) rats. In total, 12 genes were significantly elevated in the post left anterior descending artery occluded/reperfused hearts of Sh2b3(em2Mcwi) rats relative to wild-type, of which 3 (NLRP12, CCR2, and IFNγ) were significantly elevated in the left ventricle of heart failure patients carrying the MI-associated rs3184504 [T] SH2B3 risk allele.

CONCLUSIONS

These data demonstrate for the first time that SH2B3 is a crucial mediator of post-MI inflammation and fibrosis.

Shroom3 contributes to the maintenance of the glomerular filtration barrier integrity

Genome-wide association studies (GWAS) identify regions of the genome correlated with disease risk but are restricted in their ability to identify the underlying causative mechanism(s). Thus, GWAS are useful "roadmaps" that require functional analysis to establish the genetic and mechanistic structure of a particular locus. Unfortunately, direct functional testing in humans is limited, demonstrating the need for complementary approaches. Here we used an integrated approach combining zebrafish, rat, and human data to interrogate the function of an established GWAS locus (SHROOM3) lacking prior functional support for chronic kidney disease (CKD). Congenic mapping and sequence analysis in rats suggested Shroom3 was a strong positional candidate gene. Transferring a 6.1-Mb region containing the wild-type Shroom3 gene significantly improved the kidney glomerular function in FHH (fawn-hooded hypertensive) rat. The wild-type Shroom3 allele, but not the FHH Shroom3 allele, rescued glomerular defects induced by knockdown of endogenous shroom3 in zebrafish, suggesting that the FHH Shroom3 allele is defective and likely contributes to renal injury in the FHH rat. We also show for the first time that variants disrupting the actin-binding domain of SHROOM3 may cause podocyte effacement and impairment of the glomerular filtration barrier.

Abstract 025: Disruption of Shroom3 contributes to impaired glomerular podocyte function

SHROOM3 has been associated with chronic kidney disease by genome-wide association studies (GWAS), yet its role in renal function and diseases was completely unknown. Here, we leveraged the integrated physiological genomic approach in rat models and zebrafish to dissect the renal function of SHROOM3 and its contribution to renal impairments. Congenic mapping and sequence analysis in rats suggested Shroom3 was a strong positional candidate gene. Transferring a 14.50-21.40 Mbp region of the BN (Brown Norway) chromosome 14 containing the Shroom3 gene onto the FHH (Fawn-Hooded Hypertensive) background significantly lowered the level of albuminuria (31.07±3.98 vs 51.92±5.69 mg/day; P=0.001), glomerular permeability to albumin (0.48±0.02 vs 0.59±0.02; P<0.001), and glomerular sclerosis (0.72±0.03 vs 1.82±0.05; P<0.001). The FHH Shroom3 allele contains 13 amino acid variants compared to BN rat. In vivo injection of the wild-type BN Shroom3 mRNA, but not the FHH Shroom3, rescued the glomerular leakage induced by endogenous knockdown of shroom3 in zebrafish, suggesting that the FHH Shroom3 allele is defective and likely contributes to glomerular impairments seen in the FHH rat. Functional screening of the 13 FHH Shroom3 variants identified a mutation, which decreased Shroom3 binding affinity to actins (5.85±2.69 vs 52.28±16.18; P<0.05) and consequently, impaired glomerular permeability. These results demonstrate that Shroom3-mediated interaction with actin is a crucial mechanism underlying the control of glomerular filtration barrier. Furthermore, we showed that podocyte-specific shroom3 disruption in zebrafish caused glomerular impairments and podocyte effacement, indicated by reduced foot process number (2.85±0.07 vs 4.2±0.35; P<0.001) and increased foot process size (319.85±19.42 vs 164.39±8.24 nm; P<0.001). Taken together, these experiments provide the first direct demonstration that Shroom3 is functionally involved in renal pathophysiology by the control of glomerular podocyte function.

CXM: a new tool for mapping breast cancer risk in the tumor microenvironment

The majority of causative variants in familial breast cancer remain unknown. Of the known risk variants, most are tumor cell autonomous, and little attention has been paid yet to germline variants that may affect the tumor microenvironment. In this study, we developed a system called the Consomic Xenograft Model (CXM) to map germline variants that affect only the tumor microenvironment. In CXM, human breast cancer cells are orthotopically implanted into immunodeficient consomic strains and tumor metrics are quantified (e.g., growth, vasculogenesis, and metastasis). Because the strain backgrounds vary, whereas the malignant tumor cells do not, any observed changes in tumor progression are due to genetic differences in the nonmalignant microenvironment. Using CXM, we defined genetic variants on rat chromosome 3 that reduced relative tumor growth and hematogenous metastasis in the SS.BN3(IL2Rγ) consomic model compared with the SS(IL2Rγ) parental strain. Paradoxically, these effects occurred despite an increase in the density of tumor-associated blood vessels. In contrast, lymphatic vasculature and lymphogenous metastasis were unaffected by the SS.BN3(IL2Rγ) background. Through comparative mapping and whole-genome sequence analysis, we narrowed candidate variants on rat chromosome 3 to six genes with a priority for future analysis. Collectively, our results establish the utility of CXM to localize genetic variants affecting the tumor microenvironment that underlie differences in breast cancer risk.

Vascular dysfunction precedes hypertension associated with a blood pressure locus on rat chromosome 12

We previously isolated a 6.1-Mb region of SS/Mcwi (Dahl salt-sensitive) rat chromosome 12 (13.4-19.5 Mb) that significantly elevated blood pressure (BP) (Δ+34 mmHg, P < 0.001) compared with the SS-12(BN) consomic control. In the present study, we examined the role of vascular dysfunction and remodeling in hypertension risk associated with the 6.1-Mb (13.4-19.5 Mb) locus on rat chromosome 12 by reducing dietary salt, which lowered BP levels so that there were no substantial differences in BP between strains. Consequently, any observed differences in the vasculature were considered BP-independent. We also reduced the candidate region from 6.1 Mb with 133 genes to 2 Mb with 23 genes by congenic mapping. Both the 2 Mb and 6.1 Mb congenic intervals were associated with hypercontractility and decreased elasticity of resistance vasculature prior to elevations of BP, suggesting that the vascular remodeling and dysfunction likely contribute to the pathogenesis of hypertension in these congenic models. Of the 23 genes within the narrowed congenic interval, 12 were differentially expressed between the resistance vasculature of the 2 Mb congenic and SS-12(BN) consomic strains. Among these, Grifin was consistently upregulated 2.7 ± 0.6-fold (P < 0.05) and 2.0 ± 0.3-fold (P < 0.01), and Chst12 was consistently downregulated -2.8 ± 0.3-fold (P < 0.01) and -4.4 ± 0.4-fold (P < 0.00001) in the 2 Mb congenic compared with SS-12(BN) consomic under normotensive and hypertensive conditions, respectively. A syntenic region on human chromosome 7 has also been associated with BP regulation, suggesting that identification of the genetic mechanism(s) underlying cardiovascular phenotypes in this congenic strain will likely be translated to a better understanding of human hypertension.