Uracil-DNA glycosylase is involved in DNA demethylation and required for embryonic development in the zebrafish embryo

Characterization of Heat-labile Uracil-DNA Glycosylase from Oncorhynchus mykiss and its Application for Carry-over Contamination Control in RT-qPCR

BACKGROUND

Heat-labile uracil-DNA glycosylase (HL-UDG) is commonly employed to eliminate carry-over contamination in DNA amplifications. However, the prevailing HL-UDG is markedly inactivated at 50°C, rendering it unsuitable for specific one-step RT-qPCR protocols utilizing reverse transcriptase at an optimal temperature of 42°C.

OBJECTIVE

This study aimed to explore novel HL-UDG with lower inactivation temperature and for recombinant expression.

METHODS

The gene encoding an HL-UDG was cloned from the cold-water fish rainbow trout (Oncorhynchus mykiss) and expressed in Escherichia coli with high yield. The thermostability of this enzyme and other enzymatic characteristics were thoroughly examined. The novel HL-UDG was then applied for controlling carry-over contamination in one-step RT-qPCR.

RESULTS

This recombinantly expressed truncated HL-UDG of rainbow trout (OmUDG) exhibited high amino acids similarity (84.1% identity) to recombinant Atlantic cod UDG (rcUDG) and was easily denatured at 40°C. The optimal pH of OmUDG was 8.0, and the optimal concentrations of both Na+ and K+ were 10 mM. Since its inactivation temperature was lower than that of rcUDG, the OmUDG could be used to eliminate carry-over contamination in one-step RT-qPCR with moderate reverse transcription temperature.

CONCLUSION

We successfully identified and recombinantly expressed a novel HL-UDG with an inactivation temperature of 40°C. It is suitable for eliminating carry-over contamination in one-step RT-qPCR.

In vivo DNA methylation editing in zebrafish

The CRISPR/dCas9-based epigenome editing technique has driven much attention. Fused with a catalytic domain from Dnmt or Tet protein, the CRISPR/dCas9-DnmtCD or -TetCD systems possess the targeted DNA methylation editing ability and have established a series of in vitro and in vivo disease models. However, no publication has been reported on zebrafish (Danio rerio), an important animal model in biomedicine. The present study demonstrated that CRISPR/dCas9-Dnmt7 and -Tet2 catalytic domain fusions could site-specifically edit genomic DNA methylation in vivo in zebrafish and may serve as an efficient toolkit for DNA methylation editing in the zebrafish model.

DNA repair genes play a variety of roles in the development of fish embryos

Embryogenesis is one of the most important life stages because it determines an organism's healthy growth. However, embryos of externally fertilizing species, such as most fish, are directly exposed to the environment during development and may be threatened by DNA damaging factors (pollutants, UV, reactive oxygen species). To counteract the negative effects of DNA fragmentation, fish embryos evolved complex damage response pathways. DNA repair pathways have been extensively studied in some fish species, such as zebrafish (Danio rerio). Our literature review, on the other hand, revealed a paucity of knowledge about DNA damage response and repair in non-model aquaculture fish species. Further, several pieces of evidence underlie the additional role of DNA repair genes and proteins in organogenesis, spatiotemporal localization in different tissue, and its indispensability for normal embryo development. In this review, we will summarize features of different DNA repair pathways in course of fish embryo development. We describe how the expression of DNA repair genes and proteins is regulated during development, their organogenetic roles, and how the expression of DNA repair genes changes in response to genotoxic stress. This will aid in addressing the link between genotoxic stress and embryo phenotype. Furthermore, available data indicate that embryos can repair damaged DNA, but the effects of early-life stress may manifest later in life as behavioral changes, neoplasia, or neurodegeneration. Overall, we conclude that more research on DNA repair in fish embryos is needed.

Comprehensive maturity of nuclear pore complexes regulates zygotic genome activation

Nuclear pore complexes (NPCs) are channels for nucleocytoplasmic transport of proteins and RNAs. However, it remains unclear whether composition, structure, and permeability of NPCs dynamically change during the cleavage period of vertebrate embryos and affect embryonic development. Here, we report that the comprehensive NPC maturity (CNM) controls the onset of zygotic genome activation (ZGA) during zebrafish early embryogenesis. We show that more nucleoporin proteins are recruited to and assembled into NPCs with development, resulting in progressive increase of NPCs in size and complexity. Maternal transcription factors (TFs) transport into nuclei more efficiently with increasing CNM. Deficiency or dysfunction of Nup133 or Ahctf1/Elys impairs NPC assembly, maternal TFs nuclear transport, and ZGA onset, while nup133 overexpression promotes these processes. Therefore, CNM may act as a molecular timer for ZGA by controlling nuclear transport of maternal TFs that reach nuclear concentration thresholds at a given time to initiate ZGA.

