Search

Plant-specialized metabolism represents an inexhaustible source of active molecules, some of which have been used in human health for decades. Among these, monoterpene indole alkaloids (MIAs) include a wide range of valuable compounds with anticancer, antihypertensive, or neuroactive properties. This is particularly the case for the pachysiphine derivatives which show interesting antitumor and anti-Alzheimer activities but accumulate at very low levels in several Tabernaemontana species. Unfortunately, genome data in Tabernaemontanaceae are lacking and knowledge on the biogenesis of pachysiphine-related MIAs in planta remains scarce, limiting the prospects for the biotechnological supply of many pachysiphine-derived biopharmaceuticals. Here, we report a raw version of the toad tree (Tabernaemontana elegans) genome sequence. These new genomic resources led to the identification and characterization of a couple of genes encoding cytochrome P450 with pachysiphine synthase activity. Our phylogenomic and docking analyses highlight the different evolutionary processes that have been recruited to epoxidize the pachysiphine precursor tabersonine at a specific position and in a dedicated orientation, thus enriching our understanding of the diversification and speciation of the MIA metabolism in plants. These gene discoveries also allowed us to engineer the synthesis of MIAs in yeast through the combinatorial association of metabolic enzymes resulting in the tailor-made synthesis of non-natural MIAs. Overall, this work represents a step forward for the future supply of pachysiphine-derived drugs by microbial cell factories., (© 2024 The Author(s). The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Monoterpene indole alkaloids (MIAs) from Mitragyna speciosa ("kratom"), such as mitragynine and speciogynine, are promising novel scaffolds for opioid receptor ligands for treatment of pain, addiction, and depression. While kratom leaves have been used for centuries in South-East Asia as stimulant and pain management substance, the biosynthetic pathway of these psychoactives have only recently been partially elucidated. Here, we demonstrate the de novo production of mitragynine and speciogynine in Saccharomyces cerevisiae through the reconstruction of a five-step synthetic pathway from common MIA precursor strictosidine comprising fungal tryptamine 4-monooxygenase to bypass an unknown kratom hydroxylase. Upon optimizing cultivation conditions, a titer of ∼290 μg/L kratom MIAs from glucose was achieved. Untargeted metabolomics analysis of lead production strains led to the identification of numerous shunt products derived from the activity of strictosidine synthase (STR) and dihydrocorynantheine synthase (DCS), highlighting them as candidates for enzyme engineering to further improve kratom MIAs production in yeast. Finally, by feeding fluorinated tryptamine and expressing a human tailoring enzyme, we further demonstrate production of fluorinated and hydroxylated mitragynine derivatives with potential applications in drug discovery campaigns. Altogether, this study introduces a yeast cell factory platform for the biomanufacturing of complex natural and new-to-nature kratom MIAs derivatives with therapeutic potential., Competing Interests: Declaration of Competing Interest L.G.H., J.Z., J.D.K. and M.K.J. have financial interests in Biomia. J.D.K. also has financial interests in Amyris, Lygos, Demetrix, Napigen, Apertor Pharmaceuticals, Maple Bio, Ansa Biotechnologies, Berkeley Yeast and Zero Acre Farms respectively. R.P.D. and W.J.V. have financial interest in Future Genomics Technologies. All other authors have no competing interests., (Copyright © 2024. Published by Elsevier Inc.)

Zebrafish is a natural host of various Mycobacterium species and a surrogate model organism for tuberculosis research. Mycobacterium marinum is evolutionarily one of the closest non-tuberculous species related to M. tuberculosis and shares the majority of virulence genes. Although zebrafish is not a natural host of the human pathogen, we have previously demonstrated successful robotic infection of zebrafish embryos with M. tuberculosis and performed drug treatment of the infected larvae. In the present study, we examined for how long M. tuberculosis can be propagated in zebrafish larvae and tested a time series of infected larvae to study the transcriptional response via Illumina RNA deep sequencing (RNAseq). Bacterial aggregates carrying fluorescently labeled M. tuberculosis could be detected up to 9 days post-infection. The infected larvae showed a clear and specific transcriptional immune response with a high similarity to the inflammatory response of zebrafish larvae infected with the surrogate species M. marinum . We conclude that M. tuberculosis can be propagated in zebrafish larvae for at least one week after infection and provide further evidence that M. marinum is a good surrogate model for M. tuberculosis . The generated extensive transcriptome data sets will be of great use to add translational value to zebrafish as a model for infection of tuberculosis using the M. marinum infection system. In addition, we identify new marker genes such as dusp8 and CD180 that are induced by M. tuberculosis infection in zebrafish and in human macrophages at later stages of infection that can be further investigated.

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae , one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae , we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non- Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Vincent Courdavault reports financial support was provided by Horizon Europe. Michael Krogh Jensen reports a relationship with BioMIA that includes: board membership and employment. Ron Dirks reports a relationship with Future Genomics Technologies that includes: board membership and employment. Hans Jensen reports a relationship with Future Genomics Technologies that includes: board membership and employment. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

Evasion of apoptosis is critical for the development and growth of tumors. The pro-survival protein myeloid cell leukemia 1 (Mcl-1) is an antiapoptotic member of the Bcl-2 family, associated with tumor aggressiveness, poor survival, and drug resistance. Development of Mcl-1 inhibitors implies blocking of protein-protein interactions, generally requiring a lengthy optimization process of large, complex molecules. Herein, we describe the use of DNA-encoded chemical library synthesis and screening to directly generate complex, yet conformationally privileged macrocyclic hits that serve as Mcl-1 inhibitors. By applying a conceptual combination of conformational analysis and structure-based design in combination with a robust synthetic platform allowing rapid analoging, we optimized in vitro potency of a lead series into the low nanomolar regime. Additionally, we demonstrate fine-tuning of the physicochemical properties of the macrocyclic compounds, resulting in the identification of lead candidates 57 / 59 with a balanced profile, which are suitable for future development toward therapeutic use.

Monoterpene indole alkaloids (MIAs) are a structurally diverse family of specialized metabolites mainly produced in Gentianales to cope with environmental challenges. Due to their pharmacological properties, the biosynthetic modalities of several MIA types have been elucidated but not that of the yohimbanes. Here, we combine metabolomics, proteomics, transcriptomics and genome sequencing of Rauvolfia tetraphylla with machine learning to discover the unexpected multiple actors of this natural product synthesis. We identify a medium chain dehydrogenase/reductase (MDR) that produces a mixture of four diastereomers of yohimbanes including the well-known yohimbine and rauwolscine. In addition to this multifunctional yohimbane synthase (YOS), an MDR synthesizing mainly heteroyohimbanes and the short chain dehydrogenase vitrosamine synthase also display a yohimbane synthase side activity. Lastly, we establish that the combination of geissoschizine synthase with at least three other MDRs also produces a yohimbane mixture thus shedding light on the complex mechanisms evolved for the synthesis of these plant bioactives., (© 2023. The Author(s).)

Ovulation in European eel is induced by injection of 17α,20β-dihydroxy-4-pregnen-3-one (DHP) as the maturation-inducing hormone (MIH). Female eels need to ovulate within 18 h after injection to release good quality eggs. Progesterone (P), as an upstream precursor of DHP, may promote endogenous DHP production and improve egg quality. The purpose of this study was therefore to compare treatment of P with DHP on batch level, in vitro , to determine dose-response effects, and in vivo , at a single dose. For the in vitro experiment, ovarian tissue was extracted and placed in culture plates containing hormone-free medium and media supplemented with the treatment: DHP at 1, 10 and 100 ng mL -1 , or P at 10, 100 and 1,000 ng mL -1 . At the start of incubation, the folliculated oocytes were sampled for histology, microscopy and qPCR. After incubation for 12 and 18 h, the oocytes were sampled for microscopy and qPCR analysis. For the in vivo experiment, females were either injected with DHP or P at a dose of 2 mg kg -1 to assess their effects on ovulation and reproductive success. At the moment of release, eggs were sampled for RNA sequencing to compare effects of DHP and P on the expression of genes involved in egg quality aspects. Remaining eggs were fertilized and larval viability was recorded. Both DHP and P were able to induce GVBD (DHP at 10 and 100 ng mL -1 , P at 100 and 1,000 ng mL -1 ) in vitro . Expression of genes involved in oocyte maturation and ovulation was similar in vitro for both DHP and P treatments. Regarding the in vivo results, RNAseq results reflected similar DHP and P effects on the expression of genes involved in egg quality aspects. Females injected with either DHP or P ovulated, released eggs, and were equally able to produce larvae without any differences in reproductive success. Our results support the conclusion that DHP and P work equally well in vitro and in vivo . P is more attractive to apply as the price is 3,000 times lower than the price of DHP., Competing Interests: Author IG was employed by the company Rara Avis. Author RD was employed by the company Future Genomics Technologies B.V. Author WS was employed by the company Palingkwekerij Koolen B.V. Author LH was employed by the company Wageningen Eel Reproduction Experts B.V. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Jéhannet, Palstra, Meijerhof, Schipper, Giménez, Dirks, Swinkels, Heinsbroek and Komen.)

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

