Pigmentation Phenotypes Research Articles

A genetic screen of transcription factors in the Drosophila melanogaster abdomen identifies novel pigmentation genes.

Gene regulatory networks specify the gene expression patterns needed for traits to develop. Differences in these networks can result in phenotypic differences between organisms. Although loss-of-function genetic screens can identify genes necessary for trait formation, gain-of-function screens can overcome genetic redundancy and identify loci whose expression is sufficient to alter trait formation. Here, we leveraged transgenic lines from the Transgenic RNAi Project at Harvard Medical School to perform both gain- and loss-of-function CRISPR/Cas9 screens for abdominal pigmentation phenotypes. We identified measurable effects on pigmentation patterns in the Drosophila melanogaster abdomen for 21 of 55 transcription factors in gain-of-function experiments and 7 of 16 tested by loss-of-function experiments. These included well-characterized pigmentation genes, such as bab1 and dsx, and transcription factors that had no known role in pigmentation, such as slp2. Finally, this screen was partially conducted by undergraduate students in a Genetics Laboratory course during the spring semesters of 2021 and 2022. We found this screen to be a successful model for student engagement in research in an undergraduate laboratory course that can be readily adapted to evaluate the effect of hundreds of genes on many different Drosophila traits, with minimal resources.

Novel mechanisms of MITF regulation identified in a mouse suppressor screen

MITF, a basic Helix-Loop-Helix Zipper (bHLHZip) transcription factor, plays vital roles in melanocyte development and functions as an oncogene. We perform a genetic screen for suppressors of the Mitf-associated pigmentation phenotype in mice and identify an intragenic Mitf mutation that terminates MITF at the K316 SUMOylation site, leading to loss of the C-end intrinsically disordered region (IDR). The resulting protein is more nuclear but less stable than wild-type MITF and retains DNA-binding ability. As a dimer, it can translocate wild-type and mutant MITF partners into the nucleus, improving its own stability thus ensuring nuclear MITF supply. smFRET analysis shows interactions between K316 SUMOylation and S409 phosphorylation sites across monomers; these interactions largely explain the observed effects. The recurrent melanoma-associated E318K mutation in MITF, which affects K316 SUMOylation, also alters protein regulation in concert with S409. This suggests that residues K316 and S409 of MITF are impacted by SUMOylation and phosphorylation, respectively, mediating effects on nuclear localization and stability through conformational changes. Our work provides a novel mechanism of genetic suppression, and an example of how apparently deleterious mutations lead to normal phenotypes.

The chloroplast protease system degrades stromal DUF760-1 and DUF760-2 domain-containing proteins at different rates.

The chloroplast chaperone CLPC1 aids to select, unfold, and deliver hundreds of proteins to the CLP protease for degradation. Through in vivo CLPC1, trapping we previously identified dozens of proteins that are (potential) substrate adaptors or substrates for the CLP chaperone-protease system. In this study, we show that two of these highly trapped proteins, DUF760-1 and DUF760-2, are substrates for the CLP protease in Arabidopsis (Arabidopsis thaliana). Loss-of-function mutants and transgenic plants were created for phenotyping, protein expression, and localization using immunoblotting and confocal microscopy. In planta BiFC, cycloheximide chase assays, and yeast 2-hybrid analyses were conducted to determine protein interactions and protein half-life. Both DUF760 proteins directly interacted with the N-domain of CLPC1 and both were highly enriched in clpc1-1 and clpr2-1 mutants. Accordingly, in vivo cycloheximide chase assays demonstrated that both DUF760 proteins are degraded by the CLP protease. The half-life of DUF760-1 was 4 to 6 h, whereas DUF760-2 was highly unstable and difficult to detect unless CLP proteolysis was inhibited. Null mutants for DUF760-1 and DUF760-2 showed weak but differential pigment phenotypes and differential sensitivity to protein translation inhibitors. This study demonstrates that DUF760-1 and DUF760-2 are substrates of the CLP chaperone-protease system and excellent candidates for the determination of CLP substrate degrons.

The genetic basis of divergent melanic pigmentation in benthic and limnetic threespine stickleback

