Phylogenetic profiling and gene expression studies implicate a primary role of PSORS1C2 in terminal differentiation of keratinocytes

Intrapatient comparison of atopic dermatitis skin transcriptome shows differences between tape-strips and biopsies

BACKGROUND

Our knowledge of etiopathogenesis of atopic dermatitis (AD) is largely derived from skin biopsies, which are associated with pain, scarring and infection. In contrast, tape-stripping is a minimally invasive, nonscarring technique to collect skin samples.

METHODS

To construct a global AD skin transcriptomic profile comparing tape-strips to whole-skin biopsies, we performed RNA-seq on tape-strips and biopsies taken from the lesional skin of 20 moderate-to-severe AD patients and the skin of 20 controls. Differentially expressed genes (DEGs) were defined by fold-change (FCH) ≥2.0 and false discovery rate <0.05.

RESULTS

We detected 4104 (2513 Up; 1591 Down) and 1273 (546 Up; 727 Down) DEGs in AD versus controls, in tape-strips and biopsies, respectively. Although both techniques captured dysregulation of key immune genes, tape-strips showed higher FCHs for innate immunity (IL-1B, IL-8), dendritic cell (ITGAX/CD11C, FCER1A), Th2 (IL-13, CCL17, TNFRSF4/OX40), and Th17 (CCL20, CXCL1) products, while biopsies showed higher upregulation of Th22 associated genes (IL-22, S100As) and dermal cytokines (IFN-γ, CCL26). Itch-related genes (IL-31, TRPV3) were preferentially captured by tape-strips. Epidermal barrier abnormalities were detected in both techniques, with terminal differentiation defects (FLG2, PSORS1C2) better represented by tape-strips and epidermal hyperplasia changes (KRT16, MKI67) better detected by biopsies.

CONCLUSIONS

Tape-strips and biopsies capture overlapping but distinct features of the AD molecular signature, suggesting their respective utility for monitoring specific AD-related immune, itch, and barrier abnormalities in clinical trials and longitudinal studies.

Genome-wide allele and haplotype-sharing patterns suggested one unique Hmong-Mein-related lineage and biological adaptation history in Southwest China

BACKGROUND

Fine-scale genetic structure of ethnolinguistically diverse Chinese populations can fill the gap in the missing diversity and evolutionary landscape of East Asians, particularly for anthropologically informed Chinese minorities. Hmong-Mien (HM) people were one of the most significant indigenous populations in South China and Southeast Asia, which were suggested to be the descendants of the ancient Yangtze rice farmers based on linguistic and archeological evidence. However, their deep population history and biological adaptative features remained to be fully characterized.

OBJECTIVES

To explore the evolutionary and adaptive characteristics of the Miao people, we genotyped genome-wide SNP data in Guizhou HM-speaking populations and merged it with modern and ancient reference populations via a comprehensive population genetic analysis and evolutionary admixture modeling.

RESULTS

The overall genetic admixture landscape of Guizhou Miao showed genetic differentiation between them and other linguistically diverse Guizhou populations. Admixture models further confirmed that Miao people derived their primary ancestry from geographically close Guangxi Gaohuahua people. The estimated identity by descent and effective population size confirmed a plausible population bottleneck, contributing to their unique genetic diversity and population structure patterns. We finally identified several natural selection candidate genes associated with several biological pathways.

CONCLUSIONS

Guizhou Miao possessed a specific genetic structure and harbored a close genetic relationship with geographically close southern Chinese indigenous populations and Guangxi historical people. Miao people derived their major ancestry from geographically close Guangxi Gaohuahua people and experienced a plausible population bottleneck which contributed to the unique pattern of their genetic diversity and structure. Future ancient DNA from Shijiahe and Qujialing will provide new insights into the origin of the Miao people.

Comparative genomics reveals evolutionary loss of epiplakin in cetaceans

The adaptation of vertebrates to different environments was associated with changes in the molecular composition and regulation of epithelia. Whales and dolphins, together forming the clade cetaceans, have lost multiple epithelial keratins during or after their evolutionary transition from life on land to life in water. It is unknown whether the changes in keratins were accompanied by gain or loss of cytoskeletal adapter proteins of the plakin family. Here we investigated whether plakin proteins are conserved in cetaceans and other vertebrates. Comparative analysis of genome sequences showed conservation of dystonin, microtubule actin crosslinking factor 1 (MACF1), plectin, desmoplakin, periplakin and envoplakin in cetaceans. By contrast, EPPK1 (epiplakin) was disrupted by inactivating mutations in all cetaceans investigated. Orthologs of EPPK1 are present in bony and cartilaginous fishes and tetrapods, indicating an evolutionary origin of EPPK1 in a common ancestor of jawed vertebrates (Gnathostomes). In many vertebrates, EPPK1 is flanked by an as-yet uncharacterized gene that encodes protein domains homologous to the carboxy-terminal segment of MACF1. We conclude that epiplakin, unlike other plakins, was lost in cetaceans.