The Madagascar periwinkle, Catharanthus roseus , belongs to the Apocynaceae family. This medicinal plant, endemic to Madagascar, produces many important drugs including the monoterpene indole alkaloids (MIA) vincristine and vinblastine used to treat cancer worldwide. Here, we provide a new version of the C. roseus genome sequence obtained through the combination of Oxford Nanopore Technologies long-reads and Illumina short-reads. This more contiguous assembly consists of 173 scaffolds with a total length of 581.128 Mb and an N50 of 12.241 Mb. Using publicly available RNAseq data, 21,061 protein coding genes were predicted and functionally annotated. A total of 42.87% of the genome was annotated as transposable elements, most of them being long-terminal repeats. Together with the increasing access to MIA-producing plant genomes, this updated version should ease evolutionary studies leading to a better understanding of MIA biosynthetic pathway evolution., Competing Interests: Competing interests: Ron P. Dirks and Hans J. Jansen are CEO and CTO of Future Genomics Technologies, respectively., (Copyright: © 2022 Cuello C et al.)

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

The Apocynaceae tree Voacanga thouarsii, native to southern Africa and Madagascar, produces monoterpene indole alkaloids (MIA), which are specialized metabolites with a wide range of bioactive properties. Voacanga species mainly accumulates tabersonine in seeds making these species valuable medicinal plants currently used for industrial MIA production. Despite their importance, the MIA biosynthesis in Voacanga species remains poorly studied. Here, we report the first genome assembly and annotation of a Voacanga species. The combined assembly of Oxford Nanopore Technologies long-reads and Illumina short-reads resulted in 3,406 scaffolds with a total length of 1,354.26 Mb and an N50 of 3.04 Mb. A total of 33,300 protein-coding genes were predicted and functionally annotated. These genes were then used to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. A transposable element (TE) analysis showed the highest proportion of TE in Voacanga thouarsii compared with all other MIA-producing plants. In a nutshell, this first reference genome of V. thouarsii will thus contribute to strengthen future comparative and evolutionary studies in MIA-producing plants leading to a better understanding of MIA pathway evolution. This will also allow the potential identification of new MIA biosynthetic genes for metabolic engineering purposes., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Ca. Neoehrlichia mikurensis is widely prevalent in I. ricinus across Europe and has been associated with human disease. However, diagnostic modalities are limited, and much is still unknown about its biology. Here, we present the first complete Ca. Neoehrlichia mikurensis genomes directly derived from wildlife reservoir host tissues, using both long- and short-read sequencing technologies. This pragmatic approach provides an alternative to obtaining sufficient material from clinical cases, a difficult task for emerging infectious diseases, and to expensive and challenging bacterial isolation and culture methods. Both genomes exhibit a larger chromosome than the currently available Ca. Neoehrlichia mikurensis genomes and expand the ability to find new targets for the development of supportive laboratory diagnostics in the future. Moreover, this method could be utilized for other tick-borne pathogens that are difficult to culture.

Philasterides dicentrarchi is a scuticociliate that causes high mortalities in farmed fish. Although vaccination is an effective method to prevent scuticociliatosis caused by the homologous serotype, a universal vaccine has not been developed yet. Many compounds have been shown to be toxic to this ciliate species; moreover, most of them are toxic to aquatic life and cannot be used to prevent the disease. We have evaluated the toxicity to P. dicentrarchi of several compounds of natural origin to be used to reduce parasite levels in the seawater. Ciliates were exposed to several compound concentrations, and the mortality was determined at several incubation times. Tomatine, plumbagin and 2',4'-dihydroxychalcone displayed the highest anticiliate activity, with a dose-dependent response. The effects of these compounds on the EPC cell line were also evaluated, finding that 2',4'-dihydroxychalcone displayed the lowest toxicity to fish cells. At 7.54 μM, 2',4'-dihydroxychalcone inhibited 50% parasite growth but only killed about 10% of EPC cells after 24 h incubation. Finally, we evaluated the toxicity of Pseudomonas H6 surfactant (PS) to P. dicentrarchi , finding that PS was toxic to the ciliate but showed lower toxicity to EPC cells. At a concentration of 7.8 μg/mL (LC 50 for the ciliate after 3 h incubation), PS killed 14.9% of EPC cells. We conclude that 2',4'-dihydroxychalcone, and PS could be used to reduce parasite levels in seawater, thus decreasing the risk of scuticociliatosis infection in cultured fish.

Background: In physoclist fishes filling of the swimbladder requires acid secretion of gas gland cells to switch on the Root effect and subsequent countercurrent concentration of the initial gas partial pressure increase by back-diffusion of gas molecules in the rete mirabile. It is generally assumed that the rete mirabile functions as a passive exchanger, but a detailed analysis of lactate and water movements in the rete mirabile of the eel revealed that lactate is diffusing back in the rete. In the present study we therefore test the hypothesis that expression of transport proteins in rete capillaries allows for back-diffusion of ions and metabolites, which would support the countercurrent concentrating capacity of the rete mirabile. It is also assumed that in silver eels, the migratory stage of the eel, the expression of transport proteins would be enhanced., Results: Analysis of the transcriptome and of the proteome of rete mirabile tissue of the European eel revealed the expression of a large number of membrane ion and metabolite transport proteins, including monocarboxylate and glucose transport proteins. In addition, ion channel proteins, Ca 2+ -ATPase, Na + /K + -ATPase and also F 1 F 0 -ATP synthase were detected. In contrast to our expectation in silver eels the expression of these transport proteins was not elevated as compared to yellow eels. A remarkable number of enzymes degrading reactive oxygen species (ROS) was detected in rete capillaries., Conclusions: Our results reveal the expression of a large number of transport proteins in rete capillaries, so that the back diffusion of ions and metabolites, in particular lactate, may significantly enhance the countercurrent concentrating ability of the rete. Metabolic pathways allowing for aerobic generation of ATP supporting secondary active transport mechanisms are established. Rete tissue appears to be equipped with a high ROS defense capacity, preventing damage of the tissue due to the high oxygen partial pressures generated in the countercurrent system., (© 2021. The Author(s).)

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Several hypotheses explain the prevalence of undifferentiated sex chromosomes in poikilothermic vertebrates. Turnovers change the master sex determination gene, the sex chromosome or the sex determination system (e.g. XY to WZ). Jumping master genes stay main triggers but translocate to other chromosomes. Occasional recombination (e.g. in sex-reversed females) prevents sex chromosome degeneration. Recent research has uncovered conserved heteromorphic or even homomorphic sex chromosomes in several clades of non-avian and non-mammalian vertebrates. Sex determination in sturgeons (Acipenseridae) has been a long-standing basic biological question, linked to economical demands by the caviar-producing aquaculture. Here, we report the discovery of a sex-specific sequence from sterlet ( Acipenser ruthenus ). Using chromosome-scale assemblies and pool-sequencing, we first identified an approximately 16 kb female-specific region. We developed a PCR-genotyping test, yielding female-specific products in six species, spanning the entire phylogeny with the most divergent extant lineages ( A. sturio, A. oxyrinchus versus A. ruthenus, Huso huso ), stemming from an ancient tetraploidization. Similar results were obtained in two octoploid species ( A. gueldenstaedtii, A. baerii ). Conservation of a female-specific sequence for a long period, representing 180 Myr of sturgeon evolution, and across at least one polyploidization event, raises many interesting biological questions. We discuss a conserved undifferentiated sex chromosome system with a ZZ/ZW-mode of sex determination and potential alternatives. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

