# Skin Colour in Salamanders Is Modulated by Both Epitranscriptomic Methylation and Gene Expression

**Authors:** Nicholas Strowbridge, David R. Vieites, Michael G. Ritchie, Kathryn R. Elmer

PMC · DOI: 10.1111/mec.70302 · 2026-03-20

## TL;DR

This study shows that both RNA methylation and gene expression influence skin color in salamanders, revealing new molecular mechanisms behind their color diversity.

## Contribution

The study identifies RNA methylation and gene expression as key factors in salamander coloration, an underexplored area in amphibians.

## Key findings

- 129 genes were differentially expressed and 281 differentially methylated across color comparisons.
- Genes like MLANA, PMEL, TYR, and TYRP1, linked to melanin production, were identified as key players.
- A positive correlation and overlap were found between differentially methylated and expressed transcripts.

## Abstract

Animal colouration is a key trait in organismal biology, being involved in natural and sexual selection, competition, and communication. Amphibians use their highly diverse colouration in many ecological interactions, but the molecular bases of their colour variation are less well understood than in other vertebrate systems. While the genetic, structural, and cellular bases of pigmentation are increasingly understood in a range of models, potential epigenetic or epitranscriptomic effects are almost completely unexplored. The fire salamander (
Salamandra salamandra
) has striking colour patterns and polymorphisms, but the extremely large genome size of salamanders makes traditional genetic analyses infeasible. To discover loci and molecular mechanisms underlying colour differences in salamanders, we used long‐read direct RNA sequencing to test the roles of RNA methylation, gene expression and their relationship on intra‐ and inter‐individual colour variation in black, yellow and brown skinned salamanders. We found 129 differentially expressed and 281 differentially methylated genes across all pairwise comparisons. Many of the genes involved are related to pigmentation, with several directly associated with melanin production, such as Melan‐A (MLANA), Premelanosome protein (PMEL), Tyrosinase (TYR) and Tyrosinase‐related protein 1 (TYRP1). We found both a positive overall correlation and a significant overlap in transcripts that are differentially methylated and expressed. These findings suggest multiple molecular mechanisms, including gene expression and RNA methylation, contribute to amphibian colour diversity. RNA modifications as a promising area for understanding morphological variation in non‐model animals and their impact on their ecology.

## Linked entities

- **Genes:** MLANA (melan-A) [NCBI Gene 2315], PMEL (premelanosome protein) [NCBI Gene 6490], TYR (tyrosinase) [NCBI Gene 7299], TYRP1 (tyrosinase related protein 1) [NCBI Gene 7306]
- **Species:** Salamandra salamandra (taxon 57571)

## Full-text entities

