Boundary sequences flanking the mouse tyrosinase locus ensure faithful pattern of gene expression

Control of gene expression is dictated by cell-type specific regulatory sequences that physically organize the structure of chromatin, including promoters, enhancers and insulators. While promoters and enhancers convey cell-type specific activating signals, insulators prevent the cross-talk of regulatory elements within adjacent loci and safeguard the specificity of action of promoters and enhancers towards their targets in a tissue specific manner. Using the mouse tyrosinase (Tyr) locus as an experimental model, a gene whose mutations are associated with albinism, we described the chromatin structure in cells at two distinct transcriptional states. Guided by chromatin structure, through the use of Chromosome Conformation Capture (3C), we identified sequences at the 5′ and 3′ boundaries of this mammalian gene that function as enhancers and insulators. By CRISPR/Cas9-mediated chromosomal deletion, we dissected the functions of these two regulatory elements in vivo in the mouse, at the endogenous chromosomal context, and proved their mechanistic role as genomic insulators, shielding the Tyr locus from the expression patterns of adjacent genes.

Extensive profiling of genome organization and chromatin folding suggested that gene expression and chromatin topology are highly correlated 1 . Actively transcribed genes reside in large A compartments, mega-base structures characterized by the presence of active chromatin marks and high reciprocal physical interactions 2 . Conversely, inactive genetic material is confined to the B compartment. Genes can switch their respective compartment [3][4][5] , when they become active during development or differentiation, illustrating the dynamic correlation between gene expression and chromatin topology. At a smaller scale (tens to hundreds of kilobases), highly interacting chromatin regions are organized in topologically associating domains (TADs) characterized by high interdomain interactions 6 . TADs partition chromatin in independent units, delimited by boundaries enriched in CTCFoccupied regions. Disruption of boundaries results in misfolding of TADs, with spurious enhancer-promoter interactions that often lead to disease or developmental defects [7][8][9] . Within TADs, enhancers reach their target promoters through intra-TAD chromatin loops, mediated by CTCF, Cohesin, Mediator 10 and other transcription factors 11 .
The tyrosinase (TYR ) gene encodes the rate-limiting enzyme in the biosynthesis of the melanin. Mutations in TYR are the molecular cause of the most common type of albinism in Western countries (oculocutaneous albinism type 1, OCA1) a rare disease characterized by severe visual deficits and hypopigmentation in the skin, hair and eyes 12 . In the mouse, mutations at Tyr affect coat and eye pigmentation, with similar alterations in the visual pathway 13 , but without gross detrimental effects on overall physiology and viability. For this reason, Tyr has been targeted in mouse 14 and other organisms [15][16][17] to benchmark transgenesis and genome engineering strategies, from oocyte microinjection of DNA 18,19 to the most recent CRISPR nuclease 20 and base editor approaches 21,22 . Tyr is expressed in neural crest-derived melanocytes populating the skin, the hair bulbs and numerous other body locations, including the choroid, a pigmented layer found between the retina and the sclera of the eye. In addition to melanocytes, Tyr is also expressed is the retinal pigmented epithelium (RPE), a monolayer of hexagonal and binucleated pigmented cells located at the outer most layer of the retina, derived from the optic cup 23,24 . Its highly specific expression pattern makes the Tyr locus an excellent experimental model locus to explore the mechanisms of gene expression, its relation with chromatin topology and the consequences of genetic perturbations at structurally-relevant sequences. The Tyr promoter contains the minimal information necessary to correctly drive Tyr expression in vivo 18,25,26 . This transcriptional specificity is conveyed by two highly conserved motifs, the Tyr Initiator (Irn) and the M-box, both occupied by the transcription factor Mitf 27 . A regulatory element 5′ upstream in the murine Tyr locus was first identified as a DNAse I hypersensitivity site 28 . The inclusion of this element in a tyrosinase minigene resulted in copy number-dependant gene expression in transgenic mice 29 , suggesting features of a chromatin boundary 30 or a locus control region (LCR) 31,32 . Furthermore, when microinjected into oocytes from albino mice, large yeast artificial chromosomes (YACs) encompassing these Tyr upstream sequences fully rescued coat-colour pigmentation 32 , whereas shorter constructs, lacking far upstream sequences, provided weak and variegated Tyr expression 31,33 . More recently, we deleted the Tyr 5′ element at the endogenous locus using a CRISPR/Cas9 approach and generated allelic series to identify a 410 bp core enhancer sequence 34 .
Here, using in vivo mouse studies, we resolve the chromatin structure of the Tyr locus and integrate it with functional data obtained by targeting additional putative cis-regulatory sequences at the 3′ downstream of the Tyr gene.

