Transcriptome and DNA Methylation Responses in the Liver of Yellowfin Seabream Under Starvation Stress.

Fish suffer from starvation due to environmental risks such as extreme weather in the wild and due to insufficient feedings in farms. Nutrient problems from short-term or long-term starvation conditions can result in stress-related health problems for fish. Yellowfin seabream (Acanthopagrus latus) is an important marine economic fish in China. Understanding the molecular responses to starvation stress is vital for propagation and culturing yellowfin seabream. In this study, the transcriptome and genome-wide DNA methylation levels in the livers of yellowfin seabream under 14-days starvation stress were analyzed. One hundred sixty differentially expressed genes (DEGs) by RNA-Seq analysis and 737 differentially methylated-related genes by whole genome bisulfite sequencing analysis were identified. GO and KEGG pathway enrichment analysis found that energy metabolism-related pathways such as glucose metabolism and lipid metabolism were in response to starvation. Using bisulfite sequencing PCR, we confirmed the presence of CpG methylation differences within the regulatory region of a DEG ppargc1a in response to 14-days starvation stress. This study revealed the molecular responses of livers in response to starvation stress at the transcriptomic and whole genome DNA methylation levels in yellowfin seabream.

First identification of two co-existing genome-wide significant sex quantitative trait loci (QTL) in red tilapia using integrative QTL mapping.

Red tilapia ( Oreochromis spp .) is one of the most popular fish in China due to its bright red appearance, fast growth rate, and strong adaptability. Understanding the sex determination mechanisms is of vital importance for the selection of all-male lines to increase aquacultural production of red tilapia. In this research, the genetic architecture for sex from four mapping populations ( n=1 090) of red tilapia was analyzed by quantitative trait loci (QTL)-seq, linkage-based QTL mapping, and linkage disequilibrium (LD)-based genome-wide association studies. Two genome-wide significant QTL intervals associated with sex were identified on ChrLG1 (22.4-23.9 Mb) and ChrLG23 (32.0-35.9 Mb), respectively. The QTL on ChrLG1 was detected in family 1 (FAM1), FAM2, and FAM4, and the other QTL on ChrLG23 was detected in FAM3 and FAM4. Four microsatellite markers located within the QTL were successfully developed for marker-assisted selection. Interestingly, three ( lpp, sox14, and amh) of the 12 candidate genes located near or on the two QTL intervals were abundantly expressed in males, while the remaining genes were more highly expressed in females. Seven genes ( scly, ube3a, lpp, gpr17, oca2, cog4, and atp10a) were significantly differentially expressed between the male and female groups. Furthermore, LD block analysis suggested that a cluster of genes on ChrLG23 may participate in regulating sex development in red tilapia. Our study provides important information on the genetic architecture of sex in red tilapia and should facilitate further exploration of sex determination mechanisms in this species.

[Biosynthesis and metabolism regulation of terpenoids in Pogostemon cablin: a review].

The terpenoids in Pogostemon cablin have complex structures and abundant pharmacological effects. Patchouli alcohol(PA) and pogostone(PO) have a high medicinal value by virtue of anti-tumor, anti-inflammatory, antibacterial, antioxidant, and other biological activities. Due to the low content of terpenoid metabolites in P. cablin, the study of biosynthesis and metabolism regulation can provide a biosynthetic basis for obtaining high-content terpenoids. In this study, key enzyme genes in biosynthesis, transcription factors in metabolism regulation, spatio-temporal expression of terpene synthase were reviewed, aiming to provide a reference for the development, protection, and utilization of P. cablin resources.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia.

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Whole-genome sequencing of leopard coral grouper ( Plectropomus leopardus) and exploration of regulation mechanism of skin color and adaptive evolution.

Leopard coral groupers belong to the Plectropomus genus of the Epinephelidae family and are important fish for coral reef ecosystems and the marine aquaculture industry. To promote future research of this species, a high-quality chromosome-level genome was assembled using PacBio sequencing and Hi-C technology. A 787.06 Mb genome was assembled, with 99.7% (784.57 Mb) of bases anchored to 24 chromosomes. The leopard coral grouper genome size was smaller than that of other groupers, which may be related to its ancient status among grouper species. A total of 22 317 protein-coding genes were predicted. This high-quality genome of the leopard coral grouper is the first genomic resource for Plectropomus and should provide a pivotal genetic foundation for further research. Phylogenetic analysis of the leopard coral grouper and 12 other fish species showed that this fish is closely related to the brown-marbled grouper. Expanded genes in the leopard coral grouper genome were mainly associated with immune response and movement ability, which may be related to the adaptive evolution of this species to its habitat. In addition, we also identified differentially expressed genes (DEGs) associated with carotenoid metabolism between red and brown-colored leopard coral groupers. These genes may play roles in skin color decision by regulating carotenoid content in these groupers.

