Immune-Related Transcriptional Responses to Parasitic Infection in a Naturally Inbred Fish: Roles of Genotype and Individual Variation

Integrated transcriptomics and lipidomics investigation of the mechanism underlying the gastrointestinal mucosa damage of Loropetalum chinense (R.Br.) and its representative component

BACKGROUND

Loropetalum chinensis (R.Br) Oliv (Bhjm), a Chinese folk herbal medicine, was traditionally used in the treatment of wound bleeding and skin ulcers. A new drug named JIMUSAN granules used for gastrosia was developed by our group, and clinical trials have been approved. However, as the principal herb, the material basis and underlying mechanisms of Bhjm in attenuating gastrointestinal mucosa damage (GMD) remain to be systemically illuminated.

PURPOSE

An integrated strategy was used to explore the therapeutic effects and mechanisms of Bhjm and ellagic acid (EA) on GMD zebrafish, using network pharmacology, transcriptomics, lipidomics, and real-time quantitative PCR (RT-qPCR) verification.

METHODS

First, network pharmacological analysis was used to infer the major effective constituents and targets of Bhjm. Ultra high performance liquid chromatography-linear ion trap/orbitrap high resolution mass spectrometry (UHPLC-LTQ-Orbitrap HRMS) and ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) were employed to identify the chemical constituents and quantify the different types of constituents. Second, zebrafish model of GMD was established by using 2,4,6-trinitrobenzenesulfonic acid (TNBS) to evaluate the efficacy of Bhjm and EA. The potential mechanism was examined by integrated transcriptomics and lipidomics analysis. Finally, validation tests were implemented using RT-qPCR.

RESULTS

In this study, targets indentified by network pharmacology were related to inflammation and mucosal damage. Ten representative components that interacted with these targets were simultaneously determined by UHPLC-MS/MS. Sixty four compounds were identified or tentatively characterized, most of which were flavonoids and polyphenols. Bhjm and EA alleviated mucosal damage and reduced inflammation in a TNBS-induced zebrafish GMD model, indicating that EA was the main active compounds. Eight common differentially expressed genes were downregulated by Bhjm and EA, as determined by transcriptomics analysis. Lipidomics analysis confirmed 12 differential lipids, including phosphatidylcholine (PC) and triglyceride (TG). Further network enrichment analysis demonstrated that differential lipid metabolism was regulated by klf4 and hist1h2ba, and were validated by RT-qPCR.

CONCLUSION

In our study, the chemical profile of Bhjm was clarified. Moreover, the GMD repair effect and the mechanism of Bhjm and EA was comprehensively analyzed for the first time, involving inflammation and lipid metabolism. Collectively, these findings will be significantly helpful for deeply exploring the clinical application value of Bhjm.

Diversity and structure of sparids external microbiota (Teleostei) and its link with monogenean ectoparasites

Background

Animal-associated microbial communities appear to be key factors in host physiology, ecology, evolution and its interactions with the surrounding environment. Teleost fish have received relatively little attention in the study of surface-associated microbiota. Besides the important role of microbiota in homeostasis and infection prevention, a few recent studies have shown that fish mucus microbiota may interact with and attract some specific parasitic species. However, our understanding of external microbial assemblages, in particular regarding the factors that determine their composition and potential interactions with parasites, is still limited. This is the objective of the present study that focuses on a well-known fish-parasite interaction, involving the Sparidae (Teleostei), and their specific monogenean ectoparasites of the Lamellodiscus genus. We characterized the skin and gill mucus bacterial communities using a 16S rRNA amplicon sequencing, tested how fish ecological traits and host evolutionary history are related to external microbiota, and assessed if some microbial taxa are related to some Lamellodiscus species.

Results

Our results revealed significant differences between skin and gill microbiota in terms of diversity and structure, and that sparids establish and maintain tissue and species-specific bacterial communities despite continuous exposure to water. No phylosymbiosis pattern was detected for either gill or skin microbiota, suggesting that other host-related and environmental factors are a better regulator of host-microbiota interactions. Diversity and structure of external microbiota were explained by host traits: host species, diet and body part. Numerous correlations between the abundance of given bacterial genera and the abundance of given Lamellodiscus species have been found in gill mucus, including species-specific associations. We also found that the external microbiota of the only unparasitized sparid species in this study, Boops boops, harbored significantly more Fusobacteria and three genera, Shewenella, Cetobacterium and Vibrio, compared to the other sparid species, suggesting their potential involvement in preventing monogenean infection.

Conclusions

This study is the first to explore the diversity and structure of skin and gill microbiota from a wild fish family and present novel evidence on the links between gill microbiota and monogenean species in diversity and abundance, paving the way for further studies on understanding host-microbiota-parasite interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00180-1.

Molecular assessment and transcriptome profiling of wild fish populations of Oryzias mekongensis and O. songkhramensis (Adrianichthyidae: Beloniformes) from Thailand

