Next Article in Journal
Faecal Cortisol Metabolites as an Indicator of Adrenocortical Activity in Farmed Blue Foxes
Previous Article in Journal
Understanding the Benefits, Challenges, and the Role of Pet Ownership in the Daily Lives of Community-Dwelling Older Adults: A Case Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 9
10.3390/ani11092630
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis

by
Euna Jo
1,2,†,
Seung-Jae Lee
1,†,
Eunkyung Choi
1,
Jinmu Kim
1,
Jun-Hyuck Lee
3,4,* and
Hyun Park
1,*
1
Division of Biotechnology, College of Life Sciences and Biotechnology, Korea University, Seoul 02841, Korea
2
Division of Life Sciences, Korea Polar Research Institute (KOPRI), Incheon 21990, Korea
3
Research Unit of Cryogenic Novel Material, Korea Polar Research Institute (KOPRI), Incheon 21990, Korea
4
Department of Polar Sciences, University of Science and Technology, Incheon 21990, Korea
*
Authors to whom correspondence should be addressed.
†
Equal contribution.
Animals 2021, 11(9), 2630; https://doi.org/10.3390/ani11092630
Submission received: 9 July 2021 / Revised: 30 August 2021 / Accepted: 2 September 2021 / Published: 7 September 2021
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The brine shrimp Artemia is a promising model organism for ZW sex determination system, but the genes related to sex determination and differentiation of Artemia have not yet been examined in detail. In this study, the first isoform-level transcriptome sequencing was performed on female and male Artemia franciscana. By using PacBio Iso-Seq and Illumina RNA-Seq technologies, we found 39 candidate sex determination genes that showed sex-biased gene expression. The male-biased expressed genes included DMRT1 and Sad genes, which had three and seven isoforms, respectively. Among these, the Sad gene is an ecdysteroid biosynthetic pathway gene associated with arthropod molting and metamorphosis. We propose the importance and the necessity of further research on genes involved in ecdysteroid biosynthesis. These results will contribute to understand sex determination and differentiation of Artemia and other crustaceans having ZW systems.

Abstract

The brine shrimp Artemia has a ZW sex determination system with ZW chromosomes in females and ZZ chromosomes in males. Artemia has been considered a promising model organism for ZW sex-determining systems, but the genes involved in sex determination and differentiation of Artemia have not yet been identified. Here, we conducted transcriptome sequencing of female and male A. franciscana using PacBio Iso-Seq and Illumina RNA-Seq techniques to identify candidate sex determination genes. Among the 42,566 transcripts obtained from Iso-Seq, 23,514 were analyzed. Of these, 2065 (8.8%) were female specific, 2513 (10.7%) were male specific, and 18,936 (80.5%) were co-expressed in females and males. Based on GO enrichment analysis and expression values, we found 10 female-biased and 29 male-biased expressed genes, including DMRT1 and Sad genes showing male-biased expression. Our results showed that DMRT1 has three isoforms with five exons, while Sad has seven isoforms with 2–11 exons. The Sad gene is involved in ecdysteroid signaling related to molting and metamorphosis in arthropods. Further studies on ecdysteroid biosynthetic genes are needed to improve our understanding of Artemia sex determination. This study will provide a valuable resource for sex determination and differentiation studies on Artemia and other crustaceans with ZW systems.

1. Introduction

Most animals develop into two sexes, male or female, which are generally determined by a pair of sex chromosomes. In contrast to the male-heterogametic (XY) sex determination system found in most mammals, some organisms, including birds, snakes, fish, insects, and crustaceans, have a female-heterogametic (ZW) sex determination system [1]. Sex determination is an interesting theme, especially for researchers studying evolution and development. Accordingly, many studies have been conducted on sex-determining genes in several organisms, such as mammals [2,3,4], birds [5,6,7], insects [8], and crustaceans [9,10,11,12,13,14,15]. Primary sex determination depends on the sex chromosomes and is responsible for gonadal determination [16]. Well-known primary sex-determining genes are the sex-determining region of the Y chromosome (SRY) gene in humans [17] and the Z-linked doublesex and mab-3-related transcription factor 1 (DMRT1) gene in chickens [5]. Regarding secondary sex determination, sexual phenotypes are regulated by sex hormones or sex steroids secreted from gonads in vertebrates [16,18,19]. Ecdysteroids, the molting hormones in invertebrates, have been proposed to play roles as sex hormones in secondary sex determination [20,21]. They are suggested to be involved in sexual dimorphism and have much wider functions in insects [20,22,23] and crustaceans [24,25,26,27,28].
The brine shrimp Artemia is a member of the branchiopod crustaceans and is distributed worldwide in various saline habitats [29]. In the genus, seven species undergo bisexual reproduction, and numerous parthenogenetic populations, also known as A. parthenogenetica, can reproduce asexually [30,31,32]. Female and male Artemia are easily distinguished by their phenotypes. Females are usually 1–2 mm larger than males and have an ovisac. Males have a paired penis and markedly enlarged second antennae [33,34]. The genus Artemia has a ZW sex determination system, which means ZW chromosomes in females and ZZ chromosomes in males [35,36]. Since Bowen [37] first proposed the heterogametic sex chromosomes in Artemia females through sex-linked mutation experiments, the ZW system of sex determination has been supported and confirmed by amplified fragment length polymorphism (AFLP)-based linkage maps [36] and Z-linked gene expression analysis [35]. The sex-determining genes in Artemia have not yet been identified; however, some genes related to sexual differentiation have been identified. DMRT, or doublesex (Dsx), which is widely expressed in many vertebrates, insects, and crustaceans [12,15,38,39,40,41], has been suggested to be involved in sexual differentiation in Artemia [42]. The CCCH-type zinc finger gene Masculinizer (Masc), which was proposed to control the sex-specific expression of the Dsx gene of the silkworm Bombyx mori [43], has also been identified in Artemia [33]. Nevertheless, a recent study revealed that these genes are located in the autosomal region and are not primary sex determination genes [35]. It is implied that there may be different master sex determination genes unknown in other animals, so further studies on the sex determination genes of Artemia are required.
Recently, many studies on sex determination/differentiation genes through transcript expression profiles have made widespread use of long-read sequencing platforms. In particular, PacBio Iso-Seq can generate full-length transcripts without assembly, allowing splicing variants of complex isoforms to be revealed [44,45,46]. Individual genes generate various mRNA isoforms by alternative splicing, which can produce multiple functional proteins [47]. Therefore, it is necessary to conduct a transcriptome study at the isoform level. To date, sex-related transcriptome research of A. franciscana has been conducted using the 454 GS FLX titanium platform [48] and Illumina HiSeq system [35]. However, no transcriptome studies have been conducted on full-length transcripts at the isoform level in Artemia. In the present study, we sequenced A. franciscana, a bisexual Artemia species, using PacBio Iso-Seq and Illumina RNA-Seq technologies to identify genes with sex-biased expression, thereby attempting to identify candidate sex determination genes. The genus Artemia is considered to be a model organism for studying ZW sex-determining systems due to its short life cycle; ease of raising; and the appearance of two reproductive strategies, sexual and asexual [33,35,36]. Moreover, Artemia has utility as an organism for toxicological testing and as a food source for aquaculture industries [49,50,51,52]. This study provides a basis for the sex-biased gene expression and isoform-level transcriptome study of Artemia and helps to understand sex determination and differentiation in crustacean species using the ZW sex determination system.

