Background {#Sec1}
==========

Cotton linters are an important by-product of the textile industry, being the short fibre fraction that cannot be used in the textile process \[[@CR1]\]. They are obtained from cotton plant (*Gossypium* sp.), an annual shrub harvested for their high industrial interest. Cotton linters consist of high-quality cellulose fibres presenting very high cellulose content (98%) \[[@CR2]\]. They are typically used in special applications such as the production of cellulose derivatives, regenerated cellulose, or the manufacture of high added-value papers \[[@CR3]\].

In order to construct new materials and products based on renewable resources, the interest in functionalizing cellulose has gained importance in the last years. In fact, there is an increasing interest in replacing synthetic polymers by more sustainable materials to replace petroleum-based products in the packaging sector \[[@CR4]\]. Modification of cellulose by chemical or biotechnological means has been reported \[[@CR5]\]. Hydrolytic enzymes like cellulases can successfully modify cellulose, improving its reactivity and also altering fibre morphology \[[@CR6]\]. On the other hand, the oxidative enzymatic system Laccase_Tempo (2,2,6,6-tetramethyl-1-piperidinyloxy) can create new functional groups to cellulose converting primary hydroxyl groups to aldehyde or carboxyl forms \[[@CR7], [@CR8]\]. In this system, laccase, having a redox potential in the range of 0.7--0.9 V, can easily oxidize the stable oxyl-radical form of Tempo to oxoammonium ion (*E*° 0.2 V). This ion is the actual oxidant of cellulose, which can be regenerated by laccase oxidation or by acid-induced disproportionation.

A new generation of enzymes that also create functional groups in cellulose and other crystalline polysaccharides such as chitin, lytic polysaccharide monooxygenases (LPMO), have been discovered \[[@CR9], [@CR10]\]. They oxidatively cleave glycosidic linkages, leading to the formation of oxidized glucose units at different positions, resulting in the formation of aldonic acids at the C1 position and/or 4-ketoaldoses (gemdiols) at the C4 position \[[@CR11]\]. This oxidation renders the substrate more susceptible to be hydrolysed by conventional cellulases and is considered as a breakthrough in the enzymatic degradation of cellulose \[[@CR12]\]. The enzymatic effects that LPMO produce in cellulose have been mainly evaluated through their increase in cellulose degradation \[[@CR13], [@CR14]\]. However, the effect that LPMO produces on pulp fibres has been poorly investigated \[[@CR15]--[@CR17]\]. These authors demonstrated that LPMO weakens fibres cohesion, promoting their disruption during mechanical fibrillation.

The production of nanocrystalline cellulose (NCC) from cotton linters has been reported, and also the ability of cellulases to improve its yield \[[@CR2], [@CR18]\]. However, little knowledge exists about the production of nanofibrillated cellulose (NFC) from these fibres \[[@CR19]--[@CR21]\]. Interestingly, cotton provides fibres with promising interest in nanocellulose production due to its high purity and highly crystalline cellulose \[[@CR2]\]. NFC is usually produced by high-pressure homogenization, being major impediments for its commercial success the very high energy consumption of the production process and the clogging of homogenizers. Therefore, some pretreatments are needed in order to facilitate this process \[[@CR22]\]. The ability of cellulases to improve this process has been demonstrated \[[@CR23]--[@CR26]\]. Also, the improvement of mechanical fibrillation produced with Tempo--NaBr--NaClO system is well known \[[@CR20]\]. In order to replace the halide-based co-oxidizer system, laccase can be used to oxidize Tempo.

In this work, several enzymes were applied on cotton linters in order to analyse their effects on sugar release and on cellulose and fibre modifications. Four hydrolytic enzymes (cellulases) and also two oxidative enzymatic systems (a new bacterial LPMO and the Laccase_Tempo system) were used for this purpose. The effects that enzymes produced during mechanical fibrillation were also analysed, together with the optical, physical and barrier properties of the films casted from these treated linters.

Materials and methods {#Sec2}
=====================

Raw material {#Sec3}
------------

Cotton linters obtained from the second cut were supplied by CELESA (Celulosa de Levante S.A.), Tortosa, Spain. The initial fibres had an average of 0.47 mm of length, 19.67 μm of width and 38.87% fine content. Their drainability, measured as ºSR, was 12. Prior to the enzymatic pretreatments, cotton linters were refined in a valley mill for 24 h in order to reduce their average length. Refined linters, named as "R", had an average fibre length of 0.25 mm, fibre with 25.5 μm and 52.58% of fines; their drainability was increased to 77ºSR.

A commercial NFC (Com) supplied by the University of Maine, with 90% of fines, was used for comparison.

Enzymes {#Sec4}
-------

Four hydrolytic enzymes (cellulases) and two oxidative enzymes (LPMO and Laccase_Tempo) were used as pretreatments in cotton linters. Cel9B from *Paenibacillus barcinonensis* BP-23 \[[@CR27]\] was a monocomponent processive endoglucanase named as "C~9~". A commercial cellulase from Sertec20 was named as "C~50~", whereas two commercial cellulases supplied by Novozymes^®^ (Fibercare and Celluclast) were named as "C~F~" and "C~ll~", respectively. Their initial enzymatic activities were 5.5, 383, 99 and 536 U mL^−1^ for C~9~, C~50~, C~F~ and C~ll~, respectively. The commercial cellulase preparations used were not monocomponent and contained mixtures of several enzymes. Enzymatic activity was assayed by measuring the amount of reducing sugars released from carboxymethylcellulose (CMC) by the dinitrosalicylic (DNS) reagent method \[[@CR28]\]. The standard assay (100 µL reaction volume) was performed at 50 °C in 50 mM potassium acetate buffer at pH 5 for 15 min. One unit of enzymatic activity (U) was defined as the amount of enzyme that releases 1 µmol of reducing sugar equivalent per min under the assay conditions described. A standard curve of glucose was used to calculate activity units. All determinations of enzyme activity were made in triplicate.

For the oxidative treatments, an LPMO from *Streptomyces ambofaciens* (SamLPMO10C) \[[@CR29]\] and a laccase from *Trametes villosa* in combination with Tempo (2,2,6,6-tetramethyl-1-piperidinyloxy) were used. They were named as "S" and "L_Tempo", respectively. Laccase was supplied by Novozymes^®^ (Denmark) and had an activity of 746 U mL^−1^. Tempo was purchased from Sigma-Aldrich. The laccase activity was measured as the extent of oxidation of 5 mM 2,20-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) to its cation radical (*ε*~436~ = 29,300 M^−1^ cm^−1^) in 0.1 M sodium acetate buffer (pH 5) at 24 °C. One activity unit (U) was defined as the amount of enzyme converting 1 µmol of ABTS per min.

Enzymatic pretreatments on cotton linters {#Sec5}
-----------------------------------------

Pretreatments with cellulases were performed with 5 g odp (oven-dried pulp) at 10% consistency, with 10 U g^−1^ odp of enzyme in 50 mM potassium acetate buffer, pH 5, at 50 °C for 18 h. A combined treatment with C~F~ and C~ll~, named as "C~mix~", was also performed. This pretreatment was performed as described above but with 10 U g^−1^ odp of C~F~ and 10 U g^−1^ odp of C~ll~. Treatment with LPMO (S treatment) was performed with 5 g odp and 4 mg of enzyme g^−1^ odp at 5% consistency, for 72 h at 50 °C in 10 mM of ammonium acetate buffer at pH 6, with 2 mM ascorbic acid and 20 μM hydrogen peroxide. L_Tempo oxidation treatments were performed at room temperature, at 5% consistency, using 50 mM potassium acetate buffer at pH 5, 60 U g^−1^ odp laccase and 8% odp of Tempo for 18 h, according to previous works \[[@CR7], [@CR8]\].

All the enzymatic treatments were conducted in polyethylene bags that were placed in a laboratory water bath. After treatment, liquors were recovered and the resulting pulp was extensively washed as reported elsewhere for eucalyptus pulp \[[@CR30]\] in order to remove the enzymes and their degradation products. In the case of L_Tempo treatments, pulp was also washed with ethanol. Control treatments with potassium acetate buffer and ammonium acetate buffer were also performed at the same application conditions but without the addition of enzymes. They were named "C~K~" and "S~K~".

Effects on effluent properties {#Sec6}
------------------------------

Released cellooligosaccharides were quantified by the dinitrosalicylic (DNS) reagent method and analysed by thin-layer chromatography (TLC) and HPAEC-PAD (high-performance anion-exchange chromatography with pulsed amperometric detection). For reducing sugar quantification, 100 μL of DNS was added to 100 μL samples and mixtures were incubated at 100 °C for 5 min. Then, 40 μL of reaction mixtures was placed in ELISA plates, 260 μL of distilled water was added, and the absorbance at 540 nm was measured. Samples were analysed in triplicate. A standard curve of glucose was used to calculate the glucose reducing sugar equivalent of the different samples \[[@CR31]\].

For TLC analysis 10--15 μL of samples was applied on a silica gel plate (Merck, Germany) constituting the solid phase. 10 μL of an oligomer standard mixture containing cellooligosaccharides at a concentration of 20 mg mL^−1^ was applied as migration standards. The mobile phase was a mixture of chloroform, acetic acid and H~2~O in 6:7:1 ratio, respectively. The migration was repeated twice, and the silica gel plate was then sprayed (Fungilab S.A., Spain) with a developing solution, consisting of 5% H~2~SO~4~ in ethanol. Finally, the plate was heated in the oven at 100 °C for 5 min, where the spots corresponding to different cellooligosaccharides were visualized \[[@CR31]\]. For HPAEC-PAD sample preparation, after removing insoluble substrates by centrifugation, supernatants were centrifuged and diluted in water 1/20 and analysed by HPAEC-PAD using Dionex GS50, gradient pump, Dionex AS50 Autosample and electrochemical detector Waters 2465. In brief, 40-μL samples were injected on a CarboPac PA1 2 × 250 mm analytical column (Dionex). Cellooligosaccharides were eluted at 0.25 mL min^−1^ using a stepwise linear gradient from 100% eluent A (0.1 M NaOH) towards 10% eluent B (0.6 M NaOAc in 0.1 M NaOH) 10 min after injection and to 40% eluent B 15 min after injection, followed by a 5-min exponential gradient to 100% B. The column was reconditioned between each run by running initial conditions for 10 min. Standards were generated using 1, 2, 4 and 8 μg mL^−1^ cellobiose and cellobionic acid \[[@CR17]\].

Pulp characterization {#Sec7}
---------------------

The morphological properties of the fibres (viz. length and width) and the content in fines of the pulp samples were determined in accordance with TAPPI T 271 on a Metso kajaani FS300 fibre analyser. All samples were analysed in duplicate. Viscosity was determined according to ISO 5351:2010. The degree of polymerization (DP) was calculated from the intrinsic viscosity (*Ƞ*), using the equation of (SCAN-CM15:88): DP0.085 = 1.1 × \[*Ƞ*\]. Carboxyl groups were determined by measuring methylene blue adsorption onto cellulose fibres according to Davidson \[[@CR32]\]. To measure aldehyde groups samples were further oxidized with NaClO~2~ for selective conversion of aldehyde groups into carboxyl groups at room temperature for 48 h. The carboxyl content was determined with the above-described method. The carboxyl groups formed by effect of NaClO~2~ oxidation were assumed to derive from aldehyde groups originally present in the pulp. Three measures per sample were performed, and the 95% confidence interval was calculated.

High-pressure homogenization {#Sec8}
----------------------------

Prior to fibrillation, 2 g of oven-dried pulp (odp) at 1% consistency was disintegrated for 1 min at 11,200 rpm with a homogenizer (Homogenizing System UNIDRIVE X1000). Then, samples were diluted until 0.5% consistency and homogenized through the PANDA GEA 2000 homogenizer by 5 passes at 300 bar and 10 passes at 900 bar.

The yield of fibrillation (Eq. [1](#Equ1){ref-type=""}) was calculated after centrifuging 10 mL of a sample at 0.1% consistency at 2200×*g* for 20 min, removing the supernatant (containing the nanofibrillated fraction) and drying pellet (*C*) at 85 °C until constant weight.$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\text{Yield}} = \left( {1 - \frac{{C \left( {\text{g}} \right)}}{{0.01 {\text{g}}}}} \right) \times 100 \%$$\end{document}$$

Transmittance measurements were taken on samples with 0.1% of solid content. The sample was introduced in quartz cuvettes, and the transmittance was obtained with a T92 + UV spectrophotometer (PG instruments) set in the range between 400 and 800 nm. Milli-Q water was used as blank.

Fibre morphology and DP were measured as previously described in pulp samples. Electrophoretic mobility of aqueous suspensions (zeta potential) was determined using a Zetamaster model ZEM (Malvern Instruments, UK). Data were averaged over 10 measurements. All samples were analysed at room temperature.

Film characterization {#Sec9}
---------------------

After fibrillation, films with a grammage around 45--50 g m^−2^ were obtained by the film casting technique \[[@CR33]\]. Their optical and physical--mechanical properties were determined in accordance with the standards in brackets as follows: transparency (22891:2013), apparent density (ISO 534:2005), Bekk smoothness (5627:1995), and dry and wet zero-span index (ISO 15361:2000). Fibres zero-span tensile index was determined in a Zero-span 1000 Pulmac tester. For analysis of wet zero-span index, films were previously soaked in distilled water for 120 s.

Barrier properties to air and water were also analysed. Air permeance was measured with Bekk equipment. Water impermeability was measured by the water drop test (WDT) according to TAPPI standard T835 om-08. The WDT involved placing a drop of deionized water on the surface of paper and recording the time needed for complete absorption, which was signalled by vanishing of the drop specular gloss. Ten measurements per treated film sample were made and averaged. Six measures per sample were performed, and the 95% confidence interval was calculated.

The crystallinity index (CrI) of different cellulosic substrates was measured by XRD (X-ray powder diffraction). Samples were dried directly on an aluminium plate of 32 mm of diameter and 3.0 mm of thickness, that were mounted in standard sample holders for bulk samples of thickness ≤ 7 mm (PW1812/00) by means of plasticine. A PANalytical X'Pert PRO MPD Alpha1 powder diffractometer in Bragg--Brentano *θ*/2*θ* geometry of 240 mm of radius with Cu Kα1 radiation (*λ* = 1.5406 Å) at 45 kV and 40 mA, focalizing Ge (111) primary monochromator, with sample spinning at 2 revolutions per s, fixed divergence slit of 0.25º, was used. The measurement range (2*θ*) was from 2º to 50º with a step size of 0.033º and measuring time of 100 s per step. To calculate the CrI of cellulose from the XRD spectra, the peak height method used elsewhere was applied \[[@CR34]\].

Total crystallinity index (TCI) was measured using Fourier transform infrared (FTIR) spectra as previously described \[[@CR35]\].

Morphological characterization of film surface was performed by field emission scanning electron microscopy (FESEM) (JSM 7100 F) using a LED filter and a backscattered electron detector (BED).

Results and discussion {#Sec10}
======================

The initial cotton linters were long fibres, with an average length of 0.47 mm, which had been hornified (stiffened) during the drying inherent to their production. These traits made them difficult to process because they usually clog in the high-pressure homogenization apparatus. For this reason, they were mechanically refined by beating in a valley mill, which reduced fibre length to 0.25 mm and facilitated their homogenization.

The effect of enzymes on the refined cotton linters was firstly assessed on the properties of effluents released. Then, modifications produced by enzymes in fibre morphology and cellulose were analysed, together with their effect on the fibrillation improvement. Finally, the optical, physical and barrier properties of the films casted from the treated fibres were evaluated (Fig. [1](#Fig1){ref-type="fig"}) and compared with films obtained from commercial NFCs.Fig. 1General schema of the experimental work performed

Effect of enzymes on sugar release {#Sec11}
----------------------------------

Cotton linters were treated with different cellulases: endoglucanase Cel9B (C~9~) and commercial cellulases C~50~, C~F~ or C~ll~, and the amount of neutral sugars released was analysed by DNS (Table [1](#Tab1){ref-type="table"}). C~9~ and C~F~ produced a similar sugar release, much lower than that released by C~50~ and C~ll~. When C~F~ and C~ll~ were applied in the same treatment (C~mix~) the sugar release was the same as that with C~ll~ alone.Table 1Neutral sugar and cellobionic acid release produced by the enzymatic pretreatmentsSugar release (mg g^−1^ pulp)Cellobionic acid release (mg g^−1^ pulp)C~k~00C~9~22 ± 40C~50~90 ± 100C~F~34 ± 50C~ll~141 ± 110C~mix~126 ± 130S--0.8 ± 0.2SC~mix~242 ± 26.1 ± 1.1C~k~ (cellulase control treatment), C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~)

TLC analysis showed that C~9~ released mainly glucose and cellobiose (Additional file [1](#MOESM1){ref-type="media"}), being cellobiose the most abundant cellooligosaccharide released in accordance with its processive endoglucanase activity \[[@CR27]\]. Similar product pattern was reported by Garcia-Ubasart et al. \[[@CR36]\] when treating flax pulp with this enzyme. Commercial cellulases released a wider pattern of products from cotton linters, neutral sugars from glucose to cellotetraose, without noticeable differences among the enzymes (Additional file [1](#MOESM1){ref-type="media"}).

The action of LPMO, SamLPMO10C (S), was analysed determining the production of oxidized sugars in the effluents by HPAEC-PAD. S treatment released cellobionic acid and other aldonic acid oligosaccharides of higher molecular weight, together with a small fraction of neutral sugars (Fig. [2](#Fig2){ref-type="fig"}). Although the amount of cellobionic acid released to the effluents was low (Table [1](#Tab1){ref-type="table"}), the ability of SamLPMO10C to oxidize cotton linters was demonstrated. Our results are in accordance with the production of C1-oxidized oligosaccharides from phosphoric acid-swollen cellulose (PASC) by SamLPMO10C, which was also able to release aldonic acids from flax fibres \[[@CR17], [@CR29]\]. On the contrary, in other reported works, the production of aldonic acids when an LPMO belonging to the AA9 family was applied to softwood kraft pulp was not observed \[[@CR15]\].Fig. 2Sugar release produced by C~k~ (control treatment), C~9~ (Cel9B), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO) and SC~mix~ (LPMO and C~mix~) analysed by HPAEC-PAD. Enzymatic treatments with cellulases were performed at 50 °C, pH 5, for 18 h with 10 U g^−1^ odp of enzyme. (In the case of C~mix~ the enzymatic dose was 20 U g^−1^ odp.) Enzymatic treatment with LPMO (S) was performed at 50 °C, pH 6, for 72 h with 4 mg of enzyme g^−1^ odp in the presence of ascorbic acid and hydrogen peroxide

LPMO have extensively been reported to promote the efficiency of cellulases \[[@CR12]\], and in our hands, an increase in sugar release by the combined SC~mix~ treatment, that was twice more than in the single treatment with cellulases, C~mix~, confirmed this statement. Furthermore, the cellobionic acid release in the SC~mix~ treatment was eightfold higher than with S, suggesting that the oxidized fractions of cellulose created during S were cleaved and released to the effluent during the C~mix~ treatment (Table [1](#Tab1){ref-type="table"}). C~mix~ treatment was not applied after L_Tempo treatment since some authors state that cellulose oxidation produced by L_Tempo impairs the action of these enzymes \[[@CR37]\].

Effect of enzymes on fibre morphology and cellulose modification {#Sec12}
----------------------------------------------------------------

The refined cotton linter fibres used were short (0.25 mm average length) and had a large amount (more than 50%) of fibres lower than 0.2 mm (fines) (Additional file [2](#MOESM2){ref-type="media"}a). Enzyme treatment changed the morphology and size distribution of fibres. Cellulases acted on the longer fibres (around 0.2--7.6 mm) creating high amounts of fines, which showed the highest increase in their shortest fraction, fibres lower than 0.1 mm. Among the cellulases, C~9~ produced the smaller morphology change. It slightly reduced fibre length, with only 3% increase in the fine content, and it did not produce a significant effect on fibre width (Table [2](#Tab2){ref-type="table"}). Fibre degradation by C~ll~ was higher (16% fines increase), in accordance with its higher sugar release. Although C~50~ released more quantity of sugars than C~F~, they produced a similar fibre degradation (fines increased in 10%). Combined cellulase treatment (C~mix~) produced the highest increase in fines content (31%), although sugar release was not increased in the combined treatment. Fibre width was slightly reduced by C~ll~ and C~mix~ probably because of the degradation of the surface fibrillation of fibres.Table 2Effects of enzymatic pretreatments on fibre morphology and on mechanical fibrillationPulpAfter mechanical fibrillationFibre width (μm)Fibre length (mm)Fines (%)Yield (%)Transmittance at 700 nm (%)R25.50.2552.5802.4 ± 0.5C~k~24.90.2553.703.2 ± 0.7C~9~24.70.2455.4307.8 ± 0.3C~50~24.20.2259.25045.6 ± 1.2C~F~24.80.2259.01032.9 ± 0.9C~ll~23.20.2162.18027.9 ± 0.8C~mix~21.20.1870.4418 ± 152.7 ± 0.9S~k~25.10.2552.803.8 ± 0.5S23.30.2161.582.5 ± 24.3 ± 0.4SC~mix~23.20.1673.523 ± 151.8 ± 0.5L_Tempo23.00.2357.1310 ± 24.5 ± 0.6Com------11 ± 25.7 ± 0.7The confidence interval of the fibre width, fibre length and fines was less than 2% in all casesR (initial refined pulp), C~k~ (cellulase control treatment), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase) C~9~ (Cel9B), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S~k~ (LPMO control treatment), S (LPMO), SC~mix~ (LPMO and C~mix~), L_Tempo (Laccase_Tempo treatment) and Com (commercial NFC)

Whereas hydrolytic treatments with cellulases are well known to act on fibre morphology \[[@CR38]\], little knowledge exists about the fibre modification produced by oxidative treatments, particularly with LPMO enzymes. Interestingly, the two oxidative treatments performed affected fibre morphology, reducing its fibre length and width and consequently increasing fines content (Table [2](#Tab2){ref-type="table"}). The increase in fines content was more pronounced with S (16%) than with L_Tempo (6%). These results contrast with those reported by Aracri et al. \[[@CR8]\] reporting that no effect on fines content was produced by L_Tempo treatment of sisal pulps. Finally, SC~mix~ treatment produced the highest increase in fines (37%) and a large amount of fines lower than 0.1 mm (42%) (Additional file [2](#MOESM2){ref-type="media"}b), in agreement with the highest sugar release of SC~mix~ treatment, confirming that fibre degradation by cellulases was boosted by LPMO action. These results are in accordance with the proposed mechanism of LPMO that create nicking points where the cohesion of the fibres was decreased, improving the attack of cellulases \[[@CR15]\].

Changes in cellulose polymerization were assessed via intrinsic viscosity measurements (Fig. [3](#Fig3){ref-type="fig"}). Similarly to what has been reported \[[@CR26], [@CR39]\] all the tested cellulases decreased DP. In correlation with the effects of the cellulases on fibre morphology and sugar release described above, C~9~ produced lower cellulose degradation (52% decrease in DP) than the commercial cellulases applied (around 73--79%). A similar cellulose depolymerization was observed by Qing et al. in 2013 \[[@CR24]\] when C~F~ and C~ll~ were applied to a bleached eucalyptus kraft pulp at lower enzymatic doses. Contrary to our results, previous authors reported that DP of softwood and flax pulps was not significantly affected by endoglucanase C~9~ \[[@CR36], [@CR38]\], but in our case, the higher cellulose degradation produced by C~9~ in cotton linters can be due to the longer treatment applied (18 h vs. 1--2 h in previous works).Fig. 3Effect of enzymatic treatments on cellulose degree of polymerization. R (initial refined pulp), C~k~ (control treatment), C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~) and L_Tempo (Laccase_Tempo treatment). Enzymatic treatments with cellulases were performed at 50 °C, pH 5, for 18 h with 10 U g^−1^ odp of enzyme. (In the case of C~mix~ the enzymatic dose was 20 U g^−1^.) Enzymatic treatment with LPMO (S) was performed at 50 °C, pH 6, for 72 h with 4 mg of enzyme g^−1^ odp in the presence of ascorbic acid and hydrogen peroxide. Enzymatic treatment with Laccase_Tempo was performed at room temperature, pH 5, for 18 h at 60 U g^−1^ odp of laccase and 8% odp of Tempo

Concerning the oxidative treatments, cellulose DP was highly affected by L_Tempo, while S treatment produced a small decrease (Fig. [3](#Fig3){ref-type="fig"}). Depolymerization of cellulose by L_Tempo has been described to be produced by active species such as hydroxyl radicals formed in situ by side reactions of the hydroxylamine structure with oxygen during the oxidative treatment \[[@CR40]\]. Moreover, the presence of aldehyde groups produced by L_Tempo treatment can provide underestimation of viscosity values. These aldehyde groups can give depolymerization reactions through β-elimination during the viscosity determination method, with cupriethylenediamine at alkaline conditions. To avoid this problem, viscosity was also measured after treating the pulp samples with sodium borohydride (borohydride viscosity) in order to inactivate carbonyl groups by reduction to hydroxyl groups \[[@CR8]\]. Borohydride viscosity was measured in all the samples (data not shown) obtaining DP values similar to those shown in Fig. [3](#Fig3){ref-type="fig"}, with the exception of L_Tempo sample that showed an increased DP, indicating that aldehyde groups were formed in this treatment. However, even after the reductive treatment, DP of L_Tempo sample was low (572), indicating degradation of cotton linters cellulose by L_Tempo, similarly to previous results described for sisal pulps, although with less intense degradation \[[@CR8]\]. On the contrary, the low depolymerization produced by LPMO suggested that this enzyme affects fibre morphology without significantly degrading cellulose. Villares et al. also reported a slight decrease in DP by LPMO although fibre morphology was not affected \[[@CR15]\]. Interestingly, despite the high fibre modification and cellobionic acid release by the S treatment, cellulose was not significantly degraded. Subsequent treatment with the cellulase mixture, SC~mix~, did not increase cellulose depolymerization by cellulases.

Finally, the creation of functional groups on cellulose was evaluated by measuring the carboxyl and aldehyde content of fibres. Results showed a significant increase of these groups only with the oxidative treatments, where the L_Tempo treated pulps exhibited the highest content (Table [3](#Tab3){ref-type="table"}), as also appreciated by FTIR spectra (Additional file [3](#MOESM3){ref-type="media"}). A different mechanism of creating COOH groups was produced among the oxidative treatments: whereas L_Tempo oxidized cellulose as a result of the conversion of C6 primary hydroxyl groups in cellulose via an aldehyde group \[[@CR41]\], S created COOH through oxidation of alcohol in C1 position \[[@CR29]\]. A small fraction of aldehydes were also produced during L_Tempo in accordance with the previous results on borohydride viscosity. The modest increase in carboxyl group content provided by the L_Tempo system in comparison with other works under the conditions used is probably due to the absence of added oxygen during the treatment \[[@CR7]\].Table 3Effects of enzymatic pretreatments in the creation of COOH and CHO groups in cellulose fibresCOOH (µmol g^−1^ odp)CHO (µmol g^−1^ odp)R18 ± 40C~k~13 ± 10C~9~20 ± 11 ± 1C~mix~17 ± 31 ± 1S29 ± 40SC~mix~20 ± 10L_Tempo85 ± 94 ± 1R (initial refined pulp), C~k~ (cellulase control treatment), C~9~ (Cel9B), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~) and L_Tempo (Laccase_Tempo treatment)

The oxidative cleavage of cellulose by S treatment leads to the formation of a small quantity of aldonic acids (COOH groups) at the C1 position. However, a fraction of these COOH groups disappeared when C~mix~ was applied after S, suggesting they were removed. This result is in accordance with the higher cellobionic acid release in the SC~mix~ treatment previously shown, corroborating that the oxidized fractions of cellulose created by S activity were cleaved by the C~mix~ treatment.

Effect of enzymes on fibrillation improvement {#Sec13}
---------------------------------------------

Enzyme-treated samples were homogenized at high pressure, and their properties were analysed. NFC was only obtained in the hydrolytic treatment with the cellulase mixture C~mix~ (simultaneous application of C~F~ and C~ll~), while none of the other cellulase treatments gave significant amount of NFC (Table [2](#Tab2){ref-type="table"}). Nanofibrillation of this sample may have been promoted by its low fibre length (70% of fines) or by the decrease of hornification produced by these cellulases \[[@CR19], [@CR42]\]. A lower yield of NFC was obtained with oxidative S and L_Tempo treatments. In these treatments, nanofibrillation was probably stimulated by the presence of COOH groups, as it has been reported \[[@CR43], [@CR44]\]. The highest yield of NFC (23%) was produced with the SC~mix~ pretreatment, in concordance with the best performance observed in the other parameters evaluated, where a higher effect of cellulases after an LPMO treatment was achieved. This yield increase produced by LPMO can be related to the introduction of nicks in the most crystalline regions of cellulose molecules (as suggested by Villares et al. and Valenzuela et al. \[[@CR15], [@CR17]\]), rather than to the small increase of COOH that are left on fibres after cellulase treatment. Recently, it has been reported the nanofibrillation of flax pulp after a sequential pretreatment of SamLPMO10C and C~9~, obtaining a similar yield of 24% \[[@CR17]\]. Remarkably, the NFC yield obtained with C~mix~ and SC~mix~ was higher than the NFC content of a commercial nanocellulose (Table [2](#Tab2){ref-type="table"}).

Despite the fact that in some samples no NFC was obtained, other parameters were measured in order to analyse the fibrillation improvements produced by the enzymes. For example, transmittance is a simple mean to get an approximate idea about the width of the ensuing fibrils. In fact, when light passes through a medium containing randomly dispersed particles, it is scattered by the particles causing a reduction in the transparency degree, as previously reported \[[@CR24]\]. Transmittance at 700 nm was strongly improved by hydrolytic treatments in all the samples (Table [2](#Tab2){ref-type="table"}), suggesting a decrease in the amount of non-fibrillated and partially fibrillated fractions responsible for light scattering phenomenon. In accordance with the NFC yield results, the highest improvement in transmittance was produced with C~mix~ and SC~mix~ samples. Although fibres of lower length were created with C~ll~ treatment, a higher transmittance value was obtained with C~50~, followed by C~F~, C~ll~ and C~9~. Concerning the oxidative treatments, they only improved transmittance to less than 5%. The higher carboxyl content of L_Tempo sample did not produce a significant increase in transmittance, in accordance with the observations of Besbes et al., 2011, who reported that COOH content has to be higher than 300 μmol g^−1^ odp to produce a significant increase in transmittance \[[@CR44]\].

Zeta potential is a measure of the magnitude of the electrostatic or charge repulsion/attraction between particles and is one of the fundamental parameters known to affect stability. All the samples obtained after mechanical fibrillation had a Z potential around − 30 mV, that indicates that there is no agglomeration, which means sufficient mutual repulsion resulting in colloidal stability. This value was slightly increased with the oxidative treatments to − 40 mV probably due to the COOH groups (Additional file [4](#MOESM4){ref-type="media"}). However, it was reduced in the SC~mix~ treatments, correlating again the removal of the LPMO-produced COOH groups by the cellulase treatment. A similar result has been reported in NFC from flax and bleached kraft pulp \[[@CR16], [@CR17]\].

Although it has been reported that DP can be reduced during fibrillation \[[@CR24], [@CR25], [@CR45]\], in our results DP was not affected after the pass through the high-pressure homogenizer (Fig. [3](#Fig3){ref-type="fig"}). In fact, only in the S sample the DP decreased. Maybe the oxidation of glycosidic linkages during the treatment with LPMO rendered cellulose more susceptible to be degraded during fibrillation. Finally, it has to be pointed out the low DP of C~mix~ and SC~mix~ samples, indicating that cellulose chains were only formed by ≈ 300 glucose units. This value was only slightly higher than in cellulose nanocrystals (≈ 200 glucose units) obtained from cotton linters \[[@CR35]\].

Effect of enzymes on film properties {#Sec14}
------------------------------------

Films of \~ 45 μm thickness were prepared, and their optical, physico-mechanical and barrier properties were measured (Table [4](#Tab4){ref-type="table"}). Crystallinity of films was determined by XRD. It was high in all the samples (around 90%), as expected for cotton linters, although they had suffered multiple passes through the homogenizer, a process that has been reported to reduce crystallinity \[[@CR46]\]. Values obtained are similar to that reported by Hideno et al. in 2016 and higher to that obtained by Saito et al. in 2006 \[[@CR19], [@CR47]\]. Cellulases treatment slightly increased crystallinity of films, probably due to their action on the amorphous zones of cellulose more susceptible to be attacked by these enzymes \[[@CR48]\], a phenomenon observed also when commercial cellulases were applied to bleached wood pulps \[[@CR24], [@CR26]\]. Crystallinity is also an important parameter that affects the action of LPMO enzymes, where, on the contrary, higher crystalline cellulose seems to be a better substrate to be oxidized \[[@CR17], [@CR49]\]. Interestingly, in our experiments, this property was not negatively affected by the S treatment, similarly to what has been reported for NFC from flax pulps \[[@CR17]\]. The other oxidative treatment, L_Tempo, did not affect to this property neither, as previously reported \[[@CR47]\]. The lower DP produced with the enzymatic treatments did not affect cellulose crystallinity. This property was also measured from the FTIR spectra obtaining the total crystallinity index (TCI) (Additional file [3](#MOESM3){ref-type="media"}). It had a value around 1.2, without significant differences between samples, in accordance with values obtained by XRD.Table 4Effects of enzymatic pretreatments in crystallinity and physical properties of the films obtained after mechanical fibrillationCrystallinity (%)Transparency (%)Density (g cm^−3^)Bekk smoothness (s)Upper faceLower faceUpper faceR89.1 ± 0.127.8 ± 0.928.0 ± 0.10.8 ± 0.115 ± 4C~k~88.5 ± 0.535.9 ± 3.736.0 ± 3.70.8 ± 0.110 ± 1C~9~90.3 ± 0.254.3 ± 0.452.9 ± 0.21.0 ± 0.129 ± 5C~50~91.4 ± 0.183.4 ± 0.583.2 ± 1.11.1 ± 0.046 ± 6C~F~91.0 ± 0.374.9 ± 0.773.8 ± 0.51.1 ± 0.153 ± 2C~ll~91.4 ± 0.275.7 ± 0.974.7 ± 0.81.1 ± 0.034 ± 2C~mix~92.4 ± 0.286.1 ± 0.585.9 ± 1.50.9 ± 0.050 ± 6S~k~89.0 ± 0.252.9 ± 1.851.5 ± 2.30.8 ± 0.111 ± 2S90.5 ± 0.560.0 ± 0.460.9 ± 0.30.9 ± 0.010 ± 2SC~mix~92.5 ± 0.185.6 ± 0.585.3 ± 0.51.1 ± 0.153 ± 1L_Tempo89.8 ± 0.138.9 ± 3.237.8 ± 1.80.9 ± 0.114 ± 3Com81.1 ± 0.266.2 ± 0.766.7 ± 0.41.0 ± 0.019 ± 3R (initial refined pulp), C~k~ (cellulase control treatment), C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S~k~ (LPMO control treatment), S (LPMO), SC~mix~ (LPMO and C~mix~), L_Tempo (Laccase_Tempo treatment) and Com (commercial NFC)

Transparency of films was determined, showing agreement with transmittance of homogenized suspensions previously shown, and no significant differences were found between upper and lower faces of films (Table [4](#Tab4){ref-type="table"}). The highest transparency was achieved with C~mix~ and SC~mix~, where the increase in transparency was around 50 points (see Additional file [5](#MOESM5){ref-type="media"}). C~50~, C~F~, C~ll~ and C~9~ increased this property in 47, 39, 39 and 18 points, respectively. Transparency obtained with C~mix~ treatments was similar to that reported by Hideno et al. in 2016 with cotton linters and cellulase, and also to Chen et al. in 2014 in a NFC/acrylic resin composite sheet \[[@CR19], [@CR21]\]. S showed a lower increase in transparency of 7 points, being these films of higher transparency than L_Tempo films. Interestingly, crystallinity and transparency of films from cellulase treatments were higher than those of the films made from commercial NFC.

The density of the films obtained (Table [4](#Tab4){ref-type="table"}) was comparable to that of films obtained from bacterial cellulose and considerably higher than papers from wood fibres \[[@CR50]\]. Enzymatic treatment with cellulases produced the films of higher density probably due to the lower fibre size. Smoothness of non-treated films (R) was similar to films obtained from commercial NFC (Table [4](#Tab4){ref-type="table"}). Whereas smoothness was not affected by the oxidative treatments, this property was significantly increased with all the cellulases applied, particularly with C~50~, C~F~ and C~mix~ treatments. No differences were appreciated if C~mix~ was applied after S. Increased values of smoothness were obtained in the lower face of films (data not shown). The high values of smoothness and transparency of the obtained films give them the potential to be applied for printed electronics \[[@CR51]\].

Mechanical resistance of films was determined (Fig. [4](#Fig4){ref-type="fig"}). Non-treated films (R) showed dry zero-span index of 153 Nm g^−1^, similar to that of commercial NFC films (185 Nm g^−1^) and higher than that of the paper from unbleached kraft pulp reinforced with cotton linters NFC \[[@CR52]\] (8 Nm g^−1^). As previously said, DP was affected by the action of enzymes, and consequently, this could affect the physical properties of the resulting films. Interestingly, although cellulose was partially degraded with all the treatments, dry zero-span index of C~9~, S and L_Tempo samples was not significantly affected, while a significant reduction in this property was produced with all the commercial cellulases used (Fig. [4](#Fig4){ref-type="fig"}). In order to evaluate the resistance offered by a single fibre, wet zero-span index was also measured. The resistance in all the samples was reduced around 55--84%, including the one from commercial NFC. Whereas non-treated films had values around 70 Nm g^−1^, this value was reduced to 25 Nm g^−1^ in all the enzymatic treated samples, without significant differences between them. Films from commercial NFC had a slightly higher wet zero-span index (42 Nm g^−1^). These wet zero-span values obtained were significantly lower than those reported for bacterial cellulose films (around 100 Nm g^−1^) \[[@CR50]\], probably as a result of the higher crystallinity of bacterial nanocellulose.Fig. 4DP of cellulose *vs.* mechanical resistance of NFC films. R (initial refined pulp), C~k~ (control treatment), C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~) and L_Tempo (Laccase_Tempo treatment). Enzymatic treatments with cellulases were performed at 50 °C, pH 5, for 18 h with 10 U g^−1^ odp of enzyme. (In the case of C~mix~ the enzymatic dose was 20 U g^−1^.) Enzymatic treatment with LPMO (S) was performed at 50 °C, pH 6, for 72 h with 4 mg of enzyme g^−1^ odp in the presence of ascorbic acid and hydrogen peroxide. Enzymatic treatment with Laccase_Tempo was performed at room temperature, pH 5, for 18 h at 60 U g^−1^ odp of laccase and 8% odp of Tempo

Barrier properties to air and water of the obtained films were also measured. Air permeability was measured by the Bekk method (Fig. [5](#Fig5){ref-type="fig"}). The non-treated sample (R) had a similar value than films from commercial NFC. Interestingly, cellulases strongly increased the seconds that the air needed to pass through the films, i.e. decreased permeability. The most notable effect was produced with C~F~, C~50~ and SC~mix~ followed by C~ll~, C~mix~ and C~9~. On the other hand, oxidative treatments did not produce significant effects. Permeability of films from cellulase-treated samples was threefold higher than that of commercial NFC films, indicating that a strongly closed structure was formed after the enzymatic treatments. The increased fine content and fibrillation obtained with the cellulase treatments are consistent with an increased cohesion between fibre surfaces and responsible for the decreased paper permeability. These results are consistent with those of Cadena et al. who found cellulase treatments to reduce the air permeance of paper \[[@CR53]\]. Similar to smoothness, permeability was strongly decreased in the lower face of the film (data not shown).Fig. 5Barrier properties to air (Bekk permeability) and water (water drop test) of NFC films. R (initial refined pulp), C~k~ (cellulase control treatment), C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S~k~ (LPMO control treatment), S (LPMO), SC~mix~ (LPMO and C~mix~), L_Tempo (Laccase_Tempo treatment) and Com (commercial NFC). Enzymatic treatments with cellulases were performed at 50 °C, pH 5, for 18 h with 10 U g^−1^ odp of enzyme. (In the case of C~mix~ the enzymatic dose was 20 U g^−1^.) Enzymatic treatment with LPMO (S) was performed at 50 °C, pH 6, for 72 h with 4 mg of enzyme g^−1^ odp in the presence of ascorbic acid and hydrogen peroxide. Enzymatic treatment with Laccase_Tempo was performed at room temperature, pH 5, for 18 h at 60 U g^−1^ odp of laccase and 8% odp of Tempo

Film permeability was intensely related to the barrier property to water, measured by the WDT (Fig. [5](#Fig5){ref-type="fig"}). All the films from enzyme-treated samples showed an increased water impermeability although the effect was more noticeable with cellulase treatments, which showed a maximum value of 3150 s in C~F~ sample. In spite of the high impermeability compared to current cellulosic papers, it was lower than that provided by bacterial cellulose films (4000 s) \[[@CR50]\]. Interestingly, non-treated films had lower permeability to water than commercial NFC films, maybe because of the higher crystallinity of the cotton linters used.

Finally, film surface morphology was analysed by FESEM (Fig. [6](#Fig6){ref-type="fig"}). Non-treated films (R) showed fibres of different lengths and fibre widths with fibrillation. A highly entangled nano- and/or micro-fibre network was observed in enzyme-treated samples, similarly to that reported by Hu et al. in 2018 and Tarrés et al. in 2017 \[[@CR16], [@CR39]\]. Surface morphology of films demonstrates that the enzymatic treatments performed boosted mechanical delamination, since those films showed a compact structure and their structure was difficult to be visualized. Moreover, in SC~mix~ films a thin layer of nanofibres surrounding bigger fibres was appreciated.Fig. 6FESEM images of NFC films from non-treated and enzymatically treated samples. R (initial refined pulp), C~9~ (Cel9B), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~) and L_Tempo (Laccase_Tempo treatment). Enzymatic treatments with cellulases were performed at 50 °C, pH 5, for 18 h with 10 U g^−1^ odp of enzyme. (In the case of C~mix~ the enzymatic dose was 20 U g^−1^.) Enzymatic treatment with LPMO (S) was performed at 50 °C, pH 6, for 72 h with 4 mg of enzyme g^−1^ odp in the presence of ascorbic acid and hydrogen peroxide. Enzymatic treatment with Laccase_Tempo was performed at room temperature, pH 5, for 18 h at 60 U g^−1^ odp of laccase and 8% odp of Tempo

Although the presence of NFC material was not detected in films from individual cellulases and oxidative enzymes, the film properties obtained clearly show that these treatments improved fibrillation. Moreover, according to the optical, physical and barrier properties obtained, the films from enzymatically treated cotton linters seem very promising to obtain biomaterials that could replace petrol-based products.

Conclusions {#Sec15}
===========

Four hydrolytic enzymes (cellulases) were applied on cotton linters, affecting fibre morphology and degrading cellulose differently. Improved mechanical fibrillation and 18% NFC yield were obtained with a cellulase mixture (C~mix~). Application of oxidative enzymes (LPMO and L_Tempo) introduced COOH groups into cellulose. The amount of COOH groups created with L_Tempo allowed the production of NFC during mechanical fibrillation (10%). However, the smaller quantity of these groups introduced by LPMO was not enough to produce NFC. The main difference between the two oxidative treatments was that L_Tempo degraded cellulose, whereas LPMO had more effect on fibre degradation. LPMO (S) boosted the action of cellulases although the COOH groups created were released to effluents after the hydrolytic treatment. Films with high crystallinity (92%) and transparency (86%), increased smoothness, and high air and water barrier properties were obtained after cellulase treatment and mechanical fibrillation on cotton linters. The introduction of an LPMO treatment before cellulase mixture (SC~mix~ treatment) produced higher NFC yield (23%) without a further improvement in film properties.

Additional files
================

 {#Sec16}

**Additional file 1.** TLC analysis of sugars released during enzymatic treatments performed at 50 °C, pH 5, during 18 h with 10 U g^−1^ odp of enzyme (in the case of Cmix the enzymatic dose was 20 U g^−1^ odp). C~9~ (Cel9B), C~50~ (Sertec20 cellulase), C~F~ (Fibercare cellulase), C~ll~ (Celluclast cellulase), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast). M) size markers of glucose (G), cellobiose (G2), cellotriose (G3), cellotetraose (G4) and cellopentaose (G5). **Additional file 2.** Effect of the hydrolytic (a) and oxidative and C~mix~ (b) enzymatic pretreatments on fibre length distribution. The confidence interval of the length fibre distribution was less than 2% in all cases. **Additional file 3.** FTIR spectra of obtained films from control treatment (C~k~) and Laccase_Tempo treatment (L_Tempo). TCI was calculated from the ratio of the absorptions at 1372 and 2900 cm^−1^. The peak detected at 1750 cm^−1^ with L_Tempo correspond to the COOH groups. **Additional file 4.** Z potential values of the samples obtained after mechanical fibrillation. R (initial refined pulp), C~k~ (control treatment), C~9~ (Cel9B), C~mix~ (cellulase mixture consisting in Fibercare and Celluclast), S (LPMO), SC~mix~ (LPMO and C~mix~) and L_Tempo (Laccase_Tempo treatment). **Additional file 5.** Obtained films from cellulase control C~k~ (a) and C~mix~ (b) pretreatments after mechanical fibrillation.

NFC

:   nanofibrillated cellulose

R

:   refined cotton linter

C~k~

:   cellulase control treatment

LPMO

:   lytic polysaccharide monooxygenases

S

:   treatment with LPMO enzyme

S~k~

:   LPMO control treatment

C9

:   treatment with endoglucanase Cel9B

C~ll~

:   treatment with commercial cellulase Celluclast

C~F~

:   treatment with commercial cellulase Fibercare

C~50~

:   treatment with commercial cellulase from Sertec20

C~mix~

:   combined treatment with Celluclast and Fibercare

SC~mix~

:   treatment with LPMO enzyme followed by C~mix~

L_Tempo

:   treatment with laccase and Tempo

Com

:   commercial NFC

DNS

:   dinitrosalicylic

TLC

:   thin-layer chromatography

HPAEC-PAD

:   high-performance anion-exchange chromatography with pulsed amperometric detection

DP

:   degree of polymerization

WDT

:   water drop test

FESEM

:   field emission scanning electron microscopy

Odp

:   oven-dried pulp

XRD

:   X-ray powder diffraction

**Publisher\'s Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

We thank Serveis Científics Tècnics of University of Barcelona for HPAEC, XRD and FESEM analysis, and Novozymes^®^ and Sertec20 for supplying cellulases. Cotton linters were kindly supplied by CELESA.

CV conducted most of the experiments and wrote the manuscript; JP supervised the experiments and revised the manuscript; BR refined the cotton linters used and performed the analysis of fibre morphology; TV and PD took part in planning the study, checking the results and revising the work; JM performed FESEM experiments; SV obtained the LPMO and Cel9B enzymes, performed HPAEC and revised the manuscript. All authors read and approved the final manuscript.

This work was funded by the Spanish MINECO FILMBIOCEL CTQ2016-77936-R (funding also from the "Fondo Europeo de Desarrollo Regional FEDER"), MICROBIOCEL CTQ2017-84966-C2-1-R and BIOFABCEL CTQ2017-84966-C2-2-R projects. Special thanks are due to the consolidated research group AGAUR 2017 SGR 30 and to the Serra Húnter Fellow to Cristina Valls.

The data sets used and analysed during the current study are available from the corresponding author on reasonable request.

Not applicable.

Not applicable.

The authors declare that they have no competing interests.