In eels, large variations in larval mortality exist, which would impede the viable production of juvenile glass eels in captivity. The transcriptome of European eel larvae was investigated to identify physiological pathways and genes that show differential regulation between non-viable vs. viable larvae. Expression of genes involved in inflammation and host protection was higher, suggesting that non-viable larvae suffered from microbial infection. Expression of genes involved in osmoregulation was also higher, implying that non-viable larvae tried to maintain homeostasis by strong osmoregulatory adaptation. Expression of genes involved in myogenesis, neural, and sensory development was reduced in the non-viable larvae. Expression of the major histocompatibility complex class-I ( mhc1 ) gene, M-protein ( myom2 ), the dopamine 2B receptor ( d2br ), the melatonin receptor ( mtr1 ), and heat-shock protein beta-1 ( hspb1 ) showed strong differential regulation and was therefore studied in 1, 8, and 15 days post-hatch (dph) larvae by RT-PCR to comprehend the roles of these genes during ontogeny. Expression patterning of these genes indicated the start of active swimming (8 dph) and feed searching behavior (15 dph) and confirmed immunocompetence immediately after hatching. This study revealed useful insights for improving larval survival by microbial control and salinity reduction.

The ectoparasite protozoan Amyloodinium ocellatum (AO) is the causative agent of amyloodiniosis in European seabass (ESB, Dicentrarchus labrax). There is a lack of information about basic molecular immune response mechanisms of ESB during AO infestation. Therefore, to compare gene expression between experimental AO-infested ESB tissues and uninfested ESB tissues (gills and head kidney) RNA-seq was adopted. The RNA-seq revealed multiple differentially expressed genes (DEG), namely 679 upregulated genes and 360 downregulated genes in the gills, and 206 upregulated genes and 170 downregulated genes in head kidney. In gills, genes related to the immune system (perforin, CC1) and protein binding were upregulated. Several genes involved in IFN related pathways were upregulated in the head kidney. Subsequently, to validate the DEG from amyloodiniosis, 26 ESB (mean weight 14 g) per tank in triplicate were bath challenged for 2 h with AO (3.5 × 10 6 /tank; 70 dinospores/mL) under controlled conditions (26-28 °C and 34‰ salinity). As a control group (non-infested), 26 ESB per tank in triplicate were also used. Changes in the expression of innate immune genes in gills and head kidney at 2, 3, 5, 7 and 23 dpi were analysed using real-time PCR. The results indicated that the expression of cytokines (CC1, IL-8) and antimicrobial peptide (Hep) were strongly stimulated and reached a peak at 5 dpi in the early infestation stage, followed by a gradual reduction in the recovery stage (23 dpi). Noticeably, the immunoglobulin (IgM) expression was higher at 23 dpi compared to 7 dpi. Furthermore, in-situ hybridization showed positive signals of CC1 mRNA in AO infested gills compared to the control group. Altogether, chemokines were involved in the immune process under AO infestation and this evidence allows a better understanding of the immune response in European seabass during amyloodiniosis., (Copyright © 2021. Published by Elsevier B.V.)

Ceratothoa oestroides (Cymothoidea, Isopoda) is a generalist crustacean parasite that negatively affects the economic sustainability of European sea bass ( Dicentrarchus labrax ) aquaculture in the North-East Mediterranean. While mortalities are observed in fry and fingerlings, infection in juvenile and adult fish result in approximately 20% growth delay. A transcriptomic analysis (PCR array, RNA-Seq) was performed on organs (tongue, spleen, head kidney, and liver) from infected vs. Ceratothoa -free sea bass fingerlings. Activation of local and systemic immune responses was detected, particularly in the spleen, characterized by the upregulation of cytokines (also in the tongue), a general reshaping of the immunoglobulin (Ig) response and suppression of T-cell mediated responses. Interestingly, starvation and iron transport and metabolism genes were strongly downregulated, suggesting that the parasite feeding strategy is not likely hematophagous. The regulation of genes related to growth impairment and starvation supported the growth delay observed in infected animals. Most differentially expressed (DE) transcripts were exclusive of a specific organ; however, only in the tongue, the difference between infected and uninfected fish was significant. At the attachment/feeding site, the pathways involved in muscle contraction and intercellular junction were the most upregulated, whereas the pathways involved in fibrosis (extracellular matrix organization, collagen formation, and biosynthesis) were downregulated. These results suggest that parasite-inflicted damage is successfully mitigated by the host and characterized by regenerative processes that prevail over the reparative ones., Competing Interests: RD and ES are employed by the company Future Genomics Technologies. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Piazzon, Mladineo, Dirks, Santidrián Yebra-Pimentel, Hrabar and Sitjà-Bobadilla.)

A longer on-land rearing period of Gilthead seabream Sparus aurata before transfer to sea-cages would allow the farmer to benefit from exercise-enhanced growth, resilience, and robustness as induced by increasing water flow in the tanks. In this study, the physiological effects of flow-conditioning were investigated by subjecting large groups of experimental fish to minimal flow or to flow regimes inducing swimming exercise at 1 or 2 body length (BL) s -1 for a period of 8 months (February-October) in 1,500 L tanks. Fish representing the three treatment groups were then used for: (1) a stress challenge netting test and plasma cortisol measurement (baseline, peaking, and recovery levels), (2) blood plasma measurements of glucose, triglycerides, lactate, cholesterol, growth hormone (GH), and insulin-like growth factor 1 (IGF1), and (3) heart and muscle gene expression of the GH and IGF1 receptors and the muscle transcriptome by deep RNA sequencing (RNAseq). Fish size after 8 months of flow conditioning was 92 ± 27 g body weight (BW) for fish under minimal flow, 106 ± 24 g BW (+15%) at 1 BL s -1 , and 125 ± 27 g BW (+36%) at 2 BL s -1 . Flow conditioning at 1 BL s -1 provided optimal conditions for growth and uniformity, but also stress (lowest baseline plasma cortisol), robustness (higher condition factor and larger hearts), and energy mobilization (increased plasma glucose). Although flow enhanced growth linearly with swimming speed, also the percentage of lordotic fish increased with exercise, particularly high for swimming at 2 BL s -1 . The absence of important differences in plasma GH and IGF1, and expression levels of their receptors in heart and white skeletal muscle, indicated that other factors may be involved in growth enhancement. RNAseq of the white skeletal muscle showed upregulated expression of genes involved in muscle contraction, muscle development and its molecular regulation, and immune genes that may play a role in the muscle repair mechanism. An exercise regime of swimming at 1 BL s -1 can be considered as optimal for farming robust seabream although the increase of skeletal deformities should be avoided., Competing Interests: Author RD was employed by the company Future Genomics Technologies B.V. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2020 Palstra, Roque, Kruijt, Jéhannet, Pérez-Sánchez and Dirks.)

