Integration of ATAC-seq and RNA-seq Unravels Chromatin Accessibility during Sex Reversal in Orange-Spotted Grouper (Epinephelus coioides)

Comprehensive transcriptional regulatory networks in potato through chromatin accessibility and transcriptome under drought and salt stresses

Drought and high salt stress have a great negative impact on potato growth and development. However, the molecular mechanisms by which different tissues and organs of potato plants respond to drought and high salt stress at different stress times lack definition. In this study, we mapped the whole genome of THSs in potato in response to different stresses using RNA-seq and ATAC-seq technologies and constructed the unique and shared transcriptional regulatory networks (TRNs) under different stresses, stress time points, and tissues in potato. The results showed opposite trends for changes in chromatin accessibility and expression of genes under drought and salt stresses. Forty-eight hours and root-specific TRNs were more complex than those of 3 h and leaf, and there were genes with inconsistent drought- and salt-stress-regulated expression only in root-shared TRNs, demonstrating the variability of potato's response to stresses under different tissues and treatment times. These results provide a basis for elucidating the transcriptional mechanisms underlying the specific response of potato to drought or salt stress and the common response to salt and drought stress.

Meta-analysis of gonadal transcriptome provides novel insights into sex change mechanism across protogynous fishes

Protogyny, being capable of changing from female to male during their lifetime, is prevalent in 20 families of teleosts but is believed to have evolved within specific evolutionary lineages. Therefore, shared regulatory factors governing the sex change process are expected to be conserved across protogynous fishes. However, a comprehensive understanding of this mechanism remains elusive. To identify these factors, we conducted a meta-analysis using gonadal transcriptome data from seven species. We curated data pairs of ovarian tissue and transitional gonad, and employed ratios of expression level as a unified criterion for differential expression, enabling a meta-analysis across species. Our approach revealed that classical sex change-related genes exhibited differential expression levels between the ovary and transitional gonads, consistent with previous reports. These results validate our methodology's robustness. Additionally, we identified novel genes not previously linked to gonadal sex change in fish. Notably, changes in the expression levels of acetoacetyl-CoA synthetase and apolipoprotein Eb, which are involved in cholesterol synthesis and transport, respectively, suggest that the levels of cholesterol, a precursor of steroid hormones crucial for sex change, are decreased upon sex change onset in the gonads. This implies a potential universal influence of cholesterol dynamics on gonadal transformation in protogyny.

Integrative ATAC-seq and RNA-seq Analysis of the Longissimus Dorsi Muscle of Gannan Yak and Jeryak

Jeryak is the F1 generation of the cross between Gannan yak and Jersey cattle, which has the advantages of fast growth and high adaptability. The growth and development of skeletal muscle is closely linked to meat production and the quality of meat. However, the molecular regulatory mechanisms of muscle growth differences between Gannan yak and Jeryak analyzed from the perspective of chromatin opening have not been reported. In this study, ATAC-seq was used to analyze the difference of chromatin openness in longissimus muscle of Gannan yak and Jeryak. It was found that chromatin accessibility was more enriched in Jeryak compared to Gannan yak, especially in the range of the transcription start site (TSS) ± 2 kb. GO and KEGG enrichment analysis indicate that differential peak-associated genes are involved in the negative regulation of muscle adaptation and the Hippo signaling pathway. Integration analysis of ATAC-seq and RNA-seq revealed overlapping genes were significantly enriched during skeletal muscle cell differentiation and muscle organ morphogenesis. At the same time, we screened FOXO1, ZBED6, CRY2 and CFL2 for possible involvement in skeletal muscle development, constructed a genes and transcription factors network map, and found that some transcription factors (TFs), including YY1, KLF4, KLF5 and Bach1, were involved in skeletal muscle development. Overall, we have gained a comprehensive understanding of the key factors that impact skeletal muscle development in various breeds of cattle, providing new insights for future analysis of the molecular regulatory mechanisms involved in muscle growth and development.

DNA Methylation and Chromatin Accessibility Impact Subgenome Expression Dominance in the Common Carp (Cyprinus carpio)

DNA methylation and chromatin accessibility play important roles in gene expression, but their function in subgenome expression dominance remains largely unknown. We conducted comprehensive analyses of the transcriptome, DNA methylation, and chromatin accessibility in liver and muscle tissues of allotetraploid common carp, aiming to reveal the function of epigenetic modifications in subgenome expression dominance. A noteworthy overlap in differential expressed genes (DEGs) as well as their functions was observed across the two subgenomes. In the promoter and gene body, the DNA methylation level of the B subgenome was significantly different than that of the A subgenome. Nevertheless, differences in DNA methylation did not align with changes in homoeologous biased expression across liver and muscle tissues. Moreover, the B subgenome exhibited a higher prevalence of open chromatin regions and greater chromatin accessibility, in comparison to the A subgenome. The expression levels of genes located proximally to open chromatin regions were significantly higher than others. Genes with higher chromatin accessibility in the B subgenome exhibited significantly elevated expression levels compared to the A subgenome. Contrastingly, genes without accessibility exhibited similar expression levels in both subgenomes. This study contributes to understanding the regulation of subgenome expression dominance in allotetraploid common carp.

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

Characterization of chromatin accessibility and gene expression reveal the key genes involved in cotton fiber elongation

Cotton (Gossypium hirsutum L.) is an important economic crop, and cotton fiber is one of the longest plant cells, which provides an ideal model for the study of cell elongation and secondary cell wall synthesis. Cotton fiber length is regulated by a variety of transcription factors (TF) and their target genes; however, the mechanism of fiber elongation mediated by transcriptional regulatory networks is still unclear to a large extent. Here, we used a comparative assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) assay and RNA-seq analysis to identify fiber elongation transcription factors and genes using the short-fiber mutant ligon linless-2 (Li2 ) and wild type (WT). A total of 499 differential target genes were identified and GO analysis shows that differential genes are mainly involved in plant secondary wall synthesis and microtubule-binding processes. Analysis of the genomic regions preferentially accessible (Peak) has identified a number of overrepresented TF-binding motifs, highlighting sets of TFs that are important for cotton fiber development. Using ATAC-seq and RNA-seq data, we have constructed a functional regulatory network of each TF regulatory target gene and also the network pattern of TF regulating differential target genes. Further, to obtain the genes related to fiber length, the differential target genes were combined with FLGWAS data to identify the genes highly related to fiber length. Our work provides new insights into cotton fiber elongation.

Expression pattern analysis of anti-Mullerian hormone in testis development of pearlscale angelfish (Centropyge vrolikii)

In vertebrates, anti-Mullerian hormone (Amh) secreted by Sertoli cells (SC) performs a pivotal function in male sex differentiation. Compared with that of higher vertebrates, the expression pattern of Amh is more diversified in fish. In this study, the full-length complementary DNA (cDNA) of Amh in Centropyge vrolikii (Cv-Amh) was cloned and analysed, which was 2,470 bp, including a 238 bp 5'UTR, a 1,602 bp ORF and a 633 bp 3'UTR; the similarity of Amh between Cv-Amh and other fish is relatively high. The quantitative real-time PCR (qRT-PCR) results of healthy tissues and gonads at sex reversal stages in C. vrolikii showed that the expression level of Amh in the testis was significantly higher than that in other tissues (P < 0.05). Amh was weakly expressed in the vitellogenic stage ovary and perinucleolus stage ovary, but its expression significantly increased in the gonads at the hermaphroditic stage, and finally reached the highest in the pure testis after sexual reversal. The results of in situ hybridization indicated that the positive signal of Amh was strongly concentrated in SCs of testis. After Amh knockdown in the gonads, the effect on sex-related genes was tested using qRT-PCR. Among these, the expression of Dmrt1, Cyp11a, Hsd11b2, Sox8 and Sox9 significantly decreased, whereas that of Cyp19a, Sox4, Foxl2 and Sox3 increased. These results suggested that Amh could be the pivotal gene in reproductive regulation in C. vrolikii, and the data will contribute to sex-related research of C. vrolikii in the future.

ATAC-seq and RNA-seq analysis unravel the mechanism of sex differentiation and infertility in sex reversal chicken

BACKGROUND

Sex determination and differentiation are complex and delicate processes. In female chickens, the process of sex differentiation is sensitive and prone to be affected by the administration of aromatase inhibitors, which result in chicken sex reversal and infertility. However, the molecular mechanisms underlying sex differentiation and infertility in chicken sex reversal remain unclear. Therefore, we established a sex-reversed chicken flock by injecting an aromatase inhibitor, fadrozole, and constructed relatively high-resolution profiles of the gene expression and chromatin accessibility of embryonic gonads.

RESULTS

We revealed that fadrozole affected the transcriptional activities of several genes, such as DMRT1, SOX9, FOXL2, and CYP19A1, related to sex determination and differentiation, and the expression of a set of gonadal development-related genes, such as FGFR3 and TOX3, by regulating nearby open chromatin regions in sex-reversed chicken embryos. After sexual maturity, the sex-reversed chickens were confirmed to be infertile, and the possible causes of this infertility were further investigated. We found that the structure of the gonads and sperm were greatly deformed, and we identified several promising genes related to spermatogenesis and infertility, such as SPEF2, DNAI1, and TACR3, through RNA-seq.

CONCLUSIONS

This study provides clear insights into the exploration of potential molecular basis underlying sex differentiation and infertility in sex-reversed chickens and lays a foundation for further research into the sex development of birds.

Integration of ATAC-seq and RNA-seq analysis identifies key genes affecting intramuscular fat content in pigs

Meat quality is one of the most important economic traits in pig breeding and production, and intramuscular fat (IMF) content is the major factor in improving meat quality. The IMF deposition in pigs is influenced by transcriptional regulation, which is dependent on chromatin accessibility. However, how chromatin accessibility plays a regulatory role in IMF deposition in pigs has not been reported. Xidu black is a composite pig breed with excellent meat quality, which is an ideal research object of this study. In this study, we used the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) analysis to identify the accessible chromatin regions and key genes affecting IMF content in Xidu black pig breed with extremely high and low IMF content. First, we identified 21,960 differential accessible chromatin peaks and 297 differentially expressed genes. The motif analysis of differential peaks revealed several potential cis-regulatory elements containing binding sites for transcription factors with potential roles in fat deposition, including Mef2c, CEBP, Fra1, and AP-1. Then, by integrating the ATAC-seq and RNA-seq analysis results, we found 47 genes in the extremely high IMF (IMF_H) group compared with the extremely low IMF (IMF_L) group. For these genes, we observed a significant positive correlation between the differential gene expression and differential ATAC-seq signal (r² = 0.42). This suggests a causative relationship between chromatin remodeling and the resulting gene expression. We identified several candidate genes (PVALB, THRSP, HOXA9, EEPD1, HOXA10, and PDE4B) that might be associated with fat deposition. Through the PPI analysis, we found that PVALB gene was the top hub gene. In addition, some pathways that might regulate fat cell differentiation and lipid metabolism, such as the PI3K-Akt signaling pathway, MAPK signaling pathway, and calcium signaling pathway, were significantly enriched in the ATAC-seq and RNA-seq analysis. To the best of our knowledge, our study is the first to use ATAC-seq and RNA-seq to examine the mechanism of IMF deposition from a new perspective. Our results provide valuable information for understanding the regulation mechanism of IMF deposition and an important foundation for improving the quality of pork.

Establishment of a Spermatogonial Stem Cell Line with Potential of Meiosis in a Hermaphroditic Fish, Epinephelus coioides

Spermatogonial stem cells (SSCs) are unique adult stem cells capable of self-renewal and differentiation into sperm. Grouper is a protogynous hermaphroditic fish farmed widely in the tropical and subtropical seas. In this study, we established an SSC line derived from adult testis of orange-spotted grouper, Epinephelus coioides. In the presence of basic fibroblast growth factor (bFGF) and leukemia inhibitory factor (LIF), the cells could be maintained with proliferation and self-renewal over 20 months and 120 passages under in vitro culture conditions. The cells exhibited strong alkaline phosphatase activity and the characteristics of SSCs with the expression of germ cell markers, including Vasa, Dazl, and Plzf, as well as the stem cell markers Nanog, Oct4, and Ssea1. Furthermore, the cultured cells could be induced by 11-ketotestosterone treatment to highly express the meiotic markers Rec8, Sycp3, and Dmc1, and produce some spherical cells, and even sperm-like cells with a tail. The findings of this study suggested that the cultured grouper SSC line would serve as an excellent tool to study the molecular mechanisms behind SSCs self-renewal and differentiation, meiosis during spermatogenesis, and sex reversal in hermaphroditic vertebrates. Moreover, this SSC line has great application value in grouper fish aquaculture, such as germ cell transplantation, genetic manipulation, and disease research.

Pou5f1 and Nanog Are Reliable Germ Cell-Specific Genes in Gonad of a Protogynous Hermaphroditic Fish, Orange-Spotted Grouper (Epinephelus coioides)

Pluripotency markers Pou5f1 and Nanog are core transcription factors regulating early embryonic development and maintaining the pluripotency and self-renewal of stem cells. Pou5f1 and Nanog also play important roles in germ cell development and gametogenesis. In this study, Pou5f1 (EcPou5f1) and Nanog (EcNanog) were cloned from orange-spotted grouper, Epinephelus coioides. The full-length cDNAs of EcPou5f1 and EcNanog were 2790 and 1820 bp, and encoded 475 and 432 amino acids, respectively. EcPou5f1 exhibited a specific expression in gonads, whereas EcNanog was expressed highly in gonads and weakly in some somatic tissues. In situ hybridization analyses showed that the mRNA signals of EcNanog and EcPou5f1 were exclusively restricted to germ cells in gonads. Likewise, immunohistofluorescence staining revealed that EcNanog protein was limited to germ cells. Moreover, both EcPou5f1 and EcNanog mRNAs were discovered to be co-localized with Vasa mRNA, a well-known germ cell maker, in male and female germ cells. These results implied that EcPou5f1 and EcNanog could be also regarded as reliable germ cell marker genes. Therefore, the findings of this study would pave the way for elucidating the mechanism whereby EcPou5f1 and EcNanog regulate germ cell development and gametogenesis in grouper fish, and even in other protogynous hermaphroditic species.

Single-Cell Atlas of Adult Testis in Protogynous Hermaphroditic Orange-Spotted Grouper, Epinephelus coioides

Spermatogenesis is a process of self-renewal and differentiation in spermatogonial stem cells. During this process, germ cells and somatic cells interact intricately to ensure long-term fertility and accurate genome propagation. Spermatogenesis has been intensely investigated in mammals but remains poorly understood with regard to teleosts. Here, we performed single-cell RNA sequencing of ~9500 testicular cells from the male, orange-spotted grouper. In the adult testis, we divided the cells into nine clusters and defined ten cell types, as compared with human testis data, including cell populations with characteristics of male germ cells and somatic cells, each of which expressed specific marker genes. We also identified and profiled the expression patterns of four marker genes (calr, eef1a, s100a1, vasa) in both the ovary and adult testis. Our data provide a blueprint of male germ cells and supporting somatic cells. Moreover, the cell markers are candidates that could be used for further cell identification.

Comparison of differential expression genes in ovaries and testes of Pearlscale angelfish Centropyge vrolikii based on RNA-Seq analysis

Pearlscale angelfish Centropyge vrolikii is a kind of protogynous hermaphrodite fish with a natural sexual reversion. Under appropriate social conditions, a female fish can transform into a male fish spontaneously. It is an important prerequisite for artificial breeding to understand the process of its gonadal development and sexual reversion. Gonadal development is regulated by many sex-related genes. In this study, we used unreferenced RNA-Seq technology to sequence the ovary at the perinucleolus stage (O_II), ovary at the yolk vesicle stage (O_IV),IV and testis (T), respectively; screened the gonadal differential expression genes (DEGs); and analyzed the expression of these genes in different developmental stages of ovary and different sex gonads. The results showed that a total of 142,589 all-unigene samples were assembled, and gene annotation was performed by COG, GO, KEGG, KOG, Pfam, Swissprot, eggNOG, and NR functional database. Comparative analysis revealed that there were 1919 genes that were up-regulated and 1289 genes were down-regulated in comparison to O_IV vs O_II, while there were 3653 genes that were up-regulated and 2874 genes were down-regulated in comparison of O_IV vs T, there were 3345 genes that were up-regulated and 2995 genes were down-regulated in comparison of the O_II vs the T. At the same time, the results verified by RT-qPCR were consistent with the variation trend of transcriptome data. Among the results, amh, sox9b, dmrt1, dmrt2, cyp11a, cyp17a, and cyp19a were significantly expressed in the testes, while sox3, sox4, sox11, sox17, and hsd3b7 were significantly expressed in the ovaries. And, the expression of the amh, sox9b, dmrt2, and dmrt1 were low in the O_II and O_IV, while significantly increased during the ovotestis in the hermaphroditic period (O_T), and finally reached the highest level in pure testis after sex reversal. The expression of sox3, sox4, hsd3b7, sox11, and sox17 was significantly reduced during the hermaphroditic period (O_T). These results suggested that these genes may play an important role in the process of sex reversal. This study is helpful to further understand the molecular regulation mechanism of gonadal development and sexual reversion in Pearlscale angelfish and also provide important clues for future studies.

Parental selection for growth and early-life low stocking density increase the female-to-male ratio in European sea bass

In European sea bass (Dicentrarchus labrax), as in many other fish species, temperature is known to influence the sex of individuals, with more males produced at relatively high temperatures. It is however unclear to what extent growth or stress are involved in such a process, since temperature is known to influence both growth rate and cortisol production. Here, we designed an experiment aiming at reducing stress and affecting early growth rate. We exposed larvae and juveniles originating from both captive and wild parents to three different treatments: low stocking density, food supplemented with tryptophan and a control. Low stocking density and tryptophan treatment respectively increased and decreased early growth rate. Each treatment influenced the stress response depending on the developmental stage, although no clear pattern regarding the whole-body cortisol concentration was found. During sex differentiation, fish in the low-density treatment exhibited lower expression of gr1, gr2, mr, and crf in the hypothalamus when compared to the control group. Fish fed tryptophan displayed lower crf in the hypothalamus and higher level of serotonin in the telencephalon compared to controls. Overall, fish kept at low density produced significantly more females than both control and fish fed tryptophan. Parents that have been selected for growth for three generations also produced significantly more females than parents of wild origin. Our findings did not allow to detect a clear effect of stress at the group level and rather point out a key role of early sexually dimorphic growth rate in sex determination.

RNA sequencing: new technologies and applications in cancer research

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.