Human endogenous retroviruses (HERVs) and other long terminal repeat (LTR)-type retrotransposons (HERV/LTRs) have regulatory elements that possibly influence the transcription of host genes. We systematically identified and characterized these regulatory elements based on publicly available datasets of ChIP-Seq of 97 transcription factors (TFs) provided by ENCODE and Roadmap Epigenomics projects. We determined transcription factor-binding sites (TFBSs) using the ChIP-Seq datasets and identified TFBSs observed on HERV/LTR sequences (HERV-TFBSs). Overall, 794,972 HERV-TFBSs were identified. Subsequently, we identified “HERV/LTR-shared regulatory element (HSRE),” defined as a TF-binding motif in HERV-TFBSs, shared within a substantial fraction of a HERV/LTR type. HSREs could be an indication that the regulatory elements of HERV/LTRs are present before their insertions. We identified 2,201 HSREs, comprising specific associations of 354 HERV/LTRs and 84 TFs. Clustering analysis showed that HERV/LTRs can be grouped according to the TF binding patterns; HERV/LTR groups bounded to pluripotent TFs (e.g., SOX2, POU5F1, and NANOG), embryonic endoderm/mesendoderm TFs (e.g., GATA4/6, SOX17, and FOXA1/2), hematopoietic TFs (e.g., SPI1 (PU1), GATA1/2, and TAL1), and CTCF were identified. Regulatory elements of HERV/LTRs tended to locate nearby and/or interact three-dimensionally with the genes involved in immune responses, indicating that the regulatory elements play an important role in controlling the immune regulatory network. Further, we demonstrated subgroup-specific TF binding within LTR7, LTR5B, and LTR5_Hs, indicating that gains or losses of the regulatory elements occurred during genomic invasions of the HERV/LTRs. Finally, we constructed dbHERV-REs, an interactive database of HERV/LTR regulatory elements (http://herv-tfbs.com/). This study provides fundamental information in understanding the impact of HERV/LTRs on host transcription, and offers insights into the transcriptional modulation systems of HERV/LTRs and ancestral HERVs.

Animal body color is generated primarily by neural crest-derived pigment cells in the skin. Mammals and birds have only melanocytes on the surface of their bodies; however, fish have a variety of pigment cell types or chromatophores, including melanophores, xanthophores, and iridophores. The medaka has a unique chromatophore type called the leucophore. The genetic basis of chromatophore diversity remains poorly understood. Here, we report that three loci in medaka, namely, leucophore free (lf), lf-2, and white leucophore (wl), which affect leucophore and xanthophore differentiation, encode solute carrier family 2, member 15b (slc2a15b), paired box gene 7a (pax7a), and solute carrier family 2 facilitated glucose transporter, member 11b (slc2a11b), respectively. Because lf-2, a loss-of-function mutant for pax7a, causes defects in the formation of xanthophore and leucophore precursor cells, pax7a is critical for the development of the chromatophores. This genetic evidence implies that leucophores are similar to xanthophores, although it was previously thought that leucophores were related to iridophores, as these chromatophores have purine-dependent light reflection. Our identification of slc2a15b and slc2a11b as genes critical for the differentiation of leucophores and xanthophores in medaka led to a further finding that the existence of these two genes in the genome coincides with the presence of xanthophores in nonmammalian vertebrates: birds have yellow-pigmented irises with xanthophore-like intracellular organelles. Our findings provide clues for revealing diverse evolutionary mechanisms of pigment cell formation in animals.

Recent studies have revealed that a cilium-generated liquid flow in the node has a crucial role in the establishment of the left-right (LR) axis in the mouse. In fish, Kupffer's vesicle (KV), a teleost-specific spherical organ attached to the tail region, is known to have an equivalent role to the mouse node during LR axis formation. However, at present, there has been no report of an asymmetric gene expressed in KV under the control of fluid flow. Here we report the earliest asymmetric gene in teleost KV, medaka charon, and its regulation. Charon is a member of the Cerberus/DAN family of proteins, first identified in zebrafish. Although zebrafish charon was reported to be symmetrically expressed in KV, medaka charon displays asymmetric expression with more intense expression on the right side. This asymmetric expression was found to be regulated by KV flow because symmetric and up-regulated charon expression was observed in flow-defective embryos with immotile cilia or disrupted KV. Taken together, medaka charon is a reliable gene marker for LR asymmetry in KV and thus, will be useful for the analysis of the early steps downstream of the fluid flow.

Mechanisms generating diverse cell types from multipotent progenitors are crucial for normal development. Neural crest cells (NCCs) are multipotent stem cells that give rise to numerous cell-types, including pigment cells. Medaka has four types of NCC-derived pigment cells (xanthophores, leucophores, melanophores and iridophores), making medaka pigment cell development an excellent model for studying the mechanisms controlling specification of distinct cell types from a multipotent progenitor. Medaka many leucophores-3 (ml-3) mutant embryos exhibit a unique phenotype characterized by excessive formation of leucophores and absence of xanthophores. We show that ml-3 encodes sox5, which is expressed in premigratory NCCs and differentiating xanthophores. Cell transplantation studies reveal a cell-autonomous role of sox5 in the xanthophore lineage. pax7a is expressed in NCCs and required for both xanthophore and leucophore lineages; we demonstrate that Sox5 functions downstream of Pax7a. We propose a model in which multipotent NCCs first give rise to pax7a-positive partially fate-restricted intermediate progenitors for xanthophores and leucophores; some of these progenitors then express sox5, and as a result of Sox5 action develop into xanthophores. Our results provide the first demonstration that Sox5 can function as a molecular switch driving specification of a specific cell-fate (xanthophore) from a partially-restricted, but still multipotent, progenitor (the shared xanthophore-leucophore progenitor).