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.

CrusTome: a transcriptome database resource for large-scale analyses across Crustacea.

Transcriptomes from nontraditional model organisms often harbor a wealth of unexplored data. Examining these data sets can lead to clarity and novel insights in traditional systems, as well as to discoveries across a multitude of fields. Despite significant advances in DNA sequencing technologies and in their adoption, access to genomic and transcriptomic resources for nontraditional model organisms remains limited. Crustaceans, for example, being among the most numerous, diverse, and widely distributed taxa on the planet, often serve as excellent systems to address ecological, evolutionary, and organismal questions. While they are ubiquitously present across environments, and of economic and food security importance, they remain severely underrepresented in publicly available sequence databases. Here, we present CrusTome, a multispecies, multitissue, transcriptome database of 201 assembled mRNA transcriptomes (189 crustaceans, 30 of which were previously unpublished, and 12 ecdysozoans for phylogenetic context) as an evolving and publicly available resource. This database is suitable for evolutionary, ecological, and functional studies that employ genomic/transcriptomic techniques and data sets. CrusTome is presented in BLAST and DIAMOND formats, providing robust data sets for sequence similarity searches, orthology assignments, phylogenetic inference, etc. and thus allowing for straightforward incorporation into existing custom pipelines for high-throughput analyses. In addition, to illustrate the use and potential of CrusTome, we conducted phylogenetic analyses elucidating the identity and evolution of the cryptochrome/photolyase family of proteins across crustaceans.

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates.

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24ß1, -A24ß2, -A34α2, -A34ß1, and -A34ß2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded ß-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34ß1/ß2 had an additional two-stranded ß-sheet. We hypothesize that this second ß-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

Recommendations for Advancing Genome to Phenome Research in Non-Model Organisms.

The 2020 SICB Society-wide Symposium "Building Bridges from Genome to Phenome: Molecules, Methods and Models" brought together a diverse group of scientists to discuss recent progress in linking phenotype plasticity to changes at the level of the genome, epigenome, and proteome, while exploring the boundaries between variation and speciation. In a follow-up workshop, participants were asked to assess strengths and weaknesses of current approaches, to identify common barriers inhibiting their progress, and to outline the resources needed to overcome those barriers. Discussion groups generally recognized the absence of any overarching theoretical framework underlying current genome to phenome research and, therefore, called for a new emphasis on the development of conceptual models as well as the interdisciplinary collaborations needed to create and test those models. Participants also recognized a critical need for new and improved molecular and bioinformatic approaches to assist in describing function/phenotypes across phylogeny. Additionally, like all scientific endeavors, progress in genome to phenome research will be enhanced by improvements in science education and communication both within and among working groups.

Transcriptomic analysis of differentially expressed genes in the molting gland (Y-organ) of the blackback land crab, Gecarcinus lateralis, during molt-cycle stage transitions.

A transcriptome of the Gecarcinus lateralis molting gland (Y-organ or YO) contained 48,590 contiguous sequences (contigs) from intermolt (IM), early premolt (EP), mid premolt (MP), late premolt (LP), and postmolt (PM) stages. The YO is kept in the basal state in IM by molt-inhibiting hormone (MIH)/cyclic nucleotide-dependent signaling. YO activation in EP requires down-regulation of MIH signaling and activation of mechanistic target of rapamycin (mTOR)-dependent protein synthesis. Transition of the YO to the committed state in MP requires activin/transforming growth factor-beta (TGFß) signaling. YO repression occurs at the end of LP. A total of 28,179 contigs (58%) showed molt stage-specific changes in gene expression. The largest number of differentially-expressed genes (DEGs) were at the IM/EP (16,142 contigs), LP/PM (18,161 contigs), and PM/IM (8290 contigs) transitions. By contrast, the numbers of DEGs were 372 and 1502 contigs for the EP/MP and MP/LP transitions, respectively. DEG analysis of 23 signal transduction pathways showed significant changes in MIH, mTOR, activin/TGFß, Notch, MAP kinase, and Wnt signaling. Down-regulation of MIH signaling genes in premolt is consistent with reduced MIH sensitivity in MP and LP. Up-regulation of mTOR signaling genes in IM and premolt stages is consistent with its role in YO activation and sustained ecdysteroidogenesis. Up-regulation of activin/TGFß signaling genes in EP and MP is consistent with the role of a myostatin/activin-like factor in YO commitment. Notch, MAP kinase, and Wnt DEG analysis may indicate possible crosstalk with the MIH, mTOR, and activin/TGFß pathways to integrate other inputs to control YO ecdysteroidogenesis.

Analysis of Annotation and Differential Expression Methods used in RNA-seq Studies in Crustacean Systems.

