ORIGINAL RESEARCH article
Tree Plantation Systems Influence Nitrogen Retention and the Abundance of Nitrogen Functional Genes in the Solomon Islands
- 1Instituto de Ecología A.C., Red de Estudios Moleculares Avanzados, Xalapa, México
- 2Environmental Futures Research Institute, School of Natural Sciences, Griffith University, Nathan, QLD, Australia
- 3Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, QLD, Australia
- 4Environmental Futures Research Institute, School of Environment, Griffith University, Nathan, QLD, Australia
- 5Forestry College, Fujian Agriculture and Forestry University, Fuzhou, China
Tree mono-plantations are susceptible to soil nutrient impoverishment and mixed species plantations have been proposed as a way of maintaining soil fertility while enhancing biodiversity. In the Solomon Islands, mixed species plantations where teak (Tectona grandis) is inter-planted with a local tree species (Flueggea flexuosa) have been used as an alternative to teak mono-plantations and are expected to increase soil microbial diversity and modify microbial biogeochemical processes. In this study, we quantified the abundance of microbial functional genes involved in the nitrogen (N) cycle from soil samples collected in teak, flueggea, and mixed species plantations. Furthermore, we measured soil properties such as pH, total carbon (C) and total N, stable N isotope composition (δ15N), and inorganic N pools. Soil pH and δ15N were higher under teak than under flueggea, which indicates that intercropping teak with flueggea may decrease bacterial activities and potential N losses. Higher C:N ratios were found under mixed species plantations than those under teak, suggesting an enhancement of N immobilization that would help preventing fast N losses. However, inorganic N pools remained unaffected by plant cover. Inter-planting teak with flueggea in mixed species plantations generally increased the relative abundance of denitrification genes and promoted the enrichment of nosZ-harboring denitrifiers. However, it reduced the abundance of bacterial amoA (ammonia monooxygenase) genes compared to teak mono-plantations. The abundance of most denitrification genes correlated with soil total N and C:N ratio, while bacterial and archeal nitrification genes correlated positively with soil NH4+ concentrations. Altogether, these results show that the abundance of bacterial N-cycling functional guilds vary under teak and under mixed species plantations, and that inter-planting teak with flueggea may potentially alleviate N losses associated with nitrification and denitrification and favor N retention. Mixed plantations could also allow an increase in soil C and N stocks without losing the source of income that teak trees represent for local communities.
In tropical countries such as the Solomon Islands where deforestation rates are high, tree plantations are seen as a way to counteract soil degradation by restoring vegetation cover while decreasing the existing pressure on native forests (Wolfe et al., 2015). However, most established tree plantations are mono-plantations that can cause the same environmental problems as other monoculture systems, namely higher pest or disease occurrence and a modification of biogeochemical cycles and nutrient availability (Rachid et al., 2013). In order to maintain soil fertility and to enhance biodiversity, mixed plantations are being promoted and are expected to reduce plant competition for nutrients and increase soil carbon (C) and nitrogen (N) pools (Montagnini, 2000; Balieiro et al., 2008; Vigulu, 2015). Mixed-species systems are also expected to enhance soil microbial diversity by increasing the variety of carbon (C) based resources and the heterogeneity of spatial patterns in soil properties (Thevathasan and Gordon, 2004).
Many essential soil processes are primarily mediated by microbial communities (Bissett et al., 2013). Land management practices and changes in plant species composition are known to impact microbially driven processes in soils through the alteration of bacterial communities that will thereby modify nutrient availability to plants or nutrient losses from the ecosystem (Chen et al., 2003; Kourtev et al., 2003; Lindsay et al., 2010). Soil microorganisms play a central role in organic matter decomposition and in the cycling of major plant nutrients, including N (Hallin et al., 2009; Bissett et al., 2011). Yet, the importance of soil microbial diversity is often overlooked when establishing forest plantations (Levy-Booth and Winder, 2010). With the growing concerns about intensive mono-plantations comes an increasing interest in the management of soil fertility and soil bacterial communities to enhance tree growth and productivity (Lacombe et al., 2009). However, how the establishment of mixed plantations influences the functions of soil bacterial communities, as measured through microbial functional genes (MFGs), and how this relates to nutrient cycling remains to be understood.
In the Solomon Islands, teak (Tectona grandis L.f.) is often grown in smallholder plantations in order to rehabilitate the logged-over rainforests while providing a source of income to landowners (Vigulu, 2015). Teak is an economically important timber tree species grown in tropical and sub-tropical countries for its highly durable hardwood (Miranda et al., 2011). Mostly cultivated in monoculture plantations in 20–40 years rotation, its height can reach more than 20 m at maturity (Ladrach, 2009). Teak grows well on a broad range of soils but its growth has been reported to be optimal on deep and well-drained sandstones, with neutral or acid pH and high calcium, phosphorus, potassium, N, and organic matter contents (Kadambi, 1972). Currently, a new plantation system is being introduced in the Solomon Islands where teak is intercropped with a local tree species (Flueggea flexuosa Muell. Arg.) in order to overcome the reluctance of growers to thin pure teak stands. Flueggea, a small to medium tree typically 10–16 m tall, is traditionally used for house building and fencing in the Solomon Islands (Thomson, 2006). Flueggea was considered as a good candidate species for intercropping with teak as roots from both trees seem to occupy different soil depth. While teak has extensive horizontal and vertical roots and occupy a large portion of the soil volume, flueggea’s root system usually develops laterally, near the soil surface (Thomson, 2006; Vigulu, 2015).
The implementation of mixed-species systems is likely to influence nutrient cycling and the abundance of MFG associated with nutrient cycling (Rachid et al., 2013). Therefore, we aimed to determine the abundance of MFG involved in N cycling under teak mono-plantations, flueggea mono-plantations and mixed-species systems and evaluate differences in soil N pools due to tree cover. To accomplish this, we assessed the abundances of genes involved in nitrification (bacterial and archeal amoA), nitrogen fixation (nifH), and denitrification (narG, nirS, nirK, and nosZ) as indicators of microbial trait abundances (Wieder et al., 2013), and measured soil variables such as pH, total C and N content, NO3-, and NH4+. Soil N isotope composition (δ15N) was also analyzed as an indicator of N cycling rates (Högberg, 1997; Hietz et al., 2011; Reverchon et al., 2014). Finally, we discussed the existing relationships between MFG and soil characteristics in tree plantation soils.
Materials and Methods
The study site is located at Poitete, on Kolombangara Island, in the Western Province of the Solomon Islands (8° 05′ 16.33″ S and 157° 08′ 46.62″ E; Figure 1). The average annual temperature in the area is 28°C and the average rainfall is 3,600 mm, relatively evenly distributed throughout the year (Supplementary Figure S1). The soil at the study site is an Oxisol (Vigulu, 2015).
FIGURE 1. Smallholder tree plantation systems in Solomon Islands. The field site is located within the Pacific region (A), in Solomon Islands (B), at Kolombangara, Western Province (C), where teak (D), flueggea (E), and mixed-species (F) plantations were established in 2009.
Tree plantations were established in 2009 on land formerly covered with regenerated secondary tropical forests. Plantations consisted of mono-plantations of teak, mono-plantations of flueggea and mixed-species plantations (consecutive rows of teak and flueggea). All plantations were adjacent and therefore established on the same soil type, and had a planting density of 833 stems per hectare (spacing 4 m × 3 m). Three years after planting, teak trees were 16 m high with a diameter at breast height (dbh) of 19 cm, while flueggea trees were 12 m high with a dbh of 13 cm (Vigulu, 2015). The crowns and root systems of both trees were well-developed, with crown radius and root growth of teak trees being larger than those of flueggea (Vigulu, 2015).
Three replicated plots of 24 m × 24 m were established per plantation and a 24-m transect was drawn in the middle of each plot (nine transects in total). In the mixed-species plots, transects were drawn in order to cross alternatively rows of teak and rows of flueggea. Five sampling points were established 4 m apart along each transect. At each sampling point, soil samples were collected on the transect and 4 m perpendicularly on each side of the transect (three samples per sampling point), mixed, and bulked to constitute one composite sample per sampling point. All samples were collected from 0 to 15 cm soil layer, with a shovel. Soil samples were then immediately sieved (2 mm) and refrigerated until analysis (maximum of 5 days). The total number of samples per sampling was 45. Samplings were carried out in December 2012, May 2013, and December 2013.
Soil Chemical Analyses
Total C, total N, and δ15N of soil samples were determined by mass spectrometry (spectrometer GV Isoprime, Manchester, UK), following the procedure described in He et al. (2008). Soil pH was measured in water (1:5 ratio). Soil NH4+ and NO3- concentrations were determined by KCl extraction using a SmartChem 200 Discrete Chemistry Analyser as described in Bai et al. (2012). All soil NO3- concentrations were below detection level and were thus not included in the subsequent statistical analyses.
Quantification of the N-Cycle-Associated MFG
Soil DNA extractions were carried out within 1 week after sampling. DNA was extracted from 0.3 g of each soil sample using the MoBio Powersoil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA, USA) according to the manufacturer’s instructions, with a final elution in deionised water. The quantity and quality of the extracted DNA were verified with a Nanodrop spectrophometer (Thermoscientific) and diluted in deionised water (1:10).
The total bacterial community was quantified by quantitative PCR (qPCR) using the V3 hypervariable region of the 16S rRNA gene as a molecular marker. The abundances of functional genes narG, nirS, nirK, nosZ, nifH, and bacterial and archeal amoA (AOB and AOA respectively) were quantified using the primers and thermal cycling conditions described in Supplementary Table S1. Reactions were carried out in an Eppendorf Mastercycler ep realplex real-time PCR system (Eppendorf, Hamburg, Germany) in duplicate. Quantification was based on the fluorescence intensity of the SYBR Green dye (Takara) during amplification. Standard curves were obtained using 10-fold serial dilutions of plasmid DNA containing cloned narG, nirS, nirK, nosZ, nifH, amoA and 16S rRNA genes and spanning seven orders of magnitude. The 20 μL PCR mixture contained 10 μL of SYBR green PCR Master Mix [Takara SYBR Premix Ex Taq (Perfect Real Time)], 0.4 μL of each primer (10 μM) and approximately 8 ng DNA. Melting curves and agarose gels of PCR products were used at the end of each qPCR to check amplification specificity and purity of negative controls. Negative controls gave null or negligible values, and PCR efficiency for the different assays ranged from 90 to 99%. The presence of PCR inhibitors in DNA extracted from soil was estimated by a 1:10 soil DNA dilution; no inhibition was detected. All qPCR reactions were carried out immediately after DNA extraction.
The measured cycle threshold (Ct) values of standards quantification were calibrated by placing the threshold lines at the same level for each gene, to account for the different times at which samples from December 2012, May 2013, and December 2013 were processed for qPCR. Gene data were expressed in number of gene copies ng-1 DNA rather than per gram of soil to minimize any bias related to soil DNA extraction efficiency (Čuhel et al., 2010; Correa-Galeote et al., 2013; Rachid et al., 2013).
A repeated measures two-way analysis of variance (ANOVA) followed by Tukey HSD tests were conducted to detect the effects of plantation type and of sampling time on the measured soil variables and on the abundance of MFG. All data were tested for normality using Shapiro-Wilk normality test and for homogeneity of variance with Levene’s test. Gene abundance data and all soil data except inorganic N were then log-transformed to meet these assumptions. Pearson correlations were performed to analyse the relationships between gene abundances and soil chemical characteristics. A principal component analysis (PCA) was implemented to visualize how MFG abundances were distributed based on plantation type and sampling time. All statistical tests were considered significant at P < 0.05. SPSS version 22 was used for all statistical analyses except for PCA which was carried out in R (FactoMine and FactoExtra packages).
Soil chemical characteristics were significantly influenced by plantation type and sampling time, but the interaction of both factors was not significant, except for soil δ15N (Table 1). Soil pH was higher under teak plantations than under flueggea or mixed plantations, independent of sampling time (Table 2). Conversely, soil total C and N were lower under teak than under mixed plantations, except for total N in samples collected in December 2012. The lowest C:N ratios were also found under teak mono-plantations. The lowest δ15N values were found in soil under flueggea mono-plantations, regardless of the time of sampling. The highest δ15N values were found for teak plantation soil, although they did not significantly differ from mixed plantations. The effect of sampling time on soil pH, total C and total N was significant, and the highest values were generally found in the last sampling event (December 2013). No differences were found between treatments for soil NH4+ (Table 2).
TABLE 1. P-values obtained from a two-way analysis of variance (ANOVA) to detect the effects of plantation type, sampling time, and their interaction on the measured soil variables at Kolombangara, Solomon Islands.
TABLE 2. Soil chemical characteristics measured in different plantation types and at different sampling times at Kolombangara, Solomon Islands.
The abundance of MFG involved in N cycling was significantly influenced by plantation type and sampling time, and the interaction of both factors was generally significant (Table 3). The abundance of total bacteria ranged from 1.27 × 108 to 1.27 × 109 copies ng-1 DNA and was significantly affected by plantation type in the first and last samplings (Figure 2A). The denitrification genes narG and nosZ were generally most abundant in soils of mixed plantations (Figures 2B,E). The abundance of nirS was the highest under teak mono-plantations in the first sampling, and under teak and flueggea mono-plantations in the second sampling. However, in the last sampling, nirS abundance was higher under mixed plantations than under mono-plantations (Figure 2D). The abundance of AOB was higher under teak plantations than under flueggea and mixed plantations in the first and last samplings (Figure 2H). No differences were found between plantation types in nirK, nifH and AOA, except in the second sampling when AOA abundance was the lowest under mixed plantations (Figures 2C,F,G).
TABLE 3. P-values obtained from a two-way analysis of variance (ANOVA) to detect the effects of plantation type, sampling time, and their interaction on the abundances of the functional genes involved in soil N cycling at Kolombangara, Solomon Islands.
FIGURE 2. Abundance of the total bacterial community (A) and of the narG (B), nirK (C), nirS (D), nosZ (E), nifH (F), AOA (G), and AOB (H) genes in soil under plantations of teak (diamonds), flueggea (squares), and mixed species (triangles) at three sampling dates. Gene abundances are expressed in log (gene copy numbers ng-1 DNA). Asterisks indicate significant differences between plantation types (Tukey HSD test, P < 0.05).
Sampling time effects on the abundance of MFG were also detected (Table 3). The PCA biplot showed that sampling time was the main explanatory variable for gene abundance data, rather than plantation type (Figures 3A,B). Generally, the abundances of nirS and nirK genes were higher in the last sampling event (December 2013) than in the first two samplings (Figures 2C,D). On the contrary, the abundances of genes 16S, nosZ and AOA were higher in the first sampling (December 2012) than in the two following sampling events (Figures 2A,E,G), although ANOVA results showed that this difference was not significant for AOA (Supplementary Table S2). The highest abundances for AOB were found in May 2013 (Figure 2H).
FIGURE 3. Principal component analysis (PCA) biplot showing the clustering of microbial functional gene abundances per plantation type (A) and per sampling date (B). T: teak; F: flueggea; M: mixed; Dec12: December 2012; May13: May 2013; Dec13: December 2013.
The relative abundances of narG, nirK, nirS, and nosZ within the total 16S rRNA gene showed that denitrification genes were generally more abundant under mixed plantations, independently from the sampling time (Figures 4A–D). The nosZ/(nirK+nirS) ratio was also larger under mixed plantations than under mono-plantations (Figure 4E). The relative abundances of the nirK and nirS genes revealed a sharp increase in the last sampling while the relative abundance of the nosZ gene significantly decreased in December 2013. As a result, the nosZ/(nirK+nirS) ratio was the lowest in the last sampling (Figure 4E). When expressed in relation to the size of the entire bacterial community, the nifH gene abundance varied depending on plantation type, being larger under mixed plantations than under mono-plantations in the first and last samplings (Figure 4F). The AOA/AOB ratio was also influenced by the sampling date, being the largest in the first sampling. In the first two samplings, the AOA/AOB ratio was larger under flueggea than under teak and mixed plantations. However, no significant differences were detected in the AOA/AOB ratio between plantation types in December 2013 (Figure 4G).
FIGURE 4. Relative abundance of the narG (A), nirK (B), nirS (C), nosZ (D) genes, nosZ/(nirK+nirS) ratio (E), relative abundance of the nifH gene (F), and AOA/AOB ratio (G) in soil under plantations of teak (diamonds), flueggea (squares), and mixed species (triangles) at three sampling dates. Asterisks indicate significant differences between plantation types (Tukey HSD test, P < 0.05).
Pearson correlations showed that soil pH was positively related to the abundance of denitrification genes nirK and nirS, and to that of AOA (Table 4). Soil total C positively correlated with the abundance of total bacteria (16S rRNA), nirK and nirS, but negatively correlated with the abundance of AOA. Total N positively correlated with the abundance of most N-cycling genes except for nosZ and the nitrification genes (AOA and AOB). The C:N ratio strongly and positively correlated with the abundance of nirK and nirS and negatively correlated with the abundance of 16S rRNA and nosZ. Soil δ15N only correlated with the abundance of AOA while NH4+ correlated positively with that of narG and nitrification genes AOA and AOB.
TABLE 4. Pearson coefficients for correlation between functional gene abundances and soil characteristics, Kolombangara, Solomon Islands.
The objective of this study was to determine how the establishment of mixed plantations would influence soil N pools and transformations, measured through soil chemical properties and MFG abundance. Our results show that soil properties and MFG abundances were both influenced by plantation type. Inter-planting teak with flueggea decreased soil pH and may therefore reduce soil bacterial activities and nutrient cycling rates (Rousk et al., 2009), which could in turn lessen N losses from the system (Xu et al., 2013). Accelerations of N transformation rates and increases in N losses through leaching or denitrification have been associated with enriched soil 15N signals (Bai et al., 2015; Reverchon et al., 2015), due to the discrimination against the heavier N isotope during microbially mediated N transformations (Högberg, 1997). Soil δ15N values were significantly higher under teak than under flueggea plantation soils and intercropping teak with flueggea seemed to decrease soil δ15N, although not significantly, which indicates that N losses through volatilization, leaching, or denitrification could be reduced in mixed plantations.
Soil total C and N content increased when teak was inter-planted with flueggea, which translated into higher C:N ratios in soils from mixed plantations than those in teak mono-plantations. Soil C:N ratio is considered a good indicator of soil fertility as it reflects the coupling between soil organic C (SOC) and total N (Lou et al., 2012; Corral-Fernández et al., 2013). Higher C:N ratios may indicate a slowdown of SOC decomposition and N mineralisation as well as an enhancement of N immobilization (Davidson et al., 2003; Puget and Lal, 2005). Mixed plantations could therefore constitute a way to increase soil C and N stocks without losing the source of income that teak trees represent for local communities. Recently, Lang et al. (2014) showed that N acquisition and retention by subtropical trees were enhanced in mixed plantations when compared with those in mono-plantations, and that this increase could be observed as early as the sapling stage. Whether these positive effects will remain once the trees have reached their full height remains to be investigated, although evidences suggest that complementarity effects among co-occurring species, both at the above- and below-ground level, may promote resource-use efficiency in mixed species stands (Cardinale et al., 2007; Richards et al., 2010). Moreover, aging plantations have been shown to have a positive effect on SOC stocks due to an increase in tree productivity and C belowground allocation with time (Eclesia et al., 2012).
The abundance of MFG associated with N cycling was influenced by plantation type, although the sampling time effect seemed to be stronger than that of tree cover. This is consistent with findings by Yamamura et al. (2013) who reported that tree species scarcely influence N cycling genes and bacterial community structure in tropical plantations. The seasonal differences observed in the present study were unlikely due to variations in rainfall or temperature, since climatic conditions in Solomon Islands are relatively homogeneous through the year. However, seasonal variations in soil conditions, affecting plant growth and productivity, bacterial community composition, and ultimately ecosystem processes, could have altered microbial traits. Seasonal differences were detected in pH, total C and total N. Additionally, Rasche et al. (2011) showed that seasonal dynamics of microbial functional groups involved in N cycling were strongly related with seasonal changes in soil labile C and N, and in general with resource availability. This seasonality was especially apparent in the relative abundances of denitrification genes, with the relative abundance of nitrite reductase genes (nirK and nirS) showing a drastic increase in the last sampling. As a consequence, the nosZ/(nirK+nirS) ratio presented a sharp decrease at the same period, as did the AOA/AOB ratio. Denitrifier and nitrifier microorganisms (bacteria and archaea) have been shown to be particularly affected by seasonal shifts in the soil environment, such as soil temperature or moisture (Shen et al., 2008; Stres et al., 2008), which are closely related with resource availability and labile C pools (Bell et al., 2009). Seasonal patterns of tree belowground C allocation may also contribute to temporal variations in microbial dynamics (Kaiser et al., 2011; Churchland and Grayston, 2014). More recently, Graham et al. (2014) reported that temporal dynamics were critical to unraveling the relationships between soil properties and gene abundance data, and thus recommended to examine seasonal samples separately in order to develop accurate models of ecosystem functioning.
Generally, the abundance of MFG associated with denitrification (narG, nosZ, and nirS) was larger in soils from mixed plantations than in those from teak mono-plantations, especially in the last sampling. Gene copy numbers in MFG only vary from 1 to 3 in bacterial cells, while the number of 16S rRNA gene copies range from 1 to 13 per cell (Fogel et al., 1999). Therefore, whilst the observed differences in 16S rRNA gene abundance between tree plantation systems could have been due to differences in species composition of the bacterial community (Liu et al., 2015), our results indicate that denitrifier communities were more abundant under mixed plantations, which is confirmed by the relative abundances of denitrification genes. An increase in the abundance of narG (nitrate reductase) and nosZ (nitrous oxide reductase) may reduce the negative environmental impacts associated with leaching of NO3- and with N2O emissions (Henry et al., 2006). The larger nosZ/(nirK+nirS) ratio under mixed plantations indicated an enrichment in nosZ-harboring denitrifiers when teak was inter-planted with flueggea, which could in turn have implications for N2O/N2 emissions from tree plantations (Liu et al., 2015). The measurement of atypical nosZ abundance, which was not quantified in this study, would also complement the analysis of bacterial and archeal contributions to N2O emissions from plantation soils (Sanford et al., 2012).
Overall, the increase in denitrification gene abundances in mixed plantations may indicate larger denitrifying bacterial communities when teak is inter-planted with flueggea (Baudoin et al., 2009). Denitrifier bacteria have been reported to be influenced by the quantity and composition of organic compounds emitted by roots or resulting from the decomposition of organic residues (Kandeler et al., 2006; Henry et al., 2008). Mixed plantations, by promoting the diversity of bioavailable C sources associated with root exudates and litterfall decomposition, could sustain more abundant denitrifier communities than mono-plantations. The abundance of most denitrification genes correlated with soil total N content and C:N ratio, and nirK/nirS gene abundances positively correlated with total C content. These findings are consistent with results from previous studies (Lindsay et al., 2010; Ducey et al., 2013) and confirm the heterotrophic nature of denitrifier bacterial communities (Baudoin et al., 2009). Furthermore, soil total C, total N, and C:N were all influenced by plantation type, which suggests that tree species mixture is likely to be a significant driver of the denitrifier community through the modification of soil properties.
No differences in nifH abundance were detected between plantation types. This is consistent with findings from Rachid et al. (2013) who showed that soil free-living nitrogen-fixing bacteria were not altered by the establishment of mixed plantations. However, in the first and last samplings, the relative abundance of the nifH gene was larger in mixed plantations than in mono-plantations, which reflects an enrichment in N-fixing bacteria when teak was inter-planted with flueggea and could, over time, increase soil N content (Lammel et al., 2015). Nitrifier communities were also affected by plantation type, as shown by the significant decrease in AOB abundance under mixed plantations compared with teak mono-plantations. This reduction of AOB abundance in mixed plantations did not translate into alterations in soil NH4+ concentrations, although both AOA and AOB genes correlated significantly with soil NH4+. The lack of alteration in soil NH4+ may be due to the dominance of AOA in the ammonia oxidation process, as indicated by AOA:AOB ratios higher than 10 in almost all treatments (Prosser and Nicol, 2012). Furthermore, the abundance and community structure of nitrifiers have been shown to be highly dependent upon pH (Wessén et al., 2011; Levy-Booth et al., 2014), and different phylotypes of bacterial and archaeal ammonia oxidisers are selected in soils with different pH (Nicol et al., 2008), which in turn could affect nitrification rates.
Differences in N pools and transformations between different tree plantation systems may be attributed to a combination of alterations in the soil physico-chemical environment and in soil microbial communities (Anderson et al., 2014). Tree mycorrhizal symbionts could further influence N cycling by accessing different soil N sources or by enhancing the immobilization of inorganic N (Michelsen et al., 1996). While teak has been reported to form arbuscular mycorrhizal fungal associations (Rajan et al., 2000), the mycorrhizal status of flueggea is yet to be studied and may have implications for species interactions in these mixed plantations. Moreover, bacterial taxonomic and functional diversity could be affected by tree species composition directly through changes in the amount and composition of root exudates and litterfall, and indirectly through changes in soil parameters (da C Jesus et al., 2009; Rousk et al., 2010). The significant correlations found in the present study between some soil factors and MFG abundances confirmed the influence of plantation type on soil bacterial communities through the modification of the soil environment. However, there was no strong coupling between gene abundances and the measured soil parameters. This lack of correspondence may be due to different factors. First, process rates may be more readily influenced by changes in tree cover than taxonomic or functional diversity of bacterial communities (Rocca et al., 2015). Moreover, changes in soil parameters may have a stronger effect on the relative proportion of nitrifier/denitrifier bacteria than on MFG abundance (Anderson et al., 2014). Finally, the presence of functional genes does not always indicate an active bacterial community (Levy-Booth et al., 2014) and expression of MFG measured through mRNA transcripts, rather than MFG abundance, could be a better predictor of microbial functional capacity and hence N transformation rates in soil (Wood et al., 2015).
Due to the lack of proper conditions available at the sampling site, it was not possible to freeze the samples and maintain them at -80°C until DNA extraction. Sample storage at 4°C for up to 5 days prior to analysis may have induced changes in bacterial community structure and MFG abundances, as reported by Ott et al. (2004) for fecal samples. However, recent findings by Lauber et al. (2010) and Rubin et al. (2013) showed that neither storage time nor storage temperature drastically affected microbial community composition and structure from soil samples. In addition, storage conditions were consistent across our three sampling events, thereby reducing any bias resulting from sample collection and preparation.
With teak native forests disappearing, teak plantations are expanding and now cover 4.3⋅106 ha (Fernández-Moya et al., 2014). Establishing mixed plantations as alternatives to mono-species teak stands would guarantee economic returns without a subsequent impoverishment of soil fertility. Our results show that inter-planting teak with flueggea would benefit soil quality and sustainability by increasing soil C and N pools, which we hypothesized could be due to a decrease in potential losses from leaching and N2O emissions or to higher inputs from flueggea trees. In addition, soil quality under mixed plantations would be enhanced by the diversity of bioavailable C compounds from root exudates and leaf litter. Kirby and Potvin (2007) reported a positive correlation between biodiversity and organic C in soils, thereby suggesting that mixed plantations could improve soil physical structure (Xing et al., 2004). Our findings are further supported by productivity and aboveground nutrient accumulation data showing the potential for teak to be grown in mixed plantations (Vigulu, 2015). Further research should concentrate on assessing other mixed plantation models under different sets of environmental conditions, and on incorporating microbial dynamics into plantation management to improve productivity while mitigating soil fertility loss.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank the Solomon Islands Government for allowing this research and especially the Ministry of Forestry officers based at Munda for their help with soil sampling. We also thank Peter Whitehead (Kolombangara Forest Products Ltd.) and Vaeno Vigulu for access to the field site, Jupei Shen for her help with molecular analyses, Rene Diocares for his technical assistance with mass spectrometry, and Aldo Segura Cabrera for his help with the PCA analysis. We also acknowledge the constructive comments provided by three anonymous reviewers on a previous draft of the manuscript. Financial support was provided by the Australian Centre for International Agricultural Research (project FST/2012/043).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2015.01439
Figure S1 | Average rainfall data at Kolombangara, Western Province, Solomon Islands, from 1993 to 2014.
Anderson, C. R., Hamonts, K., Clough, T. J., and Condron, L. M. (2014). Biochar does not affect soil N-transformations or microbial community structure under ruminant urine patches but does alter relative proportions of nitrogen cycling bacteria. Agric. Ecosyst. Environ. 191, 63–72. doi: 10.1016/j.agee.2014.02.021
Bai, S. H., Blumfield, T. J., Xu, Z., Chen, C., and Wild, C. H. (2012). Soil organic matter dynamics and nitrogen availability in response to site preparation and management during revegetation in tropical Central Queensland, Australia. J. Soil Sediments 12, 386–395. doi: 10.1007/s11368-011-0466-9
Bai, S. H., Xu, C. Y., Xu, Z., Blumfield, T. J., Zhao, H., Wallace, H., et al. (2015). Soil and foliar nutrient and nitrogen isotope composition (δ15N) at 5 years after poultry litter and green waste biochar amendment in a macadamia orchard. Environ. Sci. Pollut. Res. Int. 22, 3803–3809. doi: 10.1007/s11356-014-3649-2
Balieiro, F. D. C., Pereira, M. G., Alves, B. J. R., Resende, A. S. D., and Franco, A. A. (2008). Soil carbon and nitrogen in pasture soil reforested with Eucalyptus and Guachapele. Rev. Bras. Ciênc. Solo 32, 1253–1260. doi: 10.1590/S0100-06832008000300033
Baudoin, E., Philippot, L., Cheneby, D., Chapuis-Lardy, L., Fromin, N., Bru, D., et al. (2009). Direct seeding mulch-based cropping increases both the activity and the abundance of denitrifier communities in a tropical soil. Soil Biol. Biochem. 41, 1703–1709. doi: 10.1016/j.soilbio.2009.05.015
Bell, C. W., Acosta-Martinez, V., McIntyre, N. E., Cox, S., Tissue, D. T., and Zak, J. C. (2009). Linking microbial community structure and function to seasonal differences in soil moisture and temperature in a Chihuahuan desert grassland. Microb. Ecol. 58, 827–842. doi: 10.1007/s00248-009-9529-5
Bissett, A., Brown, M. V., Siciliano, S. D., and Thrall, P. H. (2013). Microbial community responses to anthropogenically induced environmental change: towards a systems approach. Ecol. Lett. 16, 128–139. doi: 10.1111/ele.12109
Bissett, A., Richardson, A. E., Baker, G., and Thrall, P. H. (2011). Long-term land use effects on soil microbial community structure and function. Appl. Soil Ecol. 51, 66–78. doi: 10.1016/j.apsoil.2011.08.010
Cardinale, B. J., Wright, J. P., Cadotte, M. W., Carroll, I. T., Hector, A., Srivastava, D. S., et al. (2007). Impacts of plant diversity on biomass production increase through time because of species complementarity. Proc. Natl. Acad. Sci. U. S. A. 104, 18123–18128. doi: 10.1073/pnas.0709069104
Chen, C. R., Xu, Z. H., Blumfield, T. J., and Hughes, J. M. (2003). Soil microbial biomass during the early establishment of hoop pine plantation: seasonal variation and impacts of site preparation. For. Ecol. Manage. 186, 213–225. doi: 10.1016/S0378-1127(03)00275-5
Churchland, C., and Grayston, S. J. (2014). Specificity of plant-microbe interactions in the tree mycorrhizosphere biome and consequences for soil C cycling. Front. Microbiol. 5:261. doi: 10.3389/fmicb.2014.00261
Corral-Fernández, R., Parras-Alcántara, L., and Lozano-García, B. (2013). Stratification ratio of soil organic C, N and C: N in Mediterranean evergreen oak woodland with conventional and organic tillage. Agric. Ecosyst. Environ. 164, 252–259. doi: 10.1016/j.agee.2012.11.002
Correa-Galeote, D., Marco, D. E., Tortosa, G., Bru, D., Philippot, L., and Bedmar, E. J. (2013). Spatial distribution of N-cycling microbial communities showed complex patterns in constructed wetland sediments. FEMS Microbiol. Ecol. 83, 340–351. doi: 10.1111/j.1574-6941.2012.01479.x
Čuhel, J., Šimek, M., Laughlin, R. J., Bru, D., Chèneby, D., Watson, C. J., et al. (2010). Insights into the effect of soil pH on N2O and N2 emissions and denitrifier community size and activity. Appl. Environ. Microbiol. 76, 1870–1878. doi: 10.1128/AEM.02484-09
da C Jesus, E., Marsh, T. L., Tiedje, J. M., and de S Moreira, F. M. (2009). Changes in land use alter the structure of bacterial communities in Western Amazon soils. ISME J. 3, 1004–1011. doi: 10.1038/ismej.2009.47
Davidson, E. A., Chorover, J., and Dail, D. B. (2003). A mechanism of abiotic immobilization of nitrate in forest ecosystems: the ferrous wheel hypothesis. Glob. Change Biol. 9, 228–236. doi: 10.1046/j.1365-2486.2003.00592.x
Ducey, T. F., Ippolito, J. A., Cantrell, K. B., Novak, J. M., and Lentz, R. D. (2013). Addition of activated switchgrass biochar to an aridic subsoil increases microbial nitrogen cycling gene abundances. Appl. Soil Ecol. 65, 65–72. doi: 10.1016/j.apsoil.2013.01.006
Eclesia, R. P., Jobbagy, E. G., Jackson, R. B., Biganzoli, F., and Piñeiro, G. (2012). Shifts in soil organic carbon for plantation and pasture establishment in native forests and grasslands of South America. Glob. Change Biol. 18, 3237–3251. doi: 10.1111/j.1365-2486.2012.02761.x
Fernández-Moya, J., Alvarado, A., Forsythe, W., Ramírez, L., Algeet-Abarquero, N., and Marchamalo-Sacristán, M. (2014). Soil erosion under teak (Tectona grandis Lf) plantations: general patterns, assumptions and controversies. Catena 123, 236–242. doi: 10.1016/j.catena.2014.08.010
Fogel, G. B., Collins, C. R., Li, J., and Brunk, C. F. (1999). Prokaryotic genome size and SSU rDNA copy number: estimation of microbial relative abundance from a mixed population. Microb. Ecol. 38, 93–113. doi: 10.1007/s002489900162
Graham, E. B., Wieder, W. R., Leff, J. W., Weintraub, S. R., Townsend, A. R., Cleveland, C. C., et al. (2014). Do we need to understand microbial communities to predict ecosystem function? A comparison of statistical models of nitrogen cycling processes. Soil Biol. Biochem. 68, 279–282. doi: 10.1016/j.soilbio.2013.08.023
Hallin, S., Jones, C. M., Schloter, M., and Philippot, L. (2009). Relationship between N-cycling communities and ecosystem functioning in a 50-year-old fertilization experiment. ISME J. 3, 597–605. doi: 10.1038/ismej.2008.128
He, Y., Xu, Z., Chen, C., Burton, J., Ma, Q., Ge, Y., et al. (2008). Using light fraction and macroaggregate associated organic matters as early indicators for management-induced changes in soil chemical and biological properties in adjacent native and plantation forests of subtropical Australia. Geoderma 147, 116–125. doi: 10.1016/j.geoderma.2008.08.002
Henry, S., Baudouin, E., López-Gutiérrez, J. C., Martin-Laurent, F., Brauman, A., and Philippot, L. (2004). Quantification of denitrifying bacteria in soils by nirK gene targeted real-time PCR. J. Microbiol. Methods 59, 327–335. doi: 10.1016/j.mimet.2004.07.002
Henry, S., Bru, D., Stres, B., Hallet, S., and Philippot, L. (2006). Quantitative detection of the nosZ gene, encoding nitrous oxide reductase, and comparison of the abundances of 16S rRNA, narG, nirK, and nosZ genes in soils. Appl. Environ. Microbiol. 72, 5181–5189. doi: 10.1128/AEM.00231-06
Henry, S., Texier, S., Hallet, S., Bru, D., Dambreville, C., Chèneby, D., et al. (2008). Disentangling the rhizosphere effect on nitrate reducers and denitrifiers: insight into the role of root exudates. Environ. Microbiol. 10, 3082–3092. doi: 10.1111/j.1462-2920.2008.01599.x
Kaiser, C., Fuchslueger, L., Koranda, M., Gorfer, M., Stange, C. F., Kitzler, B., et al. (2011). Plants control the seasonal dynamics of microbial N cycling in a beech forest soil by belowground C allocation. Ecology 92, 1036–1051. doi: 10.1890/10-1011.1
Kandeler, E., Deiglmayr, K., Tscherko, D., Bru, D., and Philippot, L. (2006). Abundance of narG, nirS, nirK, and nosZ genes of denitrifying bacteria during primary successions of a glacier foreland. Appl. Environ. Microbiol. 72, 5957–5962. doi: 10.1128/AEM.00439-06
Kirby, K. R., and Potvin, C. (2007). Variation in carbon storage among tree species: implications for the management of a small-scale carbon sink project. For. Ecol. Manage. 246, 208–221. doi: 10.1016/j.foreco.2007.03.072
Kourtev, P. S., Ehrenfeld, J. G., and Haggblom, M. (2003). Experimental analysis of the effect of exotic and native plant species on the structure and function of soil microbial communities. Soil Biol. Biochem. 35, 895–905. doi: 10.1016/S0038-0717(03)00120-2
Lacombe, S., Bradley, R. L., Hamel, C., and Beaulieu, C. (2009). Do tree-based intercropping systems increase the diversity and stability of soil microbial communities? Agric. Ecosyst. Environ. 131, 25–31. doi: 10.1016/j.agee.2008.08.010
Ladrach, W. (2009). Management of Teak Plantations for Solid Wood Products. International Society of Tropical Foresters News, Special Report. Available at: http://www.istf-bethesda.org/specialreports/teca_teak/teak.pdf (accessed December 2009).
Lammel, D. R., Feigl, B. J., Cerri, C. C., and Nüsslein, K. (2015). Specific microbial gene abundances and soil parameters contribute to C, N, and greenhouse gas process rates after land use change in Southern Amazonian Soils. Front. Microbiol. 6:1057. doi: 10.3389/fmicb.2015.01057
Lang, A. C., Oheimb, G., Scherer-Lorenzen, M., Yang, B., Trogisch, S., Bruelheide, H., et al. (2014). Mixed afforestation of young subtropical trees promotes nitrogen acquisition and retention. J. Appl. Ecol. 51, 224–233. doi: 10.1111/1365-2664.12157
Lauber, C. L., Zhou, N., Gordon, J. I., Knight, R., and Fierer, N. (2010). Effect of storage conditions on the assessment of bacterial community structure in soil and human-associated samples. FEMS Microbiol. Lett. 307, 80–86. doi: 10.1111/j.1574-6968.2010.01965.x
Levy-Booth, D. J., Prescott, C. E., and Grayston, S. J. (2014). Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biol. Biochem. 75, 11–25. doi: 10.1016/j.soilbio.2014.03.021
Levy-Booth, D. J., and Winder, R. S. (2010). Quantification of nitrogen reductase and nitrite reductase genes in soil of thinned and clear-cut Douglas-fir stands by using real-time PCR. Appl. Environ. Microbiol. 76, 7116–7125. doi: 10.1128/AEM.02188-09
Lindsay, E. A., Colloff, M. J., Gibb, N. L., and Wakelin, S. A. (2010). The abundance of microbial functional genes in grassy woodlands is influenced more by soil nutrient enrichment than by recent weed invasion or livestock exclusion. Appl. Environ. Microbiol. 76, 5547–5555. doi: 10.1128/AEM.03054-09
Liu, X., Chen, C., Wang, W., Hughes, J. M., Lewis, T., Hou, E., et al. (2015). Vertical distribution of soil denitrifying communities in a wet sclerophyll forest under long-term repeated burning. Microb. Ecol. 70, 993–1003. doi: 10.1007/s00248-015-0639-y
Liu, X., Chen, C. R., Wang, W. J., Hughes, J. M., Lewis, T., Hou, E. Q., et al. (2013). Soil environmental factors rather than denitrification gene abundance control N2O fluxes in a wet sclerophyll forest with different burning frequency. Soil Biol. Biochem. 57, 292–300. doi: 10.1016/j.soilbio.2012.10.009
Lou, Y., Xu, M., Chen, X., He, X., and Zhao, K. (2012). Stratification of soil organic C, N and C: N ratio as affected by conservation tillage in two maize fields of China. Catena 95, 124–130. doi: 10.1016/j.catena.2012.02.009
Michelsen, A., Schmidt, I. K., Jonasson, S., Quarmby, C., and Sleep, D. (1996). Leaf 15N abundance of subarctic plants provides field evidence that ericoid, ectomycorrhizal and non-and arbuscular mycorrhizal species access different sources of soil nitrogen. Oecologia 105, 53–63. doi: 10.1007/BF00328791
Montagnini, F. (2000). Accumulation in above-ground biomass and soil storage of mineral nutrients in pure and mixed plantations in a humid tropical lowland. For. Ecol. Manage. 134, 257–270. doi: 10.1016/S0378-1127(99)00262-5
Nicol, G. W., Leininger, S., Schleper, C., and Prosser, J. I. (2008). The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. Environ. Microbiol. 10, 2966–2978. doi: 10.1111/j.1462-2920.2008.01701.x
Ott, S. J., Musfeldt, M., Timmis, K. N., Hampe, J., Wenderoth, D. F., and Schreiber, S. (2004). In vitro alterations of intestinal bacterial microbiota in fecal samples during storage. Diagn. Microbiol. Infect. Dis. 50, 237–245. doi: 10.1016/j.diagmicrobio.2004.08.012
Prosser, J. I., and Nicol, G. W. (2012). Archaeal and bacterial ammonia-oxidisers in soil: the quest for niche specialisation and differentiation. Trends Microbiol. 20, 523–531. doi: 10.1016/j.tim.2012.08.001
Rachid, C. T. C. C., Balieiro, F. C., Peixoto, R. S., Pinheiro, Y. A. S., Piccolo, M. C., Chaer, G. M., et al. (2013). Mixed plantations can promote microbial integration and soil nitrate increases with changes in the N cycling genes. Soil Biol. Biochem. 66, 146–153. doi: 10.1016/j.soilbio.2013.07.005
Rajan, S. K., Reddy, B. J. D., and Bagyaraj, D. J. (2000). Screening of arbuscular mycorrhizal fungi for their symbiotic efficiency with Tectona grandis. For. Ecol. Manage. 126, 91–95. doi: 10.1016/S0378-1127(99)00089-4
Rasche, F., Knapp, D., Kaiser, C., Koranda, M., Kitzler, B., Zechmeister-Boltenstern, S., et al. (2011). Seasonality and resource availability control bacterial and archaeal communities in soils of a temperate beech forest. ISME J. 5, 389–402. doi: 10.1038/ismej.2010.138
Reverchon, F., Flicker, R. C., Yang, H., Yan, G., Xu, Z., Chen, C., et al. (2014). Changes in δ15N in a soil–plant system under different biochar feedstocks and application rates. Biol. Fertil. Soils 50, 275–283. doi: 10.1007/s00374-013-0850-2
Reverchon, F., Yang, H., Ho, T. Y., Yan, G., Wang, J., Xu, Z., et al. (2015). A preliminary assessment of the potential of using an acacia—biochar system for spent mine site rehabilitation. Environ. Sci. Pollut. Res. Int. 22, 2138–2144. doi: 10.1007/s11356-014-3451-1
Richards, A. E., Forrester, D. I., Bauhus, J., and Scherer-Lorenzen, M. (2010). The influence of mixed tree plantations on the nutrition of individual species: a review. Tree Physiol. 30, 1192–1208. doi: 10.1093/treephys/tpq035
Rocca, J. D., Hall, E. K., Lennon, J. T., Evans, S. E., Waldrop, M. P., Cotner, J. B., et al. (2015). Relationships between protein-encoding gene abundance and corresponding process are commonly assumed yet rarely observed. ISME J. 9, 1693–1699. doi: 10.1038/ismej.2014.252
Rotthauwe, J., Witzel, K., and Liesack, W. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712.
Rousk, J., Bååth, E., Brookes, P. C., Lauber, C. L., Lozupone, C., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Rousk, J., Brookes, P. C., and Bååth, E. (2009). Contrasting soil pH effects on fungal and bacterial growth suggests functional redundancy in carbon mineralisation. Appl. Environ. Microbiol. 75, 1589–1596. doi: 10.1128/AEM.02775-08
Rubin, B. E., Gibbons, S. M., Kennedy, S., Hampton-Marcell, J., Owens, S., and Gilbert, J. A. (2013). Investigating the impact of storage conditions on microbial community composition in soil samples. PLoS ONE 8:e70460. doi: 10.1371/journal.pone.0070460
Sanford, R. A., Wagner, D. D., Wu, Q., Chee-Sanford, J. C., Thomas, S. H., Cruz-García, C., et al. (2012). Unexpected nondenitrifier nitrous oxide reductase gene diversity and abundance in soils. Proc. Natl. Acad. Sci. U.S.A. 109, 19709–19714. doi: 10.1073/pnas.1211238109
Severin, I., and Stal, L. J. (2010). nifH expression by five groups of phototrophs compared with nitrogenase activity in coastal microbial mats. FEMS Microbiol. Ecol. 73, 55–67. doi: 10.1111/j.1574-6941.2010.00875.x
Shen, J. P., Zhang, L. M., Zhu, Y. G., Zhang, J. B., and He, J. Z. (2008). Abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea communities of an alkaline sandy loam. Environ. Microbiol. 10, 1601–1611. doi: 10.1111/j.1462-2920.2008.01578.x
Stres, B., Danevčič, T., Pal, L., Fuka, M. M., Resman, L., Leskovec, S., et al. (2008). Influence of temperature and soil water content on bacterial, archaeal and denitrifying microbial communities in drained fen grassland soil microcosms. FEMS Microbiol. Ecol. 66, 110–122. doi: 10.1111/j.1574-6941.2008.00555.x
Thevathasan, N. V., and Gordon, A. M. (2004). Ecology of tree intercropping systems in the North temperate region: experiences from southern Ontario, Canada. Agrofor. Syst. 61, 257–268. doi: 10.1023/B:AGFO.0000029003.00933.6d
Wessén, E., Söderström, M., Stenberg, M., Bru, D., Hellman, M., Welsh, A., et al. (2011). Spatial distribution of ammonia-oxidizing bacteria and archaea across a 44-hectare farm related to ecosystem functioning. ISME J. 5, 1213–1225. doi: 10.1038/ismej.2010.206
Wieder, W. R., Cleveland, C. C., Taylor, P. G., Nemergut, D. R., Hinckley, E. L., Philippot, L., et al. (2013). Experimental removal and addition of leaf litter inputs reduces nitrate production and loss in a lowland tropical forest. Biogeochemistry 113, 629–642. doi: 10.1007/s10533-012-9793-1
Wolfe, B. T., Dent, D. H., Deago, J., and Wishnie, M. H. (2015). Forest regeneration under Tectona grandis and Terminalia amazonia plantation stands managed for biodiversity conservation in western Panama. New Forest. 46, 157–165. doi: 10.1007/s11056-014-9448-2
Wood, S. A., Almaraz, M., Bradford, M. A., McGuire, K. L., Naeem, S., Neill, C., et al. (2015). Farm management, not soil microbial diversity, controls nutrient loss from smallholder tropical agriculture. Front. Microbiol. 6:90. doi: 10.3389/fmicb.2015.00090
Xu, Y. B., Xu, Z. H., Cai, Z. C., and Reverchon, F. (2013). Review of denitrification in tropical and subtropical soils of terrestrial ecosystems. J. Soil Sediments 13, 699–710. doi: 10.1007/s11368-013-0650-1
Keywords: Tectona grandis, Flueggea flexuosa, mixed species plantations, δ15N, tropical soil, qPCR
Citation: Reverchon F, Bai SH, Liu X and Blumfield TJ (2015) Tree Plantation Systems Influence Nitrogen Retention and the Abundance of Nitrogen Functional Genes in the Solomon Islands. Front. Microbiol. 6:1439. doi: 10.3389/fmicb.2015.01439
Received: 03 August 2015; Accepted: 03 December 2015;
Published: 22 December 2015.
Edited by:Silvia Pajares Moreno, Universidad Autónoma Metropolitana, Mexico
Reviewed by:Gwen-Aelle Grelet, Landcare Research – Manaaki Whenua, New Zealand
Christopher Ryan Penton, Arizona State University, USA
Rebecca Mueller, Los Alamos National Laboratory, USA
Copyright © 2015 Reverchon, Bai, Liu and Blumfield. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Frédérique Reverchon, firstname.lastname@example.org