2. Materials and Methods

2.1. Sample Preparation and RNA Extraction

Artemia franciscana were hatched from commercial cysts (INVE Technologies NV, Dendermonde, Belgium) and cultured in 30 g/L saline water at 25 °C with aeration. The live green alga Tetraselmis sp. was used as a food source. Adult females and males of A. franciscana were collected separately, mainly based on the size of the second antennae (Figure 1). Before RNA extraction, they were not fed for 3 days to reduce the influence of their food. After washing three times with distilled water, four pooled females and males were homogenized in TRIzol Reagent (Ambion, Austin, TX, USA) followed by total RNA extraction using the Direct-zol RNA MicroPrep Kit (Zymo Research, Irvine, CA, USA). RNA quantity and quality were measured using a Qubit fluorometer (Invitrogen, Life Technologies, Carlsbad, CA, USA) and Fragment Analyzer (Agilent Technologies, Santa Clara, CA, USA), respectively.

2.2. Iso-Seq Sequencing for Full-Length Transcripts Generation

Before constructing the PacBio Iso-Seq library, RNAs from different sexes of A. franciscana were mixed in equal amounts. Complementary DNA (cDNA) was synthesized using the SMARTer PCR cDNA Synthesis Kit (Clontech, CA, USA). The Iso-Seq SMRTbell library was prepared using a DNA template prep kit 1.0 (Pacific Biosciences, Menlo Park, CA, USA) according to the manufacturer’s protocol. The library was sequenced using two SMRT cells and Sequel Sequencing Chemistry 3.0 (Pacific Biosciences, Menlo Park, CA, USA) on the Sequel platform with 20 h movie times. Raw data were processed using SMRT Link v8.0.0 with the Iso-Seq3 application (Pacific Biosciences, Menlo Park, CA, USA). Circular consensus sequences (CCSs) were built from subreads and classified into full-length non-chimeric (FLNC) reads or non-full length (NFL) reads by identifying polyA/T tails and 5′ and 3′ primers. The FLNC reads were clustered using iterative clustering for error correction (ICE) algorithm to produce consensus isoforms and then polished with NFL reads based on an Arrow accuracy of 99%.

2.3. RNA-Seq Sequencing and Mapping to Iso-Seq Data

Illumina RNA-Seq libraries were prepared individually for female and male RNAs, which were the same as those used for the Iso-Seq library construction. Following the standard protocol, paired-end libraries were constructed using a Truseq Stranded mRNA Prep kit (Illumina, San Diego, CA, USA) and sequenced on the Illumina NovaSeq 6000 platform. The quality of the raw data was assessed using the FastQC application [53]. Using CLC Genomics Workbench v12.0.3 (QIAGEN, Aarhus, Denmark), the female and male RNA-Seq data were severally mapped to the 42,566 transcripts generated from PacBio Iso-Seq. The expression levels were calculated as transcripts per million (TPM) values for each transcript. A TPM cutoff value of 1 was used to exclude transcripts with no or extremely low expression. Transcripts with a TPM ≥ 1 exclusively in one sex were considered sex-specific transcripts.

2.4. Functional Annotation and Sex-Biased Expression Analysis

The full-length transcripts generated by Iso-Seq were annotated using BLASTX [54] by aligning to the NCBI non-redundant (nr) protein database and Trinotate v2.3.0 pipeline [55]. Gene ontology (GO) term mapping and annotation were performed in terms of the biological process (BP), molecular function (MF), and cellular component (CC) using Blast2GO v5.0 [56,57]. To avoid contaminant sequences from the culture environment, the transcripts belonging to plants, algae, bacteria, and viruses were filtered out based on the BLAST best-hit species. For the remaining female and male-specific transcripts, GO enrichment analysis with a two-tailed Fisher’s exact test (p < 0.05) was conducted using OmicsBox v1.3.11 [56,58]. The heatmaps were displayed based on TPM values for transcripts assigned to the top two BP GO terms, which were enriched specifically in females or males, using MultiExperiment Viewer (MeV) v4.9.0 [59]. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway mapping was performed using the KEGG Automatic Annotation Server (KAAS) [60] and KEGG Mapper [61]. The term ‘sex-biased genes’ refers to genes that are expressed only in one sex (sex-specific) or are more highly expressed in one sex (sex-enriched) [62].

3. Results

3.1. Transcriptome Sequencing Results of Female and Male A. franciscana

To identify transcripts differentially expressed in adult females and males of A. franciscana, Iso-Seq and RNA-Seq were performed, and the data were analyzed following the bioinformatic steps shown in Figure 1. We used four adult individuals (i.e., four females or four males) for RNA extraction to minimize sample variation. The Iso-Seq library was prepared by pooling both RNAs, whereas RNA-Seq libraries were constructed independently for each RNA. Iso-Seq sequencing generated 678,253 subreads with a total read length of 2,128,734,891 bp and a mean length of 3138 bp (Table S1). A total of 42,566 high-quality isoforms were obtained after clustering and polishing the Iso-Seq data (Table S1). RNA-Seq sequencing produced 113 million reads with a total yield of ~11.4 Gb for females and 114 million reads with a total yield of ~11.5 Gb for males (Table S1).
We decided the sex-specific transcriptome sets by choosing transcripts with TPM values of one or higher in only one sex and filtering out the putative contaminant sequences. Consequently, a total of 23,514 out of 42,566 transcripts remained, among which 2065 (8.8%) were female specific, 2513 (10.7%) were male specific, and 18,936 (80.5%) were common in females and males (Figure S1, Table S2). The number of male-specific transcripts was higher than that of female-specific transcripts in A. franciscana.

