Genetic architecture of mouse skin inflammation and tumour susceptibility

Stem-cell states converge in multistage cutaneous squamous cell carcinoma development

Stem cells play a critical role in cancer development by contributing to cell heterogeneity, lineage plasticity, and drug resistance. We created gene expression networks from hundreds of mouse tissue samples (both normal and tumor) and integrated these with lineage tracing and single-cell RNA-seq, to identify convergence of cell states in premalignant tumor cells expressing markers of lineage plasticity and drug resistance. Two of these cell states representing multilineage plasticity or proliferation were inversely correlated, suggesting a mutually exclusive relationship. Treatment of carcinomas in vivo with chemotherapy repressed the proliferative state and activated multilineage plasticity whereas inhibition of differentiation repressed plasticity and potentiated responses to cell cycle inhibitors. Manipulation of this cell state transition point may provide a source of potential combinatorial targets for cancer therapy.

NCAPH drives breast cancer progression and identifies a gene signature that predicts luminal a tumour recurrence

BACKGROUND

Luminal A tumours generally have a favourable prognosis but possess the highest 10-year recurrence risk among breast cancers. Additionally, a quarter of the recurrence cases occur within 5 years post-diagnosis. Identifying such patients is crucial as long-term relapsers could benefit from extended hormone therapy, while early relapsers might require more aggressive treatment.

METHODS

We conducted a study to explore non-structural chromosome maintenance condensin I complex subunit H's (NCAPH) role in luminal A breast cancer pathogenesis, both in vitro and in vivo, aiming to identify an intratumoural gene expression signature, with a focus on elevated NCAPH levels, as a potential marker for unfavourable progression. Our analysis included transgenic mouse models overexpressing NCAPH and a genetically diverse mouse cohort generated by backcrossing. A least absolute shrinkage and selection operator (LASSO) multivariate regression analysis was performed on transcripts associated with elevated intratumoural NCAPH levels.

RESULTS

We found that NCAPH contributes to adverse luminal A breast cancer progression. The intratumoural gene expression signature associated with elevated NCAPH levels emerged as a potential risk identifier. Transgenic mice overexpressing NCAPH developed breast tumours with extended latency, and in Mouse Mammary Tumor Virus (MMTV)-NCAPHErbB2 double-transgenic mice, luminal tumours showed increased aggressiveness. High intratumoural Ncaph levels correlated with worse breast cancer outcome and subpar chemotherapy response. A 10-gene risk score, termed Gene Signature for Luminal A 10 (GSLA10), was derived from the LASSO analysis, correlating with adverse luminal A breast cancer progression.

CONCLUSIONS

The GSLA10 signature outperformed the Oncotype DX signature in discerning tumours with unfavourable outcomes, previously categorised as luminal A by Prediction Analysis of Microarray 50 (PAM50) across three independent human cohorts. This new signature holds promise for identifying luminal A tumour patients with adverse prognosis, aiding in the development of personalised treatment strategies to significantly improve patient outcomes.

NCAPH Drives Breast Cancer Progression and Identifies a Gene Signature that Predicts Luminal A Tumor Recurrence

Despite their generally favorable prognosis, luminal A tumors paradoxically pose the highest ten-year recurrence risk among breast cancers. From those that relapse, a quarter of them do it within five years after diagnosis. Identifying such patients is crucial, as long-term relapsers could benefit from extended hormone therapy, whereas early relapsers may require aggressive treatment. In this study, we demonstrate that NCAPH plays a role in the pathogenesis of luminal A breast cancer, contributing to its adverse progression in vitro and in vivo. Furthermore, we reveal that a signature of intratumoral gene expression, associated with elevated levels of NCAPH, serves as a potential marker to identify patients facing unfavorable progression of luminal A breast cancer. Indeed, transgenic mice overexpressing NCAPH generated breast tumors with long latency, and in MMTV-NCAPH/ErbB2+ double-transgenic mice, the luminal tumors formed were more aggressive. In addition, high intratumoral levels of Ncaph were associated with worse breast cancer evolution and poor response to chemotherapy in a cohort of genetically heterogeneous transgenic mice generated by backcrossing. In this cohort of mice, we identified a series of transcripts associated with elevated intratumoral levels of NCAPH, which were linked to adverse progression of breast cancer in both mice and humans. Utilizing the Least Absolute Shrinkage and Selection Operator (LASSO) multivariate regression analysis on this series of transcripts, we derived a ten-gene risk score. This score is defined by a gene signature (termed Gene Signature for Luminal A 10 or GSLA10) that correlates with unfavorable progression of luminal A breast cancer. The GSLA10 signature surpassed the Oncotype DX signature in discerning tumors with unfavorable outcomes (previously categorized as Luminal A by PAM50) across three independent human cohorts. This GSLA10 signature aids in identifying patients with Luminal A tumors displaying adverse prognosis, who could potentially benefit from personalized treatment strategies.

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.

Pan-cancer evaluation of clinical value of mitotic network activity index (MNAI) and its predictive value for immunotherapy

Increased mitotic activity is associated with the genesis and aggressiveness of many cancers. To assess the clinical value of mitotic activity as prognostic biomarker, we performed a pan-cancer study on the mitotic network activity index (MNAI) constructed based on 54-gene mitotic apparatus network. Our pan-cancer assessment on TCGA (33 tumor types, 10,061 patients) and validation on other publicly available cohorts (23 tumor types, 9,209 patients) confirmed the significant association of MNAI with overall survival, progression-free survival, and other prognostic endpoints in multiple cancer types, including lower-grade gliomas (LGG), breast invasive carcinoma (BRCA), as well as many others. We also showed significant association between MNAI and genetic instability, which provides a biological explanation of its prognostic impact at pan-cancer landscape. Our association analysis revealed that patients with high MNAI benefitted more from anti-PD-1 and Anti-CTLA-4 treatment. In addition, we demonstrated that multimodal integration of MNAI and the AI-empowered Cellular Morphometric Subtypes (CMS) significantly improved the predictive power of prognosis compared to using MNAI and CMS alone. Our results suggest that MNAI can be used as a potential prognostic biomarker for different tumor types toward different clinical endpoints, and multimodal integration of MNAI and CMS exceeds individual biomarker for precision prognosis.

Gene networks reveal stem-cell state convergence during preneoplasia and progression to malignancy in multistage skin carcinogenesis

Adult mammalian stem cells play critical roles in normal tissue homeostasis, as well as in tumor development, by contributing to cell heterogeneity, plasticity, and development of drug resistance. The relationship between different types of normal and cancer stem cells is highly controversial and poorly understood. Here, we carried out gene expression network analysis of normal and tumor samples from genetically heterogeneous mice to create network metagenes for visualization of stem-cell networks, rather than individual stem-cell markers, at the single-cell level during multistage carcinogenesis. We combined this approach with lineage tracing and single-cell RNASeq of stem cells and their progeny, identifying a previously unrecognized hierarchy in which Lgr6+ stem cells from tumors generate progeny that express a range of other stem-cell markers including Sox2, Pitx1, Foxa1, Klf5, and Cd44. Our data identify a convergence of multiple stem-cell and tumor-suppressor pathways in benign tumor cells expressing markers of lineage plasticity and oxidative stress. This same single-cell population expresses network metagenes corresponding to markers of cancer drug resistance in human tumors of the skin, lung and prostate. Treatment of mouse squamous carcinomas in vivo with the chemotherapeutic cis-platin resulted in elevated expression of the genes that mark this cell population. Our data have allowed us to create a simplified model of multistage carcinogenesis that identifies distinct stem-cell states at different stages of tumor progression, thereby identifying networks involved in lineage plasticity, drug resistance, and immune surveillance, providing a rich source of potential targets for cancer therapy.

