Detecting Embryo Developmental Potential by Single Blastomere RNA-Seq.

Recent advances in preimplantation embryo diagnostics enable a wide range of applications using single cell biopsy and molecular-based selection techniques without compromising embryo production. This study was conducted to develop a single cell embryo biopsy technique and gene expression analysis method with a very low input volume to ensure normal embryo development and to see if there are differences in gene expression profiles between day-5 biopsied bovine embryos that developed into blastocysts and embryos arrested at morula stage. Out of the 65 biopsied morulae, 32 developed to blastocysts (49.2%). Out of the 13,580 successfully annotated genes, 1204 showed a difference in mRNA expression level. Out of these, 155 genes were expressed in embryos developing to blastocysts. The pathway enrichment analysis revealed significant enrichment in "organelle biogenesis and maintenance", "mRNA splicing" and "mitochondrial translation" pathways. These findings suggest principal differences in gene expression patterns and functional networks of embryos able to reach the blastocyst stage compared to embryos arrested in development. Our preliminary data suggest that single blastomere biopsy and selected gene expression profiles at morula stage could offer additional possibilities for early preimplantation embryo selection before transfer.

Phenotypic Variation in Vietnamese Osteogenesis Imperfecta Patients Sharing a Recessive P3H1 Pathogenic Variant.

Osteogenesis imperfecta (OI) is a syndromic disorder of bone fragility with high variation in its clinical presentation. Equally variable is molecular aetiology; recessive forms are caused by approximately 20 different genes, many of which are directly implicated in collagen type I biosynthesis. Biallelic variants in prolyl 3-hydroxylase 1 (P3H1) are known to cause severe OI by affecting the competence of the prolyl 3-hydroxylation­cartilage associated protein­peptidyl-prolyl cis-trans isomerase B (P3H1-CRTAP-CyPB) complex, which acts on the Pro986 residue of collagen type I α 1 (COL1A1) and Pro707 collagen type I α 2 (COL1A2) chains. The investigation of an OI cohort of 146 patients in Vietnam identified 14 families with P3H1 variants. The c.1170+5G>C variant was found to be very prevalent (12/14) and accounted for 10.3% of the Vietnamese OI cohort. New P3H1 variants were also identified in this population. Interestingly, the c.1170+5G>C variants were found in families with the severe clinical Sillence types 2 and 3 but also the milder types 1 and 4. This is the first time that OI type 1 is reported in patients with P3H1 variants expanding the clinical spectrum. Patients with a homozygous c.1170+5G>C variant shared severe progressively deforming OI type 3: bowed long bones, deformities of ribcage, long phalanges and hands, bluish sclera, brachycephaly, and early intrauterine fractures. Although it remains unclear if the c.1170+5G>C variant constitutes a founder mutation in the Vietnamese population, its prevalence makes it valuable for the molecular diagnosis of OI in patients of the Kinh ethnicity. Our study provides insight into the clinical and genetic variation of P3H1-related OI in the Vietnamese population.

The Expression Pattern of Genes Related to Melanogenesis and Endogenous Opioids in Psoriasis.

The melanocortin system is a major regulator of stress responses in the skin and is responsible for the induction of melanin synthesis through activation of melanogenesis enzymes. The expression of both melanocortin system genes and melanogenesis enzyme genes is altered in psoriasis, and the focus here was on twelve genes related to the signal transduction between them. Additionally, five endogenous opioid system genes that are involved in cutaneous inflammation were examined. Quantitative real-time-PCR was utilized to measure mRNA expression in punch biopsies from lesional and non-lesional skin of psoriasis patients and from the skin of healthy control subjects. Most of the genes related to melanogenesis were down-regulated in patients (CREB1, MITF, LEF1, USF1, MAPK14, ICAM1, PIK3CB, RPS6KB1, KIT, and ATRN). Conversely, an up-regulation occurred in the case of opioids (PENK, PDYN, and PNOC). The suppression of genes related to melanogenesis is in agreement with the reported reduction in pigmentation signaling in psoriatic skin and potentially results from the pro-inflammatory environment. The increase in endogenous opioids can be associated with their involvement in inflammatory dysregulation in psoriasis.

