Transcriptome analysis of blood for the discovery of sex-related genes in ricefield eel Monopterus albus

The blood acts as a transfer channel for a variety of factors in the whole body. The ricefield eel (Monopterus albus) is a protogynous hermaphrodite vertebrate. Until now, no research has reported an analysis of the blood transcriptome during the process of sexual development in the ricefield eel. In this study, the transcriptome sequencing of blood samples from male and female ricefield eels was completed with a total of 34.70 Gb clean data. The clean data of each sample all reached 5.23 GB, and the percent of the Q30 basic group was 88.62% and above. A total of 106,369 unigenes were obtained after assembly, including 13,296 unigenes with a length of more than 1 kb. Further functional annotation analysis showed that there are 28,522 unigenes that can be annotated. The annotations of genes with differential expression revealed that there were 563 genes with significant differential expression in the blood of male and female ricefield eels, including 91 upregulated genes and 472 downregulated genes. Among which, 14 genes may be closely related to sex differentiation, the qPCR was used to confirmed the expression pattern of those genes and result shown that 11 genes were downregulated and 3 genes were upregulated, consistent with the results of our RNA-Seq analysis. This blood transcript dataset will open future research avenues on ricefield eel sex development and differentiation.

Expression and regulation of Smad2 by gonadotropins in the protogynous hermaphroditic ricefield eel (Monopterus albus)

Smad2, a receptor-activated Smad, plays a critical role in regulating gametogenesis. In this study, a smad2 homologue was identified and sequenced from ricefield eel ovary cDNA, and its mRNA and protein expression levels were analysed during oocyte development. The cDNA sequence of ricefield eel smad2 consisted of 1863 bp encoding a 467-amino acid protein that had high sequence homology with Smad proteins in other teleosts, especially in Poeciliopsis prolifica. The results of real-time quantitative PCR (RT-qPCR) analysis revealed that smad2 is expressed in the ovary during gonad development, increased continuously until the early vitellogenic stage in the ovaries, and then decreased with ovary maturation. Smad2 protein immunoreactivity was localized in the cytoplasm of follicular cells, oogonia, and primary growth stage oocytes. In vitro experiments revealed that follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG) promoted smad2 expression in ovary tissue in a time- and dose-dependent manner, respectively. In summary, Smad2 plays a potentially vital role in ricefield eel ovary development.

Identification and Evaluation of Long Noncoding RNAs in Response to Handling Stress in Red Cusk-Eel (Genypterus chilensis) via RNA-seq

The red cusk-eel (Genypterus chilensis) is a native species with strong potential to support Chilean aquaculture diversification. Under commercial conditions, fish are exposed to several stressors. To date, little is known about the mechanism involved in the stress response of red cusk-eel, and there is no information related to the regulation mediated by long noncoding RNAs (lncRNAs). The objective of this work was to identify for the first time the lncRNAs in the transcriptome of G. chilensis and to evaluate the differential expression levels of lncRNAs in the liver, head kidney, and skeletal muscle in response to handling stress. We used previously published transcriptome data to identify the lncRNAs by applying a series of filters based on annotation information in several databases to discard coding sequences. We identified a total of 14,614 putative lncRNAs in the transcriptome of red cusk-eel, providing a useful lncRNA reference resource to be used in future studies. We evaluated their differential expression in response to handling stress in the liver, head kidney, and skeletal muscle, identifying 112, 323, and 108 differentially expressed lncRNAs, respectively. The results suggest that handling stress in red cusk-eel generate an altered metabolic status in liver, altered immune response in head kidney, and skeletal muscle atrophy through an important coding and noncoding gene network. This is the first study that identifies lncRNAs in Genypterus genus and that evaluates the relation between handling stress and lncRNAs in teleost fish, thereby providing valuable information regarding noncoding responses to stress in Genypterus species.

Molecular diet analysis of Anguilliformes leptocephalus larvae collected in the western North Pacific

