Light-induced asymmetries in embryonic retinal gene expression are mediated by the vascular system and extracellular matrix

Left–right asymmetries in the nervous system (lateralisation) influence a broad range of behaviours, from social responses to navigation and language. The role and pathways of endogenous and environmental mechanisms in the ontogeny of lateralisation remains to be established. The domestic chick is a model of both endogenous and experience-induced lateralisation driven by light exposure. Following the endogenous rightward rotation of the embryo, the asymmetrical position in the egg results in a greater exposure of the right eye to environmental light. To identify the genetic pathways activated by asymmetric light stimulation, and their time course, we exposed embryos to different light regimes: darkness, 6 h of light and 24 h of light. We used RNA-seq to compare gene expression in the right and left retinas and telencephalon. We detected differential gene expression in right vs left retina after 6 h of light exposure. This difference was absent in the darkness condition and had already disappeared by 24 h of light exposure, suggesting that light-induced activation is a self-terminating phenomenon. This transient effect of light exposure was associated with a downregulation of the sensitive-period mediator gene DIO2 (iodothyronine deiodinase 2) in the right retina. No differences between genes expressed in the right vs. left telencephalon were detected. Gene networks associated with lateralisation were connected to vascularisation, cell motility, and the extracellular matrix. Interestingly, we know that the extracellular matrix—including the differentially expressed PDGFRB gene—is involved in morphogenesis, sensitive periods, and in the endogenous chiral mechanism of primary cilia, that drives lateralisation. Our data show a similarity between endogenous and experience-driven lateralisation, identifying functional gene networks that affect lateralisation in a specific time window.


Materials and methods
Samples and study design. Sixty-three freshly fertilized eggs of Ross 308 (Aviagen) strain were stored in darkness at 4 °C until the beginning of incubation. Having 21 eggs in each experimental group provided us with enough embryos to dissect, despite the initially unknown number of unfertilised eggs and sex composition of the sample. Eggs were incubated at a constant temperature of 37.7 °C and 40% humidity and rotated automatically in complete darkness from embryonic day 0 (E0) to embryonic day 17 (E17) at 9:00 am. Subsequently, eggs were randomly assigned to experimental groups: in the Dark condition eggs continued to stay in darkness, eggs in the 24-h light exposure condition received light from E17 9:00 am to ~ E18 12:00 pm (i.e., 24-27 consecutive hours of light), and eggs in the 6-h light exposure condition received light from E18 4:00 am to ~ E18 12:00 pm (i.e., www.nature.com/scientificreports/ 5-8 consecutive hours of light) (Fig. 1C). This timing of exposure and dissection was chosen to make sure light exposure was long enough to induce functional lateralisation 55 and at the same time target responses subsequent to those of immediate early genes [72][73][74] . Light was delivered to eggs through five white-light LEDs (18 lumens per LED) placed at 24 cm from the eggs. We conducted 44 dissections on E18 between 9 am-12 pm (unfertilised and surplus eggs were not dissected): 14 in the dark condition, 15 in the 6-h exposure condition, 15 in the 24-h exposure condition. After molecular sexing and quality control analyses on RNA samples, we chose 5 male samples for each light regime, based on the best quality control indexes. After quality control analyses on the RNA samples, we randomly choose 5 male samples for each of the light regimes based on the best quality control indexes.

