Unraveling the Molecular Mechanism of Bider Marking Formation in Dun Mongolian Horses Through Transcriptome Sequencing.

This study sequenced skin from dark and light regions of the Bider shoulder stripes, dorsal midline, and croup of Dun Mongolian horses to find genes and pathways that differ between pigmented areas. It implicates coordinated changes in pigment-related signaling (cAMP, calcium, cGMP) and tyrosine metabolism as key drivers of the stripe pattern.

What the researchers asked and why it matters

• Question: Which genes and molecular pathways underlie the formation of the Bider shoulder stripes in Dun Mongolian horses?
• Significance: Bider markings are ancestral “primitive” traits with genetic and evolutionary value. Understanding their molecular basis informs equine coat-color biology, pattern formation, and can guide conservation and selective breeding of traditional phenotypes.

Study design and methods in brief

• Subjects and tissues: Three Dun Mongolian horses; skin sampled from shoulder dark-marked areas, shoulder light-marked areas, dorsal midline, and croup, enabling within-animal regional comparisons.
• Transcriptomics: High-depth RNA sequencing produced ~893 million clean reads in total with >96% overall alignment to the horse reference, indicating high data quality.
• Differential expression: Pairwise comparisons across regions identified 140–775 differentially expressed genes (DEGs) per contrast, reflecting region-specific transcriptional programs.
• Functional analyses: DEGs underwent Gene Ontology (GO) annotation and KEGG pathway enrichment to interpret biological roles.
• Validation: Selected genes were tested by qRT-PCR, confirming the RNA-seq expression trends and supporting data reliability.

Main findings

• GO enrichment highlighted processes central to pigmentation and skin biology, including:
  • Pigment metabolism and melanogenesis
  • Skin and hair follicle development
  • Signal transduction, notably calcium and cGMP signaling
  • Immune regulation processes that can modulate melanocyte behavior
• KEGG pathways significantly enriched across multiple region comparisons included:
  • Calcium signaling pathway
  • Tyrosine metabolism (the biochemical backbone of melanin synthesis)
  • cAMP signaling pathway (a canonical driver of melanogenesis)
  • Melanoma pathway (a KEGG term that encompasses normal melanocyte regulatory networks)
• Convergence on these pigment-related pathways across different tissue contrasts suggests shared, core mechanisms shaping the darker versus lighter regions of the Bider pattern.
• qRT-PCR corroborated RNA-seq-derived expression differences, strengthening confidence in the results.

Biological interpretation: how these pathways could create a stripe

• cAMP signaling:
  • In melanocytes, activation of MC1R increases cAMP, which via PKA and CREB elevates MITF activity and expression of melanogenic enzymes (e.g., TYR, TYRP1, DCT), boosting eumelanin production.
  • Higher cAMP pathway activity in dark-marked regions would plausibly increase melanin synthesis and dendricity, deepening pigment intensity.
• Calcium signaling:
  • Intracellular Ca2+ influences MITF phosphorylation (via CaMK pathways), melanosome transport (through myosin Va/Rab27a/melanophilin complexes), and melanocyte morphology.
  • Spatial differences in Ca2+ signaling could sharpen boundaries by altering pigment production and organelle trafficking at the stripe edges.
• cGMP signaling:
  • NO–sGC–cGMP–PKG pathways modulate melanocyte differentiation and dendricity and interact with Ca2+ and cAMP signaling.
  • Regional tuning of cGMP could reinforce or counterbalance cAMP/Ca2+ effects to set light versus dark domains.
• Tyrosine metabolism:
  • Tyrosine is the substrate for melanin; upregulation of this pathway typically reflects increased activity of enzymes like tyrosinase and its partners.
  • Enrichment suggests a biochemical basis for darker pigment in the marked areas.
• “Melanoma” KEGG pathway:
  • Although named for disease, it aggregates core melanocyte signaling (MAPK, PI3K, Wnt) used in normal pigment regulation; enrichment indicates activation of these modules in patterned skin.
• Skin and hair follicle development:
  • Hair follicle stage and structure strongly influence pigment deposition; region-specific follicle programs could help establish or maintain stripes.
• Immune regulation:
  • Cytokines and paracrine mediators (e.g., endothelins, interferons, TNF) can up- or down-regulate melanogenesis; local immune tone might fine-tune pigment intensity across the pattern.

Context with prior knowledge on Dun/primitive markings

• Dun horses characteristically show dorsal and shoulder striping (“primitive markings”). Prior genetic work in dun coloration has implicated regionally restricted regulatory programs affecting pigment deposition in hair follicles.
• This study complements that body of work by identifying, at the transcriptome level in patterned skin, coordinated changes across multiple pigment-control pathways rather than a single causative gene.

Strengths and limitations

• Strengths:
  • Within-animal, region-matched sampling reduces inter-individual confounders.
  • High sequencing depth and excellent alignment rates support data quality.
  • Orthogonal qRT-PCR validation increases confidence in differential expression calls.
• Limitations:
  • Small sample size (n=3 horses) limits power and generalizability.
  • Bulk skin RNA-seq mixes cell types (melanocytes, keratinocytes, fibroblasts, immune cells), obscuring cell-specific signals.
  • Descriptive (correlative) results; causal roles of pathways/genes remain to be tested functionally.
  • Potential confounders such as hair cycle stage or sun exposure were not dissected.
  • Genotypes at known coat-color loci (e.g., MC1R, ASIP, TBX3, PMEL) were not reported, which could influence baseline expression.

Implications and applications

• Identifies candidate pathways likely orchestrating region-specific melanocyte activity that produces the Bider stripe.
• Provides a molecular foundation for future breeding, conservation of ancestral markings, and comparative studies of stripe formation across species.
• Suggests biomarkers for pigment activity that could be monitored in equine dermatology and pigmentation research.

Recommended next steps

• Single-cell or spatial transcriptomics to map melanocyte, keratinocyte, and immune-cell programs across stripe boundaries.
• Histology and immunostaining (e.g., MITF, TYR, DCT, PMEL) and direct quantification of eumelanin/pheomelanin to link expression with cellular and biochemical phenotypes.
• In vitro functional tests in equine melanocytes modulating cAMP, Ca2+, and cGMP pathways to establish causality.
• Genotyping of known coat-color loci and integration with expression to parse cis-regulatory versus signaling-driven effects.
• Longitudinal sampling across hair-cycle stages and seasons to control for follicular dynamics.

Key takeaways

• Bider stripe formation appears to result from coordinated, region-specific regulation of melanocyte pathways, especially cAMP, calcium, cGMP signaling, and tyrosine metabolism.
• High-quality transcriptomic data with qRT-PCR support provide strong evidence, but mechanistic and cell-type–resolved studies are needed to confirm causal drivers.