Natural diets of leptocephalus larvae have been enigmatic. In this study, we collected DNA samples from the gut contents and body surface of leptocephali belonging to the five Anguilliform families (Anguillidae, Chlopsidae, Congridae, Muraenidae, and Serrivomeridae) from the northwest Pacific and performed next-generation 18S rDNA sequencing. Wide variety of eukaryotes was detected in both samples, from which eight eukaryotic groups (jellyfish, conoid parasite, tunicate, copepod, krill, segmented worm, fungi, and dinoflagellate) were selected on the basis of abundance. All groups except conoid parasites were common in both the samples. Cnidarian 18S rDNA reads were the most abundant in both the samples; however, the number of samples having cnidarian reads and the read counts were significantly higher in the body surface scraping samples than in the gut content samples, regardless of careful rinsing of the body surface. These results indicate that the cnidarian DNAs are most likely found because of cross contamination from the body surface and/or environment. 18S rDNA read counts of copepod and tunicate in the gut contents were greater than or comparable with those in the body surface scraping samples, which may correspond to the previous observations of fecal pellets and larvacean houses in the leptocephali gut. Thus, the present study supports previous implications that leptocephali utilize detritus materials, so called marine snow.

Resolving deep-sea pelagic saccopharyngiform eel mysteries: Identification of Neocyema and Monognathidae leptocephali and establishment of a new fish family "Neocyematidae" based on larvae, adults and mitogenomic gene orders

Deep-sea midwater "saccopharyngiform" eels of the families Cyematidae, Monognathidae, Eurypharyngidae and Saccopharyngidae (order Anguilliformes) are extraordinary fishes having major skeletal reductions and modifications compared to the general anguilliform body structure. Little is known about most aspects of the systematics, phylogeny, and ecology of these families, and few of the approximately 30 species described from adult specimens have been matched with their leptotocephalus larvae. Based on mitogenomic sequence data from rare new specimens, we show that the long-speculated-about larval form referred to as "Leptocephalus holti", which was thought to possibly be the larva of the rare orange-colored eels of Neocyema (5 known specimens; speculated to belong to the Cyematidae) are actually the larvae of the one-jaw eels of the family Monognathidae. One of the 5 types of L. holti larvae that were collected in the Pacific is genetically matched with Monognathus jesperseni, but multiple species exist based on larval sequence data and the morphology of adult specimens. A rare leptocephalus from the Sargasso Sea, with unique morphological characteristics including many small orange spots on the gut, was found to be the larva of Neocyema, which is presently only known from the Atlantic Ocean. We demonstrate that Neocyema constitutes a separate family being most closely related to Eurypharyngidae and Saccopharyngidae based on mitogenomic DNA sequences and unique mitochondrial gene orders.

Identification of genetic linkage group 1-linked sequences in Japanese eel (Anguilla japonica) by single chromosome sorting and sequencing

Japanese eel (Anguilla japonica) constitutes one of the most important food fish in Japan; accordingly, genome sequencing and linkage mapping have been conducted for the purpose of artificial cultivation. In the next stage, integration of genomic sequences within linkage groups (LG) is required to construct higher-resolution genetic markers for quantitative trait loci mapping and selective breeding of beneficial traits in farming. In order to identify LG1-linked scaffolds from the draft genome assembly (323,776 scaffolds) reported previously, we attempted to isolate chromosomes corresponding to LG1 by flow sorting and subsequent analyses. Initially, single chromosomes were randomly collected by chromosome sorting and subjected to whole-genome amplification (WGA). A total of 60 WGA samples were screened by PCR with primers for a known LG1-linked scaffold, and five positive WGA samples were sequenced by next-generation sequencing (NGS). Following reference mapping analysis of the NGS reads, four of the five WGA samples were found to be enriched by LG1-linked sequences. These samples were cytogenetically assigned to chromosome 5 by fluorescence in situ hybridization. Using blastn searches with 82,081 contigs constructed from the NGS reads of the four WGA samples as queries, 2323 scaffolds were identified as putative LG1-linked scaffolds from the draft genome assembly. The total length of the putative LG1-linked scaffolds was 99.0 Mb, comparable to the estimated DNA amounts of chromosome 5 (91.1 Mb). These results suggest that the methodology developed herein is applicable to isolate specific chromosome DNAs and integrate unanchored scaffold sequences onto a particular LG and chromosome even in teleost fishes, in which isolation of specific chromosomes by flow sorting is generally difficult owing to similar morphologies, sizes, and GC-contents among chromosomes in the genome. The putative LG1-linked scaffolds of Japanese eel contain a total of 6833 short tandem repeats which will be available for higher-resolution linkage mapping.