Whole genome and exome sequencing reference datasets from a multi-center and cross-platform benchmark study.

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing.

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

Establishing community reference samples, data and call sets for benchmarking cancer mutation detection using whole-genome sequencing.

The lack of samples for generating standardized DNA datasets for setting up a sequencing pipeline or benchmarking the performance of different algorithms limits the implementation and uptake of cancer genomics. Here, we describe reference call sets obtained from paired tumor-normal genomic DNA (gDNA) samples derived from a breast cancer cell line-which is highly heterogeneous, with an aneuploid genome, and enriched in somatic alterations-and a matched lymphoblastoid cell line. We partially validated both somatic mutations and germline variants in these call sets via whole-exome sequencing (WES) with different sequencing platforms and targeted sequencing with >2,000-fold coverage, spanning 82% of genomic regions with high confidence. Although the gDNA reference samples are not representative of primary cancer cells from a clinical sample, when setting up a sequencing pipeline, they not only minimize potential biases from technologies, assays and informatics but also provide a unique resource for benchmarking 'tumor-only' or 'matched tumor-normal' analyses.

The Genetic Variations Associated With Time to Aseptic Loosening After Total Joint Arthroplasty.

BACKGROUND: Total joint arthroplasty (TJA) is one of the most frequent surgical procedures performed in modern hospitals, and aseptic loosening is the most common indication for revision surgeries. We conducted a systemic exploration of potential genetic determinants for early aseptic loosening. METHODS: Data from 423 patients undergoing TJA were collected and analyzed. Three analytical groups were formed based on joint arthroplasty status. Group 1 were TJA patients without symptoms of aseptic loosening of at least 1 year, group 2 were patients with primary TJA, and group 3 were patients receiving revision surgery because of aseptic loosening. Genome-wide genotyping comparing genotype frequencies between patients with and without aseptic loosening (group 3 vs groups 1 and 2) was conducted. A case-control association analysis and linear modeling were applied to identify the impact of the identified genes on implant survival with time to the revision as an outcome measure. RESULTS: We identified 52 single-nucleotide polymorphisms (SNPs) with a genome-wide suggestive P value less than 10-5 to be associated with the implant loosening. The most remarkable odds ratios (OR) were found with the variations in the IFIT2/IFIT3 (OR, 21.6), CERK (OR, 12.6), and PAPPA (OR, 14.0) genes. Variations in the genotypes of 4 SNPs-rs115871127, rs16823835, rs13275667, and rs2514486-predicted variability in the time to aseptic loosening. The time to aseptic loosening varied from 8 to 16 years depending on the genotype, indicating a substantial effect of genetic variance. CONCLUSION: Development of the aseptic loosening is associated with several genetic variations and we identified at least 4 SNPs with a significant effect on the time for loosening. These data could help to develop a personalized approach for TJA and loosening management.

Identification of an optimal method for extracting RNA from human skin biopsy, using domestic pig as a model system.

To evaluate skin tissue gene expression patterns correctly, extracting sufficient quantities of good quality RNA is essential. However, RNA extraction from skin tissue is challenging, as the hyaluronic acid-collagen matrix is extremely difficult to homogenize. Although there are multiple ways to extract RNA from skin, there are no comparative studies that identify the most critical steps, e.g. sample collection, storage and homogenization. We analysed the various steps involved in RNA extraction (i.e. biopsy collection as dry biopsy or into nucleotide stabilizing reagents, different storage conditions, enzymatic digestion, stator-rotor and bead motion-based homogenizing combined with column-based RNA purification). We hypothesised that domestic pig skin is applicable as a model for human skin studies. Altogether twenty different workflows were tested on pig skin and the four most promising workflows were tested on human skin samples. The optimal strategy for extracting human skin RNA was to collect, store and homogenize the sample in RLT lysis buffer from the RNeasy Fibrous Tissue Kit combined with beta-mercaptoethanol. Both stator-rotor and bead motion-based homogenizing were found to result in high quality and quantity of extracted RNA. Our results confirmed that domestic pig skin can be successfully used as a model for human skin RNA studies.

