Recent Positive Selection in Genes of the Mammalian Epidermal Differentiation Complex Locus

A spatial human thymus cell atlas mapped to a continuous tissue axis

T cells develop from circulating precursors, which enter the thymus and migrate throughout specialised sub-compartments to support maturation and selection. This process starts already in early fetal development and is highly active until the involution of the thymus in adolescence. To map the micro-anatomical underpinnings of this process in pre- vs. post-natal states, we undertook a spatially resolved analysis and established a new quantitative morphological framework for the thymus, the Cortico-Medullary Axis. Using this axis in conjunction with the curation of a multimodal single-cell, spatial transcriptomics and high-resolution multiplex imaging atlas, we show that canonical thymocyte trajectories and thymic epithelial cells are highly organised and fully established by post-conception week 12, pinpoint TEC progenitor states, find that TEC subsets and peripheral tissue genes are associated with Hassall's Corpuscles and uncover divergence in the pace and drivers of medullary entry between CD4 vs. CD8 T cell lineages. These findings are complemented with a holistic toolkit for spatial analysis and annotation, providing a basis for a detailed understanding of T lymphocyte development.

Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Tensor decomposition based feature extraction and classification to detect natural selection from genomic data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under non-convex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data while preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx , which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Rapid divergence of a gamete recognition gene promoted macroevolution of Eutheria

BACKGROUND

Speciation genes contribute disproportionately to species divergence, but few examples exist, especially in vertebrates. Here we test whether Zan, which encodes the sperm acrosomal protein zonadhesin that mediates species-specific adhesion to the egg's zona pellucida, is a speciation gene in placental mammals.

RESULTS

Genomic ontogeny reveals that Zan arose by repurposing of a stem vertebrate gene that was lost in multiple lineages but retained in Eutheria on acquiring a function in egg recognition. A 112-species Zan sequence phylogeny, representing 17 of 19 placental Orders, resolves all species into monophyletic groups corresponding to recognized Orders and Suborders, with <5% unsupported nodes. Three other rapidly evolving germ cell genes (Adam2, Zp2, and Prm1), a paralogous somatic cell gene (TectA), and a mitochondrial gene commonly used for phylogenetic analyses (Cytb) all yield trees with poorer resolution than the Zan tree and inferior topologies relative to a widely accepted mammalian supertree. Zan divergence by intense positive selection produces dramatic species differences in the protein's properties, with ordinal divergence rates generally reflecting species richness of placental Orders consistent with expectations for a speciation gene that acts across a wide range of taxa. Furthermore, Zan's combined phylogenetic utility and divergence exceeds those of all other genes known to have evolved in Eutheria by positive selection, including the only other mammalian speciation gene, Prdm9.

CONCLUSIONS

Species-specific egg recognition conferred by Zan's functional divergence served as a mode of prezygotic reproductive isolation that promoted the extraordinary adaptive radiation and success of Eutheria.

Survival-related indicators ALOX12B and SPRR1A are associated with DNA damage repair and tumor microenvironment status in HPV 16-negative head and neck squamous cell carcinoma patients

OBJECTIVES

To investigate prognostic-related gene signature based on DNA damage repair and tumor microenvironment statue in human papillomavirus 16 negative (HPV16-) head and neck squamous cell carcinoma (HNSCC).

METHODS

For the RNA-sequence matrix in HPV16- HNSCC in the Cancer Genome Atlas (TCGA) cohort, the DNA damage response (DDR) and tumor microenvironment (TM) status of each patient sample was estimated by using the ssGSEA algorithm. Through bioinformatics analysis in DDR_high/TM_high (n = 311) and DDR_high/TM_low (n = 53) groups, a survival-related gene signature was selected in the TCGA cohort. Two independent external validation cohorts (GSE65858 (n = 210) and GSE41613 (n = 97)) with HPV16- HNSCC patients validated the gene signature. Correlations among the clinical-related hub differentially expressed genes (DEGs) and infiltrated immunocytes were explored with the TIMER2.0 server. Drug screening based on hub DEGs was performed using the CellMiner and GSCALite databases. The loss-of-function studies were used to evaluate the effect of screened survival-related gene on the motility of HPV- HNSCC cells in vitro.

RESULTS

A high DDR level (P = 0.025) and low TM score (P = 0.012) were independent risk factors for HPV16- HNSCC. Downregulated expression of ALOX12B or SPRR1A was associated with poor survival rate and advanced cancer stages. The pathway enrichment analysis showed the DDR_high/TM_low samples were enriched in glycosphingolipid biosynthesis-lacto and neolacto series, glutathione metabolism, platinum drug resistance, and ferroptosis pathways, while the DDR_high/TM_low samples were enriched in Th17 cell differentiation, Neutrophil extracellular trap formation, PD - L1 expression and PD - 1 checkpoint pathway in cancer. Notably, the expression of ALOX12B and SPRR1A were negatively correlated with cancer-associated fibroblasts (CAFs) infiltration and CAFs downstream effectors. Sensitivity to specific chemotherapy regimens can be derived from gene expressions. In addition, ALOX12B and SPRR1A expression was associated with the mRNA expression of insulin like growth factor 1 receptor (IGF1R), AKT serine/threonine kinase 1 (AKT1), mammalian target of rapamycin (MTOR), and eukaryotic translation initiation factor 4E binding protein 1 (EIF4EBP1) in HPV negative HNSCC. Down-regulation of ALOX12B promoted HPV- HNSCC cells migration and invasion in vitro.

CONCLUSIONS

ALOX12B and SPRR1A served as a gene signature for overall survival in HPV16- HNSCC patients, and correlated with the amount of infiltrated CAFs. The specific drug pattern was determined by the gene signature.

Evolutionary context of psoriatic immune skin response

The skin is vital for protecting the body and perceiving external stimuli in the environment. Ability to adapt between environments is in part based on skin phenotypic plasticity, indicating evolved homeostasis between skin and environment. This homeostasis reflects the greater relationship between the body and the environment, and disruptions in this balance may lead to accumulation of susceptibility factors for autoimmune conditions like psoriasis. In this study, we examined the relationship between rapid, lineage-specific evolution of human skin and formation of psoriatic skin responses at the transcriptome level. We collected skin tissue biopsies from individuals with psoriasis and compared gene expression in psoriatic plaques to non-plaque psoriatic skin. We then compared these data with non-psoriatic skin transcriptome data from multiple primate species. We found 67 genes showing human-specific skin expression that are also differentially regulated in psoriatic skin; these genes are significantly enriched for skin barrier function, immunity and neuronal development. We identified six gene clusters with differential expression in the context of human evolution and psoriasis, suggesting underlying regulatory mechanisms in these loci. Human and psoriasis-specific enrichment of neuroimmune genes shows the importance of the ongoing evolved homeostatic relationship between skin and external environment. These results have implications for both evolutionary medicine and public health, using transcriptomic data to acknowledge the importance of an individual’s surroundings on their overall health.

The skin is important for protecting the body from the environment and perceiving external stimuli, creating an evolved balance between skin and the environment. We compare skin gene expression in humans with psoriasis to humans and non-human primates without psoriasis to better understand human-specific evolutionary changes in the skin. Our results suggest important evolutionary links between skin perception, human-specific skin development and immune response.

The Epidermis: Redox Governor of Health and Diseases

A functional epithelial barrier necessitates protection against dehydration, and ichthyoses are caused by defects in maintaining the permeability barrier in the stratum corneum (SC), the uppermost protective layer composed of dead cells and secretory materials from the living layer stratum granulosum (SG). We have found that loricrin (LOR) is an essential effector of cornification that occurs in the uppermost layer of SG (SG1). LOR promotes the maturation of corneocytes and extracellular adhesion structure through organizing disulfide cross-linkages, albeit being dispensable for the SC permeability barrier. This review takes psoriasis and AD as the prototype of impaired cornification. Despite exhibiting immunological traits that oppose each other, both conditions share the epidermal differentiation complex as a susceptible locus. We also review recent mechanistic insights on skin diseases, focusing on the Kelch-like erythroid cell-derived protein with the cap "n" collar homology-associated protein 1/NFE2-related factor 2 signaling pathway, as they coordinate the epidermis-intrinsic xenobiotic metabolism. Finally, we refine the theoretical framework of thiol-mediated crosstalk between keratinocytes and leukocytes in the epidermis that was put forward earlier.

Seeing Keratinocyte Proteins through the Looking Glass of Intrinsic Disorder

Epidermal keratinocyte proteins include many with an eccentric amino acid content (compositional bias), atypical ultrastructural fate (built-in protease sensitivity), or assembly visible at the light microscope level (cytoplasmic granules). However, when considered through the looking glass of intrinsic disorder (ID), these apparent oddities seem quite expected. Keratinocyte proteins with highly repetitive motifs are of low complexity but high adaptation, providing polymers (e.g., profilaggrin) for proteolysis into bioactive derivatives, or monomers (e.g., loricrin) repeatedly cross-linked to self and other proteins to shield underlying tissue. Keratohyalin granules developing from liquid–liquid phase separation (LLPS) show that unique biomolecular condensates (BMC) and proteinaceous membraneless organelles (PMLO) occur in these highly customized cells. We conducted bioinformatic and in silico assessments of representative keratinocyte differentiation-dependent proteins. This was conducted in the context of them having demonstrated potential ID with the prospect of that characteristic driving formation of distinctive keratinocyte structures. Intriguingly, while ID is characteristic of many of these proteins, it does not appear to guarantee LLPS, nor is it required for incorporation into certain keratinocyte protein condensates. Further examination of keratinocyte-specific proteins will provide variations in the theme of PMLO, possibly recognizing new BMC for advancements in understanding intrinsically disordered proteins as reflected by keratinocyte biology.

The evolution of human skin pigmentation involved the interactions of genetic, environmental, and cultural variables

The primary biological role of human skin pigmentation is as a mediator of penetration of ultraviolet radiation (UVR) into the deep layers of skin and the cutaneous circulation. Since the origin of Homo sapiens, dark, protective constitutive pigmentation and strong tanning abilities have been favored under conditions of high UVR and represent the baseline condition for modern humans. The evolution of partly depigmented skin and variable tanning abilities has occurred multiple times in prehistory, as populations have dispersed into environments with lower and more seasonal UVR regimes, with unique complements of genes and cultural practices. The evolution of extremes of dark pigmentation and depigmentation has been rare and occurred only under conditions of extremely high or low environmental UVR, promoted by positive selection on variant pigmentation genes followed by limited gene flow. Over time, the evolution of human skin pigmentation has been influenced by the nature and course of human dispersals and modifications of cultural practices, which have modified the nature and actions of skin pigmentation genes. Throughout most of prehistory and history, the evolution of human skin pigmentation has been a contingent and non-deterministic process.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

Selective sweep for an enhancer involucrin allele identifies skin barrier adaptation out of Africa

The genetic modules that contribute to human evolution are poorly understood. Here we investigate positive selection in the Epidermal Differentiation Complex locus for skin barrier adaptation in diverse HapMap human populations (CEU, JPT/CHB, and YRI). Using Composite of Multiple Signals and iSAFE, we identify selective sweeps for LCE1A-SMCP and involucrin (IVL) haplotypes associated with human migration out-of-Africa, reaching near fixation in European populations. CEU-IVL is associated with increased IVL expression and a known epidermis-specific enhancer. CRISPR/Cas9 deletion of the orthologous mouse enhancer in vivo reveals a functional requirement for the enhancer to regulate Ivl expression in cis. Reporter assays confirm increased regulatory and additive enhancer effects of CEU-specific polymorphisms identified at predicted IRF1 and NFIC binding sites in the IVL enhancer (rs4845327) and its promoter (rs1854779). Together, our results identify a selective sweep for a cis regulatory module for CEU-IVL, highlighting human skin barrier evolution for increased IVL expression out-of-Africa.

The KEAP1/NRF2 Signaling Pathway in Keratinization

Keratinization is a tissue adaptation, but aberrant keratinization is associated with skin disorders such as ichthyoses, atopic dermatitis, psoriasis, and acne. The disease phenotype stems from the interaction between genes and the environment; therefore, an understanding of the adaptation machinery may lead to a new appreciation of pathomechanisms. The KEAP1/NRF2 signaling pathway mediates the environmental responses of squamous epithelial tissue. The unpredicted outcome of the Keap1-null mutation in mice allowed us to revisit the basic principle of the biological process of keratinization: sulfur metabolism establishes unparalleled cytoprotection in the body wall of terrestrial mammals. We summarize the recent understanding of the KEAP1/NRF2 signaling pathway, which is a thiol-based sensor-effector apparatus, with particular focuses on epidermal differentiation in the context of the gene-environment interaction, the structure/function principles involved in KEAP1/NRF2 signaling, lessons from mouse models, and their pathological implications. This synthesis may provide insights into keratinization, which provides physical insulation and constitutes an essential innate integumentary defense system.

Regulation of Filaggrin, Loricrin, and Involucrin by IL-4, IL-13, IL-17A, IL-22, AHR, and NRF2: Pathogenic Implications in Atopic Dermatitis

Atopic dermatitis (AD) is an eczematous, pruritic skin disorder with extensive barrier dysfunction and elevated interleukin (IL)-4 and IL-13 signatures. The barrier dysfunction correlates with the downregulation of barrier-related molecules such as filaggrin (FLG), loricrin (LOR), and involucrin (IVL). IL-4 and IL-13 potently inhibit the expression of these molecules by activating signal transducer and activator of transcription (STAT)6 and STAT3. In addition to IL-4 and IL-13, IL-22 and IL-17A are probably involved in the barrier dysfunction by inhibiting the expression of these barrier-related molecules. In contrast, natural or medicinal ligands for aryl hydrocarbon receptor (AHR) are potent upregulators of FLG, LOR, and IVL expression. As IL-4, IL-13, IL-22, and IL-17A are all capable of inducing oxidative stress, antioxidative AHR agonists such as coal tar, glyteer, and tapinarof exert particular therapeutic efficacy for AD. These antioxidative AHR ligands are known to activate an antioxidative transcription factor, nuclear factor E2-related factor 2 (NRF2). This article focuses on the mechanisms by which FLG, LOR, and IVL expression is regulated by IL-4, IL-13, IL-22, and IL-17A. The author also summarizes how AHR and NRF2 dual activators exert their beneficial effects in the treatment of AD.

The evolution of S100A7 in primates: a model of concerted and birth-and-death evolution

The human S100A7 resides in the epidermal differentiation complex (EDC) and has been described as a key effector of innate immunity. In humans, there are five S100A7 genes located in tandem-S100A7A, S100A7P1, S100AL2, S100A7, and S100AP2. The presence of several retroelements in the S100A7A/S100A7P1 and S100A7/S100A7P2 clusters suggests that these genes were originated from a duplication around ~ 35 million years ago, during or after the divergence of Platyrrhini and Catarrhini primates. To test this hypothesis, and taking advantage of the high number of genomic sequences available in the public databases, we retrieved S100A7 gene sequences of 12 primates belonging to the Cercopithecoidea and Hominoidea (Catarrhini species). Our results support the duplication theory, with at least one gene of each cluster being identified in both Cercopithecoidea and Hominoidea species. Moreover, given the presence of an ongoing gene conversion event between S100A7 and S100A7A, a high rate of mutation in S100A7L2 and the presence of pseudogenes, we proposed a model of concerted and birth-and-death evolution to explain the evolution of S100A7 gene family. Indeed, our results suggest that S100A7L2 most likely suffered a neofunctionalization in the Catarrhini group. Being S100A7 a major protein in innate defense, we believe that our findings could open new doors in the study of this gene family in immune system.

Recent evolution of the human skin barrier

The skin is the first line of defense against the environment, with the epidermis as the outermost tissue providing much of the barrier function. Given its direct exposure to and encounters with the environment, the epidermis must evolve to provide an optimal barrier for the survival of an organism. Recent advances in genomics have identified a number of genes for the human skin barrier that have undergone evolutionary changes since humans diverged from chimpanzees. Here, we highlight a selection of key and innovative genetic findings for skin barrier evolution in our divergence from our primate ancestors and among modern human populations.

Tinkering and the Origins of Heritable Anatomical Variation in Vertebrates

Evolutionary change comes from natural and other forms of selection acting on existing anatomical and physiological variants. While much is known about selection, little is known about the details of how genetic mutation leads to the range of heritable anatomical variants that are present within any population. This paper takes a systems-based view to explore how genomic mutation in vertebrate genomes works its way upwards, though changes to proteins, protein networks, and cell phenotypes to produce variants in anatomical detail. The evidence used in this approach mainly derives from analysing anatomical change in adult vertebrates and the protein networks that drive tissue formation in embryos. The former indicate which processes drive variation-these are mainly patterning, timing, and growth-and the latter their molecular basis. The paper then examines the effects of mutation and genetic drift on these processes, the nature of the resulting heritable phenotypic variation within a population, and the experimental evidence on the speed with which new variants can appear under selection. The discussion considers whether this speed is adequate to explain the observed rate of evolutionary change or whether other non-canonical, adaptive mechanisms of heritable mutation are needed. The evidence to hand suggests that they are not, for vertebrate evolution at least.

Positive selection drives the evolution of endocrine regulatory bone morphogenetic protein system in mammals

The rapid evolution of reproductive proteins might be driven by positive Darwinian selection. The bone morphogenetic protein family is the largest within the transforming growth factor (TGF) superfamily. A little have been known about the molecular evolution of bone morphogenetic proteins exhibiting potential role in mammalian reproduction. In this study we investigated mammalian bone morphogenetic proteins using maximum likelihood approaches of codon substitutions to identify positive Darwinian selection in various species. The proportion of positively selected sites was tested by different likelihood models for individual codon, and M8 were found to be the best model. The percentage of positively elected sites under M8 are 2.20% with ω = 1.089 for BMP2, 1.6% with ω = 1.61 for BMP 4 0.53% for BMP15 with ω = 1.56 and 0.78% for GDF9 with ω = 1.93. The percentage of estimated selection sites under M8 is strong statistical confirmation that divergence of bone morphogenetic proteins is driven by Darwinian selection. For the proteins, model M8 was found significant for all proteins with ω > 1. To further test positive selection on particular amino acids, the evolutionary conservation of amino acid were measured based on phylogenetic linkage among sequences. For exploring the impact of these somatic substitution mutations in the selection region on human cancer, we identified one pathogenic mutation in human BMP4 and one in BMP15, possibly causing prostate cancer and six neutral mutations in BMPs. The comprehensive map of selection results allows the researchers to perform systematic approaches to detect the evolutionary footprints of selection on specific gene in specific species.