Expression of aryl hydrocarbon receptor-regulated genes and superoxide dismutase in the Antarctic eelpout Pachycara brachycephalum exposed to benzo[a]pyrene

The aryl hydrocarbon receptor (AhR) pathway mediates many, if not all, responses of fish to dioxin-like compounds. The Southern Ocean is progressively exposed to increasing concentrations of anthropogenic pollutants. Antarctic fish are known to accumulate those pollutants, yet nothing is known about their capability to induce chemical biotransformation via the AhR pathway. The objective of the present study was to investigate whether Antarctic eelpout, Pachycara brachycephalum, respond to anthropogenic pollutants by activation of the AhR and its target gene cytochrome P4501A (CYP1A), and of superoxide dismutase (SOD), which served as a representative for oxidative stress. We exposed P. brachycephalum to 10 and 100 mg benzo[a]pyrene (BaP)/kg body weight for 10 d and measured the expression of AhR, CYP1A, and SOD in liver tissue via quantitative polymerase chain reaction. We identified two distinct AhR isoforms in the liver of P. brachycephalum. Antarctic eelpout responded to both BaP exposures by an up-regulation of AhR and SOD, and by a particularly strong induction of CYP1A expression, which remained high until day 10 of the exposure time. Our data suggest that P. brachycephalum possesses the potential to up-regulate xenobiotic biotransformation pathways, at least at the gene expression level. The time course of the AhR and CYP1A response points to an efficient but slow xenobiotics metabolism. Moreover, BaP exposure could include adverse effects such as oxidative stress. Environ Toxicol Chem 2018;37:1487-1495. © 2018 SETAC.

Characterization, expression, and evolutionary analysis of new TLR3 and TLR5M genes cloned from the spiny eel Mastacembelus armatus

Toll-like receptors (TLRs) play an important role in innate and adaptive immunity. Here, we identify two new TLRs from the spiny eel Mastacembelus armatus (TLR3 and membrane TLR5M). Both MaTLR3 and MaTLR5M were expressed in all tested tissues; expression was highest in liver and spleen, respectively. After infection with Vibrio parahaemolyticus, expression of both TLRs fluctuated and differed significantly from controls at several time points. The predicted three-dimensional model of the MaTLR3 and MaTLR5M proteins indicates that most sites under positive selection were located in the extracellular domains of TLRs. Evolutionary analysis detected positively selected sites in the ancestral lineages of vertebrates, amphibians and reptiles. Multiple ML methods recovered 10 positively selected sites in teleost TLR3 and 24 in TLR5M, and most sites were located in leucine-rich repeat domain, possibly related to an "arms-race" co-evolution with pathogens.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads

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.

Ammonia stress under high environmental ammonia induces Hsp70 and Hsp90 in the mud eel, Monopterus cuchia

The obligatory air-breathing mud eel (Monopterus cuchia) is frequently being challenged with high environmental ammonia (HEA) exposure in its natural habitats. The present study investigated the possible induction of heat shock protein 70 and 90 (hsp70, hsc70, hsp90α and hsp90β) genes and more expression of Hsp70 and Hsp90 proteins under ammonia stress in different tissues of the mud eel after exposure to HEA (50 mM NH₄Cl) for 14 days. HEA resulted in significant accumulation of toxic ammonia in different body tissues and plasma, which was accompanied with the stimulation of oxidative stress in the mud eel as evidenced by more accumulation of malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) during exposure to HEA. Further, hyper-ammonia stress led to significant increase in the levels of mRNA transcripts for inducible hsp70 and hsp90α genes and also their translated proteins in different tissues probably as a consequence of induction of hsp70 and hsp90α genes in the mud eel. However, hyper-ammonia stress was neither associated with any significant alterations in the levels of mRNA transcripts for constitutive hsc70 and hsp90β genes nor their translated proteins in any of the tissues studied. More abundance of Hsp70 and Hsp90α proteins might be one of the strategies adopted by the mud eel to defend itself from the ammonia-induced cellular damages under ammonia stress. Further, this is the first report of ammonia-induced induction of hsp70 and hsp90α genes under hyper-ammonia stress in any freshwater air-breathing teleost.

