Identification and expression profiles of Fox transcription factors in the Yesso scallop (Patinopecten yessoensis)

A new insight into potential roles of Spfoxl-2 in the testicular development of Scylla paramamosain by RNAi and transcriptome analysis

In our previous study, we found that the Spfoxl-2 transcript was highly expressed in gonads and explored its potential target genes in the ovary of Scylla paramamosain. In the current study, we primally analyzed its potential target genes in the testis through RNAi and RNA-Seq technology and compared with that in the ovary. The results showed that a total of 7892 unigenes were differentially expressed after Spfoxl-2 silencing in the testis, including plenty of conserved genes involved in testicular development, such as Dmrt family genes, Sox family genes, Caspase family genes, Cdk family genes, Kinesin family genes, Fox family genes and other genes. Further analysis revealed that these differentially expressed genes (DEGs) were enriched in crucial pathways involved in spermatogenesis, such as DNA replication, Cell cycle, Spliceosome, Homologous recombination, Meiosis and Apoptosis. The comparison results of potential target genes in the ovary and testis reveal 135 common potential target genes, including some genes involved in the immune response. According to our knowledge, the present work was the first to disclose the functions of foxl-2 in the testis of crustacean species using transcriptome analysis. It not only identifies key genes and pathways involved in the testicular development of S. paramamosain, but also reveals a new molecular-level understanding of the function of foxl-2 in testicular development.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.

The Fox Gene Repertoire in the Annelid Owenia fusiformis Reveals Multiple Expansions of the foxQ2 Class in Spiralia

Fox genes are a large and conserved family of transcription factors involved in many key biological processes, including embryogenesis and body patterning. Although the role of Fox genes has been studied in an array of model systems, comprehensive comparative studies in Spiralia-a large clade of invertebrate animals including molluscs and annelids-are scarce but much needed to better understand the evolutionary history of this gene family. Here, we reconstruct and functionally characterize the Fox gene complement in the annelid Owenia fusiformis, a slow evolving species and member of the sister group to all remaining annelids. The genome of O. fusiformis contains at least a single ortholog for 20 of the 22 Fox gene classes that are ancestral to Bilateria, including an ortholog of the recently discovered foxT class. Temporal and spatial expression dynamics reveal a conserved role of Fox genes in gut formation, mesoderm patterning, and apical organ and cilia formation in Annelida and Spiralia. Moreover, we uncover an ancestral expansion of foxQ2 genes in Spiralia, represented by 11 paralogs in O. fusiformis. Notably, although all foxQ2 copies have apical expression in O. fusiformis, they show variable spatial domains and staggered temporal activation, which suggest cooperation and sub-functionalization among foxQ2 genes for the development of apical fates in this annelid. Altogether, our study informs the evolution and developmental roles of Fox genes in Annelida and Spiralia generally, providing the basis to explore how regulatory changes in Fox gene expression might have contributed to developmental and morphological diversification in Spiralia.

Comparative Transcriptome Reveals the Potential Modulation Mechanisms of Spfoxl-2 Affecting Ovarian Development of Scylla paramamosain

Previously, we reported the identification, tissue distribution and confirmed the roles of Spfoxl-2 in regulating vitellogenin (vtg) expression in Scylla paramamosain. Here, we primally analyzed its potential target genes in the ovary with RNAi and RNA-Seq technology. By comparing the transcriptome data of two groups (ovaries that injected with EGFP and Foxl-2 siRNA, respectively), we found 645 DEGs (differentially expressed genes), including several conserved crucial genes involved in ovarian development, such as vtg, vitellogenin receptor (vtgR), adenylate cyclase (AC), cyclinB, and cell division cycle 2 (cdc2). In addition, these DEGs were also enriched in pathways related to ovary development, including relaxin signaling pathway, ovarian steroidogenesis, and progesterone-mediated oocyte maturation. Moreover, several genes were selected for qRT-PCR to validate the accuracy of the bioinformatic result. To the best of our knowledge, the current study was the first report about foxl-2 function through comparative transcriptome analysis in crustacean species, which identified not only relevant genes and pathways involved in ovarian development of S. paramamosain, but also provided new insights into the regulatory mechanisms of foxl-2 at the molecular level in crustacean.

Comparative Transcriptome and DNA Methylation Analysis of Phenotypic Plasticity in the Pacific Abalone (Haliotis discus hannai)

Phenotypic plasticity is an adaptive mechanism used by organisms to cope with environmental fluctuations. Pacific abalone (Haliotis discus hannai) are large-scale farmed in the temperate area of northern China and in the warmer waters of southern China. RNA-seq and comparative transcriptomic analysis here were performed to determine if the northern and southern populations have evolved divergent plasticity and if functional differences are associated with protein synthesis and growth-related biological progress. The DNA methylation (5mC) landscape of H. discus hannai from the two populations using whole genomic bisulfite sequencing (WGBS), exhibited different epigenetic patterns. The southern population had significant genomic hypo-methylation that may have resulted from long-term acclimation to heat stress. Combining 790 differentially expressed genes (DEGs) and 7635 differentially methylated genes (DMGs), we found that methylation within the gene body might be important in predicting abalone gene expression. Genes related to growth, development, transduction, and apoptosis may be regulated by methylation and could explain the phenotypic divergence of H. discus hannai. Our findings not only emphasize the significant roles of adaptive plasticity in the acclimation of H. discus hannai to high temperatures but also provide a new understanding of the epigenetic mechanism underlying the phenotypic plasticity in adaptation to climate change for marine organisms.

