EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database.

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster.

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Species-specific KRAB-ZFPs function as repressors of retroviruses by targeting PBS regions.

SignificanceHosts often target the relatively conserved regions in rapidly mutating retroviruses to inhibit their replication. One of these regions is called a primer binding site (PBS), which has to be complementary to the host tRNA to initiate reverse transcription. By analyzing endogenous retroviral elements, we found that host cells use this sequence as a target in efforts to block the expression of viral elements. A specific type of zinc finger protein targets the PBS in a host genome, which not only inhibits the transcription of endogenous viruses but also inhibits the replication of exogenous retroviruses with the same PBS. Thus, our study sheds light on a strategy for searching for host restriction factors targeting retroviruses.

DSEATM: drug set enrichment analysis uncovering disease mechanisms by biomedical text mining.

Disease pathogenesis is always a major topic in biomedical research. With the exponential growth of biomedical information, drug effect analysis for specific phenotypes has shown great promise in uncovering disease-associated pathways. However, this method has only been applied to a limited number of drugs. Here, we extracted the data of 4634 diseases, 3671 drugs, 112 809 disease-drug associations and 81 527 drug-gene associations by text mining of 29 168 919 publications. On this basis, we proposed a 'Drug Set Enrichment Analysis by Text Mining (DSEATM)' pipeline and applied it to 3250 diseases, which outperformed the state-of-the-art method. Furthermore, diseases pathways enriched by DSEATM were similar to those obtained using the TCGA cancer RNA-seq differentially expressed genes. In addition, the drug number, which showed a remarkable positive correlation of 0.73 with the AUC, plays a determining role in the performance of DSEATM. Taken together, DSEATM is an auspicious and accurate disease research tool that offers fresh insights.

Gut microbiome in tumorigenesis and therapy of colorectal cancer.

Colorectal cancer (CRC) is the malignant tumor with the highest incidence in the digestive system, and the gut microbiome plays a crucial role in CRC tumorigenesis and therapy. The gastrointestinal tract is the organ harboring most of the microbiota in humans. Changes in the gut microbiome in CRC patients suggest possible host-microbe interactions, thereby hinting the potential tumorigenesis, which provides new perspective for preventing, diagnosing, or treating CRC. In this review, we discuss the effects of gut microbiome dysbiosis on CRC, and reveal the mechanisms by which gut microbiome dysbiosis leads to CRC. Gut microbiome modulation with the aim to reverse the established gut microbial dysbiosis is a novel strategy for the prevention and treatment of CRC. In addition, this review summarizes that probiotic antagonize CRC tumorigenesis by protecting intestinal barrier function, inhibiting cancer cell proliferation, resisting oxidative stress, and enhancing host immunity. Finally, we highlight clinical applications of the gut microbiome, such as gut microbiome analysis-based biomarker screening and prediction, and microbe modulation-based CRC prevention, treatment enhancement, and treatment side effect reduction. This review provides the reference for the clinical application of gut microbiome in the prevention and treatment of CRC.

SAGD: a comprehensive sex-associated gene database from transcriptomes.

Many animal species present sex differences. Sex-associated genes (SAGs), which have female-biased or male-biased expression, have major influences on the remarkable sex differences in important traits such as growth, reproduction, disease resistance and behaviors. However, the SAGs resulting in the vast majority of phenotypic sex differences are still unknown. To provide a useful resource for the functional study of SAGs, we manually curated public RNA-seq datasets with paired female and male biological replicates from the same condition and systematically re-analyzed the datasets using standardized methods. We identified 27,793 female-biased SAGs and 64,043 male-biased SAGs from 2,828 samples of 21 species, including human, chimpanzee, macaque, mouse, rat, cow, horse, chicken, zebrafish, seven fly species and five worm species. All these data were cataloged into SAGD, a user-friendly database of SAGs (http://bioinfo.life.hust.edu.cn/SAGD) where users can browse SAGs by gene, species, drug and dataset. In SAGD, the expression, annotation, targeting drugs, homologs, ontology and related RNA-seq datasets of SAGs are provided to help researchers to explore their functions and potential applications in agriculture and human health.

pyMeSHSim: an integrative python package for biomedical named entity recognition, normalization, and comparison of MeSH terms.

BACKGROUND: Many disease causing genes have been identified through different methods, but there have been no uniform annotations of biomedical named entity (bio-NE) of the disease phenotypes of these genes yet. Furthermore, semantic similarity comparison between two bio-NE annotations has become important for data integration or system genetics analysis. RESULTS: The package pyMeSHSim recognizes bio-NEs by using MetaMap which produces Unified Medical Language System (UMLS) concepts in natural language process. To map the UMLS concepts to Medical Subject Headings (MeSH), pyMeSHSim is embedded with a house-made dataset containing the main headings (MHs), supplementary concept records (SCRs), and their relations in MeSH. Based on the dataset, pyMeSHSim implemented four information content (IC)-based algorithms and one graph-based algorithm to measure the semantic similarity between two MeSH terms. To evaluate its performance, we used pyMeSHSim to parse OMIM and GWAS phenotypes. The pyMeSHSim introduced SCRs and the curation strategy of non-MeSH-synonymous UMLS concepts, which improved the performance of pyMeSHSim in the recognition of OMIM phenotypes. In the curation of 461 GWAS phenotypes, pyMeSHSim showed recall > 0.94, precision > 0.56, and F1 > 0.70, demonstrating better performance than the state-of-the-art tools DNorm and TaggerOne in recognizing MeSH terms from short biomedical phrases. The semantic similarity in MeSH terms recognized by pyMeSHSim and the previous manual work was calculated by pyMeSHSim and another semantic analysis tool meshes, respectively. The result indicated that the correlation of semantic similarity analysed by two tools reached as high as 0.89-0.99. CONCLUSIONS: The integrative MeSH tool pyMeSHSim embedded with the MeSH MHs and SCRs realized the bio-NE recognition, normalization, and comparison in biomedical text-mining.

Anti-masculinization induced by aromatase inhibitors in adult female zebrafish.

BACKGROUND: Early sex differentiation genes of zebrafish remain an unsolved mystery due to the difficulty to distinguish the sex of juvenile zebrafish. However, aromatase inhibitors (AIs) could direct juvenile zebrafish sex differentiation to male and even induce ovary-to-testis reversal in adult zebrafish. RESULTS: In order to determine the transcriptomic changes of sex differentiation in juvenile zebrafish and early sex-reversal in adult zebrafish, we sequenced the transcriptomes of juvenile and adult zebrafish treated with AI exemestane (EM) for 32 days, when juvenile zebrafish sex differentiation finished. EM treatment in females up-regulated the expression of genes involved in estrogen metabolic process, female gamete generation and oogenesis, including gsdf, macf1a and paqr5a, while down-regulated the expression of vitellogenin (vtg) genes, including vtg6, vtg2, vtg4, and vtg7 due to the lower level of Estradiol (E2). Furthermore, EM-juveniles showed up-regulation in genes related to cell death and apoptosis, such as bcl2l16 and anax1c, while the control-juveniles exhibited up-regulation of genes involved in positive regulation of reproductive process and oocyte differentiation such as zar1 and zpcx. Moreover, EM-females showed higher enrichment than control females in genes involved in VEGF signaling pathway, glycosaminoglycan degradation, hedgehog signaling pathway, GnRH signaling pathway and steroid hormone biosynthesis. CONCLUSIONS: Our study shows anti-masculinization in EM-treated adult females but not in EM-treated juveniles. This may be responsible for the lower sex plasticity in adults than juveniles.

[Genome survey and characteristic analysis of SSR in Callicarpa nudiflora].

Callicarpa nudiflora,which is a big brand of Li nationality medicine with Hainan characteristics,has the effects of dissolving stasis,hemostasis,anti-inflammatory and antibacterial. At present,there is a lack of information about the reference genome of C. nudiflora. The study of the genome size,heterozygosity rate and characteristics of SSR of C. nudiflora,can provide an effective basis for the formulation of the whole genome de novo sequencing strategy and development of SSR molecular markers of C. nudiflora. To realize this purpose,high throughput sequencing platform Illumina Hiseq was used to sequence the genome structure of C. nudiflora and K-mer analysis was applied to estimate genome size,repeat sequences and heterozygosity rate. Simple-sequence repeat( SSR) loci that are suitable as markers were identified by MISA software. The results showed the estimated genome size of C. nudiflora was 822. 43 Mb,with a 0. 85% heterozygosity rate and 71. 67% repeats,and the GC content of genome was about 49. 20%. Therefore,C. nudiflora belongs to a complex genome with high heterozygosity and repetition. SSR molecular genetic markers were analyzed in the genome sequence,and a total of 206 049 SSRs were identified,among which mono-nucleotide,di-nucleotide and tri-nucleotide repetitive motifs summed up to 198 993,accounting for 96. 57% of the total SSRs. Among the 2-6 nucleotide repeats,AT/AT,AAT/ATT,AGCC/CTGG,AAAAT/ATTTT and AGATAT/ATATCT have the largest number,respectively. This report represents the first genome-wide characterization of C. nudiflora,and provides a reference for the construction of the library for the fine sequencing of the genome,and a molecular basis for the development of SSR molecular markers as well as for the protection and utilization of gene resources.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Comparison of normalization and differential expression analyses using RNA-Seq data from 726 individual Drosophila melanogaster.

BACKGROUND: A generally accepted approach to the analysis of RNA-Seq read count data does not yet exist. We sequenced the mRNA of 726 individuals from the Drosophila Genetic Reference Panel in order to quantify differences in gene expression among single flies. One of our experimental goals was to identify the optimal analysis approach for the detection of differential gene expression among the factors we varied in the experiment: genotype, environment, sex, and their interactions. Here we evaluate three different filtering strategies, eight normalization methods, and two statistical approaches using our data set. We assessed differential gene expression among factors and performed a statistical power analysis using the eight biological replicates per genotype, environment, and sex in our data set. RESULTS: We found that the most critical considerations for the analysis of RNA-Seq read count data were the normalization method, underlying data distribution assumption, and numbers of biological replicates, an observation consistent with previous RNA-Seq and microarray analysis comparisons. Some common normalization methods, such as Total Count, Quantile, and RPKM normalization, did not align the data across samples. Furthermore, analyses using the Median, Quantile, and Trimmed Mean of M-values normalization methods were sensitive to the removal of low-expressed genes from the data set. Although it is robust in many types of analysis, the normal data distribution assumption produced results vastly different than the negative binomial distribution. In addition, at least three biological replicates per condition were required in order to have sufficient statistical power to detect expression differences among the three-way interaction of genotype, environment, and sex. CONCLUSIONS: The best analysis approach to our data was to normalize the read counts using the DESeq method and apply a generalized linear model assuming a negative binomial distribution using either edgeR or DESeq software. Genes having very low read counts were removed after normalizing the data and fitting it to the negative binomial distribution. We describe the results of this evaluation and include recommended analysis strategies for RNA-Seq read count data.

Variations in Essential Oil Yield, Composition, and Antioxidant Activity of Different Plant Organs from Blumea balsamifera (L.) DC. at Different Growth Times.

Blumea balsamifera, also named Ainaxiang, is widely used as an ancient medicinal herb in tropical and subtropical Asia. It is rich in essential oils. In this work the essential oils of B. balsamifera from different plant organs and in different months were extracted, and then analyzed by gas chromatography-mass spectrometry. The results showed that essential oil yield of young leaves was the highest (0.65 mL/100 g), followed by mature leaves (0.57 mL/100 g), and the oil yield was higher in October (0.47 mL/100 g) than other months. A total of 44 compounds were identified, representing 92.64%-96.71% of the oil. Eighteen common chemical components were found among the six plant organs, representing >80% of the oil constituents. l-borneol was the main ingredient in leaves, and its content was the highest in senescent leaves and in December. In the essential oils of young shoots and young stems, the main component was dimethoxydurene. Antioxidant activity was also determined using the 2,2-diphenyl-1-picrylhydrazyl (DPPH) and ß-carotene bleaching (BCB) assays. The results indicated that the ß-carotene bleaching activity was far stronger than the DPPH radical-scavenging capacity, and the young leaves and young shoots showed stronger antioxidant activity. Dimethoxydurene, ß-caryophyllene, and α-caryophyllene play a positive role in good antioxidant activity, while ß-eudesmol, phytol, and tetradecanal play a negative role. The antioxidant activity revealed in this study might help in developing this promising bioresource for use in the medicinal and cosmetic industries.

External application of the volatile oil from Blumea balsamifera may be safe for liver--a study on its chemical composition and hepatotoxicity.

Ainaxiang (Blumea balsamifera), also known as Sambong, is an important ancient medicinal herb in Southeast Asia. It is rich in volatile oil, and still widely used nowadays for skin wound healing and treatment of sore throats. We analyzed the volatile oil from Blumea balsamifera (BB oil) by gas chromatography-mass spectrometry (GC-MS). Forty one components, including l-borneol, were identified. Next, the damaging effects of BB oil diluted with olive oil on liver at different concentrations (100%, 50%, 20%), were evaluated, using both normal and wounded skin. Plasma ALT, AST, ALP and TBili were assessed, along with liver histopathology. The results showed that serum levels of liver toxicity markers in the high concentration groups (100% w/v) increased compared with control groups, whereas no significant changes was observed in histopathology of liver samples. In the wound groups, treatment with BB oil resulted in a decrease in serum toxicity index, compared with normal animal groups. This study confirms the safety of short term BB oil consumption, though high BB oil doses may lead to mild liver injury and this response might be weakened in the case of cutaneous wounds. These results are expected to be helpful for guiding appropriate therapeutic use of BB oil.

Black rice diet alleviates colorectal cancer development through modulating tryptophan metabolism and activating AHR pathway.

Consumption of dietary fiber and anthocyanin has been linked to a lower incidence of colorectal cancer (CRC). This study scrutinizes the potential antitumorigenic attributes of a black rice diet (BRD), abundantly rich in dietary fiber and anthocyanin. Our results demonstrate notable antitumorigenic effects in mice on BRD, indicated by a reduction in both the size and number of intestinal tumors and a consequent extension in life span, compared to control diet-fed counterparts. Furthermore, fecal transplants from BRD-fed mice to germ-free mice led to a decrease in colonic cell proliferation, coupled with maintained integrity of the intestinal barrier. The BRD was associated with significant shifts in gut microbiota composition, specifically an augmentation in probiotic strains Bacteroides uniformis and Lactobacillus. Noteworthy changes in gut metabolites were also documented, including the upregulation of indole-3-lactic acid and indole. These metabolites have been identified to stimulate the intestinal aryl hydrocarbon receptor pathway, inhibiting CRC cell proliferation and colorectal tumorigenesis. In summary, these findings propose that a BRD may modulate the progression of intestinal tumors by fostering protective gut microbiota and metabolite profiles. The study accentuates the potential health advantages of whole-grain foods, emphasizing the potential utility of black rice in promoting health.

In vivo dynamics and regulation of DNA G-quadruplex structures in mammals.

G-quadruplex (G4) is a four-stranded helical DNA secondary structure formed by guanine-rich sequence folding, and G4 has been computationally predicted to exist in a wide range of species. Substantial evidence has supported the formation of endogenous G4 (eG4) in living cells and revealed its regulatory dynamics and critical roles in several important biological processes, making eG4 a regulator of gene expression perturbation and a promising therapeutic target in disease biology. Here, we reviewed the methods for prediction of potential G4 sequences (PQS) and detection of eG4s. We also highlighted the factors affecting the dynamics of eG4s and the effects of eG4 dynamics. Finally, we discussed the future applications of eG4 dynamics in disease therapy.

Male-Biased Gut Microbiome and Metabolites Aggravate Colorectal Cancer Development.

Men demonstrate higher incidence and mortality rates of colorectal cancer (CRC) than women. This study aims to explain the potential causes of such sexual dimorphism in CRC from the perspective of sex-biased gut microbiota and metabolites. The results show that sexual dimorphism in colorectal tumorigenesis is observed in both ApcMin/ + mice and azoxymethane (AOM)/dextran sulfate sodium (DSS)-treated mice with male mice have significantly larger and more tumors, accompanied by more impaired gut barrier function. Moreover, pseudo-germ mice receiving fecal samples from male mice or patients show more severe intestinal barrier damage and higher level of inflammation. A significant change in gut microbiota composition is found with increased pathogenic bacteria Akkermansia muciniphila and deplets probiotic Parabacteroides goldsteinii in both male mice and pseudo-germ mice receiving fecal sample from male mice. Sex-biased gut metabolites in pseudo-germ mice receiving fecal sample from CRC patients or CRC mice contribute to sex dimorphism in CRC tumorigenesis through glycerophospholipids metabolism pathway. Sexual dimorphism in tumorigenesis of CRC mouse models. In conclusion, the sex-biased gut microbiome and metabolites contribute to sexual dimorphism in CRC. Modulating sex-biased gut microbiota and metabolites could be a potential sex-targeting therapeutic strategy of CRC.

Deficiency of X-linked inverted duplicates with male-biased expression and the underlying evolutionary mechanisms in the Drosophila genome.

Inverted duplicates (IDs) are pervasive in genomes and have been reported to play functional roles in various biological processes. However, the general underlying evolutionary forces that maintain IDs in genomes remain largely elusive. Through a systematic screening of the Drosophila melanogaster genome, 20,223 IDs were detected in nonrepetitive intergenic regions, far more than expectation under the neutrality model. 3,846 of these IDs were identified to have stable hairpin structure (i.e., the structural IDs). Based on whole-genome transcriptome profiling data, we found 628 unannotated expressed structural IDs, which had significantly different genomic distributions and structural properties from the unexpressed IDs. Among the expressed structural IDs, 130 exhibited higher expression in males than in females (i.e., male-biased expression). Compared with sex-unbiased ones, these male-biased IDs were significantly underrepresented on the X chromosome, similar to previously reported pattern of male-biased protein-coding genes. These analyses suggest that a selection-driven process, rather than a purely neutral mutation-driven mechanism, contributes to the maintenance of IDs in the Drosophila genome.

Evolution and function of developmentally dynamic pseudogenes in mammals.

BACKGROUND: Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored. RESULTS: To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation. CONCLUSIONS: Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

A Germline-Specific Regulator of Mitochondrial Fusion is Required for Maintenance and Differentiation of Germline Stem and Progenitor Cells.

Maintenance and differentiation of germline stem and progenitor cells (GSPCs) is important for sexual reproduction. Here, the authors identify zebrafish pld6 as a novel germline-specific gene by cross-analyzing different RNA sequencing results, and find that pld6 knockout mutants develop exclusively into infertile males. In pld6 mutants, GSPCs fail to differentiate and undergo apoptosis, leading to masculinization and infertility. Mitochondrial fusion in pld6-depleted GSPCs is severely impaired, and the mutants exhibit defects in piRNA biogenesis and transposon suppression. Overall, this work uncovers zebrafish Pld6 as a novel germline-specific regulator of mitochondrial fusion, and highlights its essential role in the maintenance and differentiation of GSPCs as well as gonadal development and gametogenesis.

Sexually Dimorphic Gene Expression in X and Y Sperms Instructs Sexual Dimorphism of Embryonic Genome Activation in Yellow Catfish (Pelteobagrus fulvidraco).

Paternal factors play an important role in embryonic morphogenesis and contribute to sexual dimorphism in development. To assess the effect of paternal DNA on sexual dimorphism of embryonic genome activation, we compared X and Y sperm and different sexes of embryos before sex determination. Through transcriptome sequencing (RNA-seq) and whole-genome bisulfite sequencing (WGBS) of X and Y sperm, we found a big proportion of upregulated genes in Y sperm, supported by the observation that genome-wide DNA methylation level is slightly lower than in X sperm. Cytokine-cytokine receptor interaction, TGF-beta, and toll-like receptor pathways play important roles in spermatogenesis. Through whole-genome re-sequencing (WGRS) of parental fish and RNA-seq of five early embryonic stages, we found the low-blastocyst time point is a key to maternal transcriptome degradation and zygotic genome activation. Generally, sexual differences emerged from the bud stage. Moreover, through integrated analysis of paternal SNPs and gene expression, we evaluated the influence of paternal inheritance on sexual dimorphism of genome activation. Besides, we screened out gata6 and ddx5 as potential instructors for early sex determination and gonad development in yellow catfish. This work is meaningful for revealing the molecular mechanisms of sex determination and sexual dimorphism of fish species.