Expanding role of CXCR2 and therapeutic potential of CXCR2 antagonists in inflammatory diseases and cancers

C-X-C motif chemokine receptor 2 (CXCR2) is G protein-coupled receptor (GPCR) and plays important roles in various inflammatory diseases and cancers, including chronic obstructive pulmonary disease (COPD), atherosclerosis, asthma, and pancreatic cancer. Upregulation of CXCR2 is closely associated with the migration of neutrophils and monocytes. To date, many small-molecule CXCR2 antagonists have entered clinical trials, showing favorable safety and therapeutic effects. Hence, we provide an overview containing the discovery history, protein structure, signaling pathways, biological functions, structure-activity relationships and clinical significance of CXCR2 antagonists in inflammatory diseases and cancers. According to the latest development and recent clinical progress of CXCR2 small molecule antagonists, we speculated that CXCR2 can be used as a biomarker and a new target for diabetes and that CXCR2 antagonists may also attenuate lung injury in coronavirus disease 2019 (COVID-19).

Vitamin D Receptor Promotes Global Nucleotide Excision Repair by Facilitating XPC Dissociation from Damaged DNA

Vitamin D receptor (VDR) is important for normal DNA repair, although the mechanism by which it acts is unclear. After focal UV irradiation to create subnuclear spots of DNA damage, epidermal keratinocytes from VDR-null mice as well as human epidermal keratinocytes depleted of VDR with small interfering RNA removed pyrimidine-pyrimidone (6-4) photoproducts more slowly than control cells. Costaining with antibodies to XPC, the DNA damage recognition sensor that initiates nucleotide excision repair, showed that XPC rapidly accumulated at spots of damage and gradually faded in control human keratinocytes. In VDR-depleted keratinocytes, XPC associated with DNA damage with comparable efficiency; however, XPC's dissociation dynamics were altered so that significantly more XPC was bound and retained over time than in control cells. The XPF endonuclease, which acts subsequently in nucleotide excision repair, bound and dissociated with comparable kinetics in control and VDR-depleted cells, but the extent of binding was reduced in the latter. These results as well as kinetic modeling of the data suggest that VDR's importance in the repair of UV-induced DNA damage is mediated in part by its ability to facilitate the dissociation of XPC from damaged DNA for the normal recruitment and assembly of other repair proteins to proceed.

Notch dimerization and gene dosage are important for normal heart development, intestinal stem cell maintenance, and splenic marginal zone B-cell homeostasis during mite infestation

Cooperative DNA binding is a key feature of transcriptional regulation. Here we examined the role of cooperativity in Notch signaling by CRISPR-mediated engineering of mice in which neither Notch1 nor Notch2 can homo- or heterodimerize, essential for cooperative binding to sequence-paired sites (SPS) located near many Notch-regulated genes. Although most known Notch-dependent phenotypes were unaffected in Notch1/2 dimer-deficient mice, a subset of tissues proved highly sensitive to loss of cooperativity. These phenotypes include heart development, compromised viability in combination with low gene dose, and the gut, developing ulcerative colitis in response to 1% dextran sulfate sodium (DSS). The most striking phenotypes-gender imbalance and splenic marginal zone B-cell lymphoma-emerged in combination with gene dose reduction or when challenged by chronic fur mite infestation. This study highlights the role of the environment in malignancy and colitis and is consistent with Notch-dependent anti-parasite immune responses being compromised in Notch dimer-deficient animals.

Pathophysiology of pachyonychia congenita-associated palmoplantar keratoderma: new insights into skin epithelial homeostasis and avenues for treatment

BACKGROUND

Pachyonychia congenita (PC), a rare genodermatosis, primarily affects ectoderm-derived epithelial appendages and typically includes oral leukokeratosis, nail dystrophy and very painful palmoplantar keratoderma (PPK). PC dramatically impacts quality of life although it does not affect lifespan. PC can arise from mutations in any of the wound-repair-associated keratin genes KRT6A, KRT6B, KRT6C, KRT16 or KRT17. There is no cure for this condition, and current treatment options for PC symptoms are limited and palliative in nature.

OBJECTIVES

This review focuses on recent progress made towards understanding the pathophysiology of PPK lesions, the most prevalent and debilitating of all PC symptoms.

METHODS

We reviewed the relevant literature with a particular focus on the Krt16 null mouse, which spontaneously develops footpad lesions that mimic several aspects of PC-associated PPK.

RESULTS

There are three main stages of progression of PPK-like lesions in Krt16 null mice. Ahead of lesion onset, keratinocytes in the palmoplantar (footpad) skin exhibit specific defects in terminal differentiation, including loss of Krt9 expression. At the time of PPK onset, there is elevated oxidative stress and hypoactive Keap1-Nrf2 signalling. During active PPK, there is a profound defect in the ability of the epidermis to maintain or return to normal homeostasis.

CONCLUSIONS

The progress made suggests new avenues to explore for the treatment of PC-based PPK and deepens our understanding of the mechanisms controlling skin tissue homeostasis. What's already known about this topic? Pachyonychia congenita (PC) is a rare genodermatosis caused by mutations in KRT6A, KRT6B, KRT6C, KRT16 and KRT17, which are normally expressed in skin appendages and induced following injury. Individuals with PC present with multiple clinical symptoms that usually include thickened and dystrophic nails, palmoplantar keratoderma (PPK), glandular cysts and oral leukokeratosis. The study of PC pathophysiology is made challenging because of its low incidence and high complexity. There is no cure or effective treatment for PC. What does this study add? This text reviews recent progress made when studying the pathophysiology of PPK associated with PC. This recent progress points to new possibilities for devising effective therapeutics that may complement current palliative strategies.

The Vitamin D Receptor as Tumor Suppressor in Skin

Cutaneous malignancies including melanomas and keratinocyte carcinomas (KC) are the most common types of cancer, occurring at a rate of over one million per year in the United States. KC, which include both basal cell carcinomas and squamous cell carcinomas, are substantially more common than melanomas and form the subject of this chapter. Ultraviolet radiation (UVR), both UVB and UVA, as occurs with sunlight exposure is generally regarded as causal for these malignancies, but UVB is also required for vitamin D synthesis in the skin. Keratinocytes are the major cell in the epidermis. These cells not only produce vitamin D but contain the enzymatic machinery to metabolize vitamin D to its active metabolite, 1,25(OH)₂D, and express the receptor for this metabolite, the vitamin D receptor (VDR). This allows the cell to respond to the 1,25(OH)₂D that it produces. Based on our own data and that reported in the literature, we conclude that vitamin D signaling in the skin suppresses UVR-induced epidermal tumor formation. In this chapter we focus on four mechanisms by which vitamin D signaling suppresses tumor formation. They are inhibition of proliferation/stimulation of differentiation with discussion of the roles of hedgehog, Wnt/β-catenin, and hyaluronan/CD44 pathways in mediating vitamin D regulation of proliferation/differentiation, regulation of the balance between oncogenic and tumor suppressor long noncoding RNAs, immune regulation, and promotion of DNA damage repair (DDR).

