Inactivation of Autophagy in Keratinocytes Reduces Tumor Growth in Mouse Models of Epithelial Skin Cancer.

Autophagy is a ubiquitous degradation mechanism, which plays a critical role in cellular homeostasis. To test whether autophagy suppresses or supports the growth of tumors in the epidermis of the skin, we inactivated the essential autophagy gene Atg7 specifically in the epidermal keratinocytes of mice (Atg7∆ep) and subjected such mutant mice and fully autophagy-competent mice to tumorigenesis. The lack of epithelial Atg7 did not prevent tumor formation in response to 7, 12-dimethylbenz(a)anthracene (DMBA) as the initiator and 12-O tetradecanoylphorbol-13-acetate (TPA) as the promoter of tumor growth. However, the number of tumors per mouse was reduced in mice with epithelial Atg7 deficiency. In the K5-SOS EGFRwa2/wa2 mouse model, epithelial tumors were initiated by Son of sevenless (SOS) in response to wounding. Within 12 weeks after tumor initiation, 60% of the autophagy-competent K5-SOS EGFRwa2/wa2 mice had tumors of 1 cm diameter and had to be sacrificed, whereas none of the Atg7∆ep K5-SOS EGFRwa2/wa2 mice formed tumors of this size. In summary, the deletion of Atg7 reduced the growth of epithelial tumors in these two mouse models of skin cancer. Thus, our data show that the inhibition of autophagy limits the growth of epithelial skin tumors.

Identification of New Biological Pathways Involved in Skin Aging From the Analysis of French Women Genome-Wide Data.

Skin aging is an ineluctable process leading to the progressive loss of tissue integrity and is characterized by various outcomes such as wrinkling and sagging. Researchers have identified impacting environmental factors (sun exposure, smoking, etc.) and several molecular mechanisms leading to skin aging. We have previously performed genome-wide association studies (GWAS) in 502 very-well characterized French women, looking for associations with four major outcomes of skin aging, namely, photoaging, solar lentigines, wrinkling, and sagging, and this has led to new insights into the molecular mechanisms of skin aging. Since individual SNP associations in GWAS explain only a small fraction of the genetic impact in complex polygenic phenotypes, we have made the integration of these genotypes into the reference Kegg biological pathways and looked for associations by the gene set enrichment analysis (GSEA) approach. 106 pathways were tested for association with the four outcomes of skin aging. This biological pathway analysis revealed new relevant pathways and genes, some likely specific of skin aging such as the WNT7B and PRKCA genes in the "melanogenesis" pathway and some likely involved in global aging such as the DDB1 gene in the "nucleotide excision repair" pathway, not picked up in the previously published GWAS. Overall, our results suggest that the four outcomes of skin aging possess specific molecular mechanisms such as the "proteasome" and "mTOR signaling pathway" but may also share common molecular mechanisms such as "nucleotide excision repair."

Single-cell transcriptomics defines keratinocyte differentiation in avian scutate scales.

The growth of skin appendages, such as hair, feathers and scales, depends on terminal differentiation of epidermal keratinocytes. Here, we investigated keratinocyte differentiation in avian scutate scales. Cells were isolated from the skin on the legs of 1-day old chicks and subjected to single-cell transcriptomics. We identified two distinct populations of differentiated keratinocytes. The first population was characterized by mRNAs encoding cysteine-rich keratins and corneous beta-proteins (CBPs), also known as beta-keratins, of the scale type, indicating that these cells form hard scales. The second population of differentiated keratinocytes contained mRNAs encoding cysteine-poor keratins and keratinocyte-type CBPs, suggesting that these cells form the soft interscale epidermis. We raised an antibody against keratin 9-like cysteine-rich 2 (KRT9LC2), which is encoded by an mRNA enriched in the first keratinocyte population. Immunostaining confirmed expression of KRT9LC2 in the suprabasal epidermal layers of scutate scales but not in interscale epidermis. Keratinocyte differentiation in chicken leg skin resembled that in human skin with regard to the transcriptional upregulation of epidermal differentiation complex genes and genes involved in lipid metabolism and transport. In conclusion, this study defines gene expression programs that build scutate scales and interscale epidermis of birds and reveals evolutionarily conserved keratinocyte differentiation genes.

