Isolation and expression analysis of an MAPKK gene from Fenneropenaeus chinensis in response to white spot syndrome virus infection

Genome-wide identification and analysis of the MAPKK gene family in Chinese mitten crab (Eriocheir sinensis) and its response to bacterial challenge

Protein kinases of the MAPK cascade family (MAPKKK-MAPKK-MAPK) play an important role in the growth and development of organisms and their response to environmental stress. The MAPKK gene families in the Chinese mitten crab Eriocheir sinensis have never been systematically analyzed. We identified four MAPKK family genes, EsMEK, EsMAPKK4, EsMAPKK6, and EsMAPKK7, in E. sinensis and analyzed their molecular features and expression patterns. All four MAPKK genes are composed of multiple exons and introns, all have a conserved domain, and all have 10 conserved motifs (except EsMEK and EsMAPKK7 which are missing motif 10). The four MAPKK genes are on four different chromosomes and have no gene duplications, and the results of phylogenetic tree analysis indicate that the ESMAPKK gene family is highly conserved evolutionarily. The EsMAPKK genes were widely expressed in all the examined tissues with higher expression in hemocytes, hepatopancreas, and gills. Notably, EsMAPKK6 was also highly expressed in the ovary. Vibrio parahaemolyticus infection significantly increased the mRNA levels of the EsMAPKK genes in hemocytes. Further disruption of the EsMAPKK gene family expression affects the expression levels of multiple antimicrobial peptides in hemocytes. Our experimental results provide a starting point for a more in-depth study of the innate immunity functional roles of members of the MAPKK gene families in E. sinensis.

MEK homologue is involved in immune response by regulating antimicrobial peptides expression in Chinese mitten crab Eriocheir sinensis

MEK activates the phosphorylation of downstream molecules involved in various immune responses. In this study, an MEK homologue gene in Chinese mitten crab Eriocheir sinensis (designated as EsMEK) was investigated. EsMEK mRNA was constitutively expressed in all tissues with higher expression in hepatopancreas, hemocytes, and gills. EsMEK protein was mainly localized in the cytoplasm. Lipopolysaccharide (LPS) and Aeromonas hydrophila challenge significantly increased the mRNA levels of EsMEK in hemocytes. In addition, the mRNA expression level of some antimicrobial peptides (AMPs), including EsWAP, EsDWD1, and EsALF decreased significantly due to the inhibition of EsMEK by specific dsRNA in LPS-challenged crabs. Downstream pathway analysis revealed that the phosphorylation of EsERK decreased prominently after EsMEK inhibition. These results suggested that EsMEK played an important role in regulating the expression of antimicrobial peptides in E. sinensis through MEK-ERK pathway.

Ovarian transcriptome and metabolic responses of RNAi-mediated farnesyl pyrophosphate synthase knockdown in Neocaridina denticulata sinensis

Methyl farnesoate (MF) is considered the equivalent of JH in crustaceans and plays an essential role in many crucial physiological processes. It is believed that farnesyl pyrophosphate synthase (FPPS) plays an essential role in the biosynthesis of mevalonate, which is a branch of the JH/MF pathway. The full-length cDNA of FPPS (NdFPPS) from Neocaridina denticulata sinensis was isolated and characterized, and the deduced amino acid of NdFPPS contained a polyprenyl_synt domain. In addition to its ubiquitous tissue expression, NdFPPS was significantly expressed in the ovary. In vivo gene silencing with dsRNA interference was performed to further investigate the function of NdFPPS. An ovarian transcriptomic analysis of dsNdFPPS experimental and control groups was used to compare, annotate, and classify differentially expressed genes (DEGs). A total of 9230 DEGs were identified in the experimental and control groups based on FPKM values, of which 5082 were up-regulated genes and 4148 were down-regulated genes. 761 GO terms and 102 KEGG pathways were enriched for the DEGs. Our results suggest that NdFPPS might play an important role in ovarian development, however, further functional study is needed to elucidate physiological role of NdFPPS in reproduction.

Genome-wide analysis of the hard clam mitogen-activated protein kinase kinase gene family and their transcriptional profiles under abiotic stress