Results
the Tyr locus is organized into two distinct chromatin loops. The mouse Tyr locus spans over 180 kb on chromosome 7 and its expression patter is restricted to melanocytes and in cells of the retinal pigmented epithelium 24 . Conversely, genes flanking the Tyr locus manifest completely different expression patters as evidenced by the ubiquitous expression of Nox4 35 and the central nervous system-specific expression of Grm5 36 gene. We used published data [37][38][39][40] to explore the chromatin landscape of the Tyr locus. In melanocytes, but not in fibroblasts, acetylation at lysine 27 of histone H3 (H3K27ac) highlights the gene promoter and regulatory elements upstream of the gene body. Upstream and downstream Tyr, chromatin carries repressive chromatin marks such as tri-methylation at lysine 27 of histone H3 (H3K27me3), covering the genes Nox4 and Grm5, inactive in melanocytes (Fig. 1a). To explore the chromatin structure of the Tyr locus in more details and relate it to known regulatory elements of the locus, we decided to analyze this chromosomal territory by performing chromosome conformation capture (3C) in both Tyr expressing cells (B16 melanoma cells) and non-expressing cells (NIH 3T3 L929 fibroblasts). We selected a 367 bp DpnII fragment containing the Irn and M-box as anchor fragment and interrogated 3C libraries. The Tyr promoter engages in multiple interactions within the Tyr locus. A strong interaction was observed within a fragment 80 kb 3′ downstream (thereafter, Tyr 3′ element) (Fig. 1b) and confirmed by Sanger sequencing (Fig. S1b). The promoter/Tyr 3′ interaction is observed in both B16 and L929 cells, indicating the constitutive nature of the chromatin loop, suggesting a structural role, independently of Tyr expression. No interactions are detected with control sequences immediately upstream and downstream the Tyr 3′ element. A cell-type specific interaction is found ~ 15 kb 5′ upstream of the Tyr gene, co-localizing with the Tyr 5′ element, analyzed previously 29,31,34 , restricted to the DpnII fragment containing the Tyr LCR A and B boxes 29,41 (Fig. 1b) and confirmed by Sanger sequencing (Fig. S1a). To confirm the Tyr 5′/Tyr promoter interaction, we performed the reciprocal analysis by using a bait primer (anchor) located at the Tyr 5′ sequence. Indeed, when this Tyr 5′ region was used as bait, we could recover signal only at the Tyr promoter confirming the specificity of the interaction (Fig. 1c). We previously showed that inactivation of this sequence in vivo by deletion results in a loss of pigmentation 34 , confirming the relevance of such interaction for gene expression. No interaction is observed between the Tyr 5′ element and far 5′ upstream elements, including the CNS-2, an element previously described to be involved in Tyr expression regulation 42,43 . No contacts are detected between the Tyr 5′ and Tyr 3′ elements, indicating that these two boundary elements probably compete with the same binding sites in the Tyr promoter and are not binding with each other (Fig. 1c). Finally, the anchor primer was positioned at the CNS-2 42 , a sequence observed to be dispensable for Tyr expression in transgenic mice generated with different YAC Tyr transgenes 31 . Using this fragment as bait we could not observe any interaction within the Tyr locus in both 3C libraries in the cell lines used (Fig. 1d).
the Tyr 5′ boundary element contains also a melanocyte-specific enhancer. The anatomical structure of the eye offers the unique opportunity to interrogate the specificity of Tyr regulatory elements with regard of the cell of origin. In fact, pigmented, neural-crest (NC) derived melanocytes in the choroid are surrounding the retinal pigmented epithelium (RPE), a single layer of cells at the basis of the neuroretina. We previously demonstrated that skin pigmentation relies on the Tyr 5′ element 34 . To address the relevance of Tyr 5′ in the eye, we compared whole-mount retinae and histological eye sections of wild-type pigmented and Tyr 5′ mice. Whole mount retina from Tyr 5′ mice display a choroid-specific loss of pigmentation. The typical distribution of hexagonal RPE cells, with unpigmented nuclei, becomes evident in the retina from Tyr 5′ mice, whereas is largely masked by underlaying pigmented choroid in the retina from wild-type animals (Fig. 2a). Histological analysis confirmed that loss of pigmentation is restricted to the melanocytes present in the choroid, whereas pigmentation of the RPE appears to be unaffected, suggesting that the Tyr 5′ element activity is specific to NC-derived cells and dispensable in the RPE (Fig. 2b). The core sequence 34 of the Tyr 5′ enhancer is enriched in Sox10 and Mitf binding motifs (Fig. S2a), and it is indeed occupied by Sox10 (Fig. 2c), a NC transcription factor that can activate Tyr 44 in melanocytes, but not in RPE, where it is not expressed 45 . These results indicates that the Tyr 5′ regulatory element, which had been shown to operate as a genomic boundary 30,34,46 , and be bound by Usf1 factor 42 , which has been associated with insulator activity 47 , also contains and functions as a transcriptional enhancer 29 to drive Tyr expression in neural crest-derived melanocytes in the skin and the choroid, but not in the RPE, in line with previous research that hypothesized the presence of additional RPE-specific regulatory elements 42,43 .

