Abstract
Every year as spring approaches and day length increases, many birds begin to reproduce, an annual expression of seasonal phenology that requires physiological preparation. In species distributed over a broad geographic range, populations that breed at higher latitudes are often migratory and delay reproduction until later in the year as compared to those breeding at lower latitudes. Dark-eyed Juncos serve as an excellent model for understanding the timing mechanisms regulating population-level variation in seasonal reproductive responses. We compared two seasonally sympatric dark-eyed junco populations in early spring. One migrates (Junco hyemalis hyemalis) and breeds in Alaska and Canada, while the other remains resident (Junco hyemalis carolinensis) and breeds in the Appalachian Mountains of Virginia USA. These populations exhibit different photoperiodic responses to the same environment with respect to activation of the HPG axis, leading to earlier gonadal recrudescence in the resident population. We caught co-wintering sympatric male migrant (n = 6) and resident (n = 7) juncos from the field in March and collected the hypothalamic tissues. We also collected blood samples to determine circulating testosterone and a wing feather to determine stable isotope ratios (δ2H) as estimate of breeding latitude. We found three differentially expressed genes, among which gonadotropin releasing hormone 1 (GnRH1) showed significantly higher expression in early breeding residents as compared to migrant juncos. The δ2H showed a positive linear correlation with testosterone levels and GnRH1 mRNA, providing strong evidence for latitudinal variation in breeding phenology. This study provides insight into the underlying neuroendocrine response giving rise to a population-level difference in the timing of reproduction observed in a seasonally sympatric (co-wintering) population of resident and migrant juncos.
Introduction
Seasonal animals have evolved different timing strategies to match their recurring physiology of migration and reproduction to the environment (). Annual changes in day length (photoperiod) are a primary entraining cue for endogenous clocks to synchronize with the environment. Avian populations that live in sympatry during the non-breeding state from winter to early spring may differ in whether or not they migrate. If they do, then they may breed at different locations in spring and summer, thus becoming allopatric (; ; ). One such species is the dark-eyed junco, which is distributed over a broad geographic range. In the junco, closely related populations winter in sympatry but differ in when they breed and whether they migrate. This pattern of sympatry in winter and allochrony in breeding makes the dark-eyed junco a powerful model to investigate the neuroendocrine basis of differences in the timing of reproduction because they pass through an overwintering period together and start to differ in their behavioral and physiological responses to similar environmental cues as spring approaches.
Use of hydrogen isotope signatures in the movement ecology of birds is a common tool in practice to predict breeding latitude (; ; ). Birds incorporate isotope signatures through their diet or nearby water bodies (; ). The local isotopic signature is incorporated into the growing feather during the post breeding molt. Because feathers are epidermal growths that are metabolically inactive, they later reflect the local isotopic signature from the geographical location where they were grown (). Stable isotopes such as hydrogen vary geographically and are used to estimate the breeding location and location of feather growth (; ; ; ). In the case of the junco, migrant and resident populations can be distinguished during winter by the lighter hydrogen stable isotope values in migrants as compared to the heavier values found in residents.
Songbirds undergo seasonal recrudescence in their reproductive processes, primarily regulated by hypothalamic GnRH neurons (; ). In the temperate zone, birds initiate their gonadal development in response to increasing day lengths during early spring. However, the time of reproductive response varies depending on breeding latitude. The seasonal reproductive response is induced primarily through the interaction of increasing day length with the hypothalamic-pituitary-gonadal (HPG) axis (; ). After the winter solstice, the photoperiod starts to increase, which stimulates the HPG on reaching a population-specific threshold level of photoperiod (; ). The annual change in day length varies with latitude, such that birds breeding at higher latitudes often require longer days to initiate gonadal recrudescence than birds breeding at lower latitudes (; ; ). The birds breeding at lower latitudes recrudesce gonads earlier, stay in the stimulatory phase for a longer period, and terminate breeding later than high latitude birds that have a short breeding period ().
In studies done on populations that overwinter in sympatry and breed at different locations, the higher latitude breeding birds have to accommodate an additional life-history event of spring migration that delays their reproduction (; ; ). This process of spring migratory preparedness involves seasonal pre-migratory fattening, increase in food intake, and an increase in metabolic activity in flight muscles and liver in phenology and behavior (; ; ). There is also a striking shift in the locomotor activity of migratory birds, switching from only day activity to nighttime activity (; ; ; ; ). These changes in behavior and physiology associated with seasonal induction of reproduction and migration are potentially an outcome of some key RNA molecules changing their expression in the hypothalamus in response to changes in spring day length and other timing cues (; ; ). The hypothalamus along with the pituitary, is a key neural center that synthesizes and releases neuropeptides and hormones critical to the regulation of seasonal reproductive and migratory processes. A key step in the induction of the gonad recrudescence is initiated by release of gonadotropin releasing hormone 1 (GnRH1) from the hypothalamus followed by its binding to the pituitary gonadotroph to stimulate release of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) (). Investigating molecules involved in gonad growth, production of sex steroids, energy homeostasis and metabolism may help advance our understanding of differences in the reproductive timing of songbirds breeding at different latitudes.
The Dark-eyed junco provides a powerful model system for investigating the neuroendocrine mechanisms that underlie variation in reproductive and migratory timing. In early spring, populations of juncos that differ in whether they migrate are found living together in the same environment (; ). As spring progresses, the locally breeding population undergoes gonadal recrudescence, while the migratory population delays recrudescence to prepare to migrate. In a common garden study on resident and migrant juncos held in captivity, neuroendocrine tissues collected during the reproductive divergence state revealed differential patterns of gene and protein expression between residents and migrants (; ; ). However, a comprehensive account of the transcriptome associated with divergent state physiology in neuroendocrine tissue is needed to better understand the underlying mechanisms regulating early spring-induced phenologies.
The objective of this study was to discover key transcripts in the hypothalamus tissue of free-living juncos on the verge of divergence in early spring. Prior studies () conducted on captives provided insight into the mechanisms through which population-level seasonal divergence is set in a controlled indoor condition, including differentially expressed genes in muscle, blood, and hypothalamus. The present study asks how the populations differ in their hypothalamic gene expression when sampled directly from their natural environment during the time of seasonal divergence.
The following hypothesis to understand how neuroendocrine response will determine the timing of reproduction in junco populations on the verge of divergence: 1) We predicted to find differences in the plasma testosterone levels as a physiological measure of gonadal recrudescence. 2) We predicted to see the differences in the key hypothalamic transcripts associated with the latitude of breeding. 3) We predicted to see differential expression of hypothalamic transcripts associated with migratory preparedness. To test these hypotheses, we sampled sympatric resident/migrant birds that breed at different latitudes and flock together during early spring with a pre-notion of individual variation in response to constantly changing natural day length and climatic conditions.
Materials and methods
The study system
The dark-eyed junco subspecies live in sympatry and flock together at Mountain Lake Biological Station (MLBS), Virginia, from early fall to spring and start to diverge as soon as the photoperiod starts increasing (Figures 1a,b). In the wintering ground, the resident and migrant juncos are identified by their beak color. The migrants are pink-billed, and residents have blue-billed resident birds that breed at MLBS (; ). The migrant juncos start to prepare for departure to the North as day length increases. Whereas resident juncos start showing gonadal recrudescence (Figure 1).
FIGURE 1
Body mass and Fat score
After taking the bird out of the mist net, we visually estimated subcutaneous fat score (0–5) deposition in the furcular and abdominal regions, as follows: no fat present (0), trace amounts of fat in the furcular region (1), trace amounts of fat in the abdominal region (2), half full in both furcular and abdomen area (3), full in both the areas (4), bulging in both the areas (5). Body mass was measured using a Pesola Spring scale (0–50 g), tared to zero with a small piece of sock used to carry the bird for body mass measurement.
Capture and tissue collection
We captured six migrants and seven resident male dark-eyed juncos using mist nets from state road 714, Giles County Virginia at Mountain Lake Biological Station (MLBS; 37.37 0N, 80.52 0W). Migrants (J. h hyemalis) and residents (J. h carolinensis) were identified using bill coloration, plumage, and wing cord (). Birds were caught from 13 to 20 March 2017 during spring initiation at MLBS (11.72 ± 0.1 h of average day length and −1 to −4 °C average daily temperature). Juncos were euthanized using isoflurane within 5 min from the time they were caught in the mist net. Before euthanizing, 100 μL of blood was collected by puncturing the wing vein for baseline testosterone hormone measurement. After the capture of the birds, the brain and testes were immediately dissected within 5 min, and tissues were flash-frozen on dry ice and stored at −80 °C until further processing of tissues. Scientific collecting permits were issued by the Virginia Department of Game and Inland Fisheries (permit # 052971). All methods were approved under a protocol (# 15-026-17) by the Indiana University Institutional Animal Care and Use Committee.
The hypothalamic area was dissected from the whole brain under the cryostat (Leica Biosystems CM1850, Buffalo Grove, IL, United States) maintained at −20 °C using a modified protocol adapted from . The frozen brain was sliced on the rostral caudal position until we saw the tractus septomesencephalicus (TrSM) as a landmark for the beginning of the hypothalamus. We sliced 30 μm thick sections and punched directly above the optic chiasma using a sterile 3 mm-diameter biopsy punch (Thermo Fisher Scientific, Integra™ 3332; Cat no. 12-460-406). All the tissue punches starting from TrSM to the end of the infundibular nuclear complex (INc) as a posterior landmark for the hypothalamus were collected in a sterile 1.5 mL eppendorf tube kept in cryostat throughout brain punching (; ; ).
Stable hydrogen isotope analyses
The distal most secondary feather from each individual was collected and stored in the coin envelopes and transported to the Indiana University, Bloomington. All feathers were cleaned using a chloroform: methanol in a 2:1 ratio to remove any external contaminants. A small fraction from the tip of the feather was cut, weighed to approximately 0.5 mg, placed in a 3 × 5 mm silver capsule, and shipped to a laboratory located in the US Geological Survey in Denver, CO, to analyze the hydrogen isotope (δ2H) values. Hydrogen isotope ratios were measured using established methods of mass spectrometry (; ; ; ). The non-exchangeable δ2H values were reported in parts per mil notation (‰) with respect to VSMOW (Vienna Standard Mean Oceanic Water) using Caribou (−157‰) and Kudu (−35.3‰) standards. We used δ2H values to estimate the latitude and photoperiod where each bird spent the breeding, molting season and incorporated local δ2H signatures in the growing feathers (). To estimate the latitude and photoperiod for each corresponding δ2H value, we obtained North America hydrogen isotope ratios for August month precipitation (δ2Hp) from the OPIC 3.2 database (). For each bird, we found all feather δ2H values in the North America compatible to δ2H precipitation values (for detailed methods see ). The δ2H values were used as a continuous variable against all the physiological, hormone measures, and expressed transcripts.
Testosterone ELISA
We measured circulating baseline testosterone levels in blood plasma by puncturing the wing vein and collecting blood into heparinized capillary tubes. The blood samples were stored at 4 °C and plasma was extracted by centrifuging blood samples at 3000 rpm for 15 min at room temperature. The plasma was separated using a Hamilton syringe and stored at −20 °C until assayed to measure testosterone. The plasma samples were thawed on ice, and 20 μL of plasma from each individual was aliquoted in glass tubes. Hormone was extracted using diethyl-ether mixed by agitating and incubation at room temperature for 20 min. After incubation, all tubes were snap frozen and the supernatant was immediately transferred to a fresh tube. This procedure was repeated three times to extract the hormone, and finally the tube was dried using high pressure liquid nitrogen. Assay buffer (250 μL) was added to each tube, and 100 μL of extracted hormone was used in duplicate wells to run the testosterone assay. Plasma testosterone was measured using high sensitivity testosterone ELISA kit (Enzo ADI-900-176) as per manufacturer’s protocol. All samples were placed on one plate, testosterone ELISA with intra-assay coefficient variability of 0.6% ± 0.06 (mean ± SE) and sensitivity of 2.6 pg/mL.
Measurement of gene expression
RNA extraction and sequencing
A total of 13 brains (6 migrant and 7 resident juncos) hypothalamii were punched and homogenized in RLT buffer containing 1% 2-Mercaptoethanol (BME) on ice. Total RNA was extracted with AllPrep DNA/ RNA Universal kit (Qiagen, Cat. # 80224) according to manufacturer protocol. RNA concentration and Integrity (RIN; ) were determined on Agilent 2200 TapeStation system. Total RNA was submitted to the Center for Genomics and Bioinformatics at Indiana University for subsequent RNA quality check, RNA library preparation and Illumina Nexseq 500 (43 nucleotides, paired end reads).
Quality filtering and mapping
The paired-end reads were de-multiplexed using bcl2fastq version v2.20.0.422 and about 28.5 million reads, with a raw read length of 43 nucleotides were assigned to each library. Raw reads were adapter trimmed and quality filtered using Trimmomatic version 0.33 () requiring a minimum base quality score of 20 averaged across a sliding window of 3 bases. Reads shorter than 20 bases post trimming were discarded. Trimmed reads were mapped to the dark-eyed junco draft genome assembly sequence using the RNA-Seq read aligner STAR version 2.5.3a ().
Differential gene expression analysis
The transcriptome was annotated against functionally annotated regions defining the potentially transcribed elements on the dark eye junco draft genome assembly homologous to protein sequences from zebra finch (Taeniopygia guttata, version 3.2.4), a Passeriformes species phylogenetically closer to junco, were identified and annotated using BLASTX (). Counts for read pairs aligning uniquely to each of the annotated transcribed elements were generated using feature Counts version 1.6.3 () and were used to identify the potentially significantly differentially expressed regions between migrant and resident birds at 5% FDR using the Bioconductor R package DESeq2 version 1.12.3 ().
Statistical Analysis
All analysis was done using R (version 3.2.0). Differences in mean hydrogen isotope ratio and T levels were calculated using a non-parametric Student’s t-test (Mann-Whitney U test) between subpopulation means. We used the Pearson correlation method to identify any significant correlation between hydrogen isotope ratios, testosterone levels, and normalized Differentially expressed gene expression values. For statistical significance, alpha was set at 0.05. The Benjamini–Hochberg procedure with a failed discovery rate of 5% was used to correct for multiple tests ().
Results
Differences in seasonal phenotypes on the verge of spring
Every year from fall to early spring, dark-eyed junco subspecies are sympatric and flock together at Mountain Lake Biological Station (MLBS), Virginia (Figures 1a,b). The pink-billed migrants arrive in the beginning of fall and depart to the north as spring approaches. Whereas blue-billed resident birds live year-round and breed at MLBS. In early spring, the migrants start preparing for migration. By the time resident juncos start showing gonadal recrudescence, the migrants are departing for their breeding ground. (Figure 1). The residents breeding and growing feathers at lower latitude had heavier δ2H than migrants (Mann-Whitney U test, p = 0.0012, Figure 2a). Also, there were significantly higher blood testosterone levels in resident as compared to migratory birds (Mann-Whitney U test, p = 0.014, Figure 2b). We did not find any significant difference in the body weight or subcutaneous fat scores of migrant and resident birds.
FIGURE 2
Differentially expressed genes
A total of 11,454 expressed sequenced tags (ESTs) were found after mapping to the dark-eyed junco draft genome assembly homologous to protein sequences from zebra finch (Taeniopygia guttata) as summarized in Supplementary Table S1. Further, we found three DEGs (FDR <0.05; Figure 3; Supplementary Table S1). Among the three DEGs, only one, gene macrophage mannose receptor 1-like (MMR1), was downregulated in resident birds (FDR = 0.0135; Figure 3a). Unlike MMR1, uncharacterized protein C21orf58 homolog (C7H21orf58; FDR = 0.0135; Figure 3b), gonadotropin releasing hormone 1 (GnRH1; FDR = 0.0216; Figure 3c) were upregulated in resident birds.
FIGURE 3
DEGs correlation with δ2H isotope and plasma testosterone levels
Macrophage mannose receptor 1 (MMR1) expression level showed a significant negative correlation to both δ2H isotope values (r2 = 0.7724, p < 0.0001) and T levels (r2 = 0.4574, p = 0.011; Figures 4a,b). In contrast, GnRH1 showed a positive correlation to both δ2H (r2 = 0.4692, p = 0.0098) and T levels (r2 = 0.3693, p = 0.0276; Figures 4c,d), indicating latitudinal differences in the HPG axis activation in anticipation to spring.
FIGURE 4
Discussion
We compared the early spring transcriptome of neuroendocrine tissues (hypothalamus) derived from two populations of free-living dark-eyed juncos that reside in the same non-breeding location but differ in timing of reproduction and whether or not they migrate. We identified only three differentially expressed genes (GnRH 1, MMR1, C7H21orf58) in the resident and migrant junco hypothalamic tissue transcriptome in early spring. The most relevant gene to this study was GnRH1 encoding transcript, which was expressed at higher levels in early breeding residents than in migrant breeding at higher latitudes. A critical step in the hypothalamic-gonadal axis (HPG) induction is the release of GnRH1 that binds to gonadotrophic cells in the anterior pituitary to release luteinizing hormone (LH) and follicle-stimulating hormone (FSH) that goes in the circulation and binds to the gonads to release sex hormones (). In addition to the differentially expressed genes (DEGs), we also found higher The hydrogen Prior research on this system using transcriptomics, proteomics, and quantification of specific target genes has shown variation in gene expression in tissues related to this phase lag in the timing of reproduction in migrants as compared to residents, including muscle, blood, and hypothalamus (; ; ). But to our knowledge, no other study has quantified gene expression in the hypothalamic-pituitary complex of two closely related sympatric populations of birds living in the wild during their divergence in reproductive timing and migratory preparedness under the same photoperiod and weather conditions in the wintering ground.
Latitude of origin is related to the timing of gonad recrudescence
Populations of dark eyed junco are distributed across different latitudes in North America and vary widely in their phenotypic (body size, feather color, eye color) and migratory behavior as well as their reproductive phenologies (). Those breeding at higher latitudes are more likely to be migratory, and in spring, these higher latitude birds increase body fattening in preparation for migration. Juncos breeding at lower latitudes are more likely to be resident, and these non -migratory birds start to grow their gonads earlier in the year in response to the increase in early spring daylengths (; ). After autumn migration, high and low latitude populations typically overwinter together in sympatry.
Feather stable hydrogen isotope ratios have been used as a proxy for breeding latitude in movement ecology studies (; ; ). Because juncos molt their feathers late in the breeding season before they migrate, analysis of feather hydrogen isotope (δ2H) ratios in sympatric dark-eyed juncos during the non-breeding season can provide an estimate of their breeding latitudes (; ; ). Birds with lighter hydrogen isotope ratios likely breed at higher latitudes and those with heavier hydrogen isotope ratios likely breed at lower latitudes (; ).
Several studies on junco populations have used feather hydrogen isotope ratios to show the latitudinal cline in spring preparedness of reproductive processes on the overwintering ground (; ; ). These studies comparing overwintering resident and migrant junco populations from two separate locations, Virginia and Colorado, indicate that lower latitude juncos begin to grow their gonads earlier, thus on shorter days, than high latitude juncos (; ). This difference is likely related to the difference in the threshold photoperiod required to stimulate the neuroendocrine system regulating reproduction ().
Temperate birds preparing to reproduce sense the change in day length and begin to rise in circulating testosterone, but activation of the hypothalamic-pituitary-gonad axis precedes that. Populations breeding at higher latitudes need relatively longer critical photoperiod than lower latitude birds to activate their HPG axis, grow gonads, and increase the circulating testosterone levels (). Several field and captive studies in junco populations show co-variation between hydrogen isotope ratios and elevation of circulating testosterone levels in response to external GnRH 1 challenge, suggesting a latitudinal cline in reproductive timing (; ; ; ). In our study, we found higher variability in circulating baseline testosterone and GnRH-1 transcript in resident juncos, suggesting a plausible effect of non-photic cues, such as food availability and exposure to predators, regulating the pace of reproductive maturity among resident birds (; ).
Latitude of origin is related to early spring hypothalamic GnRH 1 mRNA expression
In spring temperate-zone birds induce their gonads to grow in response to increasing day length. As noted, birds breeding at higher latitudes typically delay gonadal recrudescence until the days are longer because they migrate prior to reproducing (; ; ; ). This additional spring migratory phenology in migratory populations breeding at higher latitudes delays the activation of the HPG axis and uses all energy in preparing for the migratory process.
The preoptic area (POA), a key regulatory nucleus located in the hypothalamus, has neuronal populations that synthesize and release a neuropeptide, gonadotropin-releasing hormone (GnRH1). Transition from short to long day lengths is known to induce GnRH1 release from the hypothalamus that in turn stimulates the release of luteinizing hormone (LH) from gonadotroph cells in the anterior pituitary (; ; ). Increased levels of LH coming to the systemic blood circulation further stimulate gonadal recrudescence and synthesis of testosterone, a gonadal steroid (; ). Based on the studies on seasonal breeding songbirds, GnRH1 is the most critical top hierarchical molecule in stimulation of HPG axis with its elevation eventually leading to gonadal recrudescence.
In many avian species, levels of hypothalamic GnRH1 have been shown to be elevated significantly in the breeding season as compared to the non-breeding season (; ). Significant seasonal changes in the abundance of hypothalamic GnRH1 transcript/peptide have been found in several songbird species, such as house sparrow (Passer domesticus;; ), American tree sparrow (Spizella arborea; ), house finches (Carpodacus mexicanus; ); dark-eyed juncos (Junco hyemalis; ; ; ) and rufous-winged sparrows (Aimophila carpolis; ).
A prior study comparing captive resident and migrant juncos in a common garden demonstrated differential response of the HPG axis to the same day lengths despite exposure to similar photoperiod, food, and temperature in captive condition. At a given early spring day length, resident juncos showed higher GnRH mRNA expression levels than migrant juncos, providing evidence for difference in the HPG axis activation in individuals originating from different latitudes ().
In our study, the hypothalamic transcriptome data collected from early spring free-living overwintering junco populations showed higher GnRH mRNA expression in resident compared to migrant juncos. This indicates stimulation of GnRH neurons leading to earlier recrudescence in low-latitude residents as compared to high-latitude migratory juncos. The higher GnRH1 transcript levels in resident junco than migrants could be involved in earlier activation of the hypothalamic gonadal axis in residents than in migrant juncos. Thus, differential GnRH1 synthesis in junco populations experiencing similar environmental conditions suggests difference in the photic sensitivity of the neuroendocrine machinery regulating the timing of reproduction in junco populations breeding at different latitudes.
In a similar study of the neuroendocrine tissue proteome of resident and migrant junco populations collected in early spring showed key proteins involved in GnRH 1 transcription, translation, and post-translational modification, stability, and synthesis (). In addition, a proteome study on the juncos reported on here showed differential expression of Vasoactive Intestinal Peptide (VIP), a key neuropeptide molecule that has been implicated in the photoperiodic reproductive response and as a synchronizer of key circadian pacemaker cells in both birds and mammals (; ).
Research on female phenology as compared to males is lacking. A study comparing female junco populations in Virginia showed similar trends of higher expression of key reproductive genes in the ovary and ovary mass in resident juncos than high-latitude migrants (). Future work should focus on investigating the role of GnRH1 and other key genes discussed above in the seasonal divergence of reproductive timing, comparing sex and population-level differences in both free-living and captive conditions.
In addition to GnRH 1, we found two more genes differentially expressed in the hypothalamus of resident and migrant juncos. Macrophage mannose receptor 1-like (MMR1), which encodes a glycan-binding protein expressed in macrophages, was downregulated in resident juncos. The MMR1 role has been studied in several biological processes that include the regulation of circulating reproductive hormone levels, homeostasis, innate immune response, and neuroinflammation (). The protein encoded from Mannose receptor gene is involved in the clearance of circulating Luteinizing hormone (LH) from the bloodstream which could lead to a decrease in the circulating testosterone levels (). Though in our study, we found a negative correlation of MMR1 transcripts level with testosterone levels. But given the low sample size and high variability in resident testosterone levels, we would avoid deriving any mechanistic relationship between MMR1 transcript levels regulating circulating testosterone levels. Another gene, C21orf58 homolog (C7H21orf58), which is an uncharacterized protein, showed higher transcript levels in resident birds. Not much is known about the function of the C7H21orf58 gene and its role in regulating the seasonal phenologies of birds.
Conclusion
Based on transcriptomic differences in neuroendocrine tissue and circulating levels of testosterone, along with stable hydrogen isotope data revealing breeding latitude, we conclude that the HPG axis is activated earlier in the year in resident juncos that breed at lower latitudes as compared to migrant juncos that breed at higher latitudes. Future studies should focus on locating the neuronal cell bodies of GnRH 1 peptide and comparing its expression in winter non-breeding state and early spring divergent state when residents prepare to breed, while migrants prepare to migrate and delay breeding until they reach a higher latitude. In addition, future studies should extend to include other life-history states, the non-breeding wintering state, and the breeding phase at different latitudes.
The data presented in this study show the latitudinal effect on the timing of reproduction mediated via early activation of the HPG axis in birds that breed and develop at lower latitudes compared to higher latitude birds that breed in extreme long photoperiod. It will be interesting to explore the stimulation of the HPG axis at population-specific threshold photoperiod, where lower latitude resident juncos need relatively shorter photoperiod to induce gonads, and higher latitude migrants wait until they reach a critical photoperiod.
Statements
Data availability statement
The sequencing data are available in the NCBI Gene Expression Omnibus repository and are accessible through the GEO accession number # GSE305857.
Ethics statement
The animal study was approved by Virginia Department of Game and Inland Fisheries (permit # 052971)/ (# 15-026-17) Indiana University Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
DS: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review and editing. AF: Conceptualization, Data curation, Funding acquisition, Project administration, Resources, Writing – review and editing. EK: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The funds were provided by National Science Foundation (IOS- 1856423) and, Indiana University Faculty Research Support Program, (FRSP, account number # 2224640) to AMF and EDK. DS received pilot grants from Neuroscience Research priority Area (NRPA) and Advancing Innovative Research pilot grant from college of arts and sciences University of Kentucky.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1627516/full#supplementary-material
References
1
AshleyN. T.UbukaT.SchwablI.GoymannW.SalliB. M.BentleyG. E.et al (2014). Revealing a circadian clock in captive arctic-breeding songbirds, Lapland longspurs (Calcarius lapponicus), under constant illumination. J. Biol. Rhythms29, 456–469. 10.1177/0748730414552323
2
BauerC. M.FudickarA. M.Anderson-BuckinghamS.Abolins-AbolsM.AtwellJ. W.KettersonE. D.et al (2018). Seasonally sympatric but allochronic: differential expression of hypothalamic genes in a songbird during gonadal development. Proc. R. Soc. B285, 20181735. 10.1098/rspb.2018.1735
3
BertholdP. (1996). Control of bird migration. London, UK: Chapman and Hall.
4
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics, btu170.
5
BowenG. J. (2020). Gridded maps of the isotopic composition of meteoric waters. Available online at: http://www.waterisotopes.org.
6
ChmuraH. E.WingfieldJ. C.HahnT. P. (2019). Non-photic environmental cues and avian reproduction in an era of global change. J. Avian Biol., e02.
7
ChoR. N.HahnT. P.MacDougall-ShackletonS.BallG. F. (1998). Seasonal variation in brain GnRH in free-living breeding and photorefractory house finches (Carpodacus mexicanus). Gen. Comp. Endocrinol.109, 244–250. 10.1006/gcen.1997.7027
8
CormieA.ScwarczH.GaryJ. (1994). Relation between hydrogen isotopic ratios of bone collagen and rain. Geochim. Cosmochim.58, 819–830.
9
CummingsR. D. (2022). The mannose receptor ligands and the macrophage glycome. Curr. Opin. Struct. Biol.75, 102394. 10.1016/j.sbi.2022.102394
10
DawsonA. (2013). The effect of latitude on photoperiodic control of gonadal maturation, regression and molt in birds. Gen. Comp. Endocrinol.190, 129–133. 10.1016/j.ygcen.2013.05.017
11
DawsonA.GoldsmithA. R. (1983). Plasma prolactin and gonadotrophins during gonadal development and the onset of photorefractoriness in male and female starlings (Sturnus vulgaris) on artificial photoperiods. J. Endocrinol.97, 253–260. 10.1677/joe.0.0970253
12
DawsonA.GoldsmithA. R. (1997). Changes in gonadotrophin-releasing hormone (GnRH-I) in the pre-optic area and median eminence of starlings (Sturnus vulgaris) during the recovery of photosensitivity and during photostimulation. J. Reprod. Fertil.111, 1–6. 10.1530/jrf.0.1110001
13
DawsonA.FollettB. K.GoldsmithA. R.NichollsT. J. (1985). Hypothalamic gonadotrophin-releasing hormone and pituitary and plasma FSH and prolactin during photostimulation and photorefractoriness in intact and thyroidectomized starlings (Sturnus vulgaris). J. Endocrinol.105, 71–77. 10.1677/joe.0.1050071
14
DevicheP.MartinR. K.SmallT.SharpP. J. (2006). Testosterone induces testicular development but reduces GnRH-I fiber density in the brain of the house finch, Carpodacus mexicanus. General Comp. Endocrinol.147 (2), 167–174. 10.1016/j.ygcen.2005.12.024
15
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics29 (1), 15–21. 10.1093/bioinformatics/bts635
16
FudickarA. M.GreivesT. J.AtwellJ. W.StrickerC. A.KettersonE. D. (2016). Reproductive allochrony in seasonally sympatric populations maintained by differential response to photoperiod: implications for population divergence and response to climate change. Am. Nat.187, 436–446. 10.1086/685296
17
GirvanM. K.JonesJ.NorrisD. R.BargJ. J.KyserT. K.RobertsonR. J. (2007). Long-distance dispersal patterns of Male cerulean warblers (dendroica cerulea) measured by stable-hydrogen isotopes. Avian Conserv. Ecol.2, 3. 10.5751/ace-00173-020203
18
GreivesT. J.FudickarA. M.AtwellJ. W.MeddleS. L.KettersonE. D. (2016). Early spring sex differences in luteinizing hormone response to gonadotropin releasing hormone in co-occurring resident and migrant dark-eyed juncos (Junco hyemalis). General Comp. Endocrinol.236, 17–23. 10.1016/j.ygcen.2016.06.031
19
GwinnerE.CzeschlikD. (1978). On the significance of spring migratory restlessness in caged birds. Oikos30, 364–332. 10.2307/3543485
20
HahnT. P.BallG. F. (1995). Changes in brain GnRH associated with photorefractoriness in house sparrows (Passer domesticus). General Comp. Endocrinol.99 (3), 349–363. 10.1006/gcen.1995.1119
21
HobsonK. A. (1999). Tracing origins and migration of wildlife using stable isotopes: a review. Oecologia120, 314–326. 10.1007/s004420050865
22
HobsonK. A. (2003). Making migratory connections with stable isotopes. In: Avian migration. Editor BertholdP.GwinnerE.Sonnenschein.E.Berlin, Heidelberg: Springer, 379–391. 10.1007/978-3-662-05957-9_27
23
HobsonK. A. (2005). Stable isotopes and the determination of avian migratory connectivity and seasonal interactions. Auk122, 1037–1048. 10.1093/auk/122.4.1037
24
JacobsJ. D.WingfieldJ. C. (2000). Endocrine control of life-cycle stages: a constraint on response to the environment?Condor102, 35–51. 10.1093/condor/102.1.35
25
JohnstonR. A.PaxtonK. L.MooreF. R.WayneR. K.SmithT. B. (2016). Seasonal gene expression in a migratory songbird. Mol. Ecol.25, 5680–5691. 10.1111/mec.13879
26
KettersonE. D.GreivesT. J. (2025). Mechanisms matching timing to resources: comparisons of closely related seasonally sympatric, migratory and non-migratory populations. J. Avian Biol.2025, e03380. 10.1002/jav.03380
27
KettersonE. D.NolanJr V. (1976). Geographic variation and its climatic correlates in the sex ratio of eastern‐wintering dark‐eyed juncos (Junco hyemalis hyemalis). Ecology57, 679–693. 10.2307/1936182
28
KimmittA. A.HardmanJ. W.StrickerC. A.KettersonE. D. (2019). Migratory strategy explains differences in timing of female reproductive development in seasonally sympatric songbirds. Funct. Ecol.33, 1651–1662. 10.1111/1365-2435.13386
29
LanginK.ReudinkM. W.MarraP. P.NorrisD. R.KyserT. K.RatcliffeL. M. (2007). Hydrogen isotopic variation in migratory bird tissues of known origin: implications for geographic assignment. Oecologia152, 449–457. 10.1007/s00442-007-0669-3
30
LiQ.TamarkinL.LevantineP.OttingerM. (1994). Estradiol and androgen modulate chicken luteinizing hormone-releasing hormone-I release in vitro. Biol. Reprod.51, 896–903. 10.1095/biolreprod51.5.896
31
LiaoY.SmythG. K.ShiW. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics30 (7), 923–930. 10.1093/bioinformatics/btt656
32
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15, 550. 10.1186/s13059-014-0550-8
33
MagginiI.MetzgerB.VossM.VoigtC. C.BairleinF. (2016). Morphometrics and stable isotopes differentiate wintering populations of a migratory bird. Mov. Ecol.4, 20. 10.1186/s40462-016-0085-6
34
MillarR. P. (2005). GnRHs and GnRH receptors. Anim. Reprod. Sci.88, 5–28. 10.1016/j.anireprosci.2005.05.032
35
MishraI.SinghD.KumarV. (2016). Daily expression of genes coding for neurotransmitters in central and peripheral tissues of redheaded bunting: implication for circadian regulation of physiology in songbirds. Chronobiology Int.33 (3), 280–292. 10.3109/07420528.2016.1139587
36
NolanJr., V.KettersonE. D.CristolD. A.RogersC. M.ClotfelterE. D.TitusR. C.et al (2002). “Dark-eyed Junco (Junco hyemalis),” in The birds of North America online. Editor PooleA. (Ithaca, NY: Cornell Lab of Ornithology). Available online at: http://bna.birds.cornell.edu/bna.
37
NorrisD. O. (2007). Vertebrate endocrinology. 4th edn.Burlington, MA: Academic Press.
38
PowellL. A.HobsonK. A (2006). Enriched feather hydrogen isotope values for Wood Thrushes sampled in Georgia, USA, during the breeding season: implications for quantifying dispersal. Can. J. Zool.84 (9), 1331–1338. 10.1139/z06-114
39
PruittK. D.TatusovaT.KlimkeW.MaglottD. R. (2009). NCBI reference sequences: current status, policy and new initiatives. Nucleic Acids Res.37, D32–D36. 10.1093/nar/gkn721
40
RamenofskyM.CampionA. W.PérezJ. H.KrauseJ. S.NémethZ. (2017). Behavioral and physiological traits of migrant and resident white crowned sparrows: a common garden approach. J. Exp. Biol.220, 1330–1340. 10.1242/jeb.148171
41
RastogiA.KumariY.RaniS.KumarV. (2014). Neural correlates of migration: activation of hypothalamic clock (s) in and out of migratory state in the blackheaded bunting (emberiza melanocephala). PLoS One8 (10), e70065. 10.1371/journal.pone.0070065
42
ReinertB. D.WilsonF. E. (1996). Thyroid dysfunction and thyroxine dependent programming of photoinduced ovarian growth in American tree sparrows (Spizella arborea). General Comp. Endocrinol.103 (1), 71–81. 10.1006/gcen.1996.0095
43
RogersC. M.NolanV.KettersonE. D. (1993). Geographic variation in winter fat of dark-eyed juncos: displacement to a common environment. Ecology74, 1183–1190. 10.2307/1940488
44
RubensteinD. R.ChamberlainC. P.HolmesR. T.AyresM. P.WaldbauerJ. R.GravesG. R.et al (2002). Linking breeding and wintering ranges of a migratory songbird using stable isotopes. Science295, 1062–1065. 10.1126/science.1067124
45
SchroederA.MuellerO.StockerS.SalowskyR.LeiberM.GassmannM.et al (2006). The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol.31 (7), 3. 10.1186/1471-2199-7-3
46
SharmaA.DasS.SinghD.RaniS.KumarV. (2022). Differences in transcription regulation of diurnal metabolic support to physiologically contrasting seasonal life-history states in migratory songbirds. J. Ornithol.163 (1), 199–212. 10.1007/s10336-021-01926-5
47
SilverinB.MassaR.StokkanK. A. (1993). Photoperiodic adaptation to breeding at different latitudes in great tits. Gen. Comp. Endocrinol.90, 14–22. 10.1006/gcen.1993.1055
48
SinghD.KumarV. (2017). Extra-hypothalamic brain clocks in songbirds: photoperiodic state dependent clock gene oscillations in night-migratory blackheaded buntings, emberiza melanocephala. J. Photochem. Photobiol. B Biol.169, 13–20. 10.1016/j.jphotobiol.2017.02.016
49
SinghD.TrivediA. K.RaniS.PandaS.KumarV. (2015). Circadian timing in central and peripheral tissues in a migratory songbird: dependence on annual life-history states. FASEB J.29 (10), 4248–4255. 10.1096/fj.15-275339
50
SinghD.SwarupV.LeH.KumarV. (2018). Transcriptional signatures in liver reveal metabolic adaptations to seasons in migratory blackheaded buntings. Front. physiology9, 1568. 10.3389/fphys.2018.01568
51
SinghD.FudickarA.SmileyT.KettersonE. D. (2020). Comparative proteomics reveals molecular correlates of population-level variation in reproductive timing in a North American songbird. Authorea Preprints.
52
SinghD.ReedS. M.KimmittA. A.AlfordK. A.StrickerC. A.PollyP. D.et al (2021). Breeding at higher latitude is associated with higher photoperiodic threshold and delayed reproductive development in a songbird. Hormones Behav.128, 104907. 10.1016/j.yhbeh.2020.104907
53
SinghO.SinghD.MitraS.KumarA.LechanR. M.SingruP. S. (2023). TRH and NPY interact to regulate dynamic changes in energy balance in the Male zebra finch. Endocrinology164 (3), bqac195. 10.1210/endocr/bqac195
54
SmallT. W.SharpP. J.BentleyG. E.MillarR. P.TsutsuiK.MuraE.et al (2008). Photoperiod-independent hypothalamic regulation of luteinizing hormone secretion in a free-living sonoran desert bird, the rufous-winged sparrow (aimophila carpalis). Brain behav evol. Brain Behav. Evol.71 (2), 127–142. 10.1159/000111459
55
StevensonT. J.MacDougall-ShackletonS. A. (2005). Season-and age related variation in neural cGnRH-I and cGnRH-II immunoreactivity in house sparrows (Passer domesticus). General Comp.143 (1), 33–39. 10.1016/j.ygcen.2005.02.019
56
StevensonT. J.BernardD. J.BallG. F. (2009). Photoperiodic condition is associated with region-specific expression of GNRH1 mRNA in the preoptic area of the male starling (Sturnus vulgaris). Biol. Reproduction81 (4), 674–680. 10.1095/biolreprod.109.076794
57
WaiteT. A.CampbellL. G.BiologyO. (2006). Controlling the false discovery rate and increasing statistical power in ecological studies. Ecoscience13, 439–442. 10.2980/1195-6860(2006)13[439:ctfdra]2.0.co;2
58
WanamakerS. M.SinghD.ByrdA. J.SmileyT. M.KettersonE. D. (2020). Local adaptation from afar: migratory bird populations diverge in the initiation of reproductive timing while wintering in sympatry. Biol. Lett.16 (10), 20200493. 10.1098/rsbl.2020.0493
59
WingfieldJ. C.FarnerD. S. (1993). “Endocrinology of reproduction in wild species,” in Avian biology. Editors FarnerD. S.KingJ.ParkesK. (New York: Press), Vol. 9, 163–327.
60
WunderM. B.JehlJ. R.StrickerC. A. (2012). The early bird gets the shrimp: confronting assumptions of isotopic equilibrium and homogeneity in a wild bird population. J. Animal Ecol.81, 1223–1232. 10.1111/j.1365-2656.2012.01998.x
Summary
Keywords
seasonal reproduction, stable hydrogen isotope, migration, dark-eyed junco, hypothalamus, gnrh1
Citation
Singh D, Fudickar AM and Ketterson ED (2025) Seasonal divergence in reproductive timing on the verge of spring: comparing hypothalamic transcriptome of two seasonally sympatric North American songbird populations. Front. Physiol. 16:1627516. doi: 10.3389/fphys.2025.1627516
Received
12 May 2025
Accepted
12 August 2025
Published
03 September 2025
Volume
16 - 2025
Edited by
Sandra G. Velleman, The Ohio State University, United States
Reviewed by
Jamie M. Cornelius, Oregon State University, United States
Colin Guy Scanes, University of Wisconsin–Milwaukee, United States
Updates
Copyright
© 2025 Singh, Fudickar and Ketterson.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: D. Singh, devrajsingh23@uky.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.