Scp3 expression in relation to the ovarian differentiation in the protogynous hermaphroditic ricefield eel Monopterus albus

Synaptonemal complex protein 3 (Scp3), which is encoded by scp3, is a meiotic marker commonly used to trace the timing of gonadal differentiation in vertebrates. In the present study, the ricefield eel scp3 cDNA was cloned, and a fragment encoding amino acids 49 to 244 was overexpressed. The recombinant Scp3 polypeptide was purified and used to generate a rabbit anti-Scp3 polyclonal antiserum. In adult ricefield eels, scp3 mRNA was predominantly detected in the gonads and faintly detected in discrete brain areas. In the gonads, Scp3 immunoreactivities were shown to be localized to the germ cells, including meiotic primary growth oocytes, spermatocytes, and pre-meiotic spermatogonia. During early ovarian differentiation, immunoreactive Scp3 was not detected in the gonads of ricefield eels at 6 days post-hatching (dph) but was found to be abundantly localized in the cytoplasm of some oogonia at 7 dph, coinciding with the appearance of the ovarian cavity and ovarian differentiation. At 14 dph, strong Scp3 immunostaining was detected on one side of the nucleus with a distinct polarity in some germ cells, presumably at the leptotene stage. Consistent with these results, the expression of scp3 mRNA was faintly detected in the gonads of ricefield eels at 6 dph, increased at 8 dph, and then remained relatively high thereafter. Taken together, these results suggest that the appearance of immunoreactive Scp3 in oogonia could be a marker for early ovarian differentiation in ricefield eels. The translocation of the Scp3 protein from the cytoplasm to the nucleus in the oogonium of ricefield eels appears to be a controlled process that warrants further study.

Morphology and Species Composition of Southern Adriatic Sea Leptocephali Evaluated Using DNA Barcoding

Leptocephali are the characteristic larvae of the superorder Elopomorpha that are difficult to identify at the species level. In this study, we used DNA barcoding (i.e. short genetic sequences of DNA used as unique species tags) coupled with classical taxonomic methods to identify leptocephali in the southern Adriatic Sea. This information will provide an assessment of the biodiversity of the eel larvae in this region. A total of 2,785 leptocephali were collected, and using external morphology were assigned to seven morphotypes: Ariosoma balearicum, Conger conger, Gnathophis mystax, Facciolella sp., Nettastoma melanurum, Dalophis imberbis and Chlopsis bicolor. Collectively, these seven morphotypes are considered to be a good proxy for the Anguilliformes community (the main order of the Elopomorpha) in the southern Adriatic Sea (to date, seven families and sixteen species have been recorded in this region). Interestingly, the higher number of G. mystax larvae collected suggests an increased abundance of this genus. To validate the morphological identifications, we sequenced 61 leptocephali (at a 655 bp fragment from the cytochrome oxidase subunit 1 mitochondrial region) and developed barcode vouchers for the seven morphotypes. Using genetic information from reference databases, we validated three of these morphotypes. Where reference sequences were unavailable, we generated barcodes for both adult and juvenile forms to provide additional genetic information. Using this integrated approach allowed us to characterize a new species of Facciolella in the Adriatic Sea for the first time. Moreover, we also revealed a lack of differentiation, at the species level, between G. mistax and G. bathytopos, a western Atlantic Ocean species. Our morphological and barcode data have been published in the Barcoding of the Adriatic Leptocephali database. This work represents the first contribution to a wider project that aims to create a barcode database to support the assessment of leptocephali diversity in the Mediterranean Sea.

A giant anguilliform leptocephalus Thalassenchelys foliaceus Castle & Raju is a junior synonym of Congriscus maldivensis (Norman 1939)

A single specimen of giant leptocephalus Thalassenchelys foliaceus Castle & Raju 1975 was caught in subtropical waters of the western North Pacific Ocean. Mitochondrial coI gene sequence divergence between T. foliaceus and Congriscus maldivensis (Norman 1939) was 0·64 ± 0·27% (mean ± s.e.), and the myomere and vertebral counts of these species were similar, indicating T. foliaceus is a junior synonym of C. maldivensis.

Molecular phylogenetics reveals convergent evolution in lower Congo River spiny eels

BACKGROUND

