The integrated analysis of RNA-seq and microRNA-seq depicts miRNA-mRNA networks involved in Japanese flounder (Paralichthys olivaceus) albinism

Transcriptomic analyses of differentially expressed human genes, micro RNAs and long-non-coding RNAs in severe, symptomatic and asymptomatic malaria infection

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

microRNA regulation of skin pigmentation in golden-back mutant of crucian carp from a rice-fish integrated farming system

BACKGROUND

MicroRNAs (miRNAs) are endogenous small non-coding RNAs (21-25 nucleotides) that act as essential components of several biological processes. Golden-back crucian carp (GBCrC, Carassius auratus) is a naturally mutant species of carp that has two distinct body skin color types (golden and greenish-grey), making it an excellent model for research on the genetic basis of pigmentation. Here, we performed small RNA (sRNA) analysis on the two different skin colors via Illumina sequencing.

RESULTS

A total of 679 known miRNAs and 254 novel miRNAs were identified, of which 32 were detected as miRNAs with significant differential expression (DEMs). 23,577 genes were projected to be the targets of 32 DEMs, primarily those involved in melanogenesis, adrenergic signaling in cardiomyocytes, MAPK signaling pathway and wnt signaling pathway by functional enrichment. Furthermore, we built an interaction module of mRNAs, proteins and miRNAs based on 10 up-regulated and 13 down-regulated miRNAs in golden skin. In addition to transcriptional destabilization and translational suppression, we discovered that miRNAs and their target genes were expressed in the same trend at both the transcriptional and translational levels. Finally, we discovered that miR-196d could be indirectly implicated in regulating melanocyte synthesis and motility in the skin by targeting to myh7 (myosin-7) gene through the luciferase reporter assay, antagomir silencing in vivo and qRT-PCR techniques.

CONCLUSIONS

Our study gives a systematic examination of the miRNA profiles expressed in the skin of GBCrC, assisting in the comprehension of the intricate molecular regulation of body color polymorphism and providing insights for C. auratus breeding research.

Comparative transcriptome and miRNA analysis of skin pigmentation during embryonic development of Chinese soft-shelled turtle (Pelodiscus sinensis)

BACKGROUND

Aquatic animals show diverse body coloration, and the formation of animal body colour is a complicated process. Increasing evidence has shown that microRNAs (miRNAs) play important regulatory roles in many life processes. The role of miRNAs in pigmentation has been investigated in some species. However, the regulatory patterns of miRNAs in reptile pigmentation remain to be elucidated. In this study, we performed an integrated analysis of miRNA and mRNA expression profiles to explore corresponding regulatory patterns in embryonic body colour formation in the soft-shelled turtle Pelodiscus sinensis.

RESULTS

We identified 8 866 novel genes and 9 061 mature miRNAs in the skin of Chinese soft-shelled turtles in three embryonic stages (initial period: IP, middle period: MP, final period: FP). A total of 16 563 target genes of the miRNAs were identified. Furthermore, we identified 2 867, 1 840 and 4 290 different expression genes (DEGs) and 227, 158 and 678 different expression miRNAs (DEMs) in IP vs. MP, MP vs. FP, and IP vs. FP, respectively. Among which 72 genes and 25 miRNAs may be related to turtle pigmentation in embryonic development. Further analysis of the novel miRNA families revealed that some novel miRNAs related to pigmentation belong to the miR-7386, miR-138, miR-19 and miR-129 families. Novel_miR_2622 and novel_miR_2173 belong to the miR-19 family and target Kit and Gpnmb, respectively. The quantification of novel_miR_2622 and Kit revealed negative regulation, indicating that novel_miR_2622 may participate in embryonic pigmentation in P. sinensis by negatively regulating the expression of Kit.

CONCLUSIONS

miRNA act as master regulators of biological processes by controlling the expression of mRNAs. Considering their importance, the identified miRNAs and their target genes in Chinese soft-shelled turtle might be useful for investigating the molecular processes involved in pigmentation. All the results of this study may aid in the improvement of P. sinensis breeding traits for aquaculture.

Integrative analyses of mRNA and microRNA expression profiles reveal the innate immune mechanism for the resistance to Vibrio parahaemolyticus infection in Epinephelus coioides

