Toll-interacting protein is activated by the transcription factor GATA1 and Sp1 to negatively regulate NF-κB and MAPK pathways in the Japanese eel (Anguilla japonica).

Toll-interacting protein (Tollip) serves as a crucial inhibitory factor in the modulation of Toll-like receptor (TLR)-mediated innate immunological responses. The structure and function of Tollip have been well documented in mammals, yet the information in teleost remained limited. This work employed in vitro overexpression and RNA interference in vivo and in vitro to comprehensively examine the regulatory effects of AjTollip on NF-κB and MAPK signaling pathways. The levels of p65, c-Fos, c-Jun, IL-1, IL-6, and TNF-α were dramatically reduced following overexpression of AjTollip, whereas knocking down AjTollip in vivo and in vitro enhanced those genes' expression. Protein molecular docking simulations showed AjTollip interacts with AjTLR2, AjIRAK4a, and AjIRAK4b. A better understanding of the transcriptional regulation of AjTollip is crucial to elucidating the role of Tollip in fish antibacterial response. Herein, we cloned and characterized a 2.2 kb AjTollip gene promoter sequence. The transcription factors GATA1 and Sp1 were determined to be associated with the activation of AjTollip expression by using promoter truncation and targeted mutagenesis techniques. Collectively, our results indicate that AjTollip suppresses the NF-κB and MAPK signaling pathways, leading to the decreased expression of the downstream inflammatory factors, and GATA1 and Sp1 play a vital role in regulating AjTollip expression.

The expression profiles of cyp19a1, sf-1, esrs and gths in the brain-pituitary during gonadal sex differentiation in juvenile Japanese eels.

Eels are gonochoristic species whose gonadal differentiation initiates at the yellow eel stage and is influenced by environmental factors. We revealed some sex-related genes were sex dimorphically expressed in gonads during gonadal sex differentiation of Japanese eel (Anguilla japonica); however, the expression of sex-related genes in the brain-pituitary during gonadal sex differentiation in eels is still unclear. This study aimed to investigate the sex-related gene expressions in the brain-pituitary and tried to clarify their roles in the brain and gonads during gonadal sex differentiation. Based on our previous histological study, the control eels developed as males, and estradiol-17ß (E2) was used for feminization. Our results showed that during testicular differentiation, the brain cyp19a1 transcripts and aromatase proteins were increased significantly; moreover, the cyp19a1, sf-1, foxl2s, and esrs (except gperb) transcripts in the midbrain/pituitary also were increased significantly. Forebrain gnrh1 transcripts increased slightly during gonadal differentiation of both sexes, but the gnrhr1b and gnrhr2 transcripts in the midbrain/pituitary were stable during gonadal differentiation. The expression levels of gths and gh in the midbrain/pituitary were significantly increased during testicular differentiation and were much higher in males than in E2-feminized females. These results implied that endogenous estrogens might play essential roles in the brain/pituitary during testicular differentiation, sf-1, foxl2s, and esrs may have roles in cyp19a1 regulation in the midbrain/pituitary of Japanese eels. For the GnRH-GTH axis, gths, especially fshb, may be regulated by esrs and involved in regulating testicular differentiation and development in Japanese eels.

Effects of dopamine and melatonin treatment on the expression of the genes associated with artificially induced sexual maturation in Japanese eel, Anguilla japonica.

Japanese eel (Anguilla japonica) is a commercially important fish species in Asia. Understanding factors like photoperiod, temperature, and lunar cycles is crucial for successful aquaculture and managing its reproduction. Melatonin and dopamine (DA) are essential for regulating reproduction in vertebrates, including fish. This study investigated the effects of melatonin and DA on the reproductive system of mature male Japanese eels to better understand reproductive regulation in fish. To clarify the effects of these hormones on sexual maturation in eels, a critical stage in the reproductive process, sexual maturation was induced by injecting human chorionic gonadotropin, which stimulates the production of sex hormones. To check the effect of melatonin and DA on sexual maturation, DA, melatonin, and DA + domperidone were intraperitoneally injected into fish from each group (six per treatment) at a dose of 1 mg/kg body weight. The fish were then examined using quantitative RT-PCR by comparing the messenger RNA level of reproduction-related genes (gonadotropin releasing hormone 1; gnrh1, gonadotropin releasing hormone 2; gnrh2, follicle stimulating hormone; fshß, luteinizing hormone; lhß and DA receptor 2b; d2b), involved in the gonadotropic axis in eels, to those that received a control injection. The results indicate significant differences in the expression levels of gnrh1, gnrh2 and d2b in the brain and d2b, fshß, lhß in the pituitary at different stages of sexual maturation. Melatonin appears to enhance the production of sex gonadotropins, whereas DA inhibits them. These findings suggest an interaction between melatonin and DA in regulating reproduction in Japanese eels.