The lesser periwinkle Vinca minor accumulates numerous monoterpene indole alkaloids (MIAs) including the vasodilator vincamine. While the biosynthetic pathway of MIAs has been largely elucidated in other Apocynaceae such as Catharanthus roseus , the counterpart in V. minor remains mostly unknown, especially for reactions leading to MIAs specific to this plant. As a consequence, we generated a comprehensive V. minor transcriptome elaborated from eight distinct samples including roots, old and young leaves exposed to low or high light exposure conditions. This optimized resource exhibits an improved completeness compared to already published ones. Through homology-based searches using C. roseus genes as bait, we predicted candidate genes for all common steps of the MIA pathway as illustrated by the cloning of a tabersonine/vincadifformine 16- O -methyltransferase (Vm16OMT) isoform. The functional validation of this enzyme revealed its capacity of methylating 16-hydroxylated derivatives of tabersonine, vincadifformine and lochnericine with a Km 0.94 ± 0.06 µM for 16-hydroxytabersonine. Furthermore, by combining expression of fusions with yellow fluorescent proteins and interaction assays, we established that Vm16OMT is located in the cytosol and forms homodimers. Finally, a gene co-expression network was performed to identify candidate genes of the missing V. minor biosynthetic steps to guide MIA pathway elucidation.

Sturgeons are being used in aquaculture because wild populations are now endangered due to overfishing for caviar. A challenge in working with sturgeon as an aquacultured species is its long and slow reproductive development. Reproduction is a hormonally regulated process that involves hierarchical signaling between the brain, pituitary gland, and gonads. In an effort to better understand the hormonal regulation of sturgeon reproduction, we have cloned the Russian sturgeon (st), Acipenser gueldenstaedtii, luteinizing hormone receptor (stLHR) and follicle stimulating hormone receptor (stFSHR) and measured their expression from previtellogenic to mature ovarian follicles. Sturgeon LHR and FSHR expression was elevated in early-vitellogenic and mature follicles compared with pre-vitellogenic and mid-vitellogenic follicles, and only LHR expression increased during late-vitellogenesis. Recombinant sturgeon FSH and LH both activated sturgeon LHR and FSHR in a cAMP reporter assay. Further molecular characterization of these receptors was accomplished by in silico modeling and cAMP reporter assays using heterologous recombinant gonadotropins from human and piscine species. There was no apparent trend in heterologous LH and/or FSH activation of the sturgeon LHR or FSHR. These data suggest that permissive activation of LHR and FSHR are a consequence of some yet undetermined biological characteristic(s) of different piscine species., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The ectoparasite protozoan Amyloodinium ocellatum (AO) is the etiological agent of amyloodiniosis in European seabass ( Dicentrarchus labrax ) (ESB). There is a lack of information about basic molecular data on AO biology and its interaction with the host. Therefore, de novo transcriptome sequencing of AO tomonts was performed. AO trophonts were detached from infested ESB gills, and quickly becoming early tomonts were purified by Percoll ® density gradient. Tomont total RNA was processed and quality was assessed immediately. cDNA libraries were constructed using TruSeq ® Stranded mRNA kit and sequenced using Illumina sequencer. CLC assembly was used to generate the Transcriptome assembly of AO tomonts. Out of 48,188 contigs, 56.12% belong to dinophyceae wherein Symbiodinium microadriaticum had 94.61% similarity among dinophyceae. Functional annotations of contigs indicated that 12,677 had associated GO term, 9005 with KEGG term. The contigs belonging to dinophyceae resulted in the detection of several peptidases. A BLAST search for known virulent factors from the virulence database resulted in hits to Rab proteins, AP120, Ribosomal phosphoprotein, Heat-shock protein70, Casein kinases, Plasmepsin IV, and Brucipain. Hsp70 and casein kinase II alpha were characterized in-silico. Altogether, these results provide a reference database in understanding AO molecular biology, aiding to the development of novel diagnostics and future vaccines.

The present study analyses the interactions between Philasterides dicentrarchi (a ciliate parasite that causes high mortalities in cultured flatfish) and the peritoneal cells of the turbot Scophthalmus maximus during an experimental infection. The transcriptomic response was evaluated in the parasites and in the fish peritoneal cells, at 1, 2 and 4 h post-infection (hpi) in turbot injected intraperitoneally (ip) with 10 7 ciliates and at 12 and 48 hpi in turbot injected ip with 10 5 ciliates. Numerous genes were differentially expressed (DE) in P. dicentrarchi , relative to their expression in control ciliates (0 hpi): 407 (369 were up-regulated) at 1 hpi, 769 (415 were up-regulated) at 2 hpi and 507 (119 were up-regulated) at 4 hpi. Gene ontology (GO) analysis of the DE genes showed that the most representative categories of biological processes affected at 1, 2 and 4 hpi were biosynthetic processes, catabolic processes, biogenesis, proteolysis and transmembrane transport. Twelve genes of the ABC transporter family and eight genes of the leishmanolysin family were DE at 1, 2 and 4 hpi. Most of these genes were strongly up-regulated (UR), suggesting that they are involved in P. dicentrarchi infection. A third group of UR genes included several genes related to ribosome biogenesis, DNA transcription and RNA translation. However, expression of tubulins and tubulin associated proteins, such as kinesins or dyneins, which play key roles in ciliate division and movement, was down-regulated (DR). Similarly, genes that coded for lysosomal proteins or that participate in the cell cycle mitotic control, glycolysis, the Krebs cycle and/or in the electron transport chain were also DR. The transcriptomic analysis also revealed that in contrast to many parasites, which passively evade the host immune system, P. dicentrarchi strongly stimulated turbot peritoneal cells. Many genes related to inflammation were DE in peritoneal cells at 1, 2 and 4 hpi. However, the response was much lower at 12 hpi and almost disappeared completely at 48 hpi in fish that were able to kill P. dicentrarchi during the first few hpi. The genes that were DE at 1, 2 and 4 hpi were mainly related to the apoptotic process, the immune response, the Fc-epsilon receptor signalling pathway, the innate immune response, cell adhesion, cell surface receptors, the NF-kappaB signalling pathway and the MAPK cascade. Expression of toll-like receptors 2, 5 and 13 and of several components of NF-κB, MAPK and JAK/STAT signalling pathways was UR in the turbot peritoneal cells. Genes expressing chemokines and chemokine receptors, genes involved in prostaglandin and leukotriene synthesis, prostaglandins, leukotriene receptors, proinflammatory cytokines and genes involved in apoptosis were strongly UR during the first four hours of infection. However, expression of anti-inflammatory cytokines such as Il-10 and lipoxygenases with anti-inflammatory activity (i.e., arachidonate 15-lipoxygenase ) were only UR at 12 and/or 48 hpi, indicating an anti-inflammatory state in these groups of fish. In conclusion, the present study shows the regulation of several genes in P. dicentrarchi during the early stages of infection, some of which probably play important roles in this process. The infection induced a potent acute inflammatory response, and many inflammatory genes were regulated in peritoneal cells, showing that the turbot uses all the protective mechanisms it has available to prevent the entry of the parasite.

Exposure to high temperatures can lead to thermotolerance in fish, which is hypothesized to potentially improve post-release survival in species under restocking programs, like Atlantic sturgeon. The aim of this study was to determine whether Atlantic sturgeon juveniles exposed to a 4-week temperature treatment respond differently to a subsequent heat shock than juveniles exposed to heat shock for the first time (naive fish). Response to heat shock was assessed by mapping the liver transcriptome. In total, 838 unique contigs were differentially expressed between the trained and the control group (592 downregulated, 261 upregulated, and 15 down- or upregulated, depending on the condition), corresponding to genes involved in the response to heat, tissue damage, proteolysis, and metabolism. Temperature-trained fish showed 2-4-fold fewer dysregulated contigs than naive fish, indicating their ability to maintain and recover homeostasis faster. During heat shock, hspc1 was upregulated in both experimental groups, while hspa1 and dnaja4 were exclusively upregulated in the control. Overall, compensatory mechanisms were observed in addition to the heat shock response. Only two genes, fgg and apnl, were upregulated at nearly all timepoints in both groups. Peptidases were more strongly downregulated in control fish, which also showed a reduction in lipid metabolism during recovery. Keratins, pck1, gadd45ga, and gadd45gb were differentially expressed between trained and control fish, and due to their roles in tissue protection and ER stress reduction, they might be responsible for the maintenance of the transcriptional homeostasis observed in trained fish.