Eliminating predictable DNA off-target effects of cytosine base editor by using dual guiders including sgRNA and TALE

Predictable DNA off-target effect is one of the major safety concerns for the application of cytosine base editors (CBEs). To eliminate Cas9-dependent DNA off-target effects, we designed a novel effective CBE system with dual guiders by combining CRISPR with transcription activator-like effector (TALE). In this system, Cas9 nickase (nCas9) and cytosine deaminase are guided to the same target site to conduct base editing by single-guide RNA (sgRNA) and TALE, respectively. However, if nCas9 is guided to a wrong site by sgRNA, it will not generate base editing due to the absence of deaminase. Similarly, when deaminase is guided to a wrong site by TALE, base editing will not occur due to the absence of single-stranded DNA. In this way, Cas9- and TALE-dependent DNA off-target effects could be completely eliminated. Furthermore, by fusing TALE with YE1, a cytidine deaminase with minimal Cas9-independent off-target effect, we established a novel CBE that could induce efficient C-to-T conversion without detectable Cas9- or TALE-dependent DNA off-target mutations.

Transgenerational effects of parental bisphenol S exposure on zebrafish (Danio rerio) reproduction

Bisphenol S (BPS) is extensively used for production of polycarbonates and other commodities, and is often detected in environment and biota. Parental BPS exposure has been reported to interfere with reproductive development of offspring, but limited information is available on its multigenerational reproductive toxicity. In our present study, zebrafish (Danio rerio) were exposed to BPS (1 and 100 μg/L) from 3 hpf to 120 dpf, and the effects on reproduction, sex steroid hormones, DNA methylation levels and gene transcription involved in steroidogenesis and DNA methylation were investigated in unexposed F1-2 offspring. The results showed that 100 μg/L BPS exposure increased DNA methylation in F1 testes, and 1 μg/L BPS led to DNA methylation in F2 ovaries. The increased DNA methylation levels led to decreased expression of steroidogenic enzymes, including cyp11a, cyp17 and 3βhsd, which might be a main reason for the elevated plasma 17β-estradiol and decreased testosterone levels. In addition, sex ratio indicated a female dominance trend, and reproductive capacity of male fish was severely impaired. Overall, these findings suggest that parental BPS exposure impairs reproductive development of unexposed offspring via DNA methylation and BPS-induced epigenetic modification inheritance has a long-term effect on the fitness and sustainability of fish populations.

DNA methylation maintenance at the p53 locus initiates biliary-mediated liver regeneration

In cases of extensive liver injury, biliary epithelial cells (BECs) dedifferentiate into bipotential progenitor cells (BPPCs), then redifferentiate into hepatocytes and BECs to accomplish liver regeneration. Whether epigenetic regulations, particularly DNA methylation maintenance enzymes, play a role in this biliary-mediated liver regeneration remains unknown. Here we show that in response to extensive hepatocyte damages, expression of dnmt1 is upregulated in BECs to methylate DNA at the p53 locus, which represses p53 transcription, and in turn, derepresses mTORC1 signaling to activate BEC dedifferentiation. After BEC dedifferentiation and BPPC formation, DNA methylation at the p53 locus maintains in BPPCs to continue blocking p53 transcription, which derepresses Bmp signaling to induce BPPC redifferentiation. Thus, this study reveals promotive roles and mechanisms of DNA methylation at the p53 locus in both dedifferentiation and redifferentiation stages of biliary-mediated liver regeneration, implicating DNA methylation and p53 as potential targets to stimulate regeneration after extensive liver injury.

Large-scale generation and phenotypic characterization of zebrafish CRISPR mutants of DNA repair genes