3.2. Functional Annotation and GO Enrichment Analysis

To compare functional information between sex-specific transcripts in A. franciscana, GO analysis was conducted for three categories (i.e., BP, MF, and CC). First, GO-term distributions were analyzed by level 2 for 2065 female- and 2513 male-specific transcripts, respectively. In the BP category, the majority of transcripts were related to ‘cellular process’ GO terms in both females and males, followed by ‘metabolic process’ in females and ‘biological regulation’ in males (Figure 2). Additionally, the transcripts were the most distributed in the binding (MF) and cellular anatomical entity (CC) GO terms in both sexes (Figure 2). The overall GO types and proportions showed similar patterns between females and males, although there were slight differences in the BP and MF categories.
Next, we performed GO enrichment analysis for the female- and male-specific transcripts, with p-value < 0.05, to investigate which GO terms were sex specific and over-represented. After filtering to reduce to the most specific GO terms, the numbers of significantly enriched GO terms were 84 BP, 36 MF, and 19 CC in females and 209 BP, 60 MF, and 43 CC in males (Tables S3 and S4). Among these, the top 10 BP GO terms, which contain the largest number of transcripts, are presented in Figure 3. In the female set, the most abundant GO terms were ‘response to light stimulus’ (2.04%), ‘protein stabilization’ (1.55), and ‘cellular response to heat’ (1.11%) (Figure 3A). For the male set, the most enriched GO terms included ‘axon guidance’ (3.38%), ‘negative regulation of transcription by RNA polymerase II’ (2.98%), and ‘morphogenesis of embryonic epithelium’ (2.15%) (Figure 3B).

3.3. Sex-Biased Gene Expression Profile and Candidate Sex Determination Genes

In this study, we identified genes with sex-biased expression to predict candidate sex determination genes in A. franciscana. First, sex-biased transcripts assigned to the enriched GO terms were examined, especially for the major GO terms belonging to the BP category. The numbers of transcripts selected for subsequent analysis were as follows: 42 transcripts in ‘response to light stimulus’ (GO:0009416) and 32 transcripts in ‘protein stabilization’ (GO:0050821) in females and 85 transcripts in ‘axon guidance’ (GO:0007411) and 75 transcripts in ‘negative regulation of transcription by RNA polymerase II’ (GO:0000122) in males. Heatmaps for these transcripts, generated based on the TPM values and BLAST annotation, show differential expression between females and males (Figure 4).
We used two criteria to identify the sex-differentially expressed genes among the full-length transcripts: (1) the transcripts had a TPM value of at least two, which were considered ‘highly’ expressed specifically in one sex, and (2) the corresponding gene names contained only female- or male-specific isoforms based on the BLAST description. The number of identified genes with sex-differential expression was 10 in females and 29 in males (Table 1). Of the listed genes, DMRT1 and Shadow (Sad) showing male-biased expression were focused on in further isoform analyses.

3.4. Isoforms of DMRT1 and Sad Genes in A. franciscana

The mRNA isoforms of DMRT1 and Sad genes were mapped to the genome assembly (unpublished data) using Geneious v9.1.2 (Biomatters Ltd., Auckland, New Zealand) and compared. Our results showed that three male-biased DMRT1 isoforms and one Dsx isoform are expressed in both sexes (Table S5). Isoforms of the DMRT1 gene were 3191–3960 bp in length with five exons, of which only the lengths of the fifth exon were different (Figure S2, Table S5). The isoform of the Dsx gene was 1497 bp in length with six exons (Figure S2, Table S5).
Regarding Sad genes, seven male-biased isoforms were identified, out of which five isoforms showed only male-specific expression (Figure S2, Table S5). Isoforms of the Sad gene were 1865–4319 bp in length with 2 to 11 exons and showed various alternative splicing patterns (Figure S2). Ecdysteroid signaling, including the Sad gene, is associated with sexual dimorphism in addition to molting and metamorphosis in arthropods. Thus, we focused on the presence and expression of the genes that mediate ecdysteroid biosynthesis, namely Halloween genes, in A. franciscana.

3.5. KEGG Pathway for Ecdysteroid Biosynthesis

KEGG pathway analysis for insect hormone biosynthesis showed that A. franciscana harbored the spook (Spo), spookier (Spok), disembodied (Dib), and Sad genes among typical Halloween genes, whereas the phantom (Phm) gene did not appear (Figure 5). The Dib gene expressed a female-specific isoform, and the Sad gene expressed five male-specific isoforms (Figure 5). Additionally, the aldehyde dehydrogenase (ALDH) gene in the juvenile hormone (JH) biosynthetic pathway expressed one female- and one male-specific isoform, respectively (Figure 5).

4. Discussion