Differential Transcriptomic and Metabolomic Responses in the Liver of Nile Tilapia (Oreochromis niloticus) Exposed to Acute Ammonia.

Ammonia is toxic to aquatic animal. Currently, only limited works were reported on the responses of aquatic animals after ammonia exposure using "omics" technologies. Tilapia suffers from the stress of ammonia-nitrogen during intensive recirculating aquaculture. Optimizing ammonia stress tolerance has become an important issue in tilapia breeding. The molecular and biochemical mechanisms of ammonia-nitrogen toxicity have not been understood comprehensively in tilapia yet. In this study, using RNA-seq and gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS) techniques, we investigated differential expressed genes (DEGs) and metabolomes in the liver at 6 h post-challenges (6 hpc) and 24 h post-challenges (24 hpc) under high concentration of ammonia-nitrogen treatment. We detected 2258 DEGs at 6 hpc and 315 DEGs at 24 hpc. Functional enrichment analysis indicated that DEGs were significantly associated with cholesterol biosynthesis, steroid and lipid metabolism, energy conservation, and mitochondrial tissue organization. Metabolomic analysis detected 31 and 36 metabolites showing significant responses to ammonia-nitrogen stress at 6 and 24 hpc, respectively. D-(Glycerol 1-phosphate), fumaric acid, and L-malic acid were found significantly down-regulated at both 6 and 24 hpc. The integrative analysis of transcriptomics and metabolomics suggested considerable alterations and precise control of gene expression at both physiological and molecular levels in response to the stress of ammonia-nitrogen in tilapia.

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.).

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

Effect of fasting and subsequent refeeding on the transcriptional profiles of brain in juvenile Spinibarbus hollandi.

Starvation is a common stress in fish. The underlying molecular mechanisms associated with growth depression caused by feeding restriction and compensatory growth are not well understood. We investigated the effect of fasting and refeeding on the transcriptome profiles of brain in juvenile S. hollandi using RNA-seq. A total of 4.73 × 108 raw reads were obtained from nine brain samples. De novo transcriptome assembly identified 387,085 unigenes with 2.1×109 nucleotides. A total of 936 annotated unigenes showed significantly differential expression among the control, fasting, and fasting-refeeding groups. The down-regulated differentially expressed genes (DEGs) during fasting were mainly associated with cell cycle, DNA replication, and mitosis. The up-regulated DEGs were mainly related to glucose and lipid metabolism, material transportation, and transcription factors. Most decreased DEGs during fasting were restored to normal levels after refeeding. Comparing with the control group, genes associated with protein synthesis, stimulus response, and carbohydrate metabolism were significantly over-expressed and pro-opio melanocortin (POMC) was down-regulated during the refeeding period. In conclusion, fish mobilized stored energetic materials and reduced energy consumption to prolong survival during fasting. After refeeding, the down-regulation of DEGs, e.g., POMC may be associated with compensatory growth. Up-regulation of DEGs related to ribosomal protein, stimulus response, and carbohydrate metabolism may contribute to eliminate negative effect of starvation on brain. This study provided the first transcriptome data related with impact of short-time starvation and refeeding in S. hollandi brains.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq.

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

Marker-assisted selection of YY supermales from a genetically improved farmed tilapia-derived strain.

Genetically improved farmed tilapia (GIFT) and GIFT-derived strains account for the majority of farmed tilapia worldwide. As male tilapias grow much faster than females, they are often considered more desirable in the aquacultural industry. Sex reversal of females to males using the male sex hormone 17-α-methyltestosterone (MT) is generally used to induce phenotypic males during large-scale production of all male fingerlings. However, the widespread use of large quantities of sex reversal hormone in hatcheries may pose a health risk to workers and ecological threats to surrounding environments. Breeding procedures to produce genetically all-male tilapia with limited or no use of sex hormones are therefore urgently needed. In this study, by applying marker-assisted selection (MAS) for the selection of YY supermales from a GIFT-derived strain, we identified 24 XY pseudofemale and 431 YY supermale tilapias. Further performance evaluation on the progenies of the YY supermales resulted in male rates of 94.1%, 99.5% and 99.6%, respectively, in three populations, and a daily increase in body weight of 1.4 g at 3 months (n=997). Our study established a highly effective MAS procedure in the selection of YY supermales from a GIFT-derived strain. Furthermore, the development of MAS-selected YY supermales will help reduce the utilization of hormones for controlling sex in the tilapia aquaculture.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia.

