Next Article in Journal
Novel Techniques for Musculoskeletal Pain Management after Orthopedic Surgical Procedures: A Systematic Review
Previous Article in Journal
Mineral Indicators of Geologically Recent Past Habitability on Mars
Previous Article in Special Issue
Baseline Corticosterone, Stress Responses, and Leukocyte Profiles in Chicks of Precocial Birds in Rural and Urban Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphometric and Molecular Analysis of Five-Spine Epidinium Morphotypes Taken from the Rumen of European Bison, Bison bonasus

Institute of Animal Physiology, Centre of Biosciences, Slovak Academy of Sciences, 04001 Košice, Slovakia
*
Author to whom correspondence should be addressed.
Life 2023, 13(12), 2350; https://doi.org/10.3390/life13122350
Submission received: 14 November 2023 / Revised: 12 December 2023 / Accepted: 13 December 2023 / Published: 15 December 2023
(This article belongs to the Special Issue Research on Wildlife Behavior and Biodiversity)

Abstract

:
An important feature of ruminal ciliates is their phenotypic plasticity, which makes their identification difficult. The common manifestation of the phenotypic plasticity in rumen ciliates is a change in their cell size and caudal spination. We analyzed various morphotypes of Epidinium with five caudal processes (spines) taken from the rumen of European bison (Bison bonasus). In the study, the cluster analysis and K-means analysis of morphometric data could not distinguish very similar morphotypes of Epidinium with five caudal processes. However, the morphotype of E. parvicaudatum prevailed (70%). The DNA of four individual E. parvicaudatum was isolated successfully from formaldehyde-preserved samples. The partial 18S rDNA gene sequences (about 350–400 bp) were identical to Epidinium sequences in GenBank (E. caudatum, a one-spine morphotype, and E. cattanei, a five-spine morphotype). It can be assumed that these short sequences cannot distinguish the differences between the Epidinium morphospecies. Complete gene sequences from various hosts and various molecular markers are necessary to reveal the validity of the Epidinium five-spine species. In conclusion, classical morphology should be supplemented with molecular data when more morphotypes of the rumen ciliate species are present in samples.

1. Introduction

A species concept is a working definition of a species or a methodology for determining whether or not two organisms are members of the same species. Scientists commonly use three species concepts. All species concepts have limitations and work better for some organismal groups than others. The biological species concept states that a species is a group of actually or potentially interbreeding natural populations reproductively isolated from other such groups. The lineage (phylogenetic) species concept defines species as groups of organisms that share a pattern of ancestry and descent and form a single branch on the tree of life. The morphological species concept is frequently applied when there is limited or unavailable data on reproductive behavior or genetic similarity. Morphological species (morphospecies) is the historical concept of species, accepted by many generations of naturalists, zoologists, and even today by many biologists [1,2,3,4,5]. Generally, two major approaches are used for the identification of ciliates: the traditional method of live-cell morphology, fixation, staining, and morphometrics, and, more recently, DNA-based methods [6,7]. Nowadays, taxonomy is more focused on an integrative approach combining classical and molecular approaches [1,8,9]. Phenotypic plasticity is the ability of an organism to change in response to stimuli or environmental inputs. The response may involve a change in physiological state, behavior, morphology, or some combination [10]. Ciliated protozoa colonize the forestomach (rumen) of all domestic and wild ruminants and pseudoruminants (llama, camel). Rumen commensal ciliates (Protista, Ciliophora, Litostomatea, Trichostomatia) take an active part in the fermentation processes in the rumen. We know more than 225 species of rumen ophryoscolecid ciliates described morphologically [11,12,13]. Their phenotypic plasticity makes identification difficult [14,15]. The common manifestation of the phenotypic plasticity in rumen ciliates is a change in cell size and caudal spination. Entodinium, Ophryoscolex, and Epidinium genera often react in this way. Under in vitro conditions, their spinated morphotypes gradually lose their caudal projections [16,17,18,19,20,21,22]. The genus Epidinium Crawley (1923) (family Ophryoscolecidae, order Entodiniomorphida) consists of spineless and spinated morphotypes. E. capricornisi [23], E. cattanei (syn. E. ecaudatum f. fasciculus) [24,25], E. parvicaudatum [12,25,26], and E. rusa [27] represent five-spine morphospecies between which confusion may occur in ecological studies [11,28]. They differ in the size and shape of their cell body and spines. In addition, there is a discrepancy in the description and names of E. parvicaudatum, E. cattanei, and E. fasciculus among some authors [11,16,24,25,26,28,29,30]. In some studies, the descriptions and names of Dogiel monography are preferred [12,16,30]. However, E. fasciculus is not a valid name for E. cattanei [11,12,29]. There has been much discussion about whether similar epidinia, which differ only in the number of caudal spines, are different species [11,25,29]. Dogiel suggested that epidinia with varying numbers of spines are morphological variations (forms) of one species, E. ecaudatum [25], a spineless morphotype. On the other hand, the “ecaudatum group” was established with Epidinium ecaudatum as a type species characterized by an elongated body tapering posteriorly with no spine by Kofoid and McLennan [29]. The group consists of morphotypes with 0–5 caudal spines, e.g., E. caudatum, E. bicaudatum, E. tricaudatum, and E. parvicaudatum. On the contrary, E. hamatum, E. cattanei, and E. bulbiferum differ from the ecaudatum group by the body’s thick and heavy posterior part. Spinated forms also occur there [29]. E. cattanei (syn. E. ecaudatum f. fasciculus) and E. parvicaudatum (syn. E. ecaudatum f. parvicaudatum) are the relatively common five-spine morphotypes. E. cattanei, Fiorentini (1889), is characterized by its relatively short body and five long, straight caudal spines in contrast to E. parvicaudatum, Awerinzew and Mutafowa (1914), which is characterized by body tapering posteriorly with five incurved, heavy caudal spines of which the ventral is the largest [11,12,24,25,26]. In the rumen of European bison, Epidinium morphotypes with 0–5 spines occurred there [31]. Five-spine morphotypes prevailed, and they were assigned to E. parvicaudatum. However, other morphotypes were observed, and many individuals were similar to E. cattanei. In the present study, the authors hypothesized that different morphotypes might be differentiated by a multivariate analysis of morphometric data. To date, only two Epidinium morphotypes, E. caudatum (syn. E. ecaudatum f. caudatum, one-spine morphotype) and E. cattanei (syn. E. ecaudatum f. fasciculus, five-spine morphotype) have been molecularly analyzed [30,32]. E. cattanei formed a distinct group from E. caudatum [6,30,32]. We isolated the DNA from individual Epidinium five-spine morphotypes of European bison samples preserved by formaldehyde to resolve their molecular taxonomy.

