Time-course transcriptome analysis reveals regulation of Arabidopsis seed dormancy by the transcription factors WOX11/12.

The induction of seed dormancy and its release involve a finely regulated genetic program controlled by various environmental and developmental cues that are critical for plant survival and population expansion. Light plays a key role in seed dormancy and germination, but the molecular mechanisms underlying the control of dormancy are unclear. In the present study, high-resolution temporal RNA-seq in Arabidopsis identified WOX11 as encoding a hub transcription factor during the seed dormancy induction and release stages. This gene might have evolved from gymnosperms and expanded in angiosperms with highly conserved expression patterns in seeds. WOX11 and its homolog WOX12 were highly expressed from 2 d after pollination, and mRNA abundance was greatly increased during the seed dormancy induction and release stages. Further, we found that WOX11 plays a role in the regulation of seed dormancy downstream of phytochrome B (PHYB)-mediated red-light signaling during the induction stage, indicating that WOX11/12 are newly identified components of red-light signaling transduction. Taken together, our results suggest that WOX11/12-mediated PHYB signaling regulates seed dormancy in Arabidopsis, and provide insights into the developmental regulation and evolutionary adaptation of plants to changes in the light environment.

Comparative transcriptome analysis of eyes reveals the adaptive mechanism of mantis shrimp (oratosquilla oratoria) induced by a dark environment.

The light-dark cycle significantly impacts the growth and development of animals. Mantis shrimps (Oratosquilla oratoria) receive light through their complex photoreceptors. To reveal the adaptive expression mechanism of the mantis shrimp induced in a dark environment, we performed comparative transcriptome analysis with O. oratoria cultured in a light environment (Oo-L) as the control group and O. oratoria cultured in a dark environment (Oo-D) as the experimental group. In the screening of differentially expressed genes (DEGs) between the Oo-L and Oo-D groups, a total of 88 DEGs with |log2FC| > 1 and FDR < 0.05 were identified, of which 78 were upregulated and 10 were downregulated. Then, FBP1 and Pepck were downregulated in the gluconeogenesis pathway, and MKNK2 was upregulated in the MAPK classical pathway, which promoted cell proliferation and differentiation, indicating that the activity of mantis shrimp was slowed and the metabolic rate decreases in the dark environment. As a result, the energy was saved for its growth and development. At the same time, we performed gene set enrichment analysis (GSEA) on all DEGs. In the KEGG pathway analysis, each metabolic pathway in the dark environment showed a slowing trend. GO was enriched in biological processes such as eye development, sensory perception and sensory organ development. The study showed that mantis shrimp slowed down metabolism in the dark, while the role of sensory organs prominent. It provides important information for further understanding the energy metabolism and has great significance to study the physiology of mantis shrimp in dark environment.

Whole genome evaluation analysis and preliminary Assembly of Oratosquilla oratoria (Stomatopoda: Squillidae).

BACKGROUND: As the dominant species of Stomatopoda, Oratosquilla oratoria has not been fully cultivated artificially, and the fishery production mainly depends on marine fishing. Due to the lack of stomatopod genome, the development of molecular breeding of mantis shrimps still lags behind. METHODS AND RESULTS: A survey analysis was performed to obtain the genome size, GC content and heterozygosity ratio in order to provide a fundation for subsequent whole-genome sequencing. The results showed that the estimated genome size of the O. oratoria was about 2.56 G, and the heterozygosity ratio was 1.81%, indicating that it is a complex genome. Then the sequencing data was preliminarily assembled with k-mer = 51 by SOAPdenovo software to obtain a genome size of 3.01G and GC content of 40.37%. According to ReapeatMasker and RepeatModerler analysis, the percentage of repeats in O. oratoria was 45.23% in the total genome, similar to 44% in Survey analysis. The MISA tool was used to analyze the simple sequence repeat (SSR) characteristics of genome sequences including Oratosquilla oratoria, Macrobrachium nipponense, Fenneropenaeus chinensis, Eriocheir japonica sinensis, Scylla paramamosain and Paralithodes platypus. All crustacean genomes showed similar SSRs characteristics, with the highest proportion of di-nucleotide repeat sequences. And AC/GT and AGG/CCT repeats were the main types of di-nucleotide and tri-nucleotide repeats in O. oratoria. CONCLUSION: This study provided a reference for the genome assembly and annotation of the O. oratoria, and also provided a theoretical basis for the development of molecular markers of O. oratoria.

