Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

Genome-wide identification, evolution of DNA methyltransferases and their expression under salinity stress in Larimichthys crocea

DNA methyltransferases (Dnmts) are responsible for DNA methylation which influences patterns of gene expression and plays a crucial role in response to environmental changes. In this study, 7 LcDnmt genes were identified in the genome of large yellow croaker (Larimichthys crocea). The comprehensive analysis was conducted on gene structure, protein and location site of LcDnmts. LcDnmt proteins belonged to three groups (Dnmt1, Dnmt2, and Dnmt3) according to their conserved domains and phylogenetic analysis. Although Dnmt3 can be further divided into three sub groups (Dnmt3a, Dnmt3b, and Dnmt3l), there is no Dnmnt3l member in the large yellow croaker. Phylogenetic analysis revealed that the Dnmt family was highly conserved in teleosts. Expression patterns derived from the RNA-seq, qRT-PCR and Western blot analysis revealed that 2 LcDnmt genes (LcDnmt1 and LcDnmt3a2) significantly regulated under salinity stress in the liver, which was found to be dominantly expressed in the intestine and brain, respectively. These two genes may play an important role in the salinity stress of large yellow croaker and represent candidates for future functional analysis. Our results revealed the conservation of Dnmts during evolution and indicated a potential role of Dnmts in epigenetic regulation of response to salinity stress.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Gene expression profiling of DNA methyltransferase genes in Siniperca chuatsi based on transcriptome sequencing

The mandarin fish (Siniperca chuatsi) DNA methyltransferase gene 1 (dnmt1) was highly expressed in the mesonephros, head kidney and gonad, whereas dnmt2 was expressed in most tissues. dnmt3a was highly expressed in the brain and spleen, but dnmt3b was mainly expressed in the brain and head kidney. The genes dnmt1 and dnmt2 were highly expressed in the early stages of embryonic development, and dnmt3a and dnmt3b were expressed later. These genes also showed certain changes after artificial diet acclimation, salinity adaptation and immune stress.

Sexual dimorphism of DNA and histone methylation profiles in the gonads of the olive flounder Paralichthys olivaceus

DNA methylation and histone methylation are two types of the most important epigenetic modifications. However, research on their differential expression in gonads of male and female fish is limited. In this study, we examined the characteristics of DNA methylation and tri-methylation of lysine 4 of histone H3 (H3K4me3) modification profiles in the gonads of the wild-type and meio-gynogenetic olive flounders Paralichthys olivaceus. Enzyme-linked immunosorbent assay (ELISA) analysis revealed that the global DNA methylation level was higher in the testis than in the ovary. Real-time quantitative PCR (qPCR) results indicated that maintenance DNA methyltransferase gene dnmt1 and de novo DNA methyltransferase gene dnmt3a are highly expressed in the ovary, while DNA demethyltransferase genes tets are highly expressed in the testis. The inconsistency of DNA methylation and methyltransferase genes in the gonads might associate with the differential distribution in the testis. 5-mC mainly located in the spermatids of the testis was shown with immunohistochemistry (IHC). Furtherly, dnmt3a and tets are mainly located in spermatocytes and oocytes with in situ hybridization (ISH) analysis. As for H3K4me3, total level is higher in the ovary detected with western blot assay. IHC results showed that the signals of H3K4me3 in Sertoli cells of the testis were stronger than those in spermatocytes and spermatids. Methyltransferase gene kmt2b and demethylase genes kdm5a and kdm5c also exhibit much higher expression in the testis with qPCR, and ISH stronger signals of kmt2b and kdm5s were detected in spermatocytes. These results implied that DNA methylation and H3K4me3 are involved in the flounder sex differences and gametogenesis.

The DNA cytosine-5-methyltransferase 3 (DNMT3) involved in regulation of CgIL-17 expression in the immune response of oyster Crassostrea gigas

DNA methyltransferase, a key enzyme mediating DNA methylation, is involved in numerous processes including genomic imprinting, X chromosome inactivation, transposable element suppression, and immune defense in vertebrates. In the present study, a DNA cytosine-5-methyltransferase 3 was identified from oyster Crassostrea gigas (designed as CgDNMT3). There were a PWWP domain, a PHD domain and a DNA-methylase domain in the deduced amino acid sequences of CgDNMT3, and the conserved motifs I, IV, VI, Ⅷ, IX and X were identified in its C-terminal catalytic DNA-methylase domain. The mRNA transcripts of CgDNMT3 were detected in haemocytes, mantle, gill, adductor muscle, digestive gland and labial palp, with higher expression level in haemocytes (6.54 folds of those in gill, p < 0.01). The expression level of CgDNMT3 mRNA in haemocytes increased significantly after LPS primed (2.87 folds of that in control group, p < 0.05) in vitro or Vibrio splendidus challenging (1.94 folds of that in control group, p < 0.05) in vivo. Immunocytochemical analysis revealed that CgDNMT3 protein was distributed mainly in cytoplasm and partial in nucleus of oyster haemocytes. After CgDNMT3 was transfected and expressed in HEK293T cells, the DNA 5-methylcytosine (5-mc) level in the transfected group was significantly increased, which was 1.22 folds (p < 0.05) of the pcDNA-3.1 group. The expressions of oyster CgIL17-1, CgIL17-2 and CgIL17-5 in haemocytes increased (13.05 folds, 4.78 folds and 9.41 folds of that in control group, respectively) at 12 h after V. splendidus challenging, but the increase were significantly inhibited when the oysters were pre-treated with DNA methyltransferase inhibitor 5-Azacytidine, which were 9 folds, 1.93 folds and 3.22 folds of that in control group, respectively. These results collectively suggested that CgDNMT3 was a conserved member of DNA methyltransferase 3 family in oyster, and participated in regulating the expression of cytokines during immune response.

