Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation.

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

The Unreasonable Effectiveness of Reaction Diffusion in Vertebrate Skin Color Patterning.

In 1952, Alan Turing published the reaction-diffusion (RD) mathematical framework, laying the foundations of morphogenesis as a self-organized process emerging from physicochemical first principles. Regrettably, this approach has been widely doubted in the field of developmental biology. First, we summarize Turing's line of thoughts to alleviate the misconception that RD is an artificial mathematical construct. Second, we discuss why phenomenological RD models are particularly effective for understanding skin color patterning at the meso/macroscopic scales, without the need to parameterize the profusion of variables at lower scales. More specifically, we discuss how RD models (a) recapitulate the diversity of actual skin patterns, (b) capture the underlying dynamics of cellular interactions, (c) interact with tissue size and shape, (d) can lead to ordered sequential patterning, (e) generate cellular automaton dynamics in lizards and snakes, (f) predict actual patterns beyond their statistical features, and (g) are robust to model variations. Third, we discuss the utility of linear stability analysis and perform numerical simulations to demonstrate how deterministic RD emerges from the underlying chaotic microscopic agents.

NNT mediates redox-dependent pigmentation via a UVB- and MITF-independent mechanism.

Ultraviolet (UV) light and incompletely understood genetic and epigenetic variations determine skin color. Here we describe an UV- and microphthalmia-associated transcription factor (MITF)-independent mechanism of skin pigmentation. Targeting the mitochondrial redox-regulating enzyme nicotinamide nucleotide transhydrogenase (NNT) resulted in cellular redox changes that affect tyrosinase degradation. These changes regulate melanosome maturation and, consequently, eumelanin levels and pigmentation. Topical application of small-molecule inhibitors yielded skin darkening in human skin, and mice with decreased NNT function displayed increased pigmentation. Additionally, genetic modification of NNT in zebrafish alters melanocytic pigmentation. Analysis of four diverse human cohorts revealed significant associations of skin color, tanning, and sun protection use with various single-nucleotide polymorphisms within NNT. NNT levels were independent of UVB irradiation and redox modulation. Individuals with postinflammatory hyperpigmentation or lentigines displayed decreased skin NNT levels, suggesting an NNT-driven, redox-dependent pigmentation mechanism that can be targeted with NNT-modifying topical drugs for medical and cosmetic purposes.

An Unexpectedly Complex Architecture for Skin Pigmentation in Africans.

Approximately 15 genes have been directly associated with skin pigmentation variation in humans, leading to its characterization as a relatively simple trait. However, by assembling a global survey of quantitative skin pigmentation phenotypes, we demonstrate that pigmentation is more complex than previously assumed, with genetic architecture varying by latitude. We investigate polygenicity in the KhoeSan populations indigenous to southern Africa who have considerably lighter skin than equatorial Africans. We demonstrate that skin pigmentation is highly heritable, but known pigmentation loci explain only a small fraction of the variance. Rather, baseline skin pigmentation is a complex, polygenic trait in the KhoeSan. Despite this, we identify canonical and non-canonical skin pigmentation loci, including near SLC24A5, TYRP1, SMARCA2/VLDLR, and SNX13, using a genome-wide association approach complemented by targeted resequencing. By considering diverse, under-studied African populations, we show how the architecture of skin pigmentation can vary across humans subject to different local evolutionary pressures.

DNA Phenotyping: Snapshot of a Criminal.

Recent progresses in genomic technologies provide both opportunities and challenges to forensics.

The dynamics of pattern matching in camouflaging cuttlefish.

