*2.1. Data*

We used daily data from the Bitcoin Energy Consumption Index (BECI) over the period 25 February 2017 to 25 January 2022, according to the data availability from Digiconomist. BECI data have been recently used in academia [12,17]. They consist of daily data covering three series, BECI upper bound (BECI UB), BECI lower bound (BECI LB), and BECI average. Accordingly, in this work, 1825 daily observations were used for each of the three series. The plots of the three series are shown in Figure 1, in which an increase in the three series is observed from around the second quarter of 2021, which coincides with the spike in the price of Bitcoin.

**Figure 1.** Plots of BECI UB, LB, and average during the study period (25 February 2017 to 25 January 2022).

Considering both BECI upper bound (BECI UB) and BECI lower bound (BECI LB) for electricity consumptions, we argue the following. The BECI UB is defined as the breakeven point of mining revenues and electricity costs; therefore, it is more sensitive to the economic parameters. In contrast, the BECI LB is a state where all miners use the most efficient hardware, which makes it more stable and reliable for our current study [13,14] examining the long memory traits of the energy consumption series. Therefore, our analysis emphasises the lower bound results.

Our analysis involved the application of 17 windows to the three series, with the length of each sliding window being 200 days. Given that each series consists of 3400 (17\*200) daily observations, a total of 51(17\*3) windows (with 10,200 total observations, 3400\*3) were considered in our calculation for the sample period 25 February 2017 to 25 January 2022. The choice of sliding window-based estimation procedure is backed by the academic literature [18], which points to the suitability of the application of an increasing window size in a dynamic model for estimating long memory. A sliding window approach to modelling is suitable and an increasing estimation window leads to an increase in the estimation accuracy when calibrating the long memory [19].

Each series observation in the BECI indices is expressed in terawatt-hours (TWh), the standard unit of electricity consumption. The three BECI series are significantly meanreverting at first difference. Therefore, all the calculations and analyses were conducted at first difference (Δ). Long-range dependence or long memory was characterised by a slow, power law decay of the autocorrelation function (ACF).


**3. Methods and procols to study haematopoitic development during es cell differentiation**

ES cell differentiation provides a relatively easy and accessible system to study early embry‐ onic haematopoiesis. Using well-defined protocols, it is possible to effectively study the events happening *in vivo* using this *in vitro* approach. This experimental system was shown to recapitulate the early *in vivo* events of development of the haematopoietic system.

#### **3.1. Mouse embryonic fibroblasts (MEFs)**

There are several methods to keep ES cells undifferentiated. One of them consists of growing ES cells on mouse embryonic fibroblasts (MEFs). Before working with ES cells, it is recom‐ mended to prepare a good stock of MEFs to be used as feeder cell layer. For that, wild type ICR or DR4 (resistant to four drugs) [75] MEFs are harvested from E14.5 embryos and cultured in Iscove's modified Dulbecco medium (IMDM, Lonza) supplemented with 50 µg/ml penicil‐ lin-streptomicin (Gibco), 2mM L-Glutamine, 10% of FCS (PAA Laboratories) and 1,5x10-4 monothioglycerol (MTG, Sigma) under low oxygen conditions. When amplified, MEFs are harvested (TrypLE, Invitrogen) and irradiated at 30Gy to stop the cells proliferation. The cells should be frozen at around 1 million cells per ml of IMDM supplemented with 50% FCS and 10% of dimethyl sulfoxide (DMSO). Cells should be stored at-80 0 C. Thawed MEFs should be replated in one six-well plate previously coated with gelatine and let to adhere to the plastic wells overnight.  ES cell culture**

**Figure 2. Molecular regulation of early embryonic haematopoiesis.** The stages of blood development where the

**3. Methods and procols to study haematopoitic development during es cell**

ES cell differentiation provides a relatively easy and accessible system to study early embry‐ onic haematopoiesis. Using well-defined protocols, it is possible to effectively study the events happening *in vivo* using this *in vitro* approach. This experimental system was shown to

There are several methods to keep ES cells undifferentiated. One of them consists of growing ES cells on mouse embryonic fibroblasts (MEFs). Before working with ES cells, it is recom‐ mended to prepare a good stock of MEFs to be used as feeder cell layer. For that, wild type ICR or DR4 (resistant to four drugs) [75] MEFs are harvested from E14.5 embryos and cultured in Iscove's modified Dulbecco medium (IMDM, Lonza) supplemented with 50 µg/ml penicil‐

recapitulate the early *in vivo* events of development of the haematopoietic system.

function of the different genes is critical are indicated.

72 Pluripotent Stem Cell Biology - Advances in Mechanisms, Methods and Models

**3.1. Mouse embryonic fibroblasts (MEFs)**

**differentiation**

ES cells are cultured on irradiated MEFs in a media constituted of Dulbeco's modified Eagle Medium (DMEM, Gibco) supplemented with 50 µg/ml penicillin-streptomycin, 2mM L-Glutamine, 15% FCS (PAA Laboratories), 2% Leukaemia Inhibitory Factor (LIF) (conditioned medium from LIF-generating cell line, see [76]) or 50 units of recombinant ESGRO LIF/ml (Millipore) and 1,5x104 MTG (Sigma). Leukaemia inhibitory factor (LIF) – is a cytokine inhibiting differentiation. ES cells, when cultured on MEFs feeder cell layer in the presence of LIF remain undifferentiated. Upon microscopic observation they form tightly associated clusters of cells that are bright and shiny in appearance (Fig. 3A).

#### **3.3. Generation of embryoid bodies**

Embryoid bodies (EBs) are three-dimensional structures spontaneously generated by ES cells during differentiation. They contain precursors for the three primary germ layers ectoderm, endoderm and mesoderm. Two passages on gelatine are performed to remove the MEFs that would hamper ES cells differentiation. The first passage is performed in DMEM-ES media (described above), whereas for the second passage DMEM is replaced with IMDM. The ES cells are then harvested by trypsinisation and seeded into liquid cultures in non-tissue culture Petri dishes (Sterilin) in differentiation medium containing: IMDM supplemented with 15% FCS serum selected for differentiation (PAA Laboratories), 2mM L-Glutamine, 180 µg/ml transferring (Roche), 25 µg/ml Ascorbic Acid (AA, Sigma) and 4,6x10-4 MTG. The density of cell seeding should be adjusted in function of the day at which the cultures will be harvested, varying from 1,5x104 cells/ml (for day 4-6) up to 3,0 x 104 cells/ml (for days 2.5-3.5). 10-20ml of "Differentiation medium" should be used per one Petri dish.

By performing two passages on gelatine and removing feeder cell layer and LIF, ES cells become primed for differentiation and formation of three-dimensional embryoid bodies in liquid culture (Figure 3B). Early EBs contain precursors for the three primary germ layers. By day 7, hemoglobinisation can be observed as red areas present within the EBs (Figure 3C). This system is versatile and allows to access and study in details several subsequent stages of blood development such as the emergence of haemangioblast, production of blast colonies and the development of primitive and definitive blood precursors (Figure 4).

**Figure 4. The ES/EB differentiation system.** Upon differentiation of ES cells, three-dimensional structures, called em‐ bryoid bodies (EBs) are formed. Sorting EBs at day 2.5-3.5 for the expression of FLK1 enrich for BL-CFC, the *in vitro* equivalent of the haemangioblast. Upon culture, the BL-CFCs generate blast colonies that contain precursors for hae‐ matopoietic, endothelial and vascular smooth muscle cells. Day 4-6 EBs contain haematopoietic progenitors that can

Embryonic Stem Cell Differentiation – A Model System to Study Embryonic Haematopoiesis

http://dx.doi.org/10.5772/57607

75

There are two alternative approaches to study the development of blast colonies; liquid culture on gelatine or semi-solid culture in methylcellulose. Liquid culture on gelatine facilitates harvesting of the cells for flow cytometry analysis or time-lapse imaging techniques. Alterna‐ tively, the semi-solid culture is a clonogenic assay allowing the growth of of individual blast colonies and their quantification. Individual blast colonies can then further be isolated and their haematopoietic potential assessed. To perform blast colony culture, EBs are harvested between day 2.5 and 3.5 and the cells are enriched for haemangioblast by sorting for the expression of FLK1, the receptor 2 for VEGF. Isolated FLK1 positive cells are then replated on gelatinised cell culture plates at the density of 8,5 x 104 cell/10cm2 in a medium containing IMDM supplemented with 10% FCS, 2mM L-Glutamine mM, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 15% endothelial cell line-D4T conditioned medium [11], 10

give rise to various primitive and definitive haematopoietic colonies.

**3.4. Blast colonies**

**Figure 3. Morphology of ES cells and embryoid bodies (EBs).** A) Typical appearance of ES cell in culture. Cells form bright, shiny adherent colonies and are cultured on a layer of MEFs feeder cell layer. B) Typical appearance of EBs in culture. These three-dimensional structures are formed during ES cell differentiation and contain precursors for three primary germ cell layers. C) Embryoid bodies at day 7 of differentiation containing haemoglobin (indicated with aster‐ isks). Scale bar 300µm.

Embryonic Stem Cell Differentiation – A Model System to Study Embryonic Haematopoiesis http://dx.doi.org/10.5772/57607 75

**Figure 4. The ES/EB differentiation system.** Upon differentiation of ES cells, three-dimensional structures, called em‐ bryoid bodies (EBs) are formed. Sorting EBs at day 2.5-3.5 for the expression of FLK1 enrich for BL-CFC, the *in vitro* equivalent of the haemangioblast. Upon culture, the BL-CFCs generate blast colonies that contain precursors for hae‐ matopoietic, endothelial and vascular smooth muscle cells. Day 4-6 EBs contain haematopoietic progenitors that can give rise to various primitive and definitive haematopoietic colonies.

#### **3.4. Blast colonies**

**Figure 3. Morphology of ES cells and embryoid bodies (EBs).** A) Typical appearance of ES cell in culture. Cells form bright, shiny adherent colonies and are cultured on a layer of MEFs feeder cell layer. B) Typical appearance of EBs in culture. These three-dimensional structures are formed during ES cell differentiation and contain precursors for three primary germ cell layers. C) Embryoid bodies at day 7 of differentiation containing haemoglobin (indicated with aster‐

74 Pluripotent Stem Cell Biology - Advances in Mechanisms, Methods and Models

isks). Scale bar 300µm.

There are two alternative approaches to study the development of blast colonies; liquid culture on gelatine or semi-solid culture in methylcellulose. Liquid culture on gelatine facilitates harvesting of the cells for flow cytometry analysis or time-lapse imaging techniques. Alterna‐ tively, the semi-solid culture is a clonogenic assay allowing the growth of of individual blast colonies and their quantification. Individual blast colonies can then further be isolated and their haematopoietic potential assessed. To perform blast colony culture, EBs are harvested between day 2.5 and 3.5 and the cells are enriched for haemangioblast by sorting for the expression of FLK1, the receptor 2 for VEGF. Isolated FLK1 positive cells are then replated on gelatinised cell culture plates at the density of 8,5 x 104 cell/10cm2 in a medium containing IMDM supplemented with 10% FCS, 2mM L-Glutamine mM, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 15% endothelial cell line-D4T conditioned medium [11], 10 ng/ml interleukin 6 (IL-6) (Peprotech) and 5 ng/ml vascular endothelial growth factor (VEGF, Peprotech). For semi-solid cultures, FLK1 positive cells are seeded in 35 mm x 10 mm (BD Falcon) dishes at a density of 1.5 x 104 cells/ml in IMDM medium supplemented with 10% FCS, 2mM L-Glutamine mM, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 15% endothelial cell line-D4T conditioned medium [11], 10 ng/ml interleukin 6 (IL-6) (Peprotech) and 5 ng/ml vascular endothelial growth factor (VEGF, Peprotech). The medium is additionally supplemented with 10 g/L methylcellulose (dissolved in IMDM, Alfa-Aesar). Blast colonies are scored 3-4 days after replating. Plating 3.0 x 104 cells per 1 ml of semi-solid medium should result in the formation of around 300-400 of blast colonies.

It was recently shown that the haemangioblast progenitor can also generate vascular smooth muscle cells in addition to endothelial and haematopoietic cells [14, 77, 78]. Accordingly, Yamashita and co-workers also established that FLK1 positive cells can differentiate towards endothelial and mural cells [78]. Smooth muscle cells appear, under the microscope as large flat adherent cells. These cells express smooth muscle specific genes such as smooth muscle actin α, transgelin, calponin, and smooth muscle myosin heavy chain [79]. A typical immunostaining of α-smooth muscle actin is presented in figure 6A. We have recently analysed in more detail the generation of smooth muscle cells and showed that these cells are largely generated

Embryonic Stem Cell Differentiation – A Model System to Study Embryonic Haematopoiesis

http://dx.doi.org/10.5772/57607

77

**Figure 6. The differentiation of ES cells towards smooth muscle cells.** Representative image of smooth muscle cells stained with SMAα-Cy3 antibody. These cells also carry a transgenic BAC containing H2B-VENUS cDNA under the

During blast colony development *in vitro*, tight clusters of adherent cells that represent haemogenic endothelium are observed after around two days. This population of cells corresponds to an intermediate stage in the formation of blood cells from the haemangioblast. It was recently shown that these cells express TIE2 and VE-cadherin (markers of endothelial cells), c-KIT (expressed by both haematopoietic and endothelial cells) and are negative for CD41 (first haematopoietic marker) expression [20]. To isolate haemogenic endothelium, cells in day 2 liquid blast cultures are harvested and FACS sorted for the expression of these markers [20]. It is also possible to isolate a more advanced haemogenic endothelium cell population, positive for the expression of CD41 [70]. The cells are next plated onto gelatin-coated plates at

Glutamine, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 10 ng/ml Oncos‐ tatin M (R&D Systems) and 1 ng/ml bFGF (Peprotech). Cells are grown in standard culture

in IMDM medium supplemented with 10% FCS, 2 mM L-

control of the α-smooth muscle actin transcriptional regulatory elements. Scale bar 100µm.

**3.5. Culture of haemogenic endothelium**

a density of 1.2 x 105 cells/cm2

independently from the haemogenic endothelium [80].

Sorting EBs between day 2.5 and 3.5 for the expression of FLK1 enrich for blast colony forming cell (BL-CFCs) [11, 13]. During the formation of blast colonies, several distinct morphological stages can be observed. First, there is the formation of tight clusters of adherent cells and later, single round cells emerging from these tight clusters. These cells then proliferate (Figure 5A). Alongside these morphological changes, there is first an upregulation of endothelial markers, such as TIE2 and VE-Cadherin during the formation of the tight clusters. Later the expression of these endothelial markers is gradually downregulated whereas the expression of haema‐ topoietic markers (such as CD41 and then CD45) is upregulated (Figure 5B). This correlates with the emergence of round floating–haematopoietic-cells.

**Figure 5. The development of the blast colony during four days of differentiation**. A) Microscopic image showing representative blast colonies. Scale bar 100µm. B) Schematic representation of the different cell surface markers ex‐ pressed on cells during blast colony differentiation.

It was recently shown that the haemangioblast progenitor can also generate vascular smooth muscle cells in addition to endothelial and haematopoietic cells [14, 77, 78]. Accordingly, Yamashita and co-workers also established that FLK1 positive cells can differentiate towards endothelial and mural cells [78]. Smooth muscle cells appear, under the microscope as large flat adherent cells. These cells express smooth muscle specific genes such as smooth muscle actin α, transgelin, calponin, and smooth muscle myosin heavy chain [79]. A typical immunostaining of α-smooth muscle actin is presented in figure 6A. We have recently analysed in more detail the generation of smooth muscle cells and showed that these cells are largely generated independently from the haemogenic endothelium [80].

**Figure 6. The differentiation of ES cells towards smooth muscle cells.** Representative image of smooth muscle cells stained with SMAα-Cy3 antibody. These cells also carry a transgenic BAC containing H2B-VENUS cDNA under the control of the α-smooth muscle actin transcriptional regulatory elements. Scale bar 100µm.

#### **3.5. Culture of haemogenic endothelium**

ng/ml interleukin 6 (IL-6) (Peprotech) and 5 ng/ml vascular endothelial growth factor (VEGF, Peprotech). For semi-solid cultures, FLK1 positive cells are seeded in 35 mm x 10 mm (BD

2mM L-Glutamine mM, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 15% endothelial cell line-D4T conditioned medium [11], 10 ng/ml interleukin 6 (IL-6) (Peprotech) and 5 ng/ml vascular endothelial growth factor (VEGF, Peprotech). The medium is additionally supplemented with 10 g/L methylcellulose (dissolved in IMDM, Alfa-Aesar). Blast colonies

Sorting EBs between day 2.5 and 3.5 for the expression of FLK1 enrich for blast colony forming cell (BL-CFCs) [11, 13]. During the formation of blast colonies, several distinct morphological stages can be observed. First, there is the formation of tight clusters of adherent cells and later, single round cells emerging from these tight clusters. These cells then proliferate (Figure 5A). Alongside these morphological changes, there is first an upregulation of endothelial markers, such as TIE2 and VE-Cadherin during the formation of the tight clusters. Later the expression of these endothelial markers is gradually downregulated whereas the expression of haema‐ topoietic markers (such as CD41 and then CD45) is upregulated (Figure 5B). This correlates

**Figure 5. The development of the blast colony during four days of differentiation**. A) Microscopic image showing representative blast colonies. Scale bar 100µm. B) Schematic representation of the different cell surface markers ex‐

cells/ml in IMDM medium supplemented with 10% FCS,

cells per 1 ml of semi-solid medium should

Falcon) dishes at a density of 1.5 x 104

are scored 3-4 days after replating. Plating 3.0 x 104

result in the formation of around 300-400 of blast colonies.

76 Pluripotent Stem Cell Biology - Advances in Mechanisms, Methods and Models

with the emergence of round floating–haematopoietic-cells.

pressed on cells during blast colony differentiation.

During blast colony development *in vitro*, tight clusters of adherent cells that represent haemogenic endothelium are observed after around two days. This population of cells corresponds to an intermediate stage in the formation of blood cells from the haemangioblast. It was recently shown that these cells express TIE2 and VE-cadherin (markers of endothelial cells), c-KIT (expressed by both haematopoietic and endothelial cells) and are negative for CD41 (first haematopoietic marker) expression [20]. To isolate haemogenic endothelium, cells in day 2 liquid blast cultures are harvested and FACS sorted for the expression of these markers [20]. It is also possible to isolate a more advanced haemogenic endothelium cell population, positive for the expression of CD41 [70]. The cells are next plated onto gelatin-coated plates at a density of 1.2 x 105 cells/cm2 in IMDM medium supplemented with 10% FCS, 2 mM L-Glutamine, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG, 10 ng/ml Oncos‐ tatin M (R&D Systems) and 1 ng/ml bFGF (Peprotech). Cells are grown in standard culture conditions. Upon culture, at least 1-2% of these cells generate primitive and definitive haema‐ topoietic cells [20]. During this process, the haemogenic endothelial cell population (TIE2+,c-KIT+,CD41-) cells gradually acquire the expression of CD41. The cells then progress further and lose their expression of endothelial markers.

#### **3.6. Haematopoietic colonies assays**

To evaluate the presence of haematopoietic colonies, EBs should be harvested at day 4,5 or 6 and trypsinised. Cells from the EBs can be then directly used or alternatively sorted for the expression of a marker of haematopoietic progenitors, such as for example CD41. Approxi‐ mately 3.0 x 104 unsorted cells should be plated in 35 mm x 10 mm (BD Falcon) dishes in 1 ml of semisolid medium containing IMDM, 15% plasma derived serum (PDS) (Antech), 10% protein free hybridoma medium (PFHM, Gibco), 2mM L-Glutamine, 180 µg/ml transferrin, 25 µg/ml ascorbic acid, 4.6 x 10-4 M MTG and cytokines such as: 1% c-KIT ligand supernatant, 1% interleukin 3 supernatant (IL-3) (see [76]), 1µg/ml GM-CSF, 1% thrombopoietin condi‐ tioned media, 10 ng/ml IL-6 (Peprotech), 10 ng/ml macrophage colony stimulating factor (M-CSF), 5 ng/ml IL-11 (R&D Systems) and 4 U/ml of Erythropoietin (Ortho-Biotech) and 10 g/L methylcellulose (dissolved in IMDM, Alfa-Aesar). Haematopoietic colonies are assessed and scored based on their morphology. Primitive erythroid colonies are scored at day 5, whereas definitive haematopoietic colonies are usually enumerated 8 days after replating. Morphologic landmarks are used to distinguish the different types of haematopoietic colonies. Haemato‐ poietic progenitors can also be cultured in liquid conditions to allow easier access of the cells for subsequent flow cytometry analysis or cytospin assays. For that cells should be seeded at a density of 2.0 x 106 /ml in ultra low-adherence tissue culture plates (Costar) in the haemato‐ poietic medium described above with methylcellulose being replaced with IDMD medium.

The onset of emergence of primitive erythroid cells is observed within EBs by day 4 of differentiation [81]. Definitive erythroid and macrophage precursors appear shortly after and are followed by mast cells and multilineage precursors [81]. Primitive colonies appear around day 4 of culture. These colonies are round, compact and bright red in colour. By day 6-7 of culture, morphologically distinguishable definitive haematopoietic colonies are detected (Figure 7).

Although only very succinctly discussed in this chapter, this system is also very amenable to examination of the cell signalling pathways that support the development of normal haema‐ topoiesis [33, 34, 82-84]. Interestingly these pathways are also implicated in leukaemogenesis [85, 86]. Finally with the advent of novel human ES cells or human iPS cells [87], that recapit‐ ulate better the ground state and are easier to work with, and development of new methods facilitating genome editing [88-92], this experimental system is very likely to be instrumental for delivering new advances in our understanding of human haematopoietic development,

**Figure 7.** Examples of haematopoietic colonies obtained during ES cell differentiation A) Primitive erythroid colony – red in colour, compact and relatively small. B) Definitive myeloid colony – bigger in size, white and looser cells. C)

Embryonic Stem Cell Differentiation – A Model System to Study Embryonic Haematopoiesis

http://dx.doi.org/10.5772/57607

79

We apologize to the many colleagues whose work could not be cited owing to space con‐ straints. Our work is supported by Cancer Research UK (CRUK), Leukaemia and Lymphoma research (LLR) and the Biotechnology and Biological Sciences Research Council (BBSRC).

that is otherwise very difficult to study *in vivo*.

Mixed haematopoietic colony. Scale bar 50µm.

**Acknowledgements**


Author details

Klara Gadó1 \* and Gyula Domján2

\*Address all correspondence to<EMAIL_ADDRESS>
1 Faculty of Health Sciences, Department of Clinical Studies, Semmelweis University of Medicine, Budapest, Hungary

2 Faculty of Medicine, 1st Department of Medicine, Semmelweis University of Medicine, Budapest, Hungary

## References

[9] Joshi A, Hong J, Siva C. Recurrent thrombosis in patients with antiphospholipid syndrome receiving newer oral anticoagulants: A case report and review of literature. Clinical Medicine & Research. 2017;15(1–2):41-44

Riccieri V, Piccoli A, Valesini G, Doria A, Meroni PL, Tincani A. Risk factors for a first thrombotic event in antiphospholipid antibody carriers. A multicentre, retrospective

Antiphospholipid Syndrome and Thrombocytopenia http://dx.doi.org/10.5772/intechopen.72509 65

[23] Manukyan D, Rossmann H, Schulz A, Zeller T, Pfeiffer N, Binder H, Münzel T, Beutel ME, Müller-Calleja N, Wild PS, Lackner KJ. Distribution of antiphospholipid antibodies in a large population-based German cohort. Clinical Chemistry and Laboratory Medicine.

[24] Neville C, Rauch J, Kassis J, Chang ER, Joseph L, Le Comte M, Fortin PR. Thromboembolic risk in patients with high titre anticardiolipin and multiple antiphospholipid anti-

[25] Sciascia S, Sanna G, Khamashta MA, Cuadrado MJ, Erkan D, Andreoli L, et al. The estimated frequency of antiphospholipid antibodies in young adults with cerebrovascular events: A systematic review. Annals of the Rheumatic Diseases. 2015;74(11):2028-2033

[26] Sciascia S, Amigo MC, Roccatello D, Khamashta M. Diagnosing antiphospholipid syndrome: 'extra-criteria' manifestations and technical advances. Nature Reviews Rheuma-

[27] Cervera R, Piette JC, Font J, et al. Antiphospholipid syndrome: Clinical and immunologic manifestations and patterns of disease expression in a cohort of 1,000 patients. Arthritis

[28] Baroni G, Banzato A, Bison E, Denas G, Zoppellaro G, Pengo V. The role of platelets in

[29] Proulle V, Furie RA, Merrill-Skoloff G, Furie BC, Furie B. Platelets are required for enhanced activation of the endothelium and fibrinogen in a mouse thrombosis model of

[30] Vlachoyiannopoulos PG, Routsias JG. A novel mechanism of thrombosis in antiphospholipid antibody syndrome. Journal of Autoimmunity. 2010;35(3):248-255 [31] Breen KA, Sanchez K, Kirkman N, Seed PT, Parmar K, Moore GW, Hunt BJ. Endothelial and platelet microparticles in patients with antiphospholipid antibodies. Thrombosis

[32] Ali N, Auerbach HE. New-onset acute thrombocytopenia in hospitalized patients: Pathophysiology and diagnostic approach. Journal of Community Hospital Internal Medicine

[33] Buckley MF, James JW, Brown DE, et al. A novel approach to the assessment of variations in the human platelet count. Thrombosis and Haemostasis. 2000;83(3):480-484

[34] Artim-Esen B, Diz-Küçükkaya R, İnanç M. The significance and management of thrombocytopenia in antiphospholipid syndrome. Current Rheumatology Reports. 2015;17(3):14

[35] Cuadrado MJ, Mujic F, Muñoz E, Khamashta MA, Hughes GRV. Thrombocytopenia in the antiphospholipid syndrome. Annals of the Rheumatic Diseases. 1997;56:194-196

follow-up study. Annals of the Rheumatic Diseases. 2009;68(3):397-399

bodies. Thrombosis and Haemostasis. 2003;90(1):108-115

2016;54(10):1663-1670

tology. 2017;13(9):548-560

and Rheumatism. 2002;46(4):1019-1027

APS. Blood. 2014;124(4):611-622

Research. 2015;135(2):368-374

Perspectives. 2017;7(3):157-167

antiphospholipid syndrome. Platelets. 2017;7:1-5

Riccieri V, Piccoli A, Valesini G, Doria A, Meroni PL, Tincani A. Risk factors for a first thrombotic event in antiphospholipid antibody carriers. A multicentre, retrospective follow-up study. Annals of the Rheumatic Diseases. 2009;68(3):397-399

[23] Manukyan D, Rossmann H, Schulz A, Zeller T, Pfeiffer N, Binder H, Münzel T, Beutel ME, Müller-Calleja N, Wild PS, Lackner KJ. Distribution of antiphospholipid antibodies in a large population-based German cohort. Clinical Chemistry and Laboratory Medicine. 2016;54(10):1663-1670

[9] Joshi A, Hong J, Siva C. Recurrent thrombosis in patients with antiphospholipid syndrome receiving newer oral anticoagulants: A case report and review of literature. Clin-

[10] Arnaud L, Mathian A, Ruffatti A, Erkan D, Tektonidou M, Cervera R, Forastiero R, Pengo V, Lambert M, Martinez-Zamora MA, Balasch J, Zuily S, Wahl D, Amoura Z. Efficacy of aspirin for the primary prevention of thrombosis in patients with antiphospholipid antibodies: An international and collaborative meta-analysis. Autoimmunity Reviews. 2014;13(3):281-291

[11] Derksen RH, de Groot PG. Towards evidence-based treatment of thrombotic antiphospholipid

[12] Hoppe B, Burmester GR, Dörner T. Heparin or aspirin or both in the treatment of recurrent abortions in women with antiphospholipid antibody (syndrome). Current

[13] Okuma H, Kitagawa Y, Yasuda T, Tokuoka K, Takagi S. Comparison between single antiplatelet therapy and combination of antiplatelet and anticoagulation therapy for secondary prevention in ischemic stroke patients with antiphospholipid syndrome. Inter-

[14] Al Marzooqi A, Leone A, Al Saleh J, Khamashta M. Current status and future prospects for the treatment of antiphospholipid syndrome. Expert Review of Clinical Immunology.

[15] Espinosa G, Rodríguez-Pintó I, Cervera R. Catastrophic antiphospholipid syndrome: An

[16] Merashli M, Noureldine MH, Uthman I, Khamashta M. Antiphospholipid syndrome: An update. European Journal of Clinical Investigation. 2015;45(6):653-662. DOI: 10.1111/eci.12449

[17] Chaturvedi S, McCrae KR. Clinical risk assessment in the antiphospholipid syndrome: Current landscape and emerging biomarkers. Current Rheumatology Reports. 2017;19(7):43

[18] Oku K, Amengual O, Yasuda S, Atsumi T. How to identify high-risk APS patients: Clinical utility and predictive values of validated scores. Current Rheumatology Reports. 2017;19(8):51

[19] Miyakis S, Lockshin MD, Atsumi T, et al. International consensus statement on an update of the classification criteria for definite antiphospholipid syndrome (APS). Journal of

[20] Mulla MJ, Brosens JJ, Chamley LW, Giles I, Pericleous C, Rahman A, et al. Antiphospholipid antibodies induce a pro-inflammatory response in first trimester trophoblast via the TLR4/

[21] Mustonen P, Lehtonen KV, Javela K, Puurunen M. Persistent antiphospholipid antibody (aPL) in asymptomatic carriers as a risk factor for future thrombotic events: A nationwide

[22] Ruffatti A, Del Ross T, Ciprian M, Nuzzo M, Rampudda M, Bertero MT, Bergia R, Caramaschi P, Biasi D, Capsoni F, Montaguti L, Ruffini R, Brucato A, Picillo U, Fanelli V,

MyD88pathway. American Journal of Reproductive Immunology. 2009;62:96-111

ical Medicine & Research. 2017;15(1–2):41-44

Opinion in Rheumatology. 2011;23(3):299-304

national Journal of Medical Sciences. 2009;7(1):15-18

update. Panminerva Medica. 2017;59(3):254-268

Thrombosis and Haemostasis. 2006;4:295-306

prospective study. Lupus. 2014;23(14):1468-1476

syndrome. Lupus. 2010;19(4):470-474

2016;12(9):927-935

64 Thrombocytopenia

[36] Stasi R, Stipa E, Masi M, Oliva F, Sciarra A, Perroti A, et al. Prevalence and clinical significance of elevated antiphospholipid antibodies in patients with idiopathic thrombocytopenic purpura. Blood. 1994;84:4203-4208

[49] Martirosyan A, Petrek M, Kishore A, Manukyan G. Immunomodulatory effects of therapeutic plasma exchange on monocytes in antiphospholipid syndrome. Experimental and

Antiphospholipid Syndrome and Thrombocytopenia http://dx.doi.org/10.5772/intechopen.72509 67

[50] Gustavsen A, Skattum L, Bergseth G, Lorentzen B, Floisand Y, Bosnes V, Mollnes TE, Barratt-Due A. Effect on mother and child of eculizumab given before caesarean section in a patient with severe antiphospholipid syndrome: A case report. Medicine (Baltimore).

Therapeutic Medicine. 2016;12(2):1189-1195

2017;96(11):e6338

[49] Martirosyan A, Petrek M, Kishore A, Manukyan G. Immunomodulatory effects of therapeutic plasma exchange on monocytes in antiphospholipid syndrome. Experimental and Therapeutic Medicine. 2016;12(2):1189-1195

[36] Stasi R, Stipa E, Masi M, Oliva F, Sciarra A, Perroti A, et al. Prevalence and clinical significance of elevated antiphospholipid antibodies in patients with idiopathic thrombo-

[37] Harris EN, Gharavi AE, Hedge U, Derue G, Morgan SH, Englert H, et al. Anticardiolipin antibodies in autoimmune thrombocytopenic purpura. British Journal of Haematology.

[38] Abreu MM, Danowski A, Wahl DG, Amigo MC, Tektonidou M, Pacheco MS, Fleming N, Domingues V, Sciascia S, Lyra JO, Petri M, Khamashta M, Levy RA. The relevance of "noncriteria" clinical manifestations of antiphospholipid syndrome: 14th International Congress on Antiphospholipid Antibodies Technical Task Force Report on Antiphospholipid Syn-

[39] Demetrio Pablo R, Muñoz P, López-Hoyos M, Calvo V, Riancho L, Martínez-Taboada VM. Thrombocytopenia as a thrombotic risk factor in patients with antiphospholipid antibodies

[40] Palomo I, Toro C, Alarcon M. The role of platelets in the pathophysiology of atheroscle-

[41] Kamath S, Blann AD, Lip GY. Platelet activation: Assessment and quantification. Euro-

[42] Korkmaz S, Uslu AU, Sahin S, Senel S, Sencan MI. There a link between mean platelet volume and thrombotic events in antiphospholipid syndrome? Platelets. 2014;25(5):343-347

[43] Rupa-Matysek J, Gil L, Wojtasinska E, Ciepluch K, Lewandowska M, Komarnicki M. The relationship between mean platelet volume and thrombosis recurrence in patients diagnosed with antiphospholipid syndrome. Rheumatology International. 2014;34(11):1599-1605

[44] Myers B. Diagnosis and management of maternal thrombocytopenia in pregnancy. Brit-

[45] Gernsheimer T, James AH, Stasi R. How I treat thrombocytopenia in pregnancy. Blood.

[46] Sun D, Shehata N, Ye XY, Gregorovich S, De France B, Arnold DM, Shah PS, Malinowski AK. Corticosteroids compared with intravenous immunoglobulin for the treatment of immune

[47] Gamoudi D, Cutajar M, Gamoudi N, Camilleri DJ, Gatt A. Achieving a satisfactory clinical and biochemical response in antiphospholipid syndrome and severe thrombocytopenia with rituximab: Two case reports. Clinical Case Reports. 2017;5(6):845-848 [48] Pons-Estel GJ, Andreoli L, Scanzi F, Cervera R, Tincani A. The antiphospholipid syndrome in patients with systemic lupus erythematosus. Journal of Autoimmunity.

drome Clinical Features. Autoimmunity Reviews. 2015;14(5):401-414

rosis (review). Molecular Medicine Reports. 2008;1(2):179-184

thrombocytopenia in pregnancy. Blood. 2016;128(10):1329-1335

pean Heart Journal. 2001;22(17):1561-1571

ish Journal of Haematology. 2012;158(1):3-15

2013;121(1):38-47

2017;76:10-20

without disease criteria. Medicina Clínica (Barcelona). 2017;148(9):394-400

cytopenic purpura. Blood. 1994;84:4203-4208

1985;59:231-234

66 Thrombocytopenia

[50] Gustavsen A, Skattum L, Bergseth G, Lorentzen B, Floisand Y, Bosnes V, Mollnes TE, Barratt-Due A. Effect on mother and child of eculizumab given before caesarean section in a patient with severe antiphospholipid syndrome: A case report. Medicine (Baltimore). 2017;96(11):e6338

**Chapter 5**

**Provisional chapter**

**Interferon-Induced Thrombotic Microangiopathy**

**Interferon-Induced Thrombotic Microangiopathy**

DOI: 10.5772/intechopen.74226

Interferon (IFN) is an effective therapy for multiple disorders. An infrequently reported side effect is thrombotic microangiopathy (TMA): thrombotic thrombocytopenic purpura (TTP) and hemolytic uremic syndrome (HUS). We published the first comprehensive review analyzing this association with the following observations: (1) there was a higher incidence of IFN-induced TMA in myeloproliferative disorders (chronic myelogenous leukemia (CML)) than that in nonmalignant disorders (multiple sclerosis (MS), chronic hepatitis C virus infection (HCV)); (2) mean age at diagnosis was 47 years; (3) there was rare association with hairy cell leukemia (HCL), Sezary syndrome (one case each) and no cases reported for polycythemia vera (PV); (4) sex distribution was balanced (exception of higher prevalence in females for MS); (5) TMA was insidious in onset with long incubation periods (average treatment duration 40.4 months); (6) comparative analysis of mean time (months) to onset of TMA ensuing cumulative IFN exposure was: MS 68.6 vs. CML 35.5 vs. HCV 30.4; (7) confirmed TTP (low ADAMTS 13 levels) was associated with the presence of an inhibitor; (8) outcome analysis revealed complete remission in 27 (40%), persistent chronic kidney disease in 28 (42%) and fatality in 12 patients (18%); (9)

corticosteroids, plasma exchange (PEX) and rituximab are effective therapies.

**Keywords:** thrombocytopenia, interferon, thrombotic microangiopathy (TMA), thrombotic thrombocytopenia purpura (TTP), myeloproliferative disorder, HCV

Thrombocytopenia is a common side effect of IFN treatment. TMA is a distinct clinical entity with potentially fatal consequences without precise and expeditious intervention. Classic clinical presentation includes a triad of anemia, thrombocytopenia, and evidence of microangiopathic hemolysis on the peripheral blood smear. The quintessential pathologic feature is the development of microvascular thrombi affecting small or larger vessels with variable

> © 2016 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2018 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Ajay Kundra, Mirela Andrei and Jen Chin Wang

Ajay Kundra, Mirela Andrei and Jen Chin Wang

Additional information is available at the end of the chapter

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/intechopen.74226

**Abstract**

**1. Introduction**

#### **Interferon-Induced Thrombotic Microangiopathy Interferon-Induced Thrombotic Microangiopathy**

DOI: 10.5772/intechopen.74226

Ajay Kundra, Mirela Andrei and Jen Chin Wang Ajay Kundra, Mirela Andrei and Jen Chin Wang

Additional information is available at the end of the chapter Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/intechopen.74226

#### **Abstract**

Interferon (IFN) is an effective therapy for multiple disorders. An infrequently reported side effect is thrombotic microangiopathy (TMA): thrombotic thrombocytopenic purpura (TTP) and hemolytic uremic syndrome (HUS). We published the first comprehensive review analyzing this association with the following observations: (1) there was a higher incidence of IFN-induced TMA in myeloproliferative disorders (chronic myelogenous leukemia (CML)) than that in nonmalignant disorders (multiple sclerosis (MS), chronic hepatitis C virus infection (HCV)); (2) mean age at diagnosis was 47 years; (3) there was rare association with hairy cell leukemia (HCL), Sezary syndrome (one case each) and no cases reported for polycythemia vera (PV); (4) sex distribution was balanced (exception of higher prevalence in females for MS); (5) TMA was insidious in onset with long incubation periods (average treatment duration 40.4 months); (6) comparative analysis of mean time (months) to onset of TMA ensuing cumulative IFN exposure was: MS 68.6 vs. CML 35.5 vs. HCV 30.4; (7) confirmed TTP (low ADAMTS 13 levels) was associated with the presence of an inhibitor; (8) outcome analysis revealed complete remission in 27 (40%), persistent chronic kidney disease in 28 (42%) and fatality in 12 patients (18%); (9) corticosteroids, plasma exchange (PEX) and rituximab are effective therapies.

**Keywords:** thrombocytopenia, interferon, thrombotic microangiopathy (TMA), thrombotic thrombocytopenia purpura (TTP), myeloproliferative disorder, HCV


*Article* **Commensal and Opportunistic Bacteria Present in the Microbiota in Atlantic Cod (***Gadus morhua***) Larvae Differentially Alter the Hosts' Innate Immune Responses**

**Ragnhild Inderberg Vestrum 1,† , Torunn Forberg 1,† , Birgit Luef 1 , Ingrid Bakke 1 , Per Winge 2 , Yngvar Olsen 2 and Olav Vadstein 1,\***

**Abstract:** The roles of host-associated bacteria have gained attention lately, and we now recognise that the microbiota is essential in processes such as digestion, development of the immune system and gut function. In this study, Atlantic cod larvae were reared under germ-free, gnotobiotic and conventional conditions. Water and fish microbiota were characterised by 16S rRNA gene analyses. The cod larvae's transcriptional responses to the different microbial conditions were analysed by a custom Agilent 44 k oligo microarray. Gut development was assessed by transmission electron microscopy (TEM). Water and fish microbiota differed significantly in the conventional treatment and were dominated by different fast-growing bacteria. Our study indicates that components of the innate immune system of cod larvae are downregulated by the presence of non-pathogenic bacteria, and thus may be turned on by default in the early larval stages. We see indications of decreased nutrient uptake in the absence of bacteria. The bacteria also influence the gut morphology, reflected in shorter microvilli with higher density in the conventional larvae than in the germ-free larvae. The fact that the microbiota alters innate immune responses and gut morphology demonstrates its important role in marine larval development.

**Keywords:** Atlantic cod; microbiota; innate immune system; germ-free; gnotobiotic

### **1. Introduction**

The roles of the microbiota associated with vertebrate hosts, including fish, have received much attention over the last decade. Several studies have shown that the microbiota stimulates the immune system and functions as a barrier against potential pathogens [1–4], aids in epithelial development and maturation [4,5] and affects the digestion of nutrients [6,7]. There is a bias in the type of animals studied, and still relatively few studies are published on the function of microbiota in fish. Most of the bacteria associated with the fish are harmless or beneficial [8,9]. However, specific pathogens and opportunistic bacteria are also present [10,11], and bacteria present in the natural environment of the fish cause many of the infections that are associated with the mortality of marine fish larvae [12].

Germ-free animals have been popular tools used in studies of host–microbe interactions [13,14], and the use of gnotobiotic zebrafish is well-known [15,16]. Rawls et al. [16] observed that gut microbiota in zebrafish stimulated proliferation of intestinal epithelial cells, as previously seen for rodents [5,17]. They found 212 genes that were differentially regulated in germ-free fish compared to fish exposed to bacteria. Moreover, 59 of those gene-expression responses were observed in both mice and zebrafish. These genes are involved in epithelial proliferation, nutrient metabolism and innate immune responses [16].

**Citation:** Vestrum, R.I.; Forberg, T.; Luef, B.; Bakke, I.; Winge, P.; Olsen, Y.; Vadstein, O. Commensal and Opportunistic Bacteria Present in the Microbiota in Atlantic Cod (*Gadus morhua*) Larvae Differentially Alter the Hosts' Innate Immune Responses. *Microorganisms* **2022**, *10*, 24. https://doi.org/10.3390/ microorganisms10010024

Academic Editor: Francesco Marotta

Received: 15 November 2021 Accepted: 21 December 2021 Published: 24 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

Our group have developed a protocol and a cultivation system for germ-free and gnotobiotic Atlantic cod larvae [18]. Using this system, Forberg et al. [19,20] showed that bacteria regulate the transcription of genes involved in immune response and nutrient digestion in the cod larvae, and that the response differs depending on whether the bacteria are dead or alive.

However, at that time, molecular tools for Atlantic cod were poorly developed, and this allowed characterisation of only a limited part of the host transcriptome. Following the sequencing of the Atlantic cod genome [21], it was discovered that cod lack certain families of cell surface receptors, whereas others are expanded, compared with zebrafish and mice [21–24]. The functional consequences of this evolutionary divergence in immune response with regards to host–microbe interactions have, to our knowledge, not been studied. Thus, the knowledge about the early regulation of the immune system of cod and its responses to both pathogenic and non-pathogenic bacteria is still limited.

The aim of this study was to increase our knowledge about host responses induced by the early microbial colonisation of Atlantic cod larvae, focusing on the roles of microbes in regulation of the immune system and digestion. Newly hatched cod larvae were reared under three different conditions: germ-free, gnotobiotic (two probiotic candidate strains added) and conventional (uncontrolled microbial environments). The water- and fishassociated microbiota was characterised and analysed, and host responses were examined using a custom oligo microarray as well as transmission electron microscopy (TEM).

#### **2. Materials and Methods**

The experiment was carried out within the Norwegian animal welfare act guidelines, in accordance with the Animal Welfare Act of 20 December 1974, amended 19 June 2009, at a facility with permission to conduct experiments on fish (code 93) provided by the Norwegian Animal Research Authority (NARA). The experiment was approved by NARA.

#### *2.1. Cod Larval Rearing*

Atlantic cod eggs were delivered from Nofima Marin national breeding station (Havbruksstasjonen Tromsø, Norway). Upon arrival, the cod eggs (55–65-day degrees) were acclimatised in filtered (0.22 µm Micropore®) autoclaved (121 ◦C, 20 min) seawater (FASW) at 6 ± 1 ◦C, in the dark. Germ-free larvae were obtained according to the protocol of Forberg et al. [18] (information about germ-free verification in Text S1). Cod larvae (65 larvae in 2 L water) were reared under three different conditions: germ-free, conventional and gnotobiotic. For the gnotobiotic treatment, two different bacterial strains were added in equal amounts (final density of 106 cells/mL) to the rearing bottles: *Microbacterium* ND 2–7 and *Vibrio* RD 5–30, both previously isolated from cod and identified as probiotic candidates (for details, see [25]) (information about live feed and bacterial cultures in Text S2). The rearing bottles representing the conventional condition were filled with microbially matured water, from a biofilter in a seawater lab-scale aquaculture system. After stocking, the temperature was increased by 1 ◦C/day until 12 ◦C was reached.

#### *2.2. Sampling*

Each treatment had 11 replicate bottles at trial start. Cod larvae were collected at 4 (only for DNA extraction), 8, 13 and 16 dph. Fish from one bottle were sampled at 4 dph, while 3 replicate bottles were sampled at 8, 13 and 16 dph. Water was sampled from one bottle at 1 and 4 dph, and three replicate bottles at 8, 13 and 16 dph. After sampling, the bottles were taken out of the experiment, thus reducing the number of replicate bottles with time. Larvae were sacrificed with an overdose of tricaine methanesulfonate (MS-222) prior to sampling, snap-frozen in liquid nitrogen and stored at −80 ◦C for further analyses. To investigate larval growth, 10 individual larvae from each sampled bottle were freeze-dried and weighed. More details regarding sampling procedures and DNA/RNA extraction are described in Supplementary Materials (Text S3).

#### *2.3. Characterisation of Microbial Communities*

The water and fish samples from the gnotobiotic treatment were analysed by DGGE. PCR products representing the V3 region of the 16 S rRNA gene were generated using a nested PCR protocol to avoid possible amplification of eukaryotic 18S rDNA [26]. The PCR was set up and analysed as described by Bakke et al. [27].

For in-depth analysis of the microbiota in the conventional treatment, Illumina MiSeq sequencing was performed based on total DNA extracted from water and larvae sampled at 1 (only water), 4, 8, 13 and 16 dph. Larval and water samples were prepared for Illumina MiSeq sequencing by amplification of the V4 region of the 16S rRNA gene, by using the following primers (bacteria-specific V4 primer, underlined and bold) including 50 overhang, as suggested by Illumina:

#### 515F F0 TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGNNNN**GTGCCAGCM** -**GCCGCGGTAA** 3 0 and

#### 803 R 50 GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGNNNN**CTACVVGGG** -**TATCTAAKCCBK** 3 0 .

The amplicon library preparation and processing of the Illumina sequencing data was performed as described by Vestrum et al. [28]. In short, for the first stage of amplification, the reactions were run for 38 cycles for water samples and 40 cycles for cod larval samples (98 ◦C 15 s, 55 ◦C 20 s, 72 ◦C 20 s), with 0.3 µM of each primer, 0.25 mM of each dNTP, 2 mM of MgCl2, 12 µM of BSA, glycerol (10%), Phusion Hot Start II High-Fidelity DNA Polymerase and reaction buffer from Thermo Scientific in a total volume of 20 µL. All samples were normalised using the SequalPrepTM Normalisation Plate Kit (Invitrogen). A second PCR was performed to attach dual indices and Illumina sequencing adapters to the normalised v4 amplicons by using the Nextera XT Index Kit. The indexed PCR products were normalised as described above, pooled, and concentrated by using Amicon® Ultra-0.5 Centrifugal Filter Devices. The resulting amplicon library was sequenced on a MiSeq lane (Illumina, San Diego, CA, USA) with v4 reagents employing 260 bp pairedend reads at the Norwegian Sequencing Center at the University of Oslo, Norway. The Illumina sequencing data were processed with the high-performance USEARCH utility (version 11) (http://drive5.com/usearch/features.html (accessed on 22 March 2019). Taxonomy assignment was performed applying the Sintax script [29] with a confidence value threshold of 0.8 and the RDP reference dataset (version 16). OTUs of particular interest were further analysed with the RDP tools [30] Classifier and Sequence Match. OTUs representing algae, Archaea and Cyanobacteria/Chloroplast were removed from the OTU table. In addition, an OTU representing Propionibacterium acne, a well-known contaminant of DNA extraction kits [31], was removed. To remove biases due to variation in sequencing depth, analyses were performed on an OTU table that had been subsampled to 15500 sequencing reads for each sample. The subsampling threshold was chosen based on the sample with the lowest number of reads in order to keep all samples in the dataset and was performed to avoid bias due to differences in sequencing depth. The resulting Illumina sequencing data were deposited at the European Nucleotide Archive (accession numbers ERS8484975-ERS8484994).

#### *2.4. Microarray Design, Hybridisation and Annotation*

A custom, Agilent 44 k oligo microarray (A-MEXP-2226, ArrayExpress, EMBL-EBI) described by Kleppe et al. [32] was used and analysed as described by Vestrum et al. [33]. This microarray design is partly based on the Atlantic cod gene set described by Star et al. [21] as well as EST sequences from various cod tissues/developmental stages. The identified differentially regulated transcripts were used for biological term enrichment analysis and Gene Ontology term (GO term) annotation in DAVID (Database of Annotation, Visualisation and Integrated Discovery) [34,35] (using the official gene symbol for human homologues).

#### *2.5. Electron Microscopy Procedures*

For processing fish larvae for transmission electron microscopy (TEM), the protocol from Galloway et al. [36] was adopted. Shortly, germ-free and conventional reared fish larvae from 16 dph were fixed in a mixture of 2.5% paraformaldehyde, 2.5% glutardialdehyde, 0.5% sucrose and 0.11 M HEPES buffer (pH 7.4) and stored at 4 ◦C until further processing. Three larvae from each treatment were rinsed in 0.11 M HEPES buffer, post-fixed in 2% OsO4 in 1.5% potassium ferricyanide (final concentration), bulk contrasted in 1.5% uranyl acetate and dehydrated in ethanol before embedding in Epon. After polymerisation, 50–60 nm sections of the fish midgut were cut using a Leica UC6 Ultramicrotome. These ultrathin sections were collected on 200-mesh copper grids and contrasted with 4% uranyl acetate and 1% lead citrate. Sections were inspected with a FEI Company Tecnai 12 operated at 80 kV and imaged using a digital MORADA G3 CCD camera (EMSIS).

#### *2.6. Intestinal Morphometry and Statistical Methods*

Computerised morphometric measurements of microvilli lengths (from tip to base, l), diameter (2r) and abundance of microvilli (µm–2) in the midgut of the fish were made using the image processing program iTEM (Olympus Soft Imaging Solutions GmbH). Microvilli parameters were measured according to the criteria of Brown [37]. Three fish larvae per treatment were investigated (*n* = 3). For length measurements, at least 65 microvilli per fish were analysed; in total, 318 microvilli within the germ-free and 294 microvilli within the conventional treatment were measured. For diameter measurements, at least 58 microvilli per fish were analysed; in total, 214 microvilli within the germ-free and 287 microvilli within the conventional treatment were measured. To quantify the abundance of microvilli, microvilli within a total area of around 400 µm2 were counted. iTEM was used to adjust contrast in the images and to insert calibrated scale bars into images.

#### *2.7. Statistical Analyses*

Student's *t*-test (unpaired) was used to investigate significance in differences in Shannon indices, abundance of individual DGGE bands and larval growth measurements. Survival analysis was performed by the Kaplan–Meier method, and the Log-rank test was used for pairwise post hoc comparisons of survival across the groups. Ordination by Principal Coordinate Analysis (PCoA) based on Bray–Curtis similarities was used to visualize differences between sample groups, and one-way and two-way PERMANOVA based on Bray–Curtis similarities were used to test for statistically significant differences between sample groups. Similarity Percentage analysis (SIMPER) was used to identify OTUs responsible for differences (measured as Bray–Curtis similarities) between different sample groups. The multivariate analyses were performed using the program package PAST version 3.22 [38]. Venn diagrams were created using jvenn [39]. The Usearch commands Alpha\_div and Sintax\_summary were used to calculate alpha diversity indices and to generate taxa summary tables (at various taxonomic levels, as specified with the results), respectively.

Data from the intestinal morphometric study were statistically analysed with IBM SPSS Statistics (SPSS for Windows, version 26.0; SPSS Inc., Chicago, IL, USA). A Welch test was performed to investigate significant differences between the axenic- and conventionaltreated fish regarding microvilli length, diameter and abundance of microvilli (µm–2). Differences were considered statistically significant when *p* ≤ 0.05.

#### **3. Results**

#### *3.1. Larval Survival and Growth*

Daily counts of dead larvae in the rearing bottles were used to calculate the percent of survival. Kaplan–Meier survival curves showed a clear separation of the conventional group vs. the germ-free and gnotobiotic group cumulative mortality, and this difference was highly significant (Log rank post hoc *p*-values 0.000032 and 0.000005 for the pairwise comparisons) (Supplementary Figure S1) (84.9%, 84.8% and 76.0% survival, respectively). **3. Results**

*3.1. Larval Survival and Growth*

At 16 dph, the dry weight of germ-free larvae was significantly lower (average 75.5 µg) than the gnotobiotic and conventional larvae (average 96.9 and 110.1 µg, respectively) (*p* = 0.017 and 0.059, respectively) (Supplementary Figure S2). *3.2. Composition of Fish and Water Microbiota*  The DGGE profiles for the samples from the gnotobiotic rearing bottles

conventional-treated fish regarding microvilli length, diameter and abundance of microvilli (µm–2). Differences were considered statistically significant when *p* ≤ 0.05.

Daily counts of dead larvae in the rearing bottles were used to calculate the percent of survival. Kaplan–Meier survival curves showed a clear separation of the conventional group vs. the germ-free and gnotobiotic group cumulative mortality, and this difference was highly significant (Log rank post hoc *p*-values 0.000032 and 0.000005 for the pairwise comparisons) (Supplementary Figure S1) (84.9%, 84.8% and 76.0% survival, respectively). At 16 dph, the dry weight of germ-free larvae was significantly lower (average 75.5 µg) than the gnotobiotic and conventional larvae (average 96.9 and 110.1 µg, respectively) (*p*

#### *3.2. Composition of Fish and Water Microbiota* (Supplementary Figure S3), where only two bacterial strains were added, were

*Microorganisms* **2022**, *9*, x FOR PEER REVIEW 5 of 18

= 0.017 and 0.059, respectively) (Supplementary Figure S2).

The DGGE profiles for the samples from the gnotobiotic rearing bottles (Supplementary Figure S3), where only two bacterial strains were added, were consistently identical, except for the presence of one additional band in one fish sample at 8 dph. This suggests the presence of a contaminating bacterial strain. Even though bacteria were added to the same final cell density in the rearing water in the gnotobiotic treatment, *V. gallicus* clearly dominated both in water and fish samples. *Microbacterium* was detectable at low levels in water samples and present only in some fish samples. consistently identical, except for the presence of one additional band in one fish sample at 8 dph. This suggests the presence of a contaminating bacterial strain. Even though bacteria were added to the same final cell density in the rearing water in the gnotobiotic treatment, *V. gallicus* clearly dominated both in water and fish samples. *Microbacterium* was detectable at low levels in water samples and present only in some fish samples. Fish and water microbiota in the conventional treatment were characterised by amplicon sequencing. After quality trimming and chimera removal, 1,039,322 reads were

Fish and water microbiota in the conventional treatment were characterised by amplicon sequencing. After quality trimming and chimera removal, 1,039,322 reads were obtained. Two water samples and two fish samples from 4 dph were removed due to low number of reads. The estimated total (Chao1) and observed number of OTUs for each sample (Figure 1) indicate a sequencing depth of on average 95% and 82% in fish and water samples, respectively. The observed richness was generally higher in water than in fish. obtained. Two water samples and two fish samples from 4 dph were removed due to low number of reads. The estimated total (Chao1) and observed number of OTUs for each sample (Figure 1) indicate a sequencing depth of on average 95% and 82% in fish and water samples, respectively. The observed richness was generally higher in water than in fish.

**Figure 1.** Observed richness (number of OTUs) and Chao1 values in fish and water microbiota at 1 (only water) to 16 dph (D1–D16). Numbers indicate the rearing bottle sampled. One fish sample consisted of 5 pooled larvae from one rearing bottle. **Figure 1.** Observed richness (number of OTUs) and Chao1 values in fish and water microbiota at 1 (only water) to 16 dph (D1–D16). Numbers indicate the rearing bottle sampled. One fish sample consisted of 5 pooled larvae from one rearing bottle.

A PCoA ordination of the microbial communities in fish and water samples (Figure 2a) showed that the water and fish clustered separately. This was corroborated by Bray– A PCoA ordination of the microbial communities in fish and water samples (Figure 2a) showed that the water and fish clustered separately. This was corroborated by Bray– Curtis similarities (Figure 2b). The PCoA plot also indicated that both the water and fish microbiota changed over time. There were significant differences between the water and fish microbiota both early (1/4–8 dph) and late (13–16 dph) in the experiment (one-way PERMANOVA, *p* = 0.04 and 0.01 for early and late, respectively) and also between water samples early and late in the experiment (one-way PERMANOVA, *p* = 0.01). There were no significant differences in the microbiota of fish samples early and late in the experiment.

experiment.

**Figure 2.** (**a**) PCoA ordination plot based on Bray-Curtis similarities for comparison of water and fish microbiota from samples taken early (1/4–8 dph) and late (13–16 dph) in the experiment. One fish sample consists of five pooled cod larvae from one rearing bottle. (**b**) Bray-Curtis similarities within fish and water microbiota, and between fish and water microbiota. Bars indicate standard deviation. **Figure 2.** (**a**) PCoA ordination plot based on Bray-Curtis similarities for comparison of water and fish microbiota from samples taken early (1/4–8 dph) and late (13–16 dph) in the experiment. One fish sample consists of five pooled cod larvae from one rearing bottle. (**b**) Bray-Curtis similarities within fish and water microbiota, and between fish and water microbiota. Bars indicate standard deviation.

Curtis similarities (Figure 2b). The PCoA plot also indicated that both the water and fish microbiota changed over time. There were significant differences between the water and fish microbiota both early (1/4–8 dph) and late (13–16 dph) in the experiment (one-way PERMANOVA, *p* = 0.04 and 0.01 for early and late, respectively) and also between water samples early and late in the experiment (one-way PERMANOVA, *p* = 0.01). There were no significant differences in the microbiota of fish samples early and late in the

The results from PERMANOVA analysis are reflected in the taxonomic composition of the microbial communities at the order level (Figure 3). The relative abundance of Vibrionales was more than 14 times higher in fish (up to 44%) than in water (≤3%). The relative abundances of Flavobacteriales in the water increased with time. The results from PERMANOVA analysis are reflected in the taxonomic composition of the microbial communities at the order level (Figure 3). The relative abundance of Vibrionales was more than 14 times higher in fish (up to 44%) than in water (≤3%). The relative abundances of Flavobacteriales in the water increased with time. *Microorganisms* **2022**, *9*, x FOR PEER REVIEW 7 of 18

**Figure 3.** Taxa summary of bacterial orders detected at a relative abundance higher than 0.1% in at least one sample, at 1 to 16 dph (D1–D16). Numbers indicate the rearing bottle sampled. One fish sample consists of five pooled larvae from one rearing bottle. **Figure 3.** Taxa summary of bacterial orders detected at a relative abundance higher than 0.1% in at least one sample, at 1 to 16 dph (D1–D16). Numbers indicate the rearing bottle sampled. One fish sample consists of five pooled larvae from one rearing bottle.

conventional and gnotobiotic larvae. However, the genes *G-protein-coupled receptor family C group 6 member A* (*gprc6a*) (involved in regulation of inflammation, metabolism and endocrine functions) and *rhamnose binding lectin* (*rbl*) (involved in innate immunity) were downregulated in both conventional and gnotobiotic larvae compared with germ-free larvae. The *zg16* and *zg16-like* genes (involved in innate immunity) were also downregulated in gnotobiotic larvae compared with germ-free larvae. Only one gene,

*lect2*, was upregulated in conventional larvae compared with germ-free larvae.

At 8 dph, no genes in the cod larvae were differentially expressed between the

For samples from 16 dph, 82 genes were downregulated and 97 were upregulated in conventional larvae compared with germ-free larvae. In gnotobiotic larvae, only 23 were downregulated and none were upregulated compared with germ-free larvae (Supplementary Tables S1–S3). Gnotobiotic and germ-free larvae generally showed similar expression profiles. Of the genes that were upregulated in conventional compared with germ-free larvae, 74% were also upregulated compared with gnotobiotic larvae (Supplementary Figure S4). Many of the genes that were upregulated in conventional compared with germ-free larvae are involved in innate immune responses and linked to signalling and glucose transport. Examples for immunity are, e.g., *interleukin 8* (*cxcl8*), *leukocyte cell-derived chemotaxin 2* (*lect2*) and *interleukin-1 receptor-activated kinase* (*irak1*), and for signalling/transport are, e.g., *solute carrier family 2 facilitated glucose transporter member*

*3.3. Gene Expression in Atlantic Cod Larvae*

OTU 1 (*Polaribacter*, Flavobacteriales order) and OTU 2 (*Vibrio*) were the most abundant OTUs in the dataset. Both OTUs were present in both fish and water samples, but OTU 1 was far more abundant in the water than in the fish (on average 31% and 4% of the reads, respectively). OTU 2 had a higher relative abundance in the fish than in the water (average 21% and 0.4%, respectively). The third most abundant OTU was OTU 4 (*Colwellia*), which was more abundant early in the experiment (1–8 dph), both for fish (average 9.2%) and water (average 14.2%) samples, than later in the experiment (13–16 dph) (average 1.1% and 2.8% in fish and water samples, respectively).

#### *3.3. Gene Expression in Atlantic Cod Larvae*

At 8 dph, no genes in the cod larvae were differentially expressed between the treatments. For 13 dph samples, there were still no differentially expressed genes between conventional and gnotobiotic larvae. However, the genes *G-protein-coupled receptor family C group 6 member A* (*gprc6a*) (involved in regulation of inflammation, metabolism and endocrine functions) and *rhamnose binding lectin* (*rbl*) (involved in innate immunity) were downregulated in both conventional and gnotobiotic larvae compared with germ-free larvae. The *zg16* and *zg16-like* genes (involved in innate immunity) were also downregulated in gnotobiotic larvae compared with germ-free larvae. Only one gene, *lect2*, was upregulated in conventional larvae compared with germ-free larvae.

For samples from 16 dph, 82 genes were downregulated and 97 were upregulated in conventional larvae compared with germ-free larvae. In gnotobiotic larvae, only 23 were downregulated and none were upregulated compared with germ-free larvae (Supplementary Tables S1–S3). Gnotobiotic and germ-free larvae generally showed similar expression profiles. Of the genes that were upregulated in conventional compared with germ-free larvae, 74% were also upregulated compared with gnotobiotic larvae (Supplementary Figure S4). Many of the genes that were upregulated in conventional compared with germ-free larvae are involved in innate immune responses and linked to signalling and glucose transport. Examples for immunity are, e.g., *interleukin 8* (*cxcl8*), *leukocyte cell-derived chemotaxin 2* (*lect2*) and *interleukin-1 receptor-activated kinase* (*irak1*), and for signalling/transport are, e.g., *solute carrier family 2 facilitated glucose transporter member 11-like* (*slc2a11*) and *solute carrier family 2 facilitated glucose transporter member 4-like* (*slc2a4*). Biological term enrichment analysis and Gene Ontology term (GO term) annotation in DAVID showed that 15 GO terms were enriched in conventional fish compared to germfree fish (Figure 4). Several of the enriched GO terms, including the one with the highest number of genes ("regulation of nucleobase-containing metabolic process"), were related to growth and cell division. Other enriched GO terms were related to signalling and cell communication. Most of these GO terms were also enriched compared with gnotobiotic larvae (Supplementary Table S4). The KEGG pathway "bacterial invasion of epithelial cells" was enriched in conventional larvae compared with both gnotobiotic and germ-free fish.

However, 15 annotated genes had significantly lower expression in both conventional and gnotobiotic larvae compared with germ-free larvae. Interestingly, nine of these genes were involved in innate immune responses: *eosinophil peroxidase* (*epx*), *rbl*, *zymogen granule membrane protein 16* (*zgp16*), *myeloperoxidase precursor* (*mpo*), *Cytochrome b-245 heavy chain-like* (*cybb*), *immune-responsive gene 1 protein-like* (*Irg1*), *fish egg lectin* (*fel*), *N-acetylmuramoyl-Lalanine amidase-like* (*pglyrp1*) and *transmembrane protease serine 9-like* (*tmprss9*). Thus, the presence of bacteria, both as complex communities and simple gnotobiotic associations, downregulated some of the innate immune responses in the cod larvae.

Far more GO terms differed between germ-free and conventional larvae than between germ-free and gnotobiotic larvae (95 and 11, respectively). The most enriched GO terms in the germ-free larvae were "proteolysis", "negative regulation and regulation of metabolic processes", "signal transduction" and "adhesion and cell death" (Figure 5).

fish.

*11-like* (*slc2a11*) and *solute carrier family 2 facilitated glucose transporter member 4-like* (*slc2a4*). Biological term enrichment analysis and Gene Ontology term (GO term) annotation in DAVID showed that 15 GO terms were enriched in conventional fish compared to germfree fish (Figure 4). Several of the enriched GO terms, including the one with the highest number of genes ("regulation of nucleobase-containing metabolic process"), were related to growth and cell division. Other enriched GO terms were related to signalling and cell communication. Most of these GO terms were also enriched compared with gnotobiotic larvae (Supplementary Table S4). The KEGG pathway "bacterial invasion of epithelial cells" was enriched in conventional larvae compared with both gnotobiotic and germ-free

**Figure 4.** Significantly enriched Gene Ontology (GO) terms in conventional fish compared with germ-free fish. **Figure 4.** Significantly enriched Gene Ontology (GO) terms in conventional fish compared with germ-free fish. GO terms in the germ-free larvae were "proteolysis", "negative regulation and regulation of metabolic processes", "signal transduction" and "adhesion and cell death" (Figure 5).

**Figure 5.** Enriched Gene Ontology (GO) terms (including 10 or more of the 52 input genes) in germfree fish compared to conventional fish. **Figure 5.** Enriched Gene Ontology (GO) terms (including 10 or more of the 52 input genes) in germ-free fish compared to conventional fish.

Of the 21 genes included in the GO term "proteolysis", 5 were involved in the KEGG

Comparison of the midgut ultrastructure, including tight junctions, microvilli disruption/damage, intercellular space and vacuoles, showed no significant differences between germ-free and conventional cod larvae. However, the mitochondria in germ-free cod larvae were distorted. The outer membrane showed discontinuities or was missing, and structures of the Christae were reduced or hardly visible (Figure 6). In contrast, the mitochondria in conventional cod larvae showed clear Christae and a clear double

membrane.

*3.4. Ultrastructure and Morphometric Analysis of the Intestinal Tissue*

Of the 21 genes included in the GO term "proteolysis", 5 were involved in the KEGG pathway of "protein digestion and absorption". This indicates that a large fraction of the "proteolysis" GO term is related to the cod larvae's digestion.

#### *3.4. Ultrastructure and Morphometric Analysis of the Intestinal Tissue*

Comparison of the midgut ultrastructure, including tight junctions, microvilli disruption/damage, intercellular space and vacuoles, showed no significant differences between germ-free and conventional cod larvae. However, the mitochondria in germ-free cod larvae were distorted. The outer membrane showed discontinuities or was missing, and structures of the Christae were reduced or hardly visible (Figure 6). In contrast, the mitochondria in conventional cod larvae showed clear Christae and a clear double membrane. *Microorganisms* **2022**, *9*, x FOR PEER REVIEW 10 of 18

larvae showed clear Christae and a double membrane. Scale bars are 1 µm.

**Figure 6.** Transmission electron microscopy micrographs of mitochondria from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. (**a**) Mitochondria in germ-free cod larvae were distorted, outer membrane showed discontinuities or was often missing and structures of Christae were reduced or hardly visible. (**b**) Mitochondria in conventional cod larvae showed clear Christae and a double membrane. Scale bars are 1 µm. **Figure 6.** Transmission electron microscopy micrographs of mitochondria from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. (**a**) Mitochondria in germ-free cod larvae were distorted, outer membrane showed discontinuities or was often missing and structures of Christae were reduced or hardly visible. (**b**) Mitochondria in conventional cod larvae showed clear Christae and a double membrane. Scale bars are 1 µm. (**a**) (**b**) **Figure 6.** Transmission electron microscopy micrographs of mitochondria from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. (**a**) Mitochondria in germ-free cod larvae were distorted, outer membrane showed discontinuities or was often missing and structures of Christae were reduced or hardly visible. (**b**) Mitochondria in conventional cod

The microvillous brush borders in the midgut of cod larvae at 16 dph, reared under

germ-free as well as conventional conditions, were well-defined and regular. Interestingly, germ-free cod larvae had significantly longer and significantly thicker microvilli than conventional cod larvae (for both analyses, Welch test, *p* ≤ 0.001, Figure 7). Moreover, the abundance of microvilli (µm–2) in the midgut was significantly lower in the germ-free larvae than in the conventional ones (Welch test, *p* ≤ 0.001, insets in Figure 7). The microvillous brush borders in the midgut of cod larvae at 16 dph, reared under germ-free as well as conventional conditions, were well-defined and regular. Interestingly, germ-free cod larvae had significantly longer and significantly thicker microvilli than conventional cod larvae (for both analyses, Welch test, *p* ≤ 0.001, Figure 7). Moreover, the abundance of microvilli (µm–2) in the midgut was significantly lower in the germ-free larvae than in the conventional ones (Welch test, *p* ≤ 0.001, insets in Figure 7). The microvillous brush borders in the midgut of cod larvae at 16 dph, reared under germ-free as well as conventional conditions, were well-defined and regular. Interestingly, germ-free cod larvae had significantly longer and significantly thicker microvilli than conventional cod larvae (for both analyses, Welch test, *p* ≤ 0.001, Figure 7). Moreover, the abundance of microvilli (µm–2) in the midgut was significantly lower in the germ-free larvae than in the conventional ones (Welch test, *p* ≤ 0.001, insets in Figure 7).

(**a**) (**b**)

**Figure 7.** Transmission electron microscopy micrographs from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. Microvilli in conventional cod larvae were shorter and thinner compared to the microvilli in germ-free cod larvae. Insets in (**a**,**b**): Transverse section of brush border—microvilli were closer to each other in the conventional than in the germ-free cod larvae. Scale bars are 1 µm and of insets 500 nm. **Figure 7.** Transmission electron microscopy micrographs from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. Microvilli in conventional cod larvae were shorter and thinner compared to the microvilli in germ-free cod larvae. Insets in (**a**,**b**): Transverse section of brush border—microvilli were closer to each other in the conventional than in the germ-free cod larvae. Scale bars are 1 µm and of insets 500 nm. Microvilli were significantly closer to each other in the conventional than in the germ-**Figure 7.** Transmission electron microscopy micrographs from the midgut of cod larvae at 16 dph reared under (**a**) germ-free and (**b**) conventional conditions. Microvilli in conventional cod larvae were shorter and thinner compared to the microvilli in germ-free cod larvae. Insets in (**a**,**b**): Transverse section of brush border—microvilli were closer to each other in the conventional than in the germ-free cod larvae. Scale bars are 1 µm and of insets 500 nm.

Microvilli were significantly closer to each other in the conventional than in the germ-

free cod larvae. The morphometric measurements for the microvillous length, abundance

and diameter are summarised in Table 1.

and diameter are summarised in Table 1.

Microvilli were significantly closer to each other in the conventional than in the germfree cod larvae. The morphometric measurements for the microvillous length, abundance and diameter are summarised in Table 1.

**Table 1.** Average length, diameter and abundance of microvilli in germ-free and conventional cod larvae at 16 dph. *n* = 3 cod larvae per treatment were analysed.

#### **4. Discussion**

In this study, gnotobiotic husbandry of Atlantic cod larvae was used as a tool to study host–microbe interactions. Cod larvae were reared under three different treatments: germ-free, gnotobiotic (*Vibrio gallicus* and *Microbacterium* added) and conventional. The only difference between the cod larvae in the conventional treatment and those in the gnotobiotic and germ-free treatments was that the microbial conditions in the conventional treatment were uncontrolled.

For the transcriptomic analysis, the major differences between the treatments were observed at 16 dph, and therefore the discussion is focused on this time point. At 16 dph, the gut of the larvae is larger and more developed, and thus the number of niches available for bacteria may be higher than at earlier life stages, and this may allow more bacteria to coexist through selection [40]. This was also reflected in our data, as the richness increased at this time point. Importantly, our analyses were based on pooled, homogenised whole fish, leading to conservative conclusions. Analyses at the organ level could possibly discover more differences in the gene expression patterns between the fish from the different treatments, also at earlier time points.

#### *4.1. Microbial Environments*

DGGE analyses of the microbiota from the gnotobiotic treatment showed that *V. gallicus* dominated both in water and fish samples, whereas *Microbacterium* was detectable at low levels in water samples but present in only some of the fish samples. This corroborated previous studies, where *V. gallicus* was found to adhere to and grow fast in mucus [25]. This could explain the higher abundance of *V. gallicus* compared to *Microbacterium* in the cod larvae in this experiment.

For the conventional treatment, we used 16S rDNA Illumina amplicon sequencing to characterise fish and water microbiota. The results show that the fish microbiota differed significantly from the water microbiota, corroborating earlier studies [33,41–43]. The water and larval microbiota were dominated by bacterial taxa considered to represent opportunistic, rapid-growing bacteria, such as Vibrionales, Actinomycetales, Alteromonadales and Flavobacteriales. This may be a result of r-selection in the water [44]. Pulses of organic matter originating from the addition of feed and fish defecation will create a high carrying capacity in the rearing bottles. When adding new water with lower carrying capacity to the rearing bottles, the microbe–microbe competition in the water is reduced, and this favoured growth of fast-growing, opportunistic species (r-strategists) [44]. Opportunistic bacteria could potentially be detrimental for the cod larvae, and this might be the reason for the lower survival observed in the conventional treatment than in the gnotobiotic and germ-free treatment. However, overall, the survival was good, and comparable to what is typically seen in first feeding experiments with cod larvae. In our dataset, 6 OTUs were classified as Vibrio. Using the Ribosomal database project (RDP) SeqMatch tool [30] to identify the most closely related type strains for each of them (Supplementary Table S5) showed that one OTU matched *V. campbellii* and one *V. anguillarum.* Both species are known fish pathogens [45].

Our results confirm that we had distinct microbial environments in the gnotobiotic and conventional treatments. The gnotobiotic treatment represents environmental conditions with non-detrimental host–microbe interactions. The conventional treatment, on the other hand, represents an environment characterised by the presence of opportunistic bacteria, which might have had detrimental effects on the fish. Thus, the conventional treatment does not represent a "natural" microbial environment, but rather a suboptimal microbial environment, including detrimental host–microbe interactions. This is reflected in the gene expression of the fish, and likely the reason why we see (1) increased expression of some genes related to inflammatory responses and oxidative stress, and (2) lower survival of conventionally reared fish than for the germ-free and gnotobiotic cod larvae.

#### *4.2. Presence of Bacteria Downregulates Host Responses Related to Nutrient Utilisation and Innate Immune Responses*

Members of the gut microbiota in other species are well-known to aid in the digestion of, e.g., complex carbohydrates [46] and proteins [47]. Since the gut of cod larvae is functionally immature at hatching [48,49], the enzymatic activity of bacteria may aid in digestion of the live feed organisms. Even if survival was very high for the germ-free cod, and the intake of feed appeared similar in all treatments, they gained less weight than the conventional larvae. The transcriptional responses in the germ-free larvae support that they had difficulties in digesting feed, as the most enriched GO term in germ-free larvae compared with conventional larvae was proteolysis. Of the 21 genes included in the GO term "proteolysis", 5 were involved in the KEGG pathway of "protein digestion and absorption". This indicates that a large fraction of the "proteolysis" GO term is related to the host's digestion. Interestingly, germ-free rats also have higher activities of digestive enzymes, such as amylase and lipase, in their intestine than conventional rats [50]. Transcriptional responses related to fasting were downregulated in the conventional larvae compared with germ-free larvae. Host responses related to fasting were also found in the zebrafish study by Rawls et al. [16].

Similar to findings in germ-free rodents [51,52], the germ-free cod larvae showed significantly longer microvilli in the midgut than the conventional larvae (Figure 7, Table 1). This might be related to reduced renewal of the intestinal epithelium in these animals [52]. Similarly, albeit at a larger physiological scale, Willing et al. [53] found that germ-free pigs have longer villi and shorter crypts in their distal intestine, and that the shortening observed after colonisation was associated with increased cell turnover. They hypothesise that commensal bacteria contribute to enterocyte turnover through induction of inflammatory responses and cell apoptosis. In addition to the increased microvilli length, germ-free cod showed less microvilli per µm2 in the midgut than conventional reared larvae. Microvilli not only increase the cellular surface area for absorption of nutrients, they also increase the number of digestive enzymes present on the cell surface. Thus, the lower density of microvilli in the germ-free larvae could explain the apparent reduced nutrient uptake in the larvae. The distorted mitochondria observed in the germ-free cod larvae (Figure 6) also support the hypothesis of a physiological starvation state in these larvae. Hailey et al. [54] demonstrated how lipids from the mitochondrial membrane are utilised in the biogenesis of autophagosomes under starvation conditions.

The gene expression analysis indicates that certain elements of the innate immune system of cod larvae are "turned on" in the larval stage of the fish but are subsequently regulated by host–microbiota interactions. Solbakken et al. [22] described how infection by a pathogenic bacterium, *Francisella noatunensis*, dampens the intracellular immune response to allow intracellular persistence of the pathogen. Our results indicate that generation of reactive oxygen species (ROS) was lower in conventional and gnotobiotic larvae than in germ-free larvae. ROS generation is tightly linked to mitochondrial metabolism, and thus could be elevated in the disintegrating mitochondria, as observed by electron microscopy

analysis (Figure 6). ROS is also produced during recognition of non-self-substances and the immune response processes [55]. The main enriched GO terms in germ-free larvae were related to defence responses and responses to reactive oxygen species. This presumably higher ROS activity despite no bacteria present might be a consequence of the unique immune system of cod, or because phagocytes of germ-free larvae are activated by components of the live feed (such as algae). During infection by *F. noatunensis,* downregulation of ROS production in phagosomes is hypothesised to contribute to this pathogen's survival [22]. Another indication of high ROS generation in the germ-free cod was high expression of *immunoresponsive gene 1* (*irg1*). Irg1 controls macrophage function, by regulating metabolic pathways leading to increased mitochondrial ROS production that aids bacterial killing. In zebrafish, this gene was upregulated by bacterial infection [56], whereas in our experiment, the expression was lower in the fish exposed to bacteria.

Transcripts for key proteins involved in recognition of both peptidoglycans (PGLYRP) and lectins were also downregulated by the bacteria present in this study. Pglyrp proteins have been found to be expressed in zebrafish eggs, developing embryos and adult tissues that are in contact with the environment, and there are indications that they have an important role in the defence against bacteria in young fish [57,58]. Both rbl (Rhamnose-binding lectin) and fel (fish egg lectin) are known to enhance phagocytosis, and the expression of the genes is normally upregulated when the host is exposed to potentially harmful bacteria [59,60]. However, Thongda et al. [59] point to studies in catfish where *rbl* is highly upregulated by short-term fasting, indicating a link between the feeding status and the immune function. Thus, the downregulation of *rbl* in the gnotobiotic and conventional cod might be coupled to the observed transcriptional responses to fasting. Thus, enforcing the indications that the germ-free fish do not digest their food as well as the gnotobiotic and conventional fish.

It has previously been shown that the bacteria used in the gnotobiotic treatment are non-detrimental and may improve survival of cod larvae [25]. Thus, the downregulation of immune responses due to these bacteria may be a way of inducing tolerance to the colonising microbiota. However, similar downregulation observed in conventional cod larvae, where the microbiota was shown to have a slight but significantly negative effect on the survival, suggests that the immune system of the larvae is not capable of distinguishing friend from foe at this early life stage.

Several of the immune-related genes found to be downregulated in this experiment are involved in processes that have been reported to be upregulated by the presence of bacteria in gnotobiotic zebrafish and stickleback [16,61,62]. This illustrates how early host– microbiota responses differ across teleost species. In contrast to stickleback and zebrafish, Atlantic cod are heavily reliant on the innate immune response due to the loss of MHC-2, which is critical for initiation of antigen-specific immune response. In addition, cod toll-like receptor families (TLRs) have undergone genetic deletions and subsequent diversifications, possibly to compensate for the lack of the classical adaptive immunity [24,63]. This means that the type of receptors and downstream immune pathways will differ to some extent between cod and other vertebrates such as zebrafish and mice [64]. As an example, whereas LPS in high doses is lethal to mammals and zebrafish [61], Atlantic cod has a lower LPS response, with much higher LD50 values [64,65]. Analysis of the full genome revealed that Atlantic cod lacks TLR4, the mammalian LPS receptor that has a functional ortholog in zebrafish [21,66]. Despite the lack of a TLR4 ortholog, cod head kidney cells still reacted to LPS exposure by upregulation of certain immune and xenobiotic pathways [66].

#### *4.3. Presence of Bateria Induces Host Responses Related to Inflammatory Responses and Signalling*

Even though several genes related to the innate immune system seem to be downregulated by the presence of bacteria, other innate immune system responses were induced by the presence of bacteria. *Cxcl8* (*interleukin 8*) and *lect2* (*leukocyte-cell derived chemotaxin 2*) (upregulated in both 13 and 16 dph larvae) were upregulated in both conventional and gnotobiotic larvae (*cxcl8* just below the cut-off value of log2 0.8-fold change in gnotobiotic larvae). These genes that encode chemokines that attract neutrophils by chemotactic activity are induced by inflammatory stimuli caused by, e.g., microbial stress, and have been identified in teleosts earlier [67–70]. Both genes were also upregulated by in vitro LPS stimulation of cod cell cultures [66]. In fish, *lect2* is assumed to have an important role in the inflammatory response, promoting phagocytic activity of macrophages [71,72]. Since *lect2* was upregulated at both 13 and 16 dph, we suspect that this gene may be an early and important bacterial response in cod.

The KEGG pathway "bacterial invasion of epithelial cells" and GO terms related to signalling and signal transduction were enriched in conventional larvae compared with both gnotobiotic and germ-free larvae. This seems plausible, as compared to a gnotobiotic community consisting of only two probiotic candidates, the microbiota in the conventional treatment could give a higher invasion pressure, more host–microbe interactions due to higher species richness and thus more signalling and host responses. For example, a gene linked to the innate immune system, *irak1* (*interleukin-1 receptor-associated kinase 1*), was upregulated in conventional larvae compared with germ-free larvae. This gene plays a critical role in initiating the innate immune response against pathogens and has been shown to be upregulated in fish after pathogen challenge [73]. This implies that *irak1* might participate in antibacterial immunity. In gnotobiotic fish, the expression of *irak1* was at the same level as in germ-free fish, indicating that the probiotic candidates used in the gnotobiotic treatment were detected as non-pathogens by the fish.

*Slc2a11* was the most upregulated gene (log2 fold change of 2.3) in conventional compared with germ-free larvae, and it was also upregulated in conventional compared with gnotobiotic larvae. This gene encodes a glucose transporter protein GLUT11, belonging to class II of these proteins [74]. In fish, glucose transport is important for several reasons: blood glucose levels change rapidly in response to environmental disturbances, increased plasma glucose levels may be an indicator of stress and glucose intolerance has been documented in fish [75]. The expression of another GLUT protein, GLUT4, belonging to class I of these proteins and encoded by the *slc2a4* gene, was also upregulated in conventional fish compared to the germ-free fish. This is the only insulin-sensitive member of class I, and it is expressed in insulin-sensitive tissues, such as heart, muscle and adipose tissue in cod [76]. These findings indicate that the microbial community in the conventional rearing bottles induced more stress on the fish than in the gnotobiotic treatment. This is in line with the fact that potential detrimental bacteria were found in the conventional rearing bottles, and that the bacterial strains added to the gnotobiotic flasks had been characterised as probiotic candidates in an earlier study [25].

### **5. Conclusions**

To conclude, our results indicate that bacteria actively downregulate certain cod larvae immune responses, facilitating bacterial colonisation of mucosal surfaces. This concept of "downregulation" contrasts previous findings in zebrafish and stickleback and emphasizes the role of evolutionary history, highlighting the need to study host–microbe interactions in several teleost species. Similar to what has been shown in studies with other vertebrates, bacterial colonisation improved the nutritional state of the cod larvae, evident at both transcriptional and micromorphological levels and materialised as differences in growth rate. This study illustrates the dynamics between water- and host-associated microbiota and increases our insight into how Atlantic cod larvae respond physiologically and transcriptionally to bacterial colonisation.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10 .3390/microorganisms10010024/s1, Figure S1: Kaplan–Meier survival curves, Figure S2: Dry weight of cod larvae, Figure S3: DGGE gel, Figure S4: Venn diagram, Table S1: List of down-regulated genes (cut-off value of log2 fold change > 0.8) in conventional cod larvae compared with germ-free cod larvae at 16 dph. Table S2: List of up-regulated genes (cut-off value of log2 fold change > 0.8) in conventional cod larvae compared with germ-free cod larvae at 16 dph. Table S3: List of downregulated genes (cut-off value of log2 fold change > 0.8) in gnotobiotic cod larvae compared with germ-free cod larvae at 16 dph. Table S4: Enriched GO terms in conventional larvae compared with both germ-free and gnotobiotic cod larvae. Table S5: Taxonomy and most similar "type strains" for the six OTUs classified as Vibrio as inferred from the Ribosomal Database Project (RDP). Text S1: Verification of germ-free conditions, Text S2: Live feed and bacterial cultures, Text S3: Sampling procedures, Text S4: Statistical analyses.

**Author Contributions:** Conceptualisation, T.F. and O.V.; formal analysis, T.F., R.I.V., B.L. and P.W.; investigation, T.F., R.I.V. and B.L.; writing—original draft preparation, R.I.V., T.F. and B.L.; writing review and editing, R.I.V., T.F., B.L., O.V., I.B., Y.O. and P.W.; visualisation, R.I.V.; supervision, O.V. and I.B.; project administration, O.V., T.F. and R.I.V.; funding acquisition, O.V., T.F. and B.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was performed with funding from the Faculty of Natural Sciences at the NTNU— Norwegian University of Science and Technology, the Norwegian Research Council (HAVBRUK, 233865), and a People Marie Curie Actions International Incoming Fellowship FP7-PEOPLE-2013- IIF (#625655).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Norwegian Animal Facility Establishment 154 – NTNU – Sealab (protocol code 4174, 5 March 2012).

**Data Availability Statement:** Sequencing data are stored at European Nucleotide Archive, accession numbers ERS8484975-ERS8484994.

**Acknowledgments:** We thank Mari-Ann Østensen for invaluable help in conducting this experiment. Thanks to Tora Bardal and Rolf Erik Olsen for valuable training and input. The TEM work was performed at the Cellular and Molecular Imaging Core Facility (CMIC), Norwegian University of Science and Technology (NTNU). CMIC is funded by the Faculty of Medicine at NTNU and the Central Norway Regional Health Authority.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

#### **References**