Research Article
Print
Research Article
Morphological, molecular, and life cycle study of a new species of Oligogonotylus Watson, 1976 (Digenea, Cryptogonimidae) from Colombia
expand article infoVerónica Vélez-Sampedro§, Mónica Uruburu, Carolina Lenis
‡ Universidad de Antioquia, Medellín, Colombia
§ Universidad EAFIT, Medellín, Colombia
Open Access

Abstract

The present study describes Oligogonotylus andinus sp. nov. and its life cycle from a rural fish farm in Sopetrán, Antioquia, Colombia. The endemic species of snail Aroapyrgus colombiensis and the fishes Poecilia caucana and Andinoacara latifrons are identified as the first intermediate host, the second intermediate host and the definitive host, respectively. The new species was defined through an integrative approach, combining the traditional morphology of its developmental stages with molecular analyses of the markers ITS2 from ribosomal DNA and COI from mitochondrial DNA. This new species can be distinguished from its congeners by genetic divergence, the position of the vitelline fields, and the number of gonotyls. This work represents the first report of a species of this genus in South America.

Keywords

Blue mojarra, Cauca molly, Oligogonotylus, South America, taxonomy, Trematoda

Introduction

The genus Oligogonotylus Watson, 1976 (Digenea: Cryptogonimidae) includes two species, both occurring in Middle America (Choudhury et al. 2016, 2017). These infect the gastrointestinal tract of freshwater cichlid and characid fishes and have a distinctive morphological characteristic of a longitudinal row of gonotyls; these small sucker-like structures are present in a number of flukes in this family. The type species, Oligogonotylus manteri Watson, 1976, is distributed between southeastern Mexico and Costa Rica (Watson 1976; Pérez-Ponce de León et al. 2007; Sandlund et al. 2010). A second species, O. mayae Razo-Mendivil, Rosas-Valdez & Pérez-Ponce de León, 2008, is restricted to the Ría Celestun, Estero Progreso at Corchito, and the Ría Lagartos, Yucatán, Mexico (Razo-Mendivil et al. 2008). Both species utilize three hosts in their life cycle; the operculate micromollusk Pyrgophorus coronatus Pfeiffer, 1840 (Cochliopidae) as the first intermediate host (Ditrich et al. 1997; Jiménez-García and Vidal-Martínez 2005; May-Tec et al. 2020), several species of fishes as second intermediate hosts, and the native cichlid fish Mayaheros urophthalmus (Günther, 1862) as intermediate and the main definitive host among other 11 species of cichlids (Scholz et al. 1994, 1995; Pérez-Ponce de León et al. 2007; Sosa-Medina et al. 2015; May-Tec et al. 2020). In the present study, we describe a new species of Oligogonotylus and its life cycle from a rural fish farm in Colombia using an integrative approach combining ecological, molecular, and morphological data.

Methods

Sampling

We conducted fieldwork in La Miranda village, municipality of Sopetrán, Antioquia, Colombia (6°30'42"N, 75°45'22"W; Fig. 1) in March 2019 and February 2020. The sampling was performed in two fishponds supplied by La Mirandita brook in the Middle Cauca river basin. The ponds are approximately 5 m deep, 14 m long, 13 m wide for the first and 15 m wide for the second, with a total volume of 910 m3 and 1050 m3, respectively. Snails Aroapyrgus colombiensis Malek & Little, 1971 (n = 451), and fishes of the species Andinoacara latifrons (Steindachner, 1878) (n = 22), Colossoma macropomum (Cuvier, 1816) (n = 6), Oreochromis aureus (Steindachner, 1864) (n = 6), O. mossambicus (Peters, 1852) (n = 9), and Poecilia caucana (Steindachner, 1880) (n = 105) were collected from the ponds. All animals were transported live to the Laboratory of Helminthology, Programa de Estudio y Control de Enfermedades Tropicales, Facultad de Medicina, Universidad de Antioquia (Medellín, Colombia) and were placed in aquaria prepared with dechlorinated water, artificial aeration, an approximate temperature of 22–24 °C and a controlled photoperiod of 12 h of light and 12 h of darkness.

Figure 1. 

Type locality of the new species of Oligogonotylus from La Miranda village, Sopetrán, Antioquia, Colombia. Map made in ArcGIS 18.1.

Collection of the cercariae, metacercariae and adults

The search for the parasite specimens was based on standard techniques, cercarial emission in snails (Vergara and Velásquez 2009; Anucherngchai et al. 2016) and organ dissection of fishes to search for parasite infections under a stereo-dissecting microscope. To induce the cercarial emission, the snails were separated and examined individually under a microscope in 1.5 mL Eppendorf tubes with 0.1 mL of distilled water and stimulated under a cold light source for 4 h. To determine the prevalence of infection, non-emitting snails were analysed in triplicate for 21 d. The collected fishes were sacrificed with an overdose of isoeugenol according to Erikson (2011). The muscle of P. caucana was examined for the presence of metacercariae, and the digestive tract of A. latifrons, C. macropomum, O. aureus, O. mossambicus, and P. caucana were examined for the presence of juvenile and adult parasites. The trematodes were rinsed in saline solution, counted, and preserved in 1% glutaraldehyde (cercariae), alcohol-formalin-acetic acid (AFA), and 70% ethanol for morphological identification, or fixed directly in 96% molecular-grade ethanol. The ecological parameters of infections (prevalence, mean intensity and abundance of infection, and intensity range) were calculated according to Bush et al. (1997).

Morphological analyses

Cercariae were studied as temporary mounts stained with methylene blue. Metacercariae and adults were dehydrated in an alcohol series, stained with Meyer’s carmine or Borax carmine, cleared in methyl salicylate, and mounted on permanent slides in Canada balsam. A Nikon Alphaphot YS-2 with a Nikon 1.25× drawing tube was used to illustrate the parasites. A Leica DM500 (ICC50 W) microscope was used to analyse and photograph the stained specimens. Morphometric data are given as the mean and range in micrometers (µm). Line drawings and photographic images were prepared using Inkscape 0.92. Type specimens were deposited in the Colección Colombiana de Helmintos, CCH.116, Universidad de Antioquia, Medellín, Colombia.

DNA extraction, PCR, and sequencing

The DNA was extracted from individual specimens following the protocol of the genELUTE Mammalian genomic DNA miniprep kit. For barcoding, the nuclear ribosomal second internal transcribed spacer region (ITS2) was amplified using the primer pair: ITS-F 5'-CGGTGGATCACTCGGCTCGT-3' and ITSR 5'-CCTGGTTAGTTTCTTTTCCTCCGC-3' (Iwagami et al. 2000), and the partial cytochrome c oxidase subunit 1 mitochondrial gene (COI) was amplified using the primer pair: JB3 5'-TTTTTTGGGCATCCTGAGGTTTAT-3', and JB4 5'-TAAAGAAAGAACATAATGAAAATG-3' (Bowles et al. 1993). All amplification reactions were as follows: 5 min at 95 °C, followed by 35 cycles each of denaturing (94 °C, 30 s), annealing (50 °C, 45 s), and extension (72 °C, 1 min) for ITS2; annealing (53 °C, 45 s), and extension (72 °C, 1 min) for COI, and 5 min at 95 °C, followed by 35 cycles each for denaturing (94 °C, 30 s). The final extension was done at 72 °C for 5 min. PCR products were purified and sequenced by the Sanger method at Macrogen Inc., Seoul, Korea. The accession numbers in GenBank are given for each barcoded specimen.

Phylogenetic tree construction

The ITS2 and partial COI sequences of individual Oligogonotylus (ITS2: n = 4 adults, n = 1 metacercaria, n = 1 cercaria; COI: n = 6 adults) were edited and assembled in GENEIOUS 2022.0.1 (http://www.geneious.com; Kearse et al. 2012). All sequences were aligned with MUSCLE (Edgar 2004) including sequences of ITS2 and COI of other Oligogonotylus species and representative species of other genera of opisthorchids available in GenBank. The best-fitted nucleotide substitution models based on the Bayesian information criterion were estimated in JMODELTEST (Darriba et al. 2012). The genetic divergence (p-value) was calculated for each data set using MEGA X (Kumar et al. 2018) and the resulting alignments were analysed using IQ-TREE 2.1.3 (Nguyen et al. 2015) where all phylogenetic trees were constructed using the maximum-likelihood method. The ITS2 phylogenetic tree was based on the symmetrical model with gamma distribution (SYM+G); the COI phylogenetic tree was based on the general time-reversible model with gamma distribution (GTR+G); and COIITS2 concatenated tree was based on the general time-reversible model with proportion of invariable sites and gamma distribution (GTR+I+G). Downloaded sequences of some species of the same family (Cryptogonimidae) and another close family (Heterophyidae) were used in all cases as outgroup taxa to ensure the phylogenetic position of the new species. Trees were drawn to scale, and branch-support values were obtained with the ultrafast bootstrap approximation (UFBoot) (Hoang et al. 2018) implemented in IQ-TREE 2.1.3 (Nguyen et al. 2015) based on 1000 bootstrap replicates.

Results

We report a new species of Oligogonotylus from Colombia based on new morphological, molecular, and ecological evidence from a total of 1,306 worms. This taxon has a three-host life cycle, with the snail A. colombiensis as the first intermediate host, P. caucana as the second intermediate host, and A. latifrons as the definitive host (Fig. 2).

Figure 2. 

Micrographs of different stages of Oligogonotylus andinus sp. nov. A, B cercariae emerged from Aroapyrgus colombiensis, as temporary mounts C encysted metacercaria from Poecilia caucana, as a temporary mount D unencysted metacercaria isolated from P. caucana, as a permanent slide E, F paratype specimens, as temporary mounts G holotype adult specimen, isolated from A. latifrons, as a permanent slide. Scale bars: 20 μm (A, B); 20 μm (C, D); 100 µm (E–G).

Taxonomy

Family Cryptogonimidae Ward, 1917

Genus Oligogonotylus Watson, 1976

Oligogonotylus andinus sp. nov.

Diagnosis

Adult with longitudinal row of three or four sucker-like gonotyls. Vitelline follicles extending between level of pharynx and anterior margin of ovary and occupying 25–38% of body length.

Cercaria description

(Figs 2A, 3A). Pleurolophocercous type. Cercaria elongate, 149 (114–190) long by 82 (72–100) wide. Oral sucker 35 (18–46) long by 30 (24–34) wide, with three or four anterior rows of hook-like spines. Eyespots anterior to penetration glands. Pharynx located between eyespots. Esophagus and ceca undifferentiated. Penetration glands in middle body; seven pairs lead to exterior in four ducts running parallel to oral sucker in two bundles on each side. Rudimentary ventral sucker in middle body. Cystogenic glands on lateral margins of body. Excretory vesicle bilobed delimited anteriorly by primordial ventral sucker and posteriorly by insertion of tail. Tail longer than body with lateral fin, 294 (274–316) long by 22 (18–26) wide.

Figure 3. 

Developmental stages of Oligogonotylus andinus sp. nov. A pleurolophocercous cercaria emerged from Aroapyrgus colombiensis B encysted metacercaria isolated from the muscle of Poecilia caucana C unencysted metacercaria isolated from P. caucana D holotype, adult isolated from the large intestine of the trans-Andean cichlid Andinoacara latifrons. Scale bars: 20 μm (A); 20 μm (B, C); 100 µm (D).

Metacercaria description

(Figs 2B, 3B, C). Encysted form spherical with thin cyst wall of variable thickness occupying entire cyst cavity (Fig. 3B). Cyst 281 (247–327) long by 193 (136–223) wide. Unencysted forms covered with numerous, simple, single-pointed tegumental spines which are smaller and less dense posteriorly. Oral sucker (Os) large, subterminal, 79 (58–98) long by 67 (50–84) wide. Ventral sucker (Vs) just pre-equatorial to equatorial, 53 (34–70) long by 49 (34–66) wide, smaller than oral sucker; Os/Vs (ratio) 1:0.7 (1:0.6–1:0.8) long by 1:0.7 (1:0.6–1:0.9) wide. Pharynx large, 43 (36–48) long by 38 (30–46) wide. Esophagus short; ceca wide, reaching middle hindbody. Remnants of eyespots at pharyngeal level. Genital primordium indistinct, posterior to ventral sucker. Excretory vesicle large, Y-shaped, with large lateral branches reaching anteriorly to level of ventral sucker and with short posterior duct opening at excretory pore.

Adult description

(Figs 3D, 4A–F). Body oval to elongate-oval, widest at ovarian level, 517 (271–821, n = 21) long by 248 (167–335, n = 21) wide. Anterior and posterior ends of body rounded. Tegument spinose, with spines becoming shorter and less dense in posterior region. Remnants of eyespots diffuse and scattered at pharyngeal level. Oral sucker rounded, subterminal, without circumoral crown of spines, 106 (84–128, n = 21) long by 110 (66–144, n = 21) wide. Prepharynx short or inconspicuous, 18 (12–28, n = 9) long by 24 (16–28, n = 6) wide, with small extramural glands. Pharynx large, 59 (30–86, n = 18) long by 50 (26–70, n = 18) wide. Esophagus 40 (10–70, n = 13) long. Intestinal bifurcation immediately anterior to ventral sucker. Ceca reach just into post-testicular zone; distance from cecal extremities to posterior end of body 86 (60–114, n = 10). Longitudinal row of 4 (3–4, n = 17) sucker-like gonotyls located between mid-level of pharynx and anterior margin of ventral sucker, increasing in size from anterior to posterior. Ventral sucker rounded, pre-equatorial, 79 (56–100, n = 21) long by 78 (52–104, n = 21) wide, enclosed in ventrogenital sac, smaller than oral sucker; Os/Vs (ratio) 1:0.74 (1:0.61–1:0.88, n = 21) long by 1:0.72 (1:0.60–1:0.83, n = 21) wide; OsVs distance 112 (40–178, n = 15). Genital pore close to anterior margin of ventral sucker. Testes two, rounded to oval, in tandem, intercecal, in mid hindbody: anterior testis 72 (20–120, n = 20) long by 82 (50–112, n = 20) wide; posterior testis 80 (50–140, n = 21) long by 81 (46–140, n = 21) wide. Seminal vesicle tubulosaccular, undivided, mostly posterolateral to ventral sucker, 95 (30–128, n = 10) long by 82 (48–120, n = 11) wide. Cirrus sac and cirrus absent. Ovary lobed, median, post-equatorial, intercecal, between ventral sucker and posterior testis, partially overlaps anterior testis, 74 (38–108, n = 17) long by 90 (30–120, n = 17) wide. Seminal receptacle median, ovoid, dorsal, anterolateral to ovary, and ventral to vitelline duct, 65 (40–94, n = 11) long by 59 (40–82, n = 11) wide. Ootype posterolateral to ovary (on right side in type specimen and three paratypes, on left side in three other paratypes). Vitelline follicles small, irregular in shape to oval; arranged in lateral fields partly overlapping ceca or extracecal. Vitelline fields not confluent, extending from mid of pharynx to anterior region of ovary, occupying 32% (0.25–0.38, n = 15) of body length. Uterus in hindbody with several loops, passes medially between ovary and ventral sucker and forward to genital pore. Eggs small, 21 (18–24, n = 57) long by 10 (10–12, n = 57) wide. Excretory pore terminal. Excretory vesicle Y-shaped.

Figure 4. 

Paratypes of Oligogonotylus andinus sp. nov.; adults at different degrees of maturity A specimen in partly lateral orientation with the first eggs in utero; gonotyls were not observed B, D–F specimens with the seminal receptacle visible (left on B and F, right on E and D) and four gonotyls A, C, D, F specimens with the vitelline duct visible A, D, E specimens with the seminal vesicle dorsal to the ventral sucker. Scale bars: 100 μm (A–F).

Taxonomic summary

Type specimens. Holotype: permanent slide; CCH.116: 170. Paratypes: permanent slides; CCH.116: 171 to 178.

Type locality

Colombia – Antioquia • Sopetrán, La Mirandita village; 6°30'42.0"N, 75°45'22.6"W; fishponds at 750 m a.s.l., February 2020, V. Vélez-Sampedro and C. Lenis leg.

Type hosts

Andinoacara latifrons (Steindachner, 1878) (Actinopterygii: Cichlidae).

Other hosts

Cercariae in Aroapyrgus colombiensis Malek & Little, 1971 (Gastropoda, Cochliopidae); metacercariae in Poecilia caucana (Steindachner, 1880) (Actinopterygii, Poeciliidae).

Prevalence

[N (%)]. Cercariae 451 (2.0%); metacercariae 105 (38.1%); juveniles/adults 22 (100%).

Intensity

[mean (range)]. Cercariae 8.6 (2–19); metacercariae 8.3 (5–14); juveniles/adults 41.1 (5–108).

Abundance

[mean]. Cercariae (0.2); metacercariae (3.1); juveniles/adults (39.3).

Site in hosts

Cercariae in hepatopancreas; encysted metacercariae in muscle; adults in medial region of large intestine.

Etymology

The specific epithet (andinus) refers to the geographic distribution of the hosts, which are endemic to the Andes.

Key to Oligogonotylus species

1 Vitelline fields restricted to hindbody or barely extending forebody 2
Vitelline fields extended into forebody; between mid-region of pharynx and anterior margin of ovary; 3 or 4 gonotyls O. andinus sp. nov.
2 Vitelline fields between posterior margin of ventral sucker and posterior margin of posterior testis; 5–8 gonotyls O. manteri Watson, 1976
Vitelline fields between anterior margin of posterior testis and region of esophagus and pharynx; 6–8 gonotyls O. mayae Razo-Mendivil et al., 2008

Molecular phylogenetic analyses

For the analysis of the ITS2 rDNA gene region, a total of four sequences from Colombia were obtained (GenBank MW621150, MW621151, MW621152, and MW621153). The final alignment includes 18 sequences from GenBank and consists of 22 sequences of 220 bp in total. The phylogenetic tree reconstructed by ML based on the SYM+G model (-lnL = 1232.8754) shows two clades within Oligogonotylus (Fig. 5). Clade І includes five sequences of O. manteri from Mexico, Belize, and Guatemala. Clade ІІ includes four sequences of O. mayae from Yucatán. Between these clades are the four sequences that represent the new species from Colombia forming a polytomy. The sequences of the cercaria (n = 1), metacercaria (n = 1, not included in the trees), and adults (n = 4, one not included in the trees) are identical, indicating that they belong to the same species. The bootstrap values are strong to include O. andinus sp. nov. in this genus and in the family Cryptogonimidae. The intraspecific genetic distances range from 0% (O. andinus sp. nov. and O. manteri) to 0.45% (O. mayae). The interspecific genetic distances (p-distances) are 0.45–0.46% for O. andinus sp. nov. versus O. manteri, 0.91–1.39% for O. andinus sp. nov. versus O. mayae, and 1.36–1.82% for O. manteri versus O. mayae.

Figure 5. 

Phylogenetic relationships of ITS2 sequences of Oligogonotylus andinus sp. nov., Oligogonotylus manteri, and Oligogonotylus mayae. Phylogenetic tree reconstructed using the maximum-likelihood method based on the symmetrical model with gamma distribution (SYM+G). Numbers on the branches correspond to the ultrafast bootstrap approximation support values. Each terminal is identified by the alphanumeric accession number of GenBank. Terminals in orange correspond to the sequences obtained in this report. Diagram made with IQ-TREE 2.1.3.

For the analysis of the COI mtDNA gene region, five sequences from Colombia were obtained (GenBank MW658570, MW658572, MW658573, MW658574, and MW658575). The final alignment includes 22 sequences from GenBank and consists of 27 sequences of 330 bp in total. The phylogenetic tree reconstructed by ML based on the GTR+G model (-lnL = 2541.9707) shows two major monophyletic clades within Oligogonotylus (Fig. 6). Clade І includes five sequences of O. manteri from Mexico, Belize, and Guatemala. Clade ІІ includes five sequences representing the new species, all adult specimens. Clade ІІІ includes eight sequences of O. mayae from Yucatán that form a polytomy. The bootstrap values are strong (>85%) in each case supporting the position of the new species in Oligogonotylus and the family Cryptogonimidae. This topology identifies O. andinus sp. nov. and O. manteri as sister taxa with bootstrap support of 77%. The intraspecific genetic distances are 0% for O. andinus sp. nov., 0.28–3.66% for O. manteri, and 0.28–0.85% for O. mayae. The interspecific genetic distances (p-distances) are 10.27–11.28% for O. andinus sp. nov. versus O. manteri, 8.46–9.06% for O. andinus sp. nov. versus O. mayae, and 9.30–10.42% for O. manteri versus O. mayae.

Figure 6. 

Phylogenetic relationships of COI sequences of Oligogonotylus andinus sp. nov., Oligogonotylus manteri, and Oligogonotylus mayae. Phylogenetic tree reconstructed using the maximum-likelihood method based on the general time reversible with gamma distribution model (GTR+G). Numbers on the branches correspond to the ultrafast bootstrap approximation support values. Each terminal is identified by the alphanumeric accession number of GenBank. Terminals in blue correspond to the sequences obtained in this report. Diagram made with IQ-TREE 2.1.3.

The extended analysis of concatenated COI and ITS2 sequences includes four sequences from Colombia. The final alignment includes 18 sequences from GenBank and consists of 22 sequences of 330 bp in total. The phylogenetic tree reconstructed by ML based on the GTR+I+G model (-lnL = 7111.7722) shows three major clades (Fig. 7) within Oligogonotylus. Clade І includes five sequences of O. manteri from Mexico, Belize, and Guatemala; clade ІІ includes the four sequences of O. andinus sp. nov.; and clade ІІІ includes four sequences of O. mayae from Yucatán. The three clades are monophyletic within a monophyletic genus, supporting the results of the COI topology for the new species and also giving more resolution to the phylogenetic relationships within the genus. The bootstrap values are strong (>90%) for clades I, II, and III, and this topology also identifies O. andinus sp. nov. and O. manteri as sister taxa. However, the bootstrap support is low (57%). The intraspecific genetic distances are 0% for O. andinus sp. nov., 0–2.18% for O. manteri, and 0.18–0.73% for O. mayae. The interspecific genetic distances (p-distances) are 6.36–6.96% for O. andinus sp. nov. versus O. manteri, 5.64–6.04% for O. andinus sp. nov. versus O. mayae, and 6.55–7.09% for O. manteri versus O. mayae.

Figure 7. 

Phylogenetic relationships of COI and ITS2 concatenated sequences of Oligogonotylus andinus sp. nov., Oligogonotylus manteri, and Oligogonotylus mayae. Phylogenetic tree reconstructed using the maximum-likelihood method based on the general time reversible with proportion of invariable sites and gamma distribution model (GTR+I+G) model. Numbers on the branches correspond to the ultrafast bootstrap approximation support values. Terminals in red correspond to the sequences obtained in this report. Diagram made with Garli IQ-TREE 2.1.3.

The interspecific distances obtained from both the internal transcribed spacer (ITS2) and the mitochondrial gene (COI) together with the concatenated sequences provide clear support to the distinction of the new species. Our phylogenetic analysis using mitochondrial data shows that consistent with the wide distribution, the Middle American species have intraspecific mitochondrial variation —0.28–3.66% (p-distance)—over a wide latitudinal gradient from Costa Rica to Mexico. The intraspecific distance for the new species is null, consistent with all samples coming from a single Colombian population.

Remarks

Oligogonotylus andinus sp. nov. represents the first species described in the Colombian cichlid A. latifrons and the third for the genus Oligogonotylus, with O. manteri and O. mayae as being parasites of Middle American cichlids. Morphological data from cercariae and metacercariae of O. andinus are compatible with those previously reported for O. manteri (see Scholz et al. 1994). The cercaria of the new species also corresponds to the pleurolophocercous type for having a lateral fin on the tail (Schell 1970; Yamaguti 1975; Scholz et al. 1994; Bray et al. 2008) and differs from that of O. manteri in having the tegument without hair-like sensory structures (Fig. 3A); hair-like sensory structures were reported on the body margins in O. manteri.

As adults, O. andinus sp. nov. is easily distinguishable from the other species in this group based on differences in the distribution of the vitelline fields as well as the number of gonotyls. The new species has vitelline fields between the mid-level of the pharynx and the anterior margin of the ovary and three or four gonotyls. Oligogonotylus manteri has vitelline fields between the posterior margin of the ventral sucker and the posterior-end of the posterior testis and has five to eight gonotyls (Watson 1976), whereas O. mayae has vitelline fields from a level between the pharynx and the esophagus and the anterior margin of the posterior testis and six to eight gonotyls (Razo-Mendivil et al. 2008).

The analysis of the morphology of the reproductive system of O. andinus sp. nov. at different stages of maturity (Figs 3D, 4A–F) reveals that the seminal receptacle is ventral to the vitelline duct; the ootype complex overlaps the lateral region of the ovary dorsally (on right side in the type specimen and three paratypes but on left side in three other paratypes, on the same side of the seminal receptacle); and the distal region of the uterus and the seminal vesicle reach the anterior margin of the ventral sucker and converge at the genital pore. These features are not detailed in the description of O. manteri. In O. mayae, as in the new species, the seminal vesicle is tubulo-saccular and posterolateral to the ventral sucker.

Morphology-based taxonomic conclusions are supported by the analysis of molecular data from the ITS2, COI and COIITS2 barcoding. Our analyses show that the specimens from Colombia represent a new species. The polytomy formed by the sequences of the new species in the phylogenetic tree of ITS2 is due to the low rate of mutation that this molecular marker has. In conjunction, the phylogenetic trees of COI and the concatenated sequences of COI and ITS2 show that all isolates from Sopetrán represent a reciprocally monophyletic assemblage, distinct from isolates of O. manteri and O. mayae, which are shown in both topologies to be sister taxa.

Discussion

In the present study, O. andinus sp. nov. is described and associated with three endemic Andean hosts from Sopetrán, Antioquia, Colombia. Initially, molecular analyses of cercariae which had emerged from A. colombiensis enabled identification to the generic level. Subsequently, the location of the focus in the fishponds facilitated the life cycle study using five species of fishes present in the fish farm; two are endemic species (P. caucana and A. latifrons), and three are exotic cultured species (O. aureus, C. macropomum, and O. mossambicus). Of these, the two endemic fishes act as intermediate and definitive hosts of the parasite, respectively. The high prevalence of O. andinus in P. caucana (38.1%) and A. latifrons (100%) are a result of at least three factors: (a) the confinement of the snails and fishes in small semi-closed ponds, which increases the rate of parasite-host encounters; (b) the establishment of a large population of P. caucana in the fishpond, with few predators (A. latifrons and birds such as kingfishers); and (c) a voracious diet of A. latifrons that included snails, poecilids, and insects (judging from the stomach content; present study). Despite their close confinement, the farmed fish were negative for the new species because their diets are based on plants and fish feed.

Previous studies have shown that both O. manteri and O. mayae use the cochliopid snail P. coronatus (misidentified as Benthonella gaza), and that O. manteri also uses the cochliopid snail Aroapyrgus alleei Morrison, 1946 as the first intermediate host and the cichlid fish M. urophthalmus either as the second intermediate or the main definitive host among other 11 species of cichlids (Scholz et al. 1994; Jiménez-García and Vidal-Martínez 2005; Pérez-Ponce de León et al. 2007; May-Tec et al. 2020; Velázquez-Urrieta and Pérez-Ponce de León 2021). Oligogonotylus andinus sp. nov. is related to hosts whose history has been shaped by the Andes Mountains and whose life cycle occurs in lentic systems. A cochliopid snail (A. colombiensis) is the first intermediate host, a poecilid fish (P. caucana) is the second intermediate host, and a trans-Andean cichlid fish (A. latifrons) is the definitive host; thus, the present work extends the OligogonotylusCochliopidaeCichlidae association. Therefore, species of Oligogonotylus are still considered as members of the core helminth fauna of cichlids (Pérez-Ponce De León and Choudhury 2005), but also as members of the core helminth fauna of cochliopids (May-Tec et al. 2020; Velázquez-Urrieta and Pérez-Ponce de León 2021; present study).

Aroapyrgus colombiensis is distributed in the Magdalena–Cauca basin and inhabits small lotic systems with abundant stretches of slow water and decomposing dead leaves (Linares et al. 2018). It has been reported as a host for Paragonimus caliensis Little, 1968 in Valle del Cauca (Malek and Little 1971) and Antioquia (Lenis et al. 2018). Poecilia caucana, which is commonly called Pipón, Piponcita, or Cauca Molly, is distributed in the hydrographic slopes of the Pacific and the Caribbean, especially in the Magdalena–Cauca basin (Castellanos-Morales et al. 2011; Martínez et al. 2016; Mojica et al. 2018; Ruiz-Guerra and Echeverry-Galvis 2019). Normally it inhabits rivers and brooks (Roman-Valencia 1990), but small streams also lead it to lagoons and natural or artificial ponds where it is reported in Magangué, Bolívar, Los Palmitos, Sucre (Navarro et al. 2019), and Sopetrán (present study). Andinoacara latifrons is colloquially named Blue mojarra and it is grouped in two lineages in Colombia, (1) the upper Magdalena and Catatumbo clade, and (2) the upper Cauca and upper Magdalena clade (De la Ossa-Guerra et al. 2020), where our host population belongs. It inhabits muddy waters, clean-water streams, and lentic environments. It feeds on worms, insects, plant materials, fish remains, debris, and other live prey, such as crustaceans (Maldonado-Ocampo et al. 2006; Olaya-Nieto et al. 2016), making it susceptible to infection by trematodes.

This is the first study of the helminth fauna associated with A. colombiensis and A. latifrons and extends the known geographic distribution of Oligogonotylus spp. from Middle America to South America (Choudhury et al. 2016, 2017) where O. andinus sp. nov. was found at 750 m above sea level. Further field-based studies must include an increased sampling size to determine whether major intraspecific distances exist or if there are other populations in the Cauca river basin, and to reveal the ecogeographic limits of Oligogonotylus spp. in Colombia, particularly in areas with a likely suitable habitat where the species has yet to be found. Future taxonomic and ecological host-parasite research must lead to an understanding of fish parasite diversity and the risk of the transmission of zoonotic diseases, considering the dependence on fish as a food resource in rural areas of Colombia.

Acknowledgements

We acknowledge Carlos Muskus (PECET), Iván Darío Vélez (PECET), and Javier Correa Álvarez (EAFIT) for their support in conducting this research. We also thank Enderson Murillo Ramos (UBMC-PECET), Juan Manuel Martínez Cerón (EAFIT), and Jorge Luis Escobar (UdeA) for their assistance, recommendations, and contributions regarding the methodology used; we highly appreciate it. Special thanks go to Gildardo Méndez (Yayo) and his family for allowing us to carry out this research in their fish farm. This work is funded by the Helminthology Unit-PECET and the Ministry of Science, Technology and Innovation-MINCIENCIAS (grant 80740-571-2020). The authors declare that they have no competing interests exist.

References

  • Anucherngchai S, Tejangkura T, Chontananarth T (2016) Epidemiological situation and molecular identification of cercarial stage in freshwater snails in Chao-Phraya Basin, Central Thailand. Asian Pacific Journal of Tropical Biomedicine 6(6): 539–545. https://doi.org/10.1016/j.apjtb.2016.01.015
  • Bowles J, Hope M, Tiu WU, Liu X, McManus DP (1993) Nuclear and mitochondrial genetic markers highly conserved between Chinese and Philippine Schistosoma japonicum. Acta Tropica 55(4): 217–229. https://doi.org/10.1016/0001-706X(93)90079-Q
  • Bush AO, Lafferty K, Lotz J, Shostak W (1997) Parasitology meets ecology on its own terms: Margolis et al. revisited. The Journal of Parasitology 83(4): 575–583. https://doi.org/10.2307/3284227
  • Castellanos-Morales CA, Marino-Zamudio LL, Guerrero-V L, Maldonado-Ocampo JA (2011) Peces del departamento de Santander, Colombia. Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales 35: 189–212.
  • Choudhury A, Aguirre-Macedo ML, Curran SS, Ostrowski de Núñez M, Overstreet RM, Pérez-Ponce de León G, Portes Santos C (2016) Trematode diversity in freshwater fishes of the globe II: ‘New World.’. Systematic Parasitology 93(3): 271–282. https://doi.org/10.1007/s11230-016-9632-1
  • Choudhury A, García-Varela M, Pérez-Ponce De León G (2017) Parasites of freshwater fishes and the Great American Biotic Interchange: A bridge too far? Journal of Helminthology 91(2): 174–196. https://doi.org/10.1017/S0022149X16000407
  • Darriba D, Taboada G, Doallo R, Posada D (2012) jModelTest 2: More models, new heuristics and parallel computing. Nature Methods 9(8): 772–772. https://doi.org/10.1038/nmeth.2109
  • De la Ossa-Guerra LE, Santos MH, Artoni RF (2020) Genetic diversity of the cichlid Andinoacara latifrons (Steindachner, 1878) as a conservation strategy in different Colombian basins. Frontiers in Genetics 11: 1–9. https://doi.org/10.3389/fgene.2020.00815
  • Ditrich O, Scholz T, Aguirre-Macedo L, Vargas-Vázquez J (1997) Larval stages of trematodes from freshwater molluscs of the Yucatan Peninsula, Mexico. Folia Parasitologica 44: 109–127.
  • Edgar R (2004) MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32(5): 1792–1797. https://doi.org/10.1093/nar/gkh340
  • Erikson U (2011) Assessment of different stunning methods and recovery of farmed Atlantic salmon (Salmo salar): isoeugenol, nitrogen and three levels of carbon dioxide. Animal Welfare, the UFAW Journal 20: 365–375.
  • Hoang DT, Chernomor O, Von Haeseler A, Minh BQ, Vinh LS (2018) UFBoot2: Improving the ultrafast bootstrap approximation. Molecular Biology and Evolution 35(2): 518–522. https://doi.org/10.1093/molbev/msx281
  • Iwagami M, Ho LY, Su K, Lai PF, Fukushima M, Nakano M, Blair D, Kawashima K, Agatsuma T (2000) Molecular phylogeographic studies on Paragonimus westermani in Asia. Journal of Helminthology 74(4): 315–322. https://doi.org/10.1017/S0022149X00000469
  • Jiménez-García MI, Vidal-Martínez VM (2005) Temporal variation in the infection dynamics and maturation cycle of Oligogonotylus manteri (Digenea) in the cichlid fish, “Cichlasomaurophthalmus, from Yucatán, México. The Journal of Parasitology 91(5): 1008–1014. https://doi.org/10.1645/GE-380R2.1
  • Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A (2012) Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics (Oxford, England) 28(12): 1647–1649. https://doi.org/10.1093/bioinformatics/bts199
  • Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018) MX (2018) MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology and Evolution 35(6): 1547–1549. https://doi.org/10.1093/molbev/msy096
  • Lenis C, Galiano A, Vélez I, Vélez ID, Muskus C, Marcilla A (2018) Morphological and molecular characterization of Paragonimus caliensis Little, 1968 (Trematoda: Paragonimidae) from Medellin and Pichinde, Colombia. Acta Tropica 183: 95–102. https://doi.org/10.1016/j.actatropica.2018.03.024
  • Linares EL, Lasso CA, Vera-Ardila LM, Morales-Betancourt AM (2018) XVII. Moluscos dulceacuícolas de Colombia. Serie Editorial Recursos Hidrobiológicos y Pesqueros Continentales de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. Bogotá, 326 pp. https://doi.org/10.21068/a2018n03
  • Maldonado-Ocampo JA, Ortega-Lara A, Usma Oviedo JS Germán, Galvis V, Villa-Navarro FA, Vásquez Gamboa L, Prada-Pedreros S, Ardila Rodríguez C (2006) Peces de Los Andes de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt, Bogotá, 346 pp.
  • Martínez JD, Cadena CD, Torres M (2016) Critical thermal limits of Poecilia caucana (Steindachner, 1880) (Cyprinodontiformes: Poeciliidae). Neotropical Ichthyology 14(1): 1–7. https://doi.org/10.1590/1982-0224-20150171
  • May-Tec AL, Herrera-Castillo NA, Vidal-Martínez VM, Aguirre-Macedo ML (2020) Following the infection dynamics of the tropical trematode Oligogonotylus mayae in its intermediate and definitive hosts for 13 years. Journal of Helminthology 94: e208. https://doi.org/10.1017/S0022149X20000875
  • Mojica JI, Martínez-González CC, Acosta-Vela AG, Larrate-River E, González-Daza W, Ávila-Rojas FL, Martínez-Aguirre E, Forero-Cano JD (2018) Lista de los peces de la cuenca del río Mira, vertiente Pacífico, Colombia. Biota Colombiana 18(2): 191–199. https://doi.org/10.21068/c2017.v18n02a12
  • Navarro B, Tovar H, Caraballo P (2019) Composición y distribución de la ictiofauna asociada a jagüeyes, en la región Caribe colombiana. Intropica 14: 120–126. https://doi.org/10.21676/23897864.3277
  • Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evolution 32(1): 268–274. https://doi.org/10.1093/molbev/msu300
  • Olaya-Nieto CW, Camargo-Herrera L, Díaz-Sajonero V, Segura-Guevara FF (2016) Feeding habits of Cocobolo Andinoacara pulcher in the cienaga Grande de Lorica, Colombia. Revista Mvz Cordoba 21: 5189–5197. https://doi.org/10.21897/rmvz.29
  • Pérez-Ponce De León G, Choudhury A (2005) Biogeography of helminth parasites of freshwater fishes in Mexico: The search for patterns and processes. Journal of Biogeography 32(4): 645–659. https://doi.org/10.1111/j.1365-2699.2005.01218.x
  • Razo-Mendivil U, Rosas-Valdez R, Pérez-Ponce de León G (2008) A new cryptogonimid (Digenea) from the Mayan cichlid, Cichlasoma urophthalmus (Osteichthyes: Cichlidae), in several localities of the Yucatan Peninsula, Mexico. The Journal of Parasitology 94(6): 1371–1378. https://doi.org/10.1645/GE-1546.1
  • Roman-Valencia C (1990) Lista y distribución de peces en la cuenca medio del río Atrato, Chocó, Colombia. Caldasia 16: 201–208.
  • Ruiz-Guerra C, Echeverry-Galvis MÁ (2019) Prey consumed by wading birds in mangrove swamps of the Caribbean coast of Colombia. Journal of Natural History 53(29–30): 1823–1836. https://doi.org/10.1080/00222933.2019.1667037
  • Sandlund OT, Daverdin RH, Choudhury A, Brooks DR, Diserud OH (2010) A survey of freshwater fishes and their macroparasites in the Guanacaste Conservation Area (ACG), Costa Rica. NINA Report 935: 1–45. http://hdl.handle.net/11606/639
  • Schell S (1970) How to Know the Trematodes. Wm. C. Brown Company Publishers, Dubuque, Iowa, 355 pp.
  • Scholz T, Lavadores JIP, Vargas VJ, Mendoza FEF, Rodriguez CR, Vivas RC (1994) Life cycle of Oligogonotylus manteri (Digenea: Cryptogonimidae), a parasite of cichlid fishes in southern Mexico. Journal of the Helminthological Society of Washington 61: 190–199. http://bionames.org/bionames-archive/issn/1049-233X/61/190.pdf
  • Scholz T, Vargas-Vázquez J, Moravec F, Vivas-Rodríguez C, Franco-Mendoza E (1995) Cenotes (sinkholes) of the Yucatan Peninsula, Mexico, as a habitat of adult trematodes of fish. Folia Parasitologica 42: 37–47.
  • Sosa-Medina T, Vidal-Martínez VM, Aguirre-Macedo ML (2015) Metazoan parasites of fishes from the Celestun coastal lagoon, Yucatan, Mexico. Zootaxa 4007(4): 529–544. https://doi.org/10.11646/zootaxa.4007.4.4
  • Velázquez-Urrieta Y, Pérez-Ponce de León G (2021) Morphological and molecular assessment of the diversity of trematode communities in freshwater gastropods and bivalves in Los Tuxtlas tropical rainforest. Journal of Helminthology 95: E44. https://doi.org/10.1017/S0022149X21000407
  • Vergara D, Velásquez LE (2009) Larval stages of digenea from Melanoides tuberculata (Gastropoda: Thiaridae) in Medellín, Colombia. Acta Biologica Colombiana 14: 135–142.
  • Yamaguti S (1975) A Synoptical Review of the Life Histories of Digenetic Trematodes of Vertebrates, with Special Reference to the Morphology of Their Larval Forms. Keigaku Publishing Co., Tokyo, 575 pp.
login to comment