Abstract
Color molts from summer brown to winter white coats have evolved in several species to maintain camouflage year-round in environments with seasonal snow. Despite the eco-evolutionary relevance of this key phenological adaptation, its molecular regulation has only recently begun to be addressed. Here, we analyze skin transcription changes during the autumn molt of the mountain hare (Lepus timidus) and integrate the results with an established model of gene regulation across the spring molt of the closely related snowshoe hare (L. americanus). We quantified differences in gene expression among three stages of molt progression—“brown” (early molt), “intermediate,” and “white” (late molt). We found 632 differentially expressed genes, with a major pulse of expression early in the molt, followed by a milder one in late molt. The functional makeup of differentially expressed genes anchored the sampled molt stages to the developmental timeline of the hair growth cycle, associating anagen to early molt and the transition to catagen to late molt. The progression of color change was characterized by differential expression of genes involved in pigmentation, circadian, and behavioral regulation. We found significant overlap between differentially expressed genes across the seasonal molts of mountain and snowshoe hares, particularly at molt onset, suggesting conservatism of gene regulation across species and seasons. However, some discrepancies suggest seasonal differences in melanocyte differentiation and the integration of nutritional cues. Our established regulatory model of seasonal coat color molt provides an important mechanistic context to study the functional architecture and evolution of this crucial seasonal adaptation.
| Original language | English |
|---|---|
| Pages (from-to) | 1180-1192 |
| Number of pages | 13 |
| Journal | Ecology and Evolution |
| Volume | 10 |
| Issue number | 3 |
| DOIs | |
| State | Published - Feb 1 2020 |
Funding
Fundação para a Ciência e a Tecnologia (Portuguese National Funds), Grant/Award Number: PTDC/BIAEVF/1624/2014; National Science Foundation, Grant/Award Number: OIA‐1736249, DEB‐1907022 and DEB‐1743871; Fundação para a Ciência e a Tecnologia (POPH‐QREN funds from the European Social Fund and Portuguese MCTES), Grant/Award Number: PD/BD/108131/2015 and IF/00033/2014/CP1256/CT0005; NORTE2020 and PORTUGAL2020 through the European Regional Development Fund, Grant/Award Number: NORTE‐01‐0145‐FEDER‐000007; COMPETE2020 and PORTUGAL2020 through the European Regional Development Fund, Grant/Award Number: POCI‐01‐0145‐FEDER‐022184; Fundação Luso‐Americana para o Desenvolvimento; European Union's Seventh Framework Program, Grant/Award Number: 286431; M.J. Murdock Charitable Trust, Grant/Award: University of Montana Genomic Core; NIH Instrumentation, Grant/Award Number: S10RR029668 and S10RR027303; LabEx CeMEB MBB (ANR “Investissements d'avenir” program), Grant/Award Number: ANR‐10‐LABX‐04‐01. Financial support was obtained from Fundação para a Ciência e a Tecnologia (FCT), project grant "CHANGE" (reference PTDC/BIA‐EVF/1624/2014, funded by Portuguese National Funds through FCT) to JM‐F, and from the National Science Foundation (grant OIA‐1736249) to JMG. MSF was supported by POPH‐QREN funds from the European Social Fund (ESF) and Portuguese MCTES/FCT (PD/BD/108131/2015 PhD grant in the scope of BIODIV PhD program at Faculty of Sciences, University of Porto). JM‐F was supported by FCT Investigator grant IF/00033/2014/CP1256/CT0005 (POPH‐QREN funds from ESF and Portuguese MCTES/FCT). Additional support was provided by the National Science Foundation (DEB‐1907022) to JMG and LSM, by the National Science Foundation (DEB‐1743871) to LSM, by NORTE2020, PORTUGAL2020, through the European Regional Development Fund—ERDF (NORTE‐01‐0145‐FEDER‐000007), by COMPETE2020, PORTUGAL2020, and ERDF (POCI‐01‐0145‐FEDER‐022184), and by Portugal—United States of America Research Networks Program funds from Fundação Luso‐Americana para o Desenvolvimento (FLAD) to PCA and MSF. Instrumentation, laboratory, and computational support was provided by CIBIO NEW‐GEN sequencing platform, supported by European Union's Seventh Framework Program for research, technological development and demonstration under grant agreement no. 286431, the University of Montana Genomics Core, supported by a grant from the M.J. Murdock Charitable Trust, and the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley, supported by NIH S10 Instrumentation Grants S10RR029668 and S10RR027303. We thank Pierre Boursot for providing access to the Montpellier Bioinformatics Biodiversity platform supported by the LabEx CeMEB, an ANR "Investissements d'avenir" program (ANR‐10‐LABX‐04‐01) for the de novo assembly of the transcriptome, and João Pedro Marques for bioinformatic support. We thank Fabiana Neves, from CIBIO‐InBIO's IMED group for primers to amplify ACTB gene. We are grateful to René Gadient and local hunters in Grisons for their tremendous effort in collecting samples under Alpine conditions. Funding information Funda??o para a Ci?ncia e a Tecnologia (Portuguese National Funds), Grant/Award Number: PTDC/BIAEVF/1624/2014; National Science Foundation, Grant/Award Number: OIA-1736249, DEB-1907022 and DEB-1743871; Funda??o para a Ci?ncia e a Tecnologia (POPH-QREN funds from the European Social Fund and Portuguese MCTES), Grant/Award Number: PD/BD/108131/2015 and IF/00033/2014/CP1256/CT0005; NORTE2020 and PORTUGAL2020 through the European Regional Development Fund, Grant/Award Number: NORTE-01-0145-FEDER-000007; COMPETE2020 and PORTUGAL2020 through the European Regional Development Fund, Grant/Award Number: POCI-01-0145-FEDER-022184; Funda??o Luso-Americana para o Desenvolvimento; European Union's Seventh Framework Program, Grant/Award Number: 286431; M.J. Murdock Charitable Trust, Grant/Award: University of Montana Genomic Core; NIH Instrumentation, Grant/Award Number: S10RR029668 and S10RR027303; LabEx CeMEB MBB (ANR ?Investissements d'avenir? program), Grant/Award Number: ANR-10-LABX-04-01. Financial support was obtained from Funda??o para a Ci?ncia e a Tecnologia (FCT), project grant "CHANGE" (reference PTDC/BIA-EVF/1624/2014, funded by Portuguese National Funds through FCT) to JM-F, and from the National Science Foundation (grant OIA-1736249) to JMG. MSF was supported by POPH-QREN funds from the European Social Fund (ESF) and Portuguese MCTES/FCT (PD/BD/108131/2015 PhD grant in the scope of BIODIV PhD program at Faculty of Sciences, University of Porto). JM-F was supported by FCT Investigator grant IF/00033/2014/CP1256/CT0005 (POPH-QREN funds from ESF and Portuguese MCTES/FCT). Additional support was provided by the National Science Foundation (DEB-1907022) to JMG and LSM, by the National Science Foundation (DEB-1743871) to LSM, by NORTE2020, PORTUGAL2020, through the European Regional Development Fund?ERDF (NORTE-01-0145-FEDER-000007), by COMPETE2020, PORTUGAL2020, and ERDF (POCI-01-0145-FEDER-022184), and by Portugal?United States of America Research Networks Program funds from Funda??o Luso-Americana para o Desenvolvimento (FLAD) to PCA and MSF. Instrumentation, laboratory, and computational support was provided by CIBIO NEW-GEN sequencing platform, supported by European Union's Seventh Framework Program for research, technological development and demonstration under grant agreement no. 286431, the University of Montana Genomics Core, supported by a grant from the M.J. Murdock Charitable Trust, and the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley, supported by NIH S10 Instrumentation Grants S10RR029668 and S10RR027303. We thank Pierre Boursot for providing access to the Montpellier Bioinformatics Biodiversity platform supported by the LabEx CeMEB, an ANR "Investissements d'avenir" program (ANR-10-LABX-04-01) for the de novo assembly of the transcriptome, and Jo?o Pedro Marques for bioinformatic support. We thank Fabiana Neves, from CIBIO-InBIO's IMED group for primers to amplify ACTB gene. We are grateful to Ren? Gadient and local hunters in Grisons for their tremendous effort in collecting samples under Alpine conditions.
| Funder number |
|---|
| POCI‐01‐0145‐FEDER‐022184 |
| PORTUGAL2020 |
| PD/BD/108131/2015 PhD |
| DEB‐1743871, OIA‐1736249, DEB‐1907022 |
| S10RR027303, ANR‐10‐LABX‐04‐01, S10RR029668 |
| PTDC/BIA‐EVF/1624/2014 |
| 286431 |
| PD/BD/108131/2015, IF/00033/2014/CP1256/CT0005 |
| NORTE‐01‐0145‐FEDER‐000007 |
Keywords
- developmental timeline
- gene expression
- molt cycle
- seasonal coat color change
- transcriptomics