The genome and transcriptome of Japanese flounder provide insights into flatfish asymmetry

Songlin Chen, Manfred Schartl, Qingyin Wang, Deborah M. Power and colleagues analyze the genome of the Japanese flounder and its transcriptome dynamics during metamorphosis. They report a role for thyroid hormone and retinoic acid signaling, as well as phototransduction pathways, in the regulation of craniofacial asymmetry. Flatfish have the most extreme asymmetric body morphology of vertebrates. During metamorphosis, one eye migrates to the contralateral side of the skull, and this migration is accompanied by extensive craniofacial transformations and simultaneous development of lopsided body pigmentation1,2,3,4,5. The evolution of this developmental and physiological innovation remains enigmatic. Comparative genomics of two flatfish and transcriptomic analyses during metamorphosis point to a role for thyroid hormone and retinoic acid signaling, as well as phototransduction pathways. We demonstrate that retinoic acid is critical in establishing asymmetric pigmentation and, via cross-talk with thyroid hormones, in modulating eye migration. The unexpected expression of the visual opsins from the phototransduction pathway in the skin translates illumination differences and generates retinoic acid gradients that underlie the generation of asymmetry. Identifying the genetic underpinning of this unique developmental process answers long-standing questions about the evolutionary origin of asymmetry, but it also provides insight into the mechanisms that control body shape in vertebrates.

). This is further emphasized by the link with the expansion of the homeodomain-interacting protein kinase (Hipk) gene family, which is involved in p53-mediated cellular apoptosis and regulation of eye size, lens formation, and retinal lamination during late embryogenesis 9 (Fig. 1c and Supplementary Note), as well as in the dramatic changes in the eyes during metamorphosis 10 .
We reasoned that comparative transcriptomics of well-defined developmental stages of Japanese flounder and Chinese tongue sole ( Fig. 1a and Supplementary Table 5), cross-referenced with genes identified in flatfish but absent from other teleosts or PSGs, should yield the master genes, namely, genes responsive to THs that occupy a central position in the gene networks underpinning the tissue transformations in flatfish metamorphosis. We identified 2,307 differentially expressed orthologous genes (DEGs) that were common in both species when DEGs identified by pairwise comparison of each metamorphic stage were analyzed (Supplementary Table 20). This large gene set was enriched for several pathways that can be linked to recognized morphological or physiological changes during metamorphosis (for example, TH action, photo-and olfactory transduction, and nutrition and metabolism), and it contains many genes that are involved in neural development and remodeling, anterior-posterior patterning, cell proliferation, and photoreceptor development (Supplementary A role for THs in flatfish metamorphosis was identified in the early 1980s 11 . THs are obligatory for vertebrate metamorphosis 12 and function by binding to nuclear-localized receptors (Thra and Thrb), or they may have a nongenomic cellular role and activate a complex hierarchical cascade of target genes [13][14][15] (Fig. 3a). The essential role of THs in flatfish metamorphosis 16 was reiterated in the current study by the detection of 21 differentially expressed orthologs of the  Figure 1 Flatfish metamorphosis and phylogenomic evolution. (a) Gross morphology of metamorphosis in Japanese flounder. DAH, days after hatching. (b) Maximum-likelihood tree from eight teleost species. The number of gene-family expansion events is indicated in turquoise, and that of gene-family contraction events in pink. (c) Flatfish-specific expansion of the hipk gene families. Phylogenetic analysis of deduced Hipk proteins from teleost fish generated two major clusters. One cluster (left side) comprised the three subgroups of Hipk1, Hipk2, and Hipk3, which are distinct from one another. Gene numbers in all analyzed species were very similar in the subgroups. The second cluster (right side) of Hipk proteins was enriched for flatfishlineage-specific expansions, in particular one expanded Hipk2-like (Hipk2*) clade in Chinese tongue sole and one expanded Hipk1-like (Hipk1*) clade in both Chinese tongue sole and Japanese flounder. l e t t e r s TH signaling pathway in Chinese tongue sole and Japanese flounder (Supplementary Table 28). Distribution or expression of thyroxine (T4), tri-iodothyronine (T3), deiodinase 2 (dio2), dio3, thraa, and thrb1 colocalized around the eyes during Japanese flounder metamorphosis, and this pattern became progressively asymmetric with the advance toward climax (Supplementary Figs. 11-18). Treatment with methimazole (MMI), an inhibitor of TH synthesis that blocks flatfish metamorphosis 4 , inhibited eye migration (Supplementary Fig. 19 and Supplementary Table 29) and reduced the asymmetric expression of dio3, thraa, and thrb1 (Supplementary Figs. 20 and 21).
RAs are ligands of the nuclear receptors Rars (the ligand is alltrans RA, or ATRA) and Rxrs (the ligand is 9-cis-RA), which form heterodimers with TH receptors 17 and are well-recognized morphogens involved in lateralization 18 . The strong positive selection of genes involved in the RA system in Japanese flounder (aldh1l2 and rbp4b) and Chinese tongue sole (aldh1l1 and rdh12) (Supplementary Tables 17 and 18) and the enrichment of genes of the RA pathway in the DEGs during metamorphosis (rxra, rxrg, rarg, aldh1l2, rbp2, raet1, rdh12, rbp3, and lrat) (Supplementary Table 30) suggest their involvement in asymmetry formation. We hypothesized that DEGs of the RA signaling pathway might explain the acquisition of TH signaling pathway asymmetry during metamorphosis, and that this cross-talk could be the basis for the generation of asymmetry in flatfish. Genes encoding isoforms of RA receptors (Rxra and Rxrb) that heterodimerize in yeast two-hybrid assays (Supplementary Figs. 22 and 23) with Thraa and Thrb1, respectively, showed marked left-right asymmetric gene expression in the epithelial tissue surrounding the eyes in the pro-metamorphic stage, when the eye starts to migrate (Supplementary Fig. 24 and Supplementary Note), thereby preceding the pronounced asymmetry that develops at the morphological level. Other genes encoding proteins involved in RA metabolism (such as cyp26, rbp2, and aldh1-3) showed variable degrees of asymmetric expression at the metamorphic climax (Supplementary Figs. 25  and 26). MMI-the inhibitor of TH synthesis that reduced the asymmetric expression of dio3, thraa, and thrb1-had the same effect on expression of rxra and rara (Supplementary Fig. 27). In summary, genes of the TH and RA signaling pathways had coincident expression patterns in the region of the head most profoundly modified during metamorphosis. Injection of 9-cis-RA or ATRA into the suborbital area of the prospective blind side eye of pre-metamorphic Japanese l e t t e r s flounder larvae inhibited eye migration in 95% (chi-square tests, P = 1.28 × 10 −14 ) and 97% (chi-square tests, P = 2.2 × 10 −16 ) of the larvae, respectively. Injection of a noncompeting Rxra ligand, 13-cis-RA, inhibited eye migration in 42% of the larvae (versus 37% for the DMSO control) ( Fig. 3b and Supplementary Fig. 28). As expected, the proliferation of cells in the suborbital area of the blind side eye was notably reduced in 9-cis-RA-and ATRA-treated larvae as compared to that in the DMSO-treated controls ( Fig. 3c and Supplementary  Fig. 29). ATRA treatment led to upregulation of rara expression and downregulation of thraa and thrb1 expression, but it had no effect on the other retinoic acid receptors (Supplementary Figs. 30 and 31). Taken together, these results suggest that during metamorphosis, eye migration in Japanese flounder is directly regulated by the interaction of RA-and TH-induced signaling pathways.
Pigmentation also becomes asymmetric during flatfish metamorphosis, and an association exists between impaired eye migration and the failure to develop an asymmetric pigmentation pattern in skin 19 . We therefore tested the hypothesis that the left-right asymmetric RA gradients important for eye migration might also establish asymmetric pigmentation 20 . We did this by determining the concentration of RA in the pigmented (ocular) and the unpigmented (blind) skin and by comparing the two skin transcriptomes (ocular and blind) of Japanese flounder and Chinese tongue sole (Supplementary Note). The concentrations of ATRA and 9-cis-RA in the ocular side skin  Figure 3 Molecular signatures in flatfish eye migration. (a) Putative cross-talk between TH signaling and RA signaling in eye migration. The cross-talk between RA and T3 in flatfish metamorphosis is hypothesized to consist of competition between the liganded Thrb and Rar for heterodimerization with the liganded Rxr. The heterodimers bind to TRE and RARE, respectively, to regulate gene transcription 29 . 9cRA, 9-cis-RA; 13cRA, 13-cis-RA.
(b) Inhibition of eye migration at 33 DAH in Japanese flounder at different metamorphic stages after injection of different RA isoforms. DMSO was used as the control. 9-cis-RA and ATRA, which have the highest affinity for Rxrs and Rars, blocked eye migration most effectively, whereas 13-cis-RA, which has low affinity for the receptors, was less effective. P values, as determined by chi-square test, indicate significant differences between the blind and ocular sides in the different treatments relative to the control. l e t t e r s (0.0396 ng/mg and 0.1430 ng/mg, respectively) were substantially (~10×) and significantly higher than in the blind side skin (0.0036 ng/ mg ATRA, one-way analysis of variance (ANOVA), P = 0.0294; 0.0172 ng/mg 9-cis-RA , one-way ANOVA, P = 0.0452) ( Supplementary  Fig. 32). Genes involved in pigment cell development and differentiation (such as mlana, tyrp1, and dct, among others) ( Supplementary  Fig. 33) and in the pigmentation pathway (Supplementary Table 31) were upregulated in the pigmented skin. In contrast, genes involved in the RA pathway (rara, rxrb, and aldh1) and TH signaling (thrb1 and dio2) had similar abundance in pigmented and unpigmented integument (Supplementary Fig. 34), although in situ hybridization analyses indicated that rxra and aldh3 expression were asymmetric ( Fig. 4a and Supplementary Figs. 35-38). Injection of 9-cis-RA and ATRA into the dermis of the blind side of pre-metamorphic larvae (Supplementary Fig. 39) increased adult chromatophores in postmetamorphic juveniles, whereas injection of retinol or DMSO had only a minor effect (Supplementary Figs. 40 and 41). In summary, we demonstrate that an RA gradient may be crucial for the generation of asymmetric pigmentation in post-metamorphic flatfish.
Recently, opsins in the skin of the adult octopus have been related to eye-independent color changes 21 . During flatfish metamorphosis there was ubiquitous expression in skin of genes of the phototransduction pathway (rh1, rh2, lws, sws1, and sws2) ( Fig. 4b and Supplementary  Figs. 42-46), although they were also in the skin of non-flatfish teleosts (Supplementary Figs. 47 and 48). Because opsins can regulate the synthesis of RA 22 , we hypothesized that the left-right RA gradient during asymmetric pigmentation and eye movement might result from increased light exposure on the ocular side due to the body tilt. This would then lead to an exclusive distribution of adult chromatophores on the ocular side at the latest stage of metamorphosis. To test this hypothesis, we exposed the blind and ocular sides of pre-metamorphic larvae to white light or light of different wavelengths. When both sides were exposed simultaneously to white light, the larval chromatophores distributed sporadically on both sides at 38 d after hatching (DAH) and were substantially increased in number in subsequent stages. Adult-type chromatophores started to proliferate at 50 DAH, achieving symmetrical pigmentation at 100 DAH (Supplementary Fig. 49) 23 . Light of mid-length and long wavelengths (red, yellow, and green) caused a few chromatophore spots to appear on the blind side, whereas light of short wavelengths (blue and violet) caused a large number of chromatophores to appear (Fig. 4c, Supplementary Fig. 50, and Supplementary Table 32). Adult-type chromatophores developed evenly on both sides, resembling the situation after injection of RA. This suggests that light-sensing opsins under blue-light illumination may stimulate the production of RA. Most notably, the migration of the eye to the ocular side did not occur under blue-light treatment ( Fig. 4c and Supplementary Fig. 50). The gene kit, which is related to pigment-cell development 24 , was co-expressed with the genes for blue-sensitive rho protein (sws2) and rxra in pigment cells (Fig. 4d). We did not observe co-expression with kit of rh1 or opn4, at the wavelengths that overlap with sws2, or of other RA-related genes. This suggests that Sws2 promotes the synthesis of RA, which in turn regulates Rxra during establishment of the left-right asymmetric pigmentation in flatfish. The involvement in this process of opsins in the eye remains an open question. In summary, the establishment of a high-quality genome of the Japanese flounder and its comparison to the genome of another flatfish, the Chinese tongue sole, led to the identification of genes and developmental pathways that regulate metamorphosis and establish body asymmetry. Notably, the results uncover a role for the lightsensing opsins in establishing an RA gradient in the skin during metamorphosis. This functions as a developmental trigger after the dorso-ventral axis skew and augments asymmetric pigmentation. The identified cooperation between TH and RA signaling molecules in the regulation of eye migration and skin pigmentation is reminiscent of the situation in insect metamorphosis whereby ecdysteroids (Ecd) and juvenile hormone (JH) have a critical joint role 25 . In common with THs, Ecd functions via a heterodimeric complex of two nuclear receptors, the ecdysone receptor and the JH receptor ultraspiracle (USP, the ortholog of Rxr) 26 , and 9-cis-RA is structurally similar to the terpenoid JH 27,28 . We propose that nuclear receptor transcription factors and terpenoids form the basis of the molecular mechanisms that trigger the unique metamorphosis found from invertebrates to vertebrates. Our detailed study of the molecular basis of flatfish metamorphosis reinforces the central role of endocrine factors in this late developmental event and a hitherto unexpected recapitulation of pivotal mechanisms regulating early development.

MeThoDs
Methods, including statements of data availability and any associated accession codes and references, are available in the online version of the paper.

oNLINe MeThoDs
All samples used in this study were collected according to the guidelines for the care and use of experimental animals established by the Yellow Sea Fisheries Research Institute, CAFS. We estimated the number of samples necessary from the number of replicates required for sequencing and experiments. We did not use a randomized method and did not use a blinded methodology for group allocation during the experiments. We used Fisher's exact test, t-tests, one-way ANOVA and chi-square tests 31 in the analysis of gene family evolution, positively selected genes, transcriptome data, K-mer distribution, GO and KEGG enrichment, RT-PCR, detection of RA concentration, and the effect of RA injection. When required, P values were accordingly adjusted for multiple testing after considering the Benjamini-Hochberg false discovery rate (FDR) 32 .
In the experiments involving immunohistochemistry and immunofluorescence, in situ RNA hybridization, in situ cell proliferation assay, irradiation treatment, yeast two-hybrid assay, and UPLC-MS/MS, we performed at least three independent or parallel experiments, and at least three individual larval specimens were treated in each group (at least three for positive and three for negative control) correspondingly.
Fish collection and maintenance. Gynogenetic Japanese flounder were temporarily maintained at the Key Laboratory of Sustainable Development of Marine Fisheries in Qingdao. The fish used for DNA extraction was anesthetized with MS-222 (Argent Chemical Laboratories, Redmond, WA, USA), and blood was collected and stored at 4 °C for later DNA isolation. The Japanese flounder and Chinese tongue sole larvae were produced from single batches of fertilized eggs that were maintained in single incubators at Huanghai Aquaculture Co. Ltd., Yantai, China. The larvae were fed live brine shrimp (Artemia nauplii) until the end of metamorphosis. The normal metamorphic stages were defined on the basis of eye migration: pre-metamorphosis (stage before eye migration; 16 d after hatching (DAH) in Japanese flounder, and 15 DAH in tongue sole); pro-metamorphosis (eye begins to migrate; 21 DAH in Japanese flounder, and 18 DAH in tongue sole); metamorphic climax (upper edge of the migrating eye beyond the dorsal margin; 25 DAH in Japanese flounder, and 19 DAH in tongue sole), and post-metamorphosis (entire eye migrates past the dorsal midline; 29 DAH in Japanese flounder, and 22 DAH in tongue sole) 33 . We collected skin from the right (blind) and left (ocular) sides of Japanese flounder and tongue sole as asymmetrical pigmentation samples when the fish were about 3 months old. 2-mercapto-1-methylimidazole (0.3-0.34 g/l; Sigma-Aldrich, St. Louis, MO, USA) was dissolved directly in water and administered for 4 d to inhibit endogenous TH synthesis by metamorphic Japanese flounder. Larvae for whole-mount in situ RNA hybridization, immunohistochemistry, and cell proliferation analyses were anesthetized with 0.1 µg/µl MS-222, fixed in 4% paraformaldehyde (PFA) overnight, and stored in methanol at 4 °C. Larvae for RNA and DNA isolation were preserved at −80 °C.
Genome sequencing and assembly. We applied a whole-genome shotgun strategy and next-generation sequencing technology on an Illumina HiSeq 2000 platform to sequence the Japanese flounder genome. DNA was extracted from one gynogenetic female Japanese flounder according to the traditional phenol-chloroform extraction protocol 8 . To decrease the risk of nonrandomness, we constructed eight paired-end sequencing libraries with insert sizes of 170 bp, 500 bp, 800 bp, 2 kb, 5 kb, 10 kb, 20 kb, and 40 kb and sequenced them (ten lanes). The Japanese flounder genome was assembled via the de novo method using SOAPdenovo 34 . We used a K-mer-based method to estimate genome size, as described previously 35 .
Anchoring the scaffolds with genetic maps. We constructed high-resolution genetic maps for Japanese flounder using simple-sequence repeats and singlenucleotide polymorphisms, respectively 6,7 . The scaffolds, with markers located on different chromosomes, were filtered.
Whole-genome synteny analysis among different teleost species. Wholegenome alignments of Japanese flounder and tongue sole, Japanese flounder and zebrafish, and tongue sole and zebrafish were carried out on their repeat-masked genomes using LASTZ (version 1.02.00) (T = 2, C = 2, H = 2,000, Y = 3,400, L = 5,500, and K = 2,000), and chainNetsoftware 36 was used to concatenate aligned fragments to longer fragments in a particular format, such as chainnet 36 .
Repeat annotation and comparison. We detected transposable elements (TEs) in the Japanese flounder genome by combining homology-based and ab initio-based approaches. We used RepeatMasker (version open-3.2.9) against the RepBase TE library (RepBase 16.02) for the homology-based method and executed RepeatProteinMask. We constructed de novo repeat libraries by running PILER (version 1.0), RepeatScout (version 1.0.5), and LTR-Finder 37 for the ab initio-based approach and filtered the results to obtain nonredundant consensus repeat libraries. Finally, RepeatMasker was run on the genomes using these de novo libraries as reference libraries. Tandem repeats were detected in the genomic sequence data using Tandem Repeats Finder software 38 .
Transcriptome sequencing and analysis. Total RNA was isolated and purified from all of the samples by a TRIzol extraction method 8 . RNA concentration was measured using Nanodrop technology. RNAs of three samples were pooled in equal quantities to construct the Illumina sequencing libraries, and then paired-end sequences (90 bp at each end) were generated using an Illumina HiSeq 2000.
The TopHat v1.2.0 package was used to map transcriptome reads to the genome and to identify DEGs. High-quality splice junctions were also predicted by TopHat. Gene expression was measured as reads per kilobase of gene per million mapped reads (RPKM) and adjusted using a scaling normalization method 39 . Only genes with RPKM > 1 in at least one sequenced sample were considered. DEGs were detected using DEseq 40 and Cuffdiff 41 . P values were adjusted by the FDR 32 . Only genes with an adjusted P < 0.05, obtained with any method, and fold change > 2 were considered true DEGs. Genes were annotated to GO categories on the basis of the orthologous relationship between the Japanese flounder gene set and the tongue sole gene set, which had perfect GO annotation. Fisher's exact and chi-square tests were used to identify whether a list of genes (foreground genes) was enriched in a specific GO category as compared to the background genes. The KEGG automatic annotation server annotated the genes to KEGG pathways, using zebrafish and human as references 42 . The K-means algorithm was used to produce groups of DEGs using the calculated means mode with Euclidean distance based on the log fold-change data 43 .

Structural prediction and annotation of genes.
We used a homology-based method, a de novo method, and an RNA-seq-based method to predict the genes in the Japanese flounder genome. Oreochromis niloticus, Cynoglossus semilaevis, Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes, Tetraodon nigroviridis, and Homo sapiens proteins were used for the homology-based prediction, and the coding sequences were mapped onto the genome using TblastN (1 × 10 −5 ). Then, homologous genome sequences were aligned against matching proteins using Genewise to define the gene models 44 . Augustus (version 2.5.5) and Genscan (version 1.0) were used to predict coding genes de novo in the repeat-masked genome. The RNA-seq data were mapped to the genome using TopHat (version 1.2.0), and transcriptomebased gene structures were obtained with cufflinks (http://cufflinks.cbcb.umd. edu/). We did some manual checking to obtain the final reference gene set. We annotated the motifs and domains for the reference Japanese flounder gene set using InterPro against publicly available databases, including Pfam, ProDom, SMART, PRINTS, SUPERFAMILY, and PROSITE. The GO annotation information was abstracted with Iprscan (version 4.7). The Japanese flounder protein sequences were mapped to KEGG pathways to find the best hit for each gene. We also searched the Swiss-Prot and Trembl databases.
Gene family evolution. The coding and protein sequences from eight teleost species with sequenced genomes (O. niloticus, C. semilaevis, D. rerio, G. aculeatus, O. latipes, T. rubripes, T. nigroviridis, and P. olivaceus) and from human (as an outgroup species) were downloaded from the Ensemble genome site (release 60). We constructed a phylogenetic tree for the two flatfish and six other sequenced fish genomes (G. aculeatus, O. latipes, T. rubripes, T. nigroviridis, O. niloticus, and D. rerio) using single-copy gene families. The coding sequences were extracted from each single-copy family and concatenated to one super-gene in the PHYLIP format for each species. Then, PhyML (v3.0) was used to construct the phylogenetic tree 45 . The Bayesian relaxed molecular clock (BRMC) approach 46 was used to estimate species divergence time using the MCMCTREE program in the PAML package. CAFE-2.1 (ref. 47) was used to analyze changes in gene family size. A GO enrichment (MetaGO_ 20120303.RData object from GO.db package) analysis with a P value < 0.05 was performed on genes from the Japanese flounder expanding and contracting families.

Positive selection.
All related genes of other species were downloaded from Ensembl release 64 for the positive selection analysis. MUSCLE was used for multiprotein sequence alignment among the flatfish genes and their orthologs using human proteins as reference. G blocks 48 was used to eliminate poorly aligned positions and divergent regions in the DNA alignments, so they were suitable for the phylogenetic analysis. The filtered sequences were transformed to phylip format, and PhyML 45 was used to reconstruct the phylogenetic tree with fourfold degenerate sites of singlecopy orthologous genes using the HKY85 model. A branch model was used to detect the average ω across the tree (ω0), ω of the appointed branch to test (ω2), and ω of all other branches (ω1) using the following parameter settings: Codonfreq = 2, κ = 2.5, and initial ω = 0.2. Then, the chi-square test was used to check whether ω2 was significantly higher than ω1 and ω0, which would imply that the genes in the appointed branch had changed or evolved quickly. After obtaining positively selected genes (ω2 > ω1 and P < 0.05), we performed a final manual check for genes of the phototransduction and the TH and RA signaling pathways.
Whole-mount immunohistochemistry and immunofluorescence. The negative controls of larvae used for immunohistochemistry were treated with a secondary antibody but not with the primary antibody. Whole-mount preparations were analyzed for triiodothyronine (T3) or thyroxine (T4) immunohistochemistry using horseradish peroxidase (HRP)-conjugated goat antimouse secondary antibody (Immunology Consultants Laboratory, Newberg, OR, USA) or fluorescein isothiocyanate (FITC)-conjugated goat anti-mouse secondary antibody (Fitzgerald) 49 (Supplementary Note). Whole-mount specimens were observed and photographed under a dissecting microscope (Olympus, Tokyo, Japan) and documented using Image-Pro Plus6.0 software (Media Cybernetics, Bethesda, MD, USA). After the specimens had been rinsed in blocking solution (Roche), PBSTX (phosphate-buffered saline with Tween 20 and 5% Triton X-100), and PBS, larval immunofluorescence was observed and photographed under an LSM710 NL0 laser scanning microscope (Carl Zeiss, Oberkochen, Germany).
To study the effect of RA microinjection on asymmetrical body coloration in Japanese flounder, we divided the larvae into control (DMSO, RA vehicle), retinol, ATRA, and 9-cis-RA groups. The larvae were maintained at 16-18 °C in 10-L seawater containers kept in a controlled indoor rearing environment. The larvae (n = 100 per group) were microinjected at 28,31,34,37,40,43,47 and 50 DAH with ATRA, 9-cis-RA, and retinol (2 mg/ml in DMSO) under the skin on the blind side (Supplementary Fig. 38).
Whole-mount in situ RNA hybridization. Total RNA was extracted from Japanese flounder larvae and reverse-transcribed using M-MLV reverse transcriptase (Promega, Madison, WI, USA) according to the manufacturer's instructions. The DNA was subjected to PCR using the primer combinations in Supplementary Table 33 and the following PCR reaction conditions: 95 °C for 5 min followed by 40 cycles of 95 °C for 1 min, annealing temperature for 1 min, 72 °C for 1 min, and a final extension at 72 °C for 5 min. The PCR products were purified using a DNA gel extraction kit (Axygen, Tewksbury, MA, USA) and cloned into the pGEM-T Easy Vector system (Promega) according to the manufacturer's instructions. The validity and orientation of the cloned products within the plasmids were confirmed by sequencing. Both plasmids contained inserts in the same orientation in-frame with the T7 promoter and complementary to the SP6 promoter. Both sense and antisense RNA probes were prepared by in vitro transcription from the T7 and SP6 promoters, incorporating digoxigenin UTP (Roche). The plasmids were linearized with NotI or NcoI restriction enzymes to produce sense or antisense transcripts, respectively. The concentration of each labeled probe was determined by agarose gel electrophoresis, as described in the DIG RNA labeling kit (SP6/T7; Roche). The whole-mount in situ RNA hybridization was performed as previously described 49 and is detailed in the Supplementary Note. Two-and three-color fluorescence in situ RNA hybridization. The retinoid x receptor (rxrs) and retinoic acid receptor (rars) RNA probes were prepared for two-color FISH by in vitro transcription from the T7 and SP6 promoters, incorporating the Digoxigenin RNA labeling mix (Roche). We labeled related mRNAs with digoxigenin, biotin, and fluorescence for three-color FISH. rxra, rarb, rh1, and opn4 were labeled with DIG RNA labeling mix (Roche), and sws2 and kit were labeled with Biotin RNA labeling mix (Roche). sws2, kit, and dct were labeled with Fluorescence RNA labeling mix (Roche). Whole-mount in situ hybridization was carried out and processed using FISH methodology 50 (Supplementary Note). The signal from the antisense probe was photographed using an LSM 710 confocal microscope (Carl Zeiss, Jena, Germany), whereas no signal was detected with the sense probe.
Whole-mount in situ cell proliferation assay. Larvae were treated with BrdU (Sigma-Aldrich) for 8 h, anesthetized in 0.1 µg/µl MS-222, fixed in 4% PFA overnight, and stored in methanol at 4 °C. Larvae that were not labeled with BrdU were treated as negative controls. Whole-mount preparations were analyzed for BrdU immunostaining as previously described 30 (Supplementary Note). Whole-mount specimens were observed and photographed using an SZX7 dissecting microscope (Olympus) and documented using Image-Pro Plus 6.0 software. Irradiation treatment. Rectangular transparent plastic tanks (35 × 25 × 15 cm) were used to irradiate pre-metamorphic larvae. Light-emitting diodes (white, red, orange, yellow, green, blue, and violet) were set above and below the tanks to irradiate the larvae from two directions. We maintained the light intensity at the bottom of the tanks at 2,000 ± 100 Lux by changing the distance between the lights and the tanks. Irradiation time was 12 h/d, and water temperature was 24 ± 2 °C. Three parallel experiments were set up for each light color, and each experiment included 16 DAH individuals (n = 100). The irradiation treatment was carried out in an 80-m 2 dark room.

Exogenous THs rescue the nearly lethal effect of MMI treatment.
Pre-metamorphic Japanese flounder larvae at 20 DAH, which had been reared in 10 liters of water at a constant temperature (19 °C), were separated into four groups (three treatment groups and one control group) containing 200 larvae per group. Another treatment group (group B), containing 800 larvae, was cultured in 30 liters of water under the same conditions. The larvae of the treatment groups were treated with 0.3 g/l MMI for 4 d, with one-third of the water replaced daily, and keeping the MMI concentration constant. Because of the high mortality of the Japanese flounder in the MMI treatment groups, we combined the three treatment groups containing 200 larvae into a single group (group A). Thereafter, the larvae were rescued by treatment with 10 p.p.b. T3 for 5 d. In all the experiments the larvae were fed brine shrimp nauplii and rotifers.
Yeast two-hybrid assay. The coding regions of thraa and thrb1 of Japanese flounder were amplified using RT-PCR and ligated into pGBKT7 to generate pGBKT7-thraa and pGBKT7-thrb1. Similarly, the coding regions of rxra and rxrb of Japanese flounder were amplified and ligated into pGADT7 to generate pGADT7-rxra and pGADT7-rxrb. In the self-activation assay, Y190 yeast cells