In the present study, we conducted transcriptome sequencing at the isoform level and compared sex-biased gene expression of brine shrimp A. franciscana by combining PacBio IsoSeq and Illumina RNA-Seq. As a result of A. franciscana transcriptome sequencing, the number of male-specific transcripts (2513) was higher than that of female-specific transcripts (2065) (Figure S1). Sex-biased genes with male-biased expression tend to evolve faster, resulting in more functional divergence than those with female-biased expression [62,63]. This tendency, thought to contribute to sexual selection, is observed in various organisms with the XY system, such as Drosophila and Caenorhabditis elegans, as well as in organisms with the ZW system, such as birds [62,63,64]. In case of crustaceans, transcriptomic analysis of the Norway lobster Nephrops norvegicus in multiple tissues showed a higher number of male-specific transcripts compared with females [14]. However, the opposite results appeared in Pacific white shrimp Litopenaeus vannamei and water flea Daphnia species [65,66,67]. Thus, further discussion on the evolutionary rates of sex-biased genes in crustaceans is indicated.
Based on GO enrichment analysis and TPM expression values, we found 39 candidate sex-determining genes (Table 1). Of the candidate genes, 29 genes showed male-biased expression, representatively DMRT1 and Sad genes. The DMRT or Dsx gene has been reported to be involved in sex determination and differentiation in diverse species, including Artemia, and is well established as dominant in male sexual differentiation [42]. A novel DMRT gene in Chinese mitten crab Eriocheir sinensis played an essential role in testicular differentiation and was only expressed in testis [41]. Dsx homologous genes of freshwater prawn Macrobrachium rosenbergii and D. magna were suggested to play roles in male gonad development, although some of the transcripts were expressed in both testis and ovaries [12,39,40]. Farazmand et al. [42], using bisexual A. urmiana and asexual A. parthenogenetica, demonstrated that DMRT genes are involved in female ovarian development but not in male testicular differentiation. This was an unusual result, suggesting the need for research on alternative spliced forms of the DMRT family. In this study, one Dsx and three DMRT1 isoforms were identified with higher or exclusive expression in males without female-specific isoforms (Figure S2, Table S5), which seems to be related to male differentiation as in other organisms. Therefore, further studies of DMRT1 isoforms are needed to identify whether there is a difference in expression depending on the various Artemia species.
The Sad gene is a member of the Halloween gene family, which encodes the cytochrome P450 enzymes that mediate the ecdysteroidogenic pathway (steroid hormone biosynthesis) [68]. Ecdysteroid signaling, including the Sad gene, is associated with sexual dimorphism in addition to molting and metamorphosis in insects. For example, ecdysteroids contribute to extreme sexual dimorphisms in mealybugs, whose mRNA expression levels from the second nymphal instar to adult stages were higher in males than in females [69]. In beetles, ecdysteroids have been reported to control horn length dimorphisms in male heads [70], and Halloween genes, including Shade (Shd) and Phm, have been suggested as candidate target genes that are downregulated by the Dsx gene [71]. These Halloween genes are found not only in insects but also in crustaceans [72,73]. In D. pulex, which belongs to Branchiopoda along with Artemia, orthologs of typical Halloween genes (i.e., Spo/Spok, Phm, Dib, Sad, and Shd) were presented and conserved [72]. In our results, most of the Halloween genes were retained in A. franciscana transcriptome except for the Phm gene (Figure 5). Among them, Dib and Sad genes had one female- and five male-specifically expressed isoforms, respectively (Figure S2, Table S5). Taken together, there may be differences in the expression pattern of Halloween genes among arthropods. The functions and roles of Halloween genes in Artemia have not yet been studied extensively; hence, additional validation studies on ecdysteroid biosynthetic pathway genes should be conducted.
Although many studies on sex determination/differentiation on crustaceans have been conducted, compared to active studies on Decapoda (e.g., crabs, lobsters, and shrimps), research on Branchiopoda is still limited to Daphnia. Species that are easy to grow and have availability of genome sequences and embryonic developmental data are useful as new model organisms for crustaceans [74]. As mentioned in the introduction, Artemia has several advantages as a model organism for studying sex determination systems. If chromosomal-level genome sequencing of Artemia proceeds in the future, it would be possible to reveal the primary sex determination gene located on the sex chromosome. Furthermore, subsequent transcriptome research according to the embryonic and larval developmental stages would help understand the sex determination system and sexual differentiation of Artemia more clearly.

5. Conclusions

Here, we aimed to investigate genes with sex-biased expression and to suggest candidate sex determination genes in A. franciscana by transcriptomic analysis using Iso-Seq and RNA-Seq data. We identified DMRT1 and Sad as candidate genes with sex-biased expression, which were related to male sex determination/differentiation in many animals. The Sad gene regulates ecdysteroid biosynthesis and is involved in secondary sex determination and sexual dimorphism. Further research on ecdysteroidogenic genes will improve our understanding of sex determination in Artemia. To the best of our knowledge, this is the first isoform-level transcriptome study of Artemia. These results will be a valuable resource for the study of sex-biased gene expression in Artemia and other crustaceans with ZW sex determination systems.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11092630/s1, Figure S1: Venn diagram showing the number of transcripts in female and male Artemia franciscana; Figure S2: Illustration of isoforms of DMRT1, Dsx, and Sad genes in Artemia franciscana; Table S1: Statistics of Iso-Seq and RNA-Seq sequencing results for Artemia franciscana; Table S2: Summary of annotations of Artemia franciscana transcriptome; Table S3: Female-specific enriched Gene Ontology (GO) terms; Table S4: Male-specific enriched Gene Ontology (GO) terms; Table S5: Information of isoforms of DMRT1, Dsx, and Sad genes in Artemia franciscana.

Author Contributions

Conceptualization, H.P.; methodology, E.J. and H.P.; software, S.-J.L.; validation, E.J.; formal analysis, E.J., S.-J.L., E.C. and J.K.; investigation, E.J.; resources, J.-H.L. and H.P.; data curation, E.J.; writing—original draft preparation, E.J., S.-J.L., J.-H.L. and H.P.; writing—review and editing, all authors; visualization, E.J.; supervision, H.P.; project administration, H.P.; funding acquisition, J.-H.L. and H.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Development of Potential Antibiotic Compounds Using Polar Organisms Resources (grant numbers 15250103, KOPRI Grant PM21030) funded by the Ministry of Oceans and Fisheries, Korea and a Korea University Grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequencing data produced in this study have been deposited in the NCBI Sequence Read Archive (SRA) under accession numbers SRR14598203–SRR14598205 and BioProject accession number PRJNA449186.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hake, L.; O’Connor, C. Genetic mechanisms of sex determination. Nat. Educ. 2008, 1, 25. [Google Scholar]
  2. Huang, S.; Ye, L.; Chen, H. Sex determination and maintenance: The role of DMRT1 and FOXL2. Asian J. Androl. 2017, 19, 619–624. [Google Scholar] [PubMed]
  3. Greenfield, A. Understanding sex determination in the mouse: Genetics, epigenetics and the story of mutual antagonisms. J. Genet. 2015, 94, 585–590. [Google Scholar] [CrossRef]
  4. Schafer, A.J.; Goodfellow, P.N. Sex determination in humans. Bioessays 1996, 18, 955–963. [Google Scholar] [CrossRef] [PubMed]
  5. Smith, C.A.; Roeszler, K.N.; Ohnesorg, T.; Cummins, D.M.; Farlie, P.G.; Doran, T.J.; Sinclair, A.H. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature 2009, 461, 267–271. [Google Scholar] [CrossRef]
  6. Smith, C.A. Sex determination in birds: HINTs from the W sex chromosome? Sex. Dev. 2007, 1, 279–285. [Google Scholar] [CrossRef]
  7. Smith, C.A.; Sinclair, A.H. Sex determination: Insights from the chicken. Bioessays 2004, 26, 120–132. [Google Scholar] [CrossRef]
  8. Clough, E.; Oliver, B. Genomics of sex determination in Drosophila. Brief. Funct. Genom. 2012, 11, 387–394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Waiho, K.; Fazhan, H.; Zhang, Y.; Li, S.; Zhang, Y.; Zheng, H.; Ikhwanuddin, M.; Ma, H. Comparative profiling of ovarian and testicular piRNAs in the mud crab Scylla paramamosain. Genomics 2020, 112, 323–331. [Google Scholar] [CrossRef] [PubMed]
  10. Levy, T.; Sagi, A. The “IAG-switch”—A key controlling element in decapod crustacean sex differentiation. Front. Endocrinol. 2020, 11. [Google Scholar] [CrossRef] [PubMed]
  11. Song, C.; Liu, L.; Hui, M.; Liu, Y.; Liu, H.; Cui, Z. Primary molecular basis of androgenic gland endocrine sex regulation revealed by transcriptome analysis in Eriocheir sinensis. J. Oceanol. Limnol. 2019, 37, 223–234. [Google Scholar] [CrossRef]
  12. Abu Abayed, F.A.; Manor, R.; Aflalo, E.D.; Sagi, A. Screening for Dmrt genes from embryo to mature Macrobrachium rosenbergii prawns. Gen. Comp. Endocrinol. 2019, 282, 113205. [Google Scholar] [CrossRef] [PubMed]
  13. Guo, Q.; Li, S.; Lv, X.; Xiang, J.; Sagi, A.; Manor, R.; Li, F. A putative insulin-like androgenic gland hormone receptor gene specifically expressed in male Chinese shrimp. Endocrinology 2018, 159, 2173–2185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rotllant, G.; Nguyen, T.V.; Sbragaglia, V.; Rahi, L.; Dudley, K.J.; Hurwood, D.; Ventura, T.; Company, J.B.; Chand, V.; Aguzzi, J.; et al. Sex and tissue specific gene expression patterns identified following de novo transcriptomic analysis of the Norway lobster, Nephrops norvegicus. BMC Genom. 2017, 18, 622. [Google Scholar] [CrossRef] [Green Version]
  15. Chebbi, M.A.; Becking, T.; Moumen, B.; Giraud, I.; Gilbert, C.; Peccoud, J.; Cordaux, R. The genome of Armadillidium vulgare (Crustacea, Isopoda) provides insights into sex chromosome evolution in the context of cytoplasmic sex determination. Mol. Biol. Evol. 2019, 36, 727–741. [Google Scholar] [CrossRef] [Green Version]
  16. Gilbert, S.F. Chromosomal Sex Determination in Mammals. In Developmental Biology, 6th ed.; Sinauer Associates: Sunderland, MA, USA, 2000. [Google Scholar]
  17. Sinclair, A.H.; Berta, P.; Palmer, M.S.; Hawkins, J.R.; Griffiths, B.L.; Smith, M.J.; Foster, J.W.; Frischauf, A.-M.; Lovell-Badge, R.; Goodfellow, P.N. A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif. Nature 1990, 346, 240–244. [Google Scholar] [CrossRef] [Green Version]
  18. Negri, I.; Pellecchia, M. Sex steroids in insects and the role of the endosymbiont Wolbachia: A new perspective. In Sex Hormones; Dubey, R., Ed.; InTech: Rijeka, Croatia, 2012; pp. 353–374. [Google Scholar]
  19. Nakamura, M. The mechanism of sex determination in vertebrates—Are sex steroids the key-factor? J. Exp. Zool. A Ecol. Genet. Physiol. 2010, 313, 381–398. [Google Scholar] [CrossRef]
  20. De Loof, A. Ecdysteroids: The overlooked sex steroids of insects? Males: The black box. Insect Sci. 2006, 13, 325–338. [Google Scholar] [CrossRef]
  21. Chang, E.S.; Chang, S.A.; Mulder, E.P. Hormones in the lives of crustaceans: An overview. Am. Zool. 2001, 41, 1090–1097. [Google Scholar] [CrossRef] [Green Version]
  22. Brown, M.R.; Sieglaff, D.H.; Rees, H.H. Gonadal ecdysteroidogenesis in Arthropoda: Occurrence and regulation. Annu. Rev. Entomol. 2009, 54, 105–125. [Google Scholar] [CrossRef]
  23. Gilbert, L.I.; Rybczynski, R.; Warren, J.T. Control and biochemical nature of the ecdysteroidogenic pathway. Annu. Rev. Entomol. 2002, 47, 883–916. [Google Scholar] [CrossRef]
  24. Rotllant, G.; Takac, P.; Liu, L.; Scott, G.L.; Laufer, H. Role of ecdysteroids and methyl farnesoate in morphogenesis and terminal moult in polymorphic males of the spider crab Libinia emarginata. Aquaculture 2000, 190, 103–118. [Google Scholar] [CrossRef]
  25. Laufer, H.; Ahl, J.; Rotllant, G.; Baclaski, B. Evidence that ecdysteroids and methyl farnesoate control allometric growth and differentiation in a crustacean. Insect Biochem. Mol. Biol. 2002, 32, 205–210. [Google Scholar] [CrossRef]
  26. Zaleski, M.A.F.; Tamone, S.L. Relationship of molting, gonadosomatic index, and methyl farnesoate in male snow crab (Chionoecetes opilio) from the eastern Bering Sea. J. Crustac. Biol. 2014, 34, 764–772. [Google Scholar] [CrossRef] [Green Version]
  27. Ayanath, A.; Arath Raghavan, S.D. Profiling of methyl farnesoate in relation to female reproductive cycle in the freshwater crab, Travancoriana schirnerae Bott, 1969 (Crustacea: Gecarcinucidae). Invertebr. Reprod. Dev. 2020, 64, 281–290. [Google Scholar] [CrossRef]
  28. Nan, F.-H.; Wu, Y.-S.; Chang, N.-C. The effect of steroid hormone feeds on the reproductive biology of the spiny lobster, Panulirus interruptus (J.W. Randall, 1840) (Decapoda, Palinura). Crustaceana 2015, 88, 1367–1386. [Google Scholar] [CrossRef]
  29. Van Stappen, G. Zoogeography. In Artemia: Basic and Applied Biology; Abatzopoulos, T.J., Beardmore, J.A., Clegg, J.S., Sorgeloos, P., Eds.; Springer: Dordrecht, The Netherlands, 2002; pp. 171–224. [Google Scholar]
  30. Asem, A.; Eimanifar, A.; Rastegar-Pouyani, N.; Hontoria, F.; De Vos, S.; Van Stappen, G.; Sun, S.-C. An overview on the nomenclatural and phylogenetic problems of native Asian brine shrimps of the genus Artemia Leach, 1819 (Crustacea, Anostraca). ZooKeys 2020, 902, 1–15. [Google Scholar] [CrossRef] [PubMed]
  31. Asem, A.; Eimanifar, A.; Sun, S.C. Genetic variation and evolutionary origins of parthenogenetic Artemia (Crustacea: Anostraca) with different ploidies. Zool. Scr. 2016, 45, 421–436. [Google Scholar] [CrossRef]
  32. Maniatsi, S.; Baxevanis, A.D.; Kappas, I.; Deligiannidis, P.; Triantafyllidis, A.; Papakostas, S.; Bougiouklis, D.; Abatzopoulos, T.J. Is polyploidy a persevering accident or an adaptive evolutionary pattern? The case of the brine shrimp Artemia. Mol. Phylogenet. Evol. 2011, 58, 353–364. [Google Scholar] [CrossRef]
  33. Li, D.-R.; Ye, H.-L.; Yang, J.-S.; Yang, F.; Wang, M.-R.; De Vos, S.; Vuylsteke, M.; Sorgeloos, P.; Van Stappen, G.; Bossier, P.; et al. Identification and characterization of a Masculinizer (Masc) gene involved in sex differentiation in Artemia. Gene 2017, 614, 56–64. [Google Scholar] [CrossRef]
  34. Van Stappen, G. Introduction, biology and ecology of Artemia. In Manual on the Production and Use of LIVE food for Aquaculture; Lavens, P., Sorgeloos, P., Eds.; FAO: Rome, Italy, 1996; pp. 79–106. [Google Scholar]
  35. Huylmans, A.K.; Toups, M.A.; Macon, A.; Gammerdinger, W.J.; Vicoso, B. Sex-biased gene expression and dosage compensation on the Artemia franciscana Z-chromosome. Genome Biol. Evol. 2019, 11, 1033–1044. [Google Scholar] [CrossRef] [Green Version]
  36. De Vos, S.; Bossier, P.; Van Stappen, G.; Vercauteren, I.; Sorgeloos, P.; Vuylsteke, M. A first AFLP-based genetic linkage map for brine shrimp Artemia franciscana and its application in mapping the sex locus. PLoS ONE 2013, 8, e57585. [Google Scholar]
  37. Bowen, S.T. The genetics of Artemia salina. II. White eye, a sex-linked mutation. Biol. Bull. 1963, 124, 17–23. [Google Scholar] [CrossRef]
  38. Mawaribuchi, S.; Ito, Y.; Ito, M. Independent evolution for sex determination and differentiation in the DMRT family in animals. Biol. Open 2019, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Kato, Y.; Kobayashi, K.; Watanabe, H.; Iguchi, T. Environmental sex determination in the branchiopod crustacean Daphnia magna: Deep conservation of a Doublesex gene in the sex-determining pathway. PLoS Genet. 2011, 7, e1001345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Zhong, P.; Zhou, T.; Zhang, Y.; Chen, Y.; Yi, J.; Lin, W.; Guo, Z.; Xu, A.; Yang, S.; Chan, S. Potential involvement of a DMRT family member (Mr—Dsx) in the regulation of sexual differentiation and moulting in the giant river prawn Macrobrachium rosenbergii. Aquac. Res. 2019, 50, 3037–3049. [Google Scholar] [CrossRef]
  41. Zhang, E.-F.; Qiu, G.-F. A novel Dmrt gene is specifically expressed in the testis of Chinese mitten crab, Eriocheir sinensis. Dev. Genes Evol. 2010, 220, 151–159. [Google Scholar] [CrossRef] [PubMed]
  42. Farazmand, A.; Inanloo, K.; Agh, N. Expression of DMRT family genes during gonadal differentiation in two species of Artemia (Branchiopoda, Anostraca) from Urmia Lake (Iran). Crustaceana 2010, 83, 1153–1165. [Google Scholar] [CrossRef]
  43. Sakai, H.; Sakaguchi, H.; Aoki, F.; Suzuki, M.G. Functional analysis of sex-determination genes by gene silencing with LNA–DNA gapmers in the silkworm, Bombyx mori. Mech. Dev. 2015, 137, 45–52. [Google Scholar] [CrossRef] [PubMed]
  44. Minio, A.; Massonnet, M.; Figueroa-Balderas, R.; Vondras, A.M.; Blanco-Ulate, B.; Cantu, D. Iso-Seq allows genome-independent transcriptome profiling of grape berry development. G3 Genes Genomes Genet. 2019, 9, 755–767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Conesa, A.; Madrigal, P.; Tarazona, S.; Gomez-Cabrero, D.; Cervera, A.; McPherson, A.; Szcześniak, M.W.; Gaffney, D.J.; Elo, L.L.; Zhang, X.; et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016, 17, 1–19. [Google Scholar] [CrossRef] [Green Version]
  46. Rhoads, A.; Au, K.F. PacBio sequencing and its applications. Genom. Proteom. Bioinform. 2015, 13, 278–289. [Google Scholar] [CrossRef] [Green Version]
  47. Matlin, A.J.; Clark, F.; Smith, C.W.J. Understanding alternative splicing: Towards a cellular code. Nat. Rev. Mol. Cell Biol. 2005, 6, 386–398. [Google Scholar] [CrossRef]
  48. Valenzuela-Miranda, D.; Gallardo-Escárate, C.; Valenzuela-Muñoz, V.; Farlora, R.; Gajardo, G. Sex-dependent transcriptome analysis and single nucleotide polymorphism (SNP) discovery in the brine shrimp Artemia franciscana. Mar. Genomics 2014, 18, 151–154. [Google Scholar] [CrossRef]
  49. Jo, E.; Lee, S.J.; Choi, E.; Kim, J.; Lee, S.G.; Lee, J.H.; Kim, J.-H.; Park, H. Whole genome survey and microsatellite motif identification of Artemia franciscana. Biosci. Rep. 2021, 41, BSR20203868. [Google Scholar] [CrossRef] [PubMed]
  50. Bengtson, D.A.; Léger, P.; Sorgeloos, P. Use of Artemia as a food source for aquaculture. In Artemia Biology; Browne, R.A., Sorgeloos, P., Trotman, C.N.A., Eds.; CRC Press: Boca Raton, FL, USA, 2018; pp. 255–286. [Google Scholar]
  51. Nunes, B.S.; Carvalho, F.D.; Guilhermino, L.M.; Van Stappen, G. Use of the genus Artemia in ecotoxicity testing. Environ. Pollut. 2006, 144, 453–462. [Google Scholar] [CrossRef]
  52. Yu, J.; Lu, Y. Artemia spp. model-a well-established method for rapidly assessing the toxicity on an environmental perspective. Med. Res. Arch. 2018, 6, 1–15. [Google Scholar]
  53. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data, Version 0.11.9; Babraham Bioinformatics: Cambridge, UK, 2019. [Google Scholar]
  54. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007, 35, D61–D65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Bryant, D.M.; Johnson, K.; DiTommaso, T.; Tickle, T.; Couger, M.B.; Payzin-Dogru, D.; Lee, T.J.; Leigh, N.D.; Kuo, T.-H.; Davis, F.G. A tissue-mapped axolotl de novo transcriptome enables identification of limb regeneration factors. Cell Rep. 2017, 18, 762–776. [Google Scholar] [CrossRef] [Green Version]
  56. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef]
  57. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [Green Version]
  58. OmicsBox-Bioinformatics made easy. Version 1.3.11. BioBam Bioinformatics: Valencia, Spain, 2019.
  59. Saeed, A.; Sharov, V.; White, J.; Li, J.; Liang, W.; Bhagabati, N.; Braisted, J.; Klapa, M.; Currier, T.; Thiagarajan, M. TM4: A free, open-source system for microarray data management and analysis. Biotechniques 2003, 34, 374–378. [Google Scholar] [CrossRef] [Green Version]
  60. Moriya, Y.; Itoh, M.; Okuda, S.; Yoshizawa, A.C.; Kanehisa, M. KAAS: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35, W182–W185. [Google Scholar] [CrossRef] [Green Version]
  61. Kanehisa, M.; Sato, Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 2020, 29, 28–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Ellegren, H.; Parsch, J. The evolution of sex-biased genes and sex-biased gene expression. Nat. Rev. Genet. 2007, 8, 689–698. [Google Scholar] [CrossRef] [PubMed]
  63. Parsch, J.; Ellegren, H. The evolutionary causes and consequences of sex-biased gene expression. Nat. Rev. Genet. 2013, 14, 83–87. [Google Scholar] [CrossRef] [PubMed]
  64. Naurin, S.; Hansson, B.; Bensch, S.; Hasselquist, D. Why does dosage compensation differ between XY and ZW taxa? Trends Genet. 2010, 26, 15–20. [Google Scholar] [CrossRef] [PubMed]
  65. Molinier, C.; Reisser, C.M.; Fields, P.; Ségard, A.; Galimov, Y.; Haag, C.R. Identification of general patterns of sex-biased expression in Daphnia, a genus with environmental sex determination. G3 Genes Genomes Genet. 2018, 8, 1523–1533. [Google Scholar]
  66. Wang, Y.; Yu, Y.; Li, S.; Zhang, X.; Xiang, J.; Li, F. Sex-Specific transcriptome sequencing of zoea I larvae and identification of sex-linked genes using bulked segregant analysis in Pacific white shrimp Litopenaeus vannamei. Mar. Biotechnol. 2020, 22, 423–432. [Google Scholar] [CrossRef]
  67. Huylmans, A.K.; López Ezquerra, A.; Parsch, J.; Cordellier, M. De novo transcriptome assembly and sex-biased gene expression in the cyclical parthenogenetic Daphnia galeata. Genome Biol. Evol. 2016, 8, 3120–3139. [Google Scholar] [CrossRef] [Green Version]
  68. Rewitz, K.; Rybczynski, R.; Warren, J.T.; Gilbert, L.I. The Halloween genes code for cytochrome P450 enzymes mediating synthesis of the insect moulting hormone. Biochem. Soc. Trans. 2006, 34, 1256–1260. [Google Scholar] [CrossRef]
  69. Muramatsu, M.; Tsuji, T.; Tanaka, S.; Shiotsuki, T.; Jouraku, A.; Miura, K.; Vea, I.M.; Minakuchi, C. Sex-specific expression profiles of ecdysteroid biosynthesis and ecdysone response genes in extreme sexual dimorphism of the mealybug Planococcus kraunhiae (Kuwana). PLoS ONE 2020, 15, e0231451. [Google Scholar]
  70. Emlen, D.; Nijhout, H. Hormonal control of male horn length dimorphism in Onthophagus taurus (Coleoptera: Scarabaeidae): A second critical period of sensitivity to juvenile hormone. J. Insect Physiol. 2001, 47, 1045–1054. [Google Scholar] [CrossRef]
  71. Ledón-Rettig, C.C.; Zattara, E.; Moczek, A.P. Asymmetric interactions between doublesex and tissue-and sex-specific target genes mediate sexual dimorphism in beetles. Nat. Commun. 2017, 8, 1–10. [Google Scholar] [CrossRef] [PubMed]
  72. Rewitz, K.F.; Gilbert, L.I. Daphnia Halloween genes that encode cytochrome P450s mediating the synthesis of the arthropod molting hormone: Evolutionary implications. BMC Evol. Biol. 2008, 8, 60. [Google Scholar] [CrossRef] [Green Version]
  73. Xie, X.; Liu, Z.; Liu, M.; Tao, T.; Shen, X.; Zhu, D. Role of Halloween genes in ecdysteroids biosynthesis of the swimming crab (Portunus trituberculatus): Implications from RNA interference and eyestalk ablation. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2016, 199, 105–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Toyota, K.; Miyakawa, H.; Hiruta, C.; Sato, T.; Katayama, H.; Ohira, T.; Iguchi, T. Sex Determination and Differentiation in Decapod and Cladoceran Crustaceans: An Overview of Endocrine Regulation. Genes 2021, 12, 305. [Google Scholar] [CrossRef] [PubMed]
