Divergent marine anaerobic ciliates harbor closely related Methanocorpusculum endosymbionts.

Ciliates are a diverse group of protists known for their ability to establish various partnerships and thrive in a wide variety of oxygen-depleted environments. Most anaerobic ciliates harbor methanogens, one of the few known archaea living intracellularly. These methanogens increase the metabolic efficiency of host fermentation via syntrophic use of host end-product in methanogenesis. Despite the ubiquity of these symbioses in anoxic habitats, patterns of symbiont specificity and fidelity are not well known. We surveyed two unrelated, commonly found groups of anaerobic ciliates, the Plagiopylea and Metopida, isolated from anoxic marine sediments. We sequenced host 18S rRNA and symbiont 16S rRNA marker genes as well as the symbiont internal transcribed spacer region from our cultured ciliates to identify hosts and their associated methanogenic symbionts. We found that marine ciliates from both of these co-occurring, divergent groups harbor closely related yet distinct intracellular archaea within the Methanocorpusculum genus. The symbionts appear to be stable at the host species level, but at higher taxonomic levels, there is evidence that symbiont replacements have occurred. Gaining insight into this unique association will deepen our understanding of the complex transmission modes of marine microbial symbionts, and the mutualistic microbial interactions occurring across domains of life.

Role of turbulence in Dinophysis spp. growth, feeding, and toxin leakage in culture.

Dinophysis, a mixotrophic dinoflagellate that is known to prey on the ciliate Mesodinium rubrum, and retain its chloroplasts, is responsible for diarrhetic shellfish poisoning (DSP) in humans and has been identified on all U.S. coasts. Monocultures of Dinophysis have been used to investigate the growth of Dinophysis species in response to variations in environmental conditions, however, little is known about the roles of system stability (turbulence) and mixotrophy in the growth and toxicity of Dinophysis species in the U.S.. To begin to address this gap in knowledge, culturing experiments were conducted with three species (four strains) of Dinophysis, that included predator-prey co-incubation (Dinophysis spp.+ M. rubrum) and prey-only (M. rubrum) flasks. Cultures were investigated for effects of low or high turbulence on Dinophysis spp. growth, feeding, and amounts of intra- and extracellular toxins: okadaic acid and derivatives (diarrhetic shellfish toxins, DSTs) and pectenotoxins (PTXs). Turbulence did not have a measurable effect on the rates of ingestion of M. rubrum prey by Dinophysis spp. for any of the four strains, however, effects on growth and particulate and dissolved toxins were observed. High turbulence (ε = 10-2 m2s-3) significantly slowed growth of both D. acuminata and D. ovum relative to still controls, but significantly stimulated growth of the D. caudata strain. Increasing turbulence also resulted in significantly higher intracellular toxin content in D. acuminata cultures (DSTs and PTXs), but significantly reduced intracellular toxin content (PTXs) in those of D. caudata. An increase in turbulence appeared to promote toxin leakage, as D. ovum had significantly more extracellular DSTs found in the medium under high turbulence when compared to the still control. Overall, significant responses to turbulence were observed, whereby the three strains from the "Dinophysis acuminata complex" displayed a stress response to turbulence, i.e., decreasing growth, increasing intracellular toxin content and/or increasing toxin leakage, while the D. caudata strain had an opposite response, appearing stimulated by, or more tolerant of, high turbulence.

Numerical and Thermal Response of the Bacterivorous Ciliate Colpidium kleini, a Species Potentially at Risk of Extinction by Rising Water Temperatures.

We investigated the food-dependent growth and thermal response of the freshwater ciliate Colpidium kleini using numerical response (NR) experiments. This bacterivorous ciliate occurs in lotic water and the pelagial of lakes and ponds. The C. kleini strain used in this work was isolated from a small alpine lake and identified by combining detailed morphological inspections with molecular phylogeny. Specific growth rates (rmax) were measured from 5 to 21 °C. The ciliate did not survive at 22 °C. The threshold bacterial food levels (0.3 - 2.2 × 106 bacterial cells mL-1) matched the bacterial abundance in the alpine lake from which C. kleini was isolated. The food threshold was notably lower than previously reported for C. kleini and two other Colpidium species. The threshold was similar to levels reported for oligotrich and choreotrich ciliates if expressed in terms of bacterial biomass (0.05 - 0.43 mg C L-1). From the NR results, we calculated physiological mortality rates at zero food concentration. The mean mortality (0.55 ± 0.17 d-1) of C. kleini was close to the mean estimate obtained for other planktonic ciliates that do not encyst. We used the data obtained by the NR experiments to fit a thermal performance curve (TPC). The TPC yielded a temperature optimum at 17.3 °C for C. kleini, a maximum upper thermal tolerance limit of 21.9 °C, and a thermal safety margin of 4.6 °C. We demonstrated that combining NR with TPC analysis is a powerful tool to predict better a species' fitness in response to temperature and food.