In the field of crustacean biology, usage of RNA-seq to study gene expression is rapidly growing. Major advances in sequencing technology have contributed to the ability to examine complex patterns of genome activity in a wide range of organisms that are extensively used for comparative physiology, ecology and evolution, environmental monitoring, and commercial aquaculture. Relative to insect and vertebrate model organisms, however, information on the organization of crustacean genomes is virtually nonexistent, making de novo transcriptome assembly, annotation and quantification problematic and challenging. We present here a summary of the methodologies and software analyses employed in 23 recent publications, which describe de novo transcriptome assembly, annotation, and differential gene expression in a variety of crustacean experimental systems. We focus on establishing a series of best practices that will allow for investigators to produce datasets that are understandable, reproducible, and of general utility for related analyses and cross-study comparisons.

Resources and Recommendations for Using Transcriptomics to Address Grand Challenges in Comparative Biology.

High-throughput RNA sequencing (RNA-seq) technology has become an important tool for studying physiological responses of organisms to changes in their environment. De novo assembly of RNA-seq data has allowed researchers to create a comprehensive catalog of genes expressed in a tissue and to quantify their expression without a complete genome sequence. The contributions from the "Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology" symposium in this issue show the successes and limitations of using RNA-seq in the study of crustaceans. In conjunction with the symposium, the Animal Genome to Phenome Research Coordination Network collated comments from participants at the meeting regarding the challenges encountered when using transcriptomics in their research. Input came from novices and experts ranging from graduate students to principal investigators. Many were unaware of the bioinformatics analysis resources currently available on the CyVerse platform. Our analysis of community responses led to three recommendations for advancing the field: (1) integration of genomic and RNA-seq sequence assemblies for crustacean gene annotation and comparative expression; (2) development of methodologies for the functional analysis of genes; and (3) information and training exchange among laboratories for transmission of best practices. The field lacks the methods for manipulating tissue-specific gene expression. The decapod crustacean research community should consider the cherry shrimp, Neocaridina denticulata, as a decapod model for the application of transgenic tools for functional genomics. This would require a multi-investigator effort.

Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology: An Introduction to the Symposium.

Crustaceans, and decapods in particular (i.e., crabs, shrimp, and lobsters), are a diverse and ecologically and commercially important group of organisms. Understanding responses to abiotic and biotic factors is critical for developing best practices in aquaculture and assessing the effects of changing environments on the biology of these important animals. A relatively small number of decapod crustacean species have been intensively studied at the molecular level; the availability, experimental tractability, and economic relevance factor into the selection of a particular species as a model. Transcriptomics, using high-throughput next generation sequencing (NGS, coupled with RNA sequencing or RNA-seq) is revolutionizing crustacean biology. The 11 symposium papers in this volume illustrate how RNA-seq is being used to study stress response, molting and limb regeneration, immunity and disease, reproduction and development, neurobiology, and ecology and evolution. This symposium occurred on the 10th anniversary of the symposium, "Genomic and Proteomic Approaches to Crustacean Biology", held at the Society for Integrative and Comparative Biology 2006 meeting. Two participants in the 2006 symposium, the late Paul Gross and David Towle, were recognized as leaders who pioneered the use of molecular techniques that would ultimately foster the transcriptomics research reviewed in this volume. RNA-seq is a powerful tool for hypothesis-driven research, as well as an engine for discovery. It has eclipsed the technologies available in 2006, such as microarrays, expressed sequence tags, and subtractive hybridization screening, as the millions of "reads" from NGS enable researchers to de novo assemble a comprehensive transcriptome without a complete genome sequence. The symposium series concludes with a policy paper that gives an overview of the resources available and makes recommendations for developing better tools for functional annotation and pathway and network analysis in organisms in which the genome is not available or is incomplete.

Transcriptome analysis of the molting gland (Y-organ) from the blackback land crab, Gecarcinus lateralis.

In decapod crustaceans, arthropod steroid hormones or ecdysteroids regulate molting. These hormones are synthesized and released from a pair of molting glands called the Y-organs (YO). Cyclic nucleotide, mTOR, and TGFß/Smad signaling pathways mediate molt cycle-dependent phase transitions in the YO. To further identify the genes involved in the regulation of molting, a YO transcriptome was generated from three biological replicates of intermolt blackback land crab, Gecarcinus lateralis. Illumina sequencing of cDNA libraries generated 227,811,829 100-base pair (bp) paired-end reads; following trimming, 90% of the reads were used for further analyses. The trimmed reads were assembled de novo using Trinity software to generate 288,673 contigs with a mean length of 872 bp and a median length of 1842 bp. Redundancy among contig sequences was reduced by CD-HIT-EST, and the output constituted the baseline transcriptome database. Using Bowtie2, 92% to 93% of the reads were mapped back to the transcriptome. Individual contigs were annotated using BLAST, HMMER, TMHMM, SignalP, and Trinotate, resulting in assignments of 20% of the contigs. Functional and pathway annotations were carried out via gene ontology (GO) and KEGG orthology (KO) analyses; 58% and 44% of the contigs with BLASTx hits were assigned to GO and KO terms, respectively. The gene expression profile was similar to a crayfish YO transcriptome database, and the relative abundance of each contig was highly correlated among the three G. lateralis replicates. Signal transduction pathway orthologs were well represented, including those in the mTOR, TGFß, cyclic nucleotide, MAP kinase, calcium, VEGF, phosphatidylinositol, ErbB, Wnt, Hedgehog, Jak-STAT, and Notch pathways.