Genome-Wide Identification, Comparison, and Expression Analysis of Transcription Factors in Ascidian Styela clava

Tunicates include diverse species, as they are model animals for evolutionary developmental biology study. The embryonic development of tunicates is known to be extensively regulated by transcription factors (TFs). Styela clava, the globally distributed invasive tunicate, exhibits a strong capacity for environmental adaptation. However, the TFs were not systematically identified and analyzed. In this study, we reported 553 TFs categorized into 60 families from S. clava, based on the whole genome data. Comparison of TFs analysis among the tunicate species revealed that the gene number in the zinc finger superfamily displayed the most significant discrepancy, indicating this family was under the highly evolutionary selection and might be related to species differentiation and environmental adaptation. The greatest number of TFs was discovered in the Cys2His2-type zinc finger protein (zf-C2H2) family in S. clava. From the point of temporal view, more than half the TFs were expressed at the early embryonic stage. The expression correlation analysis revealed the existence of a transition for TFs expression from early embryogenesis to the later larval development in S. clava. Eight Hox genes were identified to be located on one chromosome, exhibiting different arrangement and expression patterns, compared to Ciona robusta (C. intestinalis type A). In addition, a total of 23 forkhead box (fox) genes were identified in S. clava, and their expression profiles referred to their potential roles in neurodevelopment and sensory organ development. Our data, thus, provides crucial clues to the potential functions of TFs in development and environmental adaptation in the leathery sea squirt.

Integrated mRNA and miRNA expression profile analysis of female and male gonads in Hyriopsis cumingii

Hyriopsis cumingii is an important species for freshwater pearl cultivation in China. In terms of pearl production, males have larger pearls and better glossiness than females, but there are few reports focusing on the sex of H. cumingii. In this study, six mRNA and six microRNA (miRNA) libraries were prepared from ovaries and testes. Additionally, 28,502 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Compared with testis, 14,360 mRNAs and 20 miRNAs were up-regulated in ovary, 14,142 mRNAs and 12 miRNAs were down-regulated. In DEGs, the known genes related to sex determinism and/or differentiation were also identified, such as DMRT1, SOX9, SF1 for males, FOXL2 for females, and other potentially significant candidate genes. Three sex-related pathways have also been identified, which are Wnt, Notch, and TGF-beta. In 32 DEMs, the three miRNAs (miR-9-5p, miR-92, miR-184) were paid more attention, they predicted 28 target genes, which may also be candidates for sex-related miRNAs and genes. Differential miRNAs target genes analysis reveals the pathway associated with oocyte meiosis and spermatogenesis. Overall, the findings of the study provide significant insights to enhance our understanding of sex differentiation and/or sex determination mechanisms for H. cumingii.

Circ0061052 regulation of FoxC1/Snail pathway via miR-515-5p is involved in the epithelial-mesenchymal transition of epithelial cells during cigarette smoke-induced airway remodeling

Circular RNA (circRNA) has been shown to be widely involved in a variety of lung diseases. Cigarette smoke (CS) may induce epithelial-mesenchymal transition (EMT) of airway remodeling in chronic obstructive pulmonary disease (COPD), however, in which the roles and mechanisms of circRNA have not been elucidated. In this study, we aimed to determine whether circ0061052 is involved in the EMT of human bronchial epithelial (HBE) cells and its potential mechanism for playing a biological role. Cigarette smoke extract (CSE) caused elevated EMT indicators and the increases of circ0061052 in HBE cells. Circ0061052 has a ring structure and is mainly present in the cytoplasm of HBE cells. We analyzed the regulatory relationship between circ0061052 and miR-515-5p using bioinformatics, a luciferase reporter gene, and qRT-PCR. We found that circ0061052 is mainly distributed in the cytoplasm and competitively binds to miR-515-5p, acting as a sponge for miR-515-5p. The luciferase reporter gene showed that miR-515-5p binds to the 3'UTR region of FoxC1 mRNA to inhibit its transcription. For HBE cells, overexpression of miR-515-5p antagonized the CSE-induced EMT. In addition, circ0061052 acts by binding miR-515-5p competitively to regulate the expression of FoxC1/Snail. When circ0061052 siRNA and miR-515-5p inhibitor were co-transfected into HBE cells, the inhibitor reversed the effect of circ0061052 siRNA on reducing EMT. Chronic exposure of mice to CS induced increases of circ0061052 levels, decreases of miR-515-5p levels, and the EMT in lung tissue, which caused dysfunction and airway obstruction. Overall, the results show that, by regulating miR-515-5p through a FoxC1/Snail regulatory axis, circ0061052 is involved in the CS-induced EMT and airway remodeling in COPD.