Pigmentation is an excellent trait to examine patterns of evolutionary change because it is often under natural selection. Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) exhibit distinct pigmentation phenotypes, likely an adaptation to occupation of divergent niches. The genetic architecture of pigmentation in vertebrates appears to be complex. Prior QTL mapping of threespine stickleback pigmentation phenotypes has identified several candidate loci. However—relative to other morphological phenotypes (e.g., spines or lateral plates)—the genetic architecture of threespine stickleback pigmentation remains understudied. Here, we performed QTL mapping for two melanic pigmentation traits (melanophore density and lateral barring) using benthic-limnetic F2 crosses. The two traits mapped to different chromosomes, suggesting a distinct genetic basis. The resulting QTLs were additive, but explained a relatively small fraction of the total variance (~6%). QTLs maps differed by F1 family, suggesting variation in genetic architecture or ability to detect loci of small effect. Functional analysis identified enriched pathways for candidate loci. Several of the resulting candidate loci for pigmentation, including three loci in enriched pathways (bco1, sulf1, and tyms) have been previously indicated to affect pigmentation in other vertebrates. These findings add to a growing body of evidence suggesting pigmentation is often polygenic.

Journey through the spectacular landscape of melanocortin 1 receptor.

The physiological role of α-melanocyte stimulating hormone in regulating integumental pigmentation of many vertebrate species has been recognized since the 1960's. However, its physiological significance for human pigmentation remained enigmatic until the 1990's. α-Melanocyte stimulating hormone and related melanocortins are synthesized locally in the skin, primarily by keratinocytes, in addition to the pituitary gland, and therefore act as paracrine factors for melanocytes. Human melanocytes express the melanocortin 1 receptor, which recognizes α-melanocyte stimulating hormone and the related adrenocorticotropic hormone as agonists. This review summarizes the current knowledge of the pleotropic effects of the activated melanocortin 1 receptor that maintain human melanocyte homeostasis by regulating melanogenesis and the response to environmental stressors, mainly solar radiation. Certain allelic variants of the melanocortin 1 receptor gene are associated with specific pigmentary phenotypes in various human populations. Variants associated with red hair phenotype compromise the function of the encoded receptor. Activation of the human melanocortin 1 receptor regulates eumelanin synthesis and enhances DNA damage response of melanocytes to solar radiation and oxidative stressors. We describe how synthetic selective melanocortin 1 receptor agonists can be efficacious as sunless tanning agents, for treatment of vitiligo and photosensitivity disorders, and for prevention of skin cancer, including melanoma.

Mutation of hmgA, encoding homogentisate 1,2-dioxygenase, is responsible for pyomelanin production but does not impact the virulence of Burkholderia cenocepacia in a chronic granulomatous disease mouse lung infection.

The Burkholderia cepacia complex (Bcc) is a group of Gram-negative opportunistic bacteria often associated with fatal pulmonary infections in patients with impaired immunity, particularly those with cystic fibrosis (CF) and chronic granulomatous disease (CGD). Some Bcc strains are known to naturally produce pyomelanin, a brown melanin-like pigment known for scavenging free radicals; pigment production has been reported to enable Bcc strains to overcome the host cell oxidative burst. In this work, we investigated the role of pyomelanin in resistance to oxidative stress and virulence in strains J2315 and K56-2, two epidemic CF isolates belonging to the Burkholderia cenocepacia ET-12 lineage. We previously reported that a single amino acid change from glycine to arginine at residue 378 in homogentisate 1,2-dioxygenase (HmgA) affects the pigment production phenotype: pigmented J2315 has an arginine at position 378, while non-pigmented K56-2 has a glycine at this position. Herein, we performed allelic exchange to generate isogenic non-pigmented and pigmented strains of J2315 and K56-2, respectively, and tested these to determine whether pyomelanin contributes to the protection against oxidative stress in vitro as well as in a respiratory infection in CGD mice in vivo. Our results indicate that the altered pigment phenotype does not significantly impact these strains' ability to resist oxidative stress with H2O2 and NO in vitro and did not change the virulence and infection outcome in CGD mice in vivo suggesting that other factors besides pyomelanin are contributing to the pathophysiology of these strains.IMPORTANCEThe Burkholderia cepacia complex (Bcc) is a group of Gram-negative opportunistic bacteria that are often associated with fatal pulmonary infections in patients with impaired immunity, particularly those with cystic fibrosis and chronic granulomatous disease (CGD). Some Bcc strains are known to naturally produce pyomelanin, a brown melanin-like pigment known for scavenging free radicals and overcoming the host cell oxidative burst. We investigated the role of pyomelanin in Burkholderia cenocepacia strains J2315 (pigmented) and K56-2 (non-pigmented) and performed allelic exchange to generate isogenic non-pigmented and pigmented strains, respectively. Our results indicate that the altered pigment phenotype does not significantly impact these strains' ability to resist H2O2 or NO in vitro and did not alter the outcome of a respiratory infection in CGD mice in vivo. These results suggest that pyomelanin may not always constitute a virulence factor and suggest that other features are contributing to the pathophysiology of these strains.