The lower Congo River (LCR) is a region of exceptional species diversity and endemism in the Congo basin, including numerous species of spiny eels (genus Mastacembelus). Four of these exhibit distinctive phenotypes characterized by greatly reduced optic globes deeply embedded into the head (cryptophthalmia) and reduced (or absent) melanin pigmentation, among other characteristics. A strikingly similar cryptophthalmic phenotype is also found in members of a number of unrelated fish families, strongly suggesting the possibility of convergent evolution. However, little is known about the evolutionary processes that shaped diversification in LCR Mastacembelus, their biogeographic origins, or when colonization of the LCR occurred.

METHODS

We sequenced mitochondrial and nuclear genes from Mastacembelus species collected in the lower Congo River, and compared them with other African species and Asian representatives as outgroups. We analyzed the sequence data using Maximum Likelihood and Bayesian phylogenetic inference.

RESULTS

Bayesian and Maximum Likelihood phylogenetic analyses, and Bayesian coalescent methods for species tree reconstruction, reveal that endemic LCR spiny eels derive from two independent origins, clearly demonstrating convergent evolution of the cryptophthalmic phenotype. Mastacembelus crassus, M. aviceps, and M. simbi form a clade, allied to species found in southern, eastern and central Africa. Unexpectedly, M. brichardi and brachyrhinus fall within a clade otherwise endemic to Lake Tanganikya (LT) ca. 1500 km east of the LCR. Divergence dating suggests the ages of these two clades of LCR endemics differ markedly. The age of the crassus group is estimated at ~4 Myr while colonization of the LCR by the brichardi-brachyrhinus progenitor was considerably more recent, dated at ~0.5 Myr.

CONCLUSIONS

The phylogenetic framework of spiny eels presented here, the first to include LCR species, demonstrates that cryptophthalmia and associated traits evolved at least twice in Mastacembelus: once in M. brichardi and at least once in the M. crassus clade. Timing of diversification is broadly consistent with the onset of modern high-energy flow conditions in the LCR and with previous studies of endemic cichlids. The close genetic relationship between M. brichardi and M. brachyrhinus is particularly notable given the extreme difference in phenotype between these species, and additional work is needed to better understand the evolutionary history of diversification in this clade. The findings presented here demonstrate strong, multi-trait convergence in LCR spiny eels, suggesting that extreme selective pressures have shaped numerous phenotypic attributes of the endemic species of this region.

Comparative cytogenetics of six Indo-Pacific moray eels (Anguilliformes: Muraenidae) by chromosomal banding and fluorescence in situ hybridization

A comparative cytogenetic analysis, using both conventional staining techniques and fluorescence in situ hybridization, of six Indo-Pacific moray eels from three different genera (Gymnothorax fimbriatus, Gymnothorax flavimarginatus, Gymnothorax javanicus, Gymnothorax undulatus, Echidna nebulosa and Gymnomuraena zebra), was carried out to investigate the chromosomal differentiation in the family Muraenidae. Four species displayed a diploid chromosome number 2n = 42, which is common among the Muraenidae. Two other species, G. javanicus and G. flavimarginatus, were characterized by different chromosome numbers (2n = 40 and 2n = 36). For most species, a large amount of constitutive heterochromatin was detected in the chromosomes, with species-specific C-banding patterns that enabled pairing of the homologous chromosomes. In all species, the major ribosomal genes were localized in the guanine-cytosine-rich region of one chromosome pair, but in different chromosomal locations. The (TTAGGG)n telomeric sequences were mapped onto chromosomal ends in all muraenid species studied. The comparison of the results derived from this study with those available in the literature confirms a substantial conservation of the diploid chromosome number in the Muraenidae and supports the hypothesis that rearrangements have occurred that have diversified their karyotypes. Furthermore, the finding of two species with different diploid chromosome numbers suggests that additional chromosomal rearrangements, such as Robertsonian fusions, have occurred in the karyotype evolution of the Muraenidae.

An Enriched European Eel Transcriptome Sheds Light upon Host-Pathogen Interactions with Vibrio vulnificus

