De novo assembly of the blunt snout bream (Megalobrama amblycephala) gill transcriptome to identify ammonia exposure associated microRNAs and their targets

Gill Transcriptome Sequencing and De Novo Annotation of Acanthogobius ommaturus in Response to Salinity Stress

Acanthogobius ommaturus is a euryhaline fish widely distributed in coastal, bay and estuarine areas, showing a strong tolerance to salinity. In order to understand the mechanism of adaptation to salinity stress, RNA-seq was used to compare the transcriptome responses of Acanthogobius ommaturus to the changes of salinity. Four salinity gradients, 0 psu, 15 psu (control), 30 psu and 45 psu were set to conduct the experiment. In total, 131,225 unigenes were obtained from the gill tissue of A. ommaturus using the Illumina HiSeq 2000 platform (San Diego, USA). Compared with the gene expression profile of the control group, 572 differentially expressed genes (DEGs) were screened, with 150 at 0 psu, 170 at 30 psu, and 252 at 45 psu. Additionally, among these DEGs, Gene Ontology (GO) analysis indicated that binding, metabolic processes and cellular processes were significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis detected 3, 5 and 8 pathways related to signal transduction, metabolism, digestive and endocrine systems at 0 psu, 30 psu and 45 psu, respectively. Based on GO enrichment analysis and manual literature searches, the results of the present study indicated that A. ommaturus mainly responded to energy metabolism, ion transport and signal transduction to resist the damage caused by salinity stress. Eight DEGs were randomly selected for further validation by quantitative real-time PCR (qRT-PCR) and the results were consistent with the RNA-seq data.

Molecular characterization of p38 MAPK from blunt snout bream (Megalobrama amblycephala) and its expression after ammonia stress, and lipopolysaccharide and bacterial challenge

p38 mitogen-activated protein kinase (MAPK) is an important protein which plays a key role in regulating the innate immunity, so exploring its molecular characterization is helpful in understanding the resistance against microbial infections in cultured fish. Here, a full-length cDNA of p38 MAPK was cloned from liver of blunt snout bream (Megalobrama amblycephala) which covered 2419 bp with an open reading frame of 1086 bp encoding 361 amino acids. p38 MAPK contained the characteristic structures of Thr-Gly-Tyr (TGY) motif and substrate binding site Ala-Thr-Arg-Trp (ATRW), which are conserved in MAPK family. To investigate p38 MAPK functions, two in vivo experiments were carried out to examine its expression following ammonia exposure and bacterial challenge. Also, an in vitro experiment was conducted to assess the role of p38 MAPK in inflammation of primary hepatocytes induced by lipopolysaccharide (LPS). The results showed the ubiquitous expression of p38 MAPK in all the tested tissues with varying levels. p38 MAPK mRNA expression was significantly up-regulated by ammonia stress and Aeromonas hydrophila challenge, and altered in a time-dependent manner. Moreover, the results indicated that the inflammatory response induced by LPS in hepatocytes is p38 MAPK dependent as knockdown of p38 MAPK using siRNA technology depressed the expression of IL-1β and IL-6. The findings in this study showed that p38 MAPK has anti-stress property, and plays key role in protection against bacterial infection and inflammation in blunt snout bream.

The Change of Teleost Skin Commensal Microbiota Is Associated With Skin Mucosal Transcriptomic Responses During Parasitic Infection by Ichthyophthirius multifillis