Vibrio parahaemolyticus, as one of the main pathogens of marine vibriosis, has brought huge losses to aquaculture. However, the interaction mechanism between V. parahaemolyticus and Epinephelus coioides remains unclear. Moreover, there is a lack of comprehensive multi-omics analysis of the immune response of grouper spleen to V. parahaemolyticus. Herein, E. coioides was artificially injected with V. parahaemolyticus, and it was found that the mortality was 16.7% in the early stage of infection, and accompanied by obvious histopathological lesions in the spleen. Furthermore, 1586 differentially expressed genes were screened by mRNA-seq. KEGG analysis showed that genes were significantly enriched in immune-related pathways, Acute-phase immune response, Apoptosis, Complement system and Cytokine-cytokine receptor interaction. As for miRNA-seq analysis, a total of 55 significantly different miRNAs were identified. Further functional annotation analysis indicated that the target genes of differentially expressed miRNAs were enriched in three important pathways (Phosphatidylinositol signaling system, Lysosome and Focal adhesions). Through mRNA-miRNA integrated analysis, 1427 significant miRNA–mRNA pairs were obtained and “p53 signaling pathway”, “Intestinal immune network for IgA production” were considered as two crucial pathways. Finally, miR-144-y, miR-497-x, novel-m0459-5p, miR-7133-y, miR-378-y, novel-m0440-5p and novel-m0084-3p may be as key miRNAs to regulate immune signaling pathways via the miRNA-mRNA interaction network. The above results suggest that the mRNA-miRNA integrated analysis not only sheds new light on the molecular mechanisms underlying the interaction between host and V. parahaemolyticus but also provides valuable and new insights into resistance to vibrio infection.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Identification and Characterization of Long Non-coding RNAs in the Intestine of Olive Flounder (Paralichthys olivaceus) During Edwardsiella tarda Infection

Long non-coding RNAs (lncRNAs) play widespread roles in fundamental biological processes, including immune responses. The olive flounder (Paralichthys olivaceus), an important economical flatfish widely cultured in Japan, Korea, and China, is threatened by infectious pathogens, including bacteria, viruses, and parasites. However, the role of lncRNAs in the immune responses of this species against pathogen infections is not well-understood. Therefore, in this study, we aimed to identify lncRNAs in the intestine of olive flounder and evaluate their differential expression profiles during Edwardsiella tarda infection, which is an important zoonotic and intestinal pathogen. A total of 4,445 putative lncRNAs were identified, including 3,975 novel lncRNAs and 470 annotated lncRNAs. These lncRNAs had shorter lengths and fewer exons compared with mRNAs. In total, 115 differentially expressed lncRNAs (DE-lncRNAs) were identified during E. tarda infection. To validate the expression pattern of lncRNAs, six DE-lncRNAs were randomly selected for quantitative real-time PCR. The co-located and co-expressed mRNAs of DE-lncRNAs were predicted, which were used to conduct the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The target genes of DE-lncRNAs enriched numerous immune-related processes and exhibited a strong correlation with immune-related signaling pathways. To better understand the extensive regulatory functions of lncRNAs, the lncRNA–miRNA–mRNA regulatory networks were constructed, and two potential competing endogenous RNA (ceRNA) networks, LNC_001979-novel_171-Potusc2 and LNC_001979-novel_171-Podad1, were preliminarily identified from the intestine of olive flounders for the first time. In conclusion, this study provides an invaluable annotation and expression profile of lncRNAs in the intestine of olive flounder infected with E. tarda; this forms a basis for further studies on the regulatory function of lncRNAs in the intestinal mucosal immune responses of olive flounder.

Insights Into the Albinism Mechanism for Two Distinct Color Morphs of Northern Snakehead, Channa argus Through Histological and Transcriptome Analyses

The great northern snakehead (Channa argus) is one of the most important economic and conservational fish in China. In this study, the melanocytes in the skin of two distinct color morphs C. argus were investigated and compared through employment of the microscopic analysis, hematoxylin and eosin (H&E) and Masson Fontana staining. Our results demonstrated the uneven distribution of melanocytes with extremely low density and most of them were in the state of aging or death. Meanwhile, there was no obvious pigment layer and melanocytes distribution pattern found in the albino-type (AT), while the melanocytes were evenly distributed with abundance in the bicolor-type (BT). The transcriptome analysis through Illumina HiSeq sequencing showed that a total of 34.93 Gb Clean Data was obtained, and Q30 base percentage reached 92.66%. The BT and AT northern snakeheads transcriptome data included a total of 56,039,701 and 60,410,063 clean reads (n = 3), respectively. In gene expression analyses, the sample correlation coefficients (r) were ranged between 0.92 and 1.00; the contribution of PC1 and PC2 were 50.25 and 13.73% by using PCA cluster analysis, the total number of DEGs were 1024 (559 up-regulated and 465 down-regulated), and the number of annotated DEGs was 767 (COG 172, KEGG 262, GO 288, SwissProt 548, Pfam 579 and NR 765). Additionally, 46,363 ± 873 and 44,947 ± 392 single nucleotide polymorphisms (SNPs) were compiled via genetic structure analysis, respectively. Ten key pigment-related genes were screened using qRT-PCR. And all of them revealed extremely higher expression levels in the skin of BT than those of AT. This is the first study to analyze the mechanism of albino characteristics of Channa via histology and transcriptomics, and also provide the oretical and practical support for the protection and development of germplasm resources for C. argus.

RNA-seq and microRNA-seq analysis of Japanese flounder (Paralichthys olivaceus) larvae treated by thyroid hormones

