Systems Level Analysis of Systemic Sclerosis Shows a Network of Immune and Profibrotic Pathways Connected with Genetic Polymorphisms

Systemic sclerosis (SSc) is a rare systemic autoimmune disease characterized by skin and organ fibrosis. The pathogenesis of SSc and its progression are poorly understood. The SSc intrinsic gene expression subsets (inflammatory, fibroproliferative, normal-like, and limited) are observed in multiple clinical cohorts of patients with SSc. Analysis of longitudinal skin biopsies suggests that a patient''s subset assignment is stable over 6–12 months. Genetically, SSc is multi-factorial with many genetic risk loci for SSc generally and for specific clinical manifestations. Here we identify the genes consistently associated with the intrinsic subsets across three independent cohorts, show the relationship between these genes using a gene-gene interaction network, and place the genetic risk loci in the context of the intrinsic subsets. To identify gene expression modules common to three independent datasets from three different clinical centers, we developed a consensus clustering procedure based on mutual information of partitions, an information theory concept, and performed a meta-analysis of these genome-wide gene expression datasets. We created a gene-gene interaction network of the conserved molecular features across the intrinsic subsets and analyzed their connections with SSc-associated genetic polymorphisms. The network is composed of distinct, but interconnected, components related to interferon activation, M2 macrophages, adaptive immunity, extracellular matrix remodeling, and cell proliferation. The network shows extensive connections between the inflammatory- and fibroproliferative-specific genes. The network also shows connections between these subset-specific genes and 30 SSc-associated polymorphic genes including STAT4, BLK, IRF7, NOTCH4, PLAUR, CSK, IRAK1, and several human leukocyte antigen (HLA) genes. Our analyses suggest that the gene expression changes underlying the SSc subsets may be long-lived, but mechanistically interconnected and related to a patients underlying genetic risk.     

Expression profiling of major histocompatibility and natural killer complex genes reveals candidates for controlling risk of graft versus host disease

Background

The major histocompatibility complex (MHC) is the most important genomic region that contributes to the risk of graft versus host disease (GVHD) after haematopoietic stem cell transplantation. Matching of MHC class I and II genes is essential for the success of transplantation. However, the MHC contains additional genes that also contribute to the risk of developing acute GVHD. It is difficult to identify these genes by genetic association studies alone due to linkage disequilibrium in this region. Therefore, we aimed to identify MHC genes and other genes involved in the pathophysiology of GVHD by mRNA expression profiling.

Methodology/Principal Findings

To reduce the complexity of the task, we used genetically well-defined rat inbred strains and a rat skin explant assay, an in-vitro-model of the graft versus host reaction (GVHR), to analyze the expression of MHC, natural killer complex (NKC), and other genes in cutaneous GVHR. We observed a statistically significant and strong up or down regulation of 11 MHC, 6 NKC, and 168 genes encoded in other genomic regions, i.e. 4.9%, 14.0%, and 2.6% of the tested genes respectively. The regulation of 7 selected MHC and 3 NKC genes was confirmed by quantitative real-time PCR and in independent skin explant assays. In addition, similar regulations of most of the selected genes were observed in GVHD-affected skin lesions of transplanted rats and in human skin explant assays.

Conclusions/Significance

We identified rat and human MHC and NKC genes that are regulated during GVHR in skin explant assays and could therefore serve as biomarkers for GVHD. Several of the respective human genes, including HLA-DMB, C2, AIF1, SPR1, UBD, and OLR1, are polymorphic. These candidates may therefore contribute to the genetic risk of GVHD in patients.     

Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.     

Biological processes, properties and molecular wiring diagrams of candidate low-penetrance breast cancer susceptibility genes

Background

The process of malignant transformation, progression and metastasis of melanoma is poorly understood. Gene expression profiling of human cancer has allowed for a unique insight into the genes that are involved in these processes. Thus, we have attempted to utilize this approach through the analysis of a series of primary, non-metastatic cutaneous tumors and metastatic melanoma samples.

Methods

We have utilized gene microarray analysis and a variety of molecular techniques to compare 40 metastatic melanoma (MM) samples, composed of 22 bulky, macroscopic (replaced) lymph node metastases, 16 subcutaneous and 2 distant metastases (adrenal and brain), to 42 primary cutaneous cancers, comprised of 16 melanoma, 11 squamous cell, 15 basal cell skin cancers. A Human Genome U133 Plus 2.0 array from Affymetrix, Inc. was utilized for each sample. A variety of statistical software, including the Affymetrix MAS 5.0 analysis software, was utilized to compare primary cancers to metastatic melanomas. Separate analyses were performed to directly compare only primary melanoma to metastatic melanoma samples. The expression levels of putative oncogenes and tumor suppressor genes were analyzed by semi- and real-time quantitative RT-PCR (qPCR) and Western blot analysis was performed on select genes.

Results

We find that primary basal cell carcinomas, squamous cell carcinomas and thin melanomas express dramatically higher levels of many genes, including SPRR1A/B, KRT16/17, CD24, LOR, GATA3, MUC15, and TMPRSS4, than metastatic melanoma. In contrast, the metastatic melanomas express higher levels of genes such as MAGE, GPR19, BCL2A1, MMP14, SOX5, BUB1, RGS20, and more. The transition from non-metastatic expression levels to metastatic expression levels occurs as melanoma tumors thicken. We further evaluated primary melanomas of varying Breslow's tumor thickness to determine that the transition in expression occurs at different thicknesses for different genes suggesting that the "transition zone" represents a critical time for the emergence of the metastatic phenotype. Several putative tumor oncogenes (SPP-1, MITF, CITED-1, GDF-15, c-Met, HOX loci) and suppressor genes (PITX-1, CST-6, PDGFRL, DSC-3, POU2F3, CLCA2, ST7L), were identified and validated by quantitative PCR as changing expression during this transition period. These are strong candidates for genes involved in the progression or suppression of the metastatic phenotype.

Conclusion

The gene expression profiling of primary, non-metastatic cutaneous tumors and metastatic melanoma has resulted in the identification of several genes that may be centrally involved in the progression and metastatic potential of melanoma. This has very important implications as we continue to develop an improved understanding of the metastatic process, allowing us to identify specific genes for prognostic markers and possibly for targeted therapeutic approaches.     

Transposon mutagenesis identifies cooperating genetic drivers during keratinocyte transformation and cutaneous squamous cell carcinoma progression

The systematic identification of genetic events driving cellular transformation and tumor progression in the absence of a highly recurrent oncogenic driver mutation is a challenge in cutaneous oncology. In cutaneous squamous cell carcinoma (cuSCC), the high UV-induced mutational burden poses a hurdle to achieve a complete molecular landscape of this disease. Here, we utilized the Sleeping Beauty transposon mutagenesis system to statistically define drivers of keratinocyte transformation and cuSCC progression in vivo in the absence of UV-IR, and identified both known tumor suppressor genes and novel oncogenic drivers of cuSCC. Functional analysis confirms an oncogenic role for the ZMIZ genes, and tumor suppressive roles for KMT2C, CREBBP and NCOA2, in the initiation or progression of human cuSCC. Taken together, our in vivo screen demonstrates an extremely heterogeneous genetic landscape of cuSCC initiation and progression, which can be harnessed to better understand skin oncogenic etiology and prioritize therapeutic candidates.     

Regulation of cell cycle and differentiation markers by pathogenic,non-pathogenic and opportunistic skin bacteria

Skin is the first line of defense against the physical, chemical and the biological environment. It is an ideal organ for studying molecular responses to biological infections through a variety of skin cells that specialize in immune responses. Comparative analysis of skin response to pathogenic, non-pathogenic, and commensal bacteria would help in the identification of disease specific pathways for drug targets. In this study, we investigated human breast reduction skin responses to Cutibacterium acnes (C. acnes), Staphylococcus aureus (S. aureus), Staphylococcus epidermidis (S. epidermidis), and TLR1/2 agonist using Affymetrix microarray chips. The Pam3CSK4 solution and bacterial cultures were prepared and inoculated in steel rings, that were placed on the acetone treated epidermis in a petri dish. After 24 h incubation, 8 mm punch biopsies were taken from the center of the ring, and RNA was extracted. The genome-wide expression was then analyzed using Affymetrix HG-133A gene chip microarray. We found that the C. acnes and S. aureus boosted the production of extracellular matrix components and attenuated the expression of differentiation markers. The above responses were mediated through the TLR2 pathway. Skin also responded to S. aureus and C. acnes by inducing the genes of the cell cycle machinery; this response was not TLR2-dependent. S. aureus induced, whereas C. acnes suppressed the genes associated with apoptosis; this was also not TLR2-dependent. Moreover, S. epidermis apparently did not lead to changes in gene expression. We conclude that the breast reduction skin is a very useful model to study the global gene expression in response to bacterial treatments.     

Genome-Wide Analysis Using Exon Arrays Demonstrates an Important Role for Expression of Extra-Cellular Matrix,Fibrotic Control and Tissue Remodelling Genes in Dupuytren's Disease

Dupuytren''s disease (DD) is a classic example of pathological fibrosis which results in a debilitating disorder affecting a large sector of the human population. It is characterized by excessive local proliferation of fibroblasts and over-production of collagen and other components of extracellular matrix (ECM) in the palmar fascia. The fibrosis progressively results in contracture of elements between the palmar fascia and skin causing flexion deformity or clawing of the fingers and a severe reduction in hand function. While much is known about the pathogenesis and surgical treatment of DD, little is known about the factors that cause its onset and progression, despite many years of research. Gene expression patterns in DD patients now offers the potential to identify genes that direct the pathogenesis of DD. In this study we used primary cultures of fibroblasts derived from excisional biopsies of fibrotic tissue from DD patients to compare the gene expression profiles on a genome-wide basis with normal control fibroblasts. Our investigations have identified genes that may be involved with DD pathogenesis including some which are directly relevant to fibrosis. In particular, these include significantly reduced expression levels of three matrix metallopeptidases (MMP1, MMP3, MMP16), follistatin, and STAT1, and significantly increased expression levels of fibroblast growth factors (FGF9, FGF11), a number of collagen genes and other ECM genes in DD patient samples. Many of these gene products are known to be involved in fibrosis, tumour formation and in the normal processes of tissue remodelling. In addition, alternative splicing was identified in some DD associated genes. These highly sensitive genomic investigations provide new insight into the molecular mechanisms that may underpin the development and progression of DD.     

Batch effect correction for genome-wide methylation data with Illumina Infinium platform

Background

Paclitaxel is a microtubule-stabilizing drug that has been commonly used in treating cancer. Due to genetic heterogeneity within patient populations, therapeutic response rates often vary. Here we used the NCI60 panel to identify SNPs associated with paclitaxel sensitivity. Using the panel's GI50 response data available from Developmental Therapeutics Program, cell lines were categorized as either sensitive or resistant. PLINK software was used to perform a genome-wide association analysis of the cellular response to paclitaxel with the panel's SNP-genotype data on the Affymetrix 125 k SNP array. FastSNP software helped predict each SNP's potential impact on their gene product. mRNA expression differences between sensitive and resistant cell lines was examined using data from BioGPS. Using Haploview software, we investigated for haplotypes that were more strongly associated with the cellular response to paclitaxel. Ingenuity Pathway Analysis software helped us understand how our identified genes may alter the cellular response to paclitaxel.

Results

43 SNPs were found significantly associated (FDR < 0.005) with paclitaxel response, with 10 belonging to protein-coding genes (CFTR, ROBO1, PTPRD, BTBD12, DCT, SNTG1, SGCD, LPHN2, GRIK1, ZNF607). SNPs in GRIK1, DCT, SGCD and CFTR were predicted to be intronic enhancers, altering gene expression, while SNPs in ZNF607 and BTBD12 cause conservative missense mutations. mRNA expression analysis supported these findings as GRIK1, DCT, SNTG1, SGCD and CFTR showed significantly (p < 0.05) increased expression among sensitive cell lines. Haplotypes found in GRIK1, SGCD, ROBO1, LPHN2, and PTPRD were more strongly associated with response than their individual SNPs.

Conclusions

Our study has taken advantage of available genotypic data and its integration with drug response data obtained from the NCI60 panel. We identified 10 SNPs located within protein-coding genes that were not previously shown to be associated with paclitaxel response. As only five genes showed differential mRNA expression, the remainder would not have been detected solely based on expression data. The identified haplotypes highlight the role of utilizing SNP combinations within genomic loci of interest to improve the risk determination associated with drug response. These genetic variants represent promising biomarkers for predicting paclitaxel response and may play a significant role in the cellular response to paclitaxel.     

Development of Oral Epithelial Cell Line ROE2 with Differentiation Potential from Transgenic Rats Harboring Temperature-Sensitive Simian Virus40 Large T-Antigen Gene

We have developed an immortalized oral epithelial cell line, ROE2, from fetal transgenic rats harboring temperature-sensitive simian virus 40 large T-antigen gene. The cells grew continuously at either a permissive temperature of 33°C or an intermediate temperature of 37°C. At the nonpermissive temperature of 39°C, on the other hand, growth decreased significantly, and the Sub-G1 phase of the cell cycle increased, indicating that the cells undergo apoptosis at a nonpermissive temperature. Histological and immunocytochemical analyses revealed that ROE2 cells at 37°C had a stratified epithelial-like morphology and expressed cytokeratins Krt4 and Krt13, marker proteins for oral nonkeratinized epithelial cells. Global-scale comprehensive microarray analysis, coupled with bioinformatics tools, demonstrated a significant gene network that was obtained from the upregulated genes. The gene network contained 16 genes, including Cdkn1a, Fos, Krt13, and Prdm1, and was associated mainly with the biological process of skin development in the category of biological functions, organ development. These four genes were validated by quantitative real-time polymerase chain reaction, and the results were nearly consistent with the microarray data. It is therefore anticipated that this cell line will be useful as an in vitro model for studies such as physiological functions, as well as for gene expression in oral epithelial cells.     

Expression of Novel Alzheimer’s Disease Risk Genes in Control and Alzheimer’s Disease Brains

Late onset Alzheimer’s disease (LOAD) etiology is influenced by complex interactions between genetic and environmental risk factors. Large-scale genome wide association studies (GWAS) for LOAD have identified 10 novel risk genes: ABCA7, BIN1, CD2AP, CD33, CLU, CR1, EPHA1, MS4A6A, MS4A6E, and PICALM. We sought to measure the influence of GWAS single nucleotide polymorphisms (SNPs) and gene expression levels on clinical and pathological measures of AD in brain tissue from the parietal lobe of AD cases and age-matched, cognitively normal controls. We found that ABCA7, CD33, and CR1 expression levels were associated with clinical dementia rating (CDR), with higher expression being associated with more advanced cognitive decline. BIN1 expression levels were associated with disease progression, where higher expression was associated with a delayed age at onset. CD33, CLU, and CR1 expression levels were associated with disease status, where elevated expression levels were associated with AD. Additionally, MS4A6A expression levels were associated with Braak tangle and Braak plaque scores, with elevated expression levels being associated with more advanced brain pathology. We failed to detect an association between GWAS SNPs and gene expression levels in our brain series. The minor allele of rs3764650 in ABCA7 is associated with age at onset and disease duration, and the minor allele of rs670139 in MS4A6E was associated with Braak tangle and Braak plaque score. These findings suggest that expression of some GWAS genes, namely ABCA7, BIN1, CD33, CLU, CR1 and the MS4A family, are altered in AD brains.     

Expression levels of the tetratricopeptide repeat protein gene ttc39b covary with carotenoid-based skin colour in cichlid fish

Carotenoid pigments play a major role in animal body colouration, generating strong interest in the genes involved in the metabolic processes that lead from their dietary uptake to their storage in the integument. Here, we used RNA sequencing (RNA-Seq) to test for differentially expressed genes in a taxonomically replicated design using three pairs of related cichlid fish taxa from the genera Tropheus and Aulonocara. Within each pair, taxa differed in terms of red and yellow body colouration, and high‐performance liquid chromatography (HPLC) analyses of skin extracts revealed different carotenoid profiles and concentrations across the studied taxa. Five genes were differentially expressed in all three yellow–red skin contrasts (dhrsx, nlrc3, tcaf2, urah and ttc39b), but only the tetratricopeptide repeat protein-coding gene ttc39b, whose gene product is linked to mammalian lipid metabolism, was consistently expressed more highly in the red skin samples. The RNA-Seq results were confirmed by quantitative PCR. We propose ttc39b as a compelling candidate gene for variation in animal carotenoid colouration. Since differential expression of ttc39b was correlated with the presence/absence of yellow carotenoids in a previous study, we suggest that ttc39b is more likely associated with the concentration of total carotenoids than with the metabolic formation of red carotenoids.     

Chronic Toll-like receptor 4 stimulation in skin induces inflammation,macrophage activation,transforming growth factor beta signature gene expression,and fibrosis

Introduction

The crucial role of innate immunity in the pathogenesis of systemic sclerosis (SSc) is well established, and in the past few years the hypothesis that Toll-like receptor 4 (TLR4) activation induced by endogenous ligands is involved in fibrogenesis has been supported by several studies on skin, liver, and kidney fibrosis. These findings suggest that TLR4 activation can enhance transforming growth factor beta (TGF-β) signaling, providing a potential mechanism for TLR4/Myeloid differentiation factor 88 (MyD88)-dependent fibrosis.

Methods

The expression of TLR4, CD14 and MD2 genes was analyzed by real-time polymerase chain reaction from skin biopsies of 24 patients with diffuse cutaneous SSc. In order to investigate the effects of the chronic skin exposure to endotoxin (Lipopolysaccharide (LPS)) in vivo we examined the expression of inflammation, TGF-β signaling and cellular markers genes by nanostring. We also identified cellular subsets by immunohistochemistry and flow cytometry.

Results

We found that TLR4 and its co-receptors, MD2 and CD14, are over-expressed in lesional skin from patients with diffuse cutaneous SSc, and correlate significantly with progressive or regressive skin disease as assessed by the Delta Modified Rodnan Skin Score. In vivo, a model of chronic dermal LPS exposure showed overexpression of proinflammatory chemokines, recruitment and activation of macrophages, and upregulation of TGF-β signature genes.

Conclusions

We delineated the role of MyD88 as necessary for the induction not only for the early phase of inflammation, but also for pro-fibrotic gene expression via activation of macrophages. Chronic LPS exposure might be a model of early stage of SSc when inflammation and macrophage activation are important pathological features of the disease, supporting a role for innate immune activation in SSc skin fibrosis.     

Changes in the Gene Expression Profiles of the Hypopharyngeal Gland of Worker Honeybees in Association with Worker Behavior and Hormonal Factors

The hypopharyngeal glands (HPGs) of worker honeybees undergo physiological changes along with the age-dependent role change from nursing to foraging: nurse bee HPGs secrete mainly major royal jelly proteins, whereas forager HPGs secrete mainly α-glucosidase III, which converts the sucrose in the nectar into glucose and fructose. We previously identified two other genes, Apis mellifera buffy (Ambuffy) and Apis mellifera matrix metalloproteinase 1 (AmMMP1), with enriched expression in nurse bee and forager HPGs, respectively. In the present study, to clarify the molecular mechanisms that coordinate HPG physiology with worker behavior, we first analyzed whether Ambuffy, AmMMP1, mrjp2 (a gene encoding one of major royal jelly protein isoforms), and Hbg3 (a gene encoding α-glucosidase III) expression, is associated with worker behavior in ''single-cohort colonies'' where workers of almost the same age perform different tasks. Expression of these genes correlated with the worker’s role, while controlling for age, indicating their regulation associated with the worker’s behavior. Associated gene expression suggested the possible involvement of some hormonal factors in its regulation. We therefore examined the relationship between ecdysone- and juvenile hormone (JH)-signaling, and the expression profiles of these ‘indicator’ genes (nurse bee HPG-selective genes: mrjp2 and Ambuffy, and forager HPG-selective genes: Hbg3 and AmMMP1). Expression of both ecdysone-regulated genes (ecdysone receptor, mushroom body large type Kenyon cell specific protein-1, and E74) and JH-regulated genes (Methoprene tolerant and Krüppel homolog 1) was higher in the forager HPGs than in the nurse bee HPGs, suggesting the possible roles of ecdysone- and JH-regulated genes in worker HPGs. Furthermore, 20-hydroxyecdysone-treatment repressed both nurse bee- and forager-selective gene expression, whereas methoprene-treatment enhanced the expression of forager-selective genes and repressed nurse bee-selective genes in the HPGs. Our findings suggest that both ecdysone- and JH-signaling cooperatively regulate the physiological state of the HPGs in association with the worker’s behavior.