Teleost skin serves as the first line of defense against invading pathogens, and contain a skin-associated lymphoid tissue (SALT) that elicit gut-like immune responses against antigen stimulation. Moreover, exposed to the water environment and the pathogens therein, teleost skin is also known to be colonized by diverse microbial communities. However, little is known about the interactions between microbiota and the teleost skin mucosal immune system, especially dynamic changes about the interactions under pathogen infection. We hypothesized that dramatic changes of microbial communities and strong mucosal immune response would be present in the skin of aquatic vertebrate under parasite infection. To confirm this hypothesis, we construct an infected model with rainbow trout (Oncorhynchus mykiss), which was experimentally challenged by Ichthyophthirius multifiliis (Ich). H & E staining of trout skin indicates the successful invasion of Ich and shows the morphological changes caused by Ich infection. Critically, increased mRNA expression levels of immune-related genes were detected in trout skin from experimental groups using qRT-PCR, which were further studied by RNA-Seq analysis. Here, through transcriptomics, we detected that complement factors, pro-inflammatory cytokines, and antimicrobial genes were strikingly induced in the skin of infected fish. Moreover, high alpha diversity values of microbiota in trout skin from the experimental groups were discovered. Interestingly, we found that Ich infection led to a decreased abundance of skin commensals and increased colonization of opportunistic bacteria through 16S rRNA pyrosequencing, which were mainly characterized by lose of Proteobacteria and increased intensity of Flavobacteriaceae. To our knowledge, our results suggest for the first time that parasitic infection could inhibit symbionts and offer opportunities for other pathogens' secondary infection in teleost skin.

Transcriptome profiling and histology changes in juvenile blunt snout bream (Megalobrama amblycephala) liver tissue in response to acute thermal stress

To understand the precise mechanism and the pathways activated by thermal stress in fish, we sampled livers from juvenile Megalobrama amblycephala exposed to control (25 °C) and test (35 °C) conditions, and performed short read (100 bp) next-generation RNA sequencing (RNA-seq). Using reads from different temperature, expression analysis identified a total of 440 differentially-expressed genes. These genes were related to oxidative stress, apoptosis, immune responses and so on. We used quantitative real-time reverse transcriptase PCR to assess the differential mRNA expression of selected genes that encode antioxidant enzymes and heat shock proteins in response to thermal stress. Fish exposed to thermal stress also showed liver damage associated with serum biochemical parameter changes. The set of genes identified showed regulatory modulation at different temperatures, and therefore could be further studied to determine how thermal stress damages M. amblycephala livers and the possible roles of reactive oxygen species in this process.

MiR-5702 suppresses proliferation and invasion in non-small-cell lung cancer cells via posttranscriptional suppression of ZEB1

MiRNAs have emerged as important players in tumorigenesis and progression. MiR-5702 is a newly identified miRNA; the exact role of which has not been reported. Here, we found that miR-5702 was significantly decreased in the carcinoma tissues of non-small cell lung cancer (NSCLC) patients and NSCLC cell lines. Then, our results showed that the miR-5702 mimic induced apoptosis and inhibited proliferation and invasion in A549 cells. In contrast, the miR-5702 inhibitor reduced apoptosis and increased proliferation and invasion in A549 cells. Furthermore, bioinformatics and 3'-UTR luciferase reporter assays identified that oncogene zinc finger E-box-binding homeobox 1 (ZEB1) is a target gene of miR-5702. Western blotting analysis showed that miR-5702 overexpression suppressed, and miR-5702 knockdown promoted the expression of ZEB1 protein. Finally, the ZEB1 siRNA exhibited a similar effect to the miR-5702 mimic on expression of ZEB1 and its downstream genes, cell apoptosis, cell proliferation, and cell invasion, and it could antagonize the alternations in ZEB1 expression and cell behaviors by the miR-5702 inhibitor. In conclusion, miR-5702 may function as a tumor suppressor in NSCLC, which suppresses proliferation and invasion NSCLC cells via posttranscriptional suppression of ZEB1.

Suppression of miR-26a attenuates physiological disturbances arising from exposure of Nile tilapia (Oreochromis niloticus) to ammonia

MicroRNAs may affect stress responses because they act as rapid responders at the post-translation level. In this study, we found that miR-26a is abundantly expressed in the brain and gill tissues of tilapia. Expression of miR-26a in the brain decreased significantly with increasing ammonia concentrations using stem-loop qPCR. To analyze the function of miRNA in vivo, miR-26a was stably knocked down with an antagomir in tilapia. Following ammonia challenge, miR-26a antagomir treatment significantly suppressed blood ammonia/[Cl⁻]/[K⁺] concentration and the reactive oxygen species production, while it markedly enhanced glutamine accumulation and antioxidant enzyme activity in the brain of tilapia, indicating that miR-26a may be involved in the remission of physiological disturbances resulting from ammonia stress. We strongly conclude that there is a direct link between miR-26a and the responses to ammonia in tilapia. Furthermore, bioinformatics analysis and luciferase assays demonstrated that miR-26a regulates HSP70 (heat shock protein 70) and GS (glutamine synthetase) expression by targeting their 3′-UTR and that the suppression of miR-26a could increase the intracellular level of HSP70 and GS in vivo.

Summary: Our work increases the available information about the regulation of miR-26a and indicates that miR-26a may be involved in the remission of physiological disorders upon ammonia stress in tilapia.

Discovery and characterization of conserved and novel microRNAs from blunt snout bream (Megalobrama amblycephala) by deep sequencing

MicroRNAs (miRNAs) are short, single stranded RNA molecules with approximately 22 nts in length, which regulate the stability and translation of messenger RNAs in several organisms. To increase the repertoire of miRNAs characterized in M. amblycephala, we used the deep sequencing technology to sequence a small RNA library using pooled RNA sample isolated from the 4 different tissues of M. amblycephala. A total of 309 conserved miRNAs that originated from 131 miRNA families were detected. 15 novel candidates miRNA were identified. Randomly selected 6 miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in 6 different tissues of M. amblycephala. Furthermore, the potential targets were predicted. GO analysis showed that most of the targets were involved in a broad range of physiological functions including fish growth, development, metabolism, stress responses and so on. Overall, our results significantly increased the number of novel miRNAs in M. amblycephala, which should be useful for further investigation into the role of miRNAs in regulating diverse biological processes.

Dynamic mRNA and miRNA expression analysis in response to hypoxia and reoxygenation in the blunt snout bream (Megalobrama amblycephala)

Adaptation to hypoxia is a complex process involving various pathways and regulation mechanisms. A better understanding of the genetic influence on these mechanisms could permit selection for hypoxia-sensitive fish. To aid this understanding, an integrated analysis of miRNA and mRNA expression was performed in Megalobrama amblycephala under four acute hypoxia and reoxygenation stages. A number of significantly differentially-expressed miRNAs and genes associated with oxidative stress were identified, and their functional characteristics were revealed by GO function and KEGG pathway analysis. They were found to be involved in HIF-1 pathways known to affect energy metabolism and apoptosis. MiRNA-mRNA interaction pairs were detected from comparison of expression between the four different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating hypoxia stress. As a unique resource for gene expression and regulation during hypoxia and reoxygenation, this study could provide a starting point for further studies to better understand the genetic background of hypoxia stress.

Genomic and transcriptomic approaches to study immunology in cyprinids: What is next?

Accelerated by the introduction of Next-Generation Sequencing (NGS), a number of genomes of cyprinid fish species have been drafted, leading to a highly valuable collective resource of comparative genome information on cyprinids (Cyprinidae). In addition, NGS-based transcriptome analyses of different developmental stages, organs, or cell types, increasingly contribute to the understanding of complex physiological processes, including immune responses. Cyprinids are a highly interesting family because they comprise one of the most-diversified families of teleosts and because of their variation in ploidy level, with diploid, triploid, tetraploid, hexaploid and sometimes even octoploid species. The wealth of data obtained from NGS technologies provides both challenges and opportunities for immunological research, which will be discussed here. Correct interpretation of ploidy effects on immune responses requires knowledge of the degree of functional divergence between duplicated genes, which can differ even between closely-related cyprinid fish species. We summarize NGS-based progress in analysing immune responses and discuss the importance of respecting the presence of (multiple) duplicated gene sequences when performing transcriptome analyses for detailed understanding of complex physiological processes. Progressively, advances in NGS technology are providing workable methods to further elucidate the implications of gene duplication events and functional divergence of duplicates genes and proteins involved in immune responses in cyprinids. We conclude with discussing how future applications of NGS technologies and analysis methods could enhance immunological research and understanding.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.