Abstract 7327: Investigating the cell-type specific regulation of IRF4 and its role in lung cancer via a lung cancer risk-associated pleiotropic variant

Abstract Although smoking is a main risk factor for lung cancer, genetic factors also contribute to the risk, as shown by over 50 genomic loci identified by genome-wide association studies (GWAS). A recent cross-ancestry lung cancer GWAS identified a common germline variant, rs12203592 (C/T), in chr6p25.3 as a new signal, and quantitative trait loci (eQTL) colocalization analysis using lung and blood tissues (GTEx v8) identified IRF4 as a likely target of rs12203592. rs12203592 has previously been linked with diverse traits, including pigmentation phenotypes, smoking cessation, lymphoblastic leukemia, melanoma, and other skin cancers. rs12203592 also appears to have varied allelic effects in different cell types. Namely, the lung cancer risk-associated T allele is correlated with higher IRF4 expression in lung tissue and blood cells (GTEx v8) but lower expression in primary melanocytes (n = 106; Zhang et al). Since previous studies established the melanocyte-specific allelic function of this SNP through lineage-specific transcription factors (TFs), we hypothesized that cell-type-specific TFs mediate the observed allelic effect in different tissue/cell types and further contribute to the unique role of IRF4 in lung tumorigenesis. To identify these potential cell-type-specific regulators, we performed mass spectrometry using the nuclear extracts from three cell lines representing lung, blood, and melanocyte lineage. The results nominated the NF-kB and ZEB protein families as the most prominent T-specific binders in the A549 lung cancer cell line, which was validated by electromobility shift assays. Consistent with these findings, tissue-based SNP-gene interaction analysis demonstrated that the mRNA levels of NF-kB and ZEB family members (4 of 7 genes) were correlated with IRF4 expression with a significant interaction with rs12203592 in GTEx lung tissue (P < 0.017), where the correlation is mainly observed with the T allele. Notably, the SNP interaction was less pronounced in whole blood or melanocytes, and the correlation of these TFs and IRF4 expression was mainly observed with the C allele in melanocytes. To begin to establish the roles of these TFs and IRF4 in lung tumorigenesis and given that IRF4 was shown to promote endogenous DNA damage and proliferation of lung-lineage cell lines, we will investigate the roles of NF-kB and ZEB in IRF4 regulation in these contexts using knockdown and overexpression of these TFs. To further investigate lung cancer-specific IRF4 function, we will identify potential targets of IRF4 among lung cancer-associated functional variants established through our previous studies. This data will help elucidate the context-specific roles of IRF4 in lung cancer development and provide clues to gene regulation in a pleiotropic GWAS locus. Citation Format: Alexander Kane, Jinhu Yin, Erping Long, James Feng, Cathrin Gräwe, Harsh Patel, Jinyoung Byun, Christopher Amos, Kevin Brown, Jun Xia, Michiel Vermeulen, Jiyeon Choi. Investigating the cell-type specific regulation of IRF4 and its role in lung cancer via a lung cancer risk-associated pleiotropic variant [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7327.

AcMYB10 Involved in Anthocyanin Regulation of 'Hongyang' Kiwifruit Induced via Fruit Bagging and High-Postharvest-Temperature Treatments.

Light and temperature are key factors influencing the accumulation of anthocyanin in fruit crops. To assess the effects of fruit bagging during development and high post-ripening temperature on 'Hongyang' kiwifruit, we compared the pigmentation phenotypes and expression levels of anthocyanin-related genes between bagged and unbagged treatments, and between 25 °C and 37 °C postharvest storage temperatures. Both the bagging and 25 °C treatments showed better pigmentation phenotypes with higher anthocyanin concentrations. The results of the qRT-PCR analysis revealed that the gene expression levels of LDOX (leucoanthocyanidin dioxygenase), F3GT (UDP-flavonoid 3-O-glycosyltransferase ), AcMYB10, and AcbHLH42 were strongly correlated and upregulated by both the bagging treatment and 25 °C storage. The results of bimolecular fluorescence complementation and luciferase complementation imaging assays indicated an interaction between AcMYB10 and AcbHLH42 in plant cells, whereas the results of a yeast one-hybrid assay further demonstrated that AcMYB10 activated the promoters of AcLODX and AcF3GT. These results strongly suggest that enhanced anthocyanin synthesis is caused by the promoted expression of AcLODX and AcF3GT, regulated by the complex formed by AcMYB10-AcbHLH42.

Application of a 24-SNP Multiplex Genotyping Assay System for Phenotypic Identification of Fujian Han Population.

This study aimed to evaluate the utility of 24 single nucleotide polymorphism (SNP) loci associated with iris color and hair color in phenotypic identification of the Han Chinese population in Fujian Province. The selected SNPs, known for their strong correlation with specific human phenotypic features, provide valuable reference data for developing a molecular phenotypic identification system. A multiplex genotyping assay system was established with primers for the 24 SNPs linked to iris color and hair color synthesized based on existing literature. In total, 235 unrelated individuals of Han Chinese ethnicity in Fujian Province were included in this study. PowerStats v12 was employed to calculate forensic parameters associated with the 24 SNP loci, including gene frequencies, genotype frequencies, minor allele frequencies, discrimination power (DP), polymorphism information content (PIC), and observed heterozygosity (Ho). Hardy-Weinberg equilibrium tests were conducted for each locus. The SNP genotyping results were uploaded to the HIrisPlex model (https://HIrisPlex.erasmusmc.nl/) to predict iris and hair colors, and the inferred results were compared with manually assessed images. The accuracy of pigment phenotype inference was evaluated by using ROC curves in SPSS 26.0 software. The accuracy rates of inferring brown iris and black hair phenotypes were 99.6% and 99.5%. The area under the curve (AUC) values were 0.923 and 0.980, respectively. The 24 SNP loci demonstrated high accuracy in inferring iris color and hair color; it seems to be a useful tool for forensic phenotypic identification and anthropological or evolutionary applications. Establishment of suitable pigment classification criteria and optimized prediction models is based on revealing more phenotypic genetic markers.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon.

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

Disrupting the Repeat Domain of Premelanosome Protein (PMEL) Produces Dysamyloidosis and Dystrophic Ocular Pigment Reflective of Pigmentary Glaucoma.

Pigmentary glaucoma has recently been associated with missense mutations in PMEL that are dominantly inherited and enriched in the protein's fascinating repeat domain. PMEL pathobiology is intriguing because PMEL forms functional amyloid in healthy eyes, and this PMEL amyloid acts to scaffold melanin deposition. This is an informative contradistinction to prominent neurodegenerative diseases where amyloid formation is neurotoxic and mutations cause a toxic gain of function called "amyloidosis". Preclinical animal models have failed to model this PMEL "dysamyloidosis" pathomechanism and instead cause recessively inherited ocular pigment defects via PMEL loss of function; they have not addressed the consequences of disrupting PMEL's repetitive region. Here, we use CRISPR to engineer a small in-frame mutation in the zebrafish homolog of PMEL that is predicted to subtly disrupt the protein's repetitive region. Homozygous mutant larvae displayed pigmentation phenotypes and altered eye morphogenesis similar to presumptive null larvae. Heterozygous mutants had disrupted eye morphogenesis and disrupted pigment deposition in their retinal melanosomes. The deficits in the pigment deposition of these young adult fish were not accompanied by any detectable glaucomatous changes in intraocular pressure or retinal morphology. Overall, the data provide important in vivo validation that subtle PMEL mutations can cause a dominantly inherited pigment pathology that aligns with the inheritance of pigmentary glaucoma patient pedigrees. These in vivo observations help to resolve controversy regarding the necessity of PMEL's repeat domain in pigmentation. The data foster an ongoing interest in an antithetical dysamyloidosis mechanism that, akin to the amyloidosis of devastating dementias, manifests as a slow progressive neurodegenerative disease.

Rete ridges are decreased in dyschromic burn hypertrophic scar: A histological study

Dyschromic hypertrophic scar (HTS) is a common sequelae of burn injury, however, its mechanism has not been elucidated. This work is a histological study of these scars with a focus on rete ridges. Rete ridges are important for normal skin physiology, and their absence or presence may hold mechanistic significance in post-burn HTS dyschromia. It was posited that hyper-, and hypo-pigmented areas of scars have different numbers of rete ridges. Subjects with dyschromic burn hypertrophic scar were prospectively enrolled (n = 44). Punch biopsies of hyper-, hypo-, and normally pigmented scar and skin were collected. Biopsies were paraffin embedded, sectioned, stained with H&E, and imaged. The number of rete ridges were investigated. Burn hypertrophic scars that healed without autografts were first investigated. The number of rete ridges was higher in normal skin compared to HTS that was either hypo- (p < 0.01) or hyper-pigmented (p < 0.001). This difference was similar despite scar pigmentation phenotype (p = 0.8687). Autografted hyper-pigmented scars had higher rete ridge ratio compared to non-autografted hyper-pigmented HTS (p < 0.0001). Burn hypertrophihc scars have fewer rete ridges than normal skin. This finding may explain the decreased epidermal adherence to underlying dermis associated with hypertrophic scars. Though, contrary to our hypothesis, no direct link between the extent of dyschromia and rete ridge quantity was observed, the differences in normal skin and hypertrophic scar may lead to further understanding of dyschromic scars.

Diverse and variable community structure of picophytoplankton across the Laurentian Great Lakes

AbstractThe Laurentian Great Lakes provide economic support to millions of people, drive biogeochemical cycling, and are an important natural laboratory for characterizing the fundamental components of aquatic ecosystems. Small phytoplankton are important contributors to the food web in much of the Laurentian Great Lakes. Here, for the first time, we reveal and quantify eight phenotypically distinct picophytoplankton populations across the Lakes using a multilaser flow cytometry approach, which distinguishes cells based on their pigment phenotype. The distributions and diversity of picophytoplankton flow populations varied across lakes and depths, with Lake Erie standing out with the highest diversity. By sequencing sorted cells, we identified several distinct lineages of Synechococcales spanning Subclusters 5.2 and 5.3. Distinct genotypic clusters mapped to phenotypically similar flow populations, suggesting that there may not be a clear one‐to‐one mapping between genotypes and phenotypes. This suggests genome‐level differentiation between lakes but some degree of phenotypic convergence in pigment characteristics. Our results demonstrate that ecological selection for locally adapted populations may outpace homogenization by physical transport in this interconnected system. Given the reliance of the Lakes on in situ primary production as a source for organic carbon, this work sets the foundation to test how the community structure of small primary producers corresponds to biogeochemical and food web functions of the Great Lakes and other freshwater systems.

Simultaneous CAS9 editing of cpSRP43, LHCA6, and LHCA7 in Picochlorum celeri lowers chlorophyll levels and improves biomass productivity.

High cellular pigment levels in dense microalgal cultures contribute to excess light absorption. To improve photosynthetic yields in the marine microalga Picochlorum celeri, CAS9 gene editing was used to target the molecular chaperone cpSRP43. Depigmented strains (>50% lower chlorophyll) were generated, with proteomics showing attenuated levels of most light harvesting complex (LHC) proteins. Gene editing generated two types of cpSRP43 transformants with distinct lower pigment phenotypes: (i) a transformant (Δsrp43) with both cpSRP43 diploid alleles modified to encode non-functional polypeptides and (ii) a transformant (STR30309) with a 3nt in-frame insertion in one allele at the CAS9 cut site (non-functional second allele), leading to expression of a modified cpSRP43. STR30309 has more chlorophyll than Δsrp43 but substantially less than wild type. To further decrease light absorption by photosystem I in STR30309, CAS9 editing was used to stack in disruptions of both LHCA6 and LHCA7 to generate STR30843, which has higher (5-24%) productivities relative to wild type in solar-simulating bioreactors. Maximal productivities required frequent partial harvests throughout the day. For STR30843, exemplary diel bioreactor yields of ~50 g m-2day-1 were attained. Our results demonstrate diel productivity gains in P.celeri by lowering pigment levels.

Bacterial-fungal interactions promote parallel evolution of global transcriptional regulators in a widespread Staphylococcus species.

Experimental studies of microbial evolution have largely focused on monocultures of model organisms, but most microbes live in communities where interactions with other species may impact rates and modes of evolution. Using the cheese rind model microbial community, we determined how species interactions shape the evolution of the widespread food- and animal-associated bacterium Staphylococcus xylosus. We evolved S. xylosus for 450 generations alone or in co-culture with one of three microbes: the yeast Debaryomyces hansenii, the bacterium Brevibacterium aurantiacum, and the mold Penicillium solitum. We used the frequency of colony morphology mutants (pigment and colony texture phenotypes) and whole-genome sequencing of isolates to quantify phenotypic and genomic evolution. The yeast D. hansenii strongly promoted diversification of S. xylosus. By the end of the experiment, all populations co-cultured with the yeast were dominated by pigment and colony morphology mutant phenotypes. Populations of S. xylosus grown alone, with B. aurantiacum, or with P. solitum did not evolve novel phenotypic diversity. Whole-genome sequencing of individual mutant isolates across all four treatments identified numerous unique mutations in the operons for the SigB, Agr, and WalRK global regulators, but only in the D. hansenii treatment. Phenotyping and RNA-seq experiments highlighted altered pigment and biofilm production, spreading, stress tolerance, and metabolism ofS. xylosus mutants.Fitness experiments revealed antagonistic pleiotropy, where beneficial mutations that evolved in the presence of the yeast had strong negative fitness effects in other biotic environments. This work demonstrates that bacterial-fungal interactions can have long-term evolutionary consequences within multispecies microbiomes by facilitating the evolution of strain diversity.

Genetic variants in melanogenesis proteins TYRP1 and TYR are associated with the golden rhesus macaque phenotype.

Nonhuman primates (NHPs) are vital translational research models due to their high genetic, physiological, and anatomical homology with humans. The "golden" rhesus macaque (Macaca mulatta) phenotype is a naturally occurring, inherited trait with a visually distinct pigmentation pattern resulting in light blonde colored fur. Retinal imaging also reveals consistent hypopigmentation and occasional foveal hypoplasia. Here, we describe the use of genome-wide association in 2 distinct NHP populations to identify candidate variants in genes linked to the golden phenotype. Two missense variants were identified in the Tyrosinase-related protein 1 gene (Asp343Gly and Leu415Pro) that segregate with the phenotype. An additional and distinct association was also found with a Tyrosinase variant (His256Gln), indicating the light-colored fur phenotype can result from multiple genetic mechanisms. The implicated genes are related through their contribution to the melanogenesis pathway. Variants in these 2 genes are known to cause pigmentation phenotypes in other species and to be associated with oculocutaneous albinism in humans. The novel associations presented in this study will permit further investigations into the role these proteins and variants play in the melanogenesis pathway and model the effects of genetic hypopigmentation and altered melanogenesis in a naturally occurring nonhuman primate model.

Association of functional, inherited vitamin D-binding protein variants with melanoma-specific death.

It is unclear whether genetic variants affecting vitamin D metabolism are associated with melanoma prognosis. Two functional missense variants in the vitamin D-binding protein gene (GC), rs7041 and rs4588, determine 3 common haplotypes, Gc1s, Gc1f, and Gc2, of which Gc1f may be associated with decreased all-cause death among melanoma patients based on results of a prior study, but the association of Gc1f with melanoma-specific death is unclear. We investigated the association of the Gc1s, Gc1f, and Gc2 haplotypes with melanoma-specific and all-cause death among 4490 individuals with incident, invasive primary melanoma in 2 population-based studies using multivariable Cox-proportional hazards regression. In the pooled analysis of both datasets, the patients with the Gc1f haplotype had a 37% lower risk of melanoma-specific death than the patients without Gc1f (hazard ratio [HR] = 0.63, 95% confidence interval [CI] = 0.47 to 0.83, P = .001), with adjustments for age, sex, study center, first- or higher-order primary melanoma, tumor site, pigmentary phenotypes, and Breslow thickness. Associations were similar in both studies. In pooled analyses stratified by Breslow thickness, the corresponding melanoma-specific death HRs for those patients with the Gc1f haplotype compared with those without Gc1f were 0.89 (95% CI = 0.63 to 1.27) among participants with tumor Breslow thickness equal to or less than 2.0 mm and 0.40 (95% CI = 0.25 to 0.63) among participants with tumor Breslow thickness greater than 2.0 mm (Pinteraction = .003). Our findings suggest that individuals with the GC haplotype Gc1f may have a lower risk of dying from melanoma-specifically from thicker, higher-risk melanoma-than individuals without this Gc1f haplotype.

Evolutionary walks through flower colour space driven by gene expression in Petunia and allies (Petunieae).

The structure and function of biochemical and developmental pathways determine the range of accessible phenotypes, which are the substrate for evolutionary change. Accordingly, we expect that observed phenotypic variation across species is strongly influenced by pathway structure, with different phenotypes arising due to changes in activity along pathway branches. Here, we use flower colour as a model to investigate how the structure of pigment pathways shapes the evolution of phenotypic diversity. We focus on the phenotypically diverse Petunieae clade in the nightshade family, which contains ca 180 species of Petunia and related genera, as a model to understand how flavonoid pathway gene expression maps onto pigment production. We use multivariate comparative methods to estimate co-expression relationships between pathway enzymes and transcriptional regulators, and then assess how expression of these genes relates to the major axes of variation in floral pigmentation. Our results indicate that coordinated shifts in gene expression predict transitions in both total anthocyanin levels and pigment type, which, in turn, incur trade-offs with the production of UV-absorbing flavonol compounds. These findings demonstrate that the intrinsic structure of the flavonoid pathway and its regulatory architecture underlies the accessibility of pigment phenotypes and shapes evolutionary outcomes for floral pigment production.

Leukocyte Tyrosine Kinase (Ltk) Is the Mendelian Determinant of the Axolotl Melanoid Color Variant.

The great diversity of color patterns observed among amphibians is largely explained by the differentiation of relatively few pigment cell types during development. Mexican axolotls present a variety of color phenotypes that span the continuum from leucistic to highly melanistic. The melanoid axolotl is a Mendelian variant characterized by large numbers of melanophores, proportionally fewer xanthophores, and no iridophores. Early studies of melanoid were influential in developing the single-origin hypothesis of pigment cell development, wherein it has been proposed that all three pigment cell types derive from a common progenitor cell, with pigment metabolites playing potential roles in directing the development of organelles that define different pigment cell types. Specifically, these studies identified xanthine dehydrogenase (XDH) activity as a mechanism for the permissive differentiation of melanophores at the expense of xanthophores and iridophores. We used bulked segregant RNA-Seq to screen the axolotl genome for melanoid candidate genes and identify the associated locus. Dissimilar frequencies of single-nucleotide polymorphisms were identified between pooled RNA samples of wild-type and melanoid siblings for a region on chromosome 14q. This region contains gephyrin (Gphn), an enzyme that catalyzes the synthesis of the molybdenum cofactor that is required for XDH activity, and leukocyte tyrosine kinase (Ltk), a cell surface signaling receptor that is required for iridophore differentiation in zebrafish. Wild-type Ltk crispants present similar pigment phenotypes to melanoid, strongly implicating Ltk as the melanoid locus. In concert with recent findings in zebrafish, our results support the idea of direct fate specification of pigment cells and, more generally, the single-origin hypothesis of pigment cell development.

Analysis of the genetic loci of pigment pattern evolution in vertebrates.

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.