The differential expression analysis of Brain Tissue
To analyze the role of photoperiods and circadian rhythms on gene expression patterns in the brain, we focused on samples collected at ZT0 (0600 hours local time) (Fig. 2A). We analyzed 12 RNAseq libraries from brain tissue, with six samples from each photoperiod treatment (N = 6). We generated an average of 17.5M reads per library, totaling 237933593.3 reads (Table 1). When aligned against the Galgal6 reference genome [ Ensembl 99 release version],11,867 genes showed detectable expression levels, and the dataset showed a common dispersion estimate of 0.105. Of the expressed genes, 607 genes were differentially expressed between the normal and extended photoperiod treatments at the ZT0 timepoint (6 a.m.). Among these differentially expressed genes (DEGs), 451 genes were downregulated, and 156 genes were upregulated in the normal photoperiods (Fig. 2B). The IPA analysis showed that Dopamine degradation, S100g family signaling pathway, Mitochondrial dysfunction, G-protein coupled receptor signaling, and Serotonin degradation were downregulated in extended photoperiods. PathfindR results showed the most enriched pathways such as Aldosterone-regulated sodium reabsorption, Endocrine and other factor-regulated calcium reabsorption, GABAergic synapse, Oxidative phosphorylation, Serotonergic synapse, Dopaminergic synapse and Circadian entrainment etc. with the Fold enrichment from range 1.5–6.6 (Fig. 3). The cluster enrichment analysis done with hierarchical clustering method grouped the pathways based on molecular function (Fig. 4). The master circadian clock is governed by the pituitary gland in birds and regulates the host’s physiological functions 43. The recent study reveals that the enzyme monoamine oxidase A (MAOa) plays a crucial role in dopamine metabolism, is regulated by the circadian clock, and impacts mood-related behaviors such as depression and addiction 44. In addition to synthesis and metabolism, dopamine release contributes significantly to circadian rhythms in behavior and physiology. For example, melatonin release from the pineal gland exhibits robust rhythmicity and is commonly used as a circadian marker in human studies. Melatonin release depends on the heteromerization of adrenergic receptors with dopamine D4 receptors, underscoring the pivotal role of dopamine in regulating pineal function 45. Also, dopamine neuroendocrine neurons express circadian clock genes, contributing to the circadian regulation in the host organism. This propounds that dopamine plays a dual role as a noteworthy input to the daily rhythms of circadian clock gene expression in specific parts of the brain 46.
Additionally, mitochondrial functioning regulated under circadian activity is a dynamic process crucial for maintaining cellular homeostasis and overall health 47. Mitochondria plays a role in regulating cellular redox balance, and circadian disruption impacts mitochondrial respiration, which may induce hypoxia regulated by HIF-1 signaling pathway 48. The perturbations in circadian rhythms can trade-off mitochondrial integrity and function, exacerbating oxidative stress and impairing cellular bioenergetic, which may drive a surge of reactive oxygen species (ROS) 49. These interactions underscore the interconnectedness of circadian-regulated processes, including mitochondrial function, and ROS signaling in orchestrating cellular physiology.
Table 1
provides a Quality control summary of RNASeq Reads obtained from samples of both normal photoperiod (NP) and extended photoperiods (EP).
QC | NP_ZT0 | NP_ZT24 | NP_ZT48 | EP_ZT0 | EP_ZT24 | EP_ZT48 |
Total M Seqs | 21.48 | 16.43 | 15.31 | 16.25 | 14.66 | 17.84 |
Length | 73 | 73 | 73 | 73 | 73 | 73 |
M Aligned | 18.81 | 14.8 | 13.24 | 14.32 | 12.78 | 14.55 |
% Aligned | 87% | 90% | 86% | 88% | 87% | 84% |
M Assigned | 10.54 | 8.29 | 6.93 | 6.96 | 6.34 | 6.99 |
% Assigned | 47% | 50% | 44% | 43% | 42% | 41% |
The genes following circadian rhythm
We analyzed expression intensities (CPM) of the samples of normal and extended photoperiod in DiscoRhythm to identify genes following diurnal oscillations(Fig. 5B). We found 577 oscillating genes with Cosinor and 417 with the JTK-cycle algorithm showing significant 24-hour rhythms (p < 0.05) (Fig. 5A). These two sets accounted for 723 unique genes, of which 269 were shared between the two algorithms (Supplementary file 1 & 2).
Further analyzing the 723 genes, these were found to be enriched for molecular function (MF), biological process (BP), and cellular component (CC). “The organic substance metabolic process”, “primary metabolic process,” and “cellular metabolic process” are the top three enriched biological processes with the highest gene counts involved. The top three enriched cellular components with the highest gene counts included “intracellular anatomical structure”, “organelle,” and “cytoplasm.”. The top enriched molecular functions were “heterocyclic compound binding," “organic cyclic compound binding," and “ protein binding”(Fig. 5C). To investigate differences between the NP and EP treatments, we pulled out the data for ZT0, as this time point represents the initiation of the light phase. These genes were predicted to be enriched in the canonical pathways FAK signaling, Phagosome formation, and Caveolar-mediated endocytosis signaling (activation Z-score > 2). The primary clock genes such as CLOCK, CRY1/2, BMAL1, ROR and TIPIN and four GPCR receptors FGFR, SUCNR1, HTR5, CXCL14, TACR3 showed robust 24-hr rhythms. Additionally, these receptor genes mediate G-protein coupled receptor signaling, which includes further tachykinin signaling, hydroxytryptamine receptor signaling, cytokine signaling and phagosome formation, which also showed 24-hr rhythms. These are essential genes for the overall immune health of the host 50–53. Activation of these genes shows that 24-hr rhythms are key for normal immune signaling. The activation of SUCNR1 suggests the invocation of succinate, a versatile compound that functions like hormones and cytokines54. SUCNR1 regulates metabolism by modulating circadian signaling and leptin expression, and SUCNR1 deficient adipocytes disrupts the body's ability to respond to leptin after feeding55. Succinate, a significant component within the tricarboxylic acid cycle also triggers the the production of inflammatory cytokines such as IL-1β, IL-6, IL-8, and TNF-α 56. Secondly, HTR signaling is responsible for serotonin release, and its expression is clock- dependent57. Serotonin released by the gut comprises about 95% of total serotonin and is linked with the microbiota composition 58. Serotonin also immunomodulates immune cells by secreting cytokines. Serotonin suppresses pro-inflammatory cytokines such as TNF-A and IL-1b 59. Another chemokine under clock control is CXCL14, which is secreted by a broad range of cells (immune and non-immune) and modulates inflammatory responses through GPCRs 52. Our findings of these interactions in the early life a chick, and during the window of microbiome assembly implies a close coordination of microbiota-mediated processes with circadian-regulated genes.
Correlation of gene expression and gut microbiota
To further investigate the specific interactions between circadian-regulated gene expression and the microbiota, we performed co-expression analysis using the tool WGCNA integrated the RNAseq data generated here and microbiota data from our previous study 30 .The hierarchical clustering (using ‘hclust’) showed the normal and photoperiod treatments formed two distinct clusters (Fig. 6A). The pickSoftThreshold was used for network topology analysis and set as power 5 (Fig. 6B ). Based on this power, the gene expression modules for each sample were clustered into four modules such as MEbrown, MEturquoise, MEblue, and MEgrey (Fig. 6C, and Fig. 7). Each module defines the correlation between each OTU and signifies the module-microbiome relationship (Fig. 6D) (Supplementary table 3 )
The relative abundance of specific OTUs was correlated with the Turquoise module, enriched for Cholesterol biosynthesis and Mitochondrial dysfunction. The relative abundance of OTU069 (Lactobacillus_72.9%. P = 0.007), OTU144 (Lachnospiraceae_NC2004, 71.9%, P = 0.008), Otu0134 (Defluviitaleaceae_UCG-011, 67.3%, P = 0.016), Otu0102 (Lachnospiraceae_FCS020_group, 62.1%, P = 0.03), Otu0113 (Ruminococcaceae_UCG-002, 61.9%, P = 0.03), Otu0054 (Clostridiales_vadinBB60_group, 60.4%, P = 0.037) was negatively correlated with the turquoise module. On the other hand MEturquoise is positively correlated with relative abundance of Otu0055 (Gastranaerophilales_unclassified, 58.6%, P = 0.04), OTU0012 (Alistipes, 76.4%, P = 0.003), Otu0105 (Ruminococcaceae_UCG-014, 58.6%, P = 0.04), Otu0152 (Ruminiclostridium_9, 60.8%, P = 0.03), Otu0063 (Ruminococcaceae_UCG-014, 68.1%, P = 0.014).
The correlations between the expressed genes and taxa in the Turquoise module highlight the intricate relationship between microbiome composition, cholesterol biosynthesis, and mitochondrial function. The implications of these correlations can be understood through the importance of cholesterol biogenesis for normal function. Altered cholesterol is associated with mitochondrial dysfunction and subsequent ATP deficiency, impacting neuronal function60. Additionally, mitochondrial fusion and ERK activity are crucial in cholesterol transport, affecting steroidogenesis 61 .
Mitochondrial dysfunction is also sensitive to cholesterol levels, as studies have shown that external cholesterol levels impact mitochondrial function and inflammatory responses 62. Furthermore, mitochondrial membrane transporter deficiencies lead to iron imbalance, affecting cholesterol biosynthesis and metabolic pathways 63. The circadian regulation of this crucial bidirectional interaction between cholesterol metabolism and mitochondrial function points to the important interactions indicated by this module. The correlation of these functions with specific microbial taxa in the turquoise module further suggests the circadian-regulated microbiota-homeostatic signaling crosstalk. In our previous studies 30, Rikenellaceae (Alistipes), Lachnospiraceae, and Ruminococcaceae were dominant in EP. Alistipes aids in producing saturated fatty acid 64 and the turquoise module's most enriched pathway was cholesterol biosynthesis and Mitochondrial dysfunction65. Lactobacillus abundance was correlated with the expression levels of NPC1L1, CYP7A1, and ABCG5, which regulate cholesterol in humans. These results suggest that the abundance of Lactobacillus might play a role in Cholesterol biosynthesis and Mitochondrial dysfunction pathways during microbiome assembly in chicken. The differential enrichment of Lactobacillus and its correlation with cholesterol biosynthesis can be interpreted through the lens of probiotic effects on lipid metabolism. Certain Lactobacillus strains can influence lipid metabolism and cholesterol levels. For example, abundant probiotic Lactobacillus strains can impact cholesterol assimilation and bile acid levels 66,67. Additionally, Lactobacillus strains lower total bile acids in serum, indicating a potential role in cholesterol metabolism 68. Furthermore, Lactobacillus strains like Lactobacillus fermentum and Lactobacillus plantarum have probiotic effects on hypercholesterolemia, affecting cholesterol levels in animal models 69 70. These strains have been shown to lower cholesterol levels in vitro and in vivo, suggesting underlying mechanisms for these taxa to modulate cholesterol metabolism 70.
The blue module was negatively correlated with Otu0116 (Ruminococcaceae_UCG-004, 67.5%, P = 0.016). However, it is positively correlated with Otu0046 (Ruminococcaceae_UCG-014, 58.1%, P = 0.047), Otu0142 (Lachnospiraceae, 59.01%, P = 0.043), Otu0033 (Flavonifractor, 63.5%, P = 0.03), OTU0039 (Clostridiales vadinBB60, 86%, P = 5.68E-09) and OTU0 049 (Eubacterium coprostanoligenes, 98%, P = 0.0002).
The most enriched canonical pathways in blue modules were Dopamine receptor signaling, Melatonin degradation, Myelination signaling pathway, and WNT/𝛃-catenin signaling pathway analyzed in IPA (Fig. 6E).
The enrichment of canonical pathways such as Dopamine receptor signaling, Melatonin degradation, Myelination signaling pathway, and WNT/𝛃-catenin signaling pathway in the blue module points to the functional implications of these microbiota members (Fig. 6E). These pathways play crucial roles in various physiological processes, including neurotransmission, circadian rhythm regulation, myelin formation, and cell signaling. The negative correlation with Otu0116 (Ruminococcaceae_UCG-004) may suggest a potential regulatory role of this taxon in modulating the pathways associated with Dopamine receptor signaling, Melatonin degradation, Myelination signaling pathway, and WNT/𝛃-catenin signaling pathways. Conversely, the positive correlations with Otu0046, Otu0142, Otu0033, Otu0039, and Otu0049 indicate the potential contribution of these microbial taxa to the activation or enhancement of these pathways. The presence of Ruminococcaceae and Lachnospiraceae taxa in the positive correlations is noteworthy, as these bacterial families are known for their roles in promoting gut health and metabolism. Flavonifractor, Clostridiales vadinBB60, and Eubacterium coprostanoligenes are also important gut microbiota members with potential implications for host physiology and metabolism. In chickens, there is also evidence for the the influence of these taxa on behavior: cecal microbiota transplantation with Lachnospiraceae and Ruminococcaceae UCG-005 in chickens at an early promoted aggressive behavior in recipient chickens by controlling the functions of the catecholaminergic and serotonergic systems in the brain71.
The grey module was negatively correlated Otu0181 (Ruminococcaceae, 66.%2, P = 0.0188), Otu0080 (Oscillibacter, 64.8%, P = 0.022), Otu0068 (Lachnospiraceae, 60.9%, P = 0.03) significant P < 0.05. The grey module is positively correlated with Otu0029 (Ruminococcaceae_UCG-014, 58.2%, P = 0.04), Otu0151 (Christensenellaceae_R-7_group, 60.4%, P = 0.03), Otu0100 (Lachnospiraceae, 66.2%, P = 0.01) significant. The enriched pathways in this module are the Visual Cycle, ABRA signaling pathway, Apelin cardiomyocyte signaling pathway, Regulation of actin-based motility by Rho (Fig. 6E).
The negative correlations with Otu0181 (Ruminococcaceae), Otu0080 (Oscillibacter), and Otu0068 (Lachnospiraceae) suggest a potential regulatory role of these taxa in modulating the pathways associated with the Visual Cycle, ABRA signaling pathway, Apelin cardiomyocyte signaling pathway, and Regulation of actin-based motility by Rho within the grey module. On the other hand, the positive correlations with Otu0029 (Ruminococcaceae_UCG-014), Otu0151 (Christensenellaceae_R-7_group), and Otu0100 (Lachnospiraceae) indicate a potential contribution of these microbial taxa to the activation or enhancement of these pathways. The apelin-APJ pathway can directly antagonize vascular disease-related Ang II actions 72, providing insights into the regulatory mechanisms of the apelin signaling pathway in cardiovascular health. The microbiota associations suggest intriguing hypotheses about how early-life microbiome interactions with host homeostatic pathways can affect multiple major functions.
The brown module is negatively correlated significantly (P < 0.05) with Otu0039 (Clostridiales vadinBB60_group, 60.3%, P = 0.03), Otu0037 (Clostridiales_vadinBB60_group, 58.9%, P = 0.04), Otu0074 (Ruminococcaceae_UCG-014, 58.07%, P = 0.04. However is positively correlated with Otu0118 (Anaerotruncus, 57.7%, P = 0.03), Otu0183 (Ruminiclostridium_9, 61.26%, P = 0.03), Otu0066 ([Eubacterium]_hallii_group, 64.8%, P = 0.022). Top canonical pathways in this module are the Virus Entry via Endocytic Pathways, GABA Receptor Signaling, Adipogenesis pathway, Semaphorin Neuronal Repulsive Signaling Pathway, HMGB1 Signaling (Fig. 6E).
The interactions observed within the brown module, as indicated by the significant correlations with specific OTUs, suggest a complex interplay. The negative correlations with OTUs such as Otu0039 (Clostridiales vadinBB60_group), Otu0037 (Clostridiales_vadinBB60_group), and Otu0074 (Ruminococcaceae_UCG-014) may indicate a potential regulatory role of these taxa in modulating pathways associated with the Virus Entry via Endocytic Pathways, GABA Receptor Signaling, Adipogenesis pathway, Semaphorin Neuronal Repulsive Signaling Pathway, and HMGB1 Signaling within the brown module. These interactions highlight links between gut dysbiosis, neurotransmitter disturbance, and inflammatory responses, shedding light on the potential impact of microbial taxa on behavior and neurotransmitter metabolism. Altogether, these modules indicate intriguing correlations between the microbiota and gene expression data in the context of circadian rhythms, and point to the foundational role these interactions may have on later life metabolic and immune health.
KEGG pathway analysis
The most enriched pathways found in our IPA and WGCNA analysis were GPCR signaling, Dopamine signaling, Serotonin signaling, Melatonin degradation, Myelination signaling, tachykinin signaling, and Mitochondrial dysfunction. We used the KEGG terms to understand the key molecules in these pathways to identify the potential mechanisms of importance in the microbiota-host crosstalk. Across these pathways, most enriched genes were linked with calcium signaling, MAPK Signaling, circadian entrainment and chemokine signaling (Fig. 8). The clustering in Fig. 3 shows that Enhanced functional processes in GPCR signaling, Dopamine signaling, Serotonin signaling, Melatonin degradation, Myelination signaling, tachykinin signaling, and Mitochondrial dysfunction collectively invoke a complex and intricate interplay of various signaling pathways. While many pathways could be impacted, calcium signaling is one of the key pathways. Calcium signaling is fundamental in various cellular responses to environmental stimuli, including plant defense systems 73. It also regulates cell cycle progression in response to abiotic stress 74. Furthermore, reactive oxygen species (ROS) in mitochondria trigger monoamine-induced calcium signals, influencing physiological and pathophysiological responses to dopamine 75. Calcium signaling is versatile in numerous cellular functions and crucial for translocating nuclear proteins like PKC-γ 76. The calcium signaling is a critically important phenomenon that regulates synaptic activity in the nervous system. Also, synaptic activity is the key feature of circadian oscillations, and changes in circadian rhythms alter synapse numbers. Calcium is the most important signaling molecule that plays a role in nervous excitability and regulation of biological clock 77. Calcium signaling is intricately linked to circadian rhythms, playing a significant role in regulating the molecular clock and coordinating various physiological processes throughout the day. Circadian rhythms in the SCN neurons are synchronized with calcium rhythms, indicating a causal relationship between intracellular calcium dynamics and the generation of circadian rhythms 78. Furthermore, regulating voltage-dependent calcium channels (VDCCs) by circadian mechanisms is crucial for maintaining rhythmic clock gene expression in the SCN, highlighting the importance of calcium signaling in the molecular machinery of circadian clocks 79. Additionally, mitochondrial calcium signaling has been implicated in mediating rhythmic extracellular ATP accumulation in SCN astrocytes, further emphasizing the role of calcium in coordinating circadian processes at the cellular level 80. These enriched pathways point towards key homeostatic processes. These pathways play a role in intracellular calcium signaling, which is crucial for activating lymphocytes and neurological functions (KEGG hsa04020). The gene regulatory network of differentially expressed genes involved in these pathways are shown in Fig. 9. The genes in the maroon-colored PPI network (Fig. 9) were found to be upregulated and blue-colored were downregulated in extended photoperiods.
The interaction between circadian rhythms and the gut microbiome may include calcium signaling as a key mediator. Calcium signaling regulates neurotransmitter release and cellular functions essential for circadian rhythm generation 81. The influx of calcium ions into the presynaptic terminal triggers the release of neurotransmitters into the synaptic cleft, facilitating signal transmission between neurons 82. This process is tightly regulated and involves various proteins and signaling pathways. For example, synaptotagmins, calcium sensors on synaptic vesicles, are vital for vesicle fusion with the presynaptic membrane and subsequent neurotransmitter release 83. Calcium signaling also modulates synaptic plasticity, where changes in calcium levels can impact the strength of synaptic connections and neuronal communication 84. On the other hand, the gut microbiome can be influenced by factors such as diet, disease, and probiotics, impacting cognitive function and overall health85. Given the bidirectional communication between the gut microbiome and the host, calcium signaling may serve as a signaling pathway through which circadian rhythms interact with the gut microbiome. Furthermore, the gut microbiome's role in modulating host metabolism and immune function suggests a potential link between calcium signaling, circadian rhythms, and gut microbiome-mediated effects on health and disease.