Scientific RepoRtS
| (2020) 10:15494 | https://doi.org/10.1038/s41598-020-72543-0 www.nature.com/scientificreports/ www.nature.com/scientificreports/ the Tyr locus is flanked by chromatin insulators. We previously described and dissected the Tyr 5′ element, that coincides with a peak in our 3C assay and whose inactivation in vivo results in cell-type specific loss of pigmentation 34 . Experiments in transgenic mice and in drosophila 30 and in zebrafish embryos 46 suggested that the Tyr 5′ element has insulator activity, in line with previous studies reporting binding of the boundaryassociated Usf1 factor 47 at the Tyr 5′ element 42 . Analysis of the nucleotide Tyr 5′ sequence highlighted a cluster of transcription factor binding motifs including Sox10, Mitf and the boundary-associated Usf-1 (Fig. S2a). Indeed, both the full-length Tyr 5′ element and its synthetic GAB core display enhancer-blocking activity in human cells ( Fig. 3c). At this point, we reasoned that the sequence corresponding to strongest chromatin interaction emerged in our 3C assay-located 3′ upstream the Tyr gene-could also act as a chromatin boundary as well. We hypothesized that two boundary elements, at each end of the locus, would guarantee the strict cell-type specificity of Tyr expression and shield the locus from undesired interactions with neighbouring Nox4 and Grm5 genes (Fig. 3a). Supporting our hypothesis, we observed that the Tyr 3′ region is occupied by CTCF, a protein that exert insulator activity, in several mouse tissues and cell lines (Fig. S3a). Furthermore, a SINEB2 retrotransposon, a class of repeated mobile elements associated with insulator activity 48 , was also found in this Tyr 3′ region (Fig. 3a, Fig. S3b), carrying a CTCF-binding motif (Fig. S3c). We therefore interrogated the Tyr 3′ element for its enhancer-blocking activity using an enhancer-blocking assay (EBA). Indeed, the Tyr 3′ element functions as an insulator in vitro: the full-length, 2.5 kb Tyr 3′ element displayed the highest enhancer-blocking activity. Its core component, a 241 bp element containing the SINEB2 retrotransposon including a CTCF-binding site, retains more than 50% of the activity of the full-length counterpart in less than one tenth of the size. Upon mutating the CTCF binding site we could reveal the contribution of the SINEB2 element to the boundary activity, which   www.nature.com/scientificreports/ was reduced, but not fully abolished (Fig. 3c). This was expected, in the absence of the consensus CTCF binding motif 49 .
In vivo, using an EBA in zebrafish embryos, where we tested the ability of a sequence to interrupt promoterenhancer interactions, we confirmed the statistically significant insulator property of the Tyr 3′ core element (Fig. 3d), that is capable of blocking enhancers in vivo. Likewise, upon introducing mutations at the CTCF site (Fig. S3d), the resulting enhancer-blocking activity is reduced ( Fig. 3e-g) thereby proving that the boundary activity is, at least partially, CTCF-dependent. However, the inactivation of the CTCF-binding motif does not result in a statistically significant decrease in insulator activity. We speculate that the residual activity depends on the SINEB2 sequence, and the residual CTCF binding that may occur associated to the partially mutagenized CTCF motif. These results suggested that the Tyr expression domain is flanked by two elements capable of blocking enhancer activity in vitro and in vivo, by two chromosomal insulators.
the Tyr 3′ element is dispensable for Tyr expression. In order to functionally assess the role, in vivo, at the endogenous location, of the Tyr 3′ element, we used the CRISPR/Cas9 system to produce targeted chromosomal deletions involving this Tyr 3′ site. Guide RNAs were picked flanking a 2.8 kb region encompassing the Tyr 3′ insulator, including the DpnII fragment highlighted by the 3C assay and the sequences tested in the enhancer-blocking assays. We established two mouse lines carrying deletions of different sizes in homozygosis (TYRINS3#16 and TYRINS3#26; Fig. 4a,b). Coat colour of TYRINS3 mice appeared indistinguishable from that of wild type pigmented mice (Fig. 4c). The melanin contents from skin and eye extracts from both lines were comparable to those of wild-type animals (Fig. 4d,e). Whole-mount retinae of Tyr 3′ animals were also indistinguishable from those of wild-type littermates (Fig. 4f,g), and we did not observe delay in the appearance of eye pigmentation in developing embryos (Fig. S4a,b) suggesting that perhaps the Tyr 3′ element might not be required for Tyr expression in melanocytes or RPE. Insula-tors have been classically studied by means of reporter assays or by testing for copy-number dependent activity of integrated reporter genes 50 . Such approaches limited our understanding of the actual impact of insulators on local gene expression and that of neighboring genes. Only recently, with the advent of genome engineering techniques 51,52 , the study of boundary elements at their endogenous context became accessible. This can be done by direct editing of CTCF-binding sites 53 or by introducing large rearrangements that disrupt topological chromatin domains 7 . In this context, our mouse lines with genome-edited insulator sequences represent a www.nature.com/scientificreports/ unique opportunity to investigate the functional role of these two genomic insulators and their impact on the expression of Tyr and the adjacent Nox4 and Grm5 genes at their endogenous locus. For the purpose of gaining a mechanistic insight in the functional organization of the locus, we generated homozygous animal lacking the Tyr 3′ element (this work, TYRINS3 mice) or the Tyr 5′ 34 (TYRINS5 mice), and quantified gene expression levels of Tyr and of the two flanking genes Nox4 and Grm5 in different tissues.

Inactivation of the Tyr boundaries perturbs the transcription of flanking genes in vivo.
In the eye, loss of the Tyr 5′ element in TYRINS5 34 homozygous mice results in 50% reduction of Tyr expression. This finding is in line with previous histological analyses 34 and with its enhancer activity in melanocytes of the choroid (Fig. 2). Deletion of the Tyr 3′ element in TYRINS3 homozygous mice does not affect Tyr expression in the eye. Once again, this finding is in agreement with eye melanin measurements (Fig. 4e) and whole mount retina, that are indistinguishable from wild-type mice (Fig. 4f,g). The expression of Nox4 and Grm5 remains unaltered in the eye of TYRINS3 and TYRINS5 homozygous mutant animals (Fig. 5a), respectively.
Loss of Tyr 5′ in TYRINS5 homozygous animals almost abolished Tyr expression in the skin, in agreement with their lighter coat colour and previous melanin measurements 34 . Interestingly, deletion of Tyr 3′, the element which separates Tyr and Nox4, leads to perturbation in Nox4 expression as evidenced by significant upregulation of Nox4 transcripts in the skin of TYRINS3 homozygous mice (Fig. 5b). This perturbation in the pattern and strength of Nox4 expression in the absence of the Tyr 3′ element, further confirms its function as insulator (Fig. 3B-D).
As previously reported 24 , Tyr is not expressed in the the brain of wild-type animals and this pattern does not change upon deletion of the Tyr 3′ or Tyr 5′ elements. Conversely to Tyr, Grm5 is a neuro-specific gene and is expressed in the brain of wild type mice and both Tyr 5′ and Tyr 3′ mutants. Once again, our findings clearly show that the loss of the boundary region can perturb expression of genes flanking the boundary. This is evidenced by the significant increase in the brain expression levels of Grm5 upon deleting the Tyr 5′ element that separates the Tyr and Grm5 loci, in TYRINS5 mice. Consistent with our model, deletion of the Tyr 3′ element separating Tyr and Nox4 genes in TYRINS3 mice has no bearing on Grm5 transcript levels (Fig. 5c).

A model of the Tyr locus organization.
In this study we have been able to compare, for the first time, the phenotype of various genetically modified mice produced over 20 years to understand the function of the murine Tyr locus. These mice include the classical YAC Tyr transgenic mice 31,32 , encompassing the entire Tyr expression domain; the parental NMRI albino mice, carrying the known C103S albino mutation; and the two CRISPR-mediated deletions of the Tyr 5′ element (TYRINS5 mice 34 ) and the Tyr 3′ element (TYRINS3 mice), described here (Fig. 6a,b). The YAC Tyr transgenic mice are undistinguishable from wild-type pigmented agouti mice and from TYRINS3 mice. On the contrary, TYRINS5 mice 34 display a faint gray coat colour indicating a significantly reduced Tyr expression in the skin.
Considering the mouse phenotypes observed, the 3C biochemical analyses conducted and the presence of binding sites for known nuclear factors in these two key regulatory elements at either side of the murine Tyr locus, we propose a model that could explain the altered behaviour of Tyr, Nox4 and Grm5 gene expression levels in a wild-type animal (Fig. 6c), in TYRINS3 mice (Fig. 6d) and in TYRINS5 mice 34 (Fig. 6e).

Discussion
In this work we describe the chromatin structure of the mouse Tyr locus by 3C in mouse melanoma (expressing Tyr) and fibroblast (non-expressing Tyr) cells. We describe mechanistically how the murine Tyr locus functions in these two types of cells. Through our biochemical and molecular studies, we have identified two chromatin loops: one involves the interaction between Tyr 5′ enhancer and the promoter, and correlates with Tyr expression; the other, detected in Tyr both type of cells, involves the structural protein CTCF, bound to a 3′ far upstream element.
The Tyr 5′ enhancer was initially identified in the classical "chinchilla mottled" Tyr mouse mutation (Tyr c-m ) 28 and, thereafter, its pivotal role confirmed in mouse transgenesis experiments, through the rescue of pigmentation using both small 29 and large DNA constructs 31 , or using reporter lines 44 in oocytes from albino mice. These experiments highlighted with high precision the location and the mechanism of the enhancer, that is transactivated by Sox10 and Mitf 44 . USF1 transcription factor, usually found in chromatin insulators 47 , was also reported to interact with this Tyr 5′ element 42 , likely explaining the boundary activities described in mice 30,34 , flies 30 and zebrafish 46 . Due to confounding factors, including overall transgene integrity, copy number and integration site, it remains challenging to assess up to what extent the enhancer is required to drive Tyr expression using classical transgenic mouse models, since even the smallest Tyr promoter tested in transgenic mice (270 base pairs in length) was sufficient to drive Tyr expression in both melanocytes and RPE cells 25 .
More recently, thanks to advances in genome editing techniques 52,54 , we targeted the Tyr 5′ enhancer and boundary elements at the endogenous locus (TYRINS5 mice 34 ). The phenotype of mice lacking the endogenous enhancer is less severe compared with what we had observed using engineered YACs 31,33 , where variegation was regularly reported, however highly reproducible across TYRINS5 mouse lines with similar deletions 34 . Using the TYRINS5 mice, we determined that the Tyr 5′ element is absolutely required to drive Tyr expression in neuralcrest derived melanocytes, where Tyr gene expression is almost undetectable upon deleting the Tyr 5′ element (Figs. 2b, 5b), but dispensable in the RPE (Fig. 2). Mechanistically, this correlates well with differences in gene expression between melanocytes, that express Sox10, and cells of the RPE, that are Sox10 negative 44,45 . In fact, in melanocytes, Sox10 binding overlaps the core Tyr 5′ enhancer 38 (Fig. 2c) that is enriched in Sox10 binding motifs (Fig. S2a). A RPE-specific enhancer was identified using BAC transgenesis and transient LacZ reporter mouse lines 42,43 ; however, the specificity and relevance of such distant enhancer (absent in some of the YAC Tyr transgenic mouse lines that were indistinguishable from wild-type pigmented mice, where DNA sequences upstream of the 5′ element had been deleted 31 ) remain to be validated at the native chromatin context, and in relation with the Tyr promoter. These are ongoing experiments whose results will be discussed in a following study.

Scientific RepoRtS
| (2020) 10:15494 | https://doi.org/10.1038/s41598-020-72543-0 www.nature.com/scientificreports/ At the 3′ end of the locus, at the intergenic space between Tyr and Nox4, we identified a sequence that interacts strongly with the Tyr promoter in melanoma cells and fibroblasts (Fig. 1b). This sequence contains a SINEB2 retrotransposon and it is bound by CTCF (Fig. S3). These features are typically associated with insulators 48,55 , a www.nature.com/scientificreports/ class of regulatory elements able to block enhancer-promoter interactions and to partition chromatin in transcriptionally independent domains 56 . Indeed, we could detect a strong insulator activity, in mammalian cells and in zebrafish embryos at this Tyr 3′ element comparable to that of the archetypic chicken cHS4 insulator, both CTCF-and SINEB2-dependant (Fig. 3). This further highlights the role of retrotransposons in the gene insulation 48,57 , in control of chromatin topology 58 , and in genome evolution 59,60 . Deletion of this sequence in the mouse does not affect pigmentation or melanin levels in the skin or the eye of two independent mouse lines (Fig. 4); Tyr mRNA level in the eye and skin does not depend of the Tyr 3′ insulator (Fig. 5). However, when we measured the expression level of Nox4, located further upstream and beyond the Tyr 3′ boundary, we found that Nox4 was overexpressed in the skin, suggesting that loss of the Tyr 3′ allowed Tyr regulatory elements to abnormally and ectopically transactivate Nox4. Similarly, we profiled gene expression in tissues of mice lacking the Tyr 5′ element. Grm5, the gene found 5′ upstream of Tyr and beyond the Tyr 5′, is overexpressed in the brain of animals lacking the Tyr 5′. We speculate that rearrangements in chromatin topology in the absence of the Tyr 5′ element expose the Grm5 promoter to distal regulatory sequences. Several mouse models (Fig. 6a,b) illustrated the relevance of regulatory regions at the mouse Tyr locus. Combining profiling of chromatin topology and functional study in the mouse, we devised a model for the locus: in melanocytes, the Tyr promoter engages two regulatory elements occupied by Sox10 and CTCF (Fig. 6c). Disruption of the Tyr 3′ element, that normally interacts with the Tyr promoter and exert insulator activity, results in overexpression of the flanking gene Nox4, in the skin (Fig. 6d). Perturbation of the Tyr 5′ element, that normally transactivate the Tyr promoter and functions as a boundary, results in loss of Tyr expression in gain in expression of the flanking Grm5 (Fig. 6e).
These data illustrate the role of enhancers and chromatin insulators in gene expression at the scale of a given locus and its adjacent genes. They also highlight how mutations at regulatory elements could alter the expression of genes in multiple tissues. For more than 20 years (1990 18 -2012) our ability to functionally and structurally interrogate the non-coding DNA sequences within the mouse Tyr locus was extremely limited. The extraordinarily abundance of DNA repetitive elements, including LINE1 retrotransposons in the vicinity of Tyr 30 , rendered www.nature.com/scientificreports/ impossible any attempt to apply standard homologous recombination techniques in mouse embryonic stem (ES) cells. With the advent of the CRISPR/Cas9 techniques applied to mouse functional genomics 61 , it has become possible dissecting the mechanistic role of selected regulatory elements at their endogenous chromosomal location. This study highlights the relevance of investigating the non-coding genome in their native chromosomal context, and the advantage of using mouse models over cellular systems in the study of regulatory sequences.
Briefly, 1 × 10 7 B16F0 mouse melanoma or NIH3T3-L929 mouse fibroblast cells were fixed in 2% PFA/PBS. After quenching, cells were resuspended in 500 μl lysis buffer ( Supplementary Table S1. Luminescence was measured for each construct in triplicate and used to quantify enhancer-blocking activity. The cHS4 chicken beta-globin insulator and its minimal core motif (II/III) were used as positive control; a mutated version of the II/III sequence was used as negative control. One-way ANOVA with Bonferroni post-hoc correction test was used.
Putative insulator sequences were cloned in between a hindbrain-specific enhancer and a somite-specific promoter. Single-copy integration in the zebrafish genome was obtained by Tol2-mediated transgenesis in fertilized zebrafish embryos. GFP fluorescence was acquired at 36 hpf and quantified using Laser Pix (BioRad). Enhancer-blocking activity was measured based on relative levels of somite/CNS GFP expression. Statistical analyses (median test) were calculated with IBM-SPSS v.21.
CRISPR/Cas9-mediated chromosomal deletion. A pair of RNA guides flanking the Tyr 3′ downstream element was designed using the crispr.mit.edu online tool and cloned by Golden Gate Cloning into the Esp3I sites of MLM3636 plasmid (Addgene #43860). Cas9 mRNA was prepared from hCas9 plasmid (Addgene #41815) as described 34,64 and injected into B6D2F2 (Harlan; originated from B6CBAF1/OlaHsd) mouse fertilized eggs. Founder animals were bred to albino outbred HsdWin:NMRI (Harlan) mice, the reference albino genetic background used in all previous Tyr transgenic mouse studies 29,31,33,34 . Tyr 3′ mouse lines were maintained in albino outbred HsdWin:NMRI background. Primers used for cloning and genotyping are available in Supplementary Table S1. Whole-mount retinae and Hematoxylin-Eosin stained cortical sections were prepared as described 34 . Melanin content of skin and eye biopsies was estimated by optical density 34 .
In this study, all genome-edited mouse lines were named as TYRINS3#, followed by a number corresponding to the founder they were derived from, as done for the TYRINS5# lines described before 34 . The recommended nomenclature for these newly generated mouse lines is stock-Tyr emXLmon , where 'X' is the corresponding ordinal number for each genome-edited (endonuclease mediated) mouse line 65 . For brevity, TYRINS5# and TYRINS3# mouse gene-edited lines are also referred to as Tyr 5′ and Tyr 3′, respectively, throughout this manuscript.
All animal procedures with mice reported in this work were first validated by the National Centre for Biotechnology Ethics Committee on Animal Experimentation (OEBA), thereafter approved by the National CSIC Ethics Committee and eventually authorized by the Autonomous Community of Madrid, acting as the competent authority, according to the Spanish legislation and the European Directive 2010/63/EU. All mice were housed at the registered CNB animal facility, fed and provided water and regular rodent chow ad libitum, with a light/dark cycle 08:00-20:00, according to the European and Spanish norms, and the animal welfare recommendations. Both male and female individuals were used indistinctly. mRNA expression analysis. Total RNA was isolated from cultured cells and tissues using the RNeasy Mini Kit (Qiagen) with on-column DNAse treatment. 500 ng of RNA were retrotranscribed using Super-Script III Reverse Transcriptase (Thermo Fisher). qPCR reactions were performed with TaqMan Universal PCR Master Mix (Thermo Fisher) using probes for Tyr (Mm00495817_m1), Nox4 (Mm00479246_m1), Grm5 (Mm00690332_m1) and Tbp (Mm00446973_m1) as described 66 .
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.