In spite of many decades of research, the spawning migration of the European eel Anguilla anguilla from the European coast to the Sargasso Sea remains a mystery. In particular, the role of the swimbladder as a buoyancy regulating structure is not yet understood. In this study, we exercised silver eels in a swim tunnel under elevated hydrostatic pressure. The transcriptome of gas gland tissue of these exercised eels was then compared to the known transcriptome of not exercised (control) silver eel gas gland cells. Due to the high infection rate of the eel population with the swimbladder parasite Anguillicola crassus, the comparison also included an exercised group of silver eels with a heavily damaged swimbladder, and we compared the previously published transcriptome of not exercised silver eels with a highly damaged swimbladder with the exercised group of silver eels with a heavily damaged swimbladder. The comparisons of unexercised (control) silver eels with exercised silver eels with functional swimbladder (EF), as well as with exercised silver eels with damaged swimbladder (ED), both showed a significant elevation in transcripts related to glycolytic enzymes. This could also be observed within the comparison of unexercised silver eels with a highly infected swimbladder with exercised eels with a damaged swimbladder (DED). In contrast to EF, in ED a significant elevation in transcript numbers of mitochondrial NADH dehydrogenase was observed. While in EF the transcriptional changes suggested that acid production and secretion was enhanced, in ED these changes appeared to be related to thickened tissue and thus elevated diffusion distances. The remarkable number of differentially expressed transcripts coding for proteins connected to cAMP-dependent signaling pathways indicated that metabolic control in gas gland cells includes cAMP-dependent pathways. In contrast to ED, in EF significant transcriptional changes could be related to the reconstruction of the extracellular matrix, while in ED tissue repair and inflammation was more pronounced. Surprisingly, in exercised eels hypoxia inducible transcription factor expression was elevated. In EF, a large number of genes related to the circadian clock were transcriptionally modified, which may be connected to the circadian vertical migrations observed during the spawning migration., Competing Interests: R.D. and E.S. are employees of Future Genomics Technologies, a commercial company. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

The fish parasites Saprolegnia spp. (Oomycota) and Amyloodinium ocellatum (Dinophyceae) cause important losses in freshwater and marine aquaculture industry, respectively. The possible adverse effects of compounds used to control these parasites in aquaculture resulted in increased interest on the search for natural products with antiparasitic activity. In this work, eighteen plant-derived compounds (2',4'-Dihydroxychalcone; 7-Hydroxyflavone; Artemisinin; Camphor (1R); Diallyl sulfide; Esculetin; Eucalyptol; Garlicin 80%; Harmalol hydrochloride dihydrate; Palmatine chloride; Piperine; Plumbagin; Resveratrol; Rosmarinic acid; Sclareolide; Tomatine, Umbelliferone, and Usnic Acid) have been tested in vitro . Sixteen of these were used to determine their effects on the gill cell line G1B (ATCC®CRL-2536™) and on the motility of viable dinospores of Amyloodinium ocellatum , and thirteen were screened for inhibitory activity against Saprolegnia spp. The cytotoxicity results on G1B cells determined that only two compounds (2',4'-Dihydroxychalcone and Tomatine) exhibited dose-dependent toxic effects. The highest surveyed concentrations (0.1 and 0.01 mM) reduced cell viability by 80%. Upon lowering the compound concentration the percentage of dead cells was lower than 20%. The same two compounds revealed to be potential antiparasitics by reducing in a dose-dependent manner the motility of A. ocellatum dinospores up to 100%. With respect to Saprolegnia , a Minimum Inhibitory Concentration was found for Tomatine (0.1 mM), Piperine and Plumbagin (0.25 mM), while 2',4'-Dihydroxychalcone considerably slowed down mycelial growth for 24 h at a concentration of 0.1 mM. Therefore, this research allowed to identify two compounds, Tomatine and 2',4'-Dihydroxychalcone, effective against both parasites. These compounds could represent promising candidates for the treatment of amyloodiniosis and saprolegniosis in aquaculture. Nevertheless, further in vitro and in vivo tests are required in order to determine concentrations that are effective against the considered pathogens but at the same time safe for hosts, environment and consumers., (Copyright © 2020 Tedesco, Beraldo, Massimo, Fioravanti, Volpatti, Dirks and Galuppi.)

Warfarin is the most worldwide used anticoagulant drug and rodenticide. Since it crosses placental barrier it can induce warfarin embryopathy (WE), a fetal mortality in neonates characterized by skeletal deformities in addition to brain hemorrhages. Although the effects of warfarin exposure in aquatic off target species were already described, the particular molecular toxicological mechanisms during early development are still unclear. Here, we used zebrafish (Danio rerio) to describe and compare the developmental effects of warfarin exposure (0, 15.13, 75.68 and 378.43 mM) on two distinct early developmental phases (embryos and eleuthero-embryos). Although exposure to both developmental phases induced fish mortality, only embryos exposed to the highest warfarin level exhibited features mimicking mammalian WE, e.g. high mortality, higher incidence of hemorrhages and altered skeletal development, among other effects. To gain insights into the toxic mechanisms underlying warfarin exposure, the transcriptome of embryos exposed to warfarin was explored through RNA-Seq and compared to that of control embryos. 766 differentially expressed (564 up- and 202 down-regulated) genes were identified. Gene Ontology analysis revealed particular cellular components (cytoplasm, extracellular matrix, lysosome and vacuole), biological processes (mainly amino acid and lipid metabolism and response to stimulus) and pathways (oxidative stress response and apoptosis signaling pathways) being significantly overrepresented in zebrafish embryos upon warfarin exposure. Protein-protein interaction further evidenced an altered redox system, blood coagulation and vasculogenesis, visual phototransduction and collagen formation upon warfarin exposure. The present study not only describes for the first time the WE in zebrafish, it provides new insights for a better risk assessment, and highlights the need for programming the rat eradication actions outside the fish spawning season to avoid an impact on off target fish community. The urge for the development of more species-specific anticoagulants for rodent pest control is also highlighted., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Despite efforts to restore Atlantic sturgeon in European rivers, aquaculture techniques result in animals with high post-release mortality due to, among other reasons, their low tolerance to increasing water temperature. Marker genes to monitor heat stress are needed in order to identify heat-resistant fish. Therefore, an Atlantic sturgeon cell line was exposed to different heat shock protocols (30 °C and 35 °C) and differences in gene expression were investigated. In total 3020 contigs (∼1.5%) were differentially expressed. As the core of the upregulated contigs corresponded to heat shock proteins (HSP), the heat shock factor (HSF) and the HSP gene families were annotated in Atlantic sturgeon and mapped via Illumina RNA sequencing to identify heat-inducible family members. Up to 6 hsf and 76 hsp genes were identified in the Atlantic sturgeon transcriptome resources, 16 of which were significantly responsive to the applied heat shock. The previously studied hspa1 (hsp70) gene was only significantly upregulated at the highest heat shock (35 °C), while a set of 5 genes (hspc1, hsph3a, hspb1b, hspb11a, and hspb11b) was upregulated at all conditions. Although the hspc1 (hsp90a) gene was previously used as heat shock-marker in sturgeons, we found that hspb11a is the most heat-inducible gene, with up to 3296-fold higher expression in the treated cells, constituting the candidate gene markers for in vivo trials., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Smoltification is a metamorphic event in salmon life history, which initiates downstream migration and pre-adapts juvenile salmon for seawater entry. While a number of reports concern thyroid hormones and smoltification, few and inconclusive studies have addressed the potential role of thyrotropin (TSH). TSH is composed of a α-subunit common to gonadotropins, and a β-subunit conferring hormone specificity. We report the presence and functional divergence of duplicated TSH β-subunit paralogs (tshβa and tshβb) in Atlantic salmon. Phylogeny and synteny analyses allowed us to infer that they originated from teleost-specific whole genome duplication. Expression profiles of both paralogs in the pituitary were measured by qPCR throughout smoltification in Atlantic salmon from the endangered Loire-Allier population raised in a conservation hatchery. This revealed a striking peak of tshβb expression in April, concomitant with downstream migration initiation, while tshβa expression remained relatively constant. In situ hybridization showed two distinct pituitary cell populations, tshβa cells in the anterior adenohypophysis, and tshβb cells near to the pituitary stalk, a location comparable to the pars tuberalis TSH cells involved in seasonal physiology and behaviour in birds and mammals. Functional divergence of tshβ paralogs in Atlantic salmon supports a specific role of tshβb in smoltification.

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.