Animals 11 02630 g001 550
Figure 1. Workflow of the bioinformatic steps. The photographs of female and male Artemia franciscana (scale bar = 2 mm), and the flowchart showing analysis tools used in this study.
Figure 1. Workflow of the bioinformatic steps. The photographs of female and male Artemia franciscana (scale bar = 2 mm), and the flowchart showing analysis tools used in this study.
Animals 11 02630 g001
Animals 11 02630 g002 550
Figure 2. Gene Ontology (GO) distribution of sex-specific transcripts in Artemia franciscana transcriptome. Functional annotation was conducted for female-specific transcripts (A) and male-specific transcripts (B). GO terms were divided into three major categories: biological process (BP, green), molecular function (MF, blue), and cellular component (CC, yellow). The x-axis indicates the number of transcripts, and the y-axis indicates level 2 GO terms.
Figure 2. Gene Ontology (GO) distribution of sex-specific transcripts in Artemia franciscana transcriptome. Functional annotation was conducted for female-specific transcripts (A) and male-specific transcripts (B). GO terms were divided into three major categories: biological process (BP, green), molecular function (MF, blue), and cellular component (CC, yellow). The x-axis indicates the number of transcripts, and the y-axis indicates level 2 GO terms.
Animals 11 02630 g002
Animals 11 02630 g003 550
Figure 3. Top 10 enriched Gene Ontology (GO) terms of sex-specific transcripts in Artemia franciscana transcriptome. Enrichment analysis was performed for female-specific transcripts (A) and male-specific transcripts (B). Biological process GO terms that exhibit statistically significant differences are shown in the figure (Fisher’s exact test, p < 0.05). The X-axis represents the relative frequencies, and the Y-axis represents significant enriched GO terms.
Figure 3. Top 10 enriched Gene Ontology (GO) terms of sex-specific transcripts in Artemia franciscana transcriptome. Enrichment analysis was performed for female-specific transcripts (A) and male-specific transcripts (B). Biological process GO terms that exhibit statistically significant differences are shown in the figure (Fisher’s exact test, p < 0.05). The X-axis represents the relative frequencies, and the Y-axis represents significant enriched GO terms.
Animals 11 02630 g003
Animals 11 02630 g004 550
Figure 4. Heatmap of sex-biased expressed genes involved in top two enriched Gene Ontology (GO) terms. Heatmaps were generated for female-specific transcripts (A) and male-specific transcripts (B). Each line represents a transcript annotated by BLAST, and the color scale indicates relative expression level by transcripts per million (TPM) values.
Figure 4. Heatmap of sex-biased expressed genes involved in top two enriched Gene Ontology (GO) terms. Heatmaps were generated for female-specific transcripts (A) and male-specific transcripts (B). Each line represents a transcript annotated by BLAST, and the color scale indicates relative expression level by transcripts per million (TPM) values.
Animals 11 02630 g004
Animals 11 02630 g005 550
Figure 5. Expression differences in female and male Artemia franciscana transcriptome for the insect hormone biosynthesis KEGG pathway. Colored boxes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway indicate the presence of female-specific isoforms (yellow), male-specific isoforms (purple), and isoforms expressed in both sexes (green). Heatmaps represent the expression patterns of ALDH (1.2.1.3, green), Dib (yellow), and Sad (purple) genes. The color scales indicate relative expression levels by transcripts per million (TPM) values. ALDH, aldehyde dehydrogenase; Dib, disembodied; Sad, shadow.
Figure 5. Expression differences in female and male Artemia franciscana transcriptome for the insect hormone biosynthesis KEGG pathway. Colored boxes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway indicate the presence of female-specific isoforms (yellow), male-specific isoforms (purple), and isoforms expressed in both sexes (green). Heatmaps represent the expression patterns of ALDH (1.2.1.3, green), Dib (yellow), and Sad (purple) genes. The color scales indicate relative expression levels by transcripts per million (TPM) values. ALDH, aldehyde dehydrogenase; Dib, disembodied; Sad, shadow.
Animals 11 02630 g005
Table
Table 1. Selected candidate genes with sex-biased expression in the Artemia franciscana transcriptome.
Table 1. Selected candidate genes with sex-biased expression in the Artemia franciscana transcriptome.
F/MDescriptionSex-Biased IsoformTotal Isoform
FemaleCalnexin37
DNA damage-binding protein 115
Ferrochelatase, mitochondrial25
Heat shock protein 901525
HLA-B associated transcript 313
Innexin inx1-like14
Protein ROP-like isoform 117
Ras-related protein Rab-6A isoform X127
Synapsin-like12
TFIIH basal transcription factor complex helicase XPD Subunit-like24
Male5’-AMP-activated protein kinase catalytic subunit alpha-225
7SK snRNA methylphosphate capping enzyme48
Armadillo segment polarity protein, putative15
cGMP-dependent protein kinase110
Cullin 324
Cytoplasmic FMR1-interacting protein-like22
Disheveled associated activator of morphogenesis14
Doublesex- and mab-3-related transcription factor 123
Exportin 114
Glutamine–tRNA ligase-like13
Homeobox protein not211
Immunodeficiency virus type I enhancer binding protein, putative25
Integrin alpha-PS1-like12
Integrin alpha-PS2-like111
Integrin beta-PS215
Kinesin-like protein KIF3B13
Kruppel-like factor 5-like13
Leucine-rich repeat-containing protein 4B13
Microtubule-actin cross-linking factor 1-like629
Optomotor-blind protein-like12
Ras-related C3 botulinum toxin substrate 1 (rho family, Small GTP binding protein Rac1)13
RNA polymerase-associated protein CTR9 homolog13
Serine/threonine protein phosphatase 2A catalytic subunit alpha isoform28
Shadow57
Transcription factor GATA-314
Triple functional domain protein-like712
Tyrosine protein phosphatase Lar-like13
Zinc finger protein1330
Zinc finger transcription factor15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jo, E.; Lee, S.-J.; Choi, E.; Kim, J.; Lee, J.-H.; Park, H. Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis. Animals 2021, 11, 2630. https://doi.org/10.3390/ani11092630

AMA Style

Jo E, Lee S-J, Choi E, Kim J, Lee J-H, Park H. Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis. Animals. 2021; 11(9):2630. https://doi.org/10.3390/ani11092630

Chicago/Turabian Style

Jo, Euna, Seung-Jae Lee, Eunkyung Choi, Jinmu Kim, Jun-Hyuck Lee, and Hyun Park. 2021. "Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis" Animals 11, no. 9: 2630. https://doi.org/10.3390/ani11092630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop