Genome-wide identification, evolution of histone lysine demethylases (KDM) genes and their expression during gonadal development in Nile tilapia.

Histone lysine demethylases (KDM) are responsible for histone demethylation and are involved in gene expression regulation. Previous studies have shown that histone lysine demethylation plays an important role in gonadal development of vertebrates. The KDM family consists of eight subfamilies, i.e., kdm1, kdm2, kdm3, kdm4, kdm5, kdm6, kdm7 and JmjC-only subfamily. In this study, 13 to 63 KDM genes in 23 representative species were identified based on the available version of genome assembly. Phylogenetic relationships, domain architecture, and synteny of these genes were comprehensively analyzed and the results suggested KDM genes probably originated from the early diverging metazoan and significantly expanded in vertebrates with multiple whole genome duplication, especially in the third-round whole genome duplication (3R-WGD) and polyploidization of teleosts. The subfamilies of kdm2, kdm3, kdm4, kdm5, kdm6 and kdm7 were duplicated with 1R-2R events, and duplicates of kdm2a, kdm4a, kdm5b and kdm6b were resulted from 3R-WGD. Based on transcriptome data, the KDM genes were found to be dominantly expressed in the ovary and testis. More than 80% of KDM genes displayed sexual dimorphic expression, with 15 genes dominantly expressed in ovaries, and 12 genes dominantly expressed in testes. Importantly, from transcriptome data, qRT-PCR and fluorescence in situ hybridization during sex reversal, genes with higher expression in ovary than testis, such as kdm1b and two JmjC-only subfamily members hspbap1 and riox1, were downregulated, while other genes, such as kdm3c, kdm5bb, kdm6ba, kdm6bb and kdm7b, with higher expression in testis than ovary, were upregulated in ovotestis, indicating these genes play critical roles in the gonadal development and sex reversal. This study provided new insights into the evolution of the KDM genes and a fundamental clue for understanding their important roles in sex differentiation and gonadal development in teleosts.

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.

Chromosome-level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing, Bionano and Hi-C technology.

Onychostoma macrolepis is an emerging commercial cyprinid fish species. It is a model system for studies of sexual dimorphism and genome evolution. Here, we report the chromosome-level assembly of the O.macrolepis genome obtained from the integration of nanopore long-read sequencing with physical maps produced using Bionano and Hi-C technology. A total of 87.9 Gb of nanopore sequence provided approximately 100-fold coverage of the genome. The preliminary genome assembly was 883.2 Mb in size with a contig N50 size of 11.2 Mb. The 969 corrected contigs obtained from Bionano optical mapping were assembled into 853 scaffolds and produced an assembly of 886.5 Mb with a scaffold N50 of 16.5 Mb. Finally, using the Hi-C data, 881.3 Mb (99.4% of genome) in 526 scaffolds were anchored and oriented in 25 chromosomes ranging in size from 25.27 to 56.49 Mb. In total, 24,770 protein-coding genes were predicted in the genome, and ~96.85% of the genes were functionally annotated. The annotated assembly contains 93.3% complete genes from the BUSCO reference set. In addition, we identified 409 Mb (46.23% of the genome) of repetitive sequence, and 11,213 non-coding RNAs, in the genome. Evolutionary analysis revealed that O. macrolepis diverged from common carp approximately 24.25 million years ago. The chromosomes of O. macrolepis showed an unambiguous correspondence to the chromosomes of zebrafish. The high-quality genome assembled in this work provides a valuable genomic resource for further biological and evolutionary studies of O. macrolepis.