Gene duplications and gene loss in the epidermal differentiation complex during the evolutionary land-to-water transition of cetaceans

Major protein components of the mammalian skin barrier are encoded by genes clustered in the Epidermal Differentiation Complex (EDC). The skin of cetaceans, i.e. whales, porpoises and dolphins, differs histologically from that of terrestrial mammals. However, the genetic regulation of their epidermal barrier is only incompletely known. Here, we investigated the EDC of cetaceans by comparative genomics. We found that important epidermal cornification proteins, such as loricrin and involucrin are conserved and subtypes of small proline-rich proteins (SPRRs) are even expanded in numbers in cetaceans. By contrast, keratinocyte proline rich protein (KPRP), skin-specific protein 32 (XP32) and late-cornified envelope (LCE) genes with the notable exception of LCE7A have been lost in cetaceans. Genes encoding proline rich 9 (PRR9) and late cornified envelope like proline rich 1 (LELP1) have degenerated in subgroups of cetaceans. These data suggest that the evolution of an aquatic lifestyle was accompanied by amplification of SPRR genes and loss of specific other epidermal differentiation genes in the phylogenetic lineage leading to cetaceans.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

The molecular features of normal and atopic dermatitis skin in infants, children, adolescents, and adults

BACKGROUND

Although atopic dermatitis (AD) often presents in infancy and persists into adulthood, comparative characterization of AD skin among different pediatric age groups is lacking.

OBJECTIVE

We sought to define skin biopsy profiles of lesional and nonlesional AD across different age groups (0-5-year-old infants with disease duration <6 months, 6-11-year-old children, 12-17-year-old adolescents, ≥18-year-old adults) versus age-appropriate controls.

METHODS

We performed gene expression analyses by RNA-sequencing and real-time PCR (RT-PCR) and protein expression analysis using immunohistochemistry.

RESULTS

T_H2/T_H22 skewing, including IL-13, CCL17/thymus and activation-regulated chemokine, IL-22, and S100As, characterized the common AD signature, with a global pathway-level enrichment across all ages. Nevertheless, specific cytokines varied widely. For example, IL-33, IL-1RL1/IL-33R, and IL-9, often associated with early atopic sensitization, showed greatest upregulations in infants. T_H17 inflammation presented a 2-peak curve, with highest increases in infants (including IL-17A and IL-17F), followed by adults. T_H1 polarization was uniquely detected in adults, even when compared with adolescents, with significant upregulation in adults of IFN-γ and CXCL9/CXCL10/CXCL11. Although all AD age groups had barrier abnormalities, only adults had significant decreases in filaggrin expression. Despite the short duration of the disease, infant AD presented robust downregulations of multiple barrier-related genes in both lesional and nonlesional skin. Clinical severity scores significantly correlated with T_H2/T_H22-related markers in all pediatric age groups.

CONCLUSIONS

The shared signature of AD across ages is T_H2/T_H22-skewed, yet differential expression of specific T_H2/T_H22-related genes, other T_H pathways, and barrier-related genes portray heterogenetic, age-specific molecular fingerprints.

Heritability Enrichment of Immunoglobulin G N-Glycosylation in Specific Tissues

Genome-wide association studies (GWAS) have identified over 60 genetic loci associated with immunoglobulin G (IgG) N-glycosylation; however, the causal genes and their abundance in relevant tissues are uncertain. Leveraging data from GWAS summary statistics for 8,090 Europeans, and large-scale expression quantitative trait loci (eQTL) data from the genotype-tissue expression of 53 types of tissues (GTEx v7), we derived a linkage disequilibrium score for the specific expression of genes (LDSC-SEG) and conducted a transcriptome-wide association study (TWAS). We identified 55 gene associations whose predicted levels of expression were significantly associated with IgG N-glycosylation in 14 tissues. Three working scenarios, i.e., tissue-specific, pleiotropic, and coassociated, were observed for candidate genetic predisposition affecting IgG N-glycosylation traits. Furthermore, pathway enrichment showed several IgG N-glycosylation-related pathways, such as asparagine N-linked glycosylation, N-glycan biosynthesis and transport to the Golgi and subsequent modification. Through phenome-wide association studies (PheWAS), most genetic variants underlying TWAS hits were found to be correlated with health measures (height, waist-hip ratio, systolic blood pressure) and diseases, such as systemic lupus erythematosus, inflammatory bowel disease, and Parkinson's disease, which are related to IgG N-glycosylation. Our study provides an atlas of genetic regulatory loci and their target genes within functionally relevant tissues, for further studies on the mechanisms of IgG N-glycosylation and its related diseases.