Mitogen-activated protein kinase kinase (MAPKK) was the hub component of the Mitogen-activated protein kinase (MAPK) signaling pathway and played an important role in the cellular response to environmental stress. In this study, we identified five MmMAPKK genes in hard clam Mercenaria mercenaria and found that all MmMAPKK genes contain a conserved protein kinase domain. The MmMAPKK genes derived from dispersed duplication were unevenly distributed in three chromosomes. Although the genome size was highly variable among different bivalve mollusks, the number of MAPKK genes was relatively stable. Phylogenetic analysis showed that bivalve MAPKK was divided into five clades, and amino acid sequences of MAPKK from the same clade consisted of similar conserved motifs. The syntenic analysis demonstrated that MmMAPKKs had the highest number of homologous gene pairs with Cyclina sinensis. MmMAPKKs were ubiquitously expressed in all examined tissues, and all MmMAPKK genes were highly expressed in the ovary. MmMAPKK genes showed stress-specific expression under envirionmental stress. MmMAPKK7 showed an upregulated in heat and heat plus hypoxia stress while MmMAPKK1 showed an upregulated in hypoxic stress groups. Dynamic changes of MmMAPKK7, MmMAPKK6 and MmMAPKK1 in hemocytes were observed in response to air exposure. MmMAPKK4 significantly downregulated after air exposure for five days. MmMAPKK7 and MmMAPKK6 might participate in adaptation to low salinity stress. Our results provided useful information about MAPKK and laid a foundation for further studies on MAPKK evolution in the bivalve.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

c-Raf participates in adaptive immune response of Nile tilapia via regulating lymphocyte activation

RAF proto-oncogene serine/threonine-protein kinase (c-Raf) is a MAP kinase kinase kinase (MAPKKK) that participates in the Erk1/2 pathway and plays an important role in lymphocyte activation. However, the study on how c-Raf regulates adaptive immunity in non-mammal is still limited. In present study, based on analysis of sequence characteristics of c-Raf from Oreochromis niloticus (On-c-Raf), we investigated its regulation roles on teleost lymphocyte activation. The On-c-Raf was highly conserved during evolution, which was composed of a Raf-like Ras-binding domain (RBD), a protein kinase C conserved region 1 (C1) domain and a serine/threonine protein kinase catalytic (S_TKc) domain. Its mRNA showed a wide distribution in tissues of O. niloticus and with the highest expression in gill. After Aeromonas hydrophila infection, during the adaptive immune stage transcription level of On-c-Raf was significantly upregulated on day 8, but came back to original level on day 16 and 30, suggesting the potential involvement of On-c-Raf in primary response but not memory formation. Furthermore, On-c-Raf mRNA in leukocytes of Nile tilapias was obviously induced by in vitro stimulation of T cell mitogen PHA. More importantly, in vitro stimulation of lymphocytes agonist PMA augmented phosphorylation level of On-c-Raf in leukocytes detected by western-blot and immunofluorescent. Thus, c-Raf regulated lymphocyte activation of Nile tilapia on both mRNA and phosphorylation level. Together, our results revealed that the c-Raf from teleost Nile tilapia engaged in adaptive immune response by regulating lymphocytes activation. Since the regulatory mechanism of lymphocyte-mediated adaptive immunity is largely unknown in teleost, our study provided important evidences to understand teleost adaptive immunity, and also shed a novel perspective for the evolution of adaptive immune system.

Identification of MicroRNAs and Their Target Genes Associated with Ovarian Development in Black Tiger Shrimp (Penaeus monodon) Using High-Throughput Sequencing

Plenty of evidence showing that microRNAs (miRNAs) post-transcriptionally regulate gene expression and are involved in a wide range of biological processes. However, the roles of miRNAs in ovarian development process remain largely unknown in shrimp. In the present study, high-throughput sequencing of small RNAs was performed to find specific miRNAs that are involved in ovarian development process in Penaeus monodon. Two small RNA libraries were constructed from undeveloped (UNDEV group) and developed (DEV group) ovarian tissues in P. monodon. In total, 43 differentially expressed miRNAs were identified between the two groups (P ≤ 0.05, |log2 ratio| ≥1), and their expression profiles were validated by qRT-PCR. In order to further clarify the functional roles of these differentially expressed miRNAs during ovarian development process, target gene prediction was performed. In total, 4,102 target genes of 43 miRNAs were predicted, then clustered by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; only four specific pathways related to ovarian development were obtained (P < 0.05). Dual-luciferase reporter assays and integrated expression analysis were also conducted to further clarify the interaction between the miRNAs and their target mRNAs. This study provides important information about the function of miRNAs involved in ovarian developmental stages in P. monodon.