# Screening New Xylanase Biocatalysts from the Mangrove Soil Diversity

Corinne Ivaldi, Mariane Daou, Laurent Vallon, Alexandra Bisotto, Mireille Haon, Sona Garajova, Emmanuel Bertrand, Craig B Faulds, Giuliano Sciara, Adrien Jacotot, et al.

► To cite this version:

Corinne Ivaldi, Mariane Daou, Laurent Vallon, Alexandra Bisotto, Mireille Haon, et al. Screening New Xylanase Biocatalysts from the Mangrove Soil Diversity. Microorganisms, 2021, 9 (7), pp.1-16. 10.3390/microorganisms9071484. insu-03293424

HAL Id: insu-03293424

https://insu.hal.science/insu-03293424

Submitted on 21 Jul 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Screening New Xylanase Biocatalysts from the Mangrove Soil Diversity

Corinne Ivaldi<fnref n="1" />, Mariane Daou <fnref n="2" />,<fnref n="3" /> ID, Laurent Vallon <fnref n="4" />, Alexandra Bisotto <fnref n="2" />, Mireille Haon <fnref n="2" />, Sona Garajova <fnref n="2" />, Emmanuel Bertrand <fnref n="2" /> ID, Craig B. Faulds <fnref n="2" /> ID, Giuliano Sciara <fnref n="2" />, Adrien Jacotot <fnref n="5" />,<fnref n="6" />,<fnref n="7" />, Cyril Marchand <fnref n="5" />,<fnref n="6" />, Mylène Hugoni <fnref n="4" />, Harivony Rakotoarivonina <fnref n="1" />, Marie-Noëlle Rosso <fnref n="2" /> ID, Caroline Rémond <fnref n="1" /> ID, Patricia Luis <fnref n="4" /> ID and Eric Record <fnref n="2" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Ivaldi, C.; Daou, M.; Vallon, L.; Bisotto, A.; Haon, M.; Garajova, S.; Bertrand, E.; Faulds, C.B.; Sciara, G.; Jacotot, A.; et al. Screening New Xylanase Biocatalysts from the Mangrove Soil Diversity.
*Microorganisms* **2021**, *9*, 1484. https://doi.org/10.3390/microorganisms9071484

Academic Editors: Dietmar Haltrich and Daniel Kracher

Received: 17 June 2021

Accepted: 9 July 2021

Published: 12 July 2021

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

<figure><img src="image_4.png" /></figure>

**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/).

**Abstract:** Mangrove sediments from New Caledonia were screened for xylanase sequences. One enzyme was selected and characterized both biochemically and for its industrial potential. Using a specific cDNA amplification method coupled with a MiSeq sequencing approach, the diversity of expressed genes encoding GH11 xylanases was investigated beneath *Avicenia marina* and *Rhizophora styloosa* trees during the wet and dry seasons and at two different sediment depths. GH11 xylanase diversity varied more according to tree species and season, than with respect to depth. One complete cDNA was selected (OFU29) and expressed in *Pichia pastoris*. The corresponding enzyme (called Xyn11-29) was biochemically characterized, revealing an optimal activity at 40–50 °C and at a pH of 5.5. Xyn11-29 was stable for 48 h at 35 °C, with a half-life of 1 h at 40 °C and in the pH range of 5.5–6. Xyn11-29 exhibited a high hydrolysis capacity on destarched wheat bran, with 40% and 16% of xylose and arabinose released after 24 h hydrolysis. Its activity on wheat straw was lower, with a release of 2.8% and 6.9% of xylose and arabinose, respectively. As the protein was isolated from mangrove sediments, the effect of sea salt on its activity was studied and discussed.

**Keywords:** lignocellulose degrading enzymes; xylanases; heterologous expression; biomass degradation; marine fungus; mangrove; salt adaptation

## 1. Introduction

The mangrove is a unique ecosystem located at an interface between terrestrial, estuarine, and near-shore marine ecosystems. It is exposed to tides and has to live with high concentrations of salts due to the evaporation of saline water. Mangrove trees have developed unique adaptations to survive high levels of salinity either by excreting salts through their leaves or by excluding them at the root surface with a filtration system. One of the main environmental factors controlling the distribution of vegetation appears to be sediment salinity, which is partly controlled by the position of the stand in the intertidal

<footnote n="1">¹ INRAE, FARE, UMR A 614, Chaire AFERE, Université de Reims Champagne Ardenne, 51097 Reims, France; corinne.ivaldi@univ-reims.fr (C.I.); Harivony.rakotoarivonina@univ-reims.fr (H.R.); caroline.remond@univ-reims.fr (C.R.)</footnote>

<footnote n="2">² INRAE, UMR1163, Biodiversité et Biotechnologie Fongiques, Aix-Marseille Université, 13288 Marseille, France; mariane.daou@ku.ac.ae (M.D.); alexandra.bisotto@inrae.fr (A.B.); mireille.haon@inrae.fr (M.H.); sona.garajova@inrae.fr (S.G.); emmanuel.bertrand@univ-amu.fr (E.B.); craig.faulds@univ-amu.fr (C.B.F.); giuliano.sciara@inrae.fr (G.S.); marie-noelle.rosso@inrae.fr (M.-N.R.)</footnote>

<footnote n="3">³ Department of Chemistry, Khalifa University, Abu Dhabi 127788, United Arab Emirates</footnote>

<footnote n="4">⁴ CNRS, INRAE, VetAgro Sup, UMR Ecologie Microbienne, Université Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne, France; laurent.vallon@univ-lyon1.fr (L.V.); mylene.hugoni@univ-lyon1.fr (M.H.); patricia.luis@univ-lyon1.fr (P.L.)</footnote>

<footnote n="5">⁵ Institut de Recherche pour le Développement (IRD), IMPMC, UPMC, CNRS, MNHN, 98851 Noumea, New Caledonia, France; adrien.jacotot@unc.nc (A.J.); cyril.marchand@unc.nc (C.M.)</footnote>

<footnote n="6">⁶ ISEA, Université de la Nouvelle-Calédonie, EA 7484, 8 BPR4, 98851 Noumea, New Caledonia, France</footnote>

<footnote n="7">⁷ CNRS, BRGM, ISTO, UMR 7327, Université d'Orléans, 45071 Orléans, France</footnote>

<footnote>* Correspondence: eric.record@inrae.fr</footnote>
zone [1]. In addition, mangrove trees in areas with high salinity are usually characterized
by a lower productivity and, thus, lower organic accumulation in their soils, modifying
the sediment biogeochemistry. In this hostile environment, mangrove fungi contribute
to the degradation of particulate organic matter into dissolved organic matter, including
the lignocellulosic biomass [2]. In New Caledonia, mangrove zonation is typical of the
semi-arid climate, with *Rhizophora stylosa* colonizing the lower intertidal zone and *Avicennia marina* growing at higher elevation with higher pore-water salinity and lower organic
accumulation in sediments [3]. The fungi abundantly present within these mangroves
constitute a very diverse ecological group (~100 to 1000 different taxa depending on the
sample, many of them undescribed) characterized by a strong host specificity, the fungal
communities associated with *A. marina* being very different from those associated with
*R. stylosa* [4]. In a recent study, the prokaryotic and fungal communities in New Caledonian
mangrove sediments were analyzed with respect to depth, vegetation cover, and season
using a metabarcoding approach [5]. The prokaryotic community appeared to be exclu-
sively shaped by sediment depth, while the fungal community was evenly distributed
according to all the above criteria. In parallel, the functional diversity of the same New
Caledonia mangrove sediments was examined by observing the distribution of fungal
dye-decolorizing peroxidases (DyPs) as model lignin-modifying enzymes, thought to play
a role in lignocellulose biomass degradation [6,7]. Using a functional metabarcoding ap-
proach, the diversity of expressed genes encoding fungal DyPs was investigated in surface
and deeper sediments. The highest DyP diversity was observed for the surface layers
beneath the *Rhizophora stylosa* area during the wet season. With regard to enzyme diversity,
fungal adaptation to marine conditions, especially to high salt concentrations, is of great
interest and was characterized for the mangrove fungus *Pestalotiopsis* sp. using a proteomic
approach. This study revealed that the composition of secreted lignocellulolytic enzymes
was strongly modified in the presence of sea salt, with increased secretion of xylanases and
cellulases and lower production of oxidoreductases [8].

Xylanases (EC 3.2.1.8) are glycosidases that degrade xylans into xylooligosaccharide and xylose and thus help break down hemicellulose, one of the major components of plant cell walls [9]. Xylanases are classified in the CAZy database in GH families 5, 7, 8, 10, 11, and 43, with a predominance of GH10 and 11 [10]. Microbial xylanases are used in various industrial applications such as in the food industry (production of wine and beer, fruit juice, bread), for animal feed, biofuel production, and in the pulp and paper industry [9]. These industrial processes are conducted in a range of temperatures and pHs. In this context, xylanases were screened from marine environments, as they operate at different pH and temperature ranges from their terrestrial enzyme homologues and are also active in saline conditions [11]. Several studies have focused on xylanase activity screening to find new microbial activities and enzymes from seawater and mangrove forests, sediments, sponges and algae [11]. In these marine environments, fungal marine-derived xylanases were found to act at neutral to alkaline pH [9]. In saline conditions, xylanase XynMF13A from the mangrove fungus *Phoma* sp. MF13, heterologously produced in *Pichia pastoris*, was shown to have an optimal activity in 0.5 M NaCl and to keep 53% of its initial activity in 4 M NaCl [12].

In the present work, our objective was to evaluate the diversity of expressed genes encoding xylanases from mangrove sediments and to obtain the corresponding full-length cDNAs. From this data, one recombinant xylanase, Xyn11-29, was heterologously produced in *P. pastoris*, biochemically characterized, and its biochemical properties compared with those reported in the literature. This new xylanase was tested for the hydrolysis of lignocellulosic co-products, and its hydrolytic efficiency was compared with that of another xylanase, GH11.
## 2. Materials and Methods

### 2.1. Sediment Sampling, RNA Extraction and cDNA Synthesis

Sediment samples were collected from a mangrove wetland located in Saint Vincent Bay (21°55′58″ S, 166°4′30″ E) on the west coast of New Caledonia. As previously described [3], sediment samples were collected in three independent 10 m² plots (A, B, C) located 50 m apart and defined in *Avicennia marina* (A) and *Rhizophora stylosa* (R) pristine areas. Three sediment cores (50 cm deep) were collected in 2016 at low tide with a stainless-steel corer (diameter 8 cm) in each 10 m² plot during the wet (March) and dry (November) seasons. Oxic (0–10 cm deep) and anoxic (40–50 cm deep) fractions were collected from each core and a single composite sample per fraction and per plot (A, B or C) for each tree area (*A. marina* and *R. stylosa*) was prepared by mixing equal amounts of sediments. Per season (March or November), a total of 12 different composite samples were thus obtained and separately analyzed: (i) R1A, R1B, R1C and R2A, R2B, R2C corresponding, respectively, to the oxic and anoxic fractions from the three plots (A, B, C) designated in the *R. stylosa* area (R) and (ii) A1A, A1B, A1C and A2A, A2B, A2C corresponding, respectively, to oxic and anoxic fractions of the three plots (A, B or C) localized in the *A. marina* area (A). All composite samples were frozen and kept at –70 °C until use.

Total RNA was extracted from 8–12 g of each composite sample using the RNeasy PowerSoil Total RNA kit, as previously described in Ben-Ayed et al. [6]. After a quality and quantity check, cDNA synthesis was conducted on 500 ng of total RNA using the Mint-2 cDNA synthesis kit (Evrogen, Moscow, Russia). Resulting cDNAs were then purified and used as templates to specifically capture full-length transcripts encoding fungal GH11 xylanases as previously described [13].

### 2.2. High-Throughput Sequencing of Fungal GH11 Expressed Genes and Full-Length Xyn11-29 Sequence Retrieval

Captured cDNAs from each sediment sample (12 per season) were first used as templates to specifically amplify fragments of expressed fungal GH11 xylanase-encoding genes using tagged degenerate primers fungGH11-F (5'-Tag-GGVAAGGGITGGAAYCCNGG-3') and fungGH11-R (5'-Tag-TGKCGRACIGACCARTAYTG-3') (Supplementary Table S1, [14]). All PCR amplifications were performed in triplicate in 25 µL reaction mixtures containing 2.5 µL of 10X buffer (Invitrogen, ThermoFisher Scientific, Waltham, MA, USA), 0.75 µL of MgCl₂ (50 mM), 2.5 µL of dNTPs (2 mM each), 1 µL of each primer (20 µM, Invitrogen), 0.3 µL of bovine serum albumin (20 mg mL⁻¹), 0.1 µL of Taq DNA polymerase (5U µL⁻¹, Invitrogen), and 20 ng of cDNA. The cycling conditions were 3 min at 94 °C followed by 35 cycles (45 s at 94 °C, 45 s at 50 °C and 45 s at 72 °C) and 10 min at 72 °C. Control reactions without nucleic acid were systematically run in parallel. Amplicons from the three independent PCR reactions were pooled and purified using the Agencourt AMPure XP Kit (Beckman Coulter Diagnostics, Brea, CA, USA) and quantified by fluorometry using a Qubit 2.0 fluorimeter and the Qubit dsDNA H assay kit (Invitrogen). Per season (March or November), an equimolar mix of the tagged PCR products obtained for the 12 different sediment samples was prepared and sequenced by FASTERIS (FASTERIS, Plan-les-Ouates, Switzerland) on an Illumina MiSeq sequencer (2 × 250 bp).

In parallel, a plasmid library containing full-length GH11 cDNAs captured from the 24 sediment samples was constructed as previously described [13], using the pDR196 E. coli/S. cerevisiae shuttle expression vector. The full-length cDNA sequence encoding the Xyn11-29 was retrieved from this plasmid library by PCR using two primers pairs (PMA1: 5'-CTCTCTTTTATACACACATTC-3'/Xyn11-29-F: 5'-TCCGTCAGCTGCAACGGAGCG-3 and ADH: CATAAATCATAAGAAATTCGC/Xyn11-29-R: 5'-CGCTCCGTTGCAGCTGACGGA-3') specifically targeting both plasmid and OFU29 sequences.
## 2.3. Bioinformatic Analysis and Statistics

As described in Ben-Ayed et al. [6], we used the pipeline Frogs to analyze Miseq sequences. Briefly, paired-end sequences were merged and demultiplexed to perform denoising procedures (presenting no ambiguous bases and in the range 230–350 bp). Trimmed sequences were clustered using Swar. Chimeras and scarce sequences (less than 0.005% of total abundance) were discarded. A normalization threshold of 10,311 sequences was applied to compare samples. The effects of environmental factors (tree, season, and depth) on the composition of expressed fungal genes encoding GH11 xylanases were tested using non-parametric permutation-based multi-variate analysis of variance (PERMANOVA, adonis function; [15] based on abundance dissimilarity (Bray–Curtis) matrices). This analysis was performed using the VEGAN package (http://cran.r-project.org, accessed on 19 April 2021) in R. The sequence data generated in this study were deposited in the EMBL-ENA public database (PRJEB43346 for the first sampling campaign (March) dataset and PRJEB43343 for the second sampling campaign (November) dataset).

The ProtParam tool (http://web.expasy.org/protparam/, accessed on 5 May 2021) was used to predict the theoretical molecular mass, and molar extinction coefficient of Xyn11-29. The signal peptide sequence was predicted using the tool available at http://www.cbs.dtu.dk/services/SignalP-4.1/ (accessed on 5 May 2021). N-Glycosylation sites were predicted via the NetNGlyc 1.0 server (http://www.cbs.dtu.dk/services/NetNGlyc/, accessed on 5 May 2021). For sequence comparison, BlastP was used to search for sequences with similarity to Xyn11-29 at NCBI portal (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 5 May 2021) using the non-redundant protein sequence data. For the phylogenetic analysis, Muscle alignment was used at MEGA X (https://www.megasoftware.net/, accessed on 5 May 2021). The evolutionary history was inferred by using the maximum likelihood method and JTT matrix-based model under MEGA X [16,17]. The tree with the highest log likelihood (−5562.80) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial trees for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model and then selecting the topology with the highest log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. This analysis involved 51 amino acid sequences. There were, in all, 218 positions in the final dataset.

## 2.4. Strains for Cloning and Heterologous Expression

*Escherichia coli* strain DH5α (Promega, Charbonnières, France) was used for vector storage and propagation. *P. pastoris* strain X33 (Invitrogen) was used for the heterologous expression of the Xyn11-29 encoding synthetic gene after the optimization of codons (GenScript, Piscataway, NJ, USA).

## 2.5. Cloning of the Xyn11-29 cDNA in Pichia Pastoris and Screening of Transformants

Xyn11-29 was produced using the 3PE Platform (*P. pastoris* Protein Express: www.platform3pe.com/, accessed on 5 May 2021). The cDNA encoding Xyn11-29 was synthesized after codon optimization for *P. pastoris* (GeneArt, Regensburg, Germany) and inserted into the vector pPICZαA (Invitrogen) using *Xho*I and *Xba*I restriction sites in frame with the α-factor and the (His)6 tag at the C terminus of the recombinant protein.

The *Pme*I-linearized pPICZαA recombinant plasmid was inserted into *P. pastoris*-competent cells by electroporation. Zeocin-resistant transformants were then screened for protein production. *P. pastoris* strain X33 and the pPICZαA vector are components of the *P. pastoris* Easy Select Expression System (Invitrogen). Electrocompetent cell preparation, electroporation, and screening were carried out as previously described [18]. Protein production was confirmed by SDS-PAGE.
## 2.6. Production and Purification of Recombinant Xyn11-29

The best producing transformant was grown in 2 L of BMGY (10 g L⁻¹ glycerol, 10 g L⁻¹ yeast extract, 20 g L⁻¹ peptone, 3.4 g L⁻¹ YNB, 10 g L⁻¹ ammonium sulfate, 100 mM phosphate buffer pH 6 and 0.2 g L⁻¹ of biotin) in flasks shaken at 30 °C in an orbital shaker (200 g) for 16 h to an OD600 of 2–6. Cells were then pelleted by centrifuging at 3500× g for 5 min and resuspended in 100 mL of a BMMY medium (10 g L⁻¹ yeast extract, 20 g L⁻¹ peptone, 3.4 g L⁻¹ YNB, 10 g L⁻¹ ammonium sulfate, 100 mM phosphate buffer pH 6 and 0.2 g L⁻¹ of biotin). Xyn11-29 production was induced at 20 °C in an orbital shaker (200 g) by adding 3% (v/v) methanol daily for three days.

The supernatant (500 mL) was collected after harvesting cells by centrifugation at 3500× g for 5 min at 4 °C. After adjusting the pH to 7.8, the supernatant was filtered on 0.45 µm filters (Millipore, Molsheim, France) and loaded onto a 5 mL HisTrap HP column (GE healthcare, Buc, France) connected to an Akta Xpress system (GE healthcare). Prior to loading, the column was equilibrated with buffer A Tris-HCl 50 mM pH 7.8, NaCl 150 mM and imidazole 10 mM. The (His)6-tagged recombinant enzyme was eluted with buffer B Tris-HCl 50 mM pH 7.8, NaCl 150 mM and imidazole 500 mM. The fractions eluted containing the purified protein were pooled, concentrated with a 10 kDa Vivaspin concentrator unit (Sartorius, Palaiseau, France) and dialyzed against 50 mM sodium acetate buffer pH 5.2.

Protein concentration was determined using a Nanodrop ND-2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) by adsorption at 280 nm with theoretical molecular masses (22,807 Da) and molar extinction coefficients (66,935 M⁻¹ cm⁻¹) calculated from protein sequence using Expasy tools. A fraction of eluate was loaded onto 10% Tris-glycine precast SDS-PAGE (Bio-Rad, Marnes-la-Coquette, France) to check protein purity and integrity. The molecular mass under denaturing conditions was determined with a PageRuler Prestained Protein Ladder (Thermo Fisher Scientific).

## 2.7. Standard Conditions of Xylanase Activity

Xylanase activities were measured by the quantification of reducing sugar release during incubation in the presence of beechwood xylan (Roth, Karlsruhe, Germany) (0.5%, w/v) in 50 mM citrate phosphate buffer at pH 5.5 and a temperature of 40 °C. Reducing sugars were quantified using the ferricyanide-based method described by Kidby and Davidson [19]. Measurements were performed in triplicate.

## 2.8. Influence of Temperature and pH on Xyn11-29 Activity and Enzyme Stability

To determine the influence of temperature and pH, the purified enzyme at 0.0541 mM was used. To determine the optimal pH, Xyn11-29 activity was quantified in the presence of beechwood xylans in a 50 mM sodium tartrate buffer from pH 2–4 and citrate-phosphate buffer from pH 5 to 7 at 30 °C. For the temperature profiles, the recombinant xylanase was incubated over the temperature range 20–70 °C in the presence of beechwood xylans in 50 mM citrate-phosphate buffer at pH 5.5.

To analyze thermal stability, Xyn11-29 was incubated at different temperatures (35, 40 and 50 °C) for 48 h in 50 mM citrate-phosphate buffer at pH 5.5, with regular sampling. For pH stability, the recombinant enzyme was incubated in 50 mM citrate-phosphate buffer at pH 5; 5.5 and 6 at 35 °C for 48 h. After incubation, activity was measured in standard conditions in the presence of beechwood xylans.

## 2.9. Influence of Sea Salt on Xylanase Activity

To determine the effect of sea salt on enzyme activity, activity was measured in standard conditions after adding sea salt (1–5% wt/vol) in 50 mM citrate-phosphate buffer at pH 5.5. To compare the effect of sea salt, a similar experiment was also performed with the purified GH11 xylanase Tx-Xyn11 from the bacterium *Thermobacillus xylanilyticus* in 50 mM sodium acetate buffer at 60 °C and at pH 5.8 [20].
## 2.10. Kinetic Properties

Kinetic parameters were determined in standard conditions with beechwood as substrate at concentrations of 1, 2, 5, 8, and 10 g/L. Experiments were performed in triplicate. Lineweaver-Burk plots were used to calculate the Michaelis constant $K_M$ and the enzyme turnover value $k_{cat}$. The normalized equation $v = (k_{cat}/K_M)[S]/(1 + [S]/K_M)$ yielded the catalytic efficiency values ($k_{cat}/K_M$) with their corresponding standard errors.

## 2.11. Lignocellulosic Biomass Hydrolysis

For the study of xylanase hydrolysis capacity, destarched wheat bran and wheat straw (Apache variety, 2 mm) provided by Agro Industrie Recherches Developpement (ARD, Pomacle, France) were used. Lignocellulose substrates (500 mg) were hydrated for 1 h in 10 mL of 50 mM citrate-phosphate buffer at pH 5.5 at 35 °C and in 10 mL of 50 mM sodium acetate buffer pH 5.8 at 60 °C for hydrolysis with the Xyn11-29 and Tx-Xyn11 experiment, respectively. Xyn11-29 and Tx-Xyn11 were then added at 100 IU per gram of substrate. Experiments were performed in duplicate under constant magnetic stirring (350 rpm) for 48 h. Hydrolysate samplings were performed after 0, 6, 24, and 48 h of incubation and the reaction was stopped by boiling at 100 °C for 10 min. Centrifugation (3000 × g, 10 min) separated supernatant and residual biomass. The activity of Tx-Xyn11 was also determined to compare the hydrolytic efficiency of the xylanase Xyn11-29. Tx-Xyn11 was selected for its high hydrolytic potential for lignocellulosic biomass [21].

The identification and quantification of the sugars released were carried out by HPAEC-PAD using a CarboPac PA-1 anion exchange column (4 × 250 mm, Dionex, Thermo Fisher Scientific, Waltham, MA, USA). All samples were hydrolyzed with 1M H₂SO₄ for 1 h at 100 °C. L-Fucose was used as internal standard and was added to the samples at a final concentration of 100 µg mL⁻¹. Each sample was then filtered (0.22 µm) before injection into the column. L-Arabinose and D-xylose were used as standards.

# 3. Results

## 3.1. Diversity of Fungal Expressed Genes Encoding GH11 Xylanases

The diversity of expressed genes encoding fungal GH11 xylanases was investigated in surface and deeper mangrove sediments beneath *A. marina* (A) and *R. stylosa* (R) trees during the wet (March) and dry (November) seasons in three independent plots. This diversity analysis was performed using a classical functional metabarcoding approach based on cDNAs previously enriched in fungal GH11 sequences. The normalized GH11 xylanase dataset consisted of 247,464 sequences distributed among only 15 different operational functional units (OFUs) corresponding to 15 putative GH11 xylanase encoding cDNAs (Figure 1, Supplementary Table S2). The number of OFUs per sediment sample ranged from one to seven 7. GH11 diversity, estimated with the Shannon index, was systematically higher in deep layers during the wet season (Figure 1A, Supplementary Table S3). The highest GH11 diversity was observed for the deep layers beneath the *R. stylosa* area during the wet season. By contrast, the lowest GH11 diversity was observed in the deep layers beneath the *A. marina* area during the dry season (Figure 1A, Supplementary Table S3).

Two OFUs (OFU1 and OFU2) largely dominated the GH11 expressed genes, representing 52% and 33% of the total number of sequences, respectively. Regarding their distribution, the OFUs were detected in almost all sediment samples independently of the season. However, OFU2 seemed more strongly expressed during the wet season (March) than the dry one (November) and not in all sediment samples (Figure 1B).

The difference in the composition of expressed GH11 xylanase encoding genes assessed through a non-parametric multivariate analysis of variance (PERMANOVA) highlighted strong effects of the season ($p = 0.01, R^2 = 0.158$) and of the tree species ($p = 0.02, R^2 = 0.104$). The sediment depth had a lower impact ($p = 0.20, R^2 = 0.039$) on the composition of the expressed fungal genes (Table 1).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Diversity of expressed genes encoding xylanases in surface and deeper sediments collected beneath two tree species (*Avicennia marina* and *Rhizophora stylosa*) during the wet (March) and dry (November) seasons. (A) Diversity estimated with the Shannon index and (B) distribution of the different operational functional units (OFUs).</figcaption></figure>

<table><caption>**Table 1.** Effect of environmental factors on the composition of expressed fungal genes encoding GH 11 xylanases.</caption><thead><tr><th></th><th>Df</th><th>F</th><th>p</th><th>R<sup>2</sup></th></tr></thead><tbody><tr><td><strong>Variable</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Tree species</td><td>1</td><td>4.369</td><td>0.02</td><td>0.104</td></tr><tr><td>Sediment depth</td><td>1</td><td>1.629</td><td>0.20</td><td>0.039</td></tr><tr><td>Season</td><td>1</td><td>6.641</td><td>0.01</td><td>0.158</td></tr><tr><td><strong>Interaction</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Tree × Depth</td><td>1</td><td>5.376</td><td>0.01</td><td>0.128</td></tr><tr><td>Tree × Season</td><td>1</td><td>4.755</td><td>0.02</td><td>0.113</td></tr><tr><td>Depth × Season</td><td>1</td><td>1.674</td><td>0.18</td><td>0.040</td></tr><tr><td>Tree × Depth × Season</td><td>1</td><td>1.579</td><td>0.19</td><td>0.038</td></tr><tr><td><strong>Residuals</strong></td><td>16</td><td></td><td></td><td></td></tr></tbody></table>

The differences between groups were tested using PERMANOVA analysis on Bray–Curtis dissimilarity matrices. Abbreviations: Df, degrees of freedom; F, F-test statistic.

## 3.2. Phylogeny of the Xylanase Sequences

The selection criteria for the enzyme target were: (i) representativity of the OFU in the mangrove sediments; (ii) obtainment of the full-length coding sequence by PCR; and (iii) production of the recombinant xylanase in *P. pastoris*. From the most abundant OFUs, OFU 2 and OFU 5 were amplified in our standard protocol, but not OFU1. The two amplified OFUs were not produced in *P. pastoris*, so we decided to study OFU29, a sequence close to OFU1 (99% identity 80/81 amino acid identity) that could be produced in the heterologous host strain. A BlastP was conducted at NCBI (non-redundant protein sequences) using Xyn11-29 sequence as a query. Xylanase sequences of the family GH11 from the orders Agaricales and Sebacinales emerged as proteins with the highest amino acid sequence identity (here, give the range of % identity for the first 10 best hits) (Table 2). Xyn11-29 underwent a phylogenetic analysis with the 50 best hit sequences obtained from the BlastP analysis. However, the origin of the sequence could not be inferred by this approach (not shown), probably owing to strong conservation of xylanase amino acid sequences among Agaricomycetes. In the current state of the available data, we can only suggest that Xyn11-29 has a fungal origin and belongs to a basidiomycete of the Agaricomycetes class. Alignment with the biochemically characterized mangrove-derived
xylanase XynMF13A confirmed that Xyn11-29 contained the two highly conserved catalytic residues characteristic of GH11 xylanases (Glu120 and Glu217) (Figure 2) and an analysis of the domain architecture showed that the 231 amino acid sequence was composed of a catalytic domain without a carbohydrate binding module.

<table><caption>**Table 2.** Putative xylanases presenting the highest amino acid identities with Xyn11-29. The sequences (GenBank accession number) and the fungal species are listed with their corresponding order in brackets, together with the percentage of identity and the number of residue pairs considered for each comparison.</caption><thead><tr><th></th><th>Xyn11-29</th><th></th></tr></thead><tbody><tr><td>KAF9530414.1 Crepidotus variabilis (Agaricales)</td><td>195/231(84%)</td><td></td></tr><tr><td>Lepista nuda (Agaricales)</td><td>193/231(84%)</td><td></td></tr><tr><td>PVF97526.1 Serendipita vermifera 'subsp. bescii' (Sebacinales)</td><td>195/235(83%)</td><td></td></tr><tr><td>ESK83196.1 Moniliophthora roreri (Agaricales)</td><td>187/231(81%)</td><td></td></tr><tr><td>KAF9530413.1 Crepidotus variabilis (Agaricales)</td><td>189/231(82%)</td><td></td></tr><tr><td>PVF91255.1 Serendipita vermifera 'subsp. bescii' (Sebacinales)</td><td>188/233(81%)</td><td></td></tr><tr><td>KAF5345258.1 Leucoagaricus leucothites (Agaricales)</td><td>193/232(83%)</td><td></td></tr><tr><td>CDS82539.1 uncultured eukaryote</td><td>188/233(81%)</td><td></td></tr><tr><td>KTB33854.1 Moniliophthora roreri (Agaricales)</td><td>182/231(79%)</td><td></td></tr><tr><td>KIM21487.1 Sebacina vermifera (Sebacinales)</td><td>182/231(79%)</td><td></td></tr></tbody></table><table><tbody><tr><td>MVSFTALVLTAVTAAGALAMPTN------ATEGVGLLQARGGTPSSTGTNNGYYYSWWTD</td><td>54</td></tr><tr><td>MVSFTSIITAAVAATGALAAPVADSSTDLAARSLGALTARAGTPGSQGTHNGCFYSWWTD</td><td>60</td></tr><tr><td>********: :**:*:****** *. *:..:*: * **.**.**.** * **:**** :********</td><td></td></tr><tr><td>GAAQATYTNNGGGSYSLQWSGNNGNLVGGKGWNPGYSGRVINYSGTYQPNGNSYLSVYGW</td><td>114</td></tr><tr><td>GGAQATYTNNGGGSYSVNWKT-GGNIVGGKGWNPGAA-RTTTYSGTYSPQGNSYLAIYGW</td><td>118</td></tr><tr><td>*********.********:*. **:************ : * . .*********.********:****</td><td></td></tr><tr><td>TRSPLIEYYIVESYGSYNPSSAAQRKGSVSCNGASYDILQTTRYNQPSIDGTQTFFQFWS</td><td>174</td></tr><tr><td>TRNPLVEYYVVENFGTYDPSSQASVKGSVTADGSSYKVAQTQRVNQPSIDGTQTFNQYWS</td><td>178</td></tr><tr><td>**.**:**:**:**.**:.*:**:** * . .**:**.**:**.** * * .**************:**:**</td><td></td></tr><tr><td>VRNPKKQPGGSISGTVDTACHFNAWRNYGMNLGSSHNYQIVATEGYQSSGSATINVS--</td><td>231</td></tr><tr><td>VRQNKRS-----SGSVNMKTHFDAWASKGMKLGTH-NYQIVATEGYFSSGSASITVNCA</td><td>231</td></tr><tr><td>**: *:. **:**: **:**. **:**:**: ************** .**:**.**</td><td></td></tr></tbody></table>

**Figure 2.** Sequence alignment of Xyn11-29 with the biochemically characterized xylanase XynMF13A isolated the mangrove fungus *Phoma* sp. (GenBank accession number AVV62005.1). The highly conserved glutamate catalytic residues are marked with blue arrows. Alignment was produced with Clustal Omega, and symbols below the sequences indicate full conservation of the same (asterisk) or equivalent residues (colon) and partial residue conservation (dot).

### 3.3. Heterologous Production and Purification of the Recombinant Xyn11-29

After the transformation of *P. pastoris* with the recombinant vector pPICZα-A containing the Xyn11-29 encoding cDNA, 50 transformants were selected for their resistance to zeocin and were then screened for the presence of Xyn11-29 in the extracellular medium. The best transformant was selected based on the band intensity corresponding to the recombinant protein and visualized following SDS-PAGE. The expected molecular weight for Xyn11-29 was 22.8 kDa, and after purification, the protein ran on SDS-PAGE at approximately 20 kDa with bands of slightly higher molecular masses probably corresponding to hyperglycosylated isoforms as already observed (Supplementary Figure S1) [6]. Xyn11-29 was predicted to possess three potential N-glycosylations at positions 4, 77, and 113 (starting from the first amino acids of the mature protein M PTNATEGVGL), as predicted via the NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/, accessed on 5 May 2021). The recombinant protein was purified by affinity chromatography using an IMAC column, yielding a protein concentration of 1.33 mg mL⁻¹, with a total activity of 1170 IU mL⁻¹ (Figure S1). The specific activity of the Xyn11-29 is 877 IU mg⁻¹.
### 3.4. Effect of pH and Temperature on the Activity and Stability of Xyn11-29

Recombinant Xyn11-29 was active under acidic pH conditions with an optimum at pH 5.5 for beechwood xylan hydrolysis in the pH range tested (Figure 3A). At lower pH, the Xyn11-29 activity markedly decreased, especially between pH 2 and 4, and almost no activity was detected at pH 2. At a higher pH, activities significantly decreased to less than 50% at pH 7. In the tested temperature range, the activity was optimal between 40 °C and 50 °C, while the enzyme was less active at lower temperatures (30–70% of its activity found at 50 °C) (Figure 3B). We note that Xyn11-29 was not active at temperatures above 60 °C. In conclusion, the optimal activity of Xyn11-29 was obtained at 40–50 °C in a 50 mM citrate-phosphate buffer at pH 5.5.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Effect of pH and temperature on enzyme activity: (A) Optimal pH in 50 mM sodium tartrate buffer from pH 2 to 4 and citrate-phosphate buffer from pH 5 to 7 at 30 °C; (B) Optimal temperature for beechwood xylan hydrolysis (0.5%, *w/v*) in the range 20–70 °C at pH 5.5. Activity values were calculated as a percentage of maximum activity (set to 100%) at optimum temperature and pH. Each data point (mean +/- standard deviation) is the result of triplicate experiments.</figcaption></figure>

Secondly, pH stability of the enzyme was determined. Activities were estimated after incubation at various pH values (5, 5.5 and 6) (Figure 4A). Results show that Xyn11-29 activity was relatively stable at pH 6, exhibiting 80% of its maximal activity after 48 h incubation. By contrast, at pH 5.5, the recombinant xylanase retained only 64% of its initial activity after 48 h of incubation, and at pH 5, 60% of its activity was lost after only 24 h of incubation (Figure 4A). Xyn11-29 displayed a relatively high thermal stability at 35 °C, with a half-life longer than 48 h (Figure 4B). Half-lives decreased with temperature, 1 h at 40 °C and only 3 min at 50 °C (data not presented on the graph for 50 °C). Taking into account the temperature and pH stabilities of Xyn11-29, standard conditions (40 °C, pH 5.5) were used for determining kinetic properties, and lignocellulosic biomass hydrolyses were performed at 35 °C and pH 5.5 for 48 h.

### 3.5. Catalytic Properties

The enzyme exhibited activity against beechwood xylan, with a $K_M$ value of 3.3 g L⁻¹ and a catalytic efficiency ($k_{cat}/K_M$) of 206 g L⁻¹ s⁻¹ (Table 3). Values measured for the fungal mangrove-derived xylanase, XynMF13A [12], were close to that found for Xyn11-29, with $K_M$ and $k_{cat}/K_M$ values of 3.16 g L⁻¹ and 340 g L⁻¹ s⁻¹, respectively, using beechwood xylan as the substrate.

<table><caption>**Table 3.** Kinetic constants of the recombinant Xyn11-29 produced in *P. pastoris*. Experiments were performed in triplicate against beechwood at various concentrations in the range 1–10 g L⁻¹. Activities were measured at 40 °C pH 5.5.</caption><thead><tr><th></th><th>Specific Activity<br>(IU mg<sup>-1</sup>)</th><th><i>K</i><sub>M</sub> (g L<sup>-1</sup>)</th><th><i>k</i><sub>cat</sub> (s<sup>-1</sup>)</th><th><i>k</i><sub>cat</sub>/<i>K</i><sub>M</sub><br>(L g<sup>-1</sup> s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>Xyn11-29</td><td>877 ± 37</td><td>3.30 ± 0.43</td><td>708 ± 52</td><td>216 ± 17</td></tr><tr><td>XynMF13A</td><td>1322 ± 5</td><td>3.16 ± 0.43</td><td>1075 ± 0.75</td><td>340</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** pH and temperature stability. (A) pH stability in 50 mM citrate-phosphate buffer from pH 5 to 6 at 35 °C; (B) Temperature stability from 35 °C to 50 °C at pH 5.5. Activity values were calculated as a percentage of maximum activity (set to 100%) at optimum temperature and pH. Each data point (mean +/- standard deviation) is the result of triplicate experiments.</figcaption></figure>

### 3.6. Xylanase Activity against Lignocellulosic Biomass

Xyn11-29 xylanase activity was compared with the xylanase from *Thermobacillus xylanilyticus* (*Tx*-Xyn11), already described for its efficiency for wheat straw xylan fractionation [21]. Xyn11-29 exhibited a high hydrolysis capacity on destarched wheat bran, with 40% and 16% of xylose and arabinose released after 24 h of hydrolysis, respectively (Figure 5A,B). These values are very similar to those obtained with *Tx*-Xyn11 xylanase. Conversely, the efficiency of the Xyn11-29 xylanase was low on wheat straw, with 2.6% and 6.9% of xylose and arabinose released after 24 h, respectively, while the efficiency of the *Tx*-Xyn11 xylanase was in general higher (Figure 5C,D).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Enzymatic hydrolysis capacity on destarched wheat bran (A,B) and wheat straw (C,D) of Xyn11-29 xylanase (red bars) compared with *T. xylanolyticus* xylanase, *Tx*-Xyn11 (gray bars) after 0, 6, 24, and 48 h incubation at pH 5.5 at 35 °C and at pH 5.8 at 60°C for Xyn11-29 and *Tx*-Xyn11, respectively. Results are expressed as a percentage of released sugars (arabinose and xylose). Each data point (mean +/- standard deviation) is the result of duplicate experiments.</figcaption></figure>
### 3.7. Influence of Sea Salt on Xylanase Activity

The influence of sea salt on xylanase activity was analyzed by activity measurement in standard conditions (40 °C at pH 5.5) with 1% and 5% of salt on the buffer used (Figure 6A). In the tested conditions, sea salt had no significant effect on Xyn11-29 xylanase activity even at 5% (wt/vol). In comparison, Tx-Xyn11 showed a slight decrease in its activity, particularly at 5% sea salt, with a reduction of 20% of its initial activity. In addition, analysis of the enzyme primary structure showed that Xyn11-29 had a very low recurrence of negative surface charges, with (D + E)/(R + K) = 0.73, similar to what was calculated for XynMF13a (Figure 6B). The terrestrial Tx-Xyn11 had an even balance of negative and positive charges.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Sea salt effect on Xyn11-29 xylanase activity and ratios of negatively to positively charged amino acids for Xyn11-29, *T. xylanolyticus*, *Tx-Xyn11*, and *Phoma* sp. xylanase, XynMF13A. (A) Sea salt was used at concentrations of 1% and 5%. Activities are expressed as a percentage of the activity without sea salt (set to 100%). Each data point (mean +/- standard deviation) is the result of triplicate experiments. (B) Negatively charged (D + E) to positively charged (R + K) amino acid ratios.</figcaption></figure>

## 4. Discussion

The functional diversity of fungal xylanases was investigated during both wet and dry seasons in New Caledonian mangrove sediments, from the surface to deeper horizons under the two most common tree species in this region (*Avicennia marina* and *Rhizophora stylosa*). Xylanases are involved in breaking down hemicellulose, one of the major components of plant cell walls, and they are largely characterized in the fungal kingdom. They were accordingly selected as markers for the biomass enzymatic degradation in the mangrove soil [7]. Compared with soils from terrestrial forests, which can harbor up to 103 GH11 OFUs [14], with only 15 GH11 OFUs, mangrove sediments showed a much lower diversity. This lower diversity observed may be the result of physical and chemical gradients associated with mangrove sediments that are inimical to fungal colonization and related enzymatic activities (e.g., low oxygen and nitrogen contents) [22]. The diversity of fungal communities was investigated in a previous study and as observed here for the expressed GH11 genes, the fungal diversity was also strongly affected by season and tree species [5]. This impact of tree species on fungal communities and their expressed genes is consistent with previous work conducted in classical forest ecosystems demonstrating that tree species select diverse soil fungal communities expressing different sets of lignocellulolytic enzyme-encoding genes [14]. Interestingly, the diversity of expressed genes encoding GH11 xylanases was systematically higher in deep layers during the wet season. Such higher expressed gene diversity may result from the combination of optimal condition for fungal xylanase activities such as higher carbon content in deeper layers [22] and a wet season favorable to fungal colonization in mangrove ecosystems [23].

Xyn11-29 (OFU29) was selected from the microbial diversity of the New Caledonia sediments as it was the closest homolog to the most abundant OFU (OFU1) with 99% identity, and OFU1 could not be produced in *P. pastoris*. In the current state of the data, the origin of the protein could not be inferred by a phylogenetic approach. However,
we make the supposition that the protein is a xylanase from the family GH11, which we
call Xyn11-29. In addition, Xyn11-29 shares high similarity with GH11 enzymes from the
Agaricomycetes class, within the Agaricomycotina sub- phylum. In a previous study, the
same mangrove sediments showed a predominance of Ascomycota (76.8–94.5% of total
ITS) and Basidiomycota (5.3–18.3% of total ITS) [3]. In addition, the most abundant OFU
identified in the same mangrove sediments using DyP as an alternative enzymatic marker
was attributed to a fungus belonging to the Agaricomycetes class, in agreement with our
present results [6].

A characterization of the mangrove-derived Xyn11-29 was conducted after its pro-
duction in *P. pastoris* and its purification. Concerning its general properties, Xyn11-29 was
shown to possess only a catalytic module without CBM, as is the case for 80% of xylanases
(family GH11) from eukaryotic organisms [24]. Xyn11-29 has the average size for a xy-
lanase of around 20 kDa and possesses the two key catalytic residues characteristic of the
xylanases (Glu120 and Glu217) [25]. Kinetic parameters were determined and compared
with the biochemically characterized XynMF13A isolated from the Shankou mangrove
in Guangxi province, China [12]. These two mangrove-derived xylanases exhibited very
close kinetic parameters, with a K<sub>M</sub> of 3 g L<sup>-1</sup> against beechwood xylans and a catalytic
efficiency (k<sub>cat</sub>/K<sub>M</sub>) around 200–300 g L<sup>-1</sup> s<sup>-1</sup>. No other kinetic values could be found
in the literature concerning other marine-derived xylanases. Compared with terrestrial
xylanases, K<sub>M</sub> values were relatively close, with values of 3.42–9.96 g L<sup>-1</sup> for xylanases
from *Aspergillus oryzae* and *Phanerochaete chrysosporium* (for a review see [26]). Concern-
ing stability and activity as a function of temperature and pH, Xyn11-29 exhibited rather
generic properties, with an optimal activity at pH 5–6, and a maximum activity between
40 °C and 50 °C. In addition, the recombinant protein was stable at pH 5 and 6, with a
low stability at high temperature, especially at 50 °C, with a half-life of 3 min. However,
Xyn11-29 possesses high thermal stability at lower temperature and so can be used at
moderate temperature. The mangrove-derived XynMF13A has similar properties [12],
and among the multiple examples of fungal xylanases found in terrestrial fungi the same
properties could be mostly found [27–30]. Considering these properties, we can conclude
that Xyn11-29 shares the general characteristics of the xylanases found in terrestrial fungi
and with the sole characterized purified xylanase isolated from a marine area, namely
mangrove sediments.

As xylanases have already been gainfully used for several years and, more recently,
in lignocellulose degradation. In the context of plant biomass refining [31], we tested and
compared Xyn11-29 with xylanase Tx-Xyn11 from *T. xylanilyticus*, known to efficiently
hydrolyze lignocellulosic biomass [21,32]. Xyn11-29 exhibited similar efficiency to Tx-
xyn11 for the release of sugars (xylose and arabinose) from destarched wheat bran. The
release was high, reaching 17% and 42%, respectively, for arabinose and xylose after 48 h.
Lower xylose and arabinose release was obtained on wheat straw for Xyn11-29 with 2.8%
and 6.9%, respectively. The efficiency of Tx-xyn11 on wheat straw was higher. The high
hydrolytic efficiency of Xyn11-29 on destarched wheat bran compared with wheat straw
could be explained by the low recalcitrance of wheat bran relative to wheat straw, as
the lignin content in wheat bran is low (less than 5% dry matter). The lower ability of
Xyn11-29 to hydrolyze the xylans from wheat straw relative to the xylanase Tx-xyn11 could
be explained by the kinetic parameters of both enzymes on xylans with a low degree of
branching, as is the case for wheat straw xylans [33]. When measured on beechwood
xylans, the catalytic efficiency of Xyn11-29 reached 216 L g⁻¹ s⁻¹, whereas the k₅ₐt/Kₘ was
4684 L g⁻¹ s⁻¹ for Tx-xyn11. However, this ability to degrade lignocellulosic substrates
to some extent and at moderate temperature suggests biotechnological applications for
which Xyn11-29 could be usefully tested, such as in the pulp and paper industry and the
food sector.

Several studies have been undertaken on marine microorganisms in order to discover
new biocatalysts with activity resistant to saline conditions and, more generally, to harsh
conditions (as required in the industry). For instance, the dye-colorizing sector and the
food industry need more robust biocatalysts that could exhibit these properties, and marine
organisms that have adapted to salty environments are thought to produce and secrete
enzymes that are active, especially in saline conditions [34]. Xyn 11-29 is active at up to
5% of sea salt. Activities of the mangrove-derived XynMF13A were tested in the presence
of NaCl and were quite stable after 1 h incubation in 3 M NaCl [12]. In addition, this
xylanase retained 80% of its initial activity in up to 4 M NaCl. Research has been seeking
to determine why some of these marine enzymes and not others possess this interesting
property of acting in saline conditions. To our knowledge, only three other studies with
the objective of screening marine-derived fungal strains for xylanase activities have been
published [35–37]. In these studies, only the crude extracellular medium of the fungal
culture was characterized and not the purified xylanase, so no data were provided on
the effect of salt on the related enzymes. However, other studies are available on the
effect of the saline conditions on the activity of other enzymes. For instance, the dye-
decolorizing peroxidase DyP1, purified from a fungal strain of the same New Caledonia
mangrove, was negatively affected by sea salt, with a complete inhibition of its activity at
3% of sea salt. There was a similar effect on the terrestrial *T. versicolor* DyP1, *TvDyP1* [6].
Likewise, but with a less pronounced effect, a laccase purified from the marine-derived
*Stemphylium lucomagnoense* strain from the Tunisian coast lost 50% of its initial activity
in 4% of sea salt [38]. By contrast, two laccases obtained from *Pestalotiopsis* sp. KF079
isolated from the Baltic sea mudflats were positively activated to 200–360% of their initial
activities [34]. Some studies showed that salt-adapted enzymes were characterized by
highly negative surface charges thought to contribute to protein stability and activity in
extreme osmotic conditions [39–41]. In addition, solvation has been suggested to be a key
requirement for maintaining solubility and activity of enzymes in low water activity [42]. In
this context, hydrogen bonds between negatively charged side chains and water molecules
become critical to maintaining a stable hydration shell surrounding the enzyme. In line
with this hypothesis, the laccase *PsLac2* of *Pestalotiopsis* sp. KF079 was shown to possess
a high ratio of negative to positive amino acids (D + E)/(R + K) of 3.95. Other enzymes
with an activity that remained stable or decreased in the presence of sea salts, including
Xyn11-29, have a appreciably lower ratio with values in the range from 0.67 (mangrove-
derived XynMF13A) to 1.0 (terrestrial-derived *Tx-Xyn11*). Also, the mangrove-derived
DyP1, which was negatively affected by sea salt, presents a balanced ratio of 1.19, similar
to that of *Tx-Xyn11*. Accordingly, and as was concluded for mangrove-derived DyP1, we
hypothesized that these two enzymes were proteins reminiscent of those of terrestrial origin
in their behavior in saline conditions. A recent study has shown that leaves and aerial
bodies predominantly contain fungi derived from well-studied terrestrial environments
with functions similar to their terrestrial plant counterparts [43]. By contrast, sediments
and pneumatophores are colonized by marine fungi, especially since they are periodically
submerged by tides, requiring them to grow in an environment that experiences large
and rapid fluctuations in salinity. Further research will be needed to extend our work
on the mangrove by cloning other functional enzymes and studying their origin to gain
more general insight into this complex marine environment. Through a broader inventory
of functional enzymes, a systematic characterization of model enzymes such as DyP
and xylanase should yield interesting data on the marine versus terrestrial origin of the
producing strains.

**5. Conclusions**

Mangroves are original, environmentally essential ecosystem located at the intertidal
zone between land and sea. They are adapted to a gradient of salinity, temperature,
humidity, and tree species, resulting in a diversity of the microbiota and enzymes they
produce. In this context, we have developed an approach to study the functional diversity
of the New Caledonia mangrove, which could find use in biotechnological applications.
The selected enzyme model, Xyn11-29, was shown to possess the general characteristics
of xylanases from family GH11, and to recall terrestrial xylanases with biotechnological abilities like those of their terrestrial homologs.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9071484/s1: Table S1, List of tagged primers used to amplify expressed genes encoding fungal endo-β-1,4-xylanases (GH11); Table S2, Normalized dataset where the 247,464 fungal GH11 sequences are distributed among the 15 different Operational Functional Units (OFUs) within the sediment samples; Table S3, Mean number of fungal GH11 OFUs (Operational Function Units) detected and diversity indices calculated for oxic and anoxic sediment fractions collected during the wet (March) and dry (November) seasons within two mangrove pristine areas (Avicennia marina and Rhizophora stylosa); Figure S1, SDS-PAGE of the purified Xyn11-29 produced in *P. pastoris*.

**Author Contributions:** Conceptualization, M.H. (Mylène Hugoni), C.M., P.L., and E.R.; methodology, C.I., L.V., M.D., A.B., M.H. (Mireille Haon), S.G., A.J., M.H. (Mylène Hugoni), M.-N.R., P.L., and E.R.; software, M.-N.R., P.L.; validation, M.H. (Mylène Hugoni), M.-N.R., C.R., P.L., and E.R.; formal analysis, C.I., S.G., E.B., C.B.F., G.S., P.L., and E.R.; investigation, C.I., L.V., M.D., A.B., M.H. (Mireille Haon), S.G., A.J., M.H. (Mylène Hugoni), P.L., and E.R.; resources, P.L.; data curation, C.I., P.L., and E.R.; writing—original draft preparation, C.I., P.L., and E.R.; writing—review and editing, C.I., S.G., E.B., C.B.F., G.S., C.M., M.H. (Mylène Hugoni), H.R., C.R., P.L., and E.R.; supervision, H.R., C.R., P.L., and E.R.; project administration, P.L. and E.R.; funding acquisition, P.L. and E.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Ecologie Côtière et Continentale (EC2CO, 2017–2019).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The sequence data generated in this study were deposited in the EMBL-ENA public database (PRJEB43346 for the first sampling campaign (March) dataset and PRJEB43343 for the second sampling campaign (November) dataset).

**Conflicts of Interest:** The authors have no conflicts of interest to declare. The funders had no role in the design of the study, the collected analyses, the interpretation of data, the writing of the manuscript, or the decision to publish the results.

## References

1. Kristensen, E.; Connolly, R.M.; Otero, X.L.; Marchand, C.; Ferreira, T.O.; Rivera-Monroy, V.H. Biogeochemical Cycles: Global Approaches and Perspectives. *Mangrove Ecosyst. A Glob. Biogeogr. Perspect. Struct. Funct. Serv.* **2017**, *163–209*. [CrossRef]

2. Hyde, K.D.; Lee, S.Y. Ecology of Mangrove Fungi and Their Role in Nutrient Cycling: What Gaps Occur in Our Knowledge? *Hydrobiologia* **1995**, *295*, 107–118. [CrossRef]

3. Bourgeois, C.; Alfaro, A.C.; Bisson, E.; Alcius, S.; Marchand, C. Trace Metal Dynamics in Soils and Plants along Intertidal Gradients in Semi-Arid Mangroves (New Caledonia). *Mar. Pollut. Bull.* **2020**, *156*, 111274. [CrossRef] [PubMed]

4. Arfi, Y.; Buée, M.; Marchand, C.; Levasseur, A.; Record, E. Multiple Markers Pyrosequencing Reveals Highly Diverse and Host-Specific Fungal Communities on the Mangrove Trees *Avicennia marina* and *Rhizophora stylosa*. *FEMS Microbiol. Ecol.* **2012**, *79*, 433–444. [CrossRef]

5. Luis, P.; Saint-Genis, G.; Vallon, L.; Bourgeois, C.; Bruto, M.; Marchand, C.; Record, E.; Hugoni, M. Contrasted Ecological Niches Shape Fungal and Prokaryotic Community Structure in Mangroves Sediments. *Environ. Microbiol.* **2019**, *21*, 1407–1424. [CrossRef] [PubMed]

6. Ayed, A.B.; Saint-Genis, G.; Vallon, L.; Linde, D.; Turbé-Doan, A.; Haon, M.; Daou, M.; Bertrand, E.; Faulds, C.B.; Sciara, G.; et al. Exploring the Diversity of Fungal DyPs in Mangrove Soils to Produce and Characterize Novel Biocatalysts. *J. Fungi* **2021**, *7*, 321. [CrossRef]

7. Linde, D.; Ayuso-Fernández, I.; Laloux, M.; Aguiar-Cervera, J.E.; de Lacey, A.L.; Ruiz-Dueñas, F.J.; Martínez, A.T. Comparing Ligninolytic Capabilities of Bacterial and Fungal Dye-Decolorizing Peroxidases and Class-II Peroxidase-Catalases. *Int. J. Mol. Sci.* **2021**, *22*, 2629. [CrossRef]

8. Arfi, Y.; Chevret, D.; Henrissat, B.; Berrin, J.-G.; Levasseur, A.; Record, E. Characterization of Salt-Adapted Secreted Lignocellulolytic Enzymes from the Mangrove Fungus *Pestalotiopsis* sp. *Nat. Commun.* **2013**, *4*, 1810. [CrossRef] [PubMed]

9. Singh, B. Production, Characteristics, and Biotechnological Applications of Microbial Xylanases. *Appl. Microbiol. Biotechnol.* **2019**, *103*, 8763–8784. [CrossRef]

10. Henrissat, B.; Davies, G. Structural and Sequence-Based Classification of Glycoside Hydrolases. *Curr. Opin. Struct. Biol.* **1997**, *7*, 637–644. [CrossRef]
11. Qeshmi, F.I.; Homaei, A.; Fernandes, P.; Hemmati, R.; Dijkstra, B.W.; Khajeh, K. Xylanases from Marine Microorganisms: A Brief Overview on Scope, Sources, Features and Potential Applications. *Biochim. Biophys. Acta (BBA) Proteins Proteom.* **2020**, *1868*, 140312. [CrossRef]

12. Wu, J.; Qiu, C.; Ren, Y.; Yan, R.; Ye, X.; Wang, G. Novel Salt-Tolerant Xylanase from a Mangrove-Isolated Fungus *Phoma* sp. MF13 and Its Application in Chinese Steamed Bread. *ACS Omega* **2018**, *3*, 3708–3716. [CrossRef] [PubMed]

13. Bragalini, C.; Ribière, C.; Parisot, N.; Vallon, L.; Prudent, E.; Peyretaillade, E.; Girlanda, M.; Peyret, P.; Marmeisse, R.; Luis, P. Solution Hybrid Selection Capture for the Recovery of Functional Full-Length Eukaryotic CDNs from Complex Environmental Samples. *DNA Res.* **2014**, *21*, 685–694. [CrossRef]

14. Barbi, F.; Bragalini, C.; Vallon, L.; Prudent, E.; Dubost, A.; Fraissinet-Tachet, L.; Marmeisse, R.; Luis, P. PCR Primers to Study the Diversity of Expressed Fungal Genes Encoding Lignocellulolytic Enzymes in Soils Using High-Throughput Sequencing. *PLoS ONE* **2014**, *9*, e116264. [CrossRef]

15. Anderson, M.J. Permutational Multivariate Analysis of Variance (PERMANOVA). In *Wiley StatsRef: Statistics Reference Online*; American Cancer Society: Atlanta, GA, USA, 2017; pp. 1–15. ISBN 978-1-118-44511-2.

16. Jones, D.T.; Taylor, W.R.; Thornton, J.M. The Rapid Generation of Mutation Data Matrices from Protein Sequences. *Bioinformatics* **1992**, *8*, 275–282. [CrossRef]

17. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef]

18. Haon, M.; Grisel, S.; Navarro, D.; Gruet, A.; Berrin, J.-G.; Bignon, C. Recombinant Protein Production Facility for Fungal Biomass-Degrading Enzymes Using the Yeast *Pichia pastoris*. *Front. Microbiol.* **2015**, *6*, 1002. [CrossRef] [PubMed]

19. Kidby, D.K.; Davidson, D.J. A Convenient Ferricyanide Estimation of Reducing Sugars in the Nanomole Range. *Anal. Biochem.* **1973**, *55*, 321–325. [CrossRef]

20. Debeire-Gosselin, M.; Loonis, M.; Samain, E.; Debeire, P. Purification and Properties of a 22 kDa Endoxylanase Excreted by a New Strain of Thermophilic *Bacillus*. *Prog. Biotechnol.* **1992**, *7*, 463–466.

21. Rémond, C.; Aubry, N.; Crônier, D.; Noël, S.; Martel, F.; Roge, B.; Rakotoarivonina, H.; Debeire, P.; Chabbert, B. Combination of Ammonia and Xylanase Pretreatments: Impact on Enzymatic Xylan and Cellulose Recovery from Wheat Straw. *Bioresour. Technol.* **2010**, *101*, 6712–6717. [CrossRef]

22. Aschenbroich, A.; Marchand, C.; Molnar, N.; Deborde, J.; Hubas, C.; Rybarczyk, H.; Meziane, T. Spatio-Temporal Variations in the Composition of Organic Matter in Surface Sediments of a Mangrove Receiving Shrimp Farm Effluents (New Caledonia). *Sci. Total. Environ.* **2015**, *512–513*, 296–307. [CrossRef] [PubMed]

23. Saravanakumar, K.; Rajendran, N.; Kathiresan, K.; Chen, J. Bioprospects of Microbial Enzymes from Mangrove-Associated Fungi and Bacteria. In *Advances in Food and Nutrition Research*; Academic Press: London, UK, 2016; Volume 79, pp. 99–115. [CrossRef]

24. Talens-Perales, D.; Sánchez-Torres, P.; Marín-Navarro, J.; Polaina, J. In Silico Screening and Experimental Analysis of Family GH11 Xylanases for Applications under Conditions of Alkaline PH and High Temperature. *Biotechnol. Biofuels* **2020**, *13*, 198. [CrossRef]

25. Paës, G.; Berrin, J.-G.; Beaugrand, J. GH11 Xylanases: Structure/Function/Properties Relationships and Applications. *Biotechnol. Adv.* **2012**, *30*, 564–592. [CrossRef]

26. Ahmed, S.; Riaz, S.; Jamil, A. Molecular Cloning of Fungal Xylanases: An Overview. *Appl. Microbiol. Biotechnol.* **2009**, *84*, 19–35. [CrossRef]

27. Li, C.; Kumar, A.; Luo, X.; Shi, H.; Liu, Z.; Wu, G. Highly Alkali-Stable and Cellulase-Free Xylanases from Fusarium Sp. 21 and Their Application in Clarification of Orange Juice. *Int. J. Biol. Macromol.* **2020**, *155*, 572–580. [CrossRef]

28. Pasin, T.M.; Salgado, J.C.S.; de Almeida Scarcella, A.S.; de Oliveira, T.B.; de Lucas, R.C.; Cereia, M.; Rosa, J.C.; Ward, R.J.; Buckeridge, M.S.; de Moraes, M.D.L.T. A Halotolerant Endo-1,4-β-Xylanase from *Aspergillus Clavatus* with Potential Application for Agroindustrial Residues Saccharification. *Appl. Biochem. Biotechnol.* **2020**, *191*, 1111–1126. [CrossRef] [PubMed]

29. García-Huante, Y.; Cayetano-Cruz, M.; Santiago-Hernández, A.; Cano-Ramírez, C.; Marsch-Moreno, R.; Campos, J.E.; Aguilar-Osorio, G.; Benitez-Cardoza, C.G.; Trejo-Estrada, S.; Hidalgo-Lara, M.E. The Thermophilic Biomass-Degrading Fungus *Thielavia terrestris* Co3Bag1 Produces a Hyperthermophilic and Thermostable β-1,4-Xylanase with Exo- and Endo-Activity. *Extremophiles* **2017**, *21*, 175–186. [CrossRef]

30. Heinen, P.R.; Bauermeister, A.; Ribeiro, L.F.; Messias, J.M.; Almeida, P.Z.; Moraes, L.A.B.; Vargas-Rechia, C.G.; De Oliveira, A.H.C.; Ward, R.J.; Kadowaki, M.K.; et al. GH11 Xylanase from *Aspergillus tamarii* Kita: Purification by One-Step Chromatography and Xylooligosaccharides Hydrolysis Monitored in Real-Time by Mass Spectrometry. *Int. J. Biol. Macromol.* **2018**, *108*, 291–299. [CrossRef] [PubMed]

31. Bhardwaj, N.; Verma, P. Microbial Xylanases: A Helping Module for the Enzyme Biorefinery Platform. In *Bioenergy Research: Revisiting Latest Development*; Manish, S., Neha, S., Rajeev, S., Eds.; Springer: Singapore, 2021.

32. Rakotoarivonina, H.; Revol, P.-V.; Aubry, N.; Rémond, C. The Use of Thermostable Bacterial Hemicellulases Improves the Conversion of Lignocellulosic Biomass to Valuable Molecules. *Appl. Microbiol. Biotechnol.* **2016**, *100*, 7577–7590. [CrossRef]

33. Rakotoarivonina, H.; Hermant, B.; Aubry, N.; Rabenoelina, F.; Baillieul, F.; Rémond, C. Dynamic Study of How the Bacterial Breakdown of Plant Cell Walls Allows the Reconstitution of Efficient Hemicellulasic Cocktails. *Bioresour. Technol.* **2014**, *170*, 331–341. [CrossRef]
34. Wikee, S.; Hatton, J.; Turbé-Doan, A.; Mathieu, Y.; Daou, M.; Lomascolo, A.; Kumar, A.; Lumyong, S.; Sciara, G.; Faulds, C.B.; et al. Characterization and Dye Decolorization Potential of Two Laccases from the Marine-Derived Fungus *Pestalotiopsis* sp. *Int. J. Mol. Sci.* **2019**, *20*, 1864. [CrossRef] [PubMed]

35. Dos Santos, J.A.; Vieira, J.M.F.; Videira, A.; Meirelles, L.A.; Rodrigues, A.; Taniwaki, M.H.; Sette, L.D. Marine-Derived Fungus *Aspergillus* Cf. *tubingensis* LAMAI 31: A New Genetic Resource for Xylanase Production. *AMB Express* **2016**, *6*, 25. [CrossRef] [PubMed]

36. Torres, J.M.O.; Dela Cruz, T.E.E. Production of Xylanases by Mangrove Fungi from the Philippines and Their Application in Enzymatic Pretreatment of Recycled Paper Pulps. *World J. Microbiol. Biotechnol.* **2013**, *29*, 645–655. [CrossRef]

37. Martinho, V.; Dos Santos Lima, L.M.; Barros, C.A.; Ferrari, V.B.; Passarini, M.R.Z.; Santos, L.A.; de Souza Sebastianes, F.L.; Lacava, P.T.; de Vasconcellos, S.P. Enzymatic Potential and Biosurfactant Production by Endophytic Fungi from Mangrove Forest in Southeastern Brazil. *AMB Express* **2019**, *9*, 130. [CrossRef] [PubMed]

38. Ali, W.B.; Ayed, A.B.; Turbé-Doan, A.; Bertrand, E.; Mathieu, Y.; Faulds, C.B.; Lomascolo, A.; Sciara, G.; Record, E.; Mechichi, T. Enzyme Properties of a Laccase Obtained from the Transcriptome of the Marine-Derived Fungus *Stemphylium lucomagnoense*. *Int. J. Mol. Sci.* **2020**, *21*, 8402. [CrossRef]

39. Kern, M.; McGeehan, J.E.; Streeter, S.D.; Martin, R.N.A.; Besser, K.; Elias, L.; Eborall, W.; Malyon, G.P.; Payne, C.M.; Himmel, M.E.; et al. Structural Characterization of a Unique Marine Animal Family 7 Cellobiohydrolase Suggests a Mechanism of Cellulase Salt Tolerance. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 10189–10194. [CrossRef] [PubMed]

40. Paul, S.; Bag, S.K.; Das, S.; Harvill, E.T.; Dutta, C. Molecular Signature of Hypersaline Adaptation: Insights from Genome and Proteome Composition of Halophilic Prokaryotes. *Genome Biol.* **2008**, *9*, R70. [CrossRef] [PubMed]

41. Lanyi, J.K. Salt-Dependent Properties of Proteins from Extremely Halophilic Bacteria. *Bacteriol. Rev.* **1974**, *38*, 272–290. [CrossRef]

42. DasSarma, S.; DasSarma, P. Halophiles and Their Enzymes: Negativity Put to Good Use. *Curr. Opin. Microbiol.* **2015**, *25*, 120–126. [CrossRef]

43. Lee, N.L.Y.; Huang, D.; Quek, Z.B.R.; Lee, J.N.; Wainwright, B.J. Mangrove-Associated Fungal Communities Are Differentiated by Geographic Location and Host Structure. *Front. Microbiol.* **2019**, *10*, 2456. [CrossRef]