The parasite Ichthyophthirius multifiliis infecting skin, fins and gills of a wide range of freshwater fish species, including rainbow trout, is known to induce a protective immune response in the host. Although a number of studies have reported activation of several immune genes in infected fish host, the immune response picture is still considered incomplete. In order to address this issue, a comparative transcriptomic analysis was performed on infected versus uninfected rainbow trout gills and it showed that a total of 3352 (7.2%) out of 46,585 identified gene sequences were significantly regulated after parasite infection. Of differentially expressed gene sequences, 1796 genes were up-regulated and 1556 genes were down-regulated. These were classified into 61 Gene Ontology (GO) terms and mapped to 282 reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Infection of I. multifiliis induced a clear differential expression of immune genes, related to both innate and adaptive immunity. A total of 268 (6.86%) regulated gene sequences were known to take part in 16 immune-related pathways. These involved pathways related to the innate immunity such as the Chemokine signaling pathway, Platelet activation, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, and Leukocyte transendothelial migration. Elevated transcription of genes encoding the TLR 8 gene and chemokines (CCL4, CCL19, CCL28, CXCL8, CXCL11, CXCL13, CXCL14) was recorded indicating their roles in recognition of I. multifiliis and subsequent induction of the inflammatory response, respectively. A number of upregulated genes in infected gills were associated with antigen processing/presentation and T and B cell receptor signaling (including B cell marker CD22 involved in B cell development). Overall the analysis supports the notion that I. multifiliis induces a massive and varied innate response upon which a range of adaptive immune responses are established which may contribute to the long lasting protection of immunized rainbow trout., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Forced sustained swimming exercise at optimal speed enhances growth in many fish species, particularly through hypertrophy of the white skeletal muscle. The exact mechanism of this effect has not been resolved yet. To explore the role of cortisol, we first subjected wild-type zebrafish to an exercise protocol validated for exercise-enhanced growth, and showed that exercised zebrafish, which indeed showed enhanced growth, had higher cortisol levels than the non-exercised controls. A central role was therefore hypothesized for the steroid hormone cortisol acting through the Glucocorticoid receptor (Gr). Second, we subjected wild-type zebrafish and zebrafish with a mutant Gr to exercise at optimal, suboptimal, and super-optimal speeds and compared them with non-exercised controls. Exercised zebrafish showed growth enhancement at all speeds, with highest growth at optimal speeds. In the Gr mutant fish, exercise resulted in growth enhancement similar to wild-type zebrafish, indicating that cortisol signaling through Gr cannot be considered as a main determinant of exercise-enhanced growth. Finally, the transcriptome of white skeletal muscle tissue was analyzed by RNA sequencing. The results of this analysis showed that in the muscle tissue of Gr mutant fish a lower number of genes is regulated by exercise than in wild-type fish (183 vs. 351). A cluster of 36 genes was regulated by exercise in both wild-type and mutant fish, and in this cluster genes involved in transcriptional regulation and protein ubiquitination were overrepresented. Because these two processes appear to be regulated in both wild type and mutant fish, which both display exercise-enhanced growth, we suggest that they play an important role in the growth of muscles upon exercise.