Infectious diseases are one of the principal bottlenecks for the European eel recovery. The aim of this study was to develop a new molecular tool to be used in host-pathogen interaction experiments in the eel. To this end, we first stimulated adult eels with different pathogen-associated molecular patterns (PAMPs), extracted RNA from the immune-related tissues and sequenced the transcriptome. We obtained more than 2x10⁶ reads that were assembled and annotated into 45,067 new descriptions with a notable representation of novel transcripts related with pathogen recognition, signal transduction and the immune response. Then, we designed a DNA-microarray that was used to analyze the early immune response against Vibrio vulnificus, a septicemic pathogen that uses the gills as the portal of entry into the blood, as well as the role of the main toxin of this species (RtxA1₃) on this early interaction. The gill transcriptomic profiles obtained after bath infecting eels with the wild type strain or with a mutant deficient in rtxA1₃ were analyzed and compared. Results demonstrate that eels react rapidly and locally against the pathogen and that this immune-response is rtxA1₃-dependent as transcripts related with cell destruction were highly up-regulated only in the gills from eels infected with the wild-type strain. Furthermore, significant differences in the immune response against the wild type and the mutant strain also suggest that host survival after V. vulnificus infection could depend on an efficient local phagocytic activity. Finally, we also found evidence of the presence of an interbranchial lymphoid tissue in European eel gills although further experiments will be necessary to identify such tissue.

Transcriptome profile analysis reveals specific signatures of pollutants in Atlantic eels

Identifying specific effects of contaminants in a multi-stress field context remain a challenge in ecotoxicology. In this context, "omics" technologies, by allowing the simultaneous measurement of numerous biological endpoints, could help unravel the in situ toxicity of contaminants. In this study, wild Atlantic eels were sampled in 8 sites presenting a broad contamination gradient in France and Canada. The global hepatic transcriptome of animals was determined by RNA-Seq. In parallel, the contamination level of fish to 8 metals and 25 organic pollutants was determined. Factor analysis for multiple testing was used to identify genes that are most likely to be related to a single factor. Among the variables analyzed, arsenic (As), cadmium (Cd), lindane (γ-HCH) and the hepato-somatic index (HSI) were found to be the main factors affecting eel's transcriptome. Genes associated with As exposure were involved in the mechanisms that have been described during As vasculotoxicity in mammals. Genes correlated with Cd were involved in cell cycle and energy metabolism. For γ-HCH, genes were involved in lipolysis and cell growth. Genes associated with HSI were involved in protein, lipid and iron metabolisms. Our study proposes specific gene signatures of pollutants and their impacts in fish exposed to multi-stress conditions.

Aromatase (Cyp19a1b) in the pituitary is dynamically involved in the upregulation of lhb but not fshb in the vitellogenic female ricefield eel Monopterus albus

Aromatase, encoded by Cyp19a1, is expressed in the pituitary of vertebrates; however, its physiological relevance remains poorly defined. In teleosts, the duplicated cyp19a1b is preferentially expressed in the pituitary where LH and FSH are synthesized in distinct gonadotropes. Our present study demonstrated that Cyp19a1b is colocalized with Lhb, but not Fshb, during vitellogenesis in female ricefield eels. The immunoreactive levels of Cyp19a1b and Lhb, as well as their colocalization frequency, increased during vitellogenesis toward maturation. The expression of lhb but not fshb in the pituitary fragments of female ricefield eels was induced by both estradiol (E2) and testosterone (T). In agreement, the promoter of lhb but not fshb was activated by both E2 and T. T is more potent than E2 in inducing lhb expression, whereas E2 is much more effective in activating the lhb promoter. T-induced lhb expression in the pituitary fragments was abolished by the estrogen receptor (Esr) antagonist fulvestrant and suppressed by the aromatase inhibitor letrozole, suggesting that the effect of T on lhb expression at the pituitary is largely mediated by E2. Furthermore, Lhb was shown to colocalize with Esr1 but not Esr2a. Taken together, results of the present study suggest that Cyp19a1b in LH cells may greatly upregulate lhb expression during vitellogenesis, possibly via E2 and Esr1 in an intracrine manner. The absence of Cyp19a1b in FSH cells and the insensitivity of fshb to sex steroids may contribute to the differential expression of lhb and fshb in ricefield eels and possibly other vertebrates as well.

Phylogeography of the Chinese beard eel, Cirrhimuraena chinensis Kaup, inferred from mitochondrial DNA: a range expansion after the last glacial maximum