2. Materials and Methods

The study focused on analyzing the five-spine morphotypes of Epidinium of the rumen ciliate population of the European bison (Bison bonasus). In February and March 2007, rumen content samples were obtained from 18 culled European bison in Białowieża Forest, Poland [31]. Samples were collected immediately after death (1–3 h). Approximately 5 g of the rumen content from each animal was preserved with an equal amount of 8% formaldehyde. Preserved content was filtered through four layers of gauze and transported to the Institute of Animal Physiology of the Slovak Academy of Sciences in Slovakia, where they were kept in the fridge [31].

2.1. Morphometric Analysis

The pictures of Epidinium five-spine morphotypes were taken under bright-field illumination by a Moticam Pro CCD Camera (Motic Incorporation Ltd., Hong Kong) mounted on a BA400 microscope (Motic Incorporation Ltd., Hong Kong). The wet rumen samples on microscope slides were colored by Lugol’s iodine solution. We have analyzed the set of 272 pictures of Epidinium from 15 European bison. We manually measured the ciliate length and width on the calibrated pictures using ImageJ (1.52a, Wayne Rasband, National Institutes of Health, Kensington, MD, USA). The spines were not included in the length measurement because they had a variable length. We measured the length in the median plane from the apical pole to the caudal end. The width was measured where the cell was the widest (usually in the middle plane) and the line was perpendicular to the length line. The ratio of cell length to width was calculated. The morphometric data (length and width) were statistically evaluated in Microsoft Excel (Microsoft Office 2016 Professional Plus), IBM SPSS, version 19 (IBM Corporation, Armonk, NY, USA), and PAST, version 4.04 software [33]. IBM SPSS was used for the Kolmogorov–Smirnov normality test. The multivariate analysis of morphometric data similarity (clustering and K-means analysis) was carried out according to the PAST guide [33]. The orientation of the epidinia is described according to the position of skeletal plates. It is said that the cytostome with an adoral ciliary zone (AZM) and skeletal plates are on the ventral side of the ciliate [11,34]. The dorsal ciliary zone (DZM) and two contractile vacuoles are on the dorsal ciliate side. When the ciliate is rotated with the skeletal plates on the right side, the upper side of the image can be assigned to the right side, and the lower side of the image can be given to the left side of the ciliate. In this study, we use the original two-word nomenclature of individual morphotypes and, in parentheses, the equivalent terminology, according to Dogiel [11,25].

2.2. Molecular Analysis

DNA isolation. To resolve the taxonomy of the five-spined epidinia, the DNA was isolated from 25 individual five-spine Epidinium with a long ventral caudal spine from the formaldehyde-preserved rumen fluid of European bison. DNA from individual protozoal cells was isolated by the chelex method [35]. The individual ciliate cell was picked out from the drop of formaldehyde-preserved rumen fluid under a microscope using a micromanipulator fitted with a glass capillary. The cell was washed three times in 10 µL of distilled sterile water on a glass slide and transferred to a mini tube filled with 50 µL of 5% Chelex-100 (BIO-RAD, Hercules, CA, USA) in water and vortexed for 15 s. Pre-incubated proteinase K (Merck, Rahway, NJ, USA) was then added to the reaction solution to a final concentration of 10 µg/mL and incubated for 20 min at 98 °C). After incubation, the solution was cooled to 0 °C, centrifuged (10,000 rpm, 10 min, and 4 °C), and the supernatant was used directly as a template for PCR amplification.
Polymerase-chain-reaction. For PCR amplification of the 18S rDNA gene, two primer sets synthesized by Sigma-Aldrich, St. Louis, MO, USA) were used. The primer sets EukFor (5′-AATATGGTTGATCCTGCCAGT-3′) and EukRev (5′-TGATCCTTCTGCAGGTTCACCTAC-3′) were able to amplify the nearly full-length 18S rRNA gene (of about 1500 bp), and primer sets Cill-R-GC (5′-CGCCCGCCGCGCCCCGCGCCCGGCCCGCCGCCCCCGCCCGGGGCCAATTGCAAAGATCTATCCC-3′) and Cill-F (5′-GGTGGTGCATGGCCG-3′) were able to amplify the first 400 bp of the 18S rDNA gene (Moon van der Stay, unpublished). PCR amplification was carried out using a C1000TM Thermal Cycler (BIO-RAD, Hercules, CA, USA). For amplification of the full-length 18S rDNA gene, the conditions specified by Regensbogenova et al. were used [27]. However, the original method was unsuccessful for PCR amplification [27]. The modified method was as follows: The Cill-R_GC and Cill-F primer sets were used for amplification. The PCR reaction was performed as follows: initial denaturation (94 °C, 4 min), 35 cycles of denaturation (94 °C, 1 min), annealing (55 °C, 30 s), extension (72 °C, 1 min), denaturation (94 °C, 1 min), annealing (65 °C, 30 s), and elongation (72 °C, 10 min). Electrophoresis in 1% agarose gel and 100 bp ladder as standard of molecular size were used to determine the quality and quantity of the obtained PCR products. Wizard SV Gel and PCR Clean-Up System (Promega Corporation, Madison, WI, USA) were used to purify PCR products. The concentration of DNA was measured by a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific Inc., Colorado Springs, CO, USA). Obtained amplicons were sequenced using the same primer pair for amplification at Eurofins Genomics (Ebersberg, Bayern, Germany). Obtained sequences were assembled using DNA Baser ver. 2.0 (Heracle BioSoft SRL, Cluj-Napoca, Romania) and analyzed using the Blastn algorithm against the GenBank database. For phylogenetic comparisons, the 18S rDNA sequences of Epidinium spp. and related species were downloaded from the GenBank database and aligned using the Clustal W algorithm. The neighbor-joining tree was constructed using the MEGA software version 7 with 1000 bootstrap repetitions [36]. The list of the sequences of E. parvicaudatum of European bison are summarized in Appendix A, Table A1.