Transcriptome Analysis of Host Anti-Vibrio harveyi Infection Revealed the Pathogenicity of V. harveyi to American Eel (Anguilla rostrata).

Vibrio harveyi, a recently discovered pathogenic bacterium isolated from American eels (Anguilla rostrata), poses uncertainties regarding its pathogenesis in American eel and the molecular mechanisms underlying host defense against V. harveyi infection. This study aimed to determine the LD50 of V. harveyi in American eel and assess the bacterial load in the liver, spleen, and kidney post-infection with the LD50 dose. The results showed that the LD50 of V. harveyi via intraperitoneal injection in American eels over a 14d period was determined to be 1.24 × 103 cfu/g body weight (6.2 × 104 cfu/fish). The peak bacterial load occurred at 36 h post-infection (hpi) in all three organs examined. Histopathology analysis revealed hepatic vein congestion and thrombi, tubular vacuolar degeneration, and splenic bleeding. Moreover, quantitative reverse transcription polymerase chain reaction (qRT-PCR) results indicated significant up or downregulation of 18 host immune- or anti-infection-related genes post 12 to 60 hpi following the infection. Additionally, RNA sequencing (RNA-seq) unveiled 7 hub differentially expressed genes (DEGs) and 11 encoded proteins play crucial roles in the anti-V. harveyi response in American eels. This study firstly represents the comprehensive report on the pathogenicity of V. harveyi to American eels and RNA-seq of host's response to V. harveyi infection. These findings provide valuable insights into V. harveyi pathogenesis and the strategies employed by the host's immune system at the transcriptomic level to combat V. harveyi infection.

RNA-seq analysis revealed the pathogenicity of Vibrio vulnificus to American eel (Anguilla rostrata) and the strategy of host anti-V. vulnificus infection.

Vibrio vulnificus is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-V. vulnificus infection remains uncertain. In this study, American eels were infected with different dose of V. vulnificus to determine the LD50. Then, bacterial load in the liver and kidney histopathology were assessed post the LD50 of V. vulnificus infection. Additionally, gene expressions of 18 immune related genes in the liver, spleen and kidney were detected. Furthermore, transcriptome sequencing and enrichment of differentially expressed genes (DEGs) were analyzed in the eel spleens between pre-infection (Con_0), post-36 h (Vv_36), and post-60 h (Vv_60) infection. The results showed that LD50 of V. vulnificus to American eels was determined to be 5.0 × 105 cfu/g body weight, and the bacterial load peaked at 24 and 12 h post the infection (hpi) in the kidney and liver, respectively. The histopathology was highlighted by necrotic hepatocytes and splenic cells, congestion blood vessels in liver and spleen, atrophied glomeruli and vacuolization of renal tubular epithelial cells. The results of RT-PCR revealed that 18 host immune-related genes showed significantly up or downregulated expression post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 16 DEGs play essential role to the immunosuppression in American eels, and the protein-protein interactions shed light on the widespread upregulation GEGs related to metabolism and immune response maintained the host cell homeostasis post the V. vulnificus infection, shedding new light on our understanding of the V. vulnificus pathogenesis towards understudied American eel and the host anti-V. vulnificus infection strategies in gene transcript.

Aureispira anguillae sp. nov., isolated from Japanese eel Anguilla japonica leptocephali.

