Exploring developmental gene toolkit and associated pathways in a potential new model crustacean using transcriptomic analysis

Modulation of mitochondrial dynamics genes and mtDNA during embryonic development and under UVB exposure

Studies using the embryos of the freshwater prawn Macrobrachium olfersii have reported changes in embryonic cells after exposure to ultraviolet B (UVB) radiation, such as DNA damage and apoptosis activation. Considering the importance of mitochondria in embryonic cells, this study aimed to characterize the aspects of mitochondrial morphofunctionality in M. olfersii embryos and mitochondrial responses to UVB radiation exposure. The coding sequences of genes Tfam, Nrf1, Mfn1, and Drp1 were identified from the transcriptome of M. olfersii embryos. The phylogenetic relationship showed strong amino acid identity and a highly conserved nature of the sequences. Additionally, the number of mitochondrial DNA (mtDNA) copies were higher in the early embryonic days. The results showed that the expression of the analyzed genes was highly regulated during embryonic development, increasing their levels near hatching. Furthermore, when embryos were exposed to UVB radiation, mitochondrial biogenesis was activated, recognized by higher levels of transcripts of genes Tfam and Nrf1, accompanied by mitochondrial fission. Additionally, these mitochondrial events were supported by an increase of mtDNA copies. Our results showed that UVB radiation was able to change the mitochondrial morphofunctionality, and under the current knowledge, certainly compromise embryonic cellular integrity. Additionally, mitochondria is an important cellular target of this radiation and its responses can be used to assess environmental stress caused by UVB radiation in embryos of aquatic species.

UVB radiation exposure modulates mitophagy in embryonic cells of freshwater prawn Macrobrachium olfersii: Exploring a protective organelle quality control mechanism

Aquatic environments are subject to ultraviolet B (UVB) radiation incidence, and its effects on organisms are dose-dependent. Besides DNA, mitochondria are an important target of this radiation that causes structural damage and impairs its functional dynamics. Here, we hypothesize that mitophagy acts as an organelle quality control mechanism to mitigate UVB impacts in embryonic cells. Then, freshwater prawn Macrobrachium olfersii embryos was used as a model to investigate the effects of UVB on genes (Tomm20, Opa1, Pink, Prkn, Sqstm1, and Map1lc3) and proteins (TOM20, PINK1, p62 and LC3B) involved in mitophagy modulation. The choice of genes and proteins was based on the identification of mitochondrial membrane (Tomm20, Opa1 and TOM20), mediation of mitophagy (Pink1, Prkn and PINK1), and recognition of mitochondria by the autophagosome membrane (Sqstm1, Map1lc3, p62 and LC3B). First, the phylogeny of all genes presented bootstrap values >80 and conserved domains among crustacean species. Gene expression was inherently modulated during development, with transcripts (Tomm20, Opa1, Pink, Prkn, Sqstm1, and Map1lc3) overexpressed in the initial and final stages of development. Moreover, UVB radiation induced upregulation of Tomm20, Opa1, Pink, Prkn, Sqstm1, and Map1lc3 genes at 6 h after exposure. Interestingly, after 12 h, the protein content of PINK1, p62, and LC3B increased, while TOM20 was not responsive. Despite UVB radiation's harmful effects on embryonic cells, the chronology of gene expression and protein content indicates rapid activation of mitophagy, serving as an organelle quality control mechanism, given the analyzed cells' integrity.

Phylogenetic and transcriptomic characterization of insulin and growth factor receptor tyrosine kinases in crustaceans

Receptor tyrosine kinases (RTKs) mediate the actions of growth factors in metazoans. In decapod crustaceans, RTKs are implicated in various physiological processes, such molting and growth, limb regeneration, reproduction and sexual differentiation, and innate immunity. RTKs are organized into two main types: insulin receptors (InsRs) and growth factor receptors, which include epidermal growth factor receptor (EGFR), fibroblast growth factor receptor (FGFR), vascular endothelial growth factor receptor (VEGFR), and platelet-derived growth factor receptor (PDGFR). The identities of crustacean RTK genes are incomplete. A phylogenetic analysis of the CrusTome transcriptome database, which included all major crustacean taxa, showed that RTK sequences segregated into receptor clades representing InsR (72 sequences), EGFR (228 sequences), FGFR (129 sequences), and PDGFR/VEGFR (PVR; 235 sequences). These four receptor families were distinguished by the domain organization of the extracellular N-terminal region and motif sequences in the protein kinase catalytic domain in the C-terminus or the ligand-binding domain in the N-terminus. EGFR1 formed a single monophyletic group, while the other RTK sequences were divided into subclades, designated InsR1-3, FGFR1-3, and PVR1-2. In decapods, isoforms within the RTK subclades were common. InsRs were characterized by leucine-rich repeat, furin-like cysteine-rich, and fibronectin type 3 domains in the N-terminus. EGFRs had leucine-rich repeat, furin-like cysteine-rich, and growth factor IV domains. N-terminal regions of FGFR1 had one to three immunoglobulin-like domains, whereas FGFR2 had a cadherin tandem repeat domain. PVRs had between two and five immunoglobulin-like domains. A classification nomenclature of the four RTK classes, based on phylogenetic analysis and multiple sequence alignments, is proposed.

The organizing role of Wnt signaling pathway during arthropod posterior growth

Wnt signaling pathways are recognized for having major roles in tissue patterning and cell proliferation. In the last years, remarkable progress has been made in elucidating the molecular and cellular mechanisms that underlie sequential segmentation and axial elongation in various arthropods, and the canonical Wnt pathway has emerged as an essential factor in these processes. Here we review, with a comparative perspective, the current evidence concerning the participation of this pathway during posterior growth, its degree of conservation among the different subphyla within Arthropoda and its relationship with the rest of the gene regulatory network involved. Furthermore, we discuss how this signaling pathway could regulate segmentation to establish this repetitive pattern and, at the same time, probably modulate different cellular processes precisely coupled to axial elongation. Based on the information collected, we suggest that this pathway plays an organizing role in the formation of the body segments through the regulation of the dynamic expression of segmentation genes, via controlling the caudal gene, at the posterior region of the embryo/larva, that is necessary for the correct sequential formation of body segments in most arthropods and possibly in their common segmented ancestor. On the other hand, there is insufficient evidence to link this pathway to axial elongation by controlling its main cellular processes, such as convergent extension and cell proliferation. However, conclusions are premature until more studies incorporating diverse arthropods are carried out.

Identification of Hox genes and their expression profiles during embryonic development of the emerging model organism, Macrobrachium olfersii

Hox genes encode transcription factors that specify the body segment identity during development, including crustaceans, such as amphipods and decapods, that possess a remarkable diversity of segments and specialized appendages. In amphipods, alterations of specialized appendages have been obtained using knockout experiment of Hox genes, which suggests that these genes are involved in the evolution of morphology within crustaceans. However, studies of Hox genes in crustaceans have been limited to a few species. Here, we identified the homeodomain of nine Hox genes: labial (lab), proboscipedia (pb), Deformed (Dfd), Sex combs reduced (Scr), fushi tarazu (ftz), Antennapedia (Antp), Ultrabithorax (Ubx), abdominal-A (abdA), and Abdominal-B (AbdB), and evaluated their expression by RT-qPCR and RT-PCR in the ovary, during embryonic development, and at the first larval stage (Zoea I) of the decapod Macrobrachium olfersii. The transcript levels of lab, Dfd, and ftz decreased and transcripts of pb, Scr, Antp, Ubx, abdA, and AbdB increased during embryonic development. Hox genes were expressed in mature ovaries and Zoea I larval stages, except Scr and ftz, respectively. In addition, isoforms of Dfd, Scr, Ubx, and abdA, which have been scarcely reported in crustaceans, were described. New partial sequences of 87 Hox genes from other crustaceans were identified from the GenBank database. Our results are interesting for future studies to determine the specific function of Hox genes and their isoforms in the freshwater prawn M. olfersii and to contribute to the understanding of the diversity and evolution of body plans and appendages in Crustaceans.

Regulation of gene expression during ontogeny of physiological function in the brackishwater amphipod Gammarus chevreuxi

Embryonic development is a complex process involving the co-ordinated onset and integration of multiple morphological features and physiological functions. While the molecular basis of morphological development in embryos is relatively well known for traditional model species, the molecular underpinning of the development of physiological functions is not. Here, we used global gene expression profiling to investigate the transcriptional changes associated with the development of morphological and physiological function in the amphipod crustacean Gammarus chevreuxi. We compared the transcriptomes at three timepoints during the latter half of development, characterised by different stages of the development of heart form and function: 10 days post fertilisation (dpf, Early: no heart structure visible), 15 dpf (Middle: heart present but not fully functional), and 18 dpf (Late: regular heartbeat). Gene expression profiles differed markedly between developmental stages, likely representing a change in the activity of different processes throughout the latter period of G. chevreuxi embryonic development. Differentially expressed genes belonged to one of three distinct clusters based on their expression patterns across development. One of these clusters, which included key genes relating to cardiac contractile machinery and calcium handling, displayed a pattern of sequential up-regulation throughout the developmental period studied. Further analyses of these transcripts could reveal genes that may influence the onset of a regular heartbeat. We also identified morphological and physiological processes that may occur alongside heart development, such as development of digestive caeca and the cuticle. Elucidating the mechanisms underpinning morphological and physiological development of non-model organisms will support improved understanding of conserved mechanisms, addressing the current phylogenetic gap between relatively well known model species.