Multicomponent Biomarker Approach Improves the Accuracy of Diagnostic Biomarkers for Psoriasis Vulgaris.

Accurate biomarker-based diagnosis of psoriasis vulgaris has remained a challenge; no reliable disease-specific biomarkers have yet been identified. There are several different chronic inflammatory skin diseases that can present similar clinical and dermoscopy features to psoriasis vulgaris, making accurate diagnosis more difficult. Both literature-based and data-driven selection of biomarker was conducted to select candidates for a multicomponent biomarker for psoriasis vulgaris. Support vector machine-based classification models were trained using gene expression data from locally recruited patients and validated on 7 public datasets, which included gene expression data of other inflammatory skin diseases in addition to psoriasis vulgaris. The resulting accuracy of the best classification model based on the expression levels of 4 genes (IL36G, CCL27, NOS2 and C10orf99) was 96.4%, outperforming classification based on other marker gene combinations, which were more affected by variability in gene expression profiles between different datasets and patient groups. This approach has the potential to fill the void of clinically applicable diagnostic biomarkers for psoriasis vulgaris and other inflammatory skin diseases.

Serum Amyloid Alpha Is Downregulated in Peripheral Tissues of Parkinson's Disease Patients.

We report the changed levels of serum amyloid alpha, an immunologically active protein, in Parkinson's disease (PD) patients' peripheral tissues. We have previously shown that Saa-1 and -2 (serum amyloid alpha-1,-2, genes) were among the top downregulated genes in PD patients' skin, using whole-genome RNA sequencing. In the current study, we characterized the gene and protein expression profiles of skin and blood samples from patients with confirmed PD diagnosis and age/sex matched controls. qRT-PCR analysis of PD skin demonstrated downregulation of Saa-1 and -2 genes in PD patients. However, the lowered amount of protein could not be visualized using immunohistochemistry, due to low quantity of SAA (Serum Amyloid Alpha, protein) in skin. Saa-1 and -2 expression levels in whole blood were below detection threshold based on RNA sequencing, however significantly lowered protein levels of SAA1/2 in PD patients' serum were shown with ELISA, implying that SAA is secreted into the blood. These results show that SAA is differentially expressed in the peripheral tissues of PD patients.

Enhanced Expression of Genes Related to Xenobiotic Metabolism in the Skin of Patients with Atopic Dermatitis but Not with Ichthyosis Vulgaris.

Previous transcriptome analyses underscored the importance of immunological and skin barrier abnormalities in atopic dermatitis (AD). We sought to identify pathogenic pathways involved in AD by comparing the transcriptomes of AD patients stratified for filaggrin (FLG)-null mutations to those of both healthy donors and patients with ichthyosis vulgaris. We applied RNA sequencing to analyze the whole transcriptome of nonlesional skin. We found that 607 genes (476 up-regulated and 131 down-regulated by >2-fold) and 193 genes (172 up-regulated and 21 down-regulated by >2-fold) were differentially expressed when all AD or ichthyosis vulgaris patients were compared with healthy donors, respectively. Expression of genes involved in RNA/protein turnover and adenosine triphosphate synthesis, as well as genes involved in cell death, response to oxidative stress, DNA damage/repair, and autophagy, were significantly enriched in AD skin and, to a lesser extent, in ichthyosis vulgaris skin. FLG-null mutations appear to hardly interfere with current observations. Genes related to xenobiotic metabolism were up-regulated in AD skin only, as were genes related to arachidonic, linoleic, and α-linolenic acid metabolism. Thus, this work newly links AD pathogenesis to aberrant expression of genes related to xenobiotic metabolism.

Analysis of the Expression of Repetitive DNA Elements in Osteosarcoma.