Differential Evolution of the Epidermal Keratin Cytoskeleton in Terrestrial and Aquatic Mammals

Keratins are the main intermediate filament proteins of epithelial cells. In keratinocytes of the mammalian epidermis they form a cytoskeleton that resists mechanical stress and thereby are essential for the function of the skin as a barrier against the environment. Here, we performed a comparative genomics study of epidermal keratin genes in terrestrial and fully aquatic mammals to determine adaptations of the epidermal keratin cytoskeleton to different environments. We show that keratins K5 and K14 of the innermost (basal), proliferation-competent layer of the epidermis are conserved in all mammals investigated. In contrast, K1 and K10, which form the main part of the cytoskeleton in the outer (suprabasal) layers of the epidermis of terrestrial mammals, have been lost in whales and dolphins (cetaceans) and in the manatee. Whereas in terrestrial mammalian epidermis K6 and K17 are expressed only upon stress-induced epidermal thickening, high levels of K6 and K17 are consistently present in dolphin skin, indicating constitutive expression and substitution of K1 and K10. K2 and K9, which are expressed in a body site-restricted manner in human and mouse suprabasal epidermis, have been lost not only in cetaceans and manatee but also in some terrestrial mammals. The evolution of alternative splicing of K10 and differentiation-dependent upregulation of K23 have increased the complexity of keratin expression in the epidermis of terrestrial mammals. Taken together, these results reveal evolutionary diversification of the epidermal cytoskeleton in mammals and suggest a complete replacement of the quantitatively predominant epidermal proteins of terrestrial mammals by originally stress-inducible keratins in cetaceans.

Early-onset pediatric atopic dermatitis is characterized by TH2/TH17/TH22-centered inflammation and lipid alterations

BACKGROUND

Although atopic dermatitis (AD) often starts in early childhood, detailed tissue profiling of early-onset AD in children is lacking, hindering therapeutic development for this patient population with a particularly high unmet need for better treatments.

OBJECTIVE

We sought to globally profile the skin of infants with AD compared with that of adults with AD and healthy control subjects.

METHODS

We performed microarray, RT-PCR, and fluorescence microscopy studies in infants and young children (<5 years old) with early-onset AD (<6 months disease duration) compared with age-matched control subjects and adults with longstanding AD.

RESULTS

Transcriptomic analyses revealed profound differences between pediatric patients with early-onset versus adult patients with longstanding AD in not only lesional but also nonlesional tissues. Although both patient populations harbored TH2-centered inflammation, pediatric AD also showed significant TH17/TH22 skewing but lacked the TH1 upregulation that characterizes adult AD. Pediatric AD exhibited relatively normal expression of epidermal differentiation and cornification products, which is downregulated in adults with AD. Defects in the lipid barrier (eg, ELOVL fatty acid elongase 3 [ELOVL3] and diacylglycerol o-acyltransferase 2 [DGAT2]) and tight junction regulation (eg, claudins 8 and 23) were evident in both groups. However, some lipid-associated mediators (eg, fatty acyl-CoA reductase 2 and fatty acid 2-hydroxylase) showed preferential downregulation in pediatric AD, and lipid barrier genes (FA2H and DGAT2) showed inverse correlations with transepidermal water loss, a functional measure of the epidermal barrier.

CONCLUSIONS

Skin samples from children and adult patients with AD share lipid metabolism and tight junction alterations, but epidermal differentiation complex defects are only present in adult AD, potentially resulting from chronic immune aberration that is not yet present in early-onset disease.

Comparative Analysis of Epidermal Differentiation Genes of Crocodilians Suggests New Models for the Evolutionary Origin of Avian Feather Proteins