A systematic knowledge of the roles of DNA repair genes at the level of the organism has been limited due to the lack of appropriate experimental approaches using animal model systems. Zebrafish has become a powerful vertebrate genetic model system with availability due to the ease of genome editing and large-scale phenotype screening. Here, we generated zebrafish mutants for 32 DNA repair and replication genes through multiplexed CRISPR/Cas9-mediated mutagenesis. Large-scale phenotypic characterization of our mutant collection revealed that three genes (atad5a, ddb1, pcna) are essential for proper embryonic development and hematopoiesis; seven genes (apex1, atrip, ino80, mre11a, shfm1, telo2, wrn) are required for growth and development during juvenile stage and six genes (blm, brca2, fanci, rad51, rad54l, rtel1) play critical roles in sex development. Furthermore, mutation in six genes (atad5a, brca2, polk, rad51, shfm1, xrcc1) displayed hypersensitivity to DNA damage agents. Our zebrafish mutant collection provides a unique resource for understanding of the roles of DNA repair genes at the organismal level.

Widespread Enhancer Dememorization and Promoter Priming during Parental-to-Zygotic Transition

The epigenome plays critical roles in controlling gene expression and development. However, how the parental epigenomes transit to the zygotic epigenome in early development remains elusive. Here we show that parental-to-zygotic transition in zebrafish involves extensive erasure of parental epigenetic memory, starting with methylating gametic enhancers. Surprisingly, this occurs even prior to fertilization for sperm. Both parental enhancers lose histone marks by the 4-cell stage, and zygotic enhancers are not activated until around zygotic genome activation (ZGA). By contrast, many promoters remain hypomethylated and, unexpectedly, acquire histone acetylation before ZGA at as early as the 4-cell stage. They then resolve into either activated or repressed promoters upon ZGA. Maternal depletion of histone acetyltransferases results in aberrant ZGA and early embryonic lethality. Finally, such reprogramming is largely driven by maternal factors, with zygotic products mainly contributing to embryonic enhancer activation. These data reveal widespread enhancer dememorization and promoter priming during parental-to-zygotic transition.

AP endonuclease 1 (Apex1) influences brain development linking oxidative stress and DNA repair

Brain and neurons are particularly sensitive to reactive oxygen species (ROS). Oxidative damage from ROS results in increased 8-oxoguanine in DNA followed by repair through the base excision repair (BER) pathway. We reported earlier that AP endonuclease 1 (Apex1) not only participates directly in BER but also regulates transcription factor Creb1. Here, we investigated how Apex1 affects brain to respond effectively to oxidative damage during zebrafish development. Loss of Apex1 resulted in increased ROS, 8-oxoguanine, and abasic sites as well as loss of Ogg1, which recognizes 8-oxoguanine and is required for its repair. Moreover, knock-down of Apex1 not only resulted in reduction of expression of several major proteins in the BER pathway (Polb and Ogg1), and it also resulted in maldistribution and loss of four key brain transcription factors (fezf2, otx2, egr2a, and pax2a), leading to abnormal brain development. These results were independent of p53 protein level. In contrast, exposure to exogenous H₂O₂ resulted in increased transcription and protein of Apex1 along with other BER components, as well as Creb1. Taken together, these results indicate that oxidative stress increased when the level of Apex1 was reduced, revealing a novel pathway of how Apex1 manages oxidative stress in developing brain.

The genetic program of oocytes can be modified in vivo in the zebrafish ovary

Oocytes, the irreplaceable gametes for generating a new organism, are matured in the ovary of living female animals. It is unknown whether any genetic manipulations can be applied to immature oocytes inside the living ovaries. As a proof-of-concept, we here demonstrate genetic amendments of zebrafish immature oocytes within the ovary. Oocyte microinjection in situ (OMIS) stimulates tissue repair responses, but some of the microinjected immature oocytes are matured, ovulated and fertilizable. By OMIS-mediated Cas9 approach, ntla and gata5 loci of oocytes arrested at prophase I of meiosis are successfully edited before fertilization. Through OMIS, high efficiency of biallelic mutations in single or multiple loci using Cas9/gRNAs allows immediate manifestation of mutant phenotypes in F0 embryos and multiple transgenes can co-express the reporters in F0 embryos with patterns similar to germline transgenic embryos. Furthermore, maternal knockdown of dnmt1 by antisense morpholino via OMIS results in a dramatic decrease of global DNA methylation level at the dome stage and causes embryonic lethality prior to segmentation period. Therefore, OMIS opens a door to efficiently modify the genome and provides a possibility to repair genetically abnormal oocytes in situ.

The guanine nucleotide exchange factor Net1 facilitates the specification of dorsal cell fates in zebrafish embryos by promoting maternal β-catenin activation

Wnt/β-catenin signaling is essential for the initiation of dorsal-ventral patterning during vertebrate embryogenesis. Maternal β-catenin accumulates in dorsal marginal nuclei during cleavage stages, but its critical target genes essential for dorsalization are silent until mid-blastula transition (MBT). Here, we find that zebrafish net1, a guanine nucleotide exchange factor, is specifically expressed in dorsal marginal blastomeres after MBT, and acts as a zygotic factor to promote the specification of dorsal cell fates. Loss- and gain-of-function experiments show that the GEF activity of Net1 is required for the activation of Wnt/β-catenin signaling in zebrafish embryos and mammalian cells. Net1 dissociates and activates PAK1 dimers, and PAK1 kinase activation causes phosphorylation of S675 of β-catenin after MBT, which ultimately leads to the transcription of downstream target genes. In summary, our results reveal that Net1-regulated β-catenin activation plays a crucial role in the dorsal axis formation during zebrafish development.

Abundant DNA 6mA methylation during early embryogenesis of zebrafish and pig

DNA N⁶-methyldeoxyadenosine (6mA) is a well-known prokaryotic DNA modification that has been shown to exist and play epigenetic roles in eukaryotic DNA. Here we report that 6mA accumulates up to ∼0.1–0.2% of total deoxyadenosine during early embryogenesis of vertebrates, but diminishes to the background level with the progression of the embryo development. During this process a large fraction of 6mAs locate in repetitive regions of the genome.

DNA 6mA is a poorly understood epigenetic mark present at a low abundance in eukaryotic genomes. Here the authors observe high levels in zebrafish and pig during early embryogenesis enriched to repetitive regions of the genome and followed by attenuation during development.

Impact of sperm DNA damage and oocyte-repairing capacity on trout development

Zygotic repair of paternal DNA is essential during embryo development. In spite of the interest devoted to sperm DNA damage, its combined effect with defect-repairing oocytes has not been analyzed. Modification of the breeding season is a common practice in aquaculture. This practice reduces developmental success and could affect the both factors: sperm DNA integrity and oocyte repair capacity. To evaluate the maternal role, we analyzed the progeny outcome after fertilizing in-season trout oocytes with untreated and with UV-irradiated sperm. We also analyzed the offspring obtained out of season with untreated sperm. The analysis of the number of lesions in 4 sperm nuclear genes revealed an increase of 1.22–11.18 lesions/10 kb in out-of-season sperm, similar to that obtained after sperm UV irradiation (400 µW/cm25 min). Gene expression showed in out-of-season oocytes the overexpression of repair genes (ogg1, ung, lig3, rad1) and downregulation of tp53, indicating an enhanced repairing activity and reduced capacity to arrest development upon damage. The analysis of the progeny in out-of-season embryos revealed a similar profile tolerant to DNA damage, leading to a much lower apoptotic activity at organogenesis, lower hatching rates and increased rate of malformations. The effects were milder in descendants from in-season-irradiated sperm, showing an enhanced repairing activity at epibolia. Results point out the importance of the repairing machinery provided by the oocyte and show how susceptible it is to environmental changes. Transcripts related to DNA damage signalization and repair could be used as markers of oocyte quality.

Base excision DNA repair in the embryonic development of the sea urchin, Strongylocentrotus intermedius

In actively proliferating cells, such as the cells of the developing embryo, DNA repair is crucial for preventing the accumulation of mutations and synchronizing cell division. Sea urchin embryo growth was analyzed and extracts were prepared. The relative activity of DNA polymerase, apurinic/apyrimidinic (AP) endonuclease, uracil-DNA glycosylase, 8-oxoguanine-DNA glycosylase, and other glycosylases was analyzed using specific oligonucleotide substrates of these enzymes; the reaction products were resolved by denaturing 20% polyacrylamide gel electrophoresis. We have characterized the profile of several key base excision repair activities in the developing embryos (2 blastomers to mid-pluteus) of the grey sea urchin, Strongylocentrotus intermedius. The uracil-DNA glycosylase specific activity sharply increased after blastula hatching, whereas the specific activity of 8-oxoguanine-DNA glycosylase steadily decreased over the course of the development. The AP-endonuclease activity gradually increased but dropped at the last sampled stage (mid-pluteus 2). The DNA polymerase activity was high at the first cleavage division and then quickly decreased, showing a transient peak at blastula hatching. It seems that the developing sea urchin embryo encounters different DNA-damaging factors early in development within the protective envelope and later as a free-floating larva, with hatching necessitating adaptation to the shift in genotoxic stress conditions. No correlation was observed between the dynamics of the enzyme activities and published gene expression data from developing congeneric species, S. purpuratus. The results suggest that base excision repair enzymes may be regulated in the sea urchin embryos at the level of covalent modification or protein stability.

Active DNA demethylation by DNA repair: Facts and uncertainties

Pathways that control and modulate DNA methylation patterning in mammalian cells were poorly understood for a long time, although their importance in establishing and maintaining cell type-specific gene expression was well recognized. The discovery of proteins capable of converting 5-methylcytosine (5mC) to putative substrates for DNA repair introduced a novel and exciting conceptual framework for the investigation and ultimate discovery of molecular mechanisms of DNA demethylation. Against the prevailing notion that DNA methylation is a static epigenetic mark, it turned out to be dynamic and distinct mechanisms appear to have evolved to effect global and locus-specific DNA demethylation. There is compelling evidence that DNA repair, in particular base excision repair, contributes significantly to the turnover of 5mC in cells. By actively demethylating DNA, DNA repair supports the developmental establishment as well as the maintenance of DNA methylation landscapes and gene expression patterns. Yet, while the biochemical pathways are relatively well-established and reviewed, the biological context, function and regulation of DNA repair-mediated active DNA demethylation remains uncertain. In this review, we will thus summarize and critically discuss the evidence that associates active DNA demethylation by DNA repair with specific functional contexts including the DNA methylation erasure in the early embryo, the control of pluripotency and cellular differentiation, the maintenance of cell identity, and the nuclear reprogramming.

Detection of uracil within DNA using a sensitive labeling method for in vitro and cellular applications

The role of uracil in genomic DNA has been recently re-evaluated. It is now widely accepted to be a physiologically important DNA element in diverse systems from specific phages to antibody maturation and Drosophila development. Further relevant investigations would largely benefit from a novel reliable and fast method to gain quantitative and qualitative information on uracil levels in DNA both in vitro and in situ, especially since current techniques does not allow in situ cellular detection. Here, starting from a catalytically inactive uracil-DNA glycosylase protein, we have designed several uracil sensor fusion proteins. The designed constructs can be applied as molecular recognition tools that can be detected with conventional antibodies in dot-blot applications and may also serve as in situ uracil-DNA sensors in cellular techniques. Our method is verified on numerous prokaryotic and eukaryotic cellular systems. The method is easy to use and can be applied in a high-throughput manner. It does not require expensive equipment or complex know-how, facilitating its easy implementation in any basic molecular biology laboratory. Elevated genomic uracil levels from cells of diverse genetic backgrounds and/or treated with different drugs can be demonstrated also in situ, within the cell.

Transcription Factors and DNA Repair Enzymes Compete for Damaged Promoter Sites

Transcriptional regulation is a tightly regulated, vital process. The transcription factor cyclic AMP-response element-binding protein 1 (CREB1) controls ∼25% of the mammalian transcriptome by binding the CREB1 binding site consensus sequence (CRE) sequence (TGACGTCA). DNA lesions within CRE modulate CREB1 binding negatively and positively. Because appropriate DNA lesions also interact with base excision repair proteins, we investigated whether CREB1 and repair glycosylases compete with each other. We incubated 39-mer CRE-containing double-stranded oligonucleotides with recombinant CREB1 alone or with UNG2 or OGG1, followed by EMSA. The CpG islet within CRE was modified to contain a G/U or 8-oxoG (°G)/C mispair. OGG1 and CREB1 reversibly competed for CRE containing an °G/C pair. Also, OGG1 blocked CREB1 from dimerizing by 69%, even when total CREB1 binding was reduced only by 20-30%. In contrast, bound CREB1 completely prevented access to G/U-containing CRE by UNG2 and, therefore, to base excision repair, whereas UNG2 exposure prevented CREB1 binding. CREB1 dimerization was unaffected by UNG2 when CREB1 bound to CRE, but was greatly reduced by prior UNG2 exposure. To explore physiological relevance, we microinjected zebrafish embryos with the same oligonucleotides, as a sink for endogenous CREB1. As predicted, microinjection with unmodified or lesion-containing CRE, but not scrambled CRE or scrambled CRE with a G/U mispair, resulted in increased embryo death. However, only the G/U mispair in native CRE resulted in substantial developmental abnormalities, thus confirming the danger of unrepaired G/U mispairs in promoters. In summary, CREB1 and DNA glycosylases compete for damaged CRE in vitro and in vivo, thus blocking DNA repair and resulting in transcriptional misregulation leading to abnormal development.

Uracil-DNA Glycosylase UNG Promotes Tet-mediated DNA Demethylation

In mammals, active DNA demethylation involves oxidation of 5-methylcytosine (5mC) into 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) by Tet dioxygenases and excision of these two oxidized bases by thymine DNA glycosylase (TDG). Although TDG is essential for active demethylation in embryonic stem cells and induced pluripotent stem cells, it is hardly expressed in mouse zygotes and dispensable in pronuclear DNA demethylation. To search for other factors that might contribute to demethylation in mammalian cells, we performed a functional genomics screen based on a methylated luciferase reporter assay. UNG2, one of the glycosylases known to excise uracil residues from DNA, was found to reduce DNA methylation, thus activating transcription of a methylation-silenced reporter gene when co-transfected with Tet2 into HEK293T cells. Interestingly, UNG2 could decrease 5caC from the genomic DNA and a reporter plasmid in transfected cells, like TDG. Furthermore, deficiency in Ung partially impaired DNA demethylation in mouse zygotes. Our results suggest that UNG might be involved in Tet-mediated DNA demethylation.

Si-CSP9 regulates the integument and moulting process of larvae in the red imported fire ant, Solenopsis invicta

Chemosensory proteins (CSPs) have been predicted to be involved in development; however, direct evidence for their involvement is lacking, and genetic basis is largely unknown. To determine the function of the chemosensory protein 9 (Si-CSP9) gene in Solenopsis invicta, we used RNA interference to silence Si-CSP9 in 3rd-instar larvae. The 3rd-instar larvae failed to shed their cuticle after being fed Si-CSP9-directed siRNA, and expression profiling of RNAi-treated and untreated control larvae showed that 375 genes were differentially expressed. Pathway enrichment analysis revealed that 4 pathways associated with larval development were significantly enriched. Blast analysis revealed that one fatty acid amide hydrolase (FAAH) gene was up-regulated and 4 fatty acid synthase (FAT) genes and one protein kinase DC2 gene (PKA) were down-regulated in the enriched pathways. Significantly higher expression of these genes was found in 4th-instar larvae, and Pearson correlation analysis of the expression patterns revealed significant relationships among Si-CSP9, PKA, FAAH, and FAT1-4. Moreover, we confirmed that expression levels of Si-CSP9, FAAH, and FAT1-4 were significantly reduced and that the development of 3rd-instar larvae was halted with PKA silencing. These results suggest that Si-CSP9 and PKA may be involved in the network that contributes to development of 3rd-instar larvae.

Dynamics of DNA hydroxymethylation in zebrafish

During embryonic development in mammals, most of the methylated cytosines in the paternal genome are converted to 5-hydroxymethyldeoxycytidine (hmC), as part of DNA methylation reprogramming. Recent data also suggest tissue-specific functional roles of hmC, perhaps as an epigenetic mark. However, limited data are available on the levels and tissue distribution in zebrafish. In this study, we used high-performance liquid chromatography mass spectrometry to quantify hmC and 5-methyldeoxycytidine (mC) in zebrafish during development and in different tissues of the adult fish. Low levels of mC were found at 0.5 hours postfertilization (hpf) (1-2 cell stage) (1.9%), and increased to 8.4% by 96 hpf, with similar levels observed in different adult tissues. No hmC was detected up to 12 hpf, but levels increased during development from 24 up to 96 hpf (0.23%). In tissues, the highest levels of hmC were found in the brain (0.49%), intermediate levels in muscle (0.13%), liver (0.08%), and intestine (0.06%) and low levels in testis (0.01%), with an inverse correlation between hmC and mC. Our results indicate similar tissue distribution and levels of hmC between zebrafish and mammals, but distinct differences during embryonic development. Although more research is needed, these results support the use of zebrafish as an alternative model in the elucidation of tissue-specific functions of hmC.