Niche-Specific Factors Dynamically Regulate Sebaceous Gland Stem Cells in the Skin

Oil-secreting sebaceous glands (SGs) are critical for proper skin function; however, it remains unclear how different factors act together to modulate SG stem cells. Here, we provide functional evidence that each SG lobe is serviced by its own dedicated stem cell population. Upon ablating Notch signaling in different skin subcompartments, we find that this pathway exerts dual counteracting effects on SGs. Suppressing Notch in SG progenitors traps them in a hybrid state where stem and differentiation features become intermingled. In contrast, ablating Notch outside of the SG stem cell compartment indirectly drives SG expansion. Finally, we report that a K14:K5→K14:K79 keratin shift occurs during SG differentiation. Deleting K79 destabilizes K14 in sebocytes, and attenuates SGs and eyelid meibomian glands, leading to corneal ulceration. Altogether, our findings demonstrate that SGs integrate diverse signals from different niches and suggest that mutations incurred within one stem cell compartment can indirectly influence another.

Vitamin D receptor expression and serum 25(OH)D concentration inversely associates with burden of neurofibromas

Vitamin D and its receptor may play a role in preventing tumor development and progression. As such antineoplastic effects are expected to be weak and to act over long periods, conditions with increased tumor incidence, such as the neurofibromatosis type 1 (NF1), provide suitable study models. We previously found an inverse correlation of serum 25(OH)D concentration with number of neurofibromas in NF1. Here we aim to further explore the role of the vitamin D receptor. A total of 141 adult NF1 patients were included in the study. For 101 of them, serum vitamin 25(OH)D data were available. From 87 patients, blood samples were obtained in PaxGene tubes containing a reagent to stabilize RNA immediately. mRNA of the vitamin D receptor (VDR) gene (coding for the vitamin D receptor) was measured by means of RT-PCR. Correlation of laboratory data with NF1-related tumors was statistically evaluated. Vitamin D receptor in NF1-tumors was examined by means of immunohistochemistry using an antibody against the vitamin D1 receptor. The number of dermal neurofibromas was significantly inversely correlated with VDR mRNA level and with serum 25(OH)D concentration in NF1 patients. In contrast, plexiform neurofibroma and malignant peripheral nerve sheath tumor did not correlate with these two parameters. Immunostaining did not detect vitamin D receptor in NF1-tumors. Both vitamin D and its receptor may play a role in suppressing the development of neurofibromas. Sustaining 25(OH)D at an adequate level may contribute to controlling neurofibromas and possibly also other tumors. This is especially important for individuals with lower expression of VDR.

Transcriptome signatures of p,p´-DDE-induced liver damage in Mus spretus mice

The use of DDT (1,1,1-trichloro-2,2-bis(p-chlorophenyl) ethane) in some countries, although regulated, is contributing to an increased worldwide risk of exposure to this organochlorine pesticide or its derivative p,p'-DDE [1,1-dichloro-2,2-bis(p-chlorophenyl) ethylene]. Many studies have associated p,p'-DDE exposure to type 2 diabetes, obesity and alterations of the reproductive system, but their molecular mechanisms of toxicity remain poorly understood. We have addressed this issue by using commercial microarrays based on probes for the entire Mus musculus genome to determine the hepatic transcriptional signatures of p,p'-DDE in the phylogenetically close mouse species Mus spretus. High-stringency hybridization conditions and analysis assured reliable results, which were also verified, in part, by qRT-PCR, immunoblotting and/or enzymatic activity. Our data linked 198 deregulated genes to mitochondrial dysfunction and perturbations of central signaling pathways (kinases, lipids, and retinoic acid) leading to enhanced lipogenesis and aerobic glycolysis, inflammation, cell proliferation and testosterone catabolism and excretion. Alterations of transcript levels of genes encoding enzymes involved in testosterone catabolism and excretion would explain the relationships established between p,p´-DDE exposure and reproductive disorders, obesity and diabetes. Further studies will help to fully understand the molecular basis of p,p´-DDE molecular toxicity in liver and reproductive organs, to identify effective exposure biomarkers and perhaps to design efficient p,p'-DDE exposure counteractive strategies.

Analysis of chromatin accessibility in human epidermis identifies putative barrier dysfunction-sensing enhancers

To identify putative gene regulatory regions that respond to epidermal injury, we mapped chromatin dynamics in a stratified human epidermis during barrier maturation and disruption. Engineered skin substitutes (ESS) cultured at the air-liquid interface were used as a model of developing human epidermis with incomplete barrier formation. The epidermal barrier stabilized following engraftment onto immunocompromised mice, and was compromised again upon injury. Modified formaldehyde-assisted isolation of regulatory elements (FAIRE) was used to identify accessible genomic regions characteristic of monolayer keratinocytes, ESS in vitro, grafted ESS, and tape-stripped ESS graft. We mapped differentiation- and maturation-associated changes in transcription factor binding sites enriched at each stage and observed overrepresentation of AP-1 gene family motifs in barrier-deficient samples. Transcription of TSLP, an important effector of immunological memory in response to allergen exposure, was dramatically elevated in our barrier-deficient samples. We identified dynamic DNA elements that correlated with TSLP induction and may contain enhancers that regulate TSLP. Two dynamic regions were located near the TSLP promoter and overlapped with allergy-associated SNPs rs17551370 and rs2289877, strongly implicating these loci in the regulation of TSLP expression in allergic disease. Additional dynamic chromatin regions ~250kb upstream of the TSLP promoter were found to be in high linkage disequilibrium with allergic disease SNPs. Taken together, these results define dynamic chromatin accessibility changes during epidermal development and dysfunction.

Lgr6 is a stem cell marker in mouse skin squamous cell carcinoma

The G-protein-coupled receptors Lgr4/5/6 are Wnt signalling mediators, but their functions in squamous carcinomas (SCCs) are unclear. Using lineage tracing in Lgr5-EGFP-CreERT2- and Lgr6-EGFP-CreERT2- Rosa26/Tomato reporter mice, we demonstrate that Lgr6, but not Lgr5, acts as an epithelial stem cell marker in vivo in SCCs. We identify, by single molecule in situ hybridisation and cell sorting, rare Lgr6-positive cells in immortalised keratinocytes, and show that their frequency increases in advanced SCCs. Lgr6 expression is enriched in cells with stem cell characteristics, and Lgr6 downregulation in vivo causes increased epidermal proliferation, with expanded lineage tracing from Lgr6⁺ epidermal stem cells. Surprisingly, Lgr6 germline knockout mice are predisposed to SCC development, by a mechanism that includes compensatory upregulation of Lgr5. These data provide a model for human patients with germline loss of function mutations in WNT pathway genes RSPO1 or LGR4, who show increased susceptibility to squamous tumour development.

Molecular and genetic inflammation networks in major human diseases