- **Genes:** Tyrp1 (tyrosinase-related protein 1) [NCBI Gene 22178] {aka Oca3, TRP-1, TRP1, Tyrp, b, brown}, Tyr (tyrosinase) [NCBI Gene 22173] {aka Oca1, albino, c, skc35}, MLANA (melan-A) [NCBI Gene 2315] {aka MART-1, MART1}, Txn1 (thioredoxin 1) [NCBI Gene 22166] {aka ADF, Trx1, Txn}, Fam89a (family with sequence similarity 89, member A) [NCBI Gene 69627] {aka 2310031A18Rik}, Mlana (melan-A) [NCBI Gene 77836] {aka A930034P04Rik, Mart1}, Atp6v1h (ATPase, H+ transporting, lysosomal V1 subunit H) [NCBI Gene 108664] {aka 0710001F19Rik, CGI-11, SFD, SFDalpha, SFDbeta, VMA13}, Mapk1ip1l (mitogen-activated protein kinase 1 interacting protein 1-like) [NCBI Gene 218975] {aka 9130202B12, C130032J12Rik}, Phb2 (prohibitin 2) [NCBI Gene 12034] {aka BAP, Bap37, Bcap37, REA}, Krt12 (keratin 12) [NCBI Gene 268482] {aka K12, Krt-12, Krt1-12}, TYRP1 [NCBI Gene 100034666], Myo1b (myosin IB) [NCBI Gene 17912] {aka MIH-L, MMIa}, Qdpr (quinoid dihydropteridine reductase) [NCBI Gene 110391] {aka 2610008L04Rik, D5Ertd371e, Dhpr, PKU2}, Cxcl12 (C-X-C motif chemokine ligand 12) [NCBI Gene 20315] {aka Pbsf, Scyb12, Sdf1, Tlsf, Tpar1}, Clu (clusterin) [NCBI Gene 12759] {aka ApoJ, Cli, D14Ucla3, SP-40, Sgp-2, Sgp2}, Sfrp2 (secreted frizzled-related protein 2) [NCBI Gene 20319] {aka Sdf5}, mc1r [NCBI Gene 102099064], TYRP1 (tyrosinase related protein 1) [NCBI Gene 7306] {aka CAS2, CATB, GP75, OCA3, TRP, TRP1}, Prdx3 (peroxiredoxin 3) [NCBI Gene 11757] {aka Aop1, D0Tohi1, Ef2l, Mer5, Prx3, SP22}, Atox1 (antioxidant 1 copper chaperone) [NCBI Gene 11927] {aka Atx1}, Krt13 (keratin 13) [NCBI Gene 16663] {aka CK13, K13, Krt-1.13, Krt1-13}, Lgals3 (lectin, galactose binding, soluble 3) [NCBI Gene 16854] {aka GBP, L-34, Mac-2, gal3}, Notch1 (notch 1) [NCBI Gene 18128] {aka 9930111A19Rik, Mis6, N1, Tan1, lin-12}, TYR (tyrosinase) [NCBI Gene 7299] {aka ATN, CMM8, OCA1, OCA1A, OCAIA, SHEP3}, Fmo3 (flavin containing monooxygenase 3) [NCBI Gene 14262], Gpnmb (glycoprotein (transmembrane) nmb) [NCBI Gene 93695] {aka DC-HIL, Dchil, ipd}, PMEL [NCBI Gene 101112637], PMEL (premelanosome protein) [NCBI Gene 6490] {aka D12S53E, HMB-45, HMB45, ME20, ME20-M, ME20M}, Atp7a (ATPase, copper transporting, alpha polypeptide) [NCBI Gene 11977] {aka MNK}, Dgat2 (diacylglycerol O-acyltransferase 2) [NCBI Gene 67800] {aka 0610010B06Rik, ARAT, DGAT-2}, Pard3 (par-3 family cell polarity regulator) [NCBI Gene 93742] {aka Asip, D8Ertd580e, Par-3, Par3, Pard-3, Pard3a}, TYR [NCBI Gene 780502], Pmel (premelanosome protein) [NCBI Gene 20431] {aka D10H12S53E, D12S53Eh, Pmel17, Si, Silv, gp100}
- **Diseases:** melanism (MESH:D008548), pigmentation (MESH:D010859), cytotoxicity (MESH:D064420), DM (MESH:D012734)
- **Chemicals:** glutathione (MESH:D005978), sodium bicarbonate (MESH:D017693), BH4 (MESH:C003402), Melanin (MESH:D008543), dopaquinone (MESH:C035157), sodium acetate (MESH:D019346), copper (MESH:D003300), tyrosine (MESH:D014443), pheomelanin (MESH:C018362), H2O. (MESH:D014867), DHICA (-), eumelanin (MESH:C041877), dopachrome (MESH:C001123), chloroform (MESH:D002725), ethanol (MESH:D000431), N6-methyladenosine (MESH:C010223), Trizol (MESH:C411644), MS-222 (MESH:C003636), Carotenoids (MESH:D002338)
- **Species:** Salamandra salamandra bernardezi (subspecies) [taxon 297008], Caudata (salamanders, order) [taxon 8293], Mus musculus (house mouse, species) [taxon 10090], Bacteria Latreille et al. 1825 (Bacteria stick insect, genus) [taxon 629395], Danio rerio (leopard danio, species) [taxon 7955], Homo sapiens (human, species) [taxon 9606], Tilapia (genus) [taxon 8126], Ovis aries (domestic sheep, species) [taxon 9940], Salamandra salamandra (European fire salamander, species) [taxon 57571]
- **Mutations:** C-18 C

## Figures

4 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13003378/full.md

---
Source: https://tomesphere.com/paper/PMC13003378