The complete mitochondrial genome and description of a new cryptic species of Benedenia Diesing, 1858 (Monogenea: Capsalidae), a major pathogen infecting the yellowtail kingfish Seriola lalandi Valenciennes in the South-East Pacific
Parasites & Vectors volume 12, Article number: 490 (2019)
The monogenean Benedenia seriolae parasitizes fishes belonging to the genus Seriola, represents a species complex, and causes substantial impact on fish welfare in aquaculture systems worldwide. This study reports, for the first time, the complete mitochondrial genome of B. humboldti n. sp., a new cryptic species from the South-East Pacific (SEP).
The mitogenome of B. humboldti n. sp. was assembled from short Illumina 150 bp pair-end reads. The phylogenetic position of B. humboldti n. sp. among other closely related congeneric and confamiliar capsalids was examined using mitochondrial protein-coding genes (PCGs). Morphology of B. humboldti n. sp. was examined based on fixed and stained specimens.
The AT-rich mitochondrial genome of B. humboldti is 13,455 bp in length and comprises 12 PCGs (atp8 was absent as in other monogenean genomes), 2 ribosomal RNA genes, and 22 transfer RNA genes. All protein-coding, ribosomal RNA, and transfer RNA genes are encoded on the H-strand. The gene order observed in the mitochondrial genome of B. humboldti n. sp. was identical to that of B. seriolae from Japan but different from that of B. seriolae from Australia. The genetic distance between B. humboldti n. sp. and B. seriolae from Japan was high. Minor but reliable differences in the shape of the penis were observed between Benedenia humboldti n. sp. and congeneric species.
Phylogenetic analyses based on PCGs in association with differences in the shape of the penis permitted us to conclude that the material from the South-East Pacific represents a new species of Benedenia infecting S. lalandi off the coast of Chile. The discovery of this parasite represents the first step to improving our understanding of infestation dynamics and to develop control strategies for this pathogen infecting the farmed yellowtail kingfish, Seriola lalandi, in the South-East Pacific.
Monogeneans are a clade of hermaphroditic ectoparasitic flatworms mostly restricted to the skin, fins or gills of marine and freshwater fishes . Monogeneans exhibit direct development and do not require an intermediate host to complete their life-cycle, in contrast to that reported for other parasitic flatworms (i.e. digeneans ). Monogenean infestations in farmed fish can and do become pathogenic and outbreaks often result in substantial health issues to the fish population in aquaculture systems worldwide [2, 3]. Some implications of heavy monogenean infestations include direct fish stock loss, depressed fish growth, poor fish health and welfare, reduced value of the market product, and costs associated with monitoring and treatment programs .
Among capsalid monogeneans (family Capsalidae Baird, 1853), Benedenia spp. attach to fish via a pair of anterior attachment organs and an opisthaptor which pierce the epidermis and penetrate the dermis of the host . The presence of large numbers of Benedenia spp. parasites during outbreaks causes considerable skin irritation to fish and results in the fish ‘rubbing’ themselves along the bottoms and sides of tanks/cages. Furthermore, Benedenia spp. cause skin injuries in fish that often lead to secondary infections by opportunistic pathogens such as bacteria and/or fungi . Unfortunately, genomic resources are limited in monogenean parasites and this poor knowledge is constraining our understanding of infection mechanisms, virulence and pharmacological resistance, among others, in this and other groups of ectoparasites (but there are exceptions [6, 7]).
Benedenia seriolae (Yamaguti, 1934) Meserve, 1938, is a particularly persistent problem and a major barrier to efficient finfish production and industry growth worldwide [4, 8]. Benedenia seriolae is a well-known parasite on the epidermis of the yellowtail Seriola quinqueradiata Temminck & Schlegel, and S. dumerili (Risso), cultured in Japan  and the kingfish S. lalandi Valenciennes, in Australia [3, 10], New Zealand , Mexico  and Chile . The wide distribution of this parasite might be a consequence of the pan-Pacific distribution of the host species or alternatively, might indicate the existence of a species complex . Using a barcoding approach, it has recently been demonstrated that Seriola lalandi in the South-East Pacific (SEP) is parasitized by an entity genetically different from B. seriolae in the South-West Pacific (SWP) . Importantly, no major morphological differences were observed among B. seriolae parasitizing S. lalandi from the SEP and SWP, S. quinqueradiata and S. hippos Günther. A single trait, however, i.e. the shape of the penis, appears to be dissimilar among B. seriolae from different host species and localities (SEP vs SWP) . Despite the ecological and aquaculture/fishery importance of B. seriolae, no genomic resources exist for this species that could improve our understanding of its life-cycle and its impact on the health of its host populations.
The aim of this study was to describe the complete mitochondrial genome of B. seriolae of Sepúlveda & González  from the SEP and compare it to previously published mitogenomes of B. seriolae from the SWP (Australia and Japan). Importantly, mitochondrial sequence comparison allowed for the description of a new pathogen species, Benedenia humboldti n. sp. that infects S. lalandi in the SEP. This paper describes the mitochondrial genome of B. humboldti n. sp. from the SEP focusing on codon usage profiles and nucleotide composition of protein-coding genes (PCGs). Additionally, the secondary structure of each identified tRNA gene is described and non-coding regions are examined in more detail. Selective constraints in PCGs, including those commonly used for population genetic inference, were explored.
Field collection and sequencing
A total of 4 individuals of Benedenia humboldti n. sp. (syn. Benedenia seriolae of Sepúlveda & González ) were collected with forceps from the skin of the yellowtail kingfish Seriola lalandi previously captured by artisanal fishermen in Antofagasta, Chile (23°37′S, 70°24′W). The specimens were immediately fixed in 99% ethanol within a 5 ml centrifuge tube and transported to AUSTRAL-Omics (Valdivia, Chile).
Total genomic DNA was extracted from whole individuals using a High Pure PCR Template Preparation Kit (Roche, Penzberg, Germany), following the manufacturer’s protocol. DNA concentration and purity were assessed using a Quant-iT™ PicoGreen® dsDNA Assay Kit (Thermo Fisher Scientific, Waltham, USA) on a DQ300 Hoefer Fluorometer (Hoefer Inc., Holliston, MA, USA). An Illumina Nextera XT DNA Sample Prep Kit (Illumina, San Diego, CA, USA) was used for whole genome library construction following the manufacturer’s protocol. Briefly, 1 µg of genomic DNA was randomly sheared via nebulization, DNA fragment ends were repaired, 3’ adenylated, and ligated to Illumina adapters. The resulting adapter-ligated libraries were PCR-amplified, Illumina indexes added, and pooled for multiplexed sequencing on an Illumina MiSeq sequencer (Illumina) using a pair-end 250 bp run format.
A total of 4,684,263 reads were generated and made available in FASTQ format by the company. The totality of these reads was used for the mitochondrial genome assembly of B. humboldti n. sp. from the SEP.
Mitochondrial genome assembly of Benedenia humboldti n. sp.
Contaminants, low quality sequences (Phred scores < 30), Illumina adapters and sequences with less than 50 bp were removed using the software Trimmomatic , leaving 3,380,163 pair-end high quality reads for final mitogenome assembly. The mitogenome was assembled de novo using the NOVOPlasty pipeline v.1.2.3 . NOVOPlasty uses a seed-and-extend algorithm that assembles organelle genomes from whole genome sequencing (WGS) data, starting from a related or distant single ‘seed’ sequence and an optional ‘bait’ reference mitochondrial genome . To test the reliability of the assembly, we ran NOVOPlasty using two strategies. First, we used a single fragment of the cox1 gene available in GenBank (KC633872) as a seed. Secondly, we used the complete mitochondrial genome of B. seriolae (HM222526) as a bait reference mitogenome in addition to the same partial cox1 seed. We chose to use the mitochondrial genome of B. seriolae from the SWP as a ‘bait’ reference because it is the closely related congeneric species with a published mitochondrial genome available on GenBank . The two runs used a kmer size of 49 following the developer’s suggestions .
Annotation and analysis of the Benedenia humboldti n. sp. mitochondrial genome
The newly assembled mitochondrial genome was first annotated in the MITOS web server (http://mitos.bioinf.uni-leipzig.de)  using the echinoderm/flatworm genetic code (Translation Table 9). Annotation curation and start + stop codons corrections were performed using MEGA6  and Expasy (https://web.expasy.org/). Genome visualization was conducted with OrganellarGenomeDRAW (http://ogdraw.mpimp-golm.mpg.de/index.shtml) . The open reading frames (ORFs) and codon usage profiles of PCGs were analyzed. Codon usage for each PCG was predicted using the invertebrate echinoderm/flatworm code in the Codon Usage web server (http://www.bioinformatics.org/sms2/codon_usage.html). tRNA genes were identified in the software MITFI  as implemented in the MITOS web server and the secondary structure of each tRNA was predicted using the tRNAscan-SE v.2.0 web server (http://trna.ucsc.edu/tRNAscan-SE/) . tRNA secondary structures were visualized in the RNAfold web server (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi) .
A relatively short non-coding region in the mitochondrial genome of B. humboldti n. sp. from the SEP was examined in more detail. The number of repeats in the region was investigated with the Tandem Repeat Finder v.4.09 web server (http://tandem.bu.edu/trf/trf.html) . We also attempted to discover DNA motifs in this intergenic region using the default options in MEME . Mfold (http://unafold.rna.albany.edu/) and Quickfold (http://unafold.rna.albany.edu/?q=DINAMelt/Quickfold) web servers were used to predict the secondary structure of this region with particular attention to the presence of stem-loops.
Selective constraints in PCGs were explored. Overall values of KA (the number of non-synonymous substitutions per non-synonymous site: KA = dN = SA/LA), KS (number of synonymous substitutions per synonymous site: KS = dS = SS/LS) and ω (the ratio KA/KS) were estimated for each PCG in the software KaKs_calculator v.2.0 . The above values were based on a pairwise comparison between B. humboldti n. sp. and B. seriolae from Australia (GenBank: HM222526). Next, to identify positively selected sites along the length of each examined sequence, the values of KA, Ks and ω were also calculated while adopting a sliding window (window length = 57, step length = 12) that ‘slipped’ along each sequence. The γ-MYN model  was used during calculations to account for variable mutation rates across sequence sites . If PCGs are under no selection, selective constraint (purifying selection) or diversifying selection, the ratio ω (= KA/KS) is expected to be equal to 1, < 1 or > 1, respectively .
The phylogenetic position of B. humboldti n. sp. and B. seriolae from the SWP (Australia and Japan, see below) among other species belonging to the subclass Monopisthocotylea of monogenetic flukes (class Monogenea) was examined. The newly assembled and annotated mitochondrial genome of B. humboldti n. sp., 12 sequences for B. seriolae from the SWP (available on GenBank), and those of a total of 23 other species of monopisthocotylean flukes retrieved from the GenBank database were used for the phylogenetic analysis conducted using the MitoPhAST pipeline . We used three species of monogeneans in the subclass Polyopisthocotylea as the outgroups for the analysis. MitoPhAST extracts all PCG nucleotide sequences from species available on GenBank and others provided by the user (i.e. B. humboldti n. sp. from the SEP), translates each PCG nucleotide sequence to amino acids, conducts alignments for each PCG amino acid sequence using Clustal Omega , removes poorly aligned regions with trimAl , partitions the dataset and selects best fitting models of sequence evolution for each PCG with ProtTest , and uses the concatenated and partitioned PCG amino acid alignments to perform a maximum likelihood phylogenetic analysis in the software RaxML . The robustness of the ML tree topology was assessed by bootstrap reiterations of the observed data 100 times.
Specimens of B. humboldti n. sp. were carefully removed from the skin of freshly sacrificed S. lalandi specimens. Nine specimens were fixed and stored in 70% ethanol. Fixed specimens were stained with Ehrlich’s haematoxylin for 15 min and then unstained in 1% HCl diluted in 70% ethanol. Next, each specimen was dehydrated in an ethanol series (70% × 10 min, 80% × 10 min, 90% × 10 min, 95% × 15 min and 100% × 15 min), cleared with xilene and mounted on slides in Canada balsam. Each specimen was examined and compared with other species of Benedenia based on morphological characteristics following criteria provided by the specialized literature [31,32,33,34]. The specimens were examined under an Olympus BX41 light microscope (Olympus, Tokyo, Japan) connected to a Micrometrics camera (590CU, ACCU-SCOPE Inc., Commack, NY, USA). Images were processed with Micrometric SE Premium software (ACCU-SCOPE Inc., Commack, NY, USA). Body measurements, including total body length and width, haptor length and width, hook length, testes and germarium length and width, pharynx length and width, and penis length are given in micrometers as the range followed by the mean and the number of specimens measured in parentheses. Additionally, 10 live specimens obtained from cultured S. lalandi during January 2019 were observed under a stereomicroscope (Olympus SZX7).
Results and discussion
The two strategies employed to assemble the mitochondrial genome of B. humboldti n. sp. from the SEP in NOVOPlasty resulted in identical sequences. The complete mitochondrial genome of B. humboldti n. sp. (GenBank: MK599467) is 13,455 bp in length and comprises 12 protein-coding genes (PCGs), 2 ribosomal RNA genes [rrnS (12S ribosomal RNA) and rrnL (16S ribosomal RNA)] and 22 transfer RNA (tRNA) genes. The PCG atp8 is lacking in the mitochondrial genome of B. humboldti n. sp., in agreement with that reported for the remaining monogeneans whose mitochondrial genomes have been assembled and annotated . All of the PCGs, tRNA genes and the 2 ribosomal RNA genes were encoded on the H-strand (Fig. 1, Table 1).
The gene order observed in B. humboldti n. sp. is identical to that reported in B. seriolae from Japan (unpublished sequences retrieved from GenBank) and the congeneric Benedenia hoshinai Ogawa, 1984 . In turn, gene order of B. humboldti n. sp. is different from that of B. seriolae from Australia . In B. seriolae from Australia, the trnT gene occurs between nad4 and trnF while the same gene is located between cox1 and rrnL in B. humboldti n. sp. from the SEP, B. seriolae from Japan and B. hoshinai (Fig. 2).
Gene overlaps comprising a total of 47 bp were observed in 4 gene junctions: atp6-nad2 (overlap of 1 bp), trnA(tgc)-trnD(gtc) (1 bp), trnP(tgg)-trnI (gat) (1 bp) and cox1-trnT(tgt) (40 bp) (Fig. 1, Table 1). In turn, short intergenic spaces ranging in size between 1 and 140 bp were observed in a total of 19 gene junctions. A single relatively long intergenic space involving 254 bp was assumed to be involved in replication initiation of the mitochondrial genome of B. humboldti n. sp. as it was found to contain similar features reported for the D-loop/Control Region of other invertebrates (Fig. 1).
Eleven out of the 12 PCGs in the mitochondrial genome of B. humboldti n. sp. exhibited the conventional flatworm/echinoderm mitochondrial start codon ATG (Table 1). nad2 exhibited the conventional start codon GTG, also observed in the congeneric B. seriolae from the SWP  and in Neobenedenia melleni (MacCallum, 1927) Yamaguti, 1963 . By contrast to B. humboldti n. sp. from the SEP and B. seriolae from the SWP, the congeneric B. hoshinai features the stop codon TAA . All PCGs ended with a complete and conventional termination codon, either TAG or TAA. No PCG terminated with an incomplete stop codon T, as often observed in other monogenean mitochondrial genomes [15, 35, 36].
The PCGs in the mitochondrial genome of B. humboldti n. sp. contained an A + T bias with an overall base composition of A = 25.7%, T = 46.6%, C = 11.1% and G = 16.6%. This A + T bias is within the known range reported for mitochondrial genomes in monogenetic flukes and other flatworms and likely affects codon usage. In the PCGs of B. humboldti n. sp., the most frequently used codons were UUU (Phe, n = 354 times used, 10.95% of the total), UUA (Leu, n = 311, 9.62%), AUU (Ile, n = 172, 5.32%) and UAU (Tyr, n = 153, 4.73%). Less frequently used codons included GCG (Ala, n = 1, 0.03%), CGC (Arg, n = 1, 0.03%), CCG (Pro, n = 2, 0.06%), UCG (Ser, n = 4, 0.12%) and ACG (Thr, n = 5, 0.16%) (Table 2).
The KA/KS ratios in all mitochondrial PCGs of B. humboldti n. sp. showed values < 1, indicating that all these PCGs are evolving under purifying selection. Examination of KA/KS ratio values in sliding windows across the length of each PCG sequence further indicated that purifying selection is acting along the entire PCG (Fig. 3, Additional file 1: Table S1). Interestingly, the overall KA/KS ratios observed for cox1, cox2 and cox3 (KA/KS < 0.00492, 0.00492 and 0.00502, respectively) were an order of magnitude lower than those observed for the remaining PCGs (range: 0.01454–0.07535) suggesting strong purifying selection and evolutionary constraints in the former genes (Fig. 3). Selective pressure in mitochondrial PCG has been poorly studied in monogenetic flukes but a similar pattern of widespread purifying selection in mitochondrial PCGs has been observed in other (marine) invertebrates, including flatworms .
tRNA genes encoded in the mitochondrial genome of B. humboldti n. sp. ranged in length from 58 to 61 bp and all but one [trnS1(gct)] exhibited a standard ‘cloverleaf’ secondary structure as predicted by MITFI. For the trnS1(gct) gene, MITFI predicted a secondary structure with a missing dihydrouridine arm, a feature also observed in the mitochondrial genomes of B. seriolae from the SWP and the closely related B. hoshinai and Neobenedenia melleni [15, 35, 36]. Interestingly, the RNAfold web server was not able to enforce the secondary structure of the trnH(gtg) gene predicted by MITFI resulting in the reconstruction of a tRNA with the pseudouridine stem forming a simple loop (Fig. 4). Additionally, the RNAfold web server was not able to enforce the secondary structure of the trnY(gta) gene predicted by MITFI resulting in the reconstruction of a tRNA with a missing pseudouridine arm. The anticodon nucleotides of all the tRNA genes are consistent with those found in closely related monogenean mitochondrial genomes . Benedenia hoshinai represents an exception as it exhibits the anticodon ACG instead of TCG in the trnR(tcg) gene .
The rrnL and rrnS genes identified in the mitochondrial genome of B. humboldti n. sp. were 951 and 750 nucleotides long, respectively. The rrnL gene is located between trnT(tgt) and trnC(gca). The rrnS gene is located close to the rrnL, between trnC(gca) and cox2 (Fig. 1). The two genes were A + T biased. The overall base composition of the rrnL gene was A = 31.7%, T = 42.8%, C = 10.0% and G = 15.5%, and that of the rrnS gene was A = 34.0%, T = 39.5%, C = 10.4% and G = 16.1%.
In B. humboldti n. sp., the relatively short 254 bp long non-coding region is located between the nad4l and nad4 genes (Fig. 1). The region was heavily A + T rich with an overall base composition of A = 34.3%, T = 52.4%, C = 7.5% and G = 5.9%. Visual examination of this non-coding region permitted the discovery of a single mononucleotide cytosine repeat near its 5’ end. The Tandem Repeat Finder web server analysis detected one 19-bp-long sequence (5’-TTA TAT ATT ATT TAA ATT T-3’) repeated twice (starting in position 134 and 182 from the 5’ end) in this region. The above features and several AT tandemly repeated sequences observed are in agreement to that observed in the non-coding region of the congenerics B. seriolae from the SWP  and B. hoshinai . Secondary structure prediction analysis in Mfold and Quickfold (assuming 27 °C) resulted each in three possible folding configurations, with a change in Gibbs free energy (ΔG) ranging from − 22.75 to − 23.44 kcal/mol (Additional file 2: Figure S1). In both Mfold and Quickfold, all three reconstructions featured stem-loop structures interspersed along the length of the region (Additional file 2: Figure S1). Nothing is known about replication of the mitochondrial genome in monogeneans . All the features present in the non-coding region of B. humboldti n. sp. have been observed before in the putative mitochondrial genome control region/D-loop of other invertebrates [37,38,39]. Thus, the observed mononucleotide cytosine repeats, high A + T rich base content, tandemly repeated AT sequences and predicted secondary structure(s) suggest an involvement of this non-coding region in the initiation of replication of the mitochondrial genome of B. humboldti n. sp.
The ML phylogenetic tree confirmed the monophyly of the subclass Monopisthocotylea within the class Monogenea and placed B. humboldti n. sp. from the SEP in a monophyletic clade together with B. seriolae from the SWP (Australia and Japan), B. hoshinai, and Neobenedenia melleni, in agreement with previous phylogenetic studies using a combination of partial mitochondrial and nuclear genes  (Fig. 5). Within this clade, B. humboldti n. sp. was sister to B. seriolae from Australia, a parasite of Seriola hippos. All B. seriolae from Japan clustered together into a well-supported monophyletic clade that was sister to the clade comprising B. humboldti n. sp. from the SEP and B. seriolae from Australia. Additional well-supported clades within the Monopisthocotylea included the families Dactylogiridae, Diplectanidae and Gyrodactylidae. Several nodes located near the root of the phylogenetic tree were well supported (Fig. 5). The above suggests that mitochondrial genomes alone will likely have enough phylogenetic information to reveal relationships at higher taxonomic levels within the subclass Monopisthocotylea.
The gene order herein reported for B. humboldti n. sp. is different than that of B. seriolae from Australia but identical to that reported for B. seriolae from Japan (see above). Importantly, the genetic distance between B. humboldti n. sp. and B. seriolae from the SWP (both Australia and Japan) was large (p-distance full mitogenome = 0.16; cox1 = 0.127; cytb = 0.134; rrnL = 0.096) and comparable to that previously calculated for other pairs of morphologically dissimilar species of Benedenia . Considering the information above, we examined the morphology of our specimens (from the SEP) in more detail and found minor but reliable differences when compared with B. seriolae from the SWP. In the following, we describe B. humboldti n. sp., a pathogen infecting S. lalandi off the coast of Chile.
Family Capsalidae Baird, 1853
Genus Benedenia Diesing, 1858
Benedenia humboldti Sepúlveda, González & Baeza, n. sp.
Syn. Benedenia seriolae of Sepúlveda & González 
Type-host: Seriola lalandi Valenciennes (Perciformes: Carangidae).
Type-locality: Off Antofagasta (23°37′S, 70°24′W), Chile.
Other localities: Off northern coast of Chile, from Antofagasta to Valparaíso (24°S, 70°W to 33°S, 71°W), and Archipelago of Juan Fernández (33°S, 80°W) .
Type-material: The holotype (stained whole mount) was deposited in the Chilean Museum of Natural History, Santiago, Chile, under the accession number MNHNCL PLAT-15005. Paratypes fixed in ethanol were deposited in the Chilean Museum of Natural History (3 lots: MNHNCL PLAT-15006, MNHNCL PLAT-15007 and MNHNCL PLAT-15008).
Site on host: Body surface.
Prevalence: 16% (29 out of 180 examined fish).
ZooBank registration: To comply with the regulations set out in article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new species have been submitted to ZooBank. The Life Science Identifier (LSID) of the article is urn:lsid:zoobank.org:pub:367FDE8C-75A7-4A09-A9B4-9E848F20E0F7. The LSID for the new name Benedenia humboldti is urn:lsid:zoobank.org:act:D4E5F88F-E1C5-445A-BF69-C3D4AE79CAC2.
Etymology: The specific epithet refers to Alexander von Humboldt.
[Based on 10 live specimens and nine flattened, preserved, stained and mounted adult specimens; Figs. 6, 7, 8, 9]. Total length including haptor 5526–11,210 (8147; n = 9); maximum body width at level of testes, 2553–5045 (3791; n = 9). Haptor slightly circular, with wider anterior portion, 1537–3289 (2232; n = 9) long, 1677–3421 (2376; n = 9) wide (Fig. 6). Accessory sclerites 2, located centrally on haptor, 305–654 (430; n = 5) long (Figs. 6, 7a). Anterior hamuli 2, elongated, strongly recurved distally, 374–705 (530; n = 5) long (Figs. 6, 7b); their proximal ends just overlap with proximal ends of accessory sclerites. Posterior hamuli 2, 83–118 (104; n = 4) long (Figs. 6, 7c). Hooklets 14, along haptor periphery. Tendons of extensive body musculature passing through proximal notch of accessory sclerites. Marginal valve present, substantially wider anteriorly (Fig. 6).
Anterior attachment organs 2, approximately circular or elliptical, 600–1237 (934; n = 9) long, 593–1070 (797; n = 9) wide. Pharynx 343–762 (589; n = 9) long, 412–952 (732; n = 9) wide. Eye-spots 2 pairs, dorsal, just anterior to pharynx. Gut caeca branched, not united posteriorly (Fig. 6). Testes 2, 777–1485 (1143; n = 9) long, 637–1380 (1050; n = 9) wide; each testis with variable numbers of columnar structures. Vas deferens widens to form small seminal vesicle at level just posterior to germarium, ascends along left side of germarium, and coils extensively before penetrating lateral wall of penis-sac about halfway along its length (Fig. 8); within penis, vas deferens follows longitudinal path towards distal tip of penis. Accessory gland reservoir prominent, occupies proximal quarter of penis-sac; duct long, joins vas deferens near distal tip of penis. Penis muscular, with size directly proportional to body size, 432–935 (694; n = 9) long, with proximal third broader than distal third, protrusible via common genital duct and submarginal, dorsolateral, common genital aperture (Figs. 6, 8, 9). Prominent dorsal rounded lobe present on left of common genital aperture (lo, Fig. 9). Long duct connects accessory gland reservoir to penis-sac. Glands of Goto not observed. Germarium globular, compact, 504–1100 (796; n = 9) long, 459–1110 (865; n = 9) wide, with relatively large internal fertilization chamber from which ovovitelline duct arises up to oӧtype (Fig. 8); columnar structures similar to those in testes present (s, Fig. 8). Vaginal opening on dorsal surface, posterior to common genital aperture. Vaginal opening leads to short straight duct, 83–136 (118; n = 3) long, narrowing into constricted region (cv, Fig. 8); vaginal duct travels posteriorly to enlarged proximal storage region of vagina communicating with vitelline reservoir. Oötype with thin-walled proximal region and bulbous thick-walled distal muscular region. Uterus opens into genital atrium at level of penis base. In live but not fixed specimens, connection between vitelline reservoir and oötype was detected during egg formation. Eggs tetrahedral (Fig. 7).
Of the 28 described species of Benedenia, B. humboldti n. sp. differs from B. beverleyburtonae Whittington & Deveney, 2011, B. acanthopagri (Hussey, 1986), B. anticavaginata Byrnes, 1986, B. lutjanis Whittington & Kearn, 1993 and B. ernsti Deveney & Whittington, 2010, because the latter five species possess a vaginal opening located anteriorly to the common genital pore [32, 34] or posteriorly to the left testis . In B. ovata (Goto, 1894), the vaginal pore opens at mid-body length, between the germarium and the common genital pore, and in B. sciaenae (Van Beneden, 1852) Odhner, 1905, male and female pores are separated but very closely located .
In contrast to the species listed above, in B. humboldti n. sp. the vaginal opening is located close to the left margin of the body and is posterior to the common genital pore like in most species of Benedenia. Additionally, B. humboldti n. sp. differs from B. rohdei Whittington, Kearn & Beverley-Burton, 1994, and B. jaliscana Bravo-Hollis, 1951, because the latter two species have the distal tip of penis armed with a sclerite [32, 41]. The specimens of B. humboldti n. sp. can be differentiated from other species of Benedenia described and/or reported from biogeographical regions other than the South Pacific by a combination of the following characters: body size; position of the median haptoral sclerites; size of haptor relative to body size; shape of the accessory sclerites and hamuli; relation between size of accessory sclerites and anterior hamuli [32,33,34].
Benedenia humboldti n. sp. most closely resembles B. seriolae from Seriola spp. and B. hendorffii (Linstow, 1889) Stiles & Hassall, 1908, from Coryphaena hippurus Linnaeus. Benedenia humboldti n. sp. and B. seriolae parasitize fishes belonging to the genus Seriola . The original description of B. seriolae  from S. aureovittata (= S. lalandi) was complemented [31, 43] with specimens obtained from S. quinqueradiata off Japan. Later, Whittington et al.  added morphological and morphometric information for B. seriolae from S. lalandi collected off Australia and Chile and suggested that B. seriolae was a cosmopolitan species infecting a variety of carangid fishes. Nonetheless, molecular analyses demonstrated that the species of Benedenia (identified as B. seriolae) from S. lalandi, S. quinqueradiata and S. hippos were genetically dissimilar; genetic distance was above 13% among the species but there was no significant morphometric disparity among them . The only morphological attribute that differentiates B. humboldti n. sp. (syn. B. seriolae of Sepúlveda & González ) from the SEP and B. seriolae from the SWP is penis shape. Benedenia humboldti n. sp. has a pine-nut elongated (lanceolated) penis shape while B. seriolae from the SWP has a blunt penis tip (Fig. 9).
Benedenia hendorffii was described by von Linstow  from the body surface of Coryphaena hippurus (L.) off Chile. No type-material was deposited by von Linstow  and Price  redescribed B. hendorffii based on a single specimen from an unknown fish species captured off Spokane, Washington, USA. Whittington et al.  checked the material by Price  and confirmed, based on this unique specimen, the identity of B. hendorffii. A comparison of B. humboldti n. sp. with the description and illustrations of B. hendorffii by von Linstow  revealed important differences between the two species such as the absence of a penis-sac (or a similar muscular organ), the existence of a separated uterine duct extending the length of the penis complex that opens separately and posteriorly to the male pore, and the absence of a vagina in B. hendorffii. Additionally, B. humboldti n. sp. differs from B. hendorffii redescribed by Price  by a combination of characters such as the accessory sclerites (striated in B. hendorffii) and the curvature of the distal end of the anterior hamulus (more open in B. humboldti n. sp. than in B. hendorffii). The penis shape of B. hendorffii looks similar to that of B. seriolae. We suggest that B. hendorffii should be considered a species inquirenda given the lack of type-material in the original description by von Linstow . In his description, von Linstow commented that the host specimens of C. hippurus were captured together with Seriola sp. hosts, which raises doubts about the correct identification of the host from which B. hendorffii specimens were taken. In addition, there is a lack of information about the host species from which the specimen redescribed as B. hendorffii by Price  was obtained. Finally, B. hendorffii has been found rarely on C. hippurus, and the presence of this monogenean in C. hippurus has been considered accidental .
In conclusion, this study assembled for the first time the mitochondrial genome of Benedenia humboldti n. sp., a cryptic species of great economic interest given its parasitic association with the yellowtail kingfish, Seriola lalandi, in aquaculture facilities from the SEP [8, 47, 48]. An integrative approach that included the study of the complete mitochondrial genome of Benedenia humboldti n. sp. from the SEP and B. seriolae from the SWP plus phylogenetic analyses and interrogation of morphological traits permitted us to confirm the existence of this new cryptic species in the genus Benedenia. The correct identity of this parasite represents the first step towards improving our understanding of infestation dynamics and control strategies of this pathogen in farmed S. lalandi in the SEP.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The mitochondrial genome sequence is available in the GenBank database under the accession number MK599467.
- KA :
number of nonsynonymous substitutions per nonsynonymous site
- KS :
number of synonymous substitutions per synonymous site
maximum-likelihood phylogenetic analysis
open reading frames
protein coding genes
- rrnS :
12S ribosomal RNA
- rrnL :
16S ribosomal RNA
Gibbs free energy
Rohde K. Marine parasitology. Collingwood: Csiro Publishing; 2005.
Whittington I, Cribb B, Hamwood T, Halliday J. Host-specificity of monogenean (Platyhelminth) parasites: a role for anterior adhesive areas? Int J Parasitol. 2000;30:305–20.
Ernst I, Whittington I, Corneillie S, Talbot C. Monogenean parasites in sea-cage aquaculture. Austasia Aquacult. 2002;16:46–8.
Hutson K, Ernst I, Whittington I. Risk assessment for metazoan parasites of yellowtail kingfish Seriola lalandi (Perciformes: Carangidae) in South Australian sea-cage aquaculture. Aquaculture. 2007;271:85–99.
Buchmann K, Bresciani J. Monogenea (phylum Platyhelminthes). In: Woo PT, editor. Fish diseases and disorders. Volume 1: Protozoan and metazoan infections. 2nd ed. Wallingford: CABI Publishing; 2006. p. 391–416.
Gallardo C, Valenzuela V, Núñez G. RNA-Seq analysis using de novo transcriptome assembly as a reference for the salmon louse Caligus rogercresseyi. PLoS ONE. 2014;9:e92239.
Núñez G, Valenzuela V, Pino J, Wadsworth S, Gallardo C. Insights into the olfactory system of the ectoparasite Caligus rogercresseyi: molecular characterization and gene transcription analysis of novel ionotropic receptors. Exp Parasitol. 2014;145:99–109.
Sepúlveda F, González MT. Molecular and morphological analyses reveal that the pathogen Benedenia seriolae (Monogenea: Capsalidae) is a complex species: implications for yellowtail Seriola spp. aquaculture. Aquaculture. 2014;418(419):94–100.
Ogawa K, Yokoyama H. Parasitic diseases of cultured marine fish in Japan. Fish Pathol. 1998;33:303–9.
Fensham J, Bubner E, Antignama T, Landos M, Caraguel C. Random and systematic sampling error when hooking fish to monitor skin fluke (Benedenia seriolae) and gill fluke (Zeuxapta seriolae) burden in Australian farmed yellowtail kingfish (Seriola lalandi). Prev Vet Med. 2018;153:7–14.
Sharp N, Poortenaar C, Diggles B, Willis T. Metazoan parasites of yellowtail kingfish, Seriola lalandi lalandi, in New Zealand: prevalence, intensity, and site preference. N Z J Mar Fresh. 2003;37:273–82.
Avilés A, Castelló F. Manual para el cultivo de Seriola lalandi (Pisces: Carangidae) en baja California sur, México. Mexico City: Instituto Nacional de La Pesca, Dirección General de Investigación en Acuacultura; 2004.
Bolger A, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2016;45:e18.
Perkins E, Donnellan S, Bertozzi T, Whittington I. Closing the mitochondrial circle on paraphyly of the Monogenea (Platyhelminthes) infers evolution in the diet of parasitic flatworms. Int J Parasitol. 2010;40:1237–45.
Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, et al. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phyl Evol. 2013;69:313–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW—a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41:W575–81.
Jühling F, Pütz J, Bernt M, Donath A, Middendorf M, Florentz C, Stadler PF. Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements. Nucleic Acids Res. 2012;40:2833–45.
Lowe T, Chan P. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016;44:W54–7.
Hofacker I. Vienna RNA secondary structure server. Nucleic Acids Res. 2003;31:3429–31.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.
Bailey T, Boden M, Buske F, Frith M, Grant C, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:W202–8.
Wang D, Zhang Y, Zhang Z, Zhu J, Yu J. KaKs_calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinform. 2010;8:77–80.
Wang D, Wan H, Zhang S, Yu J. γ-MYN: a new algorithm for estimating Ka and Ks with consideration of variable substitution rates. Biol Direct. 2009;4:20.
Tan M, Gan H, Schultz M, Austin C. MitoPhAST, a new automated mitogenomic phylogeny tool in the post-genomic era with a case study of 89 decapod mitogenomes including eight new freshwater crayfish mitogenomes. Mol Phylogenet Evol. 2015;85:180–8.
Sievers F, Wilm A, Dineen D, Gibson T, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.
Capella-Gutierrez S, Silla-Martinez J, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.
Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21:2104–5.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
Kearn G. Mating in the capsalid monogenean Benedenia seriolae, a skin parasite of the yellowtail, Seriola quinqueradiata, in Japan. Publ Seto Mar Biol Lab. 1992;35:273–80.
Whittington I, Deveney M, Wyborn SA. A revision of Benedenia Diesing, 1858 including a redescription of B. sciaenae (van Beneden, Odhner, 1905 and recognition of Menziesia Gibson, 1976 (Monogenea: Capsalidae). J Nat Hist. 1856;2001(35):663–777.
Deveney M, Whittington I. Three new species of Benedenia Diesing, 1858 from the Great Barrier Reef, Australia with a key to species of the genus. Zootaxa. 2010;2348:1–22.
Whittington I, Deveney M. New Benedenia species (Monogenea: Capsalidae) from Diagramma labiosum (Perciformes: Haemulidae) on the Great Barrier Reef, Australia, with oncomiracidial descriptions and a report of egg attachment to the host. J Parasitol. 2011;97:1026–34.
Kang S, Kim J, Lee J, Kim S, Min G, Park J. The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes). Mitochondrial DNA. 2012;23:176–8.
Zhang J, Wu X, Li M, Xie M, Li A. The complete mitochondrial genome of Neobenedenia melleni (Platyhelminthes: Monogenea): mitochondrial gene content, arrangement and composition compared with two Benedenia species. Mol Biol Rep. 2014;41:6583–9.
Zhang D, Hewitt G. Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997;25:99–120.
Kuhn K, Streit B, Schwenk K. Conservation of structural elements in the mitochondrial control region of Daphnia. Gene. 2008;420:107–12.
Li T, et al. A mitochondrial genome of Rhyparochromidae (Hemiptera: Heteroptera) and a comparative analysis of related mitochondrial genomes. Sci Rep. 2016;6:351375.
International Commission on Zoological Nomenclature. Amendment of articles 8, 9, 10, 21 and 78 of the International Code of Zoological Nomenclature to expand and refine methods of publication. ZooKeys. 2012;219:1–10.
Whittington I, Kearn G, Beverley-Burton M. Benedenia rohdei n. sp. (Monogenea: Capsalidae) from the gills of Lutjanus carponotatus (Perciformes: Lutjanidae) from the Great Barrier Reef, Queensland, Australia, with a description of the oncomiracidium. Syst Parasitol. 1994;28:5–13.
Yamaguti S. Studies on the helminth fauna of Japan. Part 2. Trematodes of fishes I. Jpn J Zool. 1934;5:249–541.
Hoshina T. On the monogenetic trematode, Benedenia seriolae, parasitic on yellow-tail, Seriola quinqueradiata. Bull Office Int Epizoot. 1968;69:1179–91.
Linstow O. Beitrag zur anatomie von Phylline hendorffii. Arch Mikrosk Anat. 1889;33:163–80.
Price E. The monogenetic trematodes of Latin America. In: Silva B, Travassos L, editors. Livro Jubilar do Professor Lauro Travassos. Rio de Janeiro, Brazil: Instituto Oswaldo Cruz Publishing; 1938. p. 407–13.
Williams E, Williams L. Checklists of the parasites of dolphin, Coryphaena hippurus, and Pompano Dolphin, C. equiselis with new records, corrections, and comments on the literature. Rev Fish Sci. 2010;18:73–93.
González MT, Sepúlveda F, López Z, Montenegro D, Irribarren P. Parásitos potencialmente patógenos en el cultivo de dorado (Seriola lalandi) en la macro-zona norte de Chile. 2nd ed. Antofagasta: Universidad de Antofagasta; 2013.
Bravo S, Hurtado C, Silva M. Coinfection of Caligus lalandei and Benedenia seriolae on the yellowtail kingfish Seriola lalandi farmed in a net cage in northern Chile. Lat Ame J Aquat Res. 2017;45:852–7.
Many thanks to Dr Vince P. Richards for bioinformatics support during the development of this project. MTG thanks Zambra Lopez for her assistance in the laboratory.
This study was funded by projects Fondecyt NO. 1130649 granted to MTG and MINEDUC-UA-Antofagasta no. ANT 1855 granted to MTG and JAB.
Ethics approval and consent to participate
Not applicable. Specimens of B. humboldti were carefully removed from the skin of freshly collected S. lalandi specimens available at a private fishermen facility. No permits are needed and extended to users for the collection of invertebrates and parasites from privately-owned fishes according to Chilean laws.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Selective pressure analysis in PCGs of B. humboldti n. sp. from the SEP. KA, KS and KA/KS values were calculated using the γ-MYN model and adopting a sliding window of length = 57 and step length = 12.
Secondary structure prediction analysis of the non-coding putative D-loop/CR in the mitochondrial genome of B. humboldti n. sp. from the SEP.
About this article
Cite this article
Baeza, J.A., Sepúlveda, F.A. & González, M.T. The complete mitochondrial genome and description of a new cryptic species of Benedenia Diesing, 1858 (Monogenea: Capsalidae), a major pathogen infecting the yellowtail kingfish Seriola lalandi Valenciennes in the South-East Pacific. Parasites Vectors 12, 490 (2019). https://doi.org/10.1186/s13071-019-3711-5
- Purifying selection