The analyses of the complete mitochondrial genomes of three crabs revealed novel gene rearrangements and phylogenetic relationships of Brachyura.

BACKGROUND: Brachyura crab is the largest branch of Decapoda crustacean. Phylogenetic relationships within Brachyura remain controversial to be investigated. The mitochondrial genome (mitogenome) is an important molecular marker for studying the phylogenetic relationships of Brachyura. METHODS AND RESULTS: To understand the phylogeny of Brachyura, the three complete mitogenomes from Charybdis annulata, Leptodius exaratus, and Spider crab were sequenced and annotated. Their full length was 15,747, 15,716, and 16,608 bp long, respectively. The first two crabs both contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. However, Spider crab contained 13 PCGs, two rRNA genes, 25 tRNA genes and a control region. The mitogenomes of each of the three crabs exhibited high AT content (67.8%, 69.1%, and 70.8%), with negative AT skews (-0.014, - 0.028, and - 0.017) and GC skews (-0.269, - 0.286, and - 0.341). The gene order of C. annulata was identical to the ancestor of Brachyura. Compared with the ancestor of Brachyura, L. exaratus exhibited the gene rearrangements of Val (V)-rrnS-control region, and Spider crab had the four copies of Lys (K). Phylogenetic analyses indicated that C. annulata belonged to Portunidae family, Portunoidea superfamilies, L. exaratus belonged to Xanthidae family, Xanthoidea superfamilies, and Spider crab belonged to Mithracidae family, Majoidea superfamilies. Phylogenetic analyses showed that the two species (Somanniathelphusa boyangensis and Huananpotamon lichuanense) belonging to the Potamoidea were sister groups to the Thoracotremata, thus supporting the conclusion that Heterotremata is polyphyletic. CONCLUSION: The results of this study enriched the crab mitogenome database and enabled us to better understand the phylogenetic relationships of Brachyura.

Transcriptome analysis of immune-related genes in Sesarmops sinensis hepatopancreas in reaction to peptidoglycan challenge.

Sesarmops sinensis is a dominant omnivorous crab species, which plays an important ecological function in salt marsh ecosystems. To better understand its immune system and immune related genes under pathogen infection, the transcriptome was analyzed by comparing the data of S. sinensis hepatopancreas stimulated by PBS and PGN. A set of assembly and annotation identified 39,039 unigenes with an average length of 1105 bp, obtaining 1300 differentially expressed genes (DEGs) in all, which included 466 remarkably up-regulated unigenes and 834 remarkably down-regulated unigenes. In addition, based on mensurable real time-polymerase chain reaction and high-throughput sequencing, several immune responsive genes were found to be markedly up-regulated under PGN stimulation. In conclusion, in addition to enriching the existing transcriptome data of S. sinensis, this study also clarified the immune response of S. sinensis to PGN stimulation, which will help us to further understand the crustacean's immune system.

A complement factor I (CFI) gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco.

This study isolated CFI gene from Pelteobagrus fulvidraco and named it PfCFI. The cDNA of PfCFI is 2374 bp long, including a 52 bp 5' untranslated sequence, a 222 bp 3' untranslated sequence, and an open reading frame (ORF) of 2100 bp encoding polypeptide consisting of 699 amino acids. Phylogenetic analysis revealed that the PfCFI was closely related to CFI of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis indicate that there is the PfCFI gene which expressed in all the rest of tested tissues in varied levels, and mainly distributed in liver and least in heart. The reseachers induce the expressions level of PfCFI gene in liver, spleen, head kidney and blood at different points in time after challenged with lipopolysaccharide (LPS), and polyriboinosinic polyribocytidylic acid (poly I:C), respectively. Together these results suggested that CFI gene plays an important role in resistance to pathogens in yellow catfish immunity.