The serine proteases dipeptidyl-peptidase 4 and urokinase are key molecules in human and mouse scar formation.

Despite recent advances in understanding skin scarring, mechanisms triggering hypertrophic scar formation are still poorly understood. In the present study, we investigate mature human hypertrophic scars and developing scars in mice at single cell resolution. Compared to normal skin, we find significant differences in gene expression in most cell types present in scar tissue. Fibroblasts show the most prominent alterations in gene expression, displaying a distinct fibrotic signature. By comparing genes upregulated in murine fibroblasts during scar development with genes highly expressed in mature human hypertrophic scars, we identify a group of serine proteases, tentatively involved in scar formation. Two of them, dipeptidyl-peptidase 4 (DPP4) and urokinase (PLAU), are further analyzed in functional assays, revealing a role in TGFß1-mediated myofibroblast differentiation and over-production of components of the extracellular matrix in vitro. Topical treatment with inhibitors of DPP4 and PLAU during scar formation in vivo shows anti-fibrotic activity and improvement of scar quality, most prominently after application of the PLAU inhibitor BC-11. In this study, we delineate the genetic landscape of hypertrophic scars and present insights into mechanisms involved in hypertrophic scar formation. Our data suggest the use of serine protease inhibitors for the treatment of skin fibrosis.

Epilipidomics of Senescent Dermal Fibroblasts Identify Lysophosphatidylcholines as Pleiotropic Senescence-Associated Secretory Phenotype (SASP) Factors.

During aging, skin accumulates senescent cells. The transient presence of senescent cells, followed by their clearance by the immune system, is important in tissue repair and homeostasis. The persistence of senescent cells that evade clearance contributes to the age-related deterioration of the skin. The senescence-associated secretory phenotype of these cells contains immunomodulatory molecules that facilitate clearance but also promote chronic damage. Here, we investigated the epilipidome-the oxidative modifications of phospholipids-of senescent dermal fibroblasts, because these molecules are among the bioactive lipids that were recently identified as senescence-associated secretory phenotype factors. Using replicative- and stress- induced senescence protocols, we identified lysophosphatidylcholines as universally elevated in senescent fibroblasts, whereas other oxidized lipids displayed a pattern that was characteristic for the used senescence protocol. When we tested the lysophosphatidylcholines for senescence-associated secretory phenotype activity, we found that they elicit chemokine release in nonsenescent fibroblasts but also interfere with toll-like receptor 2 and 6/CD36 signaling and phagocytic capacity in macrophages. Using matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance mass spectrometry imaging, we localized two lysophosphatidylcholine species in aged skin. This suggests that lysophospholipids may facilitate immune evasion and low-grade chronic inflammation in skin aging.

Cerebellar Degeneration-related Antigen 1 Is Ubiquitously Expressed in Human Epidermis and Dermis.

Cerebellar degeneration-related antigen 1 (CDR1) was described to be expressed in the nervous system and in different types of cancer tissues. In the present study, we demonstrate that CDR1 is in addition ubiquitously expressed in human epidermis, dermis and isolated skin cells. Both CDR1 mRNA and protein were detected in human skin-derived mast cells, melanocytes, fibroblasts and keratinocytes, suggesting that CDR1 does not have a neuron-specific function.

Pangolins Lack IFIH1/MDA5, a Cytoplasmic RNA Sensor That Initiates Innate Immune Defense Upon Coronavirus Infection.

Zoonotic infections are an imminent threat to human health. Pangolins were recently identified as carriers and intermediate hosts of coronaviruses. Previous research has shown that infection with coronaviruses activates an innate immune response upon sensing of viral RNA by interferon-induced with helicase C domain 1 (IFIH1), also known as MDA5. Here, we performed a comparative genomics study of RNA sensor genes in three species of pangolins. DDX58/RIG-I, a sensor of cytoplasmic viral RNA and toll-like receptors (TLR) 3, 7, and 8, which bind RNA in endosomes, are conserved in pangolins. By contrast, IFIH1 a sensor of intracellular double-stranded RNA, has been inactivated by mutations in pangolins. Likewise, Z-DNA-binding protein (ZBP1), which senses both Z-DNA and Z-RNA, has been lost during the evolution of pangolins. These results suggest that the innate immune response to viruses differs significantly between pangolins and other mammals, including humans. We put forward the hypothesis that loss of IFIH1 and ZBP1 provided an evolutionary advantage by reducing inflammation-induced damage to host tissues and thereby contributed to a switch from resistance to tolerance of viral infections in pangolins.

Cytosolic DNA sensing through cGAS and STING is inactivated by gene mutations in pangolins.

The release of DNA into the cytoplasm upon damage to the nucleus or during viral infection triggers an interferon-mediated defense response, inflammation and cell death. In human cells cytoplasmic DNA is sensed by cyclic GMP-AMP Synthase (cGAS) and Absent In Melanoma 2 (AIM2). Here, we report the identification of a "natural knockout" model of cGAS. Comparative genomics of phylogenetically diverse mammalian species showed that cGAS and its interaction partner Stimulator of Interferon Genes (STING) have been inactivated by mutations in the Malayan pangolin whereas other mammals retained intact copies of these genes. The coding sequences of CGAS and STING1 are also disrupted by premature stop codons and frame-shift mutations in Chinese and tree pangolins, suggesting that expression of these genes was lost in a common ancestor of all pangolins that lived more than 20 million years ago. AIM2 is retained in a functional form in pangolins whereas it is inactivated by mutations in carnivorans, the phylogenetic sister group of pangolins. The deficiency of cGAS and STING points to the existence of alternative mechanisms of controlling cytoplasmic DNA-associated cell damage and viral infections in pangolins.

TINCR is not a non-coding RNA but encodes a protein component of cornified epidermal keratinocytes.

Long non-coding RNAs have been implicated in the regulation of a plethora of biological processes, yet it has been challenging to verify that they are truly not coding for proteins. Terminal differentiation-induced non-coding RNA (TINCR) is a 3.7-kilobase mRNA that is highly abundant in epidermal keratinocytes prior to cornification. Here, we report the presence of an evolutionarily conserved open reading frame in TINCR and the identification of peptides derived from this open reading frame in the proteome of human stratum corneum. Our results demonstrate that TINCR is a protein-coding RNA and suggest that the TINCR-encoded protein is involved in keratinocyte cornification.

Deciphering the functional heterogeneity of skin fibroblasts using single-cell RNA sequencing.

Though skin fibroblasts (FB) are the main cell population within the dermis, the different skin FB subsets are not well characterized and the traditional classification into reticular and papillary FBs has little functional relevance. To fill the gap of knowledge on FB diversity in human skin, we performed single-cell RNA sequencing. Investigation of marker genes for the different skin cell subtypes revealed a heterogeneous picture of FBs. When mapping reticular and papillary FB markers, we could not detect cluster specificity, suggesting that these two populations show a higher transcriptional heterogeneity than expected. This finding was further confirmed by in situ hybridization, showing that DPP4 was expressed in both dermal layers. Our analysis identified six FB clusters with distinct transcriptional signatures. Importantly, we could demonstrate that in human skin DPP4+ FBs are the main producers of factors involved in extracellular matrix (ECM) assembly. In conclusion, we provide evidence that hitherto considered FB markers are not ideal to characterize skin FB subpopulations in single-cell sequencing analyses. The identification of DPP4+ FBs as the main ECM-producing cells in human skin will foster the development of anti-fibrotic treatments for the skin and other organs.

ATG7 is essential for secretion of iron from ameloblasts and normal growth of murine incisors during aging.

The incisors of rodents comprise an iron-rich enamel and grow throughout adult life, making them unique models of iron metabolism and tissue homeostasis during aging. Here, we deleted Atg7 (autophagy related 7) in murine ameloblasts, i.e. the epithelial cells that produce enamel. The absence of ATG7 blocked the transport of iron from ameloblasts into the maturing enamel, leading to a white instead of yellow surface of maxillary incisors. In aging mice, lack of ATG7 was associated with the growth of ectopic incisors inside severely deformed primordial incisors. These results suggest that 2 characteristic features of rodent incisors, i.e. deposition of iron on the enamel surface and stable growth during aging, depend on autophagic activity in ameloblasts. Abbreviations: ATG5: autophagy related 5; ATG7: autophagy related 7; CMV: cytomegalovirus; Cre: Cre recombinase; CT: computed tomography; FTH1: ferritin heavy polypeptide 1; GFP: green fluorescent protein; KRT5: keratin 5; KRT14: keratin 14; LGALS3: lectin, galactose binding, soluble 3; MAP1LC3/LC3: microtubule-associated protein 1 light chain 3; NCOA4: nuclear receptor coactivator 4; NRF2: nuclear factor, erythroid 2 like 2; SQSTM1: sequestosome 1.

Convergent Evolution of Cysteine-Rich Keratins in Hard Skin Appendages of Terrestrial Vertebrates.

Terrestrial vertebrates have evolved hard skin appendages, such as scales, claws, feathers, and hair that play crucial roles in defense, predation, locomotion, and thermal insulation. The mechanical properties of these skin appendages are largely determined by cornified epithelial components. So-called "hair keratins," cysteine-rich intermediate filament proteins that undergo covalent cross-linking via disulfide bonds, are the crucial structural proteins of hair and claws in mammals and hair keratin orthologs are also present in lizard claws, indicating an evolutionary origin in a hairless common ancestor of amniotes. Here, we show that reptiles and birds have also other cysteine-rich keratins which lack cysteine-rich orthologs in mammals. In addition to hard acidic (type I) sauropsid-specific (HAS) keratins, we identified hard basic (type II) sauropsid-specific (HBS) keratins which are conserved in lepidosaurs, turtles, crocodilians, and birds. Immunohistochemical analysis with a newly made antibody revealed expression of chicken HBS1 keratin in the cornifying epithelial cells of feathers. Molecular phylogenetics suggested that the high cysteine contents of HAS and HBS keratins evolved independently from the cysteine-rich sequences of hair keratin orthologs, thus representing products of convergent evolution. In conclusion, we propose an evolutionary model in which HAS and HBS keratins evolved as structural proteins in epithelial cornification of reptiles and at least one HBS keratin was co-opted as a component of feathers after the evolutionary divergence of birds from reptiles. Thus, cytoskeletal proteins of hair and feathers are products of convergent evolution and evolutionary co-option to similar biomechanical functions in clade-specific hard skin appendages.

Tissue-regenerative potential of the secretome of γ-irradiated peripheral blood mononuclear cells is mediated via TNFRSF1B-induced necroptosis.

Peripheral blood mononuclear cells (PBMCs) have been shown to produce and release a plethora of pro-angiogenetic factors in response to γ-irradiation, partially accounting for their tissue-regenerative capacity. Here, we investigated whether a certain cell subtype of PBMCs is responsible for this effect, and whether the type of cell death affects the pro-angiogenic potential of bioactive molecules released by γ-irradiated PBMCs. PBMCs and PBMC subpopulations, including CD4+ and CD8+ T cells, B cells, monocytes, and natural killer cells, were isolated and subjected to high-dose γ-irradiation. Transcriptome analysis revealed subpopulation-specific responses to γ-irradiation with distinct activation of pro-angiogenic pathways, cytokine production, and death receptor signalling. Analysis of the proteins released showed that interactions of the subsets are important for the generation of a pro-angiogenic secretome. This result was confirmed at the functional level by the finding that the secretome of γ-irradiated PBMCs displayed higher pro-angiogenic activity in an aortic ring assay. Scanning electron microscopy and image stream analysis of γ-irradiated PBMCs revealed distinct morphological changes, indicative for apoptotic and necroptotic cell death. While inhibition of apoptosis had no effect on the pro-angiogenic activity of the secretome, inhibiting necroptosis in stressed PBMCs abolished blood vessel sprouting. Mechanistically, we identified tumor necrosis factor (TNF) receptor superfamily member 1B as the main driver of necroptosis in response to γ-irradiation in PBMCs, which was most likely mediated via membrane-bound TNF-α. In conclusion, our study demonstrates that the pro-angiogenic activity of the secretome of γ-irradiated PBMCs requires interplay of different PBMC subpopulations. Furthermore, we show that TNF-dependent necroptosis is an indispensable molecular process for conferring tissue-regenerative activity and for the pro-angiogenic potential of the PBMC secretome. These findings contribute to a better understanding of secretome-based therapies in regenerative medicine.

Immunolocalization and phylogenetic profiling of the feather protein with the highest cysteine content.

Feathers are the most complex skin appendages of vertebrates. Mature feathers consist of interconnected dead keratinocytes that are filled with heavily cross-linked proteins. Although the molecular architecture determines essential functions of feathers, only few feather proteins have been characterized with regard to their amino acid sequences and evolution. Here, we identify Epidermal Differentiation protein containing DPCC Motifs (EDDM) as a cysteine-rich protein that has co-evolved with other feather proteins. The EDDM gene is located within the avian epidermal differentiation complex (EDC), a cluster of genes that has originated and diversified in amniotes. EDDM shares the exon-intron organization with EDC genes of other amniotes, including humans, and a gene encoding an EDDM-like protein is present in crocodilians, suggesting that avian EDDM arose by sequence modification of an epidermal differentiation gene present in a common ancestor of archosaurs. The EDDM protein contains multiple sequence repeats and a higher number of cysteine residues than any other protein encoded in the EDC. Immunohistochemical analysis of chicken skin and skin appendages showed expression of EDDM in barb and barbules of feathers as well as in the subperiderm on embryonic scutate scales. These results suggest that the diversification and differential expression of EDDM, besides other EDC genes, was instrumental in facilitating the evolution of the most complex molecular architecture of feathers.

A genome wide association study identifies new genes potentially associated with eyelid sagging.

Sagging eyelid is considered as an outward of skin ageing and may cause medical issues. However, little is known about the factors involved in sagging eyelid. The study, which aims at determining genetic risk factors for eyelid sagging, was conducted in a cohort of 502 unrelated Caucasian women living in the Paris region. All included participants were aged between 44 and 70 years old (mean age, 57.6 years old). The severity of sagging eyelid was graded in 6 categories by a dermatologist using standardized photographs of the face. A genome wide association study adjusted on potential risk factors (including age and smoking habits) was conducted to identify genetic associations. Two single nucleotide polymorphisms in total linkage disequilibrium on chromosome 10, rs16927253 (P = 7.07 × 10-10 ) and rs4746957 (P = 1.06 × 10-8 ), were significantly associated with eyelid sagging severity. The rs16927253-T and rs4746957-A alleles showed a dominant protective effect towards eyelid sagging. These polymorphisms are located in intronic parts of the H2AFY2 gene which encodes a member of the H2A histone family and very close to the AIFM2 gene that induces apoptosis. Additionally, single nucleotide polymorphisms with a false discovery rate below 0.25 were located nearby the type XIII collagen COL13A1 gene on chromosome 10 and in the ADAMTS18 gene on chromosome 16. Several relevant genes were identified by the genome wide association study for their potential role in the sagging eyelid severity.

A novel role for NUPR1 in the keratinocyte stress response to UV oxidized phospholipids.

Ultraviolet light is the dominant environmental oxidative skin stressor and a major skin aging factor. We studied which oxidized phospholipid (OxPL) mediators would be generated in primary human keratinocytes (KC) upon exposure to ultraviolet A light (UVA) and investigated the contribution of OxPL to UVA responses. Mass spectrometric analysis immediately or 24 h post UV stress revealed significant changes in abundance of 173 and 84 lipid species, respectively. We identified known and novel lipid species including known bioactive and also potentially reactive carbonyl containing species. We found indication for selective metabolism and degradation of selected reactive lipids. Exposure to both UVA and to in vitro UVA - oxidized phospholipids activated, on transcriptome and proteome level, NRF2/antioxidant response signaling, lipid metabolizing enzyme expression and unfolded protein response (UPR) signaling. We identified NUPR1 as an upstream regulator of UVA/OxPL transcriptional stress responses and found this protein to be expressed in the epidermis. Silencing of NUPR1 resulted in augmented expression of antioxidant and lipid detoxification genes and disturbed the cell cycle, making it a potential key factor in skin reactive oxygen species (ROS) responses intimately involved in aging and pathology.

Establishment of keratinocyte cell lines from human hair follicles.

The advent of organotypic skin models advanced the understanding of complex mechanisms of keratinocyte differentiation. However, these models are limited by both availability of primary keratinocytes and donor variability. Keratinocytes derived from cultured hair follicles and interfollicular epidermis were immortalized by ectopic expression of SV40 and hTERT. The generated keratinocyte cell lines differentiated into stratified epidermis with well-defined stratum granulosum and stratum corneum in organotypic human skin models. They behaved comparable to primary keratinocytes regarding the expression of differentiation-associated proteins, cell junction components and proteins associated with cornification and formed a barrier against biotin diffusion. Mechanistically, we found that SV40 large T-antigen expression, accompanied by a strong p53 accumulation, was only detectable in the basal layer of the in vitro reconstructed epidermis. Inhibition of DNA-methylation resulted in expression of SV40 large T-antigen also in the suprabasal epidermal layers and led to incomplete differentiation of keratinocyte cell lines. Our study demonstrates the generation of keratinocyte cell lines which are able to fully differentiate in an organotypic skin model. Since hair follicles, as source for keratinocytes, can be obtained by minimally invasive procedures, our approach enables the generation of cell lines also from individuals not available for skin biopsies.

Control of cell death-associated danger signals during cornification prevents autoinflammation of the skin.

The function of the skin as a barrier to the environment is mainly achieved by the outermost layers of the epidermis. In the granular layer, epidermal keratinocytes undergo the last steps of their terminal differentiation program resulting in cornification. The coordinated conversion of living keratinocytes into corneocytes, the building blocks of the cornified layer, represents a unique form of programmed cell death. Recent studies have identified numerous genes that are specifically expressed in terminally differentiated keratinocytes and, surprisingly, this genetic program does not only include mediators of cornification but also suppressors of pyroptosis, another mode of programmed cell death. Pyroptosis is activated by inflammasomes, leads to the release of interleukin-1 (IL-1) family cytokines, and thereby activates inflammation. In addition, inhibitors of potentially pro-inflammatory proteases and enzymes removing danger-associated cytoplasmic DNA are expressed in differentiated keratinocytes. We propose the concept of cornification as an inherently hazardous process in which damaging side effects are actively suppressed by protective mechanisms. In support of this hypothesis, loss-of-function mutations in epidermal protease inhibitors and IL-1 family antagonists suffice to induce autoinflammation. Similarly, exogenous disturbances of either cornification or its accompanying control mechanisms may be starting points for skin inflammation. Further studies into the relationship between cornification, pyroptosis and other forms of cell death will help to define the initiation phase of inflammatory skin diseases and offer new targets for disease prevention and therapy.

Blocking negative effects of senescence in human skin fibroblasts with a plant extract.

There is increasing evidence that senescent cells are a driving force behind many age-related pathologies and that their selective elimination increases the life- and healthspan of mice. Senescent cells negatively affect their surrounding tissue by losing their cell specific functionality and by secreting a pro-tumorigenic and pro-inflammatory mixture of growth hormones, chemokines, cytokines and proteases, termed the senescence-associated secretory phenotype (SASP). Here we identified an extract from the plant Solidago virgaurea subsp. alpestris, which exhibited weak senolytic activity, delayed the acquisition of a senescent phenotype and induced a papillary phenotype with improved functionality in human dermal fibroblasts. When administered to stress-induced premature senescent fibroblasts, this extract changed their global mRNA expression profile and particularly reduced the expression of various SASP components, thereby ameliorating the negative influence on nearby cells. Thus, the investigated plant extract represents a promising possibility to block age-related loss of tissue functionality.

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.