De novo transcriptome assembly from the gonads of a scleractinian coral, Euphyllia ancora: molecular mechanisms underlying scleractinian gametogenesis

Gene expression profiles of Japanese precious coral Corallium japonicum during gametogenesis

Background

Corallium japonicum, a prized resource in Japan, plays a vital role in traditional arts and fishing industries. Because of diminished stock due to overexploitation, ongoing efforts are focused on restoration through transplantation. This study aimed to enhance our understanding of the reproductive biology of these valuable corals and find more efficient methods for sex determination, which may significantly contribute to conservation initiatives.

Methods

We used 12 three-month aquarium reared C. japonicum colony fragments, conducted histological analysis for maturity and sex verification, and performed transcriptome analysis via de novo assembly and mapping using the C. rubrum transcriptome to explore gene expression differences between female and male C. japonicum.

Results

Our histological observations enabled sex identification in 33% of incompletely mature samples. However, the sex of the remaining 67% of samples, classified as immature, could not be identified. RNA-seq yielded approximately 21-31 million short reads from 12 samples. De novo assembly yielded 404,439 highly expressed transcripts. Among them, 855 showed significant differential expression, with 786 differentially expressed transcripts between females and males. Heatmap analysis highlighted 283 female-specific and 525 male-specific upregulated transcripts. Transcriptome assembly mapped to C. rubrum yielded 28,092 contigs, leading to the identification of 190 highly differentially expressed genes, with 113 upregulated exclusively in females and 70 upregulated exclusively in males. Blastp analysis provided putative protein annotations for 83 female and 72 male transcripts. Annotation analysis revealed that female biological processes were related to oocyte proliferation and reproduction, whereas those in males were associated with cell adhesion.

Discussion

Transcriptome analysis revealed sex-specific gene upregulation in incompletely mature C. japonicum and shared transcripts with C. rubrum, providing insight into its gene expression patterns. This study highlights the importance of using both de novo and reference-based assembly methods. Functional enrichment analysis showed that females exhibited enrichment in cell proliferation and reproduction pathways, while males exhibited enrichment in cell adhesion pathways. To the best of our knowledge, this is the first report on the gene expressions of each sex during the spawning season. Our findings offer valuable insights into the physiological ecology of incompletely mature red Japanese precious corals and suggest a method for identifying sex using various genes expressed in female and male individuals. In the future, techniques such as transplantation, artificial fertilization, and larval rearing may involve sex determination methods based on differences in gene expression to help conserve precious coral resources and ecosystems.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

An External Coincidence Model for the Lunar Cycle Reveals Circadian Phase-Dependent Moonlight Effects on Coral Spawning

Many marine organisms synchronously spawn at specific times to ensure the success of external fertilization in the ocean. Corals are famous examples of synchronized spawning at specific lunar phases, and two distinct spawning patterns have been observed in two dominant taxa: merulinid corals spawn at regular lunar phases, several days after the full moon, whereas Acropora corals spawn at more irregular lunar phases around the full moon. Although it has been suggested that the two coral taxa have different responses to moonlight and seawater temperature, their spawning times have never been analyzed by integrating the two environmental factors, resulting in an incomplete understanding of the regulatory mechanisms of spawning. In this study, we developed a new predictive model of coral spawning days by integrating moonlight and temperature effects based on the external coincidence model for the lunar cycle. We performed model fitting using a 10-year monitoring record of coral spawning time in Taiwan. Our model successfully demonstrated the synergistic effects of moonlight and temperature on coral spawning time (days) and provided two testable hypotheses to explain the different spawning patterns regarding the preparation (maturation) process for spawning and the sensitivity to moonlight at different circadian phases: (1) Acropora corals may have an earlier onset and longer period of preparation for spawning than merulinid corals; and (2) merulinid corals may use moonlight signals near sunset, while Acropora corals may have a similar onset at approximately midnight. This is the first study to indicate the difference in circadian phase-dependent moonlight sensitivities between coral taxa, providing a basis for underlying coral spawning mechanisms for rhythmic studies.

Administration of olaquindox impairs spermatogenesis and sperm quality by increasing oxidative stress and early apoptosis in mice

Olaquindox (OLA), a potent antibacterial agent, has been widely used as a feed additive and growth promoter in animal husbandry. Our previous study has shown that OLA administration in female mice could markedly cause sub-fertility. Here we established the model in male mice to investigate the toxic effects of OLA on mammalian spermatozoa quality and fetal development. After continuous 45 days of OLA gavage, the dosage of 60 mg/kg/day (high dose) significantly affected body weight, organ weights and coefficients, and the morphology of the testis seminiferous tubule in male mice. Dosage of 60 mg/kg/day also reduced sperm count, motility, and viability. OLA at both low-dose (5 mg/kg/day) and high-dose induced peroxidation, early apoptosis, and abnormal mitochondrial membrane potential in sperm. Significantly, high-dose OLA impaired in vitro fertilized embryo development, indicated by the decreased percentages of 2-cell and blastocyst formation. Surprisingly, the natural fertility of males was unaffected after OLA gavage, which was indicated by the comparable litter size after mating. However, paternal gavage of OLA significantly decreased the survival rate of the offspring from the age of 4 weeks. In sum, our study showed that OLA gavage in male mice damages sperm quality and offspring survival, illustrating the use of OLA as a feed additive should be strictly restricted.

Studying of Molecular Regulation of Developmental Processes of Lower Metazoans Exemplified by Cnidaria Using High-Throughput Sequencing

A unique set of features and characteristics of species of the Cnidaria phylum is the one reason that makes them a model for a various studies. The plasticity of a life cycle and the processes of cell differentiation and development of an integral multicellular organism associated with it are of a specific scientific interest. A new stage of development of molecular genetic methods, including methods for high-throughput genome, transcriptome, and epigenome sequencing, both at the level of the whole organism and at the level of individual cells, makes it possible to obtain a detailed picture of the development of these animals. This review examines some modern approaches and advances in the reconstruction of the processes of ontogenesis of cnidarians by studying the regulatory signal transduction pathways and their interactions.

Larval transcriptomic responses of a stony coral, Acropora tenuis, during initial contact with the native symbiont, Symbiodinium microadriaticum

Although numerous dinoflagellate species (Family Symbiodiniaceae) are present in coral reef environments, Acropora corals tend to select a single species, Symbiodinium microadriaticum, in early life stages, even though this species is rarely found in mature colonies. In order to identify molecular mechanisms involved in initial contact with native symbionts, we analyzed transcriptomic responses of Acropora tenuis larvae at 1, 3, 6, 12, and 24 h after their first contact with S. microadriaticum, as well as with non-native symbionts, including the non-symbiotic S. natans and the occasional symbiont, S. tridacnidorum. Some gene expression changes were detected in larvae inoculated with non-native symbionts at 1 h post-inoculation, but those returned to baseline levels afterward. In contrast, when larvae were exposed to native symbionts, we found that the number of differentially expressed genes gradually increased in relation to inoculation time. As a specific response to native symbionts, upregulation of pattern recognition receptor-like and transporter genes, and suppression of cellular function genes related to immunity and apoptosis, were exclusively observed. These findings indicate that coral larvae recognize differences between symbionts, and when the appropriate symbionts infect, they coordinate gene expression to establish stable mutualism.

Development of an in vitro tissue culture system for hammer coral (Fimbriaphyllia ancora) ovaries

In vitro gonad culture systems have proven useful to investigate intrinsic mechanisms of sexual reproduction in animals. Here we describe development of an in vitro culture method for coral ovaries. Mesenterial tissues containing both ovaries and mesenterial filaments were microscopically isolated from the scleractinian coral, Fimbriaphyllia ancora, and culture conditions were optimized. M199 diluted 10× (10% M199, pH 8.1) and supplemented with 25 mM HEPES and the antibiotics, ampicillin, penicillin and streptomycin, supported oocyte survival and maintained the structural integrity of ovaries during short-term culture (~ 6 days). Addition of a commercial antibiotic-antimycotic solution (Anti-Anti) and fetal bovine serum adversely affected ovary maintenance and caused tissue disintegration. Characterization of cultured ovaries showed that there is no difference in cell proliferation of ovarian somatic cells between culture Days 1 and 6. Moreover, the presence of oogonia and expression of a major yolk protein, vitellogenin, were confirmed in ovaries cultured for 6 days. This system will be useful for studying effects of a wide range of substances on coral oogenesis.

Chromosome anchoring in Senegalese sole (Solea senegalensis) reveals sex-associated markers and genome rearrangements in flatfish

The integration of physical and high-density genetic maps is a very useful approach to achieve chromosome-level genome assemblies. Here, the genome of a male Senegalese sole (Solea senegalensis) was de novo assembled and the contigs were anchored to a high-quality genetic map for chromosome-level scaffolding. Hybrid assembled genome was 609.3 Mb long and contained 3403 contigs with a N50 of 513 kb. The linkage map was constructed using 16,287 informative SNPs derived from ddRAD sequencing in 327 sole individuals from five families. Markers were assigned to 21 linkage groups with an average number of 21.9 markers per megabase. The anchoring of the physical to the genetic map positioned 1563 contigs into 21 pseudo-chromosomes covering 548.6 Mb. Comparison of genetic and physical distances indicated that the average genome-wide recombination rate was 0.23 cM/Mb and the female-to-male ratio 1.49 (female map length: 2,698.4 cM, male: 2,036.6 cM). Genomic recombination landscapes were different between sexes with crossovers mainly concentrated toward the telomeres in males while they were more uniformly distributed in females. A GWAS analysis using seven families identified 30 significant sex-associated SNP markers located in linkage group 18. The follicle-stimulating hormone receptor appeared as the most promising locus associated with sex within a region with very low recombination rates. An incomplete penetrance of sex markers with males as the heterogametic sex was determined. An interspecific comparison with other Pleuronectiformes genomes identified a high sequence similarity between homologous chromosomes, and several chromosomal rearrangements including a lineage-specific Robertsonian fusion in S. senegalensis.

Transcriptomic analysis of gonadal development in parasitic and non-parasitic lampreys (Ichthyomyzon spp.), with a comparison of genomic resources in these non-model species

Lampreys are jawless fishes that diverged ∼550 million years ago from other vertebrates. Sequencing of the somatic and the germline genomes of the sea lamprey (Petromyzon marinus) in 2013 and 2018, respectively, has helped to improve our understanding of the genes and gene networks that control many aspects of lamprey development. However, little is known about the genetic basis of gonadal differentiation in lampreys, partly due to the prolonged period during which their gonads remain sexually indeterminate. We performed RNA-sequencing on gonadal samples from four chestnut lamprey (Ichthyomyzon castaneus) and six northern brook lamprey (I. fossor) to identify differentially expressed genes (DEG’s) and pathways associated with transcriptomic differences in: (1) larvae during early gonadal differentiation versus definitive females (i.e., with oocytes in the slow cytoplasmic growth phase); and (2) females versus definitive males undergoing spermatogonial proliferation. We compared the mapping percentages of these transcriptomes to the two available sea lamprey reference genomes and three annotation files (Ensembl and UCSC for the somatic genome and SIMRbase for the germline genome). We found that mapping the RNA-seq reads to the germline genome gave superior results and, using Trinotate, we provided new putative annotations for 8161 genes in the somatic assembly and 880 genes for the germline assembly. We identified >2000 DEG’s between stages and sexes, as well as biological pathways associated with each. Interestingly, some of the upregulated genes (e.g., DEG’s associated with spermiation) suggest that changes in gene expression can precede morphological changes by several months. In contrast, only 81 DEG’s were evident between the chestnut lamprey (that remains sexually immature during an extended post-metamorphic parasitic feeding phase) and the nonparasitic northern brook lamprey (that undergoes sexual maturation near the end of metamorphosis), but few replicates were available for comparable stages and sexes. This work lays the foundation for identifying and confirming the orthology and the function of genes involved in gonadal development in these and other lamprey species across more developmental stages.