The Chinese beard eel (Cirrhimuraena chinensis Kaup) is an intertidal fish and a model organism for the study of impacts caused by topological fluctuations during the Pleistocene and current intricate hydrological conditions on fauna living in the coastal areas of China. In this study, we examined the phylogeographical pattern, population genetic profile and demographical history of C. chinensis using mitochondrial DNA (cytochrome b (cyt b) and control region (CR)) from 266 individuals sampled in seven localities across the coastal area of southeastern China. The combined data indicated high levels of haplotype diversity and low levels of nucleotide diversity. Analyses of molecular variance (AMOVA) and F_ST statistics suggested the absence of a significant population structure across the Chinese coast. Neutrality tests, mismatch distributions and Bayesian skyline plots uniformly indicated a recent population expansion. The phylogeographic structure of C. chinensis may be attributed to past population expansion and long-distance pelagic larval dispersal facilitated by present-day ocean currents.

Characterization of thyroid hormone receptors during early development of the Japanese eel (Anguilla japonica)

We studied the profiles of thyroid hormone receptors (TRs) in Japanese eels (Anguilla japonica) during development from hatched larvae to juveniles. Two TRαs (TRαA and TRαB) and one TRβ (TRβA) cDNA clones were generated by RACE. The TRαA, TRαB and TRβA cDNAs encoded 416, 407 and 397 amino acid proteins with much higher homologies to the Japanese conger eel (Conger myriaster) TRs than to other fish TRs. In a transiently transfected Japanese eel cell line, Hepa-E1, the TRs showed thyroid hormone (TH)-dependent activation of transcription from the TH-responsive promoter. Four TR cDNA clones, including TRβB reported in a previous study, were analyzed by real-time RT-PCR. The TR mRNA levels in hatched larvae were determined. The two TRβ mRNAs were present at low levels but there was a peak in the TRαs during the larval stage before metamorphosis. During metamorphosis, the two TRαs both exhibited peaks and expression of the two TRβs was higher than during the early growth stage. This expression pattern is similar to that of the Japanese conger eel. It is possible that thyroid hormones control the early development of Japanese eels and Japanese conger eels through TRs. This is the first analysis of the expression sequence of TRs during early larval stages of Anguilliformes.

Intestinal expression of peptide transporter 1 (PEPT1) at different life stages of Japanese eel, Anguilla japonica

The expression of peptide transporter 1 (PEPT1) was investigated at the different life stages of Japanese eel, Anguilla japonica. The cDNA encoding Japanese eel PEPT1 was cloned and sequenced. The hydrophilicity plot analysis of its deduced amino acid sequence showed high similarities with topological features of known PEPT1 molecules in other species. Tissue distribution analysis confirmed that PEPT1 mRNA was detected specifically in the anterior and posterior intestines of adult eel. In eel larvae at 13days post hatching (dph), PEPT1 mRNA expression was mainly detected in the intestinal tract regions. The trypsinogen mRNA was only detected in the gastric region including the pancreas. Intense immunoreaction for PEPT1 was observed in the apical membrane of the intestinal epithelial cells of both larval and adult eel. These results indicated that PEPT1 was an intestine-specific transporter, which was localized at the luminal side of the epithelial cells, suggesting that di/tri-peptide absorption via PEPT1 takes place in the eel intestine. According to the ontogenetic analyses by quantitative PCR, PEPT1 and trypsinogen mRNA expressions were simultaneously increased at 5-7 dph. It is thus assumed that nutrient absorption systems in the intestinal tracts of larvae become functional at this age.

Epigenetic modifications during sex change repress gonadotropin stimulation of cyp19a1a in a teleost ricefield eel (Monopterus albus)