3. Results

3.1. Morphometric Analysis

About 84% of epidinia cells with five spines had cylindrical bodies with wide body ends (Figure 1A,C–F). About 16% of epidinia were slightly tapering posteriorly (Figure 1B,J). The cells were either straight or slightly bent on the ventral side. The dorsal side was usually more convex when viewed from the right or left side. According to the shape of the spines, four morphotypes can be distinguished in the rumen of European bison. Five incurved, heavy caudal spines, of which the ventral was the largest, as is typical for E. parvicaudatum, were observed in 70% of individuals. (Figure 1A,B,E,F). They had a specific arrangement of spines: a single largest spine (sp1) on the ventral surface, two spines on the right surface (sp2 and sp3), a single spine on the dorsal surface (sp4), and a single spine (sp5) on the left surface (Figure 1E), [18]. The second morphotype (3%) had spines incurved at an angle greater than 80 degrees, usually towards the opposite side (Figure 1D,H). About 11% of epidinia had caudal spines of approximately the same length, the morphotype similar to E. cattanei. (Figure 1C,G,K). About 16% of epidinia had five short, incurved spines similar to E. rusa (Figure 1L). No individuals with furcated spines (E. capricornisi) were observed. The two contractile vacuoles on the dorsal side were regularly visible in the opening stage (Figure 1J). Epidinia had three skeletal plates on the ventral side that were typically well visible after staining with an iodine solution (Figure 1J). When epidinia were full of starchy feed, the ciliates were too dark to distinguish the intracellular organelles after iodine staining. The morphometric descriptive parameters (length, width, and ratio of length to width) of the five-spine Epidinium cells are summarized in Table 1. The values of the parameters are as follows: The length (L) varied from 67.5 µm to 144.4 µm, and the mean value was 103.1 µm. The width (W) ranged from 35.9 µm to 89.9 µm; the mean value was 57.4 µm. The ratio of L/W varied from 1.34 to 2.36, and the mean value was 1.8.
The K-means clustering analysis revealed the two main clusters in the set (Figure 2). The groups differ by ciliate size. Smaller epidinia (60%) were clustering in cluster 1. In total, 40% of epidinia were clustering in cluster 2. Epidinia in both clusters had the same mean ratio of L/W (1.8).
All morphological features described above were present in both groups and each animal (various spines and endings of the body). There were also dividing cells in both groups. In the group of smaller ciliates (cluster 1; Figure 1A–D), 23% of ciliates were in different dividing stages. In cluster 2 (larger ciliates, Figure 1E,F,H,J,L), 37% of ciliates were in the different dividing stages. The omitting of dividing ciliates from clustering had no substantial effect on the K-means analysis. The authors observed 29% of the morphotype of E. parvicaudatum (long ventral spine) in cluster 1 and 41% in cluster 2, with a prevalence of the wide-ending ciliate bodies. The morphometric descriptive parameters of both groups are summarized in Table 2.

3.2. Molecular Analysis

To resolve the taxonomy of the five-spine epidinia, the DNA was isolated from 25 individuals of five-spine Epidinium with a long ventral caudal spine from the formaldehyde-preserved rumen fluid of European bison. Only four samples were successfully sequenced. No full-length 18S rRNA PCR amplicons from formaldehyde single-cell Epidinium sp. samples were obtained using the EukFor and EukRev primer pair. Using primer-pair targeting, the first 400 nucleotides of the 18S rDNA gene (Cill) yielded specific amplicons from four cells. The phylogenetic analysis of the sequences is presented in Figure 3. These short sequences of the 18S rDNA gene clustered together with Epidinium sequences in GenBank (E. caudatum, a one-spine morphotype, and E. cattanei, a five-spine morphotype) isolated from the ruminal content of several ruminants from geographically distant localities (Figure 3). Short sequences of unidentified Epidinium sp. (JN116217-220) clustered in different group.

4. Discussion