Transcriptome analysis of differentially expressed genes in the red swamp crayfish Procambarus clarkii challenged with Aeromonas hydrophila.

As an important economic species in China, aquaculture of the crayfish Procambarus clarkii has suffered huge losses due to infection by pathogenic bacteria, mainly by Aeromonas hydrophila, which leads to high mortality and huge economic loss. To better understand the immune response of crayfish against bacterial infection, we compared and analyzed transcriptome data of hepatopancreatic tissue from P. clarkii that were either challenged with A. hydrophila or treated with PBS. After assembly and annotation of the data, 32,041 unigenes with an average length of 1512 base pairs were identified. Compared to control group, Differential gene expression (DEG) analysis revealed 608 DEGs were obtained, of which 274 unigenes were upregulated and 334 were downregulated in the A. hydrophila group. Furthermore, the expression levels of eight selected immune-related DEGs were validated by qRT-PCR, substantiating the reliability of RNA-seq results. This study not only provides effective data support for immune defense strategies of P. clarkii in response to bacterial infections, but also provides new information about the P. clarkii immune system and defense mechanisms, and a valuable basis for further studies to elucidate the molecular immune mechanisms of this species.

Transcriptome analysis reveals antioxidant defense mechanisms in the red swamp crayfish Procambarus clarkia after exposure to chromium.

Chromium (Cr) as a chromate anion has a strong redox capacity that seriously threatens the ecological environment and human health. Cr can contaminate water and impart toxicity to aquatic species. Procambarus clarkii is an important food source that once represented a large proportion of the aquaculture industry due to its rapid reproduction and high economic value. However, there have been reports on the death of P. clarkii due to heavy metal pollution. The underlying mechanism regarding heavy metal toxicity was studied in this paper. The transcriptome data of hemocytes extracted from P. clarkii injected with Cr were analyzed by high-throughput sequencing and compared to the control group. In total, 48,128,748 clean reads were obtained in the treatment group and 56,480,556 clean reads were obtained in the control group. The reads were assembled using Trinity and the identified unigenes were then annotated. Then, 421 differentially-expressed genes (DEGs) were found, 170 of which were upregulated and 251 downregulated. Many of these genes were found to be related to glutathione metabolism and transportation. The glutathione metabolic pathway of P. clarkii was thus activated by Cr exposure to detoxify and maintain body function. Validation of DEGs with quantitative real-time PCR confirms the changes in gene expression. Thus, this study provides data supporting a glutathione-focused response of P. clarkii to exposure to heavy metals.

Characterization of the complete mitochondrial genome of Helice latimera and its phylogenetic implications in Brachyura.

Mitochondrial genomes (mitogenomes) help advance our learning of molecular evolution and phylogenetic relationships. The mitogenome of H. latimera is 16,246 bp in length, which typically contains 37 animal mitogenome genes consisting of 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes, as well as a control region. The AT content of H. latimera is 69.1%. The A + T skew of the mitogenome of H. latimera was slightly negative (-0.017). The size of Thirteen PCGs is from 162 bp to 1731 bp. Twenty-two tRNA genes ranged from 62 to 73 bp and were highly A + T biased. All tRNA genes owed a typical cloverleaf structure, not including the trnS1 gene lacking a dihydroxyuridine arm. One PCG, two rRNAs, and 12 of the tRNAs were rearranged compared to the pancrustacean gene order. Phylogenetic analysis revealed the locationt of H. latimera among the Varunidae family.

Mitochondrial genome of the yellow catfish Pelteobagrus fulvidraco and insights into Bagridae phylogenetics.