Among the fish of the genus Oryzias, two species are frequently used as model animals in biological research. In Thailand, Oryzias mekongensis is usually found in natural freshwater near the Mekong Basin in the northeast region, while O. songkhramensis inhabits the Songkhram Basin. For differential morphological identification, the coloured bands on the dorsal and ventral margins of the caudal fin are used to distinguish O. mekongensis from O. songkhramensis. However, these characteristics are insufficient to justify species differentiation, and little molecular evidence is available to supplement them. This study aimed to investigate the molecular population and transcriptome profiles of adult O. mekongensis and O. songkhramensis. In the molecular tree based on cytochrome b sequences, O. mekongensis exhibited four clades that were clearly distinguished from O. songkhramensis. Clade 1 of the O. mekongensis population was close to the Mekong River and lived in the eastern portion of the upper northeast region. Clade 2 was far from the Mekong River and inhabited the middle region of the Songkhram River. Clade 3 was positioned to the west of the Songkhram River, and clade 4 was to the south of the Songkhram River Basin. After RNA sequencing using an Illumina HiSeq 2500 platform, the gene category annotations hardly differentiated the species and were discussed in the text. Based on the present findings, population dispersal of these Oryzias species might be associated with geographic variations of the upper northeast region. Molecular genetics and transcriptome profiling might advance our understanding of the evolution of teleost fish.

Molecular analysis of population and De Novo transcriptome sequencing of Thai medaka, Oryzias minutillus (Teleostei: Adrianichthyidae)

Thai medaka (Oryzias minutillus) are alternatively known as Thai rice-fish or dwarf medaka, and they widely inhabit natural freshwater environments in all regions of Thailand. In this study, we aimed to investigate the molecular genetics of the Thai medaka population in Thailand inferred from the mitochondrial control region (D-loop) and the cytochrome c oxidase subunit 1 (coxI) sequences. Furthermore, we examined RNA sequencing (RNA-seq) of adult males and females was performed with next-generation sequencing. Together, the combination of the D-loop and coxI sequences clearly distinguished the Thai medaka populations into 2 groups, such as group 1, which generally included samples from the central, northern, western, and eastern regions of the northeastern region. In this group, the fish populations seem to be a little monophyly in which the first subpopulation comprised the main samples from the northern and central regions. The second subpopulation commonly contained fish from the eastern region and specimens from the southern part of the central region near the Gulf of Thailand. Although these subgroups related to geographical distribution, bootstrap values were low in branch considered significant for both subgroups. Group 2 consisted of almost all samples from the southern population and those from the central and southern part of the northeastern region. Group 2 was found that it was made of samples from the northeastern region and samples from the southern population. A total of 73551 unigenes were identified after gene annotation. Signal transduction was the predominant protein classification among the Thai medaka orthologous groups. A differentially expressed gene (DEG) analysis identified 6 subclusters between both sexes that were composed of 257, 131, 364, 386, 114 and 108 genes. Phototransduction was the most enriched pathway and was highly expressed in males, while viral carcinogenesis, oocyte genesis, and the complement and coagulation cascades were highly expressed in females. Further details of these DEGs are discussed below. These results suggest that Thai medaka may genetically exhibit independent populations in the geographic habitats of Thailand. Moreover, these fish also reveal the genes that are conserved in other organisms and those that may be specific to this species.

Transcriptomic analysis of Baltic cod (Gadus morhua) liver infected with Contracaecum osculatum third stage larvae indicates parasitic effects on growth and immune response

High infection levels due to third-stage larvae of the anisakid nematode Contracaecum osculatum have been documented in cod from the eastern part of the Baltic sea during the latest decades. The nematode larvae mainly infect the liver of Baltic cod and prevalence of infection has reached 100% with a mean intensity up to 80 parasites per host in certain areas and size classes. Low condition factors of the cod have been observed concomitant with the rise in parasite abundance suggesting a parasitic effect on growth parameters. To investigate any association between parasite infection and physiological status of the host we performed a comparative transcriptomic analysis of liver obtained from C. osculatum infected and non-infected cod. A total of 47,025 predicted gene models showed expression in cod liver and sequences corresponding to 2084 (4.43%) unigenes were differentially expressed in infected liver when compared to non-infected liver. Of the differentially expressed unigenes (DEGs) 1240 unigenes were up-regulated while 844 unigenes were down-regulated. The Gene Ontology (GO) enrichment analysis showed that 1304 DEGs were represented in cellular process and single-organism process, cell and cell part, binding and catalytic activity. As determined by the Kyoto Encyclopedia of Gene and Genomes (KEGG) Pathways analysis, 454 DEGs were involved in 138 pathways. Ninety-seven genes were related to metabolic pathways including carbohydrate, lipid, and amino acid metabolism. Thirteen regulated genes were playing a role in immune response such as Toll-like receptor signaling, NOD-like receptor signaling, RIG-I-like receptor signalling and thirty-six genes were associated with growth processes. This indicates that the nematode infection in Baltic cod may affect on molecular mechanisms involving metabolism, immune function and growth.

Colour plasticity in response to social context and parasitic infection in a self-fertilizing fish

Many animal species rely on changes in body coloration to signal social dominance, mating readiness and health status to conspecifics, which can in turn influence reproductive success, social dynamics and pathogen avoidance in natural populations. Such colour changes are thought to be controlled by genetic and environmental conditions, but their relative importance is difficult to measure in natural populations, where individual genetic variability complicates data interpretation. Here, we studied shifts in melanin-related body coloration in response to social context and parasitic infection in two naturally inbred lines of a self-fertilizing fish to disentangle the relative roles of genetic background and individual variation. We found that social context and parasitic infection had a significant effect on body coloration that varied between genetic lines, suggesting the existence of genotype by environment interactions. In addition, individual variation was also important for some of the colour attributes. We suggest that the genetic background drives colour plasticity and that this can maintain phenotypic variation in inbred lines, an adaptive mechanism that may be particularly important when genetic diversity is low.

Contrasting DNA methylation responses of inbred fish lines to different rearing environments

Epigenetic mechanisms generate plastic phenotypes that can become locally adapted across environments. Disentangling genomic from epigenomic variation is challenging in sexual species due to genetic variation among individuals, but it is easier in self-fertilizing species. We analysed DNA methylation patterns of two highly inbred strains of a naturally self-fertilizing fish reared in two contrasting environments to investigate the obligatory (genotype-dependent), facilitated (partially depend on the genotype) or pure (genotype-independent) nature of the epigenetic variation. We found higher methylation differentiation between genotypes than between environments. Most methylation differences between environments common to both strains followed a pattern where the two genotypes (inbred lines) responded to the same environmental context with contrasting DNA methylation levels (facilitated epialleles). Our findings suggest that, at least in part, DNA methylation could depend on the dynamic interaction between the genotype and the environment, which could explain the plasticity of epigenetically mediated phenotypes.

Acting locally - affecting globally: RNA sequencing of gilthead sea bream with a mild Sparicotyle chrysophrii infection reveals effects on apoptosis, immune and hypoxia related genes

BACKGROUND

Monogenean flatworms are the main fish ectoparasites inflicting serious economic losses in aquaculture. The polyopisthocotylean Sparicotyle chrysophrii parasitizes the gills of gilthead sea bream (GSB, Sparus aurata) causing anaemia, lamellae fusion and sloughing of epithelial cells, with the consequent hypoxia, emaciation, lethargy and mortality. Currently no preventive or curative measures against this disease exist and therefore information on the host-parasite interaction is crucial to find mitigation solutions for sparicotylosis. The knowledge about gene regulation in monogenean-host models mostly comes from freshwater monopysthocotyleans and almost nothing is known about polyopisthocotyleans. The current study aims to decipher the host response at local (gills) and systemic (spleen, liver) levels in farmed GSB with a mild natural S. chrysophrii infection by transcriptomic analysis.

RESULTS

Using Illumina RNA sequencing and transcriptomic analysis, a total of 2581 differentially expressed transcripts were identified in infected fish when compared to uninfected controls. Gill tissues in contact with the parasite (P gills) displayed regulation of fewer genes (700) than gill portions not in contact with the parasite (NP gills) (1235), most likely due to a local silencing effect of the parasite. The systemic reaction in the spleen was much higher than that at the parasite attachment site (local) (1240), and higher than in liver (334). NP gills displayed a strong enrichment of genes mainly related to immune response and apoptosis. Processes such as apoptosis, inflammation and cell proliferation dominated gills, whereas inhibition of apoptosis, autophagy, platelet activation, signalling and aggregation, and inflammasome were observed in spleen. Proteasome markers were increased in all tissues, whereas hypoxia-related genes were down-regulated in gills and spleen.

CONCLUSIONS

Contrasting forces seem to be acting at local and systemic levels. The splenic down-regulation could be part of a hypometabolic response, to counteract the hypoxia induced by the parasite damage to the gills and to concentrate the energy on defence and repair responses. Alternatively, it can be also interpreted as the often observed action of helminths to modify host immunity in its own interest. These results provide the first toolkit for future studies towards understanding and management of this parasitosis.

Smell of Infection: A Novel, Noninvasive Method for Detection of Fish Excretory-Secretory Proteins

Chemical signals are produced by aquatic organisms following predatory attacks or perturbations such as parasitic infection. Ectoparasites feeding on fish hosts are likely to cause release of similar alarm cues into the environment due to the stress, wounding, and immune response stimulated upon infection. Alarm cues are often released in the form of proteins, antimicrobial peptides, and immunoglobulins that provide important insights into bodily function and infection status. Here we outline a noninvasive method to identify potential chemical cues associated with infection in fish by extracting, purifying, and characterizing proteins from water samples from cultured fish. Gel free proteomic methods were deemed the most suitable for protein detection in saline water samples. It was confirmed that teleost proteins can be characterized from water and that variation in protein profiles could be detected between infected and uninfected individuals and fish and parasite only water samples. Our novel assay provides a noninvasive method for assessing the health condition of both wild and farmed aquatic organisms. Similar to environmental DNA monitoring methods, these proteomic techniques could provide an important tool in applied ecology and aquatic biology.