This study focused on analyzing the five-spine morphotypes of Epidinium from the rumen ciliate population of the European bison. Based on the morphometry and microscopic observation, it is hard to decide if the observed Epidinium five-spine morphotypes could be separate species. However, no individuals with furcated spines (E. capricornisi) [2,14,19] were observed in European bison. In the study, other five-spine morphospecies were observed. They resembled E. parvicaudatum (70%), E. cattanei (11%), and E. rusa (16%). There was an as-yet undescribed form with spines curved inward at an angle greater than 80 degrees (3%). These morphological variations were observed in all European bison examined. It is the first observation of the simultaneous coexistence of four morphotypes of five-spine epidinia. It is an observation of the remarkable plasticity of rumen ciliates. Similarly, Dehority observed the remarkable plasticity of Entodinium dubardi from the blue duiker [6]. He observed several forms of E. dubardi similar to other non-spine Entodinium species. Polymorphism was also observed in the species of Diplodinium and Ophryoscolex [5,7,8]. According to previous descriptions, E parvicaudatum has been detected worldwide in ruminants (cattle, gaur, sheep, goats, and giraffes) [2]. E. cattanei has been described from cattle, zebu, sheep, and goats worldwide [2]. E. rusa has been described from sambar deer in India [18]. Compared to the previous descriptions of E. parvicaudatum and E. cattanei [2], the dimensions of the five-spine epidinia were in the range of both morphotypes. The K-means analysis provided the clustering of epidinia according to size. We can speculate if these clusters represent the two different species; all morphological features used to distinguish epidinia species (various spines and endings of the body) were present in both groups. The populations of ruminal protozoa of the same morphotype differing in cell size were observed in some ecological studies [22]. The large and small Eudiplodinium maggii populations were observed in the rumen of European bison [22]. It is supposed that it could be a reaction to predation. Entodinia, epidinia, and other ciliates are objects of predation by other rumen ciliates, usually by Polyplastron, Elytroplastron, and Entodinium bursa [11,31,37,38,39,40]. The epidinia were predated by Elytroplastron bubali in European bison [22]. Therefore, the larger size may be a strategy to prevent predation. It is consistent with the assumption that larger and broader ciliates resist predation better. It can be said that the variability of spines is not surprising. It is the primary protection against predation. The ability to form different spines probably allows for better survival of the species. In addition, ciliates without a long spine may represent a population of upper daughter cells after the division of the mother cell (Figure 1I). During division, short caudal spines were visible in upper daughter cells before separation. The spines likely grow after the cells have separated. However, the presence of a long barb seems to be advantageous, as this form was the most numerous. The plasticity of the spination was proven during cultivation experiments. Clarke reported that Epidinium mixed culture (individuals with 1–5 spines) lost all caudal spines after seven months in vitro [12]. Coleman et al. reported a similar loss of the caudal spine in Epidinium caudatum culture after 4–5 months in vitro [10]. Dehority used clonal cultures of E. caudatum (one-spine form), E. tricaudatum, and E. parvicaudatum [9]. During a short period (about seven days), he reported bifurcated and trifurcated spines in epidinia and changes in the counts of spines. In the study of Miltko et al., the length of spines of E. cattanei (syn. E. ecaudatum f. fasciculus) only decreased after 12 months in vitro [8]. It seems that the five-spined form was relatively stable. The original and other descriptions of E. parvicaudatum describe long ciliates with posteriorly tapering cell endings and E. cattanei as relatively short, wide ciliates [2,15]. The morphotype description contrasts our observations, where most ciliates with a long ventral spine had wide bodies along their entire length. The tapering of the cell end does not appear to be a valid character for classifying Epidinium morphotypes, nor does the size and shape of the spines. For now, the authors are inclined to believe that the different five-spined morphotypes in European bison probably do not represent separate species.
The molecular taxonomy of European bison five-spine Epidinium was carried out to resolve their taxonomical relationship. According to the authors, this is the first report on the isolation of DNA from samples of rumen ciliates preserved by formaldehyde. The lack of amplification of the full 18S rDNA gene was probably due to extensive DNA fragmentation during formaldehyde-preserved sample storage [41]. The DNA fragmentation extracted from formalin-preserved paraffin-embedded specimens (FFPE) is well known. Formaldehyde functions as a cross-linking agent for cell (tissue) preservation. It stabilizes the tissue (cell) structure by creating a covalent linkage between macromolecules, such as DNA–DNA, DNA–protein, and protein–protein [42]. Reversing the formaldehyde-formed cross-linking during DNA extraction causes DNA fragmentation [42]. The degree of damage is affected not only by the preservation procedure but also by the storage time and storage conditions [43]. The significant degradation of DNA was observed after 4 to 6 years of sample storage [43]. On the other hand, DNA isolated from FFPE samples kept at a low temperature (4 °C) was less fragmented than DNA from samples kept at laboratory temperature [44]. However, many studies have demonstrated the feasibility of using FFPE DNA for next-generation sequencing, including gene-specific targeted sequencing, whole-exome, and whole-genome sequencing [45]. Lin et al. [46] successfully amplified DNA from FFPE samples. They assumed that DNA from FFPE samples with fragments less than 200bp may be suitable for genetic studies. Even though long-term storage of the European bison samples led to intense DNA degradation, we obtained shorter PCR amplicons. In our study, using a primer pair targeting the first 400 nucleotides of the 18S rDNA gene (primers Cill-R_GC and Cill-F) and modifying PCR reactions yielded specific amplicons from four cells. The main changes in the PCR protocol were made in the annealing and extension phases.
The sequences of Epidinium of European bison were identical to several Epidinium isolates in the GenBank. However, it is possible that the short sequences cannot distinguish the differences between the Epidinium morphospecies. It cannot be excluded that epidinia from European bison are genetically distinct from other Epidinium morphospecies. In addition, the sequences of unidentified Epidinium sp. (JN116217-220) clustered in another group. They are only short partial sequences of the 18S rDNA gene. No morphometric description can be found on the Epidinium sp. Until now, only the 18S rDNA genes have been used in the molecular analyses of the genus Epidinium and other rumen ciliates. Furthermore, few Epidinium species have been sequenced. To date, the complete sequences of 18S rDNA genes of the species E. caudatum (E. ecaudatum f. caudatum) and E. cattanei (syn. E. ecaudatum f. fasciculus) are available [21,23]. Different Epidinium morphotypes were not molecularly analyzed. Phylogenetic analysis of the complete sequences of 18S rDNA genes showed differences between E. caudatum and E. cattanei [21,23]. The full-length 18S rRNA sequences of the one-spine morphotype (E. caudatum) differ from the five-spine morphotype (E. cattanei) by 12 nucleotide substitutions (data not shown). Identifying and classifying the organisms based on conserved and variable regions of 18S rDNA is a standard procedure in taxonomy studies of eukaryotic microorganisms [47]. 18S rDNA genes are highly conserved intra-species markers among eukaryotes, and finding potential primer sites specific to certain taxa can be challenging [48,49]. The internal transcribed spacer regions (ITS1 and ITS2) that separate nuclear ribosomal genes are considered to be quite variable and have helped discriminate interspecific and intraspecific genetic variation from a wide diversity of life from bacteria to plants and vertebrates. Among rumen ciliates, the ITS1 and ITS2 regions were used to discriminate variations between eight isolates of Isotricha prostoma from two continents (North America and Australia) and different hosts (sheep and cattle). No differences were reported, indicating that the regions are 100% conserved [50]. Other rumen ciliates were not evaluated by this marker. No other markers were used in the evaluation of the phylogenetic relationship of rumen ciliates. The ITS1 and ITS2 were used in the phylogenetic analysis of other free-living litostomatean ciliates [51]. The single- and multilocus analyses of ITS1-5.8S rDNA-ITS2 were in agreement with previous 18S rDNA gene phylogenies, supporting that both 18S rDNA and ITS region sequences are effective tools for resolving phylogenetic relationships among the litostomateans [51]. The single- and multilocus analyses were also used in the phylogenetic study of endocommensal astome ciliates from earthworms. Their systematic positions were determined using the 18S rDNA gene, ITS1-5.8S-ITS2region, and the barcoding D1/D2 domains of the 28S rDNA gene [52].