A tiny, long-distance hunter.

Lacrymaria olor cytoskeleton and membrane "origami" enables rapid cell hyperextension.

Curved crease origami and topological singularities enable hyperextensibility of L. olor.

Fundamental limits of cellular deformations, such as hyperextension of a living cell, remain poorly understood. Here, we describe how the single-celled protist Lacrymaria olor, a 40-micrometer cell, is capable of reversibly and repeatably extending its necklike protrusion up to 1200 micrometers in 30 seconds. We discovered a layered cortical cytoskeleton and membrane architecture that enables hyperextensions through the folding and unfolding of cellular-scale origami. Physical models of this curved crease origami display topological singularities, including traveling developable cones and cytoskeletal twisted domain walls, which provide geometric control of hyperextension. Our work unravels how cell geometry encodes behavior in single cells and provides inspiration for geometric control in microrobotics and deployable architectures.

Cellular interactions and evolutionary origins of endosymbiotic relationships with ciliates.

As unicellular predators, ciliates engage in close associations with diverse microbes, laying the foundation for the establishment of endosymbiosis. Originally heterotrophic, ciliates demonstrate the ability to acquire phototrophy by phagocytizing unicellular algae or by sequestering algal plastids. This adaptation enables them to gain photosynthate and develop resistance to unfavorable environmental conditions. The integration of acquired phototrophy with intrinsic phagotrophy results in a trophic mode known as mixotrophy. Additionally, ciliates can harbor thousands of bacteria in various intracellular regions, including the cytoplasm and nucleus, exhibiting species specificity. Under prolonged and specific selective pressure within hosts, bacterial endosymbionts evolve unique lifestyles and undergo particular reductions in metabolic activities. Investigating the research advancements in various endosymbiotic cases within ciliates will contribute to elucidate patterns in cellular interaction and unravel the evolutionary origins of complex traits.

Genomic insights into the cellular specialization of predation in raptorial protists.

BACKGROUND: Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited. RESULTS: In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria. CONCLUSIONS: Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

Exploring the landscape of symbiotic diversity and distribution in unicellular ciliated protists.

BACKGROUND: The eukaryotic-bacterial symbiotic system plays an important role in various physiological, developmental, and evolutionary processes. However, our current understanding is largely limited to multicellular eukaryotes without adequate consideration of diverse unicellular protists, including ciliates. RESULTS: To investigate the bacterial profiles associated with unicellular organisms, we collected 246 ciliate samples spanning the entire Ciliophora phylum and conducted single-cell based metagenome sequencing. This effort has yielded the most extensive collection of bacteria linked to unicellular protists to date. From this dataset, we identified 883 bacterial species capable of cohabiting with ciliates, unveiling the genomes of 116 novel bacterial cohabitants along with 7 novel archaeal cohabitants. Highlighting the intimate relationship between ciliates and their cohabitants, our study unveiled that over 90% of ciliates coexist with bacteria, with individual hosts fostering symbiotic relationships with multiple bacteria concurrently, resulting in the observation of seven distinct symbiotic patterns among bacteria. Our exploration of symbiotic mechanisms revealed the impact of host digestion on the intracellular diversity of cohabitants. Additionally, we identified the presence of eukaryotic-like proteins in bacteria as a potential contributing factor to their resistance against host digestion, thereby expanding their potential host range. CONCLUSIONS: As the first large-scale analysis of prokaryotic associations with ciliate protists, this study provides a valuable resource for future research on eukaryotic-bacterial symbioses. Video Abstract.

Theory of Stoichiometric Intraguild Predation: Algae, Ciliate, and Daphnia.

Consumers respond differently to external nutrient changes than producers, resulting in a mismatch in elemental composition between them and potentially having a significant impact on their interactions. To explore the responses of herbivores and omnivores to changes in elemental composition in producers, we develop a novel stoichiometric model with an intraguild predation structure. The model is validated using experimental data, and the results show that our model can well capture the growth dynamics of these three species. Theoretical and numerical analyses reveal that the model exhibits complex dynamics, including chaotic-like oscillations and multiple types of bifurcations, and undergoes long transients and regime shifts. Under moderate light intensity and phosphate concentration, these three species can coexist. However, when the light intensity is high or the phosphate concentration is low, the energy enrichment paradox occurs, leading to the extinction of ciliate and Daphnia. Furthermore, if phosphate is sufficient, the competitive effect of ciliate and Daphnia on algae will be dominant, leading to competitive exclusion. Notably, when the phosphorus-to-carbon ratio of ciliate is in a suitable range, the energy enrichment paradox can be avoided, thus promoting the coexistence of species. These findings contribute to a deeper understanding of species coexistence and biodiversity.

Reductitherus cryptostomus n. gen., n. sp. (Ciliophora: Armophorea: Clevelandellida), a remarkable new nyctotherid from an Australian cockroach, Parapanesthia gigantea (Blaberidae: Panesthiinae).

Ciliates of the family Nyctotheridae (Armophorea: Clevelandellida) are frequent intestinal symbionts of various invertebrates and some poikilotherm vertebrates. Depending on the classification scheme, there are between 15 and 18 recognized genera of Nyctotheridae, the majority of which exhibit a rather uniform morphology. They have round to ellipsoidal cells with an adoral zone of membranelles that begins anteriorly in an adoral groove and continues posteriorly into the buccal cavity where it extends deep into the cell in the peristomial funnel. The taxonomy of the Nyctotheridae is primarily based on the number and location of kinetal sutures. The only known divergence from the relatively conservative nyctortherid body plan are the bizarre symbionts of Panesthiinae cockroaches, ciliates of the family Clevelandellidae, which forms a clade nested within the Nyctotheridae genus Nyctotherus. In this study we report another ciliate that diverges morphologically from the canonical Nyctotheridae body plan, and which is also found in Panesthiinae hosts. The novel ciliate Reductitherus cryptostomus n. gen., n. sp. differs from the rest of Nyctotheridae by absence of the anterior adoral groove, a shortened adoral zone completely enclosed in a notably small buccal cavity, and two strongly reduced kinetal sutures, one left anterodorsal and the other right posterodorsal.

Continuous warming drives the colonization dynamics of periphytic ciliate fauna in marine environments.

In order to evaluate the influence of global warming on the ecosystem processes in marine environments, the changes in colonization dynamics of periphytic microbiota were studied using the periphytic ciliate communities as the test organism fauna under a continuous warming gradient of 22℃ (control), 25℃, 28℃, 31℃, and 34 ℃. The results demonstrated that (1) the test ciliate communities generally showed a similar temporal pattern in within the colonization process under the water temperatures from 22 up to 28℃; however, (2) the colonization dynamics were significantly changed, and the fitness of colonization curves to the MacArthur-Wilson model equation was failed under the temperature increased by 6 ℃, and (3) the loading or assimilative capacity of the test aquatic ecosystem was decreased with the increase of water temperature. Therefore, this study suggests that continuous warming may significantly drive the colonization dynamics of periphytic ciliates in marine ecosystems.

Aggregation and proliferation of B cells and T cells in MALTs upon Cryptocaryon irritans infection in large yellow croaker Larimichthys crocea.

Mucosal immunity in mucosa-associated lymphoid tissues (MALTs) plays crucial roles in resisting infection by pathogens, including parasites, bacteria and viruses. However, the mucosal immune response in the MALTs of large yellow croaker (Larimichthys crocea) upon parasitic infection remains largely unknown. In this study, we investigated the role of B cells and T cells in the MALTs of large yellow croaker following Cryptocaryon irritans infection. Upon C. irritans infection, the total IgM and IgT antibody levels were significantly increased in the skin mucus and gill mucus. Notably, parasite-specific IgM antibody level was increased in the serum, skin and gill mucus following parasitic infection, while the level of parasite-specific IgT antibody was exclusively increased in MALTs. Moreover, parasitic infection induced both local and systemic aggregation and proliferation of IgM+ B cells, suggesting that the increased levels of IgM in mucus may be derived from both systemic and mucosal immune tissues. In addition, we observed significant aggregation and proliferation of T cells in the gill, head kidney and spleen, suggesting that T cells may also be involved in the systemic and mucosal immune responses upon parasitic infection. Overall, our findings provided further insights into the role of immunoglobulins against pathogenic infection, and the simultaneous aggregation and proliferation of both B cells and T cells at mucosal surfaces suggested potential interactions between these two major lymphocyte populations during parasitic infection.

SMORES: a simple microfluidic operating room for the examination and surgery of Stentor coeruleus.

Ciliates are powerful unicellular model organisms that have been used to elucidate fundamental biological processes. However, the high motility of ciliates presents a major challenge in studies using live-cell microscopy and microsurgery. While various immobilization methods have been developed, they are physiologically disruptive to the cell and incompatible with microscopy and/or microsurgery. Here, we describe a Simple Microfluidic Operating Room for the Examination and Surgery of Stentor coeruleus (SMORES). SMORES uses Quake valve-based microfluidics to trap, compress, and perform surgery on Stentor as our model ciliate. Compared with previous methods, immobilization by physical compression in SMORES is more effective and uniform. The mean velocity of compressed cells is 24 times less than that of uncompressed cells. The compression is minimally disruptive to the cell and is easily applied or removed using a 3D-printed pressure rig. We demonstrate cell immobilization for up to 2 h without sacrificing cell viability. SMORES is compatible with confocal microscopy and is capable of media exchange for pharmacokinetic studies. Finally, the modular design of SMORES allows laser ablation or mechanical dissection of a cell into many cell fragments at once. These capabilities are expected to enable biological studies previously impossible in ciliates and other motile species.

Morphological reconstruction during cell regeneration in the ciliate Spirostomum ambiguum.

When the ciliate Spirostomum ambiguum is transected into two pieces, both fragments regenerate and proliferate. In the anterior fragments, which have lost their contractile vacuoles due to transection, new contractile vacuoles were formed at their posterior ends in a few minutes. When the cells were cut into three pieces, new contractile vacuoles were formed in the anterior and middle fragments, both at their posterior ends. Thus, the anterior-posterior axis of S. ambiguum was maintained after transection. Morphological repair, including the formation of the contractile vacuole, was also observed when only the anteriormost portion was transected to cut out a small fragment that did not contain part of the macronucleus. Scanning electron microscopy was performed to observe changes in the shape of the cleavage surface of S. ambiguum during the wound healing process. Within minutes after cutting, the cut surface was covered with a cilia-free membrane, preventing leakage of cytoplasmic contents. The surface of the cut area then rounded with time and was covered with cilia, completing the repair of the cut area in about one day.

Diatom bloom trigger notable variations in microzooplanktonic ciliate composition, body-size spectrum and biotic-abiotic interaction in the Arctic Ocean.

How microzooplanktonic ciliate adaptative strategies differ across diatom bloom and non-diatom bloom areas in the Arctic Ocean remains poorly documented. To address this gap, two different situations were categorized in the Arctic Ocean at summer 2023: diatom bloom stations (DBS) (genus Thalassiosira, chain-like) and non-diatom bloom stations (nDBS). Total abundance of ciliate at 3 m and 25 m in DBS was 2.8 and 1.8 folds higher than in nDBS, respectively. Aloricate ciliates were singled out in both DBS and nDBS, whilst their average abundance and biomass of large size-fraction (>50 µm) in former were 4.5-5.6 folds higher than in latter. Regarding tintinnids, high abundance of Ptychocylis acuta (Bering Strait species) mainly occurred at DBS, coupled with distribution of co-occurring Pacific-origin species Salpingella sp.1, collectively suggested a strong intrusion of Pacific Inflow during summer 2023. Additionally, presence of high abundance of Acanthostomella norvegica and genus Parafavella in nDBS might indicate the trajectory of the Transpolar Drift. Alternatively, tintinnids can serve as credible bioindicators for either monitoring currents or evaluating microzooplankton Borealization. Average abundance of total ciliate within 15-135 µm body-size spectrum in DBS was higher than nDBS. Moreover, spearman's rank correlation between biotic and abiotic analysis revealed that temperature and dissolved oxygen at DBS determined tintinnid species richness and ciliate total abundance, respectively. The results clearly demonstrate that remarkable divergences in large size-fraction of ciliate abundance between DBS and nDBS validate their irreplaceable role in controlling phytoplankton outbreak and associated biological processes in polar seas.

Coinfection of slime feather duster worms (Annelida, Myxicola) by different gregarine apicomplexans (Selenidium) and astome ciliates reflects spatial niche partitioning and host specificity.

Individual organisms can host multiple species of parasites (or symbionts), and one species of parasite can infect different host species, creating complex interactions among multiple hosts and parasites. When multiple parasite species coexist in a host, they may compete or use strategies, such as spatial niche partitioning, to reduce competition. Here, we present a host­symbiont system with two species of Selenidium (Apicomplexa, Gregarinida) and one species of astome ciliate co-infecting two different species of slime feather duster worms (Annelida, Sabellidae, Myxicola) living in neighbouring habitats. We examined the morphology of the endosymbionts with light and scanning electron microscopy (SEM) and inferred their phylogenetic interrelationships using small subunit (SSU) rDNA sequences. In the host 'Myxicola sp. Quadra', we found two distinct species of Selenidium; S. cf. mesnili exclusively inhabited the foregut, and S. elongatum n. sp. inhabited the mid to hindgut, reflecting spatial niche partitioning. Selenidium elongatum n. sp. was also present in the host M. aesthetica, which harboured the astome ciliate Pennarella elegantia n. gen. et sp. Selenidium cf. mesnili and P. elegantia n. gen. et sp. were absent in the other host species, indicating host specificity. This system offers an intriguing opportunity to explore diverse aspects of host­endosymbiont interactions and competition among endosymbionts.

Microplastics drive community dynamics of periphytic protozoan fauna in marine environments.

The pollution of microplastics (MPs) to the marine environment has become a widespread focus of attention. To assess MP-induced ecotoxicity on marine ecosystems, periphytic protozoan communities were used as test organisms and exposed to five concentrations of MPs: 0, 1, 5, 25, and 125 mg l-1. Protozoan samples were collected using microscope slides from coastal waters of the Yellow Sea, northern China. A total of 13 protozoan species were identified and represented different tolerance to MP-induced ecotoxicity. Inhibition effects of MPs on the test protozoan communities were clearly shown in terms of both the species richness and individual abundance and followed linear relationships to MP concentrations. The community patterns were driven by MPs and significantly shifted at concentrations over 5 mg l-1. Our findings demonstrated that MPs may induce the community-level ecotoxic response of periphytic protozoan fauna and followed significant community dynamics. Thus, it is suggested that periphytic protozoan fauna may be used as useful community-based test model organisms for evaluating MP-induced ecotoxicity in marine environments.

Transcriptomic analysis of large yellow croaker (Larimichthys crocea) reveals the suppression of the inflammatory response from Cryptocaryon irritans infection.

Large yellow croaker (Larimichthys crocea) is the most productive marine fish in China. Cryptocaryon irritans is an extremely destructive parasite that causes great economic losses in large yellow croaker aquaculture industry. Therefore, it is very necessary to study the immune response of large yellow croaker in response to C. irritans infection. In this study, the transcriptomic profiles of large yellow croaker were sequenced and analyzed in the brain and head kidney at 72 h after C. irritans infection. Cytokines and chemokines related terms were significantly enriched based on the GO enrichment of down-regulated differentially expressed genes (DEGs) from the head kidney. Meanwhile, cytokine-cytokine receptor interaction was significantly enriched based on the KEGG enrichment of up-regulated DEGs from the brain and down-regulated DEGs from the head kidney, respectively. Moreover, the majority of inflammation-related DEGs were significantly up-regulated in the brain, but distinctly down-regulated in the head kidney. These results showed that the brain and head kidney might play different roles against C. irritans infection, and the inflammatory response of large yellow croaker may be restrained during C. irritans infection. Taken together, the transcriptomic analyses will be helpful to more comprehensively understand the immune mechanism of teleost against C. irritans infection, and provide a theoretical basis for the prevention and treatment of Cryptosporidiosis.

Does prey availability influence the detection of Dinophysis spp. by the imaging FlowCytobot?

The Imaging FlowCytobot (IFCB) is a field-deployable imaging-in-flow cytometer that is increasingly being used to monitor harmful algae. The IFCB acquires images of suspended particles based on their chlorophyll-a fluorescence and/or the amount of light they scatter (side scattering). The present study hypothesized that fluorescence-based image acquisition would undercount Dinophysis spp., a genus of non-constitutive mixotrophs, when prey is limited. This is because Dinophysis spp. acquire plastids via ingestion of their ciliate prey Mesodinium spp., and lose photosynthetic capacity and autofluorescence in the absence of prey. Even small blooms of Dinophysis spp. can be highly toxic and result in diarrhetic shellfish poisoning (DSP), highlighting the importance of accurately detecting low abundances. To explore this, laboratory experiments were conducted to determine optimal IFCB settings for a fed culture of Dinophysis acuminata, and an existing time series of IFCB observations collected in Puget Sound (Washington, U.S.A) was used to compare Dinophysis spp. abundance estimates from samples triggered via side scattering versus fluorescence in relation to Mesodinium spp. abundance. This study introduces a quantitative approach for optimizing the detection of target harmful algae which can be repeated across multiple IFCBs and demonstrates the effects of IFCB calibration on Dinophysis spp. detection. The laboratory experiments showed that IFCB settings for fluorescence-based image acquisition need to be fairly sensitive to accurately detect D. acuminata cells. A poorly calibrated IFCB can miss a significant proportion of D. acuminata abundance whatever the method used to trigger the image acquisition. Field results demonstrated that the physiological status of Dinophysis spp. can influence their detection by the IFCB when triggering on fluorescence. This was observed during a 7-day period when the IFCB failed to detect Dinophysis spp. cells when triggering on fluorescence while cells were still detected using the side scattering triggering method as well as observed by microscopy. During this period, Mesodinium spp. was not detected, IFCB-derived autofluorescence level of individual cells of Dinophysis spp. was low, and less than 50 % of Dinophysis spp. cells exhibited autofluorescence under the microscope. Together, this indicates that the unique feeding ecology of Dinophysis spp. may affect their detection by the IFCB when cells are starved.

Transcriptomic Variation of Amphiprion Percula (Lacepède, 1802) in Response to Infection with Cryptocaryon Irritans Brown, 1951.

Cryptocaryon irritans (Brown 1951) frequently infect the Pomacentridae fishes causing severe economic losses. However, the anti-C. irritans' molecular mechanism in these fishes remains largely unknown. To address this issue, we conducted RNA-Seq for C. irrtians-infected gills of the clownfish Amphiprion percula (Lacepède 1802) at the early (day 1) and late (day 3) stages of infection. A total of 1655 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs showed a vast genetic variation related to the following aspects: ECM-receptor interaction, P13K-Akt signalling, cytokine-cytokine receptor interaction, and endocytosis. During the early phase of infection, key genes involved in ATP production, energy homeostasis, and stress control were abruptly increased. In the late phase, however, acute response molecules of the peripheral nervous system (synaptic transmission and local immunity), metabolic system triggering glycogen synthesis, energy maintenance, and osmoregulation were found to be critical. The highest number of upregulated genes (URGs) recovered during the early phase was included under the 'biological process' category, which primarily functions as response to stimuli, signalling, and biological regulation. In the late phase, most of the URGs were related to gene regulation and immune system processes under 'molecular function' category. The immune-related URGs of early infection include major histocompatibility complex (MHC) class-II molecules apparently triggering CD4+ T-cell-activated Th responses, and that of late infection include MHC class-1 molecules for the possible culmination of CD8+ T-cell triggered cytotoxicity. The high level of genic single nucleotide polymorphisms (SNPs) identified during the late phase of infection is likely to influence their susceptibility to secondary infection. In summary, the identified DEGs and their related metabolic and immune-related pathways and the SNPs may provide new insights into coordinating the immunological events and improving resistance in Pomacentridae fishes against C. irritans.