A novel filamentous eel-leptocephalus pathogenic marine bacterium, designated strain EL160426T, was isolated from Japanese eel, Anguilla japonica, leptocephali reared at a laboratory in Mie, Japan. In experimental infection studies on eel larvae, the strain EL160426T caused massive larval mortality and was reisolated from moribund leptocephali. Characteristically, observations of infected larvae found that EL160426T forms columnar colonies on the cranial surface of larvae. The novel isolate exhibited growth at 15-30 °C, pH 7-9, and seawater concentrations of 60-150% (W/V). Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain EL160426T was most closely related to Aureispira maritima 59SAT with 97.7% sequence similarity. The whole genome sequence analysis of the strain EL160426T showed that the strain maintained a circular chromosome with a size of approximately 7.58 Mbp and the DNA G + C content was 36.2%. The major respiratory quinone was MK-7 and the predominant cellular fatty acids were 16:0, 20:4 w6c (arachidonic acid), 17:0 iso and 16:0 N alcohol. DNA relatedness between the closest phylogenetic neighbor strain EL160426T and A. maritima (JCM23207T) was less than 13%. On the basis of the polyphasic taxonomic data, the strain represents a novel species of the genus Aureispira, for which the name Aureispira anguillae sp. nov. is proposed. The type strain is EL160426T (= JCM 35024 T = TSD-286 T).

Exploring bacterial community composition and immune gene expression of European eel larvae (Anguilla anguilla) in relation to first-feeding diets.

European eel (Anguilla anguilla) is a commercially important species for fisheries and aquaculture in Europe and the attempt to close the lifecycle in captivity is still at pioneering stage. The first feeding stage of this species is characterized by a critical period between 20 to 24 days post hatch (dph), which is associated with mortalities, indicating the point of no return. We hypothesized that this critical period might also be associated with larvae-bacterial interactions and the larval immune status. To test this, bacterial community composition and expression of immune and stress-related genes of hatchery-produced larvae were explored from the end of endogenous feeding (9 dph) until 28 dph, in response to three experimental first-feeding diets (Diet 1, Diet 2 and Diet 3). Changes in the water bacterial community composition were also followed. Results revealed that the larval stress/repair mechanism was activated during this critical period, marked by an upregulated expression of the hsp90 gene, independent of the diet fed. At the same time, a shift towards a potentially detrimental larval bacterial community was observed in all dietary groups. Here, a significant reduction in evenness of the larval bacterial community was observed, and several amplicon sequence variants belonging to potentially harmful bacterial genera were more abundant. This indicates that detrimental larvae-bacteria interactions were likely involved in the mortality observed. Beyond the critical period, the highest survival was registered for larvae fed Diet 3. Interestingly, genes encoding for pathogen recognition receptor TLR18 and complement component C1QC were upregulated in this group, potentially indicating a higher immunocompetency that facilitated a more successful handling of the harmful bacteria that dominated the bacterial community of larvae on 22 dph, ultimately leading to better survival, compared to the other two groups.

Analysis of Alternative Splicing and Long Noncoding RNAs After the Edwardsiella anguillarum Infected the Immunized European Eels (Anguilla anguilla) Revealed the Role of Outer Membrane Protein A in OmpA Subunit Vaccine.

Edwardsiella anguillarum is a bacterium that commonly infects cultivated eels. Outer membrane protein A (OmpA) emulsified with Freund's adjuvant has been shown to be an effective fishery vaccine against this pathogen. However, the specific roles of OmpA in the vaccine have not been fully explored. In this study, we performed RNA-seq in the liver of a European eel (Anguilla anguilla) after challenge with E. anguillarum in eels previously immunized with an OmpA subunit vaccine. Our aim was to elucidate the differentially alternative splicing (DAS) and differentially expressed long noncoding RNAs (DE-lncRNAs) using a genome-wide transcriptome. The results showed after that at 28 days post-immunization, eels challenged with E. anguillarum (Con_inf) exhibited severe pathological changes in the liver. In contrast, the OmpA infused eels (OmpA_inf group) showed infiltrated lymphocytes, while Freund's adjuvant-inoculated eels (FCIA_inf group) showed edema of hepatocytes and blood coagulation. The relative percent survival (RPS) was 77.7% and 44.4% for OmpA_inf and FCIA_inf compared to the Con_inf group. We identified 37 DE-lncRNAs and 293 DAS genes between OmpA_inf and FCIA_inf. Interactions between DAS gene-expressed proteins indicated that 66 expressed proteins formed 20 networks. Additionally, 33 DE-lncRNAs interacted with 194 target genes formed 246 and 41 networks in co-expression and co-location. Taken together, our findings demonstrate that the OmpA subunit vaccine elicits a higher RPS and provides novel insights into the role of OmpA through DAS genes and DE-lncRNAs perspective. These results are significant for the development of fishery subunit vaccines.