It has been well-recognized that inflammation alongside tissue repair and damage maintaining tissue homeostasis determines the initiation and progression of complex diseases. Albeit with the accomplishment of having captured the most critical inflammation-involved molecules, genetic susceptibilities, epigenetic factors, and environmental factors, our schemata on the role of inflammation in complex diseases remain largely patchy, in part due to the success of reductionism in terms of research methodology per se. Omics data alongside the advances in data integration technologies have enabled reconstruction of molecular and genetic inflammation networks which shed light on the underlying pathophysiology of complex diseases or clinical conditions. Given the proven beneficial role of anti-inflammation in coronary heart disease as well as other complex diseases and immunotherapy as a revolutionary transition in oncology, it becomes timely to review our current understanding of the molecular and genetic inflammation networks underlying major human diseases. In this review, we first briefly discuss the complexity of infectious diseases and then highlight recently uncovered molecular and genetic inflammation networks in other major human diseases including obesity, type II diabetes, coronary heart disease, late onset Alzheimer's disease, Parkinson's disease, and sporadic cancer. The commonality and specificity of these molecular networks are addressed in the context of genetics based on genome-wide association study (GWAS). The double-sword role of inflammation, such as how the aberrant type 1 and/or type 2 immunity leads to chronic and severe clinical conditions, remains open in terms of the inflammasome and the core inflammatome network features. Increasingly available large Omics and clinical data in tandem with systems biology approaches have offered an exciting yet challenging opportunity toward reconstruction of more comprehensive and dynamic molecular and genetic inflammation networks, which hold great promise in transiting network snapshots to video-style multi-scale interplays of disease mechanisms, in turn leading to effective clinical intervention.

Chemical Carcinogenesis Models of Cancer: Back to the Future

Over a century has elapsed since the first demonstration that exposure to chemicals in coal tar can cause cancer in animals. These observations provided an essential causal mechanistic link between environmental chemicals and increased risk of cancer in human populations. Mouse models of chemical carcinogenesis have since led to the concept of multistage tumor development through distinct stages of initiation, promotion, and progression and identified many of the genetic and biological events involved in these processes. Recent breakthroughs in DNA sequencing have now given us tools to dissect complete tumor genome architectures and revealed that chemically induced cancers in the mouse carry a high point mutation load and mutation signatures that reflect the causative agent used for tumor induction. Chemical carcinogenesis models may therefore provide a route to identify the causes of mutation signatures found in human cancers and further inform studies of therapeutic drug resistance and responses to immunotherapy, which are dependent on mutation load and genetic heterogeneity.

Cichorium intybus root extract: A "vitamin D-like" active ingredient to improve skin barrier function

During the aging process, the human skin suffers many alterations including dryness, skin barrier function damage. The skin barrier function is important to the prevention of skin alterations and maintenance of homeostasis. So, the objective of this study was to assess the clinical efficacy on skin barrier function of Cichorium intybus root extract in cosmetic formulations with or without UV filters. Fifty women, aged between 45 and 60 years, were divided into two groups. One group received vehicle formulations containing UV filters, and the other group received formulations without UV filters. Both groups received a formulation containing the extract and the vehicle. The formulations were applied twice daily to the upper arms after washing with sodium lauryl sulphate. Transepidermal water loss (TEWL) and skin microrelief were evaluated before and after a 14- and 28-day period of treatment. The control regions and regions where the vehicles were applied showed an increase in the TEWL. For the formulations containing the extract, decreased TEWL and improved microrelief were observed when compared to the vehicle and control areas after a 28-day period. In conclusion, Cichorium intybus root extract showed protective and restructuring effects on the skin and stands out as an innovative ingredient to improve skin barrier function.

Panx3 links body mass index and tumorigenesis in a genetically heterogeneous mouse model of carcinogen-induced cancer

Background

Body mass index (BMI) has been implicated as a primary factor influencing cancer development. However, understanding the relationship between these two complex traits has been confounded by both environmental and genetic heterogeneity.

Methods

In order to gain insight into the genetic factors linking BMI and cancer, we performed chemical carcinogenesis on a genetically heterogeneous cohort of interspecific backcross mice ((Mus Spretus × FVB/N) F1 × FVB/N). Using this cohort, we performed quantitative trait loci (QTL) analysis to identify regions linked to BMI. We then performed an integrated analysis incorporating gene expression, sequence comparison between strains, and gene expression network analysis to identify candidate genes influencing both tumor development and BMI.

Results

Analysis of QTL linked to tumorigenesis and BMI identified several loci associated with both phenotypes. Exploring these loci in greater detail revealed a novel relationship between the Pannexin 3 gene (Panx3) and both BMI and tumorigenesis. Panx3 is positively associated with BMI and is strongly tied to a lipid metabolism gene expression network. Pre-treatment Panx3 gene expression levels in normal skin are associated with tumor susceptibility and inhibition of Panx function strongly influences inflammation.

Conclusions

These studies have identified several genetic loci that influence both BMI and carcinogenesis and implicate Panx3 as a candidate gene that links these phenotypes through its effects on inflammation and lipid metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0334-8) contains supplementary material, which is available to authorized users.

Gene Expression Architecture of Mouse Dorsal and Tail Skin Reveals Functional Differences in Inflammation and Cancer

Inherited germline polymorphisms can cause gene expression levels in normal tissues to differ substantially between individuals. We present an analysis of the genetic architecture of normal adult skin from 470 genetically unique mice, demonstrating the effect of germline variants, skin tissue location, and perturbation by exogenous inflammation or tumorigenesis on gene signaling pathways. Gene networks related to specific cell types and signaling pathways, including sonic hedgehog (Shh), Wnt, Lgr family stem cell markers, and keratins, differed at these tissue sites, suggesting mechanisms for the differential susceptibility of dorsal and tail skin to development of skin diseases and tumorigenesis. The Pten tumor suppressor gene network is rewired in premalignant tumors compared to normal tissue, but this response to perturbation is lost during malignant progression. We present a software package for expression quantitative trait loci (eQTL) network analysis and demonstrate how network analysis of whole tissues provides insights into interactions between cell compartments and signaling molecules.

Dane M, Esch A, Feiler HS, Wang NJ, Hardwicke MA, Laquerre S, Jackson J, W. Wood K, Weber B, Spellman PT, Aparicio S, Wooster R, Caldas C, Gray JW. Genome co-amplification upregulates a mitotic gene network activity that predicts outcome and response to mitotic protein inhibitors in breast cancer

BACKGROUND

High mitotic activity is associated with the genesis and progression of many cancers. Small molecule inhibitors of mitotic apparatus proteins are now being developed and evaluated clinically as anticancer agents. With clinical trials of several of these experimental compounds underway, it is important to understand the molecular mechanisms that determine high mitotic activity, identify tumor subtypes that carry molecular aberrations that confer high mitotic activity, and to develop molecular markers that distinguish which tumors will be most responsive to mitotic apparatus inhibitors.

METHODS

We identified a coordinately regulated mitotic apparatus network by analyzing gene expression profiles for 53 malignant and non-malignant human breast cancer cell lines and two separate primary breast tumor datasets. We defined the mitotic network activity index (MNAI) as the sum of the transcriptional levels of the 54 coordinately regulated mitotic apparatus genes. The effect of those genes on cell growth was evaluated by small interfering RNA (siRNA).

RESULTS