The Hedgehog pathway in penaeid shrimp: developmental expression and evolution of splice junctions in Pancrustacea

Penaeid shrimp embryos undergo holoblastic division, gastrulation by invagination, and hatching as a nauplius larva. Posterior segments form and differentiate during larval development. Hedgehog (Hh) pathway genes from penaeid shrimp and other pancrustaceans were identified by in silico analysis of genomes and transcriptomes, and mapped onto a recent pancrustacean phylogeny to determine patterns of intron gains and losses. Penaeus vannamei, P. japonicus, and P. monodon Hh proteins were encoded by four exons. Amphipod, isopod, and ostracod hh were also encoded by four exons, but hh from other arthropod groups contained three conserved exons. The novel hh intron is hypothesized to have arisen independently in the malacostracan ancestor and Ostracoda by a transposon insertion. Shared patterns of ptc, smo, and ci exon structure were found for Malacostraca, Branchiopoda + Hexapoda, Hexanauplia (Thecostraca + Copepoda), Multicrustacea (Thecostraca + Copepoda + Malacostraca), and Pancrustacea minus Oligostraca. mRNA expression of P. vannamei of hh, ptc, and ci from developmental transcriptomes of zygotes through postlarvae showed low expression from zygote to gastrula, which increased at limb bud, peaked at unhatched nauplius, and declined in nauplius and later larval stages. smo expression was found in zygotes, peaked in gastrula, and declined in limb bud and later stages. These results are consistent with a role for Hh signaling during segmentation in penaeid shrimp.

Signaling Pathways That Regulate the Crustacean Molting Gland

A pair of Y-organs (YOs) are the molting glands of decapod crustaceans. They synthesize and secrete steroid molting hormones (ecdysteroids) and their activity is controlled by external and internal signals. The YO transitions through four physiological states over the molt cycle, which are mediated by molt-inhibiting hormone (MIH; basal state), mechanistic Target of Rapamycin Complex 1 (mTORC1; activated state), Transforming Growth Factor-β (TGFβ)/Activin (committed state), and ecdysteroid (repressed state) signaling pathways. MIH, produced in the eyestalk X-organ/sinus gland complex, inhibits the synthesis of ecdysteroids. A model for MIH signaling is organized into a cAMP/Ca2+-dependent triggering phase and a nitric oxide/cGMP-dependent summation phase, which maintains the YO in the basal state during intermolt. A reduction in MIH release triggers YO activation, which requires mTORC1-dependent protein synthesis, followed by mTORC1-dependent gene expression. TGFβ/Activin signaling is required for YO commitment in mid-premolt. The YO transcriptome has 878 unique contigs assigned to 23 KEGG signaling pathways, 478 of which are differentially expressed over the molt cycle. Ninety-nine contigs encode G protein-coupled receptors (GPCRs), 65 of which bind a variety of neuropeptides and biogenic amines. Among these are putative receptors for MIH/crustacean hyperglycemic hormone neuropeptides, corazonin, relaxin, serotonin, octopamine, dopamine, allatostatins, Bursicon, ecdysis-triggering hormone (ETH), CCHamide, FMRFamide, and proctolin. Contigs encoding receptor tyrosine kinase insulin-like receptor, epidermal growth factor (EGF) receptor, and fibroblast growth factor (FGF) receptor and ligands EGF and FGF suggest that the YO is positively regulated by insulin-like peptides and growth factors. Future research should focus on the interactions of signaling pathways that integrate physiological status with environmental cues for molt control.

Ultraviolet-B radiation induces transcriptional modulation of components associated with the extracellular matrix in embryos of decapod Macrobrachium olfersii

The extracellular matrix (ECM) is a non-cellular and three-dimensional structure, constituted by a macromolecular dynamic network that involves the cells in all animal tissues, including embryonic ones. Several studies with vertebrates and cell cultures have reported deleterious effects of ultraviolet-B (UVB) radiation on the components associated with the ECM. However, studies focusing on the UVB radiation effects on ECM components of crustaceans during embryonic development are very scarce. Thus, the aim of this study was to identify the coding sequences of components associated with the ECM and to evaluate the effect of UVB radiation on embryos of the ecologically-important decapod Macrobrachium olfersii. To evaluate the modulation of these ECM components during embryonic development, the transcript levels of Col4α1, Itgβ, Lamα, Mmp1 and Timp in M. olfersii embryos were analyzed at early developmental stages (E1, E3 and E4), intermediate developmental stage (E7) and late developmental stages (E10 and E14). In addition, embryos at E7, which correspond to a landmark of crustacean development, were analyzed after 12 h of UVB exposure to verify UVB effects on the ECM components. The ECM component sequences were similar to other decapods, suggesting conservation of these genes among crustaceans. The results showed modulations of the ECM components of M. olfersii embryos that reflect the need for each component in the cellular mechanisms, necessary for normal embryonic development. After UVB exposure, embryos showed opacity of embryonic tissues and it was found the overexpression of Col4α1, Itgβ, Mmp1 and Timp transcript levels (1.82-, 1.52-, 2.34- and 6.27-fold, respectively). These impairments can compromise important events for normal embryonic development, such as growth of optic lobes, caudal papilla, ramification of appendages and differentiation of organic systems. The results presented here, together with the effects on morphology, cell proliferation, differentiation, and apoptosis demonstrated previously, strengthen the knowledge of the complex impacts of UVB radiation on freshwater embryos. Nevertheless, our results encourage further investigations focusing on the assessment of UVB effects on different organisms in order to better understand the myriad of UVB effects on ECM components.

Mutually exclusive locales for N-linked glycans and disorder in human glycoproteins

Several post-translational protein modifications lie predominantly within regions of disorder: the biased localization has been proposed to expand the binding versatility of disordered regions. However, investigating a representative dataset of 500 human N-glycoproteins, we observed the sites of N-linked glycosylations or N-glycosites, to be predominantly present in the regions of predicted order. When compared with disordered stretches, ordered regions were not found to be enriched for asparagines, serines and threonines, residues that constitute the sequon signature for conjugation of N-glycans. We then investigated the basis of mutual exclusivity between disorder and N-glycosites on the basis of amino acid distribution: when compared with control ordered residue stretches without any N-glycosites, residue neighborhoods surrounding N-glycosites showed a depletion of bulky, hydrophobic and disorder-promoting amino acids and an enrichment for flexible and accessible residues that are frequently found in coiled structures. When compared with control disordered residue stretches without any N-glycosites, N-glycosite neighborhoods were depleted of charged, polar, hydrophobic and flexible residues and enriched for aromatic, accessible and order-promoting residues with a tendency to be part of coiled and β structures. N-glycosite neighborhoods also showed greater phylogenetic conservation among amniotes, compared with control ordered regions, which in turn were more conserved than disordered control regions. Our results lead us to propose that unique primary structural compositions and differential propensities for evolvability allowed for the mutual spatial exclusion of N-glycosite neighborhoods and disordered stretches.

Identification of a transforming growth factor-β type I receptor transcript in Eriocheir sinensis and its molting-related expression in muscle tissues

The transforming growth factor-β (TGF-β) signaling pathway is conserved across animals, and knowledge of its roles during the molt cycle in crustaceans is presently very limited. This study investigates the roles of the TGF-β receptor in molting-related muscle growth in Eriocheir sinensis. Using the RT-PCR and RACE techniques, we obtained a 1722 bp cDNA sequence encoding a transforming growth factor-β type I receptor in Eriocheir sinensis, designated EsTGFBRI, which contains a 124 bp 5'-untranslated region, a 20 bp partial 3'-untranslated region and a 1578 bp open reading frame encoding 525 amino acids. The deduced EsTGFBRI contains an N-terminal 24 amino acid signal peptide, an activin type I and II receptor domain, a transmembrane helix region, a glycine-serine-rich motif, and a conserved serine/threonine kinase catalytic domain including an activation loop. The qRT-PCR results showed that EsTGFBRI gene was highly expressed in the intermolt testis and ovary in mature crabs. In juvenile crabs, the mRNA levels of EsTGFBRI in claw and abdominal muscles in the later premolt D_3-4 stage were significantly higher than those in the intermolt C and postmolt A-B stages. There was no significant change in EsTGFBRI mRNA levels in walking leg muscles during the molt cycle. The results suggest that EsTGFBRI is probably play roles in molting-related muscle growth in E. sinensis. This study provides a necessary basis for elucidating the functions of TGF-β-like signaling mediated by TGFBRI in molting-related muscle growth in crustaceans.

microRNAs in Macrobrachium olfersii embryos: Identification, their biogenesis components and potential targets

In embryonic development, microRNAs (miRNAs) regulate the complex gene expression associated with the complexity of embryogenesis. Today, few studies have been conducted on the identification of miRNAs and components of miRNA biogenesis on embryonic development in crustaceans, especially in prawns. In this context, the aim of this study was to identify in silico components of miRNA biogenesis, and miRNAs and potential target genes during embryonic development in the prawn Macrobrachium olfersii through small RNAs and transcriptome analyses. Using the miRDeep2 program, we identified 17 miRNA precursors in M. olfersii, which seven (miR-9, miR-10, miR-92, miR-125, miR-305, miR-1175, and miR-2788) were reported in the miRBase database, indicating high evolutionary conservation of these sequences among animals. The other 10 miRNAs of M. olfersii were novel miRNAs and only similar to Macrobrachium niponnense miRNAs, indicating genus-specific miRNAs. In addition, eight key components of miRNA biogenesis (DROSHA, PASHA/DGCR8, XPO5, RAN, DICER, TRBP2, AGO, and PIWI) were identified in M. olfersii embryos unigenes. In the annotation of miRNA targets, 516 genes were similar to known sequences in the GenBank database. Regarding the conserved miRNAs, we verified that they were differentially expressed during embryonic development in M. olfersii. In conclusion, this is the first study that identifies conserved and novel miRNAs in the prawn M. olfersii with some miRNA target genes involved in embryonic development. Our results will allow further studies on the function of these miRNAs and miRNA biogenesis components during embryonic development in M. olfersii and other prawns of commercial interest.

Wnt gene family members and their expression profiling in Litopenaeus vannamei

The Wnt gene family encodes secreted glycoproteins involved in a wide variety of biological processes, including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt pathway exists in all metazoan animals, however, the relevant research is rare in crustaceans. Here we described 12 Wnt genes representing 12 Wnt gene subfamilies in the Pacific white shrimp, Litopenaeus vannamei. Based on homolog annotations and phylogenetic analyses, we named these 12 Wnt genes as LvWnt1, LvWnt2, LvWnt4-11, LvWnt16, and LvWntA. All the corresponding LvWnt proteins shared a conserved Wnt1 domain and 22 conserved cysteine residues. LvWnt1 and LvWnt6 were adjacent in a scaffold in the shrimp genome. Furthermore, we performed expression analyses of LvWnt genes at different developmental stages, during the molting process, in different tissues and after different pathogenic infection. We showed that each LvWnt gene had a unique expression pattern at different developmental stages but only a few of them expressed in adult shrimp. All the investigated LvWnt genes were initially expressed at the gastrula or limb bud embryo stages. Among them, LvWnt8 was specifically high expressed only in early embryos. LvWntA and LvWnt5 displayed high and similar expression profiles during the molting process, and LvWnt6 and LvWnt16 were specifically expressed in the thoracic ganglion, ventral nerve, intestines and gill tissues, respectively. We also found the expression of LvWntA, LvWnt5, LvWnt6, LvWnt9, and LvWnt16 were varied in the different tissues after infected with Staphylococcus aureus, Vibrio parahaemolyticus and white spot syndrome virus (WSSV), which indicated that they might participate in immune response in L. vannamei. This study provided an insight into the repertoire of the Wnt gene structure and expression in shrimps, and furthermore, might promote the understanding of development, growth and immune response of shrimps and crustaceans.

Effect of UVB radiation exposure in the expression of genes and proteins related to apoptosis in freshwater prawn embryos

Our previous studies showed that embryos of the freshwater prawn Macrobrachium olfersii exposed to ultraviolet B (UVB) radiation exhibited DNA damage, excessive ROS production, mitochondrial dysfunction and increased hsp70 expression, which are able, independently or together, to induce apoptosis. Thus, we attempted to elucidate some key apoptosis-related genes (ARG) and apoptosis-related proteins (ARP) and their expression during different stages of embryonic development, as well as to characterize the chronology of ARG expression and ARP contents after UVB radiation insult. We demonstrate that p53, Bax and Caspase3 genes are active in the embryonic cells at early embryonic developmental stages, and that the Bcl2 gene is active from the mid-embryonic stage. After UVB radiation exposure, we found an increase in ARP such as p53 and Bak after 3h of exposure. Moreover, an increase in ARG transcript levels for p53, Bax, Bcl2 and Caspase3 was observed at 6h after UVB exposure. Then, after 12h of UVB radiation exposure, an increase in Caspase3 gene expression and protein was observed, concomitantly with an increased number of apoptotic cells. Our data reveal that ARG and ARP are developmentally regulated in embryonic cells of M. olfersii and that UVB radiation causes apoptosis after 12h of exposure. Overall, we demonstrate that embryonic cells of M. olfersii are able to active the cell machinery against environmental changes, such as increased incidence of UVB radiation in aquatic ecosystems.