Osteosarcoma (OS) is a rare malignant bone tumor. It affects mostly young persons and has poor outcome with the present treatment. No improvement was observed since the introduction of chemotherapy. The better understanding of osteosarcoma development could indicate better management strategy. Repetitive DNA elements were found to play a role in cancer mechanism especially in epithelial tumors but not yet analyzed in osteosarcoma. We conducted the study to analyse the expression profile of repetitive elements (RE) in osteosarcoma. Methods: Fresh bone paired (tumor and normal bone) samples were obtained from excised parts of tumors of 18 patients with osteosarcoma. We performed sequencing of RNA extracted from 36 samples (18 tumor tissues and 18 normal bone for controls), mapped raw reads to the human genome and identified the REs. EdgeR package was used to analyse the difference in expression of REs between osteosarcoma and normal bone. Results: 82 REs were found differentially expressed (FDR < 0.05) between osteosarcoma and normal bone. Out of all significantly changed REs, 35 were upregulated and 47 were downregulated. HERVs (THE1C-int, LTR5, MER57F and MER87B) and satellite elements (HSATII, ALR-alpha) were the most significantly differential expressed elements between osteosarcoma and normal tissues. These results suggest significant impact of REs in the osteosarcoma. The role of REs should be further studied to understand the mechanism they have in the genesis of osteosarcoma.

Whole transcriptome analysis identifies differentially regulated networks between osteosarcoma and normal bone samples.

We performed whole transcriptome analysis of osteosarcoma bone samples. Initially, we sequenced total RNA from 36 fresh-frozen samples (18 tumoral bone samples and 18 non-tumoral paired samples) matching in pairs for each osteosarcoma patient. We also performed independent gene expression analysis of formalin-fixed paraffin-embedded samples to verify the RNAseq results. Formalin-fixed paraffin-embedded samples allowed us to analyze the effect of chemotherapy. Data were analyzed with DESeq2, edgeR and Reactome packages of R. We found 5365 genes expressed differentially between the normal bone and osteosarcoma tissues with an FDR below 0.05, of which 3399 genes were upregulated and 1966 were downregulated. Among those genes, BTNL9, MMP14, ABCA10, ACACB, COL11A1, and PKM2 were expressed differentially with the highest significance between tumor and normal bone. Functional annotation with the reactome identified significant changes in the pathways related to the extracellular matrix degradation and collagen biosynthesis. It was suggested that chemotherapy may induce the modification of ECM with important collagen biosynthesis. Taken together, our results indicate that changes in the degradation of extracellular matrix seem to be an important mechanism of osteosarcoma and efficient chemotherapy induces the genes related to bone formation. Impact statement Osteosarcoma is a rare disease but it is of interest to many scientists all over the world because the current standard treatment still has poor results. We sequenced total RNA from 36 fresh-frozen paired samples (18 tumoral bone samples and 18 non-tumoral paired samples) from osteosarcoma patients. We found that differences in the gene expressions between the normal and affected bones reflected the changes in the regulation of the degradation of collagen and extracellular matrix. We believe that these findings contribute to the understanding of OS and suggest ideas for further studies.

Looking beyond the brain to improve the pathogenic understanding of Parkinson's disease: implications of whole transcriptome profiling of Patients' skin.

BACKGROUND: Parkinson's Disease is a progressive neurodegenerative disease, characterized by symptoms of motor impairment, resulting from the loss of dopaminergic neurons in the midbrain, however non-neuronal symptoms are also common. Although great advances have been made in the pathogenic understanding of Parkinson's Disease in the nervous system, little is known about the molecular alterations occurring in other non-neuronal organ systems. In addition, a higher rate of melanoma and non-melanoma skin cancer has been observed in the Parkinson's Disease population, indicating crosstalk between these diseases. METHODS: To understand the molecular pathogenesis and gene expression alterations of Parkinson's Disease in peripheral tissues, and in order to explore the possible link between skin cancer and neurodegeneration, whole transcriptomic profiling of patients' skin was performed. Skin biopsies from 12 patients and matched controls were collected, and processed with high-throughput RNA-sequencing analysis. RESULTS: This analysis resulted in a large collection of over 1000 differentially expressed genes, among which clear biological and functional networks could be distinguished. The central functional processes altered in patients skin can be grouped into six broad categories: impaired cellular metabolism and mitochondrial dysfunction, defective protein metabolism, disturbed skin homeostasis, dysfunctional nuclear processes, altered signalling and tumour pathways, as well as disordered immune regulation. CONCLUSIONS: These results demonstrate that the molecular alterations leading to neurodegeneration in the CNS are systemic and manifest also in peripheral tissues, thereby indicating the presence of "skin-brain" crosstalk in Parkinson's Disease. In addition, the extensive homeostatic imbalance and basal stress can lead to increased susceptibility to external and internal mutagenic hazards in these patients, and thus provide a possible molecular link for the crosstalk between skin cancer and Parkinson's Disease.