Dnmt3aa but Not Dnmt3ab Is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus)

Dnmt3a, a de novo methyltransferase, is essential for mammalian germ line DNA methylation. Only one Dnmt3a is identified in mammals, and homozygous mutants of Dnmt3a are lethal, while two Dnmt3a paralogs, dnmt3aa and dnmt3ab, are identified in teleosts due to the third round of genome duplication, and homozygous mutants of dnmt3aa and dnmt3ab are viable in zebrafish. The expression patterns and roles of dnmt3aa and dnmt3ab in gonadal development remain poorly understood in teleosts. In this study, we elucidated the precise expression patterns of dnmt3aa and dnmt3ab in tilapia gonads. Dnmt3aa was highly expressed in oogonia, phase I and II oocytes and granulosa cells in ovaries and spermatogonia and spermatocytes in testes, while dnmt3ab was mainly expressed in ovarian granulosa cells and testicular spermatocytes. The mutation of dnmt3aa and dnmt3ab was achieved by CRISPR/Cas9 in tilapia. Lower gonadosomatic index (GSI), increased apoptosis of oocytes and spermatocytes and significantly reduced sperm quality were observed in dnmt3aa^−/− mutants, while normal gonadal development was observed in dnmt3ab^−/− mutants. Consistently, the expression of apoptotic genes was significantly increased in dnmt3aa^−/− mutants. In addition, the 5-methylcytosine (5-mC) level in dnmt3aa^−/− gonads was decreased significantly, compared with that of dnmt3ab^−/− and wild type (WT) gonads. Taken together, our results suggest that dnmt3aa, not dnmt3ab, plays important roles in maintaining gametogenesis in teleosts.

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

Studying environmental sex determination (ESD) in cichlids provides a phylogenetic and comparative approach to understand the evolution of the underlying mechanisms, their impact on the evolution of the overlying systems, and the neuroethology of life history strategies. Natural selection normally favors parents who invest equally in the development of male and female offspring, but evolution may favor deviations from this 50:50 ratio when environmental conditions produce an advantage for doing so. Many species of cichlids demonstrate ESD in response to water chemistry (temperature, pH, and oxygen concentration). The relative strengths of and the exact interactions between these factors vary between congeners, demonstrating genetic variation in sensitivity. The presence of sizable proportions of the less common sex towards the environmental extremes in most species strongly suggests the presence of some genetic sex-determining loci acting in parallel with the ESD factors. Sex determination and differentiation in these species does not seem to result in the organization of a final and irreversible sexual fate, so much as a life-long ongoing battle between competing male- and female-determining genetic and hormonal networks governed by epigenetic factors. We discuss what is and is not known about the epigenetic mechanism behind the differentiation of both gonads and sex differences in the brain. Beyond the well-studied tilapia species, the 2 best-studied dwarf cichlid systems showing ESD are the South American genus Apistogramma and the West African genus Pelvicachromis. Both species demonstrate male morphs with alternative reproductive tactics. We discuss the further neuroethology opportunities such systems provide to the study of epigenetics of alternative life history strategies and other behavioral variation.

Genome-wide identification, evolution of Krüppel-like factors (klfs) and their expressions during GCRV challenge in grass carp (Ctenopharyngodonidella)

The Krüppel-like factors (KLFs) are a family of transcription factors containing three highly conserved tandem zinc finger structures, and each member participates in multiple physiological and pathological processes. The publication of genome sequences and the application of bioinformatics tools have led to the discovery of numerous gene families in fishes. Here, 24 klf genes were re-annotated in grass carp. Subsequently, the number of klf family members were investigated in some representative vertebrate species. Then, a series of bioinformatics analysis showed that grass carp klfs in the same subfamily had similar genome structure patterns and conserved distribution patterns of motifs, which supported their molecular evolutionary relationships. Furthermore, the mRNA expression profiles showed that 24 grass carp klfs were ubiquitously expressed in 11 different tissues, and some of them displayed tissue-enriched expression patterns. Finally, the expressions of the evolutionarily expanded klf members (klf2a, 2b, 2l, 5a, 5b, 5l, 6a, 6b, 7a, 7b, 11a, 11b, 12a, 12b, 15 and 15l) during GCRV infection were also analyzed. The results suggested that grass carp klf genes with common evolutionary sources may share functional diversity and conservation. In conclusion, this study provides preliminary clues for further researches on grass carp klf members and their underlying transcriptional regulatory mechanisms during GCRV infection.

The round goby genome provides insights into mechanisms that may facilitate biological invasions

Background

The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success.

Results

We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby’s capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions.

Conclusions

The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0731-8) contains supplementary material, which is available to authorized users.

Evolutionary history of DNA methylation related genes in chordates: new insights from multiple whole genome duplications

DNA methylation is an important epigenetic mechanism involved in many biological processes, i.e. gametogenesis and embryonic development. However, increased copy numbers of DNA methylation related genes (dnmt, tet and tdg) have been found during chordate evolution due to successive whole genome duplication (WGD) events. Their evolutionary history and phylogenetic relationships remain unclear. The present study is the first to clarify the evolutionary history of DNA methylation genes in chordates. In particular, our results highlight the fixation of several dnmt3-related genes following successive WGD throughout evolution. The rainbow trout genome offered a unique opportunity to study the early evolutionary fates of duplicated genes due to a recent round of WGD at the radiation of salmonids. Differences highlighted in transcriptional patterns of these genes during gametogenesis and ontogenesis in trout indicated that they might be subjected to sub- or neo-functionalisation after WDG. The fixation of multiple dnmt3 genes in genomes after WGD could contribute to the diversification and plastic adaptation of the teleost.