Tissue-specific promoters regulate the transcription of cyp19a1 in the brain-pituitary-gonad axis of Japanese eel (Anguilla japonica).

Aromatase is a key enzyme that catalyzes the biosynthesis of estrogens. Previous study indicated that putative tissue-specific promoters of the one aromatase gene (cyp19a1) may drive the differential regulatory mechanisms of cyp19a1 expression in Anguilla japonica. In the present study, for elucidating the transcription characteristics and the function of putative tissue-specific promoters of cyp19a1 in the brain-pituitary-gonad (BPG) axis during vitellogenesis, we investigated the transcriptional regulation of cyp19a1 by 17ß-estrogen (E2), testosterone (T), or human chorionic gonadotropin (HCG) in A. japonica. The expression of estrogen receptor (esra), androgen receptor (ara), or luteinizing hormone receptor (lhr) was up-regulated as cyp19a1 in response to E2, T, or HCG, respectively in the telencephalon, diencephalon, and pituitary. The expression of cyp19a1 was also upregulated in the ovary by HCG or T in a dose-dependent manner. Unlike in the brain and pituitary, the expression of esra and lhr, rather than ara, was upregulated by T in the ovary. Subsequently, four primary subtypes of 5'-untranslated terminal regions of cyp19a1 transcripts and the corresponding two 5' flanking regions (promoter P.I and P.II) were identified. The P.II existed in all BPG axis tissues, whereas the P.I with strong transcriptional activity was brain- and pituitary-specific. Furthermore, the transcriptional activity of promoters, the core promoter region, and the three putative hormone receptor response elements were validated. The transcriptional activity did not change when the HEK291T cells co-transfected with P.II and ar vector were exposed to T. These results suggested that the expression of cyp19a1 was upregulated indirectly through esra and lhr rather than ara by T in the ovary, whereas the expression of cyp19a1 was upregulated directly through androgen receptor and the downstream androgen response element of tissue-specific P.I in the brain and pituitary. The results of the study reveal the regulatory mechanisms of estrogen biosynthesis and provide a reference for optimizing the technology of artificially induced maturation in eels.

Validating the occurrence of the giant mottled eel (Anguilla marmorata) in the Galápagos Islands.

The giant mottled eel (Anguilla marmorata) is distributed mostly in the Indo-West Pacific. However, a few records indicate the presence of this eel in the Tropical Central and East Pacific. In April 2019, an eel specimen was caught in a small stream in San Cristobal Island, Galápagos. Morphological and molecular characters (16S and Cytb mtDNA sequences) confirmed the species as A. marmorata Quoy & Gaimard, 1824. The re-discovery of A. marmorata in Galápagos supports the hypothesis of an eastward range expansion from the west, probably through the North Equatorial Counter-Current.

Molecular identification of critically endangered European eels (Anguilla anguilla) in US retail outlets.

The European eel (Anguilla anguilla) has declined by over 90% since the early 1980s and has been listed as critically endangered. Yet, despite strict export bans from the European Union, the European eel is still sold illegally in many countries. Efforts to monitor the trade of European eels have been primarily concentrated in Asian markets where concerningly high rates of European eel have been reported. Comparably fewer studies have assessed the identities of eel samples from the United States (US), despite the obvious implications for eel conservation. To address this knowledge gap, we purchased 137 eel products (134 freshwater eels and three saltwater eels) from grocers, sushi bars, and restaurants in nine states across the US from 2019 to 2021. Seven samples (5.2%) labeled as freshwater eels (or "unagi") were identified as European eels using a combination of mitochondrial (cytochrome b) and nuclear (18S rRNA) restriction digestion assays, a fast and inexpensive molecular tool for seafood identification that can identify hybrids between European eels (A. anguilla) and American eels (A. rostrata). No hybrids between European and American eels were found and all seven samples identified with restriction digestion as European eels were confirmed by sequencing of cytochrome b and 18S rRNA. Frequency of European eels in US markets did not significantly correlate with state or retail type. Although illegal eel exports are likely reaching US consumers, the frequency of European eel samples in this study of the US market is much lower than found in other non-European countries.