5. Conclusions

The morphometric and molecular analysis results are inconclusive to decide whether the five-spine Epidinium from the European bison represents one or more species. For now, it can be recommended to name the Epidinium morphotypes according to Dogiel [25]. Thus, the spination will be included in the name; although, it is unclear how this feature relates to the phylogenetic relationship of these morphotypes. Complete sequences from various hosts are necessary to reveal the validity of the Epidinium five-spine morphospecies. Using other potential molecular markers (large ribosomal DNA subunit—LSU rDNA, 28S rDNA) and internal transcribed spacers (ITS1 and ITS2) can produce a more reliable view of the relationships between rumen ciliate morphospecies [1,2,8,9,53,54,55,56].

Author Contributions

Conceptualization, S.K.; Data curation, S.K.; Formal analysis, S.I., A.K., P.P. and S.K.; Funding acquisition, S.K.; Investigation, S.I., A.K., P.P. and S.K.; Methodology, A.K., P.P. and S.K.; Project administration, S.K.; Supervision, S.K.; Validation, P.P. and S.K.; Visualization, P.P. and S.K.; Writing—original draft, S.I. and S.K.; Writing—review and editing, S.I., A.K., P.P. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Grant Scientific Agency of the Ministry of Education, Science, Research and Sport of the Slovak Republic and the Slovak Academy of Sciences (VEGA 2/0009/08 and VEGA 2/0008/21).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and analyzed during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank V. Kolarčik, University of P. J. Šafárik, for his skilled advice in the multivariate analysis of morphometric data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Accession numbers of the partial sequences of the 18S rDNA gene of the rumen commensal ciliate Epidinium parvicaudatum (syn. Epidinium ecaudatum f. parvicaudatum) isolated from the formaldehyde-preserved rumen samples of European bison (Bison bonasus).
Table A1. Accession numbers of the partial sequences of the 18S rDNA gene of the rumen commensal ciliate Epidinium parvicaudatum (syn. Epidinium ecaudatum f. parvicaudatum) isolated from the formaldehyde-preserved rumen samples of European bison (Bison bonasus).
Species NameIsolate LabelGenBank Accession NumberCollection YearCountry of Origin of the HostLocalityHost OrganismSource
E. parvicaudatum (Sharp, 1914) Clone 11OR6333372007PolandBialowieza ForestBison bonasus (Linnaeus, 1758)rumen
E. parvicaudatum (Sharp, 1914)Clone 12OR6333382007PolandBialowieza ForestBison bonasus (Linnaeus, 1758)rumen
E. parvicaudatum (Sharp, 1914)Clone 18OR6333392007PolandBialowieza ForestBison bonasus (Linnaeus, 1758)rumen
E. parvicaudatum (Sharp, 1914)Clone 19OR6333402007PolandBialowieza ForestBison bonasus (Linnaeus, 1758)rumen

References

  1. Warren, A.; Patterson, D.J.; Dunthorn, M.; Clamp, J.C.; Achilles-Day, U.E.; Aescht, E.; Al-Farraj, S.A.; Al-Quraishy, S.; Al-Rasheid, K.; Carr, M.; et al. Beyond the “Code”: A Guide to the Description and Documentation of Biodiversity in Ciliated Protists (Alveolata, Ciliophora). J. Eukaryot. Microbiol. 2017, 64, 539–554. [Google Scholar] [CrossRef] [PubMed]
  2. Zhao, Y.; Yi, Z.; Gentekaki, E.; Zhan, A.; Al-Farraj, S.A.; Song, W. Utility of combining morphological characters, nuclear and mitochondrial genes: An attempt to resolve the conflicts of species identification for ciliated protists. Mol. Phylogenet. Evol. 2016, 94, 718–729. [Google Scholar] [CrossRef] [PubMed]
  3. Jia, T.Z.; Caudan, M.; Mamajanov, I. Origin of Species before Origin of Life: The Role of Speciation in Chemical Evolution. Life 2021, 11, 154. [Google Scholar] [CrossRef] [PubMed]
  4. Mallet, J. Species, Concepts of. In Encyclopedia of Biodiversity; Elsevier: Amsterdam, The Netherlands, 2013; pp. 679–691. [Google Scholar] [CrossRef]
  5. Finlay, B.J.; Corliss, J.O.; Esteban, G.; Fenchel, T. Biodiversity at the Microbial Level: The Number of Free-Living Ciliates in the Biosphere. Q. Rev. Biol. 1996, 71, 221–237. [Google Scholar] [CrossRef]
  6. Rossi, M.F.; Dias, R.J.P.; Senra, M.V.X.; Martinele, I.; Soares, C.A.G.; D’Agosto, M. Molecular Phylogeny of the Family Ophryscolecidae (Ciliophora, Litostomatea) Inferred from 18S rDNA Sequences. J. Eukaryot. Microbiol. 2015, 62, 584–590. [Google Scholar] [CrossRef] [PubMed]
  7. Cedrola, F.; Rossi, M.; Dias, R.J.P.; Martinele, I.; D’Agosto, M. Methods for Taxonomic Studies of Rumen Ciliates (Alveolata: Ciliophora): A Brief Review. Zool. Sci. 2015, 32, 8–15. [Google Scholar] [CrossRef]
  8. Kittelmann, S.; Devente, S.R.; Kirk, M.R.; Seedorf, H.; Dehority, B.A.; Janssen, P.H. Phylogeny of Intestinal Ciliates, Including Charonina ventriculi, and Comparison of Microscopy and 18S rRNA Gene Pyrosequencing for Rumen Ciliate Community Structure Analysis. Appl. Environ. Microbiol. 2015, 81, 2433–2444. [Google Scholar] [CrossRef]
  9. Clamp, J.C.; Lynn, D.H. Investigating the biodiversity of ciliates in the ‘Age of Integration’. Eur. J. Protistol. 2017, 61, 314–322. [Google Scholar] [CrossRef]
  10. West-Eberhard, M.J. Phenotypic Plasticity. Encycl. Ecol. 2008, 2701–2707. [Google Scholar] [CrossRef]
  11. Williams, A.G.; Coleman, G.S. The Rumen Protozoa; Springer: New York, NY, USA, 1992. [Google Scholar] [CrossRef]
  12. Ogimoto, K.; Imai, S. Atlas of Rumen Microbiology; Japan Scientific Societies Press: Tokyo, Japan, 1981. [Google Scholar]
  13. Cedrola, F.; Senra, M.V.X.; Fregulia, P.; D’agosto, M.; Dias, R.J.P. Insights into the systematics of the family Ophryoscolecidae (Ciliophora, Entodiniomorphida). J. Eukaryot. Microbiol. 2022, 69, e12915. [Google Scholar] [CrossRef]
  14. Cedrola, F.; Dias, R.J.P.; Martinele, I.; D’agosto, M. Polymorphism and inconsistencies in the taxonomy of Diplodinium anisacanthum da Cunha, 1914 (Ciliophora, Entodiniomorphida, Ophryoscolecidae) and taxonomic notes on the genus Diplodinium Schuberg, 1888. Zootaxa 2017, 4306, 249–260. [Google Scholar] [CrossRef]
  15. Dehority, B.A. Rumen Ciliate Protozoa of the Blue Duiker (Cephalophus monticola), with Observations on Morphological Variation Lines Within the Species Entodinium dubardi. J. Eukaryot. Microbiol. 1994, 41, 103–111. [Google Scholar] [CrossRef] [PubMed]
  16. Miltko, R.; Wereszka, K.; Kasperowicz, A.; Michalowski, T. Observation on the growth and morphological variability of the rumen ciliates Epidinium caudatum f. fasciculus cultured in vitro. Endocytobiosis Cell Res. 2006, 17, 164–170. [Google Scholar]
  17. Miltko, R.; Michałowski, T.; Pristas, P.; Javorsky, P.; Hackstein, J.H.P. Factors influencing morphological variability of rumen ciliates from the genus Ophryoscolex. J. Anim. Feed Sci. 2006, 15, 35–38. [Google Scholar] [CrossRef]
  18. Dehority, B.A. Studies on the variation in caudal spination of Epidinium. Zootaxa 2006, 1305, 33–39. [Google Scholar] [CrossRef]
  19. Coleman, G.S.; Davies, J.I.; Cash, M.A. The Cultivation of the Rumen Ciliates Epidinium ecaudatum caudatum and Polyplastron multivesiculatum in vitro. J. Gen. Microbiol. 1972, 73, 509–521. [Google Scholar] [CrossRef]
  20. Coleman, G.S.; Laurie, J.I.; Bailey, J.E. The Cultivation of the Rumen Ciliate Entodinium bursa in the Presence of Entodinium caudatum. J. Gen. Microbiol. 1977, 101, 253–258. [Google Scholar] [CrossRef]
  21. Clarke, R.T.J. The Cultivation of Some Rumen Oligotrich Protozoa. J. Gen. Microbiol. 1963, 33, 401–408. [Google Scholar] [CrossRef]
  22. Mah, R.A. Factors Influencing the in vitro Culture of the Rumen Ciliate Ophryoscolex purkynei Stein. J. Protozool. 1964, 11, 546–552. [Google Scholar] [CrossRef]
  23. Imai, S.; Abe, M.; Ogimoto, K. Ciliate protozoa from the rumen of the Japanese serow, Capricornis crispus (Temminck). Jpn. J. Vet. Sci. 1981, 43, 359–367. [Google Scholar] [CrossRef]
  24. Fiorentini, A. Intorno ai Protisti Stomaco dei Bovini; Frat Fusi: Pavia, Italy, 1889. [Google Scholar]
  25. Dogiel, V.A. Monographie der Familie Ophryoscolecidae. Teil I. Arch fur Protistenkd. 1927, 59, 1–289. [Google Scholar]
  26. Sharp, R.C. Diplodinium ecaudatum with an Account of Its Neuromotor Apparatus; University of California Press: Berkeley, CA, USA, 1914; pp. 43–122. [Google Scholar]
  27. Banerjee, A.K. Studies on parasitic ciliates from Indian ruminants. Proc. Zool. Soc. Bengal. 1955, 8, 87–101. [Google Scholar]
  28. Göçmen, B.; Atatür, M.K. Some Rumen Ciliates (Isotrichidae, Trichostomatida; Epidininae, Ophryoscolecidae) of the Domestic Goat (Capra hircus L.) in Turkey. Turkish J. Zool. 2002, 26, 15–26. [Google Scholar]
  29. Kofoid, C.A.; MacLennan, R.F. Ciliates from Bos indicus Linn. III. Epidinium Crawley, Epiplastron Gen. Nov. and Ophryoscolex Stein; University of California Press: Berkeley, CA, USA, 1933. [Google Scholar]
  30. der Staay, S.Y.M.-V.; van der Staay, G.W.; Michalowski, T.; Jouany, J.-P.; Pristas, P.; Javorský, P.; Kišidayová, S.; Varadyova, Z.; McEwan, N.R.; Newbold, C.J.; et al. The symbiotic intestinal ciliates and the evolution of their hosts. Eur. J. Protistol. 2014, 50, 166–173. [Google Scholar] [CrossRef]
  31. Kišidayová, S.; Durkaj, D.; Mihaliková, K.; Váradyová, Z.; Puchalska, J.; Szumacher-Strabel, M.; Cieślak, A.; Gizejewski, Z. Rumen Ciliated Protozoa of the Free-Living European Bison (Bison bonasus, Linnaeus). Front. Microbiol. 2021, 12, 658448. [Google Scholar] [CrossRef] [PubMed]
  32. Wright, A.-D.G.; Dehority, B.A.; Lynn, D.H. Phylogeny of the Rumen Ciliates Entodinium, Epidinium and Polyplastron (Litostomatea: Entodiniomorphida) Inferred from Small Subunit Ribosomal RNA Sequences. J. Eukaryot. Microbiol. 1997, 44, 61–67. [Google Scholar] [CrossRef] [PubMed]
  33. Hammer, D.A.T.; Ryan, P.D.; Hammer, Ø.; Harper, D.A.T. Past: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol. Electron. 2001, 4, 178. [Google Scholar]
  34. Burk, A.D. Laboratory Manual for Classification and Morphology of Rumen Ciliate Protozoa; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar] [CrossRef]
  35. Regensbogenova, M.; Kišidayová, S.; Michalowski, T.; Javorsky, P.; Moon-Van Der Staay, S.Y.; Hackstein, J.H.P.; Mcewan, N.R.; Jouany, J.-P.; Newbold, J.C.; Pristas, P. Rapid identification of rumen protozoa by restriction analysis of amplified 18S rRNA gene. Acta Protozool. 2004, 43, 219–224. [Google Scholar]
  36. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  37. Dehority, B.A. Physiological characteristics of several rumen protozoa grown in vitro with observations on within and among species variation. Eur. J. Protistol. 2010, 46, 271–279. [Google Scholar] [CrossRef]
  38. Eadie, J.M. Inter-Relationships between Certain Rumen Ciliate Protozoa. J. Gen. Microbiol. 1962, 29, 579–588. [Google Scholar] [CrossRef]
  39. Martinele, I.; D’Agosto, M. Predação e canibalismo entre protozoários ciliados (Ciliophora: Entodiniomorphida: Ophryoscolecidae) no rúmen de ovinos (Ovis aries). Rev. Bras. Zool. 2008, 25, 451–455. [Google Scholar] [CrossRef]
  40. Dehority, B.A. Microbial Interactions in the Rumen. Rev. Fac. Agron. (LUZ) 1998, 15, 69–86. Available online: https://www.revfacagronluz.org.ve/v15_1/v151z009.html (accessed on 6 November 2023).
  41. Paireder, S.; Werner, B.; Bailer, J.; Werther, W.; Schmid, E.; Patzak, B.; Cichna-Markl, M. Comparison of protocols for DNA extraction from long-term preserved formalin fixed tissues. Anal. Biochem. 2013, 439, 152–160. [Google Scholar] [CrossRef] [PubMed]
  42. Ottestad, A.L.; Emdal, E.F.; Grønberg, B.H.; Halvorsen, T.O.; Dai, H.Y. Fragmentation assessment of FFPE DNA helps in evaluating NGS library complexity and interpretation of NGS results. Exp. Mol. Pathol. 2022, 126, 104771. [Google Scholar] [CrossRef] [PubMed]
  43. Guyard, A.; Boyez, A.; Pujals, A.; Robe, C.; Van Nhieu, J.T.; Allory, Y.; Moroch, J.; Georges, O.; Fournet, J.-C.; Zafrani, E.-S.; et al. DNA degrades during storage in formalin-fixed and paraffin-embedded tissue blocks. Virchows Arch. 2017, 471, 491–500. [Google Scholar] [CrossRef]
  44. Nam, S.K.; Im, J.; Kwak, Y.; Han, N.; Nam, K.H.; Na Seo, A.; Lee, H.S. Effects of Fixation and Storage of Human Tissue Samples on Nucleic Acid Preservation. Korean J. Pathol. 2014, 48, 36–42. [Google Scholar] [CrossRef]
  45. Robbe, P.; Popitsch, N.; Knight, S.J.; Antoniou, P.; Becq, J.; He, M.; Kanapin, A.; Samsonova, A.; Vavoulis, D.V.; Ross, M.T.; et al. Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: Pilot study for the 100,000 Genomes Project. Genet Med. 2018, 20, 1196–1205. [Google Scholar] [CrossRef]
  46. Lin, J.; Kennedy, S.H.; Svarovsky, T.; Rogers, J.; Kemnitz, J.W.; Xu, A.; Zondervan, K.T. High-quality genomic DNA extraction from formalin-fixed and paraffin-embedded samples deparaffinized using mineral oil. Anal. Biochem. 2009, 395, 265–267. [Google Scholar] [CrossRef]
  47. Wang, Y.; Tian, R.M.; Gao, Z.M.; Bougouffa, S.; Qian, P.-Y. Optimal Eukaryotic 18S and Universal 16S/18S Ribosomal RNA Primers and Their Application in a Study of Symbiosis. PLoS ONE 2014, 9, e90053. [Google Scholar] [CrossRef]
  48. Ishaq, S.L.; Wright, A.-D.G. Design and Validation of Four New Primers for Next-Generation Sequencing To Target the 18S rRNA Genes of Gastrointestinal Ciliate Protozoa. Appl. Environ. Microbiol. 2014, 80, 5515–5521. [Google Scholar] [CrossRef]
  49. Karnati, S.K.R.; Yu, Z.; Sylvester, J.T.; Dehority, B.A.; Morrison, M.; Firkins, J.L. Technical note: Specific PCR amplification of protozoal 18S rDNA sequences from DNA extracted from ruminal samples of cows1. J. Anim. Sci. 2003, 81, 812–815. [Google Scholar] [CrossRef] [PubMed]
  50. Wright, A.-D.G. Analysis of intraspecific sequence variation among eight isolates of the rumen symbiont, Isotricha prostoma (Ciliophora), from two continents. J. Eukaryot. Microbiol. 1999, 46, 445–446. [Google Scholar] [CrossRef] [PubMed]
  51. Vd’ačný, P.; Bourland, W.A.; Orsi, W.; Epstein, S.S.; Foissner, W. Genealogical analyses of multiple loci of litostomatean ciliates (Protista, Ciliophora, Litostomatea). Mol. Phylogenet. Evol. 2012, 65, 397–411. [Google Scholar] [CrossRef] [PubMed]
  52. Obert, T.; Vďačný, P. Multi-gene phylogeny of the subclass Astomatia (Protista: Ciliophora) refreshed with two rare astome ciliates from the digestive tube of endogeic earthworms. Org. Divers. Evol. 2021, 21, 59–77. [Google Scholar] [CrossRef]
  53. Newbold, C.J.; De La Fuente, G.; Belanche, A.; Ramos-Morales, E.; McEwan, N.R. The Role of Ciliate Protozoa in the Rumen. Front. Microbiol. 2015, 6, 1313. [Google Scholar] [CrossRef]
  54. Santoferrara, L.F.; McManus, G.B.; Alder, V.A. Utility of Genetic Markers and Morphology for Species Discrimination within the Order Tintinnida (Ciliophora, Spirotrichea). Protist 2013, 164, 24–36. [Google Scholar] [CrossRef]
  55. Berney, C.; Ciuprina, A.; Bender, S.; Brodie, J.; Edgcomb, V.; Kim, E.; Rajan, J.; Parfrey, L.W.; Adl, S.; Audic, S.; et al. UniEuk: Time to Speak a Common Language in Protistology! J. Eukaryot. Microbiol. 2017, 64, 407–411. [Google Scholar] [CrossRef]
  56. Belanche, A.; Kingston-Smith, A.H.; Griffith, G.W.; Newbold, C.J. A multi-kingdom study reveals the plasticity of the rumen microbiota in response to a shift from non-grazing to grazing diets in sheep. Front. Microbiol. 2019, 10, 122. [Google Scholar] [CrossRef]
Figure 1. Photomicrographs of five-spine morphotypes of Epidinium from the rumen of European bison (Bison bonasus). Samples were colorized with an iodine solution. The protozoan right side (A,D,E,H,K); the protozoan left side (B,C,F,I,J); the protozoan dorsal side (G,L); a dividing ciliate (I). The typical representatives of the E. parvicaudatum morphotype (B,E,F,J). DS, dorsal side; VS, ventral side; sp1–sp5, the numbered sequence of spines; MaN, the macronucleus; AZM, adoral zone of membranels; DZM, dorsal zone of membranels; SkP, skeletal plates; PV, pulsatile vacuoles. The length of the scale bar is 20 µm in all pictures.
Figure 1. Photomicrographs of five-spine morphotypes of Epidinium from the rumen of European bison (Bison bonasus). Samples were colorized with an iodine solution. The protozoan right side (A,D,E,H,K); the protozoan left side (B,C,F,I,J); the protozoan dorsal side (G,L); a dividing ciliate (I). The typical representatives of the E. parvicaudatum morphotype (B,E,F,J). DS, dorsal side; VS, ventral side; sp1–sp5, the numbered sequence of spines; MaN, the macronucleus; AZM, adoral zone of membranels; DZM, dorsal zone of membranels; SkP, skeletal plates; PV, pulsatile vacuoles. The length of the scale bar is 20 µm in all pictures.
Life 13 02350 g001
Figure 2. The clustering of the five-spine morphotypes of Epidinium from the rumen of European bison (Bison bonasus) by K-means analysis.
Figure 2. The clustering of the five-spine morphotypes of Epidinium from the rumen of European bison (Bison bonasus) by K-means analysis.
Life 13 02350 g002
Figure 3. Rooted neighbor-joining tree showing relatedness of Epidinium parvicaudatum (syn. E. ecaudatum f. parvicaudatum) partial 18S rDNA sequences obtained through this study to the related rumen ciliate sequences in GenBank database. Numbers at nodes are bootstrap values after 1000 repetitions. As an outgroup, the 18S rDNA sequence of Entodinium sp. LDK-2011 clone 8 was used.
Figure 3. Rooted neighbor-joining tree showing relatedness of Epidinium parvicaudatum (syn. E. ecaudatum f. parvicaudatum) partial 18S rDNA sequences obtained through this study to the related rumen ciliate sequences in GenBank database. Numbers at nodes are bootstrap values after 1000 repetitions. As an outgroup, the 18S rDNA sequence of Entodinium sp. LDK-2011 clone 8 was used.
Life 13 02350 g003
Table 1. Descriptive statistics of the body length, body width, and the ratio of body length to body width of the five-spine rumen ciliate Epidinium from European bison.
Table 1. Descriptive statistics of the body length, body width, and the ratio of body length to body width of the five-spine rumen ciliate Epidinium from European bison.
EpidiniumLength (µm)Width (µm)L/W
Mean103.157.41.8
Standard Error0.960.600.01
Median101.2956.021.80
Standard Deviation15.719.700.18
Minimum67.4535.931.34
Maximum144.4289.902.36
Count266266266
Table 2. Descriptive statistics of the body length, body width, and the ratio of body length to body width of the five-spine rumen ciliate Epidinium from the European bison after K-means clustering.
Table 2. Descriptive statistics of the body length, body width, and the ratio of body length to body width of the five-spine rumen ciliate Epidinium from the European bison after K-means clustering.
Cluster 1Cluster 2
Length
(µm)
Width
(µm)
L/WLength
(µm)
Width
(µm)
L/W
Mean92.251.71.8118.765.81.8
Standard Error0.730.490.011.010.760.02
Median93.7351.181.80116.1465.161.83
St. Deviation9.396.300.1710.698.050.20
Minimum58.733.51.4101.2450.001.34
Maximum109.968.92.3157.4989.902.36
Count166166166112112112
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ivorová, S.; Kopčaková, A.; Pristaš, P.; Kišidayová, S. Morphometric and Molecular Analysis of Five-Spine Epidinium Morphotypes Taken from the Rumen of European Bison, Bison bonasus. Life 2023, 13, 2350. https://doi.org/10.3390/life13122350

AMA Style

Ivorová S, Kopčaková A, Pristaš P, Kišidayová S. Morphometric and Molecular Analysis of Five-Spine Epidinium Morphotypes Taken from the Rumen of European Bison, Bison bonasus. Life. 2023; 13(12):2350. https://doi.org/10.3390/life13122350

Chicago/Turabian Style

Ivorová, Silvia, Anna Kopčaková, Peter Pristaš, and Svetlana Kišidayová. 2023. "Morphometric and Molecular Analysis of Five-Spine Epidinium Morphotypes Taken from the Rumen of European Bison, Bison bonasus" Life 13, no. 12: 2350. https://doi.org/10.3390/life13122350

APA Style

Ivorová, S., Kopčaková, A., Pristaš, P., & Kišidayová, S. (2023). Morphometric and Molecular Analysis of Five-Spine Epidinium Morphotypes Taken from the Rumen of European Bison, Bison bonasus. Life, 13(12), 2350. https://doi.org/10.3390/life13122350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop