Guide for library design and bias correction for large-scale transcriptome studies using highly multiplexed RNAseq methods.

BACKGROUND: Standard RNAseq methods using bulk RNA and recent single-cell RNAseq methods use DNA barcodes to identify samples and cells, and the barcoded cDNAs are pooled into a library pool before high throughput sequencing. In cases of single-cell and low-input RNAseq methods, the library is further amplified by PCR after the pooling. Preparation of hundreds or more samples for a large study often requires multiple library pools. However, sometimes correlation between expression profiles among the libraries is low and batch effect biases make integration of data between library pools difficult. RESULTS: We investigated 166 technical replicates in 14 RNAseq libraries made using the STRT method. The patterns of the library biases differed by genes, and uneven library yields were associated with library biases. The former bias was corrected using the NBGLM-LBC algorithm, which we present in the current study. The latter bias could not be corrected directly, but could be solved by omitting libraries with particularly low yields. A simulation experiment suggested that the library bias correction using NBGLM-LBC requires a consistent sample layout. The NBGLM-LBC correction method was applied to an expression profile for a cohort study of childhood acute respiratory illness, and the library biases were resolved. CONCLUSIONS: The R source code for the library bias correction named NBGLM-LBC is available at https://shka.github.io/NBGLM-LBC and https://shka.bitbucket.io/NBGLM-LBC . This method is applicable to correct the library biases in various studies that use highly multiplexed sequencing-based profiling methods with a consistent sample layout with samples to be compared (e.g., "cases" and "controls") equally distributed in each library.

Intracellular signalling pathways and cytoskeletal functions converge on the psoriasis candidate gene CCHCR1 expressed at P-bodies and centrosomes.

BACKGROUND: CCHCR1 (Coiled-Coil α-Helical Rod protein 1) is a putative psoriasis candidate gene with the risk alleles CCHCR1*WWCC and *Iso3, the latter inhibiting the translation of isoform 1. CCHCR1 was recently shown to be a centrosomal protein, as well as a component of cytoplasmic processing bodies (P-bodies) that regulate mRNA turnover. The function of CCHCR1 has remained unsettled, partly because of the inconsistent findings; it has been shown to play a wide variety of roles in divergent processes, e.g., cell proliferation and steroidogenesis. Here we utilized RNA sequencing (RNAseq) using HEK293 cells overexpressing isoforms 1 or 3 (Iso1, Iso3 cells), in combination with the coding non-risk or risk (*WWCC) haplotype of CCHCR1. Our aim was to study the overall role of CCHCR1 and the effects of its variants. RESULTS: The overexpression of CCHCR1 variants in HEK293 cells resulted in cell line-specific expression profiles though several similarities were observable. Overall the Iso1 and Iso3 cells showed a clear isoform-specific clustering as two separate groups, and the Non-risk and Risk cells often exhibited opposite effects. The RNAseq supported a role for CCHCR1 in the centrosomes and P-bodies; the most highlighted pathways included regulation of cytoskeleton, adherens and tight junctions, mRNA surveillance and RNA transport. Interestingly, both the RNAseq and immunofluorescent localization revealed variant-specific differences for CCHCR1 within the P-bodies. CONCLUSIONS: CCHCR1 influenced a wide variety of signaling pathways, which could reflect its active role in the P-bodies and centrosomes that both are linked to the cytoskeleton; as a centrosomal P-body protein CCHCR1 may regulate diverse cytoskeleton-mediated functions, such as cell adhesion and -division. The present findings may explain the previous inconsistent observations about the functions of CCHCR1.

Gene expression analysis of skin grafts and cultured keratinocytes using synthetic RNA normalization reveals insights into differentiation and growth control.

BACKGROUND: Keratinocytes (KCs) are the most frequent cells in the epidermis, and they are often isolated and cultured in vitro to study the molecular biology of the skin. Cultured primary cells and various immortalized cells have been frequently used as skin models but their comparability to intact skin has been questioned. Moreover, when analyzing KC transcriptomes, fluctuation of polyA+ RNA content during the KCs' lifecycle has been omitted. RESULTS: We performed STRT RNA sequencing on 10 ng samples of total RNA from three different sample types: i) epidermal tissue (split-thickness skin grafts), ii) cultured primary KCs, and iii) HaCaT cell line. We observed significant variation in cellular polyA+ RNA content between tissue and cell culture samples of KCs. The use of synthetic RNAs and SAMstrt in normalization enabled comparison of gene expression levels in the highly heterogenous samples and facilitated discovery of differences between the tissue samples and cultured cells. The transcriptome analysis sensitively revealed genes involved in KC differentiation in skin grafts and cell cycle regulation related genes in cultured KCs and emphasized the fluctuation of transcription factors and non-coding RNAs associated to sample types. CONCLUSIONS: The epidermal keratinocytes derived from tissue and cell culture samples showed highly different polyA+ RNA contents. The use of SAMstrt and synthetic RNA based normalization allowed the comparison between tissue and cell culture samples and thus proved to be valuable tools for RNA-seq analysis with translational approach. Transciptomics revealed clear difference both between tissue and cell culture samples and between primary KCs and immortalized HaCaT cells.

Aberrant splicing of genes involved in haemoglobin synthesis and impaired terminal erythroid maturation in SF3B1 mutated refractory anaemia with ring sideroblasts.

Refractory anaemia with ring sideroblasts (RARS) is distinguished by hyperplastic inefficient erythropoiesis, aberrant mitochondrial ferritin accumulation and anaemia. Heterozygous mutations in the spliceosome gene SF3B1 are found in a majority of RARS cases. To explore the link between SF3B1 mutations and anaemia, we studied mutated RARS CD34(+) marrow cells with regard to transcriptome sequencing, splice patterns and mutational allele burden during erythroid differentiation. Transcriptome profiling during early erythroid differentiation revealed a marked up-regulation of genes involved in haemoglobin synthesis and in the oxidative phosphorylation process, and down-regulation of mitochondrial ABC transporters compared to normal bone marrow. Moreover, mis-splicing of genes involved in transcription regulation, particularly haemoglobin synthesis, was confirmed, indicating a compromised haemoglobinization during RARS erythropoiesis. In order to define the phase during which erythroid maturation of SF3B1 mutated cells is most affected, we assessed allele burden during erythroid differentiation in vitro and in vivo and found that SF3B1 mutated erythroblasts showed stable expansion until late erythroblast stage but that terminal maturation to reticulocytes was significantly reduced. In conclusion, SF3B1 mutated RARS progenitors display impaired splicing with potential downstream consequences for genes of key importance for haemoglobin synthesis and terminal erythroid differentiation.

Optimized single-cell RNA sequencing protocol to study early genome activation in mammalian preimplantation development.

Here, we present a modification of single-cell tagged reverse transcription protocol to study gene expression on a single-cell level or with limited RNA input. We describe different enzymes for reverse transcription and cDNA amplification, modified lysis buffer, and additional clean-up steps before cDNA amplification. We also detail an optimized single-cell RNA sequencing method for handpicked single cells, or tens to hundreds of cells, as input material to study mammalian preimplantation development. For complete details on the use and execution of this protocol, please refer to Ezer et al.1.

Influence of FLG loss-of-function mutations in host-microbe interactions during atopic skin inflammation.

BACKGROUND: Loss-of-function mutations in the filaggrin (FLG) gene directly alter skin barrier function and critically influence atopic inflammation. While skin barrier dysfunction, Th2-associated inflammation and bacterial dysbiosis are well-known characteristics of atopic dermatitis (AD), the mechanisms interconnecting genotype, transcriptome and microbiome remain largely elusive. OBJECTIVE: In-depth analysis of FLG genotype-associated skin gene expression alterations and host-microbe interactions in AD. METHODS: Multi-omics characterization of a cohort of AD patients carrying heterozygous loss-of-function mutations in the FLG gene (ADMut) (n = 15), along with matched wild-type (ADWt) patients and healthy controls. Detailed clinical characterization, microarray gene expression and 16 S rRNA-based microbial marker gene data were generated and analyzed. RESULTS: In the context of filaggrin dysfunction, the transcriptome was characterized by dysregulation of barrier function and water homeostasis, while the lesional skin of ADWt demonstrated the specific upregulation of pro-inflammatory cytokines and T-cell proliferation. S. aureus dominated the microbiome in both patient groups, however, shifting microbial communities could be observed when comparing healthy with non-lesional ADWt or ADMut skin, offering the opportunity to identify microbe-associated transcriptomic signatures. Moreover, an AD core signature with 28 genes, including CCL13, CCL18, BTC, SCIN, RAB31 and PCLO was identified. CONCLUSIONS: Our integrative approach provides molecular insights for the concept that FLG loss-of-function mutations are a genetic shortcut to atopic inflammation and unravels the complex interplay between genotype, transcriptome and microbiome in the human holobiont.

Biomarkers of nanomaterials hazard from multi-layer data.

There is an urgent need to apply effective, data-driven approaches to reliably predict engineered nanomaterial (ENM) toxicity. Here we introduce a predictive computational framework based on the molecular and phenotypic effects of a large panel of ENMs across multiple in vitro and in vivo models. Our methodology allows for the grouping of ENMs based on multi-omics approaches combined with robust toxicity tests. Importantly, we identify mRNA-based toxicity markers and extensively replicate them in multiple independent datasets. We find that models based on combinations of omics-derived features and material intrinsic properties display significantly improved predictive accuracy as compared to physicochemical properties alone.

Toxicogenomic Profiling of 28 Nanomaterials in Mouse Airways.

Toxicogenomics opens novel opportunities for hazard assessment by utilizing computational methods to map molecular events and biological processes. In this study, the transcriptomic and immunopathological changes associated with airway exposure to a total of 28 engineered nanomaterials (ENM) are investigated. The ENM are selected to have different core (Ag, Au, TiO2, CuO, nanodiamond, and multiwalled carbon nanotubes) and surface chemistries (COOH, NH2, or polyethylene glycosylation (PEG)). Additionally, ENM with variations in either size (Au) or shape (TiO2) are included. Mice are exposed to 10 µg of ENM by oropharyngeal aspiration for 4 consecutive days, followed by extensive histological/cytological analyses and transcriptomic characterization of lung tissue. The results demonstrate that transcriptomic alterations are correlated with the inflammatory cell infiltrate in the lungs. Surface modification has varying effects on the airways with amination rendering the strongest inflammatory response, while PEGylation suppresses toxicity. However, toxicological responses are also dependent on ENM core chemistry. In addition to ENM-specific transcriptional changes, a subset of 50 shared differentially expressed genes is also highlighted that cluster these ENM according to their toxicity. This study provides the largest in vivo data set currently available and as such provides valuable information to be utilized in developing predictive models for ENM toxicity.

Fetal HLA-G mediated immune tolerance and interferon response in preeclampsia.

BACKGROUND: Fetal immune tolerance is crucial for pregnancy success. We studied the link between preeclampsia, a severe pregnancy disorder with uncertain pathogenesis, and fetal human leukocyte antigen G (HLA-G) and other genes regulating maternal immune responses. METHODS: We assessed sex ratios and regulatory HLA-G haplotypes in population cohorts and series of preeclampsia and stillbirth. We studied placental mRNA expression of 136 genes by sequencing and HLA-G and interferon alpha (IFNα) protein expression by immunohistochemistry. FINDINGS: We found underrepresentation of males in preeclamptic births, especially those delivered preterm or small for gestational age. Balancing selection at HLA-G associated with the sex ratio, stillbirth, and preeclampsia. We observed downregulation of HLA-G, its receptors, and many other tolerogenic genes, and marked upregulation of IFNA1 in preeclamptic placentas. INTERPRETATION: These findings indicate that an evolutionary trade-off between immune tolerance and protection against infections at the maternal-fetal interface promotes genetic diversity in fetal HLA-G, thereby affecting survival, preeclampsia, and sex ratio. We highlight IFNA1 as a potential mediator of preeclampsia and a target for therapeutic trials. FUNDING: Finnish Medical Foundation, Päivikki and Sakari Sohlberg Foundation, Karolinska Institutet Research Foundation, Scandinavia-Japan Sasakawa Foundation, Japan Eye Bank Association, Astellas Foundation for Research on Metabolic Disorders, Japan Society for the Promotion of Science, Knut and Alice Wallenberg Foundation, Swedish Research Council, Medical Society Liv och Hälsa, Sigrid Jusélius Foundation, Helsinki University Hospital and University of Helsinki, Jane and Aatos Erkko Foundation, Academy of Finland, Finska Läkaresällskapet, Novo Nordisk Foundation, Finnish Foundation for Pediatric Research, and Emil Aaltonen Foundation.

Multiparametric Profiling of Engineered Nanomaterials: Unmasking the Surface Coating Effect.

Despite considerable efforts, the properties that drive the cytotoxicity of engineered nanomaterials (ENMs) remain poorly understood. Here, the authors inverstigate a panel of 31 ENMs with different core chemistries and a variety of surface modifications using conventional in vitro assays coupled with omics-based approaches. Cytotoxicity screening and multiplex-based cytokine profiling reveals a good concordance between primary human monocyte-derived macrophages and the human monocyte-like cell line THP-1. Proteomics analysis following a low-dose exposure of cells suggests a nonspecific stress response to ENMs, while microarray-based profiling reveals significant changes in gene expression as a function of both surface modification and core chemistry. Pathway analysis highlights that the ENMs with cationic surfaces that are shown to elicit cytotoxicity downregulated DNA replication and cell cycle responses, while inflammatory responses are upregulated. These findings are validated using cell-based assays. Notably, certain small, PEGylated ENMs are found to be noncytotoxic yet they induce transcriptional responses reminiscent of viruses. In sum, using a multiparametric approach, it is shown that surface chemistry is a key determinant of cellular responses to ENMs. The data also reveal that cytotoxicity, determined by conventional in vitro assays, does not necessarily correlate with transcriptional effects of ENMs.

MMP-10 (Stromelysin-2) and MMP-21 in human and murine squamous cell cancer.

The squamous cell cancers (SCC) of renal transplant recipients are more aggressive and metastasize earlier than those of the non-immunocompromised population. Matrix metalloproteinases (MMPs) have a central role in tumor initiation, invasion and metastasis. Our aim was to compare the expression of MMPs-10, -12 and -21 in SCCs from immunosuppressed (IS) and control patients and the contribution of MMPs-10 and -21 to SCC development in the FVB/N-Tg(KRT5-Nfkbia)3Rto mouse line. Immunohistochemical analysis of 25 matched pairs of SCCs, nine of Bowen's disease and timed back skin biopsies of mice with selective inhibition of Rel/NF-kappaB signalling were performed. Semiquantitatively assessed stromal MMP-10 expression was higher (P = 0.009) in the control group when compared with IS patients. Tumor cell-derived MMP-10, -12 and -21 expression did not differ between the groups but stromal fibroblasts of the control SCCs tended to express MMP-21 more abundantly. MMP-10 expression was observed already in Bowen's disease while MMP-21 was absent. MMP-10 and -21 were present in inflammatory or stromal cells in ageing mice while dysplastic keratinocytes and invasive cancer were negative. Our results suggest that MMP-10 may be important in the initial stages of SCC progression and induced in the stroma relating to the general host-response reaction to skin cancer. MMP-21 does not associate with invasion of SCC but may be involved in keratinocyte differentiation.

The human long non-coding RNA gene RMRP has pleiotropic effects and regulates cell-cycle progression at G2.

RMRP was the first non-coding nuclear RNA gene implicated in a disease. Its mutations cause cartilage-hair hypoplasia (CHH), an autosomal recessive skeletal dysplasia with growth failure, immunodeficiency, and a high risk for malignancies. This study aimed to gain further insight into the role of RNA Component of Mitochondrial RNA Processing Endoribonuclease (RMRP) in cellular physiology and disease pathogenesis. We combined transcriptome analysis with single-cell analysis using fibroblasts from CHH patients and healthy controls. To directly assess cell cycle progression, we followed CHH fibroblasts by pulse-labeling and time-lapse microscopy. Transcriptome analysis identified 35 significantly upregulated and 130 downregulated genes in CHH fibroblasts. The downregulated genes were significantly connected to the cell cycle. Multiple other pathways, involving regulation of apoptosis, bone and cartilage formation, and lymphocyte function, were also affected, as well as PI3K-Akt signaling. Cell-cycle studies indicated that the CHH cells were delayed specifically in the passage from G2 phase to mitosis. Our findings expand the mechanistic understanding of CHH, indicate possible pathways for therapeutic intervention and add to the limited understanding of the functions of RMRP.

Cationic gold nanoparticles elicit mitochondrial dysfunction: a multi-omics study.

Systems biology is increasingly being applied in nanosafety research for observing and predicting the biological perturbations inflicted by exposure to nanoparticles (NPs). In the present study, we used a combined transcriptomics and proteomics approach to assess the responses of human monocytic cells to Au-NPs of two different sizes with three different surface functional groups, i.e., alkyl ammonium bromide, alkyl sodium carboxylate, or poly(ethylene glycol) (PEG)-terminated Au-NPs. Cytotoxicity screening using THP-1 cells revealed a pronounced cytotoxicity for the ammonium-terminated Au-NPs, while no cell death was seen after exposure to the carboxylated or PEG-modified Au-NPs. Moreover, Au-NR3+ NPs, but not the Au-COOH NPs, were found to trigger dose-dependent lethality in vivo in the model organism, Caenorhabditis elegans. RNA sequencing combined with mass spectrometry-based proteomics predicted that the ammonium-modified Au-NPs elicited mitochondrial dysfunction. The latter results were validated by using an array of assays to monitor mitochondrial function. Au-NR3+ NPs were localized in mitochondria of THP-1 cells. Moreover, the cationic Au-NPs triggered autophagy in macrophage-like RFP-GFP-LC3 reporter cells, and cell death was aggravated upon inhibition of autophagy. Taken together, these studies have disclosed mitochondria-dependent effects of cationic Au-NPs resulting in the rapid demise of the cells.

Microbe-host interplay in atopic dermatitis and psoriasis.

Despite recent advances in understanding microbial diversity in skin homeostasis, the relevance of microbial dysbiosis in inflammatory disease is poorly understood. Here we perform a comparative analysis of skin microbial communities coupled to global patterns of cutaneous gene expression in patients with atopic dermatitis or psoriasis. The skin microbiota is analysed by 16S amplicon or whole genome sequencing and the skin transcriptome by microarrays, followed by integration of the data layers. We find that atopic dermatitis and psoriasis can be classified by distinct microbes, which differ from healthy volunteers microbiome composition. Atopic dermatitis is dominated by a single microbe (Staphylococcus aureus), and associated with a disease relevant host transcriptomic signature enriched for skin barrier function, tryptophan metabolism and immune activation. In contrast, psoriasis is characterized by co-occurring communities of microbes with weak associations with disease related gene expression. Our work provides a basis for biomarker discovery and targeted therapies in skin dysbiosis.

The human GIMAP5 gene has a common polyadenylation polymorphism increasing risk to systemic lupus erythematosus.

BACKGROUND: Several members of the GIMAP gene family have been suggested as being involved in different aspects of the immune system in different species. Recently, a mutation in the GIMAP5 gene was shown to cause lymphopenia in a rat model of autoimmune insulin-dependent diabetes. Thus it was hypothesised that genetic variation in GIMAP5 may be involved in susceptibility to other autoimmune disorders where lymphopenia is a key feature, such as systemic lupus erythematosus (SLE). MATERIAL AND METHODS: To investigate this, seven single nucleotide polymorphisms in GIMAP5 were analysed in five independent sets of family-based SLE collections, containing more than 2000 samples. RESULT: A significant increase in SLE risk associated with the most common GIMAP5 haplotype was found (OR 1.26, 95% CI 1.02 to 1.54, p = 0.0033). In families with probands diagnosed with trombocytopenia, the risk was increased (OR 2.11, 95% CI 1.09 to 4.09, p = 0.0153). The risk haplotype bears a polymorphic polyadenylation signal which alters the 3' part of GIMAP5 mRNA by producing an inefficient polyadenylation signal. This results in higher proportion of non-terminated mRNA for homozygous individuals (p<0.005), a mechanism shown to be causal in thalassaemias. To further assess the functional effect of the polymorphic polyadenylation signal in the risk haplotype, monocytes were treated with several cytokines affecting apoptosis. All the apoptotic cytokines induced GIMAP5 expression in two monocyte cell lines (1.5-6 times, p<0.0001 for all tests). CONCLUSION: Taken together, the data suggest the role of GIMAP5 in the pathogenesis of SLE.

Sequence analysis of pooled bacterial samples enables identification of strain variation in group A streptococcus.

Knowledge of the genomic variation among different strains of a pathogenic microbial species can help in selecting optimal candidates for diagnostic assays and vaccine development. Pooled sequencing (Pool-seq) is a cost effective approach for population level genetic studies that require large numbers of samples such as various strains of a microbe. To test the use of Pool-seq in identifying variation, we pooled DNA of 100 Streptococcus pyogenes strains of different emm types in two pools, each containing 50 strains. We used four variant calling tools (Freebayes, UnifiedGenotyper, SNVer, and SAMtools) and one emm1 strain, SF370, as a reference genome. In total 63719 SNPs and 164 INDELs were identified in the two pools concordantly by at least two of the tools. Majority of the variants (93.4%) from six individually sequenced strains used in the pools could be identified from the two pools and 72.3% and 97.4% of the variants in the pools could be mined from the analysis of the 44 complete Str. pyogenes genomes and 3407 sequence runs deposited in the European Nucleotide Archive respectively. We conclude that DNA sequencing of pooled samples of large numbers of bacterial strains is a robust, rapid and cost-efficient way to discover sequence variation.

Investigation of rare and low-frequency variants using high-throughput sequencing with pooled DNA samples.

High-throughput sequencing using pooled DNA samples can facilitate genome-wide studies on rare and low-frequency variants in a large population. Some major questions concerning the pooling sequencing strategy are whether rare and low-frequency variants can be detected reliably, and whether estimated minor allele frequencies (MAFs) can represent the actual values obtained from individually genotyped samples. In this study, we evaluated MAF estimates using three variant detection tools with two sets of pooled whole exome sequencing (WES) and one set of pooled whole genome sequencing (WGS) data. Both GATK and Freebayes displayed high sensitivity, specificity and accuracy when detecting rare or low-frequency variants. For the WGS study, 56% of the low-frequency variants in Illumina array have identical MAFs and 26% have one allele difference between sequencing and individual genotyping data. The MAF estimates from WGS correlated well (r = 0.94) with those from Illumina arrays. The MAFs from the pooled WES data also showed high concordance (r = 0.88) with those from the individual genotyping data. In conclusion, the MAFs estimated from pooled DNA sequencing data reflect the MAFs in individually genotyped samples well. The pooling strategy can thus be a rapid and cost-effective approach for the initial screening in large-scale association studies.

NOD-like receptor signaling and inflammasome-related pathways are highlighted in psoriatic epidermis.

Psoriatic skin differs distinctly from normal skin by its thickened epidermis. Most gene expression comparisons utilize full-thickness biopsies, with substantial amount of dermis. We assayed the transcriptomes of normal, lesional, and non-lesional psoriatic epidermis, sampled as split-thickness skin grafts, with 5'-end RNA sequencing. We found that psoriatic epidermis contains more mRNA per total RNA than controls, and took this into account in the bioinformatic analysis. The approach highlighted innate immunity-related pathways in psoriasis, including NOD-like receptor (NLR) signaling and inflammasome activation. We demonstrated that the NLR signaling genes NOD2, PYCARD, CARD6, and IFI16 are upregulated in psoriatic epidermis, and strengthened these findings by protein expression. Interestingly, PYCARD, the key component of the inflammasome, showed an altered expression pattern in the lesional epidermis. The profiling of non-lesional skin highlighted PSORS4 and mitochondrially encoded transcripts, suggesting that their gene expression is altered already before the development of lesions. Our data suggest that all components needed for the active inflammasome are present in the keratinocytes of psoriatic skin. The characterization of inflammasome pathways provides further opportunities for therapy. Complementing previous transcriptome studies, our approach gives deeper insight into the gene regulation in psoriatic epidermis.

Matrilysin-2 (matrix metalloproteinase-26) is upregulated in keratinocytes during wound repair and early skin carcinogenesis.

Matrilysin-2 (matrix metalloproteinase (MMP)-26) is a small protein of the MMP family expressed in some epithelial carcinomas and normal tissues. We studied its role in benign skin disorders characterized by epithelial proliferation, in wound repair, skin cancer, and regulation in keratinocyte (KC) cultures. MMP-26 is expressed by laminin-5-positive KC in the migrating area during wound repair, in benign skin disorders characterized by inflammation and microdisruptions of basement membrane, but in intact skin only in hair follicles. It was detected in occasional atypical KC in pre-malignant lesions but not in basal cell cancer islands. Although MMP-26 was expressed in grades I and II squamous cell cancers (SCC), it was not present in dedifferentiated grade III tumors. MMP-26 was neither co-expressed with its close homologue matrilysin-1 nor with the proliferation marker Ki-67. But in tissue samples it either co-localized or was detected in adjacent cells of same regions with the tumor suppressor p16. In KC and HaCaT cell cultures, 12-phorbol-13-myristate-acetate, epidermal growth factor, tumor necrosis factor-alpha, transforming growth factor-beta1, interleukin-1 (IL-1)beta, IL-6, insulin-like growth factor, gamma-IFN, retinoic acid, dexamethasone, four matrices or ras-transformation were unable to upregulate MMP-26 expression. The expression pattern of MMP-26 suggests that it may be upregulated in basal KC even without tumorigenesis because of altered cell-matrix interactions and inflammation and, unlike most MMP, becomes downregulated during histological dedifferentiation of SCC. Thus, lack of MMP-26 in SCC could be a marker of aggressive growth.