Seaweeds may represent immuno-stimulants that could be used as health-promoting fish feed components. This study was performed to gain insights into the immunomodulatory effects of dietary seaweeds in Atlantic salmon. Specifically tested were 10% inclusion levels of Laminaria digitata (SW1) and a commercial blend of seaweeds (Oceanfeed ® ) (SW2) against a fishmeal based control diet (FMC). Differences between groups were assessed in growth, feed conversion ratio and blood parameters hematocrit and hemoglobin. After a LPS challenge of fish representing each of the three groups, RNAseq was performed on the head kidney as major immune organ to determine transcriptomic differences in response to the immune activation. Atlantic salmon fed with dietary seaweeds did not show major differences in performance in comparison with fishmeal fed fish. RNAseq resulted in ∼154 million reads which were mapped against a NCBI Salmo salar reference and against a de novo assembled S. salar reference for analyses of expression of immune genes and ontology of immune processes among the 87,600 cDNA contigs. The dietary seaweeds provoked a more efficient immune response which involved more efficient identification of the infection site, and processing and presentation of antigens. More specifically, chemotaxis and the chemokine-mediated signaling were improved and therewith the defense response to Gram-positive bacterium reduced. Specific Laminaria digitata effects included reduction of the interferon-gamma-mediated signaling. Highly upregulated and specific for this diet was the expression of major histocompatibility complex class I-related gene protein . The commercial blend of seaweeds caused more differential expression than Laminaria digitata and improved immune processes such as receptor-mediated endocytosis and cell adhesion , and increased the expression of genes involved in response to lipopolysaccharide and inflammatory response . Particularly, expression of many important immune receptors was up-regulated illustrating increased responsiveness. NF-kappa-B inhibitor alpha is an important gene that marked the difference between both seaweed diets as Laminaria digitata inhibits the expression for this cytokine while the blend of seaweeds stimulates it. It can be concluded that the inclusion of seaweeds such as Laminaria digitata can have important modulatory effects on the immune capacity of Atlantic salmon resulting in a more efficient immune response.

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., (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

To enable selection of novel chemicals for new processes, there is a recognized need for alternative toxicity screening assays to assess potential risks to man and the environment. For human health hazard assessment these screening assays need to be translational to humans, have high throughput capability, and from an animal welfare perspective be harmonized with the principles of the 3Rs (Reduction, Refinement, Replacement). In the area of toxicology a number of cell culture systems are available but while these have some predictive value, they are not ideally suited for the prediction of developmental and reproductive toxicology (DART). This is because they often lack biotransformation capacity, multicellular or multi- organ complexity, for example, the hypothalamus pituitary gonad (HPG) axis and the complete life cycle of whole organisms. To try to overcome some of these limitations in this study, we have used Caenorhabditis elegans (nematode) and Danio rerio embryos (zebrafish) as alternative assays for DART hazard assessment of some candidate chemicals being considered for a new commercial application. Nematodes exposed to Piperazine and one of the analogs tested showed a slight delay in development compared to untreated animals but only at high concentrations and with Piperazine as the most sensitive compound. Total brood size of the nematodes was also reduced primarily by Piperazine and one of the analogs. In zebrafish Piperazine and analogs showed developmental delays. Malformations and mortality in individual fish were also scored. Significant malformations were most sensitively identified with Piperazine, significant mortality was only observed in Piperazine and only at the higest dose. Thus, Piperazine seemed the most toxic compound for both nematodes and zebrafish. The results of the nematode and zebrafish studies were in alignment with data obtained from conventional mammalian toxicity studies indicating that these have potential as developmental toxicity screening systems. The results of these studies also provided reassurance that none of the Piperazines tested are likely to have any significant developmental and/or reproductive toxicity issues to humans when used in their commercial applications., (Copyright © 2017. Published by Elsevier Ltd.)

Using Illumina sequencing, we investigated transcriptional changes caused by the nematode Anguillicola crassus within yellow and silver eels by comparing swimbladder samples of uninfected yellow with infected yellow eels, and uninfected silver with infected silver eels, respectively. In yellow eel gas gland, the infection caused a modification of steady state mRNA levels of 1675 genes, most of them being upregulated. Functional annotation analysis based on GO terms was used to categorize identified genes with regard to swimbladder metabolism or response to the infection. In yellow eels, the most prominent category was 'immune response', including various inflammatory components, complement proteins, and immunoglobulins. The elevated expression of several glucose and monocarboxylate transporters indicated an attempt to maintain the level of glucose metabolism, even in due to the infection thickened swimbladder tissue. In silver eel swimbladder tissue, on the contrary, the mRNA levels of only 291 genes were affected. Genes in the categories 'glucose metabolism' and 'ROS metabolism' barely responded to the infection and even the reaction of the immune system was much less pronounced compared to infected yellow eels. However, in the category 'extracellular matrix', the mRNA levels of several mucin genes were strongly elevated, suggesting increased mucus production as a defense reaction against the parasite. The present study revealed a strong reaction to an Anguillicola crassus infection on mRNA expression levels in swimbladder tissue of yellow eels, whereas in silver eels the changes ware almost negligible. A possible explanation for this difference is that the silvering process requires so much energy that there is not much scope to cope with the additional challenge of a nematode infection. Another possible explanation could be that gas-secreting activity of the silver eel swimbladder was largely reduced, which could coincide with a reduced responsiveness to other challenges, like a nematode infection.

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.

Background: The introduction of the MinION TM sequencing device by Oxford Nanopore Technologies may greatly accelerate whole genome sequencing. It has been shown that the nanopore sequence data, in combination with other sequencing technologies, is highly useful for accurate annotation of all genes in the genome. However, it also offers great potential for de novo assembly of complex genomes without using other technologies. In this manuscript we used nanopore sequencing as a tool to classify yeast strains. Methods: We compared various technical and software developments for the nanopore sequencing protocol, showing that the R9 chemistry is, as predicted, higher in quality than R7.3 chemistry. The R9 chemistry is an essential improvement for assembly of the extremely AT-rich mitochondrial genome. Results: In this study, we used this new technology to sequence and de novo assemble the genome of a recently isolated ethanologenic yeast strain, and compared the results with those obtained by classical Illumina short read sequencing. This strain was originally named Candida vartiovaarae ( Torulopsis vartiovaarae ) based on ribosomal RNA sequencing. We show that the assembly using nanopore data is much more contiguous than the assembly using short read data. Conclusions: The mitochondrial and chromosomal genome sequences showed that our strain is clearly distinct from other yeast taxons and most closely related to published Cyberlindnera species. In conclusion, MinION-mediated long read sequencing can be used for high quality de novo assembly of new eukaryotic microbial genomes., Competing Interests: Competing interests: HJJ and CVH are members of the Nanopore Community, and have previously received flow cells free of charge, as well as travel expense reimbursements from Oxford Nanopore Technologies.

Background : The introduction of the MinION sequencing device by Oxford Nanopore Technologies may greatly accelerate whole genome sequencing. Nanopore sequence data offers great potential for de novo assembly of complex genomes without using other technologies. Furthermore, Nanopore data combined with other sequencing technologies is highly useful for accurate annotation of all genes in the genome. In this manuscript we used nanopore sequencing as a tool to classify yeast strains. Methods : We compared various technical and software developments for the nanopore sequencing protocol, showing that the R9 chemistry is, as predicted, higher in quality than R7.3 chemistry. The R9 chemistry is an essential improvement for assembly of the extremely AT-rich mitochondrial genome. We double corrected assemblies from four different assemblers with PILON and assessed sequence correctness before and after PILON correction with a set of 290 Fungi genes using BUSCO. Results : In this study, we used this new technology to sequence and de novo assemble the genome of a recently isolated ethanologenic yeast strain, and compared the results with those obtained by classical Illumina short read sequencing. This strain was originally named Candida vartiovaarae ( Torulopsis vartiovaarae ) based on ribosomal RNA sequencing. We show that the assembly using nanopore data is much more contiguous than the assembly using short read data. We also compared various technical and software developments for the nanopore sequencing protocol, showing that nanopore-derived assemblies provide the highest contiguity. Conclusions : The mitochondrial and chromosomal genome sequences showed that our strain is clearly distinct from other yeast taxons and most closely related to published Cyberlindnera species. In conclusion, MinION-mediated long read sequencing can be used for high quality de novo assembly of new eukaryotic microbial genomes., Competing Interests: Competing interests: HJJ and CVH are members of the Nanopore Community, and have previously received flow cells free of charge, as well as travel expense reimbursements from Oxford Nanopore Technologies.

This study evaluates the effects of temperature on hCG-induced spermatogenesis in European eel (Anguilla anguilla), subjected to three thermal regimes: T10: 10°C (first 4weeks), 15°C (next 3weeks) and 20°C (last 6weeks); T15: 15°C (first 4weeks) and 20°C (last 9weeks); and T20: constant 20°C for the duration of the experiment. At 10°C, maturation stopped in the A spermatogonial stage (SPG1), and no further maturation was observed until the temperature was ≥15°C. With the aim of explaining these results, the influence of temperature on steroidogenic enzyme gene expression and steroid synthesis was tested. The initial synthesis of androgens (T and 11-KT) increased at SPG1, and was not influenced by temperature. Likewise, the gene expression of the steroidogenic enzymes linked to androgen synthesis (aacyp11a1, aacyp17-I and aa11βHSD) also increased at SPG1. In contrast, no correlation was seen between the increase in E2 and the aacyp19a1 gene expression peak in the testes, with E2 increasing as a consequence of the seawater acclimation carried out before hormonal treatment, and peaking the aacyp19a1 gene expression at B spermatogonial stage (SPG2). Aacyp21 gene expression was also higher at SPG2, and this stage was only reached when the rearing temperature was ≥15°C. In conclusion, androgen synthesis is not dependent on temperature, but further maturation requires higher temperatures in order to induce a change in the steroidogenic pathway towards estrogen and progestin synthesis. This study demonstrates that temperature plays a crucial role in European eel maturation, even perhaps controlling gonad development during the reproductive migration., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Using Illumina sequencing, transcriptional changes occurring during silvering in swimbladder tissue of the European eel have been analyzed by comparison of yellow and silver eel tissue samples. Functional annotation analysis based on GO terms revealed significant expression changes in a number of genes related to the extracellular matrix, important for the control of gas permeability of the swimbladder, and to reactive oxygen species (ROS) defense, important to cope with ROS generated under hyperbaric oxygen partial pressures. Focusing on swimbladder tissue metabolism, levels of several mRNA species encoding glucose transport proteins were several-fold higher in silver eels, while enzymes of the glycolytic pathway were not affected. The significantly higher steady state level of a transcript encoding for membrane bound carbonic anhydrase, however, suggested that CO2 production in the pentose phosphate shunt and diffusion of CO2 was of particular importance in silver eel swimbladder. In addition, the mRNA level of a large number of genes related to immune response and to sexual maturation was significantly modified in the silver eel swimbladder. The modification of several processes related to protein metabolism and transport, cell cycle, and apoptosis suggested that these changes in swimbladder metabolism and permeability were achieved by increasing cell turn-over. The impact of an infection of the swimbladder with the nematode Anguillicola crassus has been assessed by comparing these expression changes with expression changes observed between uninfected yellow eel swimbladder tissue and infected silver eel swimbladder tissue. In contrast to uninfected silver eel swimbladder tissue, in infected tissue the mRNA level of several glycolytic enzymes was significantly elevated, and with respect to extracellular matrix, several mucin genes were many-fold higher in their mRNA level. Modification of many immune related genes and of the functional categories "response to DNA damage stimulus" and "cellular response to stress" illustrated the damaging effect of the nematode infection. This study has identified a range of cellular processes in the swimbladder of silver eels that appear to be altered by nematode infection. These altered cellular processes could contribute to detrimental changes in swimbladder function that, in turn, may lead to impairment of spawning migration.

Complete sexual maturation of European eels (Anguilla anguilla) in captivity can only be achieved via injections with gonadotropins. For female eels this procedure takes 4-6months and the response ranges from "unresponsive" to final maturation and ovulation. Reproductive success could be significantly increased via early selection of responders based on predictive markers and minimally invasive sampling methods. To get a better understanding of the genetic background of ovarian maturation of the European eel we performed a pilot deep-sequencing transcriptome analysis of ovarian tissue derived from a yellow eel, a prepubertal silver eel and a post-spawning matured eel. Two key players in steroidogenesis were strongly correlated with advanced sexual maturation, namely P450c17 and liver receptor homolog-1, suggesting that blood plasma steroids might qualify as minimally invasive markers for early detection of responders. Since the predictive value of plasma sex steroid levels for final maturation of the European eel had not yet been carefully examined, we performed an extensive artificial maturation trial. Farmed silver eels were treated with pituitary extracts and sampled at multiple time intervals. Expression of steroidogenesis-related genes in ovarian tissue of responding and non-responding eels after four weekly injections with pituitary extract was compared using a custom-built microarray and RNAseq. Increased expression of 17β-hsd1 was strongly linked to sexual maturation. Blood plasma levels of sex steroids were measured using ELISAs. We show that a 2.5-fold increase in blood-plasma estradiol level after 4 weekly pituitary extract injections is a strong predictor of final sexual maturation of female European eel., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Reproductive homing migration of salmonids requires accurate interaction between the reception of external olfactory cues for navigation to the spawning grounds and the regulation of sexual maturation processes. This study aimed at providing insights into the hypothesized functional link between olfactory sensing of the spawning ground and final sexual maturation. We have therefore assessed the presence and expression levels of olfactory genes by RNA sequencing (RNAseq) of the olfactory rosettes in homing chum salmon Oncorhynchus keta Walbaum from the coastal sea to 75 km upstream the rivers at the pre-spawning ground. The progression of sexual maturation along the brain-pituitary-gonadal axis was assessed through determination of plasma steroid levels by time-resolved fluoroimmunoassays (TR-FIA), pituitary gonadotropin subunit expression and salmon gonadotropin-releasing hormone (sgnrh) expression in the brain by quantitative real-time PCR. RNAseq revealed the expression of 75 known and 27 unknown salmonid olfactory genes of which 13 genes were differentially expressed between fish from the pre-spawning area and from the coastal area, suggesting an important role of these genes in homing. A clear progression towards final maturation was characterised by higher plasma 17α,20β-dihydroxy-4-pregnen-3-one (DHP) levels, increased pituitary luteinizing hormone β subunit (lhβ) expression and sgnrh expression in the post brain, and lower plasma testosterone (T) and 17β-estradiol (E2) levels. Olfactomedins and ependymin are candidates among the differentially expressed genes that may connect olfactory reception to the expression of sgnrh to regulate final maturation.

Since its discovery in mammals as a key-hormone in reproduction and metabolism, leptin has been identified in an increasing number of tetrapods and teleosts. Tetrapods possess only one leptin gene, while most teleosts possess two leptin genes, as a result of the teleost third whole genome duplication event (3R). Leptin acts through a specific receptor (LEPR). In the European and Japanese eels, we identified two leptin genes, and for the first time in vertebrates, two LEPR genes. Synteny analyses indicated that eel LEPRa and LEPRb result from teleost 3R. LEPRb seems to have been lost in the teleost lineage shortly after the elopomorph divergence. Quantitative PCRs revealed a wide distribution of leptins and LEPRs in the European eel, including tissues involved in metabolism and reproduction. Noticeably, leptin1 was expressed in fat tissue, while leptin2 in the liver, reflecting subfunctionalization. Four-month fasting had no impact on the expression of leptins and LEPRs in control European eels. This might be related to the remarkable adaptation of silver eel metabolism to long-term fasting throughout the reproductive oceanic migration. In contrast, sexual maturation induced differential increases in the expression of leptins and LEPRs in the BPG-liver axis. Leptin2 was strikingly upregulated in the liver, the central organ of the reproductive metabolic challenge in teleosts. LEPRs were differentially regulated during sexual maturation, which may have contributed to the conservation of the duplicated LEPRs in this species. This suggests an ancient and positive role of the leptin system in the vertebrate reproductive function. This study brings new insights on the evolutionary history of the leptin system in vertebrates. Among extant vertebrates, the eel represents a unique case of duplicated leptins and leptin receptors as a result of 3R.

The translational value of zebrafish high-throughput screens can be improved when more knowledge is available on uptake characteristics of potential drugs. We investigated reference antibiotics and 15 preclinical compounds in a translational zebrafish-rodent screening system for tuberculosis. As a major advance, we have developed a new tool for testing drug uptake in the zebrafish model. This is important, because despite the many applications of assessing drug efficacy in zebrafish research, the current methods for measuring uptake using mass spectrometry do not take into account the possible adherence of drugs to the larval surface. Our approach combines nanoliter sampling from the yolk using a microneedle, followed by mass spectrometric analysis. To date, no single physicochemical property has been identified to accurately predict compound uptake; our method offers a great possibility to monitor how any novel compound behaves within the system. We have correlated the uptake data with high-throughput drug-screening data from Mycobacterium marinum-infected zebrafish larvae. As a result, we present an improved zebrafish larva drug-screening platform which offers new insights into drug efficacy and identifies potential false negatives and drugs that are effective in zebrafish and rodents. We demonstrate that this improved zebrafish drug-screening platform can complement conventional models of in vivo Mycobacterium tuberculosis-infected rodent assays. The detailed comparison of two vertebrate systems, fish and rodent, may give more predictive value for efficacy of drugs in humans., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

The European eel is a critically endangered species that cannot be reproduced in captivity yet. Artificial maturation of female European eels can be achieved via a laborious and expensive procedure, including weekly injections with pituitary extracts for up to 6 months. The success rate is highly variable and a minimally invasive method for early selection of responsive eels would prevent the unnecessary and lengthy treatment of non-responding individuals. Since sexual maturation of European eels is accompanied by morphological changes of the pectoral fin, we examined whether fin could be used to monitor the response to the hormone treatment. Farmed eels were subjected to weekly injections with pituitary extracts and representative groups were sampled at 0 and 14-18 weeks of hormone treatment. Responders and non-responders were identified based on the gonado-somatic index. Transcriptomes of pectoral fin samples obtained at the start and end of the trial were mapped using Illumina RNAseq. Responders showed 384 and non-responders only 54 differentially expressed genes. Highly stringent selection based on minimum expression levels and fold-changes and a manual re-annotation round yielded 23 up-regulated and 21 down-regulated maturation marker genes. The up-regulated markers belong to five categories: proteases, skin/mucus structural proteins, steroid hormone signaling, tyrosine/dopamine metabolism and lipid metabolism. The down-regulated markers are either blood markers or lectin-related genes. In conclusion, pectoral fin transcriptomes are a rich source of indicator markers for monitoring hormone induced sexual maturation of female European eels. In addition, these markers provide important new insight into several fundamental processes in eel biology., (Copyright © 2014 Elsevier Inc. All rights reserved.)