BACKGROUND: Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited. RESULTS: Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = - 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = - 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs. CONCLUSIONS: This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

The FTO Gene Is Associated with Growth and Omega-3/-6 Ratio in Asian Seabass.

Polymorphisms in the FTO gene are associated with obesity and body mass index in humans and livestock. Little information of whether FTO plays an important role in aquaculture fish species is available. We cloned and characterized the FTO gene in an economically important food fish species: Asian seabass (Lates calcarifer). The full-length cDNA of the gene is 3679 bp, containing an ORF of 1935 bp encoding 644 amino acids, a 216 bp 5' UTR and a 1538 bp 3' UTR. The gene consisted of nine exons and eight introns and was 117,679 bp in length. Phylogenetic analysis revealed that the gene in Asian seabass was closely related to those of Japanese flounder and Nile tilapia. Analysis of its expressions using qRT-PCR showed that it was expressed ubiquitously, but was higher in the liver, stomach and intestine. Comparative analysis of the genomic sequences of part of intron 1 of the gene among 10 unrelated individuals identified two SNPs. Analysis of associations between SNPs and traits (i.e. growth, oil content, omega-3 and -6 contents) in an F2 family demonstrated that the two SNPs were significantly associated with growth, oil content, omega-3 content and omega-3/-6 ratio. Altogether, our data suggest that the gene or/and its linked genes play an important role in growth and fatty acid synthesis, and that the SNPs associated with traits may be used as markers for selecting quicker growth and higher omega-3/-6 ratio at the fingerling stage.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq.

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Identifying selectively important amino acid positions associated with alternative habitat environments in fish mitochondrial genomes.

Fish species inhabitating seawater (SW) or freshwater (FW) habitats have to develop genetic adaptations to alternative environment factors, especially salinity. Functional consequences of the protein variations associated with habitat environments in fish mitochondrial genomes have not yet received much attention. We analyzed 829 complete fish mitochondrial genomes and compared the amino acid differences of 13 mitochondrial protein families between FW and SW fish groups. We identified 47 specificity determining sites (SDS) that associated with FW or SW environments from 12 mitochondrial protein families. Thirty-two (68%) of the SDS sites are hydrophobic, 13 (28%) are neutral, and the remaining sites are acidic or basic. Seven of those SDS from ND1, ND2 and ND5 were scored as probably damaging to the protein structures. Furthermore, phylogenetic tree based Bayes Empirical Bayes analysis also detected 63 positive sites associated with alternative habitat environments across ten mtDNA proteins. These signatures could be important for studying mitochondrial genetic variation relevant to fish physiology and ecology.

Acute hypoxia stress induced abundant differential expression genes and alternative splicing events in heart of tilapia.

Hypoxia is one of the critical environmental stressors for fish in aquatic environments. Although accumulating evidences indicate that gene expression is regulated by hypoxia stress in fish, how genes undergoing differential gene expression and/or alternative splicing (AS) in response to hypoxia stress in heart are not well understood. Using RNA-seq, we surveyed and detected 289 differential expressed genes (DEG) and 103 genes that undergo differential usage of exons and splice junctions events (DUES) in heart of a hypoxia tolerant fish, Nile tilapia, Oreochromis niloticus following 12h hypoxic treatment. The spatio-temporal expression analysis validated the significant association of differential exon usages in two randomly selected DUES genes (fam162a and ndrg2) in 5 tissues (heart, liver, brain, gill and spleen) sampled at three time points (6h, 12h, and 24h) under acute hypoxia treatment. Functional analysis significantly associated the differential expressed genes with the categories related to energy conservation, protein synthesis and immune response. Different enrichment categories were found between the DEG and DUES dataset. The Isomerase activity, Oxidoreductase activity, Glycolysis and Oxidative stress process were significantly enriched for the DEG gene dataset, but the Structural constituent of ribosome and Structural molecule activity, Ribosomal protein and RNA binding protein were significantly enriched only for the DUES genes. Our comparative transcriptomic analysis reveals abundant stress responsive genes and their differential regulation function in the heart tissues of Nile tilapia under acute hypoxia stress. Our findings will facilitate future investigation on transcriptome complexity and AS regulation during hypoxia stress in fish.