Many cephalopods escape detection using camouflage1. This behaviour relies on a visual assessment of the surroundings, on an interpretation of visual-texture statistics2-4 and on matching these statistics using millions of skin chromatophores that are controlled by motoneurons located in the brain5-7. Analysis of cuttlefish images proposed that camouflage patterns are low dimensional and categorizable into three pattern classes, built from a small repertoire of components8-11. Behavioural experiments also indicated that, although camouflage requires vision, its execution does not require feedback5,12,13, suggesting that motion within skin-pattern space is stereotyped and lacks the possibility of correction. Here, using quantitative methods14, we studied camouflage in the cuttlefish Sepia officinalis as behavioural motion towards background matching in skin-pattern space. An analysis of hundreds of thousands of images over natural and artificial backgrounds revealed that the space of skin patterns is high-dimensional and that pattern matching is not stereotyped-each search meanders through skin-pattern space, decelerating and accelerating repeatedly before stabilizing. Chromatophores could be grouped into pattern components on the basis of their covariation during camouflaging. These components varied in shapes and sizes, and overlay one another. However, their identities varied even across transitions between identical skin-pattern pairs, indicating flexibility of implementation and absence of stereotypy. Components could also be differentiated by their sensitivity to spatial frequency. Finally, we compared camouflage to blanching, a skin-lightening reaction to threatening stimuli. Pattern motion during blanching was direct and fast, consistent with open-loop motion in low-dimensional pattern space, in contrast to that observed during camouflage.

Distinct melanocyte subpopulations defined by stochastic expression of proliferation or maturation programs enable a rapid and sustainable pigmentation response.

The ultraviolet (UV) radiation triggers a pigmentation response in human skin, wherein, melanocytes rapidly activate divergent maturation and proliferation programs. Using single-cell sequencing, we demonstrate that these 2 programs are segregated in distinct subpopulations in melanocytes of human and zebrafish skin. The coexistence of these 2 cell states in cultured melanocytes suggests possible cell autonomy. Luria-Delbrück fluctuation test reveals that the initial establishment of these states is stochastic. Tracking of pigmenting cells ascertains that the stochastically acquired state is faithfully propagated in the progeny. A systemic approach combining single-cell multi-omics (RNA+ATAC) coupled to enhancer mapping with H3K27 acetylation successfully identified state-specific transcriptional networks. This comprehensive analysis led to the construction of a gene regulatory network (GRN) that under the influence of noise, establishes a bistable system of pigmentation and proliferation at the population level. This GRN recapitulates melanocyte behaviour in response to external cues that reinforce either of the states. Our work highlights that inherent stochasticity within melanocytes establishes dedicated states, and the mature state is sustained by selective enhancers mark through histone acetylation. While the initial cue triggers a proliferation response, the continued signal activates and maintains the pigmenting subpopulation via epigenetic imprinting. Thereby our study provides the basis of coexistence of distinct populations which ensures effective pigmentation response while preserving the self-renewal capacity.

The dynamic behavior of chromatophores marks the transition from bands to spots in leopard geckos.

Reptilian skin coloration is spectacular and diverse, yet little is known about the ontogenetic processes that govern its establishment and the molecular signaling pathways that determine it. Here, we focus on the development of the banded pattern of leopard gecko hatchlings and the transition to black spots in the adult. With our histological analyses, we show that iridophores are present in the white and yellow bands of the hatchling and they gradually perish in the adult skin. Furthermore, we demonstrate that melanophores can autonomously form spots in the absence of the other chromatophores both on the regenerated skin of the tail and on the dorsal skin of the Mack Super Snow (MSS) leopard geckos. This color morph is characterized by uniform black coloration in hatchlings and black spots in adulthood; we establish that their skin is devoid of xanthophores and iridophores at both stages. Our genetic analyses identified a 13-nucleotide deletion in the PAX7 transcription factor of MSS geckos, affecting its protein coding sequence. With our single-cell transcriptomics analysis of embryonic skin, we confirm that PAX7 is expressed in iridophores and xanthophores, suggesting that it plays a key role in the differentiation of both chromatophores. Our in situ hybridizations on whole-mount embryos document the dynamics of the skin pattern formation and how it is impacted in the PAX7 mutants. We hypothesize that the melanophores-iridophores interactions give rise to the banded pattern of the hatchlings and black spot formation is an intrinsic capacity of melanophores in the postembryonic skin.

UV-Protection Timer Controls Linkage between Stress and Pigmentation Skin Protection Systems.

Skin sun exposure induces two protection programs: stress responses and pigmentation, the former within minutes and the latter only hours afterward. Although serving the same physiological purpose, it is not known whether and how these programs are coordinated. Here, we report that UVB exposure every other day induces significantly more skin pigmentation than the higher frequency of daily exposure, without an associated increase in stress responses. Using mathematical modeling and empirical studies, we show that the melanocyte master regulator, MITF, serves to synchronize stress responses and pigmentation and, furthermore, functions as a UV-protection timer via damped oscillatory dynamics, thereby conferring a trade-off between the two programs. MITF oscillations are controlled by multiple negative regulatory loops, one at the transcriptional level involving HIF1α and another post-transcriptional loop involving microRNA-148a. These findings support trait linkage between the two skin protection programs, which, we speculate, arose during furless skin evolution to minimize skin damage.

Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation.

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

Genetic Basis and Evolutionary Forces of Sexually Dimorphic Color Variation in a Toad-Headed Agamid Lizard.

In the animal kingdom, sexually dimorphic color variation is a widespread phenomenon that significantly influences survival and reproductive success. However, the genetic underpinnings of this variation remain inadequately understood. Our investigation into sexually dimorphic color variation in the desert-dwelling Guinan population of the toad-headed agamid lizard (Phrynocephalus putjatai) utilized a multidisciplinary approach, encompassing phenotypic, ultrastructural, biochemical, genomic analyses, and behavioral experiments. Our findings unveil the association between distinct skin colorations and varying levels of carotenoid and pteridine pigments. The red coloration in males is determined by a genomic region on chromosome 14, housing four pigmentation genes: BCO2 and three 6-pyruvoyltetrahydropterin synthases. A Guinan population-specific nonsynonymous single nucleotide polymorphism in BCO2 is predicted to alter the electrostatic potential within the binding domain of the BCO2-ß-carotene complex, influencing their interaction. Additionally, the gene MAP7 on chromosome 2 emerges as a potential contributor to the blue coloration in subadults and adult females. Sex-specific expression patterns point to steroid hormone-associated genes (SULT2B1 and SRD5A2) as potential upstream regulators influencing sexually dimorphic coloration. Visual modeling and field experiments support the potential selective advantages of vibrant coloration in desert environments. This implies that natural selection, potentially coupled with assortative mating, might have played a role in fixing color alleles, contributing to prevalence in the local desert habitat. This study provides novel insights into the genetic basis of carotenoid and pteridine-based color variation, shedding light on the evolution of sexually dimorphic coloration in animals. Moreover, it advances our understanding of the driving forces behind such intricate coloration patterns.

Direct detection of natural selection in Bronze Age Britain.

We developed a novel method for efficiently estimating time-varying selection coefficients from genome-wide ancient DNA data. In simulations, our method accurately recovers selective trajectories and is robust to misspecification of population size. We applied it to a large data set of ancient and present-day human genomes from Britain and identified seven loci with genome-wide significant evidence of selection in the past 4500 yr. Almost all of them can be related to increased vitamin D or calcium levels, suggesting strong selective pressure on these or related phenotypes. However, the strength of selection on individual loci varied substantially over time, suggesting that cultural or environmental factors moderated the genetic response. Of 28 complex anthropometric and metabolic traits, skin pigmentation was the only one with significant evidence of polygenic selection, further underscoring the importance of phenotypes related to vitamin D. Our approach illustrates the power of ancient DNA to characterize selection in human populations and illuminates the recent evolutionary history of Britain.

Lrig1- and Wnt-dependent niches dictate segregation of resident immune cells and melanocytes in murine tail epidermis.

The barrier-forming, self-renewing mammalian epidermis comprises keratinocytes, pigment-producing melanocytes and resident immune cells as first-line host defense. In murine tail skin, interfollicular epidermis patterns into pigmented 'scale' and hypopigmented 'interscale' epidermis. Why and how mature melanocytes accumulate in scale epidermis is unresolved. Here, we delineate a cellular hierarchy among epidermal cell types that determines skin patterning. Already during postnatal development, melanocytes co-segregate with newly forming scale compartments. Intriguingly, this process coincides with partitioning of both Langerhans cells and dendritic epidermal T cells to interscale epidermis, suggesting functional segregation of pigmentation and immune surveillance. Analysis of non-pigmented mice and of mice lacking melanocytes or resident immune cells revealed that immunocyte patterning is melanocyte and melanin independent and, vice versa, immune cells do not control melanocyte localization. Instead, genetically enforced progressive scale fusion upon Lrig1 deletion showed that melanocytes and immune cells dynamically follow epithelial scale:interscale patterns. Importantly, disrupting Wnt-Lef1 function in keratinocytes caused melanocyte mislocalization to interscale epidermis, implicating canonical Wnt signaling in organizing the pigmentation pattern. Together, this work uncovers cellular and molecular principles underlying the compartmentalization of tissue functions in skin.