In vertebrates, cytochrome P450 aromatase, encoded by cyp19a1, converts androgens to estrogens and plays important roles in gonadal differentiation and development. The present study examines whether epigenetic mechanisms are involved in cyp19a1a expression and subsequent gonadal development in the hermaphroditic ricefield eel. The expression of the ricefield eel cyp19a1a was stimulated by gonadotropin via the cAMP pathway in the ovary but not the ovotestis or testis. The CpG within the cAMP response element (CRE) of the cyp19a1a promoter was hypermethylated in the ovotestis and testis compared with the ovary. The methylation levels of CpG sites around CRE in the distal region (region II) and around steroidogenic factor 1/adrenal 4 binding protein sites and TATA box in the proximal region (region I) were inversely correlated with cyp19a1a expression during the natural sex change from female to male. In vitro DNA methylation decreased the basal and forskolin-induced activities of cyp19a1a promoter. Chromatin immunoprecipitation assays indicated that histone 3 (Lys9) in both regions I and II of the cyp19a1a promoter were deacetylated and trimethylated in the testis, and in contrast to the ovary, phosphorylated CRE-binding protein failed to bind to these regions. Lastly, the DNA methylation inhibitor 5-aza-2'-deoxycytidine reversed the natural sex change of ricefield eels. These results suggested that epigenetic mechanisms involving DNA methylation and histone deacetylation and methylation may abrogate the stimulation of cyp19a1a by gonadotropins in a male-specific fashion. This may be a mechanism widely used to drive natural sex change in teleosts as well as gonadal differentiation in other vertebrates.

Experimental and computational study of inter- and intra- species specificity of gonadotropins for various gonadotropin receptors

The gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH) and their receptors play critical roles in vertebrate reproduction. In order to study intra- and interspecies ligand promiscuity of gonadotropins, COS-7 cells were transiently transfected with one of the gonadotropin receptor genes, FSHR or LHR, and tested for activation by gonadotropins from representative fish orders: Aquilliformes (eel; e), Salmoniformes (trout; tr), and Perciformes (tilapia; ta), and of mammalian origin: porcine (p), bovine (b) and human (h). The study reveals complex relations between the gonadotropin hormones and their receptors. Each gonadotropin activated its own cognate receptor. However, taLHR was also activated by hCG and eLHR was activated by hFSH, hCG, and trFSH. For FSHR, the only cross-reactivity detected was for hFSHR, which was activated by pFSH and bFSH. These findings are of great interest and applicability in the context of activation of various GTHRs by their ligands and by ligands from other vertebrates. Analysis of the three-dimensional models of the structures highlights the importance of residues outside of the currently established hormone-receptor interface region. In addition, the interface residues in taFSHR and the effect of exon duplication, which causes an insert in the LRR domain, are suggested to affect the interaction and binding of taFSH.

Genomics in eels--towards aquaculture and biology

Freshwater eels (genus Anguilla), especially the species inhabiting the temperate areas such as the European, American and Japanese eels, are important aquaculture species. Although artificial reproduction has been attempted since the 1930s and large numbers of studies have been conducted, it has not yet fully succeeded. Problems in eel artificial breeding are highly diverse, for instance, lack of basic information about reproduction in nature, no appropriate food for larvae, high mortality, and high individual variation in adults in response to maturation induction. Over the last decade, genomic data have been obtained for a variety of aquatic organisms. Recent technological advances in sequencing and computation now enable the accumulation of genomic information even for non-model species. The draft genome of the European eel Anguilla anguilla has been recently determined using Illumina technology and transcriptomic data based on next generation sequencing have been emerging. Extensive genomic information will facilitate many aspects of the artificial reproduction of eels. Here, we review the progress in genome-wide studies of eels, including additional analysis of the European eel genome data, and discuss future directions and implications of genomic data for aquaculture.

Primitive duplicate Hox clusters in the European eel's genome

The enigmatic life cycle and elongated body of the European eel (Anguilla anguilla L., 1758) have long motivated scientific enquiry. Recently, eel research has gained in urgency, as the population has dwindled to the point of critical endangerment. We have assembled a draft genome in order to facilitate advances in all provinces of eel biology. Here, we use the genome to investigate the eel's complement of the Hox developmental transcription factors. We show that unlike any other teleost fish, the eel retains fully populated, duplicate Hox clusters, which originated at the teleost-specific genome duplication. Using mRNA-sequencing and in situ hybridizations, we demonstrate that all copies are expressed in early embryos. Theories of vertebrate evolution predict that the retention of functional, duplicate Hox genes can give rise to additional developmental complexity, which is not immediately apparent in the adult. However, the key morphological innovation elsewhere in the eel's life history coincides with the evolutionary origin of its Hox repertoire.