Differential Gene Expression Profiles and Alternative Isoform Regulations in Gill of Nile Tilapia in Response to Acute Hypoxia.

Fish often encounters exposures to acute environmental hypoxia either spatially or temporally. Gill organ plays important roles in response to hypoxic stress in fish. Few studies focus on the molecular regulation mechanisms of gills under hypoxic stress. In this study, we investigated the transcriptomic response to 12-h acute hypoxia in gill of a hypoxia tolerant fish, Nile tilapia Oreochromis niloticus through RNA sequencing (RNA-Seq). We sequenced messenger RNA from three control samples and three hypoxia-treated samples. Bioinformatics analysis identified 239 differentially expressed genes (DEG) and 34 genes (DUES) that had significant differential alternative isoform regulation events in at least one exonic region in gill in response to acute hypoxia. The spatiotemporal expression analysis in five tissues (heart, liver, brain, gill, and spleen) sampled at three time points (6, 12, and 24 h) under hypoxia treatment confirmed the significant association of differential exon usages in two DUES genes (TLDC2 and SSX2IPA) with hypoxia conditions. Further functional analysis suggested several energy and immune response-related pathways, e.g., metabolic pathway and antigen processing and presentation, contained the most abundant DEG genes. We found that some GO biological processes for DEG genes were significantly enriched under hypoxic stress, such as glycolysis, metabolic process, generation of precursor metabolites and energy, and cholesterol metabolic process. Our findings suggest abundant differential gene expression changes and alternative isoform regulation events in genes involved in the hypoxia response in gill. Our results provide a basis for exploring the gene regulation mechanism under hypoxic stress in fish.

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia.

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6-14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.

Characterization and functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in tilapia.

Hypoxia is a major cause of fish morbidity and mortality in the aquatic environment. Hypoxia-inducible factors are very important modulators in the transcriptional response to hypoxic stress. In this study, we characterized and conducted functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in Nile tilapia (Oreochromis niloticus). By cloning and Sanger sequencing, we obtained the full length cDNA sequences for HIF1α (2686bp) and HIF1αn (1308bp), respectively. The CDS of HIF1α includes 15 exons encoding 768 amino acid residues and the CDS of HIF1αn contains 8 exons encoding 354 amino acid residues. The complete CDS sequences of HIF1α and HIF1αn cloned from tilapia shared very high homology with known genes from other fishes. HIF1α show differentiated expression in different tissues (brain, heart, gill, spleen, liver) and at different hypoxia exposure times (6h, 12h, 24h). HIF1αn expression level under hypoxia is generally increased (6h, 12h, 24h) and shows extremely highly upregulation in brain tissue under hypoxia. A functional determination site analysis in the protein sequences between fish and land animals identified 21 amino acid sites in HIF1α and 2 sites in HIF1αn as significantly associated sites (α = 0.05). Phylogenetic tree-based positive selection analysis suggested 22 sites in HIF1α as positively selected sites with a p-value of at least 95% for fish lineages compared to the land animals. Our study could be important for clarifying the mechanism of fish adaptation to aquatic hypoxia environment.

Copy Number Variations in Tilapia Genomes.

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R 2 > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.

Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia.

Sexual dimorphism is an interesting biological phenomenon. Previous studies showed that DNA methylation might play a role in sexual dimorphism. However, the overall picture of the genome-wide methylation landscape in sexually dimorphic species remains unclear. We analyzed the DNA methylation landscape and transcriptome in hybrid tilapia (Oreochromis spp.) using whole genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq). We found 4,757 sexually dimorphic differentially methylated regions (DMRs), with significant clusters of DMRs located on chromosomal regions associated with sex determination. CpG methylation in promoter regions was negatively correlated with the gene expression level. MAPK/ERK pathway was upregulated in male tilapia. We also inferred active cis-regulatory regions (ACRs) in skeletal muscle tissues from WGBS datasets, revealing sexually dimorphic cis-regulatory regions. These results suggest that DNA methylation contribute to sex-specific phenotypes and serve as resources for further investigation to analyze the functions of these regions and their contributions towards sexual dimorphisms.