The clinical features of osteogenesis imperfecta in Vietnam.

PURPOSE: Osteogenesis imperfecta (OI) has not been studied in a Vietnamese population before. The aim of this study was to systematically collect epidemiological information, investigate clinical features and create a clinical database of OI patients in Vietnam for future research and treatment strategy development. METHOD: Participants underwent clinical and physical examinations; also medical records were reviewed. Genealogical information was collected and family members' phenotypical manifestations recorded. Cases were classified according to the Sillence classification. RESULTS: In total, 146 OI patients from 120 families were studied: 46 with OI Type I, 46 with Type III and 54 with Type IV. Almost patients had skeletal deformations. One hundred and forty-two had a history of fractures, 117 blue sclera, 89 dentinogenesis imperfecta and 26 hearing loss. The total number of fractures was 1,932. Thirty-four patients had intra-uterine fractures and nine had perinatal fractures. Surgery was performed 163 times in 58 patients; 100 osteosyntheses and 63 osteotomies. Bisphosphonate treatment was used in 37 patients. The number of affected individuals and predominance of severe forms of OI indicate that the disease is under diagnosed in Vietnam, especially in cases without a family history or with mild form of OI. Deformities appeared in all patients with different severity and localisation, affecting mostly the lower limbs. OI medical and surgical treatment rates are low and in most cases surgery was performed due to fractures. CONCLUSIONS: Compared to previous studies, our results indicate a lower OI prevalence and greater severity of symptoms in the Vietnamese population when compared with other areas. Further investigation, improved diagnosis and treatment are needed to increase the patients' quality of life.

Association analysis of class II cytokine and receptor genes in vitiligo patients.

The loss of melanocytes in vitiligo is mainly attributed to defective autoimmune mechanisms and lately autoinflammatory mediators have become more emphasized. Among these, a number of class II cytokines and their receptors have displayed altered expression patterns in vitiligo. Thus, we selected 30 SNPs from the regions of respective genes to be genotyped in Estonian case-control sample (109 and 328 individuals, respectively). For more precise analyses, patients were divided into subgroups based on vitiligo progression activity, age of onset, sex, occurrence of vitiligo among relatives, extent of depigmented areas, appearance of Köbner's phenomenon, existence of halo nevi, occurrence of spontaneous repigmentation, and amount of thyroid peroxidase antibodies. No associations appeared in whole vitiligo group. In subgroups, several allelic and haplotype associations were found. The strongest involved SNPs rs12301088 (near IL26 gene), that was associated with familial vitiligo and existence of halo nevi, and rs2257167 (IFNAR1 gene), that was associated with female vitiligo. Additionally, haplotypes consisting of rs12301088 and rs12321603 alleles (IL26-IL22 genes), that were associated with familial vitiligo and existence of halo nevi. In conclusion, several genetic associations with vitiligo subphenotypes were revealed and functional explanations to these remain to be determined in respective studies.

Smoking-induced expression of the GPR15 gene indicates its potential role in chronic inflammatory pathologies.

Despite the described clear epigenetic effects of smoking, the effect of smoking on genome-wide gene expression in the blood is obscure. We therefore studied the smoking-induced changes in the gene-expression profile of the peripheral blood. RNA was extracted from the whole blood of 48 individuals with a detailed smoking history (24 never-smokers, 16 smokers, and 8 ex-smokers). Gene-expression profiles were evaluated with RNA sequencing, and results were analyzed separately in 24 men and 24 women. In the male smokers, 13 genes were statistically significantly (false-discovery rate <0.1) differentially expressed; in female smokers, 5 genes. Although most of the differentially expressed genes were different between the male and female smokers, the G-protein-coupled receptor 15 gene (GPR15) was differentially expressed in both male and female smokers compared with never-smokers. Analysis of GPR15 methylation identified significantly greater hypomethylation in smokers compared with that in never-smokers. GPR15 is the chemoattractant receptor that regulates T-cell migration and immunity. Up-regulation of GPR15 could explain to some extent the health hazards of smoking with regard to chronic inflammatory diseases.

Melanocytes in the skin--comparative whole transcriptome analysis of main skin cell types.

Melanocytes possess several functions besides a role in pigment synthesis, but detailed characteristics of the cells are still unclear. We used whole transcriptome sequencing (RNA-Seq) to assess differential gene expression of cultivated normal human melanocytes with respect to keratinocytes, fibroblasts and whole skin. The present results reveal cultivated melanocytes as highly proliferative cells with possible stem cell-like properties. The enhanced readiness to regenerate makes melanocytes the most vulnerable cells in the skin and explains their high risk of developing into malignant melanoma.

Whole exome sequencing of a single osteosarcoma case--integrative analysis with whole transcriptome RNA-seq data.

BACKGROUND: Osteosarcoma (OS) is a prevalent primary malignant bone tumour with unknown etiology. These highly metastasizing tumours are among the most frequent causes of cancer-related deaths. Thus, there is an urgent need for different markers, and with our study, we were aiming towards finding novel biomarkers for OS. METHODS: For that, we analysed the whole exome of the tumorous and non-tumour bone tissue from the same patient with OS applying next-generation sequencing. For data analysis, we used several softwares and combined the exome data with RNA-seq data from our previous study. RESULTS: In the tumour exome, we found wide genomic rearrangements, which should qualify as chromotripsis-we detected almost 3,000 somatic single nucleotide variants (SNVs) and small indels and more than 2,000 copy number variants (CNVs) in different chromosomes. Furthermore, the somatic changes seem to be associated to bone tumours, whereas germline mutations to cancer in general. We confirmed the previous findings that the most significant pathway involved in OS pathogenesis is probably the WNT/ß-catenin signalling pathway. Also, the IGF1/IGF2 and IGF1R homodimer signalling and TP53 (including downstream tumour suppressor gene EI24) pathways may have a role. Additionally, the mucin family genes, especially MUC4 and cell cycle controlling gene CDC27 may be considered as potential biomarkers for OS. CONCLUSIONS: The genes, in which the mutations were detected, may be considered as targets for finding biomarkers for OS. As the study is based on a single case and only DNA and RNA analysis, further confirmative studies are required.

Expression of class II cytokine genes in children's skin.

Immune regulation of the skin plays an important role in susceptibility and development of illnesses. The aim of our study was to localise the interleukin (IL)-10 family of cytokines, in children's skin and to determine possible age-related differences in the expression level. The mRNA expression level of IL10, IL19, IL20, IL22, IL24, IL26, IL28B, IL29 and their receptors IL10RA, IL10RB, IL20RA, IL20RB, IL22RA1, IL22RA2, IL28RA was compared in skin biopsies of children and adults and in childrens' skin cells by quantitative real-time PCR (qRT-PCR). Immunohistochemistry was performed to confirm the qRT-PCR findings. We found age-related differences in the expression of IL10RB, IL20, IL20RA, IL22RA1, IL22RA2, IL26 and IL28RA genes. Cell type-dependent expression of IL10 family cytokines was apparent in the skin. In addition to previously known differences in systemic immunological response of adults and children, the present results reveal differences in immune profile of adult and juvenile skin.