The mitochondrial genome (mitogenome) can provide important information for understanding phylogenetic analysis and molecular evolution. Herein, we amplified the complete mitogenome sequence of Pelteobagrus fulvidraco. The mitogenome was 16,526 bp in length and included 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and a non-coding control region (D-loop). Both the organization and location of genes in the mitogenome were consistent with those from Siluriformes fishes previously published in GenBank. The phylogenetic relationships based on Bayesian inference (BI) and Maximum likelihood (ML) methods showed that P. fulvidraco has close relationships with Pelteobagrus eupogon and Tachysurus intermedius, suggesting that P. fulvidraco belongs to Tachysurus. This study provides evidence that Tachysurus, Pseudobagrus and Leiocassis do not form monophyly, but that these three genera form a monophyletic group. Our results provide reference for further phylogenetic research of the Bagridae species.

A non-mammalian Toll-like receptor 26 (TLR26)gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco.

In this study, we identified a fish-specific Toll-like receptor (TLR) in Pelteobagrus fulvidraco, an economically important freshwater fish in China. This TLR, PfTLR26, was shown to be encoded by a 3084 bp open reading frame (ORF), producing a polypeptide 1027 amino acids in length. The PfTLR26 protein contains a signal peptide, eight leucine-rich repeat (LRR) domains, two LRR_TYP domains in the extracellular region, and a Toll/interleukin (IL)-1 receptor (TIR) domain in the cytoplasmic region, consistent with the characteristic TLR domain architecture. This predicted 117.1 kDa protein was highly homologous to those of other fish, with phylogenetic analysis revealing the closest relation to TLR26 of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis showed that the PfTLR26 gene was expressed in all tissues tested, with the highest expression levels seen in the head kidney and blood, and the lowest seen in muscle. PfTLR26 exhibited significant upregulation in liver, spleen, head kidney, and blood at different time points following challenge with the common TLR agonists lipopolysaccharide (LPS) and polyriboinosinic polyribocytidylic acid (Poly I:C). Taken together, these results suggest that PfTLR26 may be an important component of the P. fulvidraco innate immune system, participating in the transduction of TLR signaling under pathogen stimulation.

Characterization and expression analysis of immune-related genes in the red swamp crayfish, Procambarus clarkii in response to lipopolysaccharide challenge.

To learn more about red swamp crayfish related genes in response to bacterial infections, we investigated immune-related genes induced by lipopolysaccharide (LPS) in the hepatopancreas using high-throughput sequencing method. In present the study, a total of 55,107 unigenes were identified, with an average length of 678 bp. A total of 2215 differentially expressed genes (DEGs) were found, including 669 up-regulated genes and 1546 down-regulated genes. The result of Gene ontology (GO) analysis revealed that 3017 DEGs were enriched in 19 biological process subcategories, 17 cellular component subcategories and 15 molecular function subcategories. The top 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that "ribosome" was the most abundant group, which had 34 DEGs. KEGG enrichment analysis identified several immune response pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results exhibited that several immune responsive genes were greatly up-regulated following LPS stimulation as observed in the results of high-throughput sequencing. Overall, this study provides new insight into the immune defense mechanisms of P. clarkii against LPS infection.

Transcriptome analysis of differential expressed genes in hepatopancreas of Procambarus clarkii challenged with peptidoglycan.

Procambarus clarkii is one of the most economically important species in Chinese aquaculture, and is widely cultured. Infection of P. clarkii populations with bacterial pathogens causes high mortality and great economic loss, therefore disease control is of significant economic importance. P. clarkii is a model system for studying immune responses in invertebrates, and its immune system consists solely of the innate response. In the present study, we examined gene expression related to immune function in P. clarkii in response to pathogen challenge. The transcriptome of hepatopancreas tissue from P. clarkii challenged with peptidoclycan (PGN) was analyzed and compared to control specimens. After assembly and annotation, 48,661 unigenes were identified with an average length of 671.54 bp. A total of 2533 differentially expressed genes (DEGs) were obtained, including 765 significantly up-regulated unigenes and 1757 significantly down-regulated unigenes. Gene ontology (GO) analysis demonstrated 19 biological process subcategories, 16 cellular component subcategories, and 17 molecular function subcategories that were enriched among these DEGs. Enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database revealed enrichment among immune responses pathways. Taken together, this study not only enriches the existing P. clarkii transcriptome database, but also elucidates immune responses of crayfish that are activated in response to PGN challenge.

Transcriptome-wide identification of differentially expressed genes in Procambarus clarkii in response to chromium challenge.

Because of the high protein content and rich meat quality of crayfish Procambarus clarkii, it has become widely popular in China in recent years and has a high economic value. When P. clarkii is stimulated by heavy metals, it reacts to oxidation. P. clarkii has evolved antioxidant defense systems, including antioxidant enzymes such as catalase (CAT). The hexavalent form of Cr (VI) is a pathogenic factor that is of particular concern in aqueous systems because of its great toxicity to living organisms. In this study, we characterized the transcriptome of P. clarkii using a RNA sequencing method and performed a comparison between K2Cr2O7-treated samples and controls. In total, 34,237 unigenes were annotated. We identified 5098 significantly differentially expressed genes (DEGs), including 2536 and 2562 were significantly up-regulated and down-regulated, respectively. In addition, quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. Our results contribute to a more comprehensive understanding of the antioxidant defense system used by P. clarkii in response to heavy metal stress.

New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis.

Iron in excess can have toxic effects on living organisms. In China, the freshwater crayfish Procambarus clarkii is a source of aquatic food with high-quality protein and has significant commercial value. P. clarkii shows oxidative stress on exposure to heavy metals, and antioxidant enzymes, such as ubiquitination enzymes and proteasomes, play important roles in oxidative stress. To understand the antioxidant defense system of P. clarkii, we analyzed the hepatopancreas transcriptomes of P. clarkii after stimulation with FeCl3. In total, 5199 differentially expressed genes (DEGs) were identified (2747 upregulated and 2452 downregulated). GO analysis revealed that these DEGs belonged to 16 cellular component, 16 molecular function, and 19 biological process subcategories. A total of 1069 DEGs were classified into 25 categories by using COG. Some antioxidant defense pathways, such as "Ubiquitin mediated proteolysis" and "Glutathione metabolism," were identified using KEGG. In addition, quantitative real time-PCR (qRT-PCR) substantiated the up-regulation of a random selection of DEGs including antioxidant and immune defense genes. We obtained information for P. clarkii transcriptome databases and new insights into the responses of P. clarkii hepatopancreas to heavy metals.

Adaptive evolution of osmoregulatory-related genes provides insight into salinity adaptation in Chinese mitten crab, Eriocheir sinensis.

Osmoregulation is an important mechanism by which euryhaline crustaceans regulate osmotic and ionic concentrations. The Chinese mitten crab (Eriocheir sinensis) is a strong osmoregulating animal model among crustacean species, as it can maintain its hemolymph composition and survives well in either seawater or freshwater. Osmoregulation by E. sinensis during physiological adaptation has been studied extensively. However, the genetic basis of osmoregulation in E. sinensis for acclimating to changing salinities remains unclear. The current study investigated five genes involved in E. sinensis osmoregulation and compared them with a representative marine crab Portunus trituberculatus to test whether adaptive evolution has occurred changing salinity conditions. The results showed that carbonic anhydrase (CA), cytochrome P450 4C (CYP4C), glutamate dehydrogenase (GDH), and the Na+/H+ exchanger (NHE) have undergone positive selection (i.e., directional selection) in E. sinensis. Thus, the positive selection in CA and NHE suggests that E. sinensis has enhanced capacity for maintaining systemic acid-base balance and ion regulation. GDH and CYP4C also demonstrated positive selection in E. sinensis, suggesting that E. sinensis might have acquired an enhanced capacity to metabolize glutamate and synthesize ecdysteroids in response to a change in osmotic concentration. The present study provides new insight into the molecular genetic basis of salinity adaption in E. sinensis.

Transcriptomic analysis of immune-related genes in the lipopolysaccharide-stimulated hepatopancreas of the mudflat crab Helice tientsinensis.