Alternative splicing in the fiddler crab cognate ecdysteroid receptor: variation in receptor isoform expression and DNA binding properties in response to hormone.

RXR cDNA cloning from three Uca species led to the identification of 4 conserved isoforms, indicative of alternative splicing in the hinge and ligand binding domains (LBD). Sequencing of overlapping clones from a Ucapugilator genomic library identified EcR isoforms matching previously identified cDNA variants; in addition, a cryptic exon in the LBD was detected and evidence for expression of this new isoform was obtained from next-generation sequencing. RNA-seq analysis also identified a new amino terminal EcR variant. EcR and RXR transcript abundance increases throughout ovarian maturation in U. pugilator, while cognate receptor transcript abundance remains constant in a related Indo-Pacific species with a different reproductive strategy. To examine if crab RXR LBD isoforms have different physical properties in vitro, electromobility shift assays were performed with different EcR isoforms. The cognate crab and fruit fly receptors differ in their responses to hormone. Ecdysteroids did not increase DNA binding for the crab heterodimers, while ecdysteroids stimulate binding for Drosophilamelanogaster EcR/USP heterodimers. In swapping experiments, UpEcR/USP heterodimers did not show ligand-responsive differences in DNA binding; both crab RXR LBD isoforms, however, conferred ligand-responsive increases in DNA binding with DmEcRs. These data indicate that both UpRXR LBD isoforms can heterodimerize with the heterologous DmEcR receptors and promote ligand and DNA binding. Unresponsiveness of the cognate receptors to ecdysteroid, however, suggest additional factors may be required to mediate endogenous, perhaps isoform-specific, differences in EcR conformation, consistent with previously reported effects of UpRXR isoforms on UpEcR ligand-binding affinities.

Ecdysteroid receptor signaling disruption obstructs blastemal cell proliferation during limb regeneration in the fiddler crab, Uca pugilator.

To study ecdysteroid signaling during limb regeneration, we have applied RNAi (dsRNA) mediated silencing to EcR/RXR, the genes encoding the ecdysteroid receptor heterodimer, in the fiddler crab Uca pugilator. We injected RNAi into the blastemal chamber during early limb regeneration. Silencing was evaluated by knockdown in receptor transcript abundance, and disruption was evaluated by changes in growth rate and morphology of limb regenerates. q-PCR results indicated a 50% drop in transcript abundance 48h post injection in both RNAi (dsEcR/dsRXR) injected ipsilateral and uninjected contralateral blastemas in experimental animals relative to controls. EcR/RXR transcript levels further decreased over time. Several phenotypes were associated with knockdown. The experimental blastema failed to develop; microscopic examination of the arrested blastema revealed an absence of the cuticular ingrowths characteristic of the beginnings of limb segmentation and cell proliferation assays revealed that the arrested blastema had few dividing cells. Ecdysteroid levels were also lowered in experimental animals; given the bilateral effects of RNAi on limb buds in experimental animals, these results suggest RNAi had a systemic effect. Although hormone titers in experimental animals rose to comparable control levels during the late proecdysial phase of limb regeneration, most experimental crabs failed to molt and died. The overall failure to molt indicates that RNAi receptor knockdown has long-term effects. The combined effects of receptor knockdown indicate that, although circulating ecdysteroid titers are normally low during basal limb bud growth, signaling via the ecdysteroid receptor pathway is necessary for establishment of blastemal cell proliferation and development in the regenerating limbs of U. pugilator.

EST library sequencing of genes expressed during early limb regeneration in the fiddler crab and transcriptional responses to ecdysteroid exposure in limb bud explants.