High MNAI was enriched in basal-like breast tumors and was associated with reduced survival duration and preferential sensitivity to inhibitors of the mitotic apparatus proteins, polo-like kinase, centromere associated protein E and aurora kinase designated GSK462364, GSK923295 and GSK1070916, respectively. Co-amplification of regions of chromosomes 8q24, 10p15-p12, 12p13, and 17q24-q25 was associated with the transcriptional upregulation of this network of 54 mitotic apparatus genes, and we identify transcription factors that localize to these regions and putatively regulate mitotic activity. Knockdown of the mitotic network by siRNA identified 22 genes that might be considered as additional therapeutic targets for this clinically relevant patient subgroup.

CONCLUSIONS

We define a molecular signature which may guide therapeutic approaches for tumors with high mitotic network activity.

Sebaceous Gland Atrophy in Psoriasis: An Explanation for Psoriatic Alopecia?

In a transcriptome study of lesional psoriatic skin (PP) versus normal skin, we found a coexpressed gene module (N5) enriched 11.5-fold for lipid biosynthetic genes. We also observed fewer visible hairs in PP skin, compared with uninvolved nonlesional psoriatic skin or normal skin (P < 0.0001). To ask whether these findings might be due to abnormalities of the pilosebaceous unit, we carried out three-dimensional morphometric analysis of paired PP and nonlesional psoriatic skin biopsies. Sebaceous glands were markedly atrophic in PP versus nonlesional psoriatic skin (91% average reduction in volume, P = 0.031). Module N5 genes were strongly downregulated in PP versus normal skin (fold change < 0.25, 44.4-fold) and strongly upregulated in sebaceous hyperplasia (fold change > 4, 54.1-fold). The intersection of PP-downregulated and sebaceous hyperplasia-upregulated gene lists generated a gene expression signature consisting solely of module N5 genes, whose expression in PP versus normal skin was inversely correlated with the signature of IL17-stimulated keratinocytes. Despite loss of visible hairs, morphometry identified elongated follicles in PP versus nonlesional psoriatic skin (average 1.7 vs. 1.2 μm, P = 0.020). These results document sebaceous gland atrophy in nonscalp psoriasis, identify a cytokine-regulated set of sebaceous gland signature genes, and suggest that loss of visible hair in PP skin may result from abnormal sebaceous gland function.

Protective Effect of Topical Vitamin D3 against Photocarcinogenesis in a Murine Model

Background

Although the incidence of non-melanoma skin cancer is increasing, there are no effective practical preventive measures other than avoiding sun exposure.

Objective

To elucidate the protective effect of topical application of biologically active vitamin D₃ (calcitriol) on skin cancer development caused by exposure to ultraviolet (UV).

Methods

Groups of hairless mice were topically treated with either calcitriol or vehicle immediately after exposure to UVB and UVA three times weekly for the initial 20 weeks, and without UV exposure in the following 6 weeks. Tumor number was counted and biopsies were done for histopathologic analysis. The changes of cyclobutane pyrimidine dimer (CPD) were evaluated 1 hour and 11 hours after short term of UV exposure and application of calcitriol. For safety evaluation, blood test and body weights were evaluated at 23rd and 25th week.

Results

Total tumor count and number of tumors less than 3 mm in size tended to be fewer in calcitriol group, and tumors more than 3 mm in size showed significantly lower tumor formation rate in calcitriol group. Single application of calcitriol reduced CPD at 1 hour and 11 hours after UV exposure. Histopathologic analysis showed tumors with lower grade malignancy in calcitriol group which suggested a delay in tumor progression. However, serum levels of calcium and phosphate in calcitriol group were above normal range, and weight loss was found.

Conclusion

Topical calcitriol may suppress the formation and progression of UV-induced non-melanoma skin cancer by enhancing the repair mechanism of UV damage.

A keratin scaffold regulates epidermal barrier formation, mitochondrial lipid composition, and activity

Epidermal keratin filaments are important components and organizers of the cornified envelope and regulate mitochondrial metabolism by modulating their membrane composition.

Keratin intermediate filaments (KIFs) protect the epidermis against mechanical force, support strong adhesion, help barrier formation, and regulate growth. The mechanisms by which type I and II keratins contribute to these functions remain incompletely understood. Here, we report that mice lacking all type I or type II keratins display severe barrier defects and fragile skin, leading to perinatal mortality with full penetrance. Comparative proteomics of cornified envelopes (CEs) from prenatal KtyI^−/− and KtyII^−/−_K8 mice demonstrates that absence of KIF causes dysregulation of many CE constituents, including downregulation of desmoglein 1. Despite persistence of loricrin expression and upregulation of many Nrf2 targets, including CE components Sprr2d and Sprr2h, extensive barrier defects persist, identifying keratins as essential CE scaffolds. Furthermore, we show that KIFs control mitochondrial lipid composition and activity in a cell-intrinsic manner. Therefore, our study explains the complexity of keratinopathies accompanied by barrier disorders by linking keratin scaffolds to mitochondria, adhesion, and CE formation.

HZE Radiation Non-Targeted Effects on the Microenvironment That Mediate Mammary Carcinogenesis

Clear mechanistic understanding of the biological processes elicited by radiation that increase cancer risk can be used to inform prediction of health consequences of medical uses, such as radiotherapy, or occupational exposures, such as those of astronauts during deep space travel. Here, we review the current concepts of carcinogenesis as a multicellular process during which transformed cells escape normal tissue controls, including the immune system, and establish a tumor microenvironment. We discuss the contribution of two broad classes of radiation effects that may increase cancer: radiation targeted effects that occur as a result of direct energy deposition, e.g., DNA damage, and non-targeted effects (NTE) that result from changes in cell signaling, e.g., genomic instability. It is unknown whether the potentially greater carcinogenic effect of high Z and energy (HZE) particle radiation is a function of the relative contribution or extent of NTE or due to unique NTE. We addressed this problem using a radiation/genetic mammary chimera mouse model of breast cancer. Our experiments suggest that NTE promote more aggressive cancers, as evidenced by increased growth rate, transcriptomic signatures, and metastasis, and that HZE particle NTE are more effective than reference γ-radiation. Emerging evidence suggest that HZE irradiation dampens antitumor immunity. These studies raise concern that HZE radiation exposure not only increases the likelihood of developing cancer but also could promote progression to more aggressive cancer with a greater risk of mortality.

Animal Models of Chemical Carcinogenesis: Driving Breakthroughs in Cancer Research for 100 Years

The identification of carcinogens in the workplace, diet, and environment through chemical carcinogenesis studies in animals has directly contributed to a reduction of cancer burden in the human population. Reduced exposure to these carcinogens through lifestyle changes, government regulation, or change in industry practices has reduced cancer incidence in exposed populations. In addition to providing the first experimental evidence for cancer's relationship to chemical and radiation exposure, animal models of environmentally induced cancer have and will continue to provide important insight into the causes, mechanisms, and conceptual frameworks of cancer. More recently, combining chemical carcinogens with genetically engineered mouse models has emerged as an invaluable approach to study the complex interaction between genotype and environment that contributes to cancer development. In the future, animal models of environmentally induced cancer are likely to provide insight into areas such as the epigenetic basis of cancer, genetic modifiers of cancer susceptibility, the systems biology of cancer, inflammation and cancer, and cancer prevention.

Equalizer reduces SNP bias in Affymetrix microarrays

BACKGROUND

Gene expression microarrays measure the levels of messenger ribonucleic acid (mRNA) in a sample using probe sequences that hybridize with transcribed regions. These probe sequences are designed using a reference genome for the relevant species. However, most model organisms and all humans have genomes that deviate from their reference. These variations, which include single nucleotide polymorphisms, insertions of additional nucleotides, and nucleotide deletions, can affect the microarray's performance. Genetic experiments comparing individuals bearing different population-associated single nucleotide polymorphisms that intersect microarray probes are therefore subject to systemic bias, as the reduction in binding efficiency due to a technical artifact is confounded with genetic differences between parental strains. This problem has been recognized for some time, and earlier methods of compensation have attempted to identify probes affected by genome variants using statistical models. These methods may require replicate microarray measurement of gene expression in the relevant tissue in inbred parental samples, which are not always available in model organisms and are never available in humans.

RESULTS

By using sequence information for the genomes of organisms under investigation, potentially problematic probes can now be identified a priori. However, there is no published software tool that makes it easy to eliminate these probes from an annotation. I present equalizer, a software package that uses genome variant data to modify annotation files for the commonly used Affymetrix IVT and Gene/Exon platforms. These files can be used by any microarray normalization method for subsequent analysis. I demonstrate how use of equalizer on experiments mapping germline influence on gene expression in a genetic cross between two divergent mouse species and in human samples significantly reduces probe hybridization-induced bias, reducing false positive and false negative findings.

CONCLUSIONS

The equalizer package reduces probe hybridization bias from experiments performed on the Affymetrix microarray platform, allowing accurate assessment of germline influence on gene expression.

The intracellular domains of Notch1 and Notch2 are functionally equivalent during development and carcinogenesis

Although Notch1 and Notch2 are closely related paralogs and function through the same canonical signaling pathway, they contribute to different outcomes in some cell and disease contexts. To understand the basis for these differences, we examined in detail mice in which the Notch intracellular domains (N1ICD and N2ICD) were swapped. Our data indicate that strength (defined here as the ultimate number of intracellular domain molecules reaching the nucleus, integrating ligand-mediated release and nuclear translocation) and duration (half-life of NICD-RBPjk-MAML-DNA complexes, integrating cooperativity and stability dependent on shared sequence elements) are the factors that underlie many of the differences between Notch1 and Notch2 in all the contexts we examined, including T-cell development, skin differentiation and carcinogenesis, the inner ear, the lung and the retina. We were able to show that phenotypes in the heart, endothelium, and marginal zone B cells are attributed to haploinsufficiency but not to intracellular domain composition. Tissue-specific differences in NICD stability were most likely caused by alternative scissile bond choices by tissue-specific γ-secretase complexes following the intracellular domain swap. Reinterpretation of clinical findings based on our analyses suggests that differences in outcome segregating with Notch1 or Notch2 are likely to reflect outcomes dependent on the overall strength of Notch signals.

The genetics of splicing in neuroblastoma

Regulation of mRNA splicing, a critical and tightly regulated cellular function, underlies the majority of proteomic diversity and is frequently disrupted in disease. Using an integrative genomics approach, we combined both genomic data and exon-level transcriptome data in two somatic tissues (cerebella and peripheral ganglia) from a transgenic mouse model of neuroblastoma, a tumor that arises from the peripheral neural crest. Here, we describe splicing quantitative trait loci associated with differential splicing across the genome that we use to identify genes with previously unknown functions within the splicing pathway and to define de novo intronic splicing motifs that influence splicing from hundreds of bases away. Our results show that these splicing motifs represent sites for functional recurrent mutations and highlight novel candidate genes in human cancers, including childhood neuroblastoma.

SIGNIFICANCE

Somatic mutations with predictable downstream effects are largely relegated to coding regions, which comprise less than 2% of the human genome. Using an unbiased in vivo analysis of a mouse model of neuroblastoma, we have identified intronic splicing motifs that translate into sites for recurrent somatic mutations in human cancers.

Genetic variation in the functional ENG allele inherited from the non-affected parent associates with presence of pulmonary arteriovenous malformation in hereditary hemorrhagic telangiectasia 1 (HHT1) and may influence expression of PTPN14

HHT shows clinical variability within and between families. Organ site and prevalence of arteriovenous malformations (AVMs) depend on the HHT causative gene and on environmental and genetic modifiers. We tested whether variation in the functional ENG allele, inherited from the unaffected parent, alters risk for pulmonary AVM in HHT1 mutation carriers who are ENG haploinsufficient. Genetic association was found between rs10987746 of the wild type ENG allele and presence of pulmonary AVM [relative risk = 1.3 (1.0018-1.7424)]. The rs10987746-C at-risk allele associated with lower expression of ENG RNA in a panel of human lymphoblastoid cell lines (P = 0.004). Moreover, in angiogenically active human lung adenocarcinoma tissue, but not in uninvolved quiescent lung, rs10987746-C was correlated with expression of PTPN14 (P = 0.004), another modifier of HHT. Quantitative TAQMAN expression analysis in a panel of normal lung tissues from 69 genetically heterogeneous inter-specific backcross mice, demonstrated strong correlation between expression levels of Eng, Acvrl1, and Ptpn14 (r2 = 0.75-0.9, P < 1 × 10(-12)), further suggesting a direct or indirect interaction between these three genes in lung in vivo. Our data indicate that genetic variation within the single functional ENG gene influences quantitative and/or qualitative differences in ENG expression that contribute to risk of pulmonary AVM in HHT1, and provide correlative support for PTPN14 involvement in endoglin/ALK1 lung biology in vivo. PTPN14 has been shown to be a negative regulator of Yap/Taz signaling, which is implicated in mechanotransduction, providing a possible molecular link between endoglin/ALK1 signaling and mechanical stress. EMILIN2, which showed suggestive genetic association with pulmonary AVM, is also reported to interact with Taz in angiogenesis. Elucidation of the molecular mechanisms regulating these interactions in endothelial cells may ultimately provide more rational choices for HHT therapy.

Chaperone Hsp47 Drives Malignant Growth and Invasion by Modulating an ECM Gene Network

The extracellular matrix (ECM) is a determining factor in the tumor microenvironment that restrains or promotes malignant growth. In this report, we show how the molecular chaperone protein Hsp47 functions as a nodal hub in regulating an ECM gene transcription network. A transcription network analysis showed that Hsp47 expression was activated during breast cancer development and progression. Hsp47 silencing reprogrammed human breast cancer cells to form growth-arrested and/or noninvasive structures in 3D cultures, and to limit tumor growth in xenograft assays by reducing deposition of collagen and fibronectin. Coexpression network analysis also showed that levels of microRNA(miR)-29b and -29c were inversely correlated with expression of Hsp47 and ECM network genes in human breast cancer tissues. We found that miR-29 repressed expression of Hsp47 along with multiple ECM network genes. Ectopic expression of miR-29b suppressed malignant phenotypes of breast cancer cells in 3D culture. Clinically, increased expression of Hsp47 and reduced levels of miR-29b and -29c were associated with poor survival outcomes in breast cancer patients. Our results show that Hsp47 is regulated by miR-29 during breast cancer development and progression, and that increased Hsp47 expression promotes cancer progression in part by enhancing deposition of ECM proteins.

Expression quantitative trait loci and receptor pharmacology implicate Arg1 and the GABA-A receptor as therapeutic targets in neuroblastoma

The development of targeted therapeutics for neuroblastoma, the third most common tumor in children, has been limited by a poor understanding of growth signaling mechanisms unique to the peripheral nerve precursors from which tumors arise. In this study, we combined genetics with gene-expression analysis in the peripheral sympathetic nervous system to implicate arginase 1 and GABA signaling in tumor formation in vivo. In human neuroblastoma cells, either blockade of ARG1 or benzodiazepine-mediated activation of GABA-A receptors induced apoptosis and inhibited mitogenic signaling through AKT and MAPK. These results suggest that ARG1 and GABA influence both neural development and neuroblastoma and that benzodiazepines in clinical use may have potential applications for neuroblastoma therapy.

Refining the role for adult stem cells as cancer cells of origin

Significant progress has been made to identify the cells at the foundation of tumorigenesis, the cancer cell of origin (CCO). The majority of data points towards resident adult stem cells (ASCs) or primitive progenitors as the CCO for those cancers studied, highlighting the importance of stem cells not only as propagators but also as initiators of cancer. Recent data suggest tumor initiation at the CCOs can be regulated through both intrinsic and extrinsic signals and that the identity of the CCOs and their propensity to initiate tumorigenesis is context dependent. In this review, we summarize some of the recent findings regarding CCOs and solid tumor initiation and highlight its relation with bona fide human cancer.

Modeling cutaneous squamous carcinoma development in the mouse

Cutaneous squamous cell carcinoma (SCC) is one of the most common cancers in Caucasian populations and is associated with a significant risk of morbidity and mortality. The classic mouse model for studying SCC involves two-stage chemical carcinogenesis, which has been instrumental in the evolution of the concept of multistage carcinogenesis, as widely applied to both human and mouse cancers. Much is now known about the sequence of biological and genetic events that occur in this skin carcinogenesis model and the factors that can influence the course of tumor development, such as perturbations in the oncogene/tumor-suppressor signaling pathways involved, the nature of the target cell that acquires the first genetic hit, and the role of inflammation. Increasingly, studies of tumor-initiating cells, malignant progression, and metastasis in mouse skin cancer models will have the potential to inform future approaches to treatment and chemoprevention of human squamous malignancies.

RNA-seq permits a closer look at normal skin and psoriasis gene networks

Psoriasis is a chronic inflammatory skin disease driven by aberrant signals from the immune system. In this issue Li et al. present the first large RNA-seq analysis of gene expression in normal skin and psoriasis lesions, providing a more comprehensive view of mRNA expression than earlier microarray studies. This study’s size enables gene co-expression analysis, a method illustrating which pathways are altered by the presence of disease.

Mouse and human BAC transgenes recapitulate tissue-specific expression of the vitamin D receptor in mice and rescue the VDR-null phenotype

The biological actions of 1,25-dihydroxyvitamin D3 (1,25(OH)2D3) are mediated by the vitamin D receptor (VDR), which is expressed in numerous target tissues in a cell type-selective manner. Recent studies using genomic analyses and recombineered bacterial artificial chromosomes (BACs) have defined the specific features of mouse and human VDR gene loci in vitro. In the current study, we introduced recombineered mouse and human VDR BACs as transgenes into mice and explored their expression capabilities in vivo. Individual transgenic mouse strains selectively expressed BAC-derived mouse or human VDR proteins in appropriate vitamin D target tissues, thereby recapitulating the tissue-specific expression of endogenous mouse VDR. The mouse VDR transgene was also regulated by 1,25(OH)2D3 and dibutyryl-cAMP. When crossed into a VDR-null mouse background, both transgenes restored wild-type basal as well as 1,25(OH)2D3-inducible gene expression patterns in the appropriate tissues. This maneuver resulted in the complete rescue of the aberrant phenotype noted in the VDR-null mouse, including systemic features associated with altered calcium and phosphorus homeostasis and disrupted production of parathyroid hormone and fibroblast growth factor 23, and abnormalities associated with the skeleton, kidney, parathyroid gland, and the skin. This study suggests that both mouse and human VDR transgenes are capable of recapitulating basal and regulated expression of the VDR in the appropriate mouse tissues and restore 1,25(OH)2D3 function. These results provide a baseline for further dissection of mechanisms integral to mouse and human VDR gene expression and offer the potential to explore the consequence of selective mutations in VDR proteins in vivo.

Principles and methods of integrative genomic analyses in cancer

Combined analyses of molecular data, such as DNA copy-number alteration, mRNA and protein expression, point to biological functions and molecular pathways being deregulated in multiple cancers. Genomic, metabolomic and clinical data from various solid cancers and model systems are emerging and can be used to identify novel patient subgroups for tailored therapy and monitoring. The integrative genomics methodologies that are used to interpret these data require expertise in different disciplines, such as biology, medicine, mathematics, statistics and bioinformatics, and they can seem daunting. The objectives, methods and computational tools of integrative genomics that are available to date are reviewed here, as is their implementation in cancer research.

Identification of Hipk2 as an essential regulator of white fat development

Homeodomain-interacting protein kinase 2 (Hipk2) has previously been implicated in the control of several transcription factors involved in embryonic development, apoptosis, cell proliferation, and tumor development, but very little is understood about the exact mechanisms through which Hipk2 influences these processes. Analysis of gene expression in normal tissues from genetically heterogeneous mouse or human populations can reveal network motifs associated with the structural or functional components of the tissue, and may predict roles for genes of unknown function. Here we have applied this network strategy to uncover a role for the Hipk2 gene in the transcriptional system controlling adipogenesis. Both in vitro and in vivo models were used to show that knockdown or loss of Hipk2 specifically inhibits white adipose cell differentiation and tissue development. In addition, loss of Hipk2 leads to induction of pockets of multilocular brown fat-like cells in remaining white adipose depots, which express markers of brown and beige fat such as uncoupling protein 1 and transmembrane protein 26. These changes are accompanied by increased insulin sensitivity in Hipk2 knockout mice and reduced high-fat diet-induced weight gain, highlighting a potential role for this kinase in diseases such as diabetes and obesity. Our study underscores the versatility and power of a readily available tissue, such as skin, for network modeling of systemic transcriptional programs involved in multiple pathways, including lipid metabolism and adipogenesis.

Emerging role for leucine-rich repeat-containing G-protein-coupled receptors LGR5 and LGR4 in cancer stem cells

The concept of cancer stem cells has gained considerable interest in the last few decades, partly because of their potential implication in therapy resistance. However, the lack of specific cellular surface markers for these cells has impeded their isolation, making the characterization of this cellular subpopulation technically challenging. Recent studies have indicated that leucine-rich repeat-containing G-protein-coupled receptor 4 and 5 (LGR4 and LGR5) expression in multiple organs may represent a global marker of adult stem cells. This review aims to give an overview of LGR4 and LGR5 as cancer stem cell markers and their function in development.

Hair follicle disruption facilitates pathogenesis to UVB-induced cutaneous inflammation and basal cell carcinoma development in Ptch(+/-) mice

Hairless mice carrying homozygous mutations in hairless gene manifest rudimentary hair follicles (HFs), epidermal cysts, hairless phenotype, and enhanced susceptibility to squamous cell carcinomas. However, their susceptibility to basal cell carcinomas (BCCs), a neoplasm considered originated from HF-localized stem cells, is unknown. To demonstrate the role of HFs in BCC development, we bred Ptch(+/-)/C57BL6 with SKH-1 hairless mice, followed by brother-sister cross to get F2 homozygous mutant (hairless) or wild-type (haired) mice. UVB-induced inflammation was less pronounced in shaved haired than in hairless mice. In hairless mice, inflammatory infiltrate was found around the rudimentary HFs and epidermal cysts. Expression of epidermal IL1f6, S100a8, vitamin D receptor, repetin, and major histocompatibility complex II, biomarkers depicting susceptibility to cutaneous inflammation, was also higher. In these animals, HF disruption altered susceptibility to UVB-induced BCCs. Tumor onset in hairless mice was 10 weeks earlier than in haired littermates. The incidence of BCCs was significantly higher in hairless than in haired animals; however, the magnitude of sonic hedgehog signaling did not differ significantly. Overall, 100% of hairless mice developed >12 tumors per mouse after 32 weeks of UVB therapy, whereas haired mice developed fewer than three tumors per mouse after 44 weeks of long-term UVB irradiation. Tumors in hairless mice were more aggressive than in haired littermates and manifested decreased E-cadherin and enhanced mesenchymal proteins. These data provide novel evidence that disruption of HFs in Ptch(+/-) mice enhances cutaneous susceptibility to inflammation and BCCs.

Systems biology perspectives on the carcinogenic potential of radiation

This review focuses on recent experimental and modeling studies that attempt to define the physiological context in which high linear energy transfer (LET) radiation increases epithelial cancer risk and the efficiency with which it does so. Radiation carcinogenesis is a two-compartment problem: ionizing radiation can alter genomic sequence as a result of damage due to targeted effects (TE) from the interaction of energy and DNA; it can also alter phenotype and multicellular interactions that contribute to cancer by poorly understood non-targeted effects (NTE). Rather than being secondary to DNA damage and mutations that can initiate cancer, radiation NTE create the critical context in which to promote cancer. Systems biology modeling using comprehensive experimental data that integrates different levels of biological organization and time-scales is a means of identifying the key processes underlying the carcinogenic potential of high-LET radiation. We hypothesize that inflammation is a key process, and thus cancer susceptibility will depend on specific genetic predisposition to the type and duration of this response. Systems genetics using novel mouse models can be used to identify such determinants of susceptibility to cancer in radiation sensitive tissues following high-LET radiation. Improved understanding of radiation carcinogenesis achieved by defining the relative contribution of NTE carcinogenic effects and identifying the genetic determinants of the high-LET cancer susceptibility will help reduce uncertainties in radiation risk assessment.

Transcriptome analysis of psoriasis in a large case-control sample: RNA-seq provides insights into disease mechanisms

To increase our understanding of psoriasis, we utilized RNA-seq to assay the transcriptomes of lesional psoriatic and normal skin. We sequenced polyadenylated RNA-derived cDNAs from 92 psoriatic and 82 normal punch biopsies, generating an average of ~38 million single-end 80-bp reads per sample. Comparison of 42 samples examined by both RNA-seq and microarray revealed marked differences in sensitivity, with transcripts identified only by RNA-seq having much lower expression than those also identified by microarray. RNA-seq identified many more differentially expressed transcripts enriched in immune system processes. Weighted gene co-expression network analysis (WGCNA) revealed multiple modules of coordinately expressed epidermal differentiation genes, overlapping significantly with genes regulated by the long non-coding RNA TINCR, its target gene, staufen-1 (STAU1), the p63 target gene ZNF750, and its target KLF4. Other coordinately expressed modules were enriched for lymphoid and/or myeloid signature transcripts and genes induced by IL-17 in keratinocytes. Dermally-expressed genes were significantly down-regulated in psoriatic biopsies, most likely due to expansion of the epidermal compartment. These results demonstrate the power of WGCNA to elucidate gene regulatory circuits in psoriasis, and emphasize the influence of tissue architecture in both differential expression and co-expression analysis.

Germline sequence variants in TGM3 and RGS22 confer risk of basal cell carcinoma

To search for new sequence variants that confer risk of cutaneous basal cell carcinoma (BCC), we conducted a genome-wide association study of 38.5 million single nucleotide polymorphisms (SNPs) and small indels identified through whole-genome sequencing of 2230 Icelanders. We imputed genotypes for 4208 BCC patients and 109 408 controls using Illumina SNP chip typing data, carried out association tests and replicated the findings in independent population samples. We found new BCC susceptibility loci at TGM3 (rs214782[G], P = 5.5 × 10⁻¹⁷, OR = 1.29) and RGS22 (rs7006527[C], P = 8.7 × 10⁻¹³, OR = 0.77). TGM3 encodes transglutaminase type 3, which plays a key role in production of the cornified envelope during epidermal differentiation.

Beyond expectations: novel insights into epidermal keratin function and regulation

The epidermis is a stratified epithelium that relies on its cytoskeleton and cell junctions to protect the body against mechanical injury, dehydration, and infections. Keratin intermediate filament proteins are involved in many of these functions by forming cell-specific cytoskeletal scaffolds crucial for the maintenance of cell and tissue integrity. In response to various stresses, the expression and organization of keratins are altered at transcriptional and posttranslational levels to restore tissue homeostasis. Failure to restore tissue homeostasis in the presence of keratin gene mutations results in acute and chronic skin disorders for which currently no rational therapies are available. Here, we review the recent progress on the role of keratins in cytoarchitecture, adhesion, signaling, and inflammation. By focusing on epidermal keratins, we illustrate the contribution of keratin isotypes to differentiated epithelial functions.

An integrated systems genetics screen reveals the transcriptional structure of inherited predisposition to metastatic disease

Metastasis is the result of stochastic genomic and epigenetic events leading to gene expression profiles that drive tumor dissemination. Here we exploit the principle that metastatic propensity is modified by the genetic background to generate prognostic gene expression signatures that illuminate regulators of metastasis. We also identify multiple microRNAs whose germline variation is causally linked to tumor progression and metastasis. We employ network analysis of global gene expression profiles in tumors derived from a panel of recombinant inbred mice to identify a network of co-expressed genes centered on Cnot2 that predicts metastasis-free survival. Modulating Cnot2 expression changes tumor cell metastatic potential in vivo, supporting a functional role for Cnot2 in metastasis. Small RNA sequencing of the same tumor set revealed a negative correlation between expression of the Mir216/217 cluster and tumor progression. Expression quantitative trait locus analysis (eQTL) identified cis-eQTLs at the Mir216/217 locus, indicating that differences in expression may be inherited. Ectopic expression of Mir216/217 in tumor cells suppressed metastasis in vivo. Finally, small RNA sequencing and mRNA expression profiling data were integrated to reveal that miR-3470a/b target a high proportion of network transcripts. In vivo analysis of Mir3470a/b demonstrated that both promote metastasis. Moreover, Mir3470b is a likely regulator of the Cnot2 network as its overexpression down-regulated expression of network hub genes and enhanced metastasis in vivo, phenocopying Cnot2 knockdown. The resulting data from this strategy identify Cnot2 as a novel regulator of metastasis and demonstrate the power of our systems-level approach in identifying modifiers of metastasis.