Dissection.
A total of 44 embryos were dissected. Experiments with bird embryos did not require specific protocol approval when these studies were carried out (i.e., 2016), as per the Italian and European Community laws for the ethical treatment of animals (Legislative decree law 26/2014). Our reporting follows ARRIVE guidelines. For each embryo, the egg was retrieved from the incubator, placed at − 20 °C for 2 min and then immediately opened by cutting along the widest part of the shell with scissors. Embryonic membranes were then disrupted and the embryo was culled by rapid decapitation followed by the removal of both eyes and then brain. We included embryos that macroscopically corresponded to developmental Stage 44 and took notes of potential Only right eye is directly exposed to light. (B) Representation of the visual system in chicks. The two arrows that depart from the retina indicate the tectofugal and thalamofugal pathways (see 55 ). The left tectofugal pathway (orange pathway via right retina, left optic tectum, left nucleus rotundus (NR) and left entopallium in the left telencephalon) and left thalamofugal pathway (red pathway via right retina, left thalamic nucleus geniculatus lateralis pars dorsalis (GLd), left Wulst in the left telencephalon and right Wulst in the right telencephalon) are shown. The right tectofugal and thalamofugal pathways (fed by the left retina) are represented in blue. The regions used in the RNA-seq analysis (retina and telencephalon) are coloured. (C) Schematic representation of the experimental phases. Eggs were incubated in the dark from embryonic day 0 (E0) to embryonic day 17 (E17). On day E17, eggs were randomly assigned to one of three different experimental treatments. One group continued incubation in darkness (Dark), one group was exposed to 6 h of light (6-h light) and one group was exposed to 24 h of light (24-h light www.nature.com/scientificreports/ differences 61 . The dissection of retinas and telencephalons were performed in parallel by two experimenters to minimize sample collection time. Brain tissue was macroscopically subdivided in two regions: telencephalon and remainder of the brain. Of these regions, only telencephalon was further subdivided into left and right samples and processed for RNA-seq. Macrodissected brain samples were frozen immediately on dry ice and stored at − 80 °C until RNA extraction. For retina dissection, first, eyelids were cut away and the whole eye was carefully removed without piercing the eye globe. Then, dissection of the retina from the enucleated eyes was performed by cutting the eyes in half to expose the vitreous and retina. The translucent gelatinous vitreous humour was carefully removed with forceps exposing the retina. Retina and pigmented epithelium were removed from the eye cup and placed immediately on dry ice and stored at − 80 °C until RNA extraction.

Molecular sexing.
Since anatomic asymmetries in the thalamofugal pathway 75 and behavioural asymmetries 76,77 are more pronounced in male chicks compared to females (reviewed in Ref. 78 ), we focused our study on males only. Sexing was performed according to Ref. 79 .
RNA extraction and quality control. Total RNA was extracted and DNAse-treated using RNeasy mini kit (Qiagen) with DNAse treatment, following the manufacturer's protocol. RNA concentrations were measured www.nature.com/scientificreports/ using the Qubit RNA Broad-Range Assay kit and RNA quality was analysed on an Agilent Bioanalyzer 2100. Only samples with RNA integrity numbers (RIN) ≥ 9 were considered for RNA-seq inclusion.
RNA-seq analysis. RNA-seq data were generated from a total of 60 samples (5 embryos × 3 experimental groups × 2 tissues × 2 hemispheres) on an Illumina HiSeq 4000 at the Genomic facilities of the Wellcome Trust Human Genetic Centre (WTHGC), University of Oxford. Five biological replicates were used for each tissue (telencephalon and retina), side (left and right) and condition (dark, 6-h and 24-h light exposure). One µg of input RNA was used for the standard TruSeq Stranded mRNA Library Preparation Kit that includes isolation of mRNA by polyA-selection. The libraries were then pooled according to their tissue, side and treatment group. The sequence libraries had a median number of 11,232,013 reads before trimming (10,261,901 reads after trimming). All data were checked for quality with FastQC (FastQC v0.11.5, 2017 https:// www. bioin forma tics. babra ham. ac. uk/ proje cts/ fastqc/) and contaminants with FastQ Screen (v0.10.0), then trimmed for adapter sequences with Trimmomatic (v0.36) 80 . The statistics of number of reads per sample are shown in Supplementary Table 4.
Differential gene expression analysis. Surviving reads were quantified by mapping against the Ensembl galGal6 chicken transcriptome (downloaded from ftp:// ftp. ncbi. nlm. nih. gov/ genom es/ all/ GCF/ 000/ 002/ 315/ GCF_ 00000 2315.6_ GRCg6a/ GCF_ 00000 2315.6_ GRCg6a_ rna. fna. gz) using Salmon (v0.12.0) 81 . Mapping efficiency ranged from 76 to 81%. Transcript read counts were collapsed to genes with tximport 82 in R version 4.0.3 (https:// www.R-proje ct/ org), then counts and fragment lengths were imported into DESeq2 (1.30.0) 83 . Principal components analysis (PCA) was applied to visualize variance in the data, and one individual in the 24-h light exposure treatment group was excluded from further analysis because it was an outlier and it was annotated as underdeveloped at the time of dissection. Sequencing replicates of each sample library were similar by PCA, so they were combined to increase statistical power. Differential gene expression analysis was conducted separately for retina and telencephalon tissues. With DESeq2, we compared gene expression in left and right sides of the same tissue in the same treatment conditions, adding a factor to our model to control for individual effects, and we also looked for aggregate differences in gene expression between treatment groups (reference code in Available data for details). An adjusted p-value cut-off of p < 0.05 was applied for all comparisons. Analysis was conducted with R (https:// www.R-proje ct. org/) and RStudio (http:// www. rstud io. com). The packages dplyR (https:// CRAN.R-proje ct. org/ packa ge= dplyr), biomaRt 84 , tximport 82 and independent filtering was employed to optimize detection of genes below this threshold. Shrinkage estimation was applied to moderate fold-changes of low abundance genes to facilitate ranking for downstream analysis of gene ontology.

Results
Differential gene expression: lateralised and overall effects. To identify genes differentially expressed in the right and left hemisphere in the absence or presence of asymmetric light exposure, we compared left vs right gene expression in both tissues (retina and telencephalon) in the Dark, 6-h light and 24-h light conditions, using a within-subjects design ( Fig. 2A). In the retina, we identified 56 differentially expressed (DE) genes (adjusted p-value < 0.05, Supplementary Table 1) in the right vs left side after 6 h of light exposure (Fig. 2C). All differentially expressed genes were downregulated in the right vs. left retina (see also Ref. 85 ). Besides genes involved in the circulatory system development and extracellular matrix/cell motility (see "Gene set enrichment analysis" section below) we identified two differentially expressed genes of interest. The first, DIO2, is involved in sensitive periods 66,67 and in the development of the retina in the chicken and other vertebrates 68 . The second, PAX2, is a transcription factor that plays an important role in the glia during retinal development in the chicken embryo 86 and in eye development in vertebrates 87 . We found no significant left-right differences in gene expression in either darkness or after 24 h of light exposure in the retina. No significant differences in gene expression were identified in the telencephalon.
To identify the overall, non-lateralised effect of light exposure, we compared gene expression of both brain sides (i.e., right plus left retina, right plus left telencephalon) between different light exposure conditions (Dark vs 6-h light, Dark vs 24-h light, 6-h light vs 24-h light) (Fig. 2B). In both tissues, the largest effect has been found in the Dark vs 6-h light exposure comparison (52 differentially expressed (DE) genes in the retina, 23 in the telencephalon). In the 6-h vs 24-h comparison we found 4 DE genes in the retina, 3 in the telencephalon; in the dark vs 24-h exposure comparison we found 13 DE genes in the retina (4 upregulated, 9 downregulated, including CAPN15, that is involved in the development of the visual system) and 0 in the telencephalon (see Supplementary Table 2A-E). The reduction of DE genes after more than 6 h of light exposure corroborates behavioural and neuroanatomic evidence that, in different species, embryonic light exposure has a transient effect (e.g. see 45 for different time windows in chicks, see 46 for zebrafish). Behavioural evidence suggests that, in domestic chick embryos (E19), this self-terminating sensitive period ends after 6 h of light exposure, during the last three days of incubation 39 . Gene set enrichment analysis. To identify biological pathways connected to significant asymmetrical light exposure, we used the differentially expressed genes identified above (Retina dark vs 6-h exposure) and performed a gene enrichment analysis via g:Profiler 88,89 , with 16,356 protein-coding genes used in the model as the background, and FDR at 1%. In the Retina dark vs 6-h exposure comparison, 32 biological process GO terms, 14 cellular component terms and 23 molecular function GO terms were enriched (Fig. 3A). A partial list of enriched terms is shown in Fig. 3B (for the complete list see Supplementary Table 3).
The enriched biological processes include two systems, namely the development of the vascular/circulatory system and the formation of the anatomical structure, which highlight an involvement of ontogenetic plasticity. www.nature.com/scientificreports/ The enriched biological processes relate to two main systems: vascular/circulatory system development and anatomical structure formation. In line with these results, in studying the functional genetics of handedness and lateralisation, Schmitz et al. 38 found enrichment on GO terms related to anatomic structure development, including cardiovascular system development, artery development and epithelial tube morphogenesis. The enriched cellular components include the extracellular matrix (that provides physical scaffolding for the cellular constituents and can initiate biochemical and biomechanical cues required for tissue morphogenesis, differentiation and homeostasis), plasma membrane (GO:0005886), haptoglobin-haemoglobin complex (GO:0031838), cell periphery (GO:0071944) and haemoglobin complex (GO:0005833) as top hits. Finally, the enriched molecular functions include metal ion binding (GO:0046872), cation binding (GO:0043169) and platelet-derived growth factor binding (GO:0048407), haptoglobin binding (GO:0031720) and oxygen carrier activity (GO:0005344) and metallocarboxypeptidase activity (GO:0004181) as top hits. Using the same approach, no functional enrichment was found in other tissues and conditions. KEGG 90 pathway analysis confirmed the significant enrichments for extracellular matrix (ECM) receptorinteraction, that is important in tissue morphogenesis and in the maintenance of tissue and function, and focal adhesion, that is involved in cell motility proliferation and differentiation. Interestingly, the membrane of primary cilia has been associated with receptors for ECM proteins 26,91 . Similarly, in studying the functional genetics involved in handedness with KEGG pathway analysis, Schmitz et al. found enrichment for ECM-receptor interaction and focal adhesion. Moreover, they found enrichment in the TGF-beta signalling pathway. This pathway is connected to PDGFRB (platelet-derived growth factor receptor β), which we found to be differentially expressed in right vs left retina at 6 h of light exposure. This gene, essential for the development of the cardiovascular system, helps in the rearrangement of the actin cytoskeleton and is involved in the dynamic primary cilia mechanism 24,92 . This further evidence points toward the implication of primary cilia in the experience-driven, light-triggered gene pathway of lateralisation initiated in the eye 13 .
Overall, gene set enrichment analysis shows an involvement of blood vessels/circulatory system development and processes related to morphogenesis in response to 6 h of light exposure in the retina. Moreover, we detected a temporal-specific effect in the retina present after 6 h of light stimulation but not after 24 h. No enrichment was present in other experimental conditions (bilateral tissue sampled at different light-exposure times) and in the telencephalon.   www.nature.com/scientificreports/ Interaction network. To investigate the protein interaction network of the differentially regulated genes involved with embryonic light exposure (left-right asymmetry after 6-h light exposure), we analysed the differentially expressed genes using STRING v11.0 93 and visualised its output with Cytoscape 94,95 . We analysed the full network (physical and non-physical interactions) with all available interaction sources and medium confidence cut-off (required interaction score ≥ 0.4). Of the 56 differentially expressed genes for Retina 6-h light exposure, 53 could be found in STRING. The 53 nodes in the overall network produce 40 interactions (vs 6 expected interactions if the nodes were selected at random), with a significant protein-protein interaction enrichment (p < 1.0e−16), involving 25 genes. Two interaction subnetworks have been identified (see Fig. 4).

A.
The main network (left in Fig. 4) includes nodes involved in the vascular/circulatory system (blue circle) and in connective tissue/cell adhesion/extracellular matrix (yellow outline); the haemoglobin-complex subnetwork (right) is linked to the vascular/circulatory system.

Discussion
The mechanisms that drive the ontogenesis of lateralisation are far from being understood. Due to ease of egg manipulation and brain architecture, birds are prominent models to investigate the role of environmental factors in lateralisation 11 and eye development 59 . We used chicken embryos to identify genetic pathways involved in lateralisation induced by asymmetrical light stimulation. We compared gene expression in the right and left retinas and telencephalon in embryos exposed to different light regimes. In natural settings, light exposure is regulated by the parents, that systematically interrupt incubation in the last days before hatching, as shown in pigeons and chickens 96,97 . In these periods, light reaches the egg and embryo, enhancing ontogenetic plasticity.
After exposing embryos to darkness, 6-h light or 24-h light, we dissected samples at the same developmental stage (embryonic day E18) and compared the effects of endogenous and stimulus-driven lateralisation induced by light in the right and left brain. We identified significantly differentially expressed genes in the retina in the 6-h light exposure condition ( Fig. 2A). Because all samples have been collected at the same developmental stage, we conclude that the observed differential gene expression between left and right retina is due to light exposure.
All differentially expressed genes were downregulated in the right vs left retina, suggesting a potential inhibitory effect of light during ontogenesis. In line with this scenario, Manns et al. 85 have previously shown that embryonic light stimulation in pigeon induces asymmetric inhibition of TrkB ⁄Ras activity (related to morphogenesis) in the tectum fed by the right eye. These cellular asymmetries emerge in the tectum only after hatching. A downregulation in the extracellular matrix (ECM), with consequent inhibition of plasticity 98 , is consistent with an inhibitory effect of light. In the present work, we exposed embryos to at least 6 h of light, thus avoiding the signature of peak expression of immediate early genes. Future work should investigate early effect of light exposure (before 6 h of exposure) on neurons and ECM, as well as the connection between retinal light exposure and inhibition in other brain areas.
Our findings support the idea that differential gene expression induced by light exposure terminates between 6 and 24 h of light exposure in the retina. The lack of differences observed after 24 h of light exposure, while differences are detected at 6 h of light exposure, confirms previous findings about the presence of transient sensitive periods in the development of the nervous system in chicks (for chicks see 62,74,99 ), and other vertebrates 66,100 (see 101 for developmental changes in human foetuses). Moreover, the timing of this effect is in line with the fact www.nature.com/scientificreports/ that inversion of lateralisation in chicks has been achieved only between 2.5 and 6 h of light exposure at E19/ E20 39 . Overall, this pattern of data suggests that, between 6 and 24 h of light exposure the sensitive period of neural plasticity that stabilises lateralisation in the retina self-terminates. This hypothesis is corroborated by the significant differential expression of DIO2 (iodothyronine deiodinase 2) in the 6-h light exposure but not in the 24-h light exposure. DIO2 is a thyroid hormone converting enzyme with an essential role in neurodevelopment and the regulation of the sensitive periods in mammals and birds 66 . DIO2 has a role in the development of the retina in different vertebrates, including the chicken, where is involved in spatiotemporally specified waves of cone differentiation 68 . In the chick retina, the expression of thyroid hormone regulators follows a series of spatiotemporal waves, moving from the central retina to the periphery. It has been suggested that DIO2 likely blocks/postpones photoreceptor differentiation, thus increasing photoreceptors in specific areas of the retina 102 . More in general, modulation of DIO2 can finely tune the availability of triiodothyronine (and hence the presence of thyroid hormones) at specific times during brain development. This gene regulates the start of the sensitive period for filial imprinting in chicks (the self-terminating mechanism that mediates affiliative learning), as well as the regulation of predisposed preferences in domestic chicks, such as the preference for biological motion 67,103 and objects that move changing in speed. Evidence that light exposure has a transient effect in a specific time window comes also by the comparison of the samples from both sides of each tissue (right plus left retina, right plus left telencephalon). In these comparisons, in both tissues we identified more differentially expressed genes in the darkness vs 6-h light exposure than in the 6-h light exposure vs 24-h light exposure conditions (23 vs 3 DE genes in the telencephalon; 52 vs 4 DE genes in the retina, see Fig. 2B).
While the retina showed right vs left differential gene expression at 6 h of light exposure, we did not detect any difference between sides in the telencephalon. A similar lack of significant left-right differences in telencephalon has been previously reported in fish 104 and in human foetuses in the perisylvian area (neocortex), an area with high interhemispheric asymmetry 105 , while only mild and age-dependent asymmetries have been recently found in human foetuses 101 ).
The telencephalon is a large portion of the chicken pallium fed by both the thalamofugal and tectofugal pathways (see Fig. 1B. These brain areas are homologous to the geniculocortical and extrageniculocortical systems in mammals 11 ). Lack of lateralised gene expression at E18 can hence be due to different factors. Since we dissected a large portion of the brain, the lack of differential expression could be due to more regionally restricted lateralised effects mediated by light. In addition, light exposure may have its maximum effect in the telencephalon more than 24 h after the beginning of light stimulation. Moreover, the stimulation of both tectofugal and thalamofugal pathways may activate both hemispheres, since the thalamofugal pathway has both ipsilateral and contralateral projections in the telencephalon, or some genes involved in ontogenetic development and plasticity might not have transcriptional modulation. Further experiments should clarify whether more anatomically restricted areas of the embryonic telencephalon are affected by light exposure and on what time course.
Gene enrichment analyses have revealed two main pathways associated with lateralised differences induced by light stimulation: the vascular/circulatory system and the extracellular matrix (ECM). An interesting overlap between gene networks involved in the vascular system and the ECM/connective tissue appears in the interaction network analysis (Fig. 4). It is worth discussing the connection between the vascular system, ECM and lateralisation in more detail. As for the vascular system, one of the first macroscopic asymmetries observed during vertebrate development is the right-side looping of the heart tube (see 70 for chicks, reviewed in Ref. 106 ). More in general, organogenesis is a highly asymmetric process and the biophysical mechanisms engaged in organogenesis involve the ECM. The ECM is engaged in organogenesis via proliferation and traction forces that are mediated for instance by the cytoskeleton (reviewed in Ref. 107 ). The ECM plays an important role in controlling experience-dependent plasticity in the visual cortex during ontogenesis, acting on spine dynamics and axonal sprouting (as reviewed here Ref. 98 ). This is confirmed by our results of significant differential expression in the ECM. Together with previous literature, our findings suggest a connection between the vascular system and extracellular matrix in establishing asymmetries. Overall, the involvement of blood vessels/vascular system and extracellular matrix in asymmetries support the hypothesis that genes associated with lateralisation might be involved in anatomic structure development guided by environmental cues, rather than in asymmetrical body formation per se (see also Ref. 38 ).
The ECM is connected to ion activity and to the asymmetric movement of primary cilia, that drives asymmetries in mice and other vertebrate species 23 . Moreover, cilia are present in the retina at the level of photoreceptors and in other retinal cell types 28,29 , acting in the detection and transduction of visual signals, macromolecules and other signals 108 . At the molecular level, we found enrichment of ion activity, including Ca 2+ ion binding, that is connected to the activity of cilia. Previous reports have highlighted the role of calcium ions in lateralisation (e.g., fish 109 , mice 31 , chicks 110 ). The membrane of primary cilia includes receptors for ECM proteins 26,91 , and extensive signalling has been established between cilia and ECM, with defects in cilia associated with developmental defects 24,26 . Enrichment for ECM-receptor interaction and focal adhesion has been found in the functional genetics of handedness 38 , as well as enrichment in the TGF-beta signalling pathway. This pathway is connected to PDGFRB (platelet-derived growth factor receptor β), one of the differentially expressed genes in the 6-h light exposure left-right comparison. This gene, essential for the development of the cardiovascular system, helps in the rearrangement of the actin cytoskeleton and is involved in the primary cilia mechanism 92 . PDGFR β receptors are primarily located and activated outside the cilium, inducing disassembly of the cilium during morphogenesis. By controlling the assembly and disassembly of their primary cilium, cells can modulate the response to extracellular morphogenetic signals, and determine their cellular architecture (see dynamics in the corneal endothelium of mice 111 ). Overall, our results suggest that there are similarities between the endogenous lateralisation mechanisms and environmentally driven lateralisation since in both cases there is an involvement of the extracellular matrix and the platelet-derived growth factor receptor PDGFRβ, that are usually connected to primary cilia. www.nature.com/scientificreports/ Among the differentially expressed genes in the right vs left retina there is PAX2, a transcription factor that plays an important role in the glia during retinal development, as previously shown in the chicken embryo 86 and in the development of the visual system in vertebrates 87 . PAX2 is involved in the generation of the optic stalk and in the first stage of chiasm asymmetric axonal migration 87,112 . PAX2 expression has also been reported in the chick retina at later developmental stages in Müller glial cells, astrocyte-like radial cells that can undergo reprogramming and repair in teleost fish and regenerate retinal neurons through a process of dedifferentiation and acquisition of multipotent progenitor identity 113 (see also Ref. 114 for a review). Our findings highlight genes that mediate the effect of environmental experience (light exposure) in the ontogeny of left-right asymmetries. The role of environmental factors in lateralisation is becoming progressively clearer 4,11,17 . Moreover, additive genetic factors explain only one quarter of the variance in human handedness, and it has been suggested that nonshared environmental factors can explain part of the remaining variance 20,115 (but see Refs. [116][117][118]. Hence, to understand the ontogeny of functional and brain lateralization we must elucidate the interaction between endogenous and environmental factors 2,3,11,14,17,66,67,[69][70][71] . It must also be considered how noncoding microRNA 72 , or epigenetic mechanisms, which affect gene activity and expression by modifying DNA accessibility or chromatin structure, mediate long-term effects of gene-environment interactions 42,71 .
Although retinal asymmetries have not been investigated in chickens yet, left-right differences in the distribution of retinal cells have been reported in other birds, such as starlings and cockatoos 119,120 . Based on the asymmetrical gene expression in the retina that we have documented, future studies should investigate the role of asymmetrical gene expression in the symmetrical/asymmetrical differentiation of retinal tissues.
To summarise, our data show a molecular association between endogenous and experience-driven lateralisation mediated by light exposure. We identified functional gene networks connected to the vascular system and extracellular matrix as mediators of light-induced asymmetries in the retina, and their connections. The evidence of a temporally defined window of these effects calls for further studies on modulating lateralisation via the molecular pathways that open and close sensitive periods in the developing brain, such as thyroid hormones 66,74,121 , with important biomedical implications regarding the mechanisms of morphogenesis, ontogenetic plasticity-such as opening/closing of sensitive periods, and for regulating left-right asymmetries.

Data availability
All of the sequencing data used are available in the NCBI Short Read Archive under the Bioproject PRJNA743180. All code used in this analysis is available online at https:// github. com/ peter thorp e5/ Gallu sggal6_ RNAseq.