We have constructed directional and randomly primed cDNA libraries from mRNAs isolated during progressive stages of fiddler crab (Uca pugilator) limb regeneration. Data from these libraries are being assembled into an on-line database (http://www.genome.ou.edu/crab.html) that is both BLAST and keyword searchable; the data set is also available through GenBank. The first characterized library was made from mRNA isolated 4 days post-autotomy, when the first sign of morphological differentiation, cuticle secretion, is observed. Analysis of 1698 cDNA clones led to assignment of 473 contigs and 417 singlets, for a total of 890 sequences. Of these, ∼86% showed no assignments to characterized genes on database searching, while 14% could be assigned to a known ortholog in the COG (Clusters of Orthologous Groups) database. BLAST searches to specific protein domains in the Gene Ontology database led to assignments for ∼40% of the assembled sequences. Sequence similarity searches of other crustacean EST databases produced hits to 13-30% of the Uca query sequences. The ESTs include several genes that may be potentially ecdysteroid-responsive, such as homologs to chaperone proteins and cuticle protein genes, as well as homologs to arthropod proteins involved in retinoid/terpenoid metabolism. We have tested 3 potential candidate genes for their ability to be induced by ecdysteroid in limb bud explants; an arthropodial cuticle protein gene, and the nuclear receptor genes EcR and RXR. A subset of early blastemal limb buds (8 days post autotomy) show a positive response to ecdysteroid by 1-1.5 h, followed by a decrease in transcript abundance at longer periods of sustained incubation. Later stage buds (12 days post autotomy-late premolt) show decreases in steady-state mRNA levels by 1.5 h, or are completely refractory to ecdysteroid exposure.

Crustacean retinoid-X receptor isoforms: distinctive DNA binding and receptor-receptor interaction with a cognate ecdysteroid receptor.

We have identified cDNA clones that encode homologs of the ecdysteroid receptor (EcR) and retinoid-X receptor (RXR)/USP classes of nuclear receptors from the fiddler crab Uca pugilator (UpEcR and UpRXR). Several UpRXR cDNA splicing variants were found in coding regions that could potentially influence function. A five-amino acid (aa) insertion/deletion is located in the "T" box in the hinge region. Another 33-aa insertion/deletion is found inside the ligand-binding domain (LBD), between helix 1 and helix 3. Ribonuclease protection assays (RPA) showed that four UpRXR transcripts [UpRXR(+5+33), UpRXR(-5+33), UpRXR(+5-33) and UpRXR(-5-33)] were present in regenerating limb buds. UpRXR(-5+33) was the most abundant transcript present in regenerating limb buds in both early blastema and late premolt growth stages. Expression vectors for these UpRXR variants and UpEcR were constructed, and the proteins expressed in E. coli and in vitro expression systems. The expressed crab nuclear receptors were then characterized by electrophoretic mobility shift assay (EMSA) and glutathione S-transferase (GST) pull down experiments. EMSA results showed that UpEcR/UpRXR(-5+33) heterocomplexes bound with a series of hormone response elements (HREs) including eip28/29, IRper-1, DR-4, and IRhsp-1 with appreciable affinity. Competition EMSA also showed that the affinity decreased as sequence composition deviated from a perfect consensus element. Binding to IRper-1 HREs occurred only if the heterodimer partner UpRXR contained the 33-aa LBD insertion. UpRXR lacking both the 5-aa and 33-aa insertion bound to a DR-1G HRE in the absence of UpEcR. The results of GST-pull down experiments showed that UpEcR interacted only with UpRXR variants containing the 33-aa insertion, and not with those lacking the 33-aa insertion. These in vitro receptor protein-DNA and receptor protein-protein interactions occurred in the absence of hormone (20-hydroxyecdysone and 9-cis retinoid acid, 9-cis RA). Transactivation studies using a hybrid UpEcR ligand-binding domain construct and UpRXR (+/-33) ligand-binding domain constructs also showed that the 33-aa insertion was indispensable in mediating ecdysteroid stimulated transactivation.

Characterization of crab EcR and RXR homologs and expression during limb regeneration and oocyte maturation.

We report here complete coding sequences for the Uca pugilator homologs of the ecdysteroid (UpEcR) and retinoid-X receptors (UpRXR). Library screenings recovered cDNA clones containing a unique amino terminal open-reading frame (A/B domain) for each gene, most similar to insect B1 EcR and USP1/RXR isoforms. Splicing variants in the UpRXR ligand-binding domain were also identified, in a region critical for folding of Drosophila and lepidopteran USP. UpEcR and UpRXR proteins were able to associate, and both are required for binding to an ecdysteroid HRE; these interactions were not hormone-dependent. Ribonuclease protection assays (RPA) were conducted using A/B domain and 'common' (C or E) domain probes on RNA isolated from various stages of regenerating limb buds and ovaries. For several of the limb bud and ovarian stages examined, the relative level of A/B domain sequence protected was significantly less than common domain suggesting alternative amino terminal isoforms other than those isolated through cloning. This is the first report of UpEcR and UpRXR transcription during ovarian maturation, implicating the ovary as a potential target for hormonal control in Crustacea.