The mudflat crab Helice tientsinensis is one of the most economically important aquaculture species in China. Nevertheless, it is susceptible to various diseases caused by viruses, bacteria and rickettsia-like organisms. A better understanding of the immune system and genes related to the responses to bacterial and viral infection is required. Herein, the hepatopancreas transcriptome of H. tientsinensis was analyzed by comparing control and lipopolysaccharide (LPS)-stimulated RNA-Seq data, yielding 91,885,038 bp and 13.78 Gb of clean reads. Following assembly and annotation, 93,207 unigenes with an average length of 883 bp were identified, of which 31,674 and 13,700 were annotated in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Following LPS, 4845 differentially expressed genes (DEGs) were identified, of which 2491 and 2354 were up- and down-regulated, respectively. To further investigate immune-related DEGs, KEGG enrichment analysis identified immune response pathways, most notably the peroxisome and Toll-like receptor signaling pathways. Quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. This systematic transcriptomic analysis of the innate immune pathway in H. tientsinensis expands our understanding of the immune system in crabs.

De novo transcriptome assembly and analysis of differential gene expression following lipopolysaccharide challenge in Pelteobagrus fulvidraco.

The yellow catfish, Pelteobagrus fulvidraco, has been recognized as an important freshwater aquaculture species in Eastern and Southeast Asia. To gain a better understanding of the immune response in P. fulvidraco, we analyzed its transcriptome following stimulation with lipopolysaccharide (LPS). Phosphate buffer saline (PBS) was used as control. Following assembly and annotation, 72,152 unigenes with an average length of 1090 bp were identified. A total of 370 differentially expressed genes (DEGs) in the P. fulvidraco were observed at 12 h post LPS treatment, including 197 up-regulated genes and 173 down-regulated genes. Clusters of Orthologous Groups of proteins (KOG/COG) annotation demonstrated that a total of 18,819 unigenes classified into 26 categories. Gene ontology (GO) analysis revealed 20 biological process subcategories, 7 cellular component subcategories and 20 molecular function subcategories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified immune responses pathways. Quantitative reverse transcription polymerase chain reaction measured the expression of 18 genes involved in the immune response. CXCL2-like chemokine (CXCL2), goose-type lysozyme (LYZ G), and cathepsin K (CTSK) were significantly up-regulated. This study enriches the P. fulvidraco transcriptome database and provides insight into the immune response of P. fulvidraco against infection.

The complete mitochondrial genome of Clostera anachoreta (Lepidoptera: Notodontidae) and phylogenetic implications for Noctuoidea species.

In this study, the complete mitochondrial genome (mitogenome) of Clostera anachoreta (Lepidoptera: Notodontidae) was sequenced. It comprises 15,456 base pairs (bp), including 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and one non-coding control region (CR), as found in other lepidopterans. Gene order is identical to that of typical lepidopterans. There are 15 intergenic spacers ranging from 2 to 49bp, and 9 overlapping regions ranging from 1 to 8bp, occurring throughout the genome. The CR is 347bp long. All PCGs are initiated by ATN codons. We found a typical gene rearrangement in C. anachoreta (tRNAMet-tRNAIle-tRNAGln), which is different from ancestral insects (tRNAIle-tRNAGln-tRNAMet). The gene rearrangement can be explained by a duplication/random loss model. Phylogenetic analyses indicate that C. anachoreta belongs to Notodontidae, and that the monophyly of Lepidopteran families is well supported.

Characterisation of the complete mitochondrial genome of Helice wuana (Grapsoidea: Varunidae) and comparison with other Brachyuran crabs.

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Herein, we sequenced, annotated, and characterised the mitogenome of the crab Helice wuana to better understand its molecular evolution and phylogeny. The 16,359bp mitogenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region. The genome composition is highly A+T biased 68.42%, and exhibits a negative AT-skew (-0.036) and GC-skew (-0.269) among Brachyura crabs. Gene rearrangements were detected, as was tandem duplication followed by random loss, which explains the translocation of mitochondrial genes. Phylogenetic analysis showed that H. wuana and H. tientsinensis clustered on one branch with high nodal support values. These results confirm that the placement of H. wuana within the Varunidae family of Thoracotrematan crabs. This study will provided a better understanding for gene rearrangements and crab evolution in the further.