The Role of DNA Methylation Reprogramming During Sex Determination and Transition in Zebrafish

DNA methylation is a prevalent epigenetic modification in vertebrates, and it has been shown to be involved the regulation of gene expression and embryo development. However, it remains unclear how DNA methylation regulates sexual development, especially in species without sex chromosomes. To determine this, we utilized zebrafish to investigate DNA methylation reprogramming during juvenile germ cell development and adult female-to-male sex transition.We reveal that primordial germ cells(PGCs) undergo significant DNA methylation reprogramming during germ cell development, and the methylome of PGCs is reset to an oocyte/ovary-like pattern at 9 days post fertilization(9 dpf). When DNA methyltransferase(DNMT) activity in juveniles was blocked after 9 dpf, the zebrafish developed into females. We also show that Tet3 is involved in PGC development. Notably, we find that DNA methylome reprogramming during adult zebrafish sex transition is similar to the reprogramming during the sex differentiation from 9 dpf PGCs to sperm. Furthermore, inhibiting DNMT activity can prevent the female-to-male sex transition, suggesting that methylation reprogramming is required for zebrafish sex transition. In summary, DNA methylation plays important roles in zebrafish germ cell development and sexual plasticity.     

Mapping the zebrafish brain methylome using reduced representation bisulfite sequencing

Reduced representation bisulfite sequencing (RRBS) has been used to profile DNA methylation patterns in mammalian genomes such as human, mouse and rat. The methylome of the zebrafish, an important animal model, has not yet been characterized at base-pair resolution using RRBS. Therefore, we evaluated the technique of RRBS in this model organism by generating four single-nucleotide resolution DNA methylomes of adult zebrafish brain. We performed several simulations to show the distribution of fragments and enrichment of CpGs in different in silico reduced representation genomes of zebrafish. Four RRBS brain libraries generated 98 million sequenced reads and had higher frequencies of multiple mapping than equivalent human RRBS libraries. The zebrafish methylome indicates there is higher global DNA methylation in the zebrafish genome compared with its equivalent human methylome. This observation was confirmed by RRBS of zebrafish liver. High coverage CpG dinucleotides are enriched in CpG island shores more than in the CpG island core. We found that 45% of the mapped CpGs reside in gene bodies, and 7% in gene promoters. This analysis provides a roadmap for generating reproducible base-pair level methylomes for zebrafish using RRBS and our results provide the first evidence that RRBS is a suitable technique for global methylation analysis in zebrafish.     

Sequence analysis and functional study of thymidylate synthase from zebrafish, Danio rerio

The thymidylate synthase (TS), an important target for many anticancer drugs, has been cloned from different species. But the cDNA property and function of TS in zebrafish are not well documented. In order to use zebrafish as an animal model for screening novel anticancer agents, we isolated TS cDNA from zebrafish and compared its sequence with those from other species. The open reading frame (ORF) of zebrafish TS cDNA sequence was 954 nucleotides, encoding a 318-amino acid protein with a calculated molecular mass of 36.15 kDa. The deduced amino acid sequence of zebrafish TS was similar to those from other organisms, including rat, mouse and humans. The zebrafish TS protein was expressed in Escherichia coli and purified to homogeneity. The purified zebrafish TS showed maximal activity at 28 degrees C with similar K(m) value to human TS. Western immunoblot assay confirmed that TS was expressed in all the developmental stages of zebrafish with a high level of expression at the 1-4 cell stages. To study the function of TS in zebrafish embryo development, a short hairpin RNA (shRNA) expression vector, pSilencer 4.1-CMV/TS, was constructed which targeted the protein-coding region of zebrafish TS mRNA. Significant change in the development of tail and epiboly was found in zebrafish embryos microinjected pSilencer4.1-CMV/TS siRNA expression vector.     

DNA甲基化与脊椎动物胚胎发育

DNA甲基化是指DNA甲基转移酶(DNMT)将DNA序列中的5′胞嘧啶转变为5′甲基胞嘧啶的化学修饰, 可以调控基因的时空特异性表达, 从而影响细胞命运决定和分化等生物学过程。近年来研究发现, DNA甲基化在脊椎动物胚胎早期发育中有重要作用, Dnmt基因的缺失会影响胚胎早期发育和多个器官的形成及分化, 如胚胎早期致死、内脏器官和神经系统终末分化缺陷以及血液发生紊乱等。文章总结了DNA甲基化转移酶在小鼠和斑马鱼发育过程中的动态变化, 并系统阐述了DNA甲基化在胚胎早期发育和器官发生中的作用, 重点揭示DNA 甲基化转移酶与组蛋白甲基化转移酶如何协同调控DNA甲基化从而影响基因转录的分子机理。DNA甲基化作为一种关键的表观遗传学因素, 全面系统地理解其在胚胎发育过程中的作用机制对靶向治疗人类相关疾病有一定的理论指导意义。     

DNA methylation regulates gabrb2 mRNA expression: developmental variations and disruptions in l‐methionine‐induced zebrafish with schizophrenia‐like symptoms

Single nucleotide polymorphisms (SNPs) in the human type A gamma‐aminobutyric acid (GABA) receptor β₂ subunit gene (GABRB2) have been associated with schizophrenia and quantitatively correlated with mRNA expression in the postmortem brain tissue of patients with schizophrenia. l ‐Methionine (MET) administration has been reported to cause a recrudescence of psychotic symptoms in patients with schizophrenia, and similar symptoms have been generated in MET‐induced mice. In this study, a zebrafish animal model was used to evaluate the relationship between the gabrb2 mRNA expression and its promoter DNA methylation in developmental and MET‐induced schizophrenia‐like zebrafish. The results indicated developmental increases in global DNA methylation and decreases in gabrb2 promoter methylation in zebrafish. A significant increase in gabrb2 mRNA levels was observed after GABA was synthesized. Additionally, the MET‐triggered schizophrenia‐like symptoms in adult zebrafish, involving social withdrawal and cognitive dysfunction analyzed with social interaction and T‐maze behavioral tests, were accompanied by significantly increased DNA methylation levels in the global genome and the gabrb2 promoter. Furthermore, the significant correlation between gabrb2 mRNA expression and gabrb2 promoter methylation observed in the developmental stages became non‐significant in MET‐triggered adult zebrafish. These findings demonstrate that gabrb2 mRNA expression is associated with DNA methylation varies by developmental stage and show that these epigenetic association mechanisms are disrupted in MET‐triggered adult zebrafish with schizophrenia‐like symptoms. In conclusion, these results provide plausible epigenetic evidence of the GABA_A receptor β₂ subunit involvement in the schizophrenia‐like behaviors and demonstrate the potential use of zebrafish models in neuropsychiatric research.     

Pseudoloma neurophilia n. g., n. sp., a new microsporidium from the central nervous system of the zebrafish (Danio rerio)

An unusual xenoma-forming microsporidium was discovered in the central nervous system of moribund zebrafish from a laboratory colony in Eugene, Oregon. Infected fish were often emaciated and lethargic, and histological examination commonly revealed severe myelitis and myositis associated with the infection. Based on its structure, development, and small subunit ribosomal DNA sequence it is unique among fish microsporidia. Spores are uninucleate, ovoid to pyriform, with a prominent posterior vacuole. Spores average 5.4 x 2.7 microm with 13-16 coils of the polar filament. The microsporidium produces xenomas within the spinal cord and hindbrain of fish, and xenomas contained sporophorous vesicles with up to 16 spores. Sporoblasts and presporoblast stages (probably sporonts) are found occasionally in small aggregates dispersed randomly throughout xenomas. It clustered in the "Ichthyosporidium group" along with other fish microsporidian genera based on rDNA sequence analysis. The rDNA sequence of the zebrafish microsporidium was most similar to that of Ichthyosporidium, but showed only 12.1% similarity and therefore this microsporidium can be considered a distinct genus and species, which we have named Pseudoloma neurophilia n. g., n. sp.     

Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.     

Nucleotide sequence of cDNA clones coding for a brain-type fatty acid binding protein and its tissue-specific expression in adult zebrafish (Danio rerio)

We have determined the nucleotide sequence for two cDNA clones coding for a fatty acid binding protein (FABP) from zebrafish (Danio rerio). Comparison of the sequence with GenBank entries revealed extensive amino acid identity between this zebrafish FABP and brain FABPs (B-FABP) from other species. The zebrafish B-FABP cDNA hybridized to single restriction fragments of total zebrafish genomic DNA digested with the restriction endonucleases BglII or EcoRI suggesting that a single copy of the B-FABP gene is present in the zebrafish genome. Northern blot analysis demonstrated that the zebrafish B-FABP mRNA is approximately 850 nucleotides in length. In situ hybridization revealed that the B-FABP mRNA was expressed in the periventricular gray zone of the optic tectum of the adult zebrafish brain.     

Impact of DNA methyltransferase inhibitor 5-azacytidine on cardiac development of zebrafish in vivo and cardiomyocyte proliferation,apoptosis, and the homeostasis of gene expression in vitro

Cardiac development is a peculiar process involving coordinated cellular differentiation, migration, proliferation, and apoptosis. DNA methylation plays a key role in genomic stability, tissue-specific gene expression, cell proliferation, and apoptosis. Hypomethylation in the global genome has been reported in cardiovascular diseases. However, little is known about the impact and specific mechanism of global hypomethylation on cardiomyocytes. In the present study, we explored the impact of DNA methyltransferase inhibitors 5-azacytidine on cardiac development. In vivo experiment showed that hypomethylation of zebrafish embryos with 5-azacytidine exposure significantly reduced survival, induced malformations, and delayed general development process. Furthermore, zebrafish embryos injected with 5-azacytidine developed pericardial edema, ventricular volume reduction, looping deformity, and reduction in heart rate and ventricular shortening fraction. Cardiomyocytes treated with 5-azacytidine in vitro decreased proliferation and induced apoptosis in a concentration-dependent manner. Furthermore, 5-azacytidine treatment in cardiomyocytes resulted in 20 downregulated genes expression and two upregulated genes expression in 45 candidate genes, which indicated that DNA methylation functions as a bidirectional modulator in regulating gene expression. In conclusion, these results show the regulative effects of the epigenetic modifier 5-azacytidine in cardiac development of zebrafish embryos in vivo and cardiomyocyte proliferation and apoptosis and the homeostasis of gene expression in vitro, which offer a novel understanding of aberrant DNA methylation in the etiology of cardiovascular disease.     

DNA barcodes from century‐old type specimens using next‐generation sequencing

Type specimens have high scientific importance because they provide the only certain connection between the application of a Linnean name and a physical specimen. Many other individuals may have been identified as a particular species, but their linkage to the taxon concept is inferential. Because type specimens are often more than a century old and have experienced conditions unfavourable for DNA preservation, success in sequence recovery has been uncertain. This study addresses this challenge by employing next‐generation sequencing (NGS) to recover sequences for the barcode region of the cytochrome c oxidase 1 gene from small amounts of template DNA. DNA quality was first screened in more than 1800 century‐old type specimens of Lepidoptera by attempting to recover 164‐bp and 94‐bp reads via Sanger sequencing. This analysis permitted the assignment of each specimen to one of three DNA quality categories – high (164‐bp sequence), medium (94‐bp sequence) or low (no sequence). Ten specimens from each category were subsequently analysed via a PCR‐based NGS protocol requiring very little template DNA. It recovered sequence information from all specimens with average read lengths ranging from 458 bp to 610 bp for the three DNA categories. By sequencing ten specimens in each NGS run, costs were similar to Sanger analysis. Future increases in the number of specimens processed in each run promise substantial reductions in cost, making it possible to anticipate a future where barcode sequences are available from most type specimens.     

The spectra of point mutations in vertebrate genomes

In spite of the importance of point mutations for evolution and human diseases, their natural spectrum of incidence in different species is not known. Here I propose to determine these spectra by comparing consecutive sequence periods in stretches of repetitive DNA. The article presents the analysis of more than 51,000 such point mutations identified by this approach in the genomes of human, chimpanzee, rat, mouse, pufferfish, zebrafish, and sea squirt. I propose to explain the observed spectra by auto‐mutagenic mechanisms of genome variation involving the inter‐conversions of nucleotides, single base‐pair inversions and their combinations.     

Short Inverted-Repeat Transposable Elements in Teleost Fish and Implications for a Mechanism of Their Amplification

Angel is the first miniature inverted-repeat transposable element (MITE) isolated from fish. Angel elements are imperfect palindromes with the potential to form stem-loop structures in vitro. Despite sequence divergence of elements of up to 55% within and between species, their inverted repeat structures have been maintained, implying functional importance. We estimate that there are about 10³–10⁴ Angels scattered throughout the zebrafish genome, evidence that this family of transposable elements has been significantly amplified over the course of evolution. Angel elements and Xenopus MITEs carry common sequence motifs at their termini, indicating common origin and/or related mechanisms of transposition. We present a model in which MITEs take advantage of the basic cellular mechanism of DNA replication for their amplification, which is dependent on the characteristic inverted repeat structures of these elements. We propose that MITEs are genomic parasites that transpose via a DNA intermediate, which forms by a folding-back of a single strand of DNA, that borrow all of the necessary factors for their amplification from products encoded in the genomes in which they reside. DNA polymorphisms in different lines of zebrafish were detected by PCR using Angel-specific primers, indicating that such elements, combined with other transposons in vertebrate genomes, will be useful molecular tools for genome mapping and genetic analyses of mutations. Received: 7 April 1998 / Accepted: 7 April 1998     

Zebrafish HSF4: a novel protein that shares features of both HSF1 and HSF4 of mammals

Heat-shock proteins (hsps) have important roles in the development of the eye lens. We previously demonstrated that knockdown of hsp70 gene expression using morpholino antisense technology resulted in an altered lens phenotype in zebrafish embryos. A less severe phenotype was seen with knockdown of heat-shock factor 1 (HSF1), suggesting that, while it likely plays a role in hsp70 regulation during lens formation, other regulatory factors are also involved. Heat-shock factor 4 plays an important role in mammalian lens development, and an expressed sequence tag encoding zebrafish HSF4 has been identified. The deduced amino acid sequence shares structural similarities with mammalian HSF4 including the lack of an HR-C domain. However, the HR-C domain is absent due to a severe C-terminal truncation within zebrafish HSF4 (zHSF4) relative to the mammalian protein. Surprisingly, the amino acid composition of the zHSF4 DNA binding domain shares a greater degree of identity with HSF1 proteins than it does with mammalian HSF4 proteins. Consistent with this, the binding affinity of in vitro synthesized zHSF4 for discontinuous heat-shock response element sequences is more limited, similar to what has been previously observed for HSF1 proteins. Hsf4 mRNA is expressed in zebrafish adult eye tissue but is only observed in developing embryonic tissue at 60 h post-fertilization or later. This, together with the lack of an observable phenotype following morpholino-based antisense knockdown of hsf4, suggests that zHSF4 is unlikely to play a role in regulating early embryonic lens development.     

Probing the role of sequence in the assembly of three‐dimensional DNA crystals

DNA is a widely used biopolymer for the construction of nanometer‐scale objects due to its programmability and structural predictability. One long‐standing goal of the DNA nanotechnology field has been the construction of three‐dimensional DNA crystals. We previously determined the X‐ray crystal structure of a DNA 13‐mer that forms a continuously hydrogen bonded three‐dimensional lattice through Watson‐Crick and non‐canonical base pairs. Our current study sets out to understand how the sequence of the Watson‐Crick duplex region influences crystallization of this 13‐mer. We screened all possible self‐complementary sequences in the hexameric duplex region and found 21 oligonucleotides that crystallized. Sequence analysis showed that one specific Watson‐Crick pair influenced the crystallization propensity and the speed of crystal self‐assembly. We determined X‐ray crystal structures for 13 of these oligonucleotides and found sequence‐specific structural changes that suggests that this base pair may serve as a structural anchor during crystal assembly. Finally, we explored the crystal self‐assembly and nucleation process. Solution studies indicated that these oligonucleotides do not form base pairs in the absence of cations, but that the addition of divalent cations leads to rapid self‐assembly to higher molecular weight complexes. We further demonstrate that crystals grown from mixtures of two different oligonucleotide sequences contain both oligonucleotides. These results suggest that crystal self‐assembly is nucleated by the formation of the Watson‐Crick duplexes initiated by a simple chemical trigger. This study provides new insight into the role of sequence for the assembly of periodic DNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 618–626, 2015.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).