Flatfish pigmentation is a complex process, affected by environmental factors including light, nutrients, and hormones. Of those, the thyroid hormone has been reported to increase the albinism rate of Japanese flounder (Paralichthys olivaceus). However, the underlying mechanism remains unclear. In the present study, triiodothyronine (T3), thyroxine, and thiourea were introduced into P. olivaceus larvae from 16 to 57 days after hatching (DAH). By comparison of albinism rate, T3 treatment and control larvae of 42 DAH were chosen for mRNA and miRNA high-throughput sequencing analyses. A total of 337 miRNAs were identified via miRNA-seq, and 12 miRNAs exhibited significantly differential expression patterns in D42_T3 versus D42_Con (TPM > 10, fold change ≥ 1.5 or ≤ 0.67 and q ≤ 0.05). These differentially expressed miRNAs targeted 3658 putative genes, which further enriched to 10 GO terms (q < 0.05). RNA-seq identified 146 differentially expressed genes (DEGs) in D42_T3 versus D42_Con (|log₂ fold change| > 1 and q < 0.005), including pigmentation-related genes such as the receptor tyrosine-protein kinase erbB-3, pro-opiomelanocortin A, and melanotransferrin, and the growth-related gene somatotropin. These DEGs were significantly enriched to 15 GO terms and 8 KEGG pathways (q < 0.05), which included several sugar metabolic pathways (glycolysis/gluconeogenesis and the pentose phosphate pathway). Integrated analysis revealed that 26 overlapping genes between DEGs and mRNAs were targeted by miRNAs. Furthermore, seven mRNA-miRNA pairs exhibited reversed regulation patterns. This provides important clues to understand the role of thyroid hormones in flatfish pigmentation.

Integrated Analysis of Transcriptomic, miRNA and Proteomic Changes of a Novel Hybrid Yellow Catfish Uncovers Key Roles for miRNAs in Heterosis

Heterosis is a complex biological phenomenon in which hybridization produces offspring that exhibit superior phenotypic characteristics compared with the parents. Heterosis is widely utilized in agriculture, for example in fish farming; however, its underlying molecular basis remains elusive. To gain a comprehensive and unbiased molecular understanding of fish heterosis, we analyzed the mRNA, miRNA, and proteomes of the livers of three catfish species, Pelteobagrus fulvidraco, P. vachelli, and their hybrid, the hybrid yellow catfish "Huangyou-1" (P. fulvidraco ♀ × P. vachelli ♂). Using next-generation sequencing and mass spectrometry, we show that the nonadditive, homoeolog expression bias and expression level dominance pattern were readily identified at the transcriptional, post-transcriptional, or protein levels, providing the evidence for the widespread presence of dominant models during hybridization. A number of predicted miRNA-mRNA-protein pairs were found and validated by qRT-PCR and PRM assays. Furthermore, several diverse key pathways were identified, including immune defense, metabolism, digestion and absorption, and cell proliferation and development, suggesting the vital mechanisms involved in the generation of the heterosis phenotype in progenies. We propose that the high parental expression of genes/proteins (growth, nutrition, feeding, and disease resistance) coupled with low parental miRNAs of the offspring, are inherited from the mother or father, thus indicating that the offspring were enriched with the advantages of the father or mother. We provide new and important information about the molecular mechanisms of heterosis, which represents a significant step toward a more complete elucidation of this phenomenon.

Comparative analysis of distinctive transcriptome profiles with biochemical evidence in bisphenol S- and benzo[a]pyrene-exposed liver tissues of the olive flounder Paralichthys olivaceus

Flounder is a promising model species for environmental monitoring of coastal regions. To assess the usefulness of liver transcriptome profiling, juvenile olive flounder Paralichthys olivaceus were exposed to two pollutants, bisphenol S (BPS) and benzo[a]pyrene (BaP), which have different chemical characteristics and have distinct modes of metabolic action in teleost. Six hours after intraperitoneal injection with BPS (50 mg/kg bw) or BaP (20 mg/kg bw), liver transcriptomes were analyzed using the Illumina Hiseq 3000 platform. Interestingly, the transcriptome was highly sensitive and was distinctively expressed in response to each chemical. The primary effect of BPS was significantly increased transcription of egg process and vitellogenesis related genes, including vitellogenins (vtg1, vtg2), zona pellucida sperm-binding proteins (zp3, zp4), and estrogen receptors (erα, erβ), with increases in plasma 17β-estradiol (E2) and vitellogenin (VTG) concentrations. Following BaP treatment, detoxification- and biotransformation-related genes such as cyp1a1 and UDP-glucuronosyltransferase (ugt1a1) were significantly increased, with an increase in EROD activity. In both transcriptomes, mRNA expression of genes involved in antioxidant defense systems was increased, while genes involved in innate immunity were decreased upon BPS or BaP exposure with a decrease in complement activity. This study provides useful insight into the chemical-specific hepatic transcriptional response of P. olivaceus and suggests a basis for further studies examining biomarker application of liver transcriptomes for environmental pollution.