The epidermis of amniotes forms a protective barrier against the environment and the differentiation program of keratinocytes, the main cell type in the epidermis, has undergone specific alterations in the course of adaptation of amniotes to a broad variety of environments and lifestyles. The epidermal differentiation complex (EDC) is a cluster of genes expressed at late stages of keratinocyte differentiation in both sauropsids and mammals. In the present study, we identified and analyzed the crocodilian equivalent of the EDC. The gene complement of the EDC of both the American alligator and the saltwater crocodile were determined by comparative genomics, de novo gene prediction and identification of EDC transcripts in published transcriptome data. We found that crocodilians have an organization of the EDC similar to that of their closest living relatives, the birds, with which they form the clade Archosauria. Notable differences include the specific expansion of a subfamily of EDC genes in crocodilians and the loss of distinct ancestral EDC genes in birds. Identification and comparative analysis of crocodilian orthologs of avian feather proteins suggest that the latter evolved by cooption and sequence modification of ancestral EDC genes, and that the amplification of an internal highly cysteine-enriched amino acid sequence motif gave rise to the feather component epidermal differentiation cysteine-rich protein in the avian lineage. Thus, sequence diversification of EDC genes contributed to the evolutionary divergence of the crocodilian and avian integuments.

Evolution of the MC5R gene in placental mammals with evidence for its inactivation in multiple lineages that lack sebaceous glands

MC5R is one of five melanocortin receptor genes found in placental mammals. MC5R plays an important role in energy homeostasis and is also expressed in the terminal differentiation of sebaceous glands. Among placental mammals there are multiple lineages that either lack or have degenerative sebaceous glands including Cetacea (whales, dolphins, and porpoises), Hippopotamidae (hippopotamuses), Sirenia (manatees and dugongs), Proboscidea (elephants), Rhinocerotidae (rhinos), and Heterocephalus glaber (naked mole rat). Given the loss or diminution of sebaceous glands in these taxa, we procured MC5R sequences from publicly available genomes and transcriptomes, supplemented by a newly generated sequence for Choeropsis liberiensis (pygmy hippopotamus), to determine if this gene remains intact or is inactivated in association with loss/reduction of sebaceous glands. Our data set includes complete MC5R sequences for 114 placental mammal species including two individuals of Mammuthus primigenius (woolly mammoth) from Oimyakon and Wrangel Island. Complete loss or inactivation of the MC5R gene occurs in multiple placental lineages that have lost sebaceous glands (Cetacea, West Indian manatee, African elephant, white rhinoceros) or are characterized by unusual skin (pangolins, aardvarks). Both M. primigenius individuals share inactivating mutations with the African elephant even though sebaceous glands have been reported in the former. MC5R remains intact in hippopotamuses and the naked mole rat, although slightly elevated dN/dS ratios in these lineages allow for the possibility that the accumulation of inactivating mutations in MC5R may lag behind the relaxation of purifying selection. For Cetacea and Hippopotamidae, the absence of shared inactivating mutations in two different skin genes (MC5R, PSORS1C2) is consistent with the hypothesis that semi-aquatic lifestyles were acquired independently in these clades following divergence from a common ancestor.

Epidermal cornification is preceded by the expression of a keratinocyte-specific set of pyroptosis-related genes

The homeostasis of the epidermis depends on keratinocyte differentiation and cornification, a mode of programmed cell death that does not elicit inflammation. Here, we report that cornification is associated with the expression of specific genes that control multiple steps of pyroptosis, another form of cell death that involves the processing and release of interleukin-1 family (IL1F) cytokines. Expression levels of pro-inflammatory IL1A and IL1B and of the pyroptotic pore-forming gasdermin (GSDM) D were downregulated during terminal differentiation of human keratinocytes in vitro. By contrast, negative regulators of IL-1 processing, including NLR family pyrin domain containing 10 (NLRP10) and pyrin domain-containing 1 (PYDC1), the anti-inflammatory IL1F members IL-37 (IL1F7) and IL-38 (IL1F10), and GSDMA, were strongly induced in differentiated keratinocytes. In human tissues, these keratinocyte differentiation-associated genes are expressed in the skin at higher levels than in any other organ, and mammalian species, that have lost the epidermal cornification program during evolution, i.e. whales and dolphins, lack homologs of these genes. Together, our results suggest that human epidermal cornification is accompanied by a tight control of pyroptosis and warrant further studies of potential defects in the balance between cornification and pyroptosis in skin pathologies.

Filaggrin has evolved from an "S100 fused-type protein" (SFTP) gene present in a common ancestor of amphibians and mammals

The expression of filaggrin in differentiated keratinocytes and the association of filaggrin mutations with ichthyosis vulgaris and atopic dermatitis suggest that this prototypical member of the S100 fused-type protein (SFTP) family plays a key role in the epidermal barrier to the environment. Here, we report that SFTP genes are present not only in amniotes but also in amphibians. Four SFTPs are expressed in the skin of the frog Xenopus laevis. The results of this study indicate that filaggrin has evolved from an ancestral SFTP that may have contributed to skin modifications during the evolutionary transition to terrestrial life.