Implementing differentially pigmented skin models for predicting drug response variability across human ancestries.

Persistent racial disparities in health outcomes have catalyzed legislative reforms and heightened scientific focus recently. However, despite the well-documented properties of skin pigments in binding drug compounds, their impact on therapeutic efficacy and adverse drug responses remains insufficiently explored. This perspective examines the intricate relationships between variation in melanin-based skin pigmentation and pharmacokinetics and -dynamics, highlighting the need for considering diversity in skin pigmentation as a variable to advance the equitability of pharmacological interventions. The article provides guidelines on the selection of New Approach Methods (NAMs) to foster inclusive study designs in preclinical drug development pipelines, leading to an improved level of translatability to the clinic.

Temporal analysis of melanogenesis identifies fatty acid metabolism as key skin pigment regulator.

Therapeutic methods to modulate skin pigmentation has important implications for skin cancer prevention and for treating cutaneous hyperpigmentary conditions. Towards defining new potential targets, we followed temporal dynamics of melanogenesis using a cell-autonomous pigmentation model. Our study elucidates 3 dominant phases of synchronized metabolic and transcriptional reprogramming. The melanogenic trigger is associated with high MITF levels along with rapid uptake of glucose. The transition to pigmented state is accompanied by increased glucose channelisation to anabolic pathways that support melanosome biogenesis. SREBF1-mediated up-regulation of fatty acid synthesis results in a transient accumulation of lipid droplets and enhancement of fatty acids oxidation through mitochondrial respiration. While this heightened bioenergetic activity is important to sustain melanogenesis, it impairs mitochondria lately, shifting the metabolism towards glycolysis. This recovery phase is accompanied by activation of the NRF2 detoxication pathway. Finally, we show that inhibitors of lipid metabolism can resolve hyperpigmentary conditions in a guinea pig UV-tanning model. Our study reveals rewiring of the metabolic circuit during melanogenesis, and fatty acid metabolism as a potential therapeutic target in a variety of cutaneous diseases manifesting hyperpigmentary phenotype.

Imprinting sets the stage for speciation.

Sexual imprinting-a phenomenon in which offspring learn parental traits and later use them as a model for their own mate preferences-can generate reproductive barriers between species1. When the target of imprinting is a mating trait that differs among young lineages, imprinted preferences may contribute to behavioural isolation and facilitate speciation1,2. However, in most models of speciation by sexual selection, divergent natural selection is also required; the latter acts to generate and maintain variation in the sexually selected trait or traits, and in the mating preferences that act upon them3. Here we demonstrate that imprinting, in addition to mediating female mate preferences, can shape biases in male-male aggression. These biases can act similarly to natural selection to maintain variation in traits and mate preferences, which facilitates reproductive isolation driven entirely by sexual selection. Using a cross-fostering study, we show that both male and female strawberry poison frogs (Oophaga pumilio) imprint on coloration, which is a mating trait that has diverged recently and rapidly in this species4. Cross-fostered females prefer to court mates of the same colour as their foster mother, and cross-fostered males are more aggressive towards rivals that share the colour of their foster mother. We also use a simple population-genetics model to demonstrate that when both male aggression biases and female mate preferences are formed through parental imprinting, sexual selection alone can (1) stabilize a sympatric polymorphism and (2) strengthen the trait-preference association that leads to behavioural reproductive isolation. Our study provides evidence of imprinting in an amphibian and suggests that this rarely considered combination of rival and sexual imprinting can reduce gene flow between individuals that bear divergent mating traits, which sets the stage for speciation by sexual selection.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice.

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Genetic adaptation of skin pigmentation in highland Tibetans.

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Rapid radiation in a highly diverse marine environment.

Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.