Lipoprotein receptors in ovary of eel, Anguilla australis: molecular characterisation of putative vitellogenin receptors.

Lipoprotein receptors, including low-density lipoprotein receptor (LDLr) relatives (Lrs) and LDLr-related proteins (Lrps), belong to the LDLr supergene family and participate in diverse physiological functions. In this study, novel sequences of lr and lrp genes expressed in the ovary of the short-finned eel, Anguilla australis, during early gonadal development are presented. The genes encoding the LDLr-like, Lrp1-like, Lrp1b-like, Lrp3, Lrp4-like, Lrp5-like, Lrp6, Lrp10, Lrp11, Lrp12-like, and Lr11-like proteins were found and identified by sequence and structure analysis, in addition to phylogenetic analysis. Genes encoding proteins previously implicated in follicle development and vitellogenin (Vtg) uptake in oviparous vertebrates were also identified, i.e. lr8 (including lr8 + and lr8- variants) and lrp13; their identification was reinforced by conserved synteny with orthologues in other teleost fish. Compared to other lr/lrp genes, the genes encoding Lr8 + , Lr8-, and Lrp13 were highly expressed in ovary during early development, decreasing as oocyte development advanced when induced by hypophysation. Furthermore, lr8 + , lr8-, and lrp13 were dominantly expressed in the ovary when compared with 17 other tissues. Finally, this study successfully detected the expression of both lr8 variants, which showed different expression patterns to those reported in other oviparous vertebrates and provided the first characterisation of Lrp13 in Anguilla sp. We propose that lr8 + , lr8-, and lrp13 encode putative Vtg receptors in anguillid eels.

Effects of the short-term fasting and refeeding on growth-related genes in Japanese eel (Anguilla japonica) larvae.

The Japanese eel (Anguilla japonica) spends a long period as the leptocephalus larval form under current rearing conditions. The duration of the larval stage until metamorphosis is influenced by body size and growth; however, little knowledge exists of the regulatory mechanism of growth in eel larvae. The present study focused on growth hormone (GH), insulin-like growth factors (IGFs), and IGF binding protein (IGFBP) as the central regulators of growth in teleost fishes and transforming growth factor-beta 3 (TGF-ß3) as a possible key modulator of muscle growth and body component synthesis. Japanese eel IGFBP-1a and TGF-ß3, comprising 264 and 411 amino acids, respectively, were cloned. Short-term (5-day) fasting and refeeding experiments were performed to understand changes in growth-related genes affected by nutritional status. The relative expression of gh increased with fasting and subsequently decreased with refeeding to the basal levels of the fed control. Relative igf-1 and igf-2 expression levels were high in the fasted group. Relative igf-1 was reduced after 2-day refeeding, whereas igf-2 decreased to the basal level after 1-day refeeding, suggesting that IGF-1 and IGF-2 might be regulated independently and contribute to postnatal growth in eel larvae. Relative igfbp-1a expression was sharply increased by fasting, whereas tgf-ß3 showed high and low values in the fed and fasted groups, respectively. Relative igfbp-1a and tgf-ß3 levels were negatively and positively correlated with body size, respectively. These results suggest that igfbp-1a and tgf-ß3 are potential indices of growth for exploring optimal rearing conditions to shorten the larval stage in Japanese eels.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of European eel, Anguilla anguilla.

Eels are important aquaculture species for which an increasing number of reference genes are being identified and applied. In this study, five housekeeping genes [RPL7 (ribosomal protein L7), 18 S (18 S ribosomal RNA), EF1A (elongation factor 1α), ACTB (ß-actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase)] were chosen to evaluate their reliability as reference genes for quantitative real-time PCR (qPCR) for the study of Anguilla anguilla. The expression of the selected genes in different eel tissues was determined using qPCR at different growth stages or upon challenge by Anguillid herpesvirus (AngHV), and the expression levels of these genes were then compared and evaluated using the geNorm and NormFinder algorithms. Then, RefFinder was used to comprehensively rank the examined housekeeping genes. Interestingly, the expression of the evaluated housekeeping genes exhibited tissue-dependent and treatment-dependent variations. In different growth periods A. anguilla tissues, the most stable genes were the following: ACTB in mucus; 18 S in skin and kidney; RPL7 in muscle, gill, intestine and brain; EF1A in heart and liver; and GAPDH in spleen. In contrast, in AngHV-challenged A. anguilla tissues, the most stable genes were the following: 18 S in mucus; RPL7 in skin, gill, heart, spleen, kidney and intestine; EF1A in muscle and liver; and ACTB in brain. Further comparison analysis indicated that the expression of RPL7 and EF1A was stable in multiple A. anguilla tissues in different growth periods and in eels challenged by AngHV. Nonetheless, the expression level of GAPDH in eel tissues was lower, and it was unstable in several tissues. These results indicated that the selection of reference genes for qPCR analysis in A. anguilla should be made in accordance with experimental parameters, and both RPL7 and EF1A could be used as reference genes for qPCR study of A. anguilla at different growth stages or upon challenge by AngHV. The reference genes identified in this study could improve the accuracy of qPCR data and facilitate further studies aimed at understanding the biology of eels.

Blood-based gene expression as non-lethal tool for inferring salinity-habitat history of European eel (Anguilla anguilla).

The European eel is a facultative catadromous species, meaning that it can skip the freshwater phase or move between marine and freshwater habitats during its continental life stage. Otolith microchemistry, used to determine the habitat use of eel or its salinity history, requires the sacrifice of animals. In this context, blood-based gene expression may represent a non-lethal alternative. In this work, we tested the ability of blood transcriptional profiling to identify the different salinity-habitat histories of European eel. Eels collected from different locations in Norway were classified through otolith microchemistry as freshwater residents (FWR), seawater residents (SWR) or inter-habitat shifters (IHS). We detected 3451 differentially expressed genes from blood by comparing FWR and SWR groups, and then used that subset of genes in a machine learning approach (i.e., random forest) to the extended FWR, SWR, and IHS group. Random forest correctly classified 100% of FWR and SWR and 83% of the IHS using a minimum of 30 genes. The implementation of this non-lethal approach may replace otolith-based microchemistry analysis for the general assessment of life-history tactics in European eels. Overall, this approach is promising for the replacement or reduction of other lethal analyses in determining certain fish traits.

Immune function modulation during artificial ovarian maturation in Japanese eel (Anguilla japonica): A transcriptome profiling approach.

The Japanese eel (Anguilla japonica) experiences dramatic internal and external environmental changes during its transoceanic reproductive migrations. Here, we assess immune function changes in the primary and secondary immune organs (head kidney and spleen) of A. japonica during artificial ovarian maturation at the previtellogenic (PV), midvitellogenic (MV), and ovulating (OV) stages by transcriptome analyses. Stress responses were also assessed by determining the serum concentrations of lysozyme, alkaline phosphatase, acid phosphatase, total antioxidant capacity, and superoxide dismutase. Our results showed that together with increased serum 17ß-estrogen and testosterone, lysozyme activity and antioxidant capacity were suppressed during artificial ovarian maturation. Comparisons across these developmental stages identified 60 (head kidney) and 36 (spleen) differentially expressed genes associated with the immune system. Genes related to the key activation markers of innate immune function, such as CXCL10, CXCL11, CCL20, HSP90B, MMP9, and MMP13, were upregulated and significantly enriched in the interleukin-17 signaling pathway. Adaptive immune function-related genes (IGM and MHC1) were upregulated in the head kidney from PV to MV, and their levels increased thereafter in the spleen. Moreover, a correlation between Pax5 expression and IGM expression in the spleen of MV (IGM+/Pax5+) and OV (IGM++/Pax5-) stage suggests that adaptive immune function was enhanced during ovarian maturation. To our knowledge, the present study is the first to describe transcriptome profiling of immune organs during ovarian maturation in teleost. Our findings suggest that the interleukin-17 pathway and IgM may play important roles in spawning.

Preliminary insight into parental contributions to Japanese eel (Anguilla japonica) preleptocephali spawned on different nights.

Parentage sibship-inference analyses were conducted using mtDNA sequencing and six microsatellite genotypes of 182 Japanese eel preleptocephali that were collected from one net-tow near the West Mariana Ridge in May 2014. At least 328 parents were involved in producing the 182 preleptocephali, and several parents may have spawned a few times during 3 days of a spawning period. Half-sibs suggested that a few parents mated with 1-3 partners, indicating that the Japanese eel can form spawning aggregations in which several parents mate with each other in the ocean.

TAC3/TACR3 System Function in the Catadromous Migration Teleost, Anguilla japonica.

Neurokinin B (NKB), a member of the tachykinin (TAC) family, plays important roles in mammalian neuropeptide secretion in related to reproduction. However, its potential role in spawning migration teleost is less clear. In the present study, Japanese eel (Anguilla japonica) was employed to study the performance of NKB in regulating reproduction. Results showed that two tac3 and one tacr3 genes were identified in Japanese eel. Sequence analysis showed that two tac3 transcripts, tac3a and tac3b, encode four NKBs: NKBa-13, NKBa-10, NKBb-13, and NKBb-10. However, compared with other species, a mutation caused early termination of TACR3 protein was confirmed, leading to the loss of the 35 amino acid (aa) C-terminal of the receptor. Expression analysis in different tissues showed that both tac3a and tac3b mRNAs were highly expressed in the brain. In situ hybridization localized both tac3a and tac3b mRNAs to several brain regions, mainly in the telencephalon and hypothalamus. Because of the mutation in TACR3 of Japanese eel, we further analyzed whether it could activate the downstream signaling pathway. Luciferase assay results showed the negative regulation of cAMP Response Element (CRE) and Sterol Response Element (SRE) signal pathways by Japanese eel NKBs. Intraperitoneal injection of four different NKB mature peptides at 100 ng/g had negative effect on either gnrh or gth gene expression. However, the high concentration of NKBa-10 and NKBb-13 (1,000 ng/g) upregulated mgnrh and fshb or lhb expression level significantly, which may be mediated by other receptors. In general, the NKBs/NK3Rs system has important functions in regulating eel puberty onset.

Comprehensive Relationship Analysis of the Long Noncoding RNAs (lncRNAs) and the Target mRNAs in Response to the Infection of Edwardsiella anguillarum in European eel (Anguilla anguilla) Inoculated with Freund's Adjuvant.

Freund's complete adjuvant (FCA) and incomplete adjuvant (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the non-specific immune enhancement. As long noncoding RNAs (lncRNAs) play vital regulating roles in various biological activities, in this study, we examined the genome-wide expression of transcripts in the liver of European eel (Anguilla anguilla, Aa) inoculated with FCA and FIA (FCIA) to elucidate the regulators of lncRNAs in the process of Edwardsiella anguillarum (Ea) infection and Aa anti-Ea infection using strand-specific RNA-seq. After eels were challenged by Ea at 28 days post the first inoculation (dpi), compared to the control uninfected eels (Li group), the control infected eels (Con_Li group) showed severe bleeding, hepatocyte atrophy, and thrombi formed in the hepatic vessels of the liver, although eels inoculated with FCIA (FCIA_Li group) also formed slight thrombi in the hepatic vessels. Compared to the FCIA_Li group, there was about 10 times colony-forming unit (cfu) in the Con_Li group per 100 µg liver tissue, and the relative percent survival (RPS) of eels was 50% in FCIA_Li vs Con_Li. Using high-throughput transcriptomics, differential expressed genes (DEGs) and transcripts were identified and the results were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). Interactions between the differential expressed lncRNAs (DE-lncRNAs) and the target DEGs were explored using Cytoscape according to their co-expression and co-location relationship. We found 13,499 lncRNAs (10,176 annotated and 3423 novel lncRNAs) between 3 comparisons of Con_Li vs Li, FCIA_Li vs Li, and FCIA_Li vs Con_Li, of which 111, 110, and 129 DE-lncRNAs were ascertained. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs targeted by DE-lncRNAs revealed these DEGs mainly involved in single-organism cellular process in BP, membrane in CC and binding in MF, and KEGG pathways showed that the target DEGs in co-expression and co-location enriched in cell adhesion molecules. Finally, 118 DE-lncRNAs target 1161 DEGs were involved in an interaction network of 8474 co-expression and 333 co-location-related links, of which 16 DE-lncRNAs play vital roles in anti-Ea infection. Taken together, the interaction networks revealed that DE-lncRNAs underlies the process of Ea infection and Aa anti-Ea infection.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla).

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with Hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, probably reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggest interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight Hox clusters, and the results suggest important roles of methylation at Hox genes for adaptive processes.