The evolutionary adaptation of shrimp hemocyanin subtypes and the consequences on their structure and functions

Hemocyanin is the main respiratory protein of arthropods and is formed by hexameric and/or oligomeric subunits. Due to changes in the living environment and gene rearrangement, various hemocyanin subtypes and subunits evolved in crustaceans. This paper reviews the various hemocyanin subtypes and isoforms in shrimp and analyses published genomic data of sixteen hemocyanin family genes from Litopenaeus vannamei to explore the evolution of hemocyanin genes, subunits, and protein structure. Analysis of hemocyanin subtypes distribution and structure in various tissues was also performed and related to multiple and tissue-specific functions, i.e., immunological activity, immune signaling, phenoloxidase activity, modulation of microbiota homeostasis, and energy metabolism. The functional diversity of shrimp hemocyanin due to molecular polymorphism, transcriptional regulation, alternative splicing, degradation into functional peptides, interaction with other proteins or genes, and structural differences will also be highlighted for future research. Inferences would be drawn from other crustaceans to explain how evolution has changed the structure-function of hemocyanin and its implication for evolutionary research into the multifunctionality of hemocyanin and other related proteins in shrimp.

The pyruvate kinase of the whiteleg shrimp Litopenaeus vannamei: Gene structure and responses to short term hypoxia

The shrimp Litopenaeus vannamei is the main farmed crustaceans worldwide. This crustacean suffers environmental changes in oxygen availability that affect its energy metabolism. Pyruvate kinase (PK) catalyzes the last reaction of glycolysis and is key for the regulation of glycolysis and gluconeogenesis. There is ample knowledge about mammalian PK, but in crustaceans, the information is very scarce. In this study, we analyzed in silico the structures of the PK gene and protein. Also, the effects of hypoxia on gene expression, enzymatic activity, glucose, and lactate in hepatopancreas and muscle were analyzed. The PK gene is 15,103 bp and contains 11 exons and 10 introns, producing four mRNA variants by alternative splicing and named PK1, PK2, PK3 and PK4, and two proteins with longer C-terminus and two with a 12 bp insertion. The promoter contains putative binding sites for transcription factors (TF) that are typically involved in stress responses. The deduced amino acid sequences contain the classic domains, binding sites for allosteric effectors and potential reversible phosphorylation residues. Protein modeling indicates a homotetramer with highly conserved structure. The effect of hypoxia for 6 and 12 h showed tissue-specific patterns, with higher expression, enzyme activity and lactate in muscle, but higher glucose in hepatopancreas. Changes in response to hypoxia were detected at 12 h in expression with induction in muscle and reduction in hepatopancreas, while enzyme activity was maintained, and glucose and lactate decreased. These results show rapid changes in expression and metabolites, while enzyme activity was maintained to cope with short-term hypoxia.

Molecular isolation and expression analysis of hemocyanin isoform 2 of Macrobrachium rosenbergii

OBJECTIVE

Hemocyanin is a copper-bearing protein in the hemolymph of many arthropods and mollusks and functions as an oxygen transport and important nonspecific immune protein.

METHODS

In this study, complementary DNA of hemocyanin isoform 2 of the prawn Macrobrachium rosenbergii (MrHc2) was isolated by rapid amplification of cDNA ends and mRNA expression was characterized to elucidate molecular basis of its function.

RESULT

With a molecular mass of 77.3 kDa, MrHc2 contained three domains: hemocyanin-all-alpha, hemocyanin-copper-containing, and hemocyanin-immunoglobulin-like domains. Molecular phylogenetic analysis revealed that MrHc2 belongs to the γ-type subunit and is closely related to hemocyanin subunit 1 of the palaemonid shrimp Macrobrachium nipponense. In addition, MrHc2 resided in a different clade relative to hemocyanin (MrHc) of M. rosenbergii (α-type subunit) and in a different subclade relative to the hemocyanin proteins of penaeid shrimp. The messenger RNA transcript of MrHc2 was highly expressed in the hepatopancreas and weakly expressed in the gills, intestine, stomach, muscle, and hemocytes. Upon challenge with M. rosenbergii nodavirus (MrNV), the expression of MrHc2 was 1.96-, 2.93-, and 1.96-fold on days 3, 4, and 5, respectively, and then gradually declined to basal levels on day 7.

CONCLUSION

This study suggests that MrHc2 plays an important role in the innate immune response of M. rosenbergii to MrNV.

Protein Diversity and Immune Specificity of Hemocyanin From Shrimp Litopenaeus vannamei

Hemocyanin is an important non-specific innate immune defense molecule with phenoloxidase, antiviral, antibacterial, hemolytic, and antitumor activities. To better understand the mechanism of functional diversity, proteomics approach was applied to characterize hemocyanin (HMC) expression profiles from Litopenaeus vannamei. At first, hemocyanin was purified by Sephadex G-100 and DEAE-cellulose (DE-52) columns from shrimp serum, and 34 protein spots were identified as HMC on the 2-DE gels. Furthermore, we found that 9 HMC spots about 75 or 77 kDa were regulated by Streptococcus agalactiae and Vibrio parahaemolyticus infection at 6, 12, and 24 h. In addition, 6 different pathogen-binding HMC fractions, viz., HMC-Mix, HMC-Vp, HMC-Va, HMC-Vf, HMC-Ec, and HMC-Sa, showed different agglutinative and antibacterial activities. Moreover, lectin-blotting analysis showed significant differences in glycosylation level among HMC isomers and bacteria-binding HMC fractions. Particularly, the agglutinative activities of the HMC fractions were almost completely abolished when HMC was deglycosylated by O-glycosidase, which suggest that O-linked sugar chains of HMC played important roles in the innate immune recognition. Our findings demonstrated for the first time that L. vannamei HMC had molecular diversity in protein level, which is closely associated with its ability to recognize diverse pathogens, whereas glycan modification probably contributed to HMC’s diversity and multiple immune activities.

Resolving hemocyanin isoform complexity in haemolymph of black tiger shrimp Penaeus monodon - implications in aquaculture, medicine and food safety

Hemocyanin (Hc) is a multifunctional macromolecule involved in oxygen transport and non-specific immunity in shrimp. Hc is crucial in physiology and nutrition linked with optimal performance in aquaculture production systems. In medicine, Hc has been approved for clinical use in humans as adjuvant and anticancer therapeutic. In contrast, Hc has also been identified as one of the proteins causing anaphylaxis following shrimp consumption. The role of individual Hc isoforms remains unknown due to a lack of resolved Hc isoforms. We successfully identified eleven different Penaeus monodon hemocyanin (PmoHc) γ isoforms including two truncated isoforms (50 and 20 kDa) and one PmoHc β isoform in haemolymph using proteomics informed by transcriptomics. Amino acid sequence homology ranged from 24 to 97% between putative PmoHc gene isoforms. Hc isoforms showed specific patterns of transcript expression in shrimp larval stages and adult hepatopancreas. These findings enable isoform level investigations aiming to define molecular mechanisms underpinning Hc functionality in shrimp physiology and immunity, as well as their individual immunogenic role in human allergy. Our research demonstrates the power of proteomics informed by transcriptomics to resolve isoform complexity in non-model organisms and lay the foundations for improved performance within the aquaculture industry and advance allergenic applications in medicine. SIGNIFICANCE: The roles of hemocyanin (Hc) in shrimp homeostasis and immunity as well as in human allergy are not well understood because the complexity of Hc isoforms has remained unresolved. Our results have confirmed the existence of at least 12 individual Hc isoforms in shrimp haemolymph and validated putative Hc gene assemblies from transcriptomics. Our findings will enable monitoring the expression of specific Hc isoforms in shrimp haemolymph during different environmental, nutritional and pathogenic conditions, thus providing insights into isoform specific functional roles. In medicine, the potential allergenicity of each Hc isoform could be determined and advance allergenic applications. Lastly, since Hc comprises up to 95% of the total protein in haemolymph, these isoforms become ideal targets for prawn provenance, traceability and food contamination studies.

Moderate reductions in dissolved oxygen may compromise performance in an ecologically-important estuarine invertebrate

Coastal ecosystems, including estuaries, are increasingly pressured by expanding hypoxic regions as a result of human activities such as increased release of nutrients and global warming. Hypoxia is often defined as oxygen concentrations below 2 mL O₂ L^-1. However, taxa vary markedly in their sensitivity to hypoxia and can be affected by a broad spectrum of low oxygen levels. To better understand how reduced oxygen availability impacts physiological and molecular processes in invertebrates, we investigated responses of an estuarine amphipod to an ecologically-relevant level of moderate hypoxia (~2.6 mL O₂ L^-1) or severe hypoxia (~1.3 mL O₂ L^-1). Moderate hypoxia elicited a reduction in aerobic scope, and widespread changes to gene expression, including upregulation of metabolic genes and stress proteins. Under severe hypoxia, a marked hyperventilatory response associated with maintenance of aerobic performance was accompanied by a muted transcriptional response. This included a return of metabolic genes to baseline levels of expression and downregulation of transcripts involved in protein synthesis, most of which indicate recourse to hypometabolism and/or physiological impairment. We conclude that adverse ecological effects may occur under moderate hypoxia through compromised individual performance and, therefore, even modest declines in future oxygen levels may pose a significant challenge to coastal ecosystems.

Comparative proteomics in the wild: Accounting for intrapopulation variability improves describing proteome response in a Gammarus pulex field population exposed to cadmium

High-throughput proteomics can be performed on animal sentinels for discovering key molecular biomarkers signing the physiological response and adaptation of organisms. Ecotoxicoproteomics is today amenable by means of proteogenomics to small arthropods such as Gammarids which are well known sentinels of aquatic environments. Here, we analysed two regional Gammarus pulex populations to characterize the potential proteome divergence induced in one site by natural bioavailable mono-metallic contamination (cadmium) compared to a non-contaminated site. Two RNAseq-derived protein sequence databases were established previously on male and female individuals sampled from the reference site. Here, individual proteomes were acquired on 10 male and 10 female paired organisms sampled from each site. Proteins involved in protein lipidation, carbohydrate metabolism, proteolysis, innate immunity, oxidative stress response and lipid transport were found more abundant in animals exposed to cadmium, while hemocyanins were found in lower abundance. The intrapopulation proteome variability of long-term exposed G. pulex was inflated relatively to the non-contaminated population. These results show that, while remaining a challenge for such organisms with not yet sequenced genomes, taking into account intrapopulation variability is important to better define the molecular players induced by toxic stress in a comparative field proteomics approach.

Evolutionary adaptation analysis of immune defense and hypoxia tolerance in two closely related Marsupenaeus species based on comparative transcriptomics

Kuruma shrimp, a major farmed shrimp species in the world, includes two cryptic or sibling species, Form I (Marsupenaeus japonicus) and Form II (Marsupenaeus pulchricaudatus). Due to the lack of genomic resources, little is known about the molecular mechanisms associated with immune defense and hypoxia tolerance. Here, we sequenced the transcriptomes of two closely related Marsupenaeus species and compared genomic divergence. This study obtained 77049 and 84561 unigenes with N50 values of 1281bp and 1244bp for M. japonicus and M. pulchricaudatus, respectively, and 5036 pairs of putative orthologs were identified between two Marsupenaeus species. Estimation of Ka/Ks ratios indicated that 165 orthologous genes may be under positive selection (Ka/Ks > 0.5), including 49 pairs with a Ka/Ks ratio >1. According to the peak of synonymous rates, the divergence time between M. japonicus and M. pulchricaudatus was about 0.26-0.69 Mya. These positively selected orthologous genes related to the immune process mainly comprised single VWC domain protein, legumain, ras-related C3 botulinum, caspase, C-type lectin and were enriched in functions related to immune (Toll-like receptor and PI3K-Akt signaling) and hypoxia signaling (HIF-1 signaling and VEGF signaling). In this study, dozens of caspase-like unigenes were screened from two Marsupenaeus transcriptomes. Among these, the PjCaspase orthologous gene was subjected to positive selection (Ka/Ks = 1.22), which had different secondary and three-dimensional structure prediction. Based on the single copy caspase gene, eight populations of Marsupenaeus species were divided into two phylogeographic lineages from the East and South China. We characterized the transcriptomes of the two Marsupenaeus species and obtained several key orthologs associated with immune defense and hypoxia tolerance, which provides new insights into the immunity and genetic divergence of the two varieties. Moreover, this study will facilitate further comparative genomic studies of the two varieties.

Shrimp humoral responses against pathogens: antimicrobial peptides and melanization

Diseases have caused tremendous economic losses and become the major problem threatening the sustainable development of shrimp aquaculture. The knowledge of host defense mechanisms against invading pathogens is essential for the implementation of efficient strategies to prevent disease outbreaks. Like other invertebrates, shrimp rely on the innate immune system to defend themselves against a range of microbes by recognizing and destroying them through cellular and humoral immune responses. Detection of microbial pathogens triggers the signal transduction pathways including the NF-κB signaling, Toll and Imd pathways, resulting in the activation of genes involved in host defense responses. In this review, we update the discovery of components of the Toll and Imd pathways in shrimp and their participation in the regulation of shrimp antimicrobial peptide (AMP) synthesis. We also focus on a recent progress on the two most powerful and the best-studied shrimp humoral responses: AMPs and melanization. Shrimp AMPs are mainly cationic peptides with sequence diversity which endues them the broad range of activities against microorganisms. Melanization, regulated by the prophenoloxidase activating cascade, also plays a crucial role in killing and sequestration of invading pathogens. The progress and emerging research on mechanisms and functional characterization of components of these two indispensable humoral responses in shrimp immunity are summarized and discussed. Interestingly, the pattern recognition protein (PRP) crosstalk is evidenced between the proPO activating cascade and the AMP synthesis pathways in shrimp, which enables the innate immune system to build up efficient immune responses.

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.