Plant n-alkane production from litterfall altered the diversity and community structure of alkane degrading bacteria in litter layer in lowland subtropical rainforest in Taiwan

Some results tend to be very general and need more details, especially when referencing n-alkane degrading bacteria. Very general results in some figures, no details and not everything shown from what is stated in the conclusions. More details need to be presented, as supplementary results? [Reply] We appreciate for the insightful criticizes. More detail information regarding n-alkane bacteria has been presented as supplementary results. The data include the phylogenetic trees of OTUs in Proteobacteria and Actinobacteria together with alkB contained bacteria. A summary table with detail information of these OTUs is also presented.


Introduction
Alkanes are saturated hydrocarbons that account for a small fraction of the total organic carbon, with a concentration range in natural habitats of not more than a couple of thousand parts per million by weight (Rojo, 2010;Gomez-Coca et al., 2016).In non-contaminated habitats, long-chain nalkanes with an odd/even predominance are produced by plants or algae to act as chemoattractants or agents that protect against microbial invasion or water loss (Koch et al., 2009;Feakins et al., 2016).Alkanes are exceptionally stable molecules that can survive in soil and sediment, which make them useful as biomarkers in paleoecology (Sachse et al., 2012).Because alkanes are the major components in petroleum, natural gas, and diesel fuel, the presence of environmental alkanes also serves as an index of oil contamination (Afzal et al., 2013).Although alkanes are very stable, most can be degraded by several microorganisms.Alkanedegrading bacteria are commonly upregulated when the concentration of alkanes is increased; it has thus been suggested that these bacteria be used as a biological indicator for oil and gas prospecting (Rasheed et al., 2012;Xu et al., 2013).
T.-Y.Huang et al.: Litterfall altered alkane-degrading-bacteria structure Extensive knowledge on the degradation of n-alkane by microorganisms has been acquired in the past decades.The first step of enzymatic reaction is oxidization of n-alkane to the corresponding primary alcohol by alkane terminal hydroxylases.This step is crucial because activation of the alkane molecule requires an enzyme system that is not widespread (Rojo, 2010).The products can be further oxidized into a corresponding aldehyde and fatty acid, which are digested by most microorganisms.Although fungi and yeasts can degrade alkanes, most research has focused on the role of bacteria in the degradation of alkanes (Rojo, 2009(Rojo, , 2010;;Singh et al., 2012).To date, alkane monooxygenase gene (alkB) hydroxylases are the most commonly found alkane hydroxylases in both Gram-negative and Grampositive bacteria.More than 60 genera of aerobic bacteria and 5 genera of anaerobic bacteria are known to degrade n-alkanes (Nie et al., 2014).Although alkanes can be metabolized anaerobically, the growth of anaerobic alkane degraders is significantly slower than that of aerobic ones, especially in the litter layer environment where the oxygen supply is adequate (Wentzel et al., 2007).The alkB-encoded membrane-bound homologous protein is an inducible integral membrane-bound alkane hydroxylase that plays a vital role in aerobic alkane degradation (Beilen et al., 2003).Evidence has revealed that the number of alkB bacteria dynamically changes in response to a fluctuation in environmental alkanes.In soil or water environments, when polluted with crude oil, the alkane-degrading bacteria were upregulated (Afzal et al., 2013).For example, a 10-fold increase in the relative abundance of oil degraders in contaminated versus clean sand was estimated for a Gulf of Mexico beach (Kostka et al., 2011).In oil and gas reservoirs, it has been shown that the abundance of alkB bacteria was substantially higher in the surface soils with no significant detection of soil alkane (Xu et al., 2013).The dynamic relationship between alkane and alkB bacteria in soil from an agriculture research farm and a laboratory maintained under controlled conditions has been previously studied.Giebler et al. (2013) reported that in a controlled environment, sizeable communities of alkane-degrading bacteria, especially Proteobacteria and Actinobacteria, were detected at the soil-litter interface when it was enriched with artificial alkane supplements.The researchers demonstrated that many of the alkane-degrading bacteria were not detected before the alkanes were supplemented, suggesting that some alkane degraders may have been present in very low copy numbers and below the limits of detection.It is reported that the substrate was discovered to modulate the biomass of microorganisms based on its initial concentration (Schmidt, 1992).Thus, different alkane yields of litterfall in various natural habitats may upregulate the growth of various alkane degraders.The relationship between alkB bacteria abundance and alkane substrate is more complicated in natural habitats.The dynamic changes in alkB-degrading bacteria populations were demonstrated to be driven by numerous factors, such as the original source of the alkane and soil type.Although knowledge of the regulatory mechanisms regarding environmental plant litter on soil bacterial levels and of the dynamics of alkane-degrading bacteria in natural habitats has been beneficial and important in many aspects, studies on their dynamic changes are currently lacking and should be conducted (Schulz et al., 2012).
The Nanjenshan Reserve, a lowland subtropical rainforest in southern Taiwan, consists of several dynamic forest plots for which forest structure, vegetation pattern, climate, and topography data have been collected at well-established permanent study sites.Four dynamic forest plots and a transect zone in the habitat have been documented since 1989 (Chen and Huang, 1986;Hsieh et al., 2000;Tsui et al., 2004;Fan et al., 2005;Chao et al., 2007Chao et al., , 2008;;Chao et al., 2010).Among them, three plots captured our attention, namely plot I (120 • 50 51 E, 22 • 04 54 N), plot II (120 • 50 36 E, 22 • 04 52 N), and the Lanjenchi plot (120 • 51 38 E, 22 • 03 23 N).A ravine habitat is located across plots I and II, whereas leeward and windward habitats are located on plot III.Surveys conducted over the past decade have demonstrated that the annual yield of litterfall in the ravine habitat was higher than that in the windward and leeward habitats.Moreover, studies have shown that in the leeward habitat, plant height was higher and plant density was lower when compared with the windward habitat (8.41 ± 1.73 vs. 4.63 ± 0.88 m; 7505 vs. 20 065 trees ha −1 ).The plant height and density in the ravine habitat were 9.45 ± 1.35 m and 4257 trees ha −1 (Chin, 2008).In addition, the litterfall plays an essential role in sustaining the microorganism food chain in the forest ecosystem.The litter-leaf, the main component of litterfall, similarly provides the major source of plantderived alkanes.The significant effects of litterfall derived alkanes on the bacterial alkane degradation potential have also been observed (Giebler et al., 2013 andSchulz et al., 2012).As the litterfall in the habitats of the Nanjenshan Reserve was different, investigation into the effects of litterfall on various habitats and microbial communities in these habitats is possible.
Environmental microbial communities are highly complex and diverse and difficult to culture under laboratory conditions.Next-generation sequencing (NGS) is a powerful tool for high throughput sequencing technologies.With this recent advent platform and computational methods, we can conduct genome studies of microbes in these habitats, which have numerous variables such as annual litterfall productivity and plant vegetation (Degnan and Ochman, 2012).We investigated the productivity of long-chain n-alkane in the pristine natural habitats of the Nanjenshan Reserve and explored their relationship with microbial communities, particularly emphasizing the correlation between n-alkane production and bacteria containing the alkB gene.During the present 1-year follow-up study, gas chromatography with a flame ionization detector (GC-FID) was used to determine n-alkane (C 14 -C 35 ) concentrations.Bioinformatics analyses and community amplicon sequencing data mining were per-formed to identify microbial communities and read numbers of bacterial lineages that contained alkB genes in the three rainforest habitats.Semi-quantitative polymerase chain reaction (PCR) was performed to test the results of the bioinformatics analysis.
2 Materials and methods

Sample description
A map of the sampling sites is provided in Fig. S1 in the Supplement.Twenty-two custom-made aluminum stands with frames (0.71 m × 0.71 m) covered by a nylon mesh (1 mm i.d.) were used to collect litterfall.Twelve stands were located in the ravine habitat, whereas 10 stands were located in the Lanjenchi plot.Of those on the Lanjenchi plot, six stands were placed in the windward habitat and four in the leeward habitat.The horizontal distance between stands was approximately 50 m.The number of stands in the three habitats of the Nanjenshan Reserve was mainly selected based on the stands available in this area.Each nylon mesh containing litterfall was collected and replaced monthly from early October 2012 to late September 2013.After their collection, the litterfall samples were sorted into litter-leaf, litter-branches, litter-flowers, litter-fruits, and miscellaneous fractions.The sorted samples were weighed after being oven-dried at 40 • C for 14 days (Conklin-Brittain et al., 2006).The litter-leaves were further separated into genus and species.Samples of the litter layer and underlying soils from the three habitats were collected in 2015 and 2016, respectively.The collected samples, including litter-leaf, litter layer, and soils, were subjected to n-alkane analysis, and litter layer samples from each habitat were further subjected to DNA extraction and consequential assays.
The annual yield (litterfall or litter-leaf) in each stand was used to calculate the productivity of each habitat.The annual litterfall and litter-leaf yields in each habitat were calculated by averaging the annual stand weights.Statistically significant differences between habitats were evaluated through Student's t test using the t test calculator (SD format).

Extraction of n-alkanes
All samples from the litter-leaf (42 plant species) and litter layer and all soils were dried and well homogenized.Approximately 0.5 g of sample powder was dissolved in dichloromethane (250 mL) and Soxhlet extracted for 16 h.A surrogate standard of 5-α-cholestane (Sigma) was added to each sample before extraction.The extracts were concentrated into 5-10 mL volumes using a rotary evaporator with a water bath and further concentrated to reduce the volume to approximately 1 mL, after which they were fractionated using silica gel chromatography (Silica gel 60, 3 %).The n-alkane fraction was eluted with 4 mL of hexane and reduced to a 1 mL volume under a stream of nitrogen.An in-ternal standard of squalene (Acros) was added to the concentrated 1 mL of hexane fraction prior to instrumental GC analysis.Quantification of the n-alkanes was performed using a PerkinElmer Clarus 500 gas chromatograph equipped with an autosampler, PerkinElmer Elite-5 CB fused silica capillary column (30 m in length, 0.32 mm i.d., film thickness 0.25 µm), and flame ionization detector (FID).The GC oven temperature was programmed to increase from 70 to 310 • C at a rate of 4 • C min −1 .The n-alkanes between C 14 and C 35 were identified and quantified using external standards that contained known concentrations of all n-alkanes of interest (Dr.Ehrenstorfer GmbH), 5-α-cholestane, and squalene.
The n-alkane concentrations in 42 litter-leaf samples were measured and then used to estimate the n-alkane concentration from the litter-leaf and to assess litter-leaf productivity in the three habitats.The n-alkane concentration in the litter-leaf from each stand (C stand , in µg g −1 ) was summed by weighting the concentration in one species in the litter-leaf by that species' proportion of the total litter-leaf collected in the stand as in the following equation: where C stand is n-alkane concentration in the litter-leaf from each stand (µg g −1 ), W i is the dry weight of litter-leaf of each species (g), W t is the dry weight of total litter-leaf (g), and C i is the n-alkane concentrations (µg g −1 ) of 42 litter-leaf samples.
Calculation of the litter-leaf n-alkane concentration in each habitat was then performed by averaging the n-alkane litter-leaf concentrations in all the stands.The annual nalkane productivity in each stand was calculated by multiplying the n-alkane litter-leaf concentration in each stand by the annual yield of that stand.The n-alkane litter-leaf productivities of the stands in a habitat were then averaged to estimate the annual n-alkane productivity in the three habitats, which was evaluated using Student's t test.

DNA extraction
For DNA extraction in the litter layer, several pieces of leaves from the litter layer were randomly selected, cut into 2 mm × 2mm chips, and well blended.The DNA extraction procedure is described in detail in the user manual, Genomic DNA from soil (Macherey-Nagel).Briefly, a bulk mixture (approximately 0.1-0.2g) was transferred to a tube containing ceramic beads.Lysis buffer (SL1, 700 µL) and enhancer SX (150 µL) were added to each sample and vortexed for 5 min at room temperature.Precipitate contaminants were removed through centrifugation for 2 min at 11 000 g. DNA supernatant was transferred to a new tube, 150 µL of SL3 was added, and the mixture was vortexed for 5 s before precipitation by centrifugation for 1 min at 11 000 g.The supernatant inhibitors were removed using an inhibitor removal www.biogeosciences.net/15/1815/2018/Biogeosciences, 15, 1815-1826, 2018 column.The column that bound the DNA was washed four times and dried before DNA elution.DNA elution was performed by centrifugation for 30 s at 11 000 g after incubation with Buffer SE (50 µL).A final DNA extraction (approximately 30-100 ng µL −1 in 100 µL of Tris buffer) was eluted.

The construction and analysis of community amplicon sequencing library
FASTA files were further processed using the sequence analysis tool, USEARCH.We merged all sequence files, removed duplicates, and clustered the sequences into operational taxonomic units (OTUs) at 97 % pairwise identity with a minimum cluster size of two to construct an OTU reference library.Then, a comparison between samples and the ref-erence library at a 97 % sequence identity level was performed to create OTU tables, which contained the number of DNA sequence reads of each OTU.A 16S UTAX reference database was employed as a BLAST (Basic Local Alignment Search Tool) library.Finally, these data were combined to determine the relative abundances and subsequent visual patterns of a heat map.Principal coordinate analysis (PCoA) was used to identify the relative distance between OTUs.To plot PCoA figures, the calc_distmx command together with a phylip_lower_triangular parameter was used to calculate the distance of the OTUs.The axes file contained plotting coordinates and was derived from mothur using the PCoA command.The PCoA figures were plotted using the R program.
2.6 Proportion of alkB-gene-containing bacteria in the habitats estimated using the NCBI database Data mining for community amplicon sequencing was performed to search for microbial communities that may have had alkB genes.Because bacteria in the same genus may have similar gene sequences and functions, nucleotide blasting was performed to search for homogeneous sequences in the OTUs.Bacteria with DNA sequences 95 % or more similar to a known gene were defined as lineage-like bacteria.The NCBI database was used to construct an alkB gene reference library.We collected representative 16s gene sequences of all known alkB-gene-containing bacteria in the phyla Actinobacteria and Proteobacteria from the current NCBI database to investigate the alkB gene family populations in different habitats.The alkB gene sequences in bacterial 16s rRNA were downloaded to a reference nucleotide file and combined with the OTU libraries of Actinobacteria and Proteobacteria.
To search the alkB family gene-like lineage in the current OTU database, sequence alignment and phylogenetic analysis of the alkB nucleotide sequences were conducted using the Molecular Evolutionary Genetics Analysis 7 (MEGA 7) program.MEGA 7 used parsimony, neighbor-joining, and maximum likelihood analyses to create a 16S rRNA phylogenetic gene tree.The NCBI nucleotide BLAST program was employed to test the similarity of DNA sequences between OTUs and adjacent references.DNA sequences with a similarity of more than 95 % were manually selected to determine the alkB-lineage OTU DNA sequence numbers read in the library.The Shannon-Wiener diversity index was used to calculate the plant diversity, overall bacteria communities, and alkB-lineage bacteria and was determined using the PRIMER-5 software (Plymouth Routines in Multivariate Ecological Research).

Semi-quantitative PCR for determining alkB gene levels in litter layer samples
Because alkB genes are highly abundant and diverse, a standard method for alkB gene quantification was not available.
We employed semi-quantitative PCR for measuring environmental alkB in a method that was modified from previous studies because the alkB-targeting primer covers Proteobacteria and Actinobacteria (Kloos et al., 2006;Jurelevicius et al., 2013).Litter layer samples from three habitats in Nanjenshan were subjected to semi-quantitative PCR to quantify the DNA levels of alkB genes in the natural habitats.Litter layer samples from different stands in each habitat were randomly selected for sample collection at three time points.After collection, litter layer samples were sent immediately to the laboratory at room temperature.Litter layer patches were randomly selected and cut into approximately 2 mm × 2 mm chips.Samples (250 mg) were taken for DNA extraction.Following DNA elution, 100 µL of buffer solution and 2 µL of eluted DNA were taken for semi-quantitative PCR.The DNA was mixed with an alkB gene primer set (forward primer: 5'-AAY ACN GCN CAY GAR CTN GGN CAY AA -3'; reverse primer: 5'-GCR TGR TGR TCN GAR TGN CGY TG -3', 1 µL, 0.4 µM), 5 µL Fast-Run Taq Master Mix with Dye, and 16 µL ddH 2 O with the final reaction volume of 25 µL.Thermal cycling conditions for PCR were an initial temperature of 95 • C for 5 min, followed by 30, 27, or 24 cycles of 95 • C for 30 s, 55 • C for 30 s, 72 • C for 1 min, and finally 72 • C for 5 min.The expected amplicons of the PCR were 548 bps.Following PCR, aliquots (5 µL) of each amplicon were confirmed using electrophoresis on 1.5 % agarose gel followed by 0.5 µg mL −1 SYBR Safe DNA staining and imaging using an ultraviolet transillumination system.The average band density levels were compared with a windward habitat density setting of 100 %.Statistically significant differences were evaluated using Student's t test.
3 Results and discussion

Quantitative estimates of annual litterfall in the Nanjenshan Reserve
The annual average litterfall was estimated from the litterfall at the sampling sites.The annual yields of litterfall in the three habitats in Nanjenshan are illustrated in Fig. 1.As shown in the records from the past decade, the litterfall was higher in the ravine habitat than in the other habitats.Our results on the litterfall of Nanjenshan indicated that the annual weight productivity in the ravine habitat was approximately 30 % higher than in the leeward or windward habitats, which is consistent with previous studies.According to records from 1999 to 2007, seasonal litterfall in Nanjenshan ranged from 2 to 7 t ha −1 ; the increases in litterfall were contributed by typhoons, precipitation, and the effects of monsoon seasons.In this study, the annual litterfall was discov- Annual productivity of litterfall in ravine habitat was higher than windward and leeward habitats (* p < 0.05).
ered to be approximately 7-10 t ha −1 , which is considered to be a high volume when compared with the records.The primary reason for the large amount of litterfall during this study was the occurrence of typhoons.Because typhoons increase litterfall in the three habitats equally, the ranking in annual litterfall volume was not altered by the advent of typhoons.One of the reasons that the ravine habitat had more annual litterfall may be the topography of the ravine habitat, which is more suitable for plant growth than that of the other habitats.The average height and diameter of cross sections of plants in the ravine habitat were the largest of those in all the habitats.Because litter-leaf was the major part of litterfall and average n-alkane resources, we focused on n-alkane levels in the litter-leaf in these habitats.

Example of n-alkane measurements in the Nanjenshan Reserve
The n-alkane concentrations in litter-leaf from 42 plants were measured using GC-FID.Figure S2  terns.This paper is the first to reveal the n-alkane levels of plant vegetation in Nanjenshan.

Annual litter-leaf productivity and n-alkane productivity rankings
A total of 42 plant species were identified in this study.Data from the plant species revealed the total amount of leaves in the annual litterfall, and the resultant ranking of leaf nalkane concentrations in these plants is presented in Table 1.The leaf n-alkane concentration and dry weight data were employed to make a gross estimate of the n-alkane productivity of the litter-leaf in the Nanjenshan Reserve.The nalkane concentration was discovered to change during leaf development and was affected by environmental parameters (Jetter and Schäffer, 2001;Kahmen et al., 2011;Hoffmann et al., 2013).For example, temperature and relative humidity affected the composition of n-alkanes in both Acacia and Eucalyptus in Australia (Hoffmann et al., 2013).In our pilot study, the n-alkane concentration variation among plant species was less than 15 % of the standard error, as determined through comparison of several samples from different stands and habitats for the same species.The ranges of leaf n-alkane concentrations between species were from 2694 (µg g −1 ) in Ilex rotunda to 123 (µg g −1 ) in Alniphyllum pterospermum.These data demonstrated that intraspecies variation in n−alkane concentration could be ignored compared with the interspecies variation.Therefore, using a representative plant n-alkane concentration was reasonable for estimating the overall litter-leaf n-alkane production in a habitat.When comparing the plant vegetation and leaf nalkane concentration, the plants with leaves that had a high n-alkane concentration tended to grow in the ravine habitat; Ilex rotunda and Celastrus kusanoi, for example.

Estimates of annual n-alkane productivity in the three habitats
Figure 2 presents the n-alkane concentrations from the litterleaf to the bulk soil and the average litter-leaf n-alkane productivity in the three habitats.The average litter-leaf nalkane concentration was highest in the ravine habitat.Moreover, the n-alkane concentration in different layers formed a steep gradient from the litterfall to the bulk soil.The n-alkane flux in the litterfall in the ravine habitat was approximately twice that in the other habitats.We concluded that the ravine habitat had higher n-alkane input than the other habitats.This characteristic offers us an example of higher n-alkane input in the natural habitat to investigate the consequence of microbial community.
The balance between litterfall yield and decomposition affects the development of organic carbon in soil layers.This study and previous reports have indicated that the litterfall in ravine habitats is higher than that in other habitats.However, a study of these habitats demonstrated that the total organic carbon in the litter layer and bulk soil in the ravine habitat was equal to or less than that in the windward and leeward habitats (Guo, 2010).Therefore, decomposition effects may be stronger in ravine habitats.Figure 2a reveals that the nalkane concentration decreased significantly from litter-leaf to surface soil, which could be due to the effects of bacteria.It is plausible that other organic compounds were degraded in addition to the n-alkanes.Several effects such as leaching, erosion, and sediment transportation could play roles in governing organic compound accumulation; therefore, we cannot rule out the diffusion effects of other physical factors (Quinton et al., 2006;Kirkels et al., 2014).We performed NGS-based analysis and bioinformatics studies to discover the roles of microbial communities in n-alkane degradation.

Bacterial community composition in the Nanjenshan Reserve
NGS-based study was performed to reveal the microbial communities in the litter layer from the three habitats.Figure 3 shows the relative abundance in community amplicon sequencing data of OTUs that were grouped by phylum.
More detailed information about the relative abundance in class, order, and family levels is presented in Fig. S5.Interestingly, a similarity was identified in microbial community structure between the windward and leeward habitats.Because the windward and leeward habitats are located on the same plot of a different orography, several parameters such as plant vegetation and soil properties were similar between the habitats.As mentioned previously, the most considerable differences between the windward and leeward habitats were the average height of plants and average number of trees per hectare.The physical parameters that affect plant growth in mountainous areas did not appear to affect the relative abundance of microbial communities.Conversely, the relative abundance of microbial communities in the ravine habitat was different from that in the windward and leeward habitats, with more members of the Actinobacteria phyla and fewer members of the Proteobacteria and Acidobacteria phyla.The more abundant microbial communities of the Actinobacteria phyla may have been because of higher n-alkane productivity in the ravine habitat.Although many phyla have been identified as carrying the alkB gene, our data revealed that the only microbial communities carrying the alkB gene in our study were Proteobacteria and Actinobacteria.

Prediction of bacteria carrying alkB genes in Nanjenshan Reserve
The relative abundance of bacteria containing the alkB gene in the Actinobacteria and Proteobacteria phyla from different habitats is illustrated in Fig. 4. The phylogenetic trees of OTUs in the Proteobacteria and Actinobacteria phyla with alkB reference strains are shown in Figs.S3 and S4.The nucleotide sequence similarities of OTUs and the adjacent alkB Table 1.The resultant ranking of n-alkane concentration of 42 litter-leaf samples with the nomenclature and the semi-quantitative results of annual litter-leaf productivity in three habitats.The plus sign (+) was used to denote the annual productivity rankings of litter-leaf in each habitat.Five plus signs (+++++) denote the litter-leaf productivity rankings from 1 to 5; four plus signs (++++) denote the productivity rankings from 6 to 10; three plus signs (+ + +) denote the rankings from 11 to 15; two plus signs (++) denote the rankings from 16 to 20; one plus sign (+) denotes the rankings after 21.A blank space denotes no litter-leaf yield for the entire year in the habitat.reference strains were manually blasted.Table S1 in the Supplement reveals the nomenclature of OTUs in Proteobacteria and Actinobacteria with nucleotide sequences that had similarity of higher than 95 % when compared with adjacent alkB reference strains.Figure 4 illustrates that a significant number of reads in the OTUs whose sequences were similar to those of alkB-gene-containing bacteria was identified for the ravine habitat.Although the relative abundance of Proteobac- teria in the ravine habitat was 20 % less than that in the other habitats (Fig. 3), the number of DNA sequence reads of alkBlineage-like bacteria in the ravine habitat was greater than 1.5 times those in the windward and leeward areas (Fig. 4).To summarize, there were twice as many DNA sequence reads of alkB-lineage-like bacteria in the ravine habitat than in the other habitats.Despite the relative abundance of DNA sequence reads in many OTUs of the alkB-lineage-like bacteria, Proteobacteria, and Actinobacteria phyla were limited in the leeward and windward habitats (Fig. 4); these two habitats had small amounts of alkB-lineage-like bacteria (approximately 0.1 % relative abundance) in most OTUs (Table S1).Our data provide evidence that there is a universal presence of alkane degraders in natural habitats.

PCoA of bacteria carrying alkB genes
PCoA was used to visualize the similarities of DNA sequences in Proteobacteria and Actinobacteria from the three habitats.The biodiversity index in OTUs of alkB-lineagelike bacteria and other sample types is shown in Table S2.
Figure 5  alkB-lineage-like bacteria OTUs and relative locations are presented in the figure.Thirty-three OTUs contained alkBlineage-like bacteria in this study.The numbers of OTUs that contained alkB-lineage-like bacteria in the ravine, windward, and leeward habitats were 30, 30, and 25, respectively.The Shannon-Wiener index of OTUs that contained alkBlineage-like bacteria in the ravine, windward, and leeward habitats was 2.51, 2.39, and 2.46, respectively.Although the relative abundance of OTUs that contained alkB-lineage-like bacteria in the ravine habitat was more than twice that in the windward and leeward habitats (Fig. 4), the Shannon-Wiener indices of the three habitats were similar, suggesting that the effective species numbers of alkB-lineage-like bacteria were almost identical.The results indicate the preexistence of diverse alkB-gene-containing bacteria in natural habitats, with some bacteria so small in number that they are undetectable using normal PCR.They could, however, proliferate in the habitats with the appropriate substrate supplements.This finding is consistent with the results from a previous study using a laboratory-controlled system and an agricultural research farm (Schulz et al., 2012;Giebler et al., 2013).

Empirical testing of bacteria carrying alkB gene
In the bacterial community amplicon sequencing study, equal amounts of DNA were amplified to estimate the bacterial communities in each habitat.To compare the productivity of the alkB-gene-containing bacteria in different habitats, we conducted semi-quantitative PCR to determine relative alkB gene concentrations.Figure 6a shows a triplicate study of DNA staining of amplicons from the three habitats after various semi-quantitative PCR cycles.No bands were detected if fewer than 24 amplification cycles were performed.Marginal PCR products could be detected by increasing the cycle number to 27. Figure 6b shows the statistical results regarding alkB gene levels in the three Nanjenshan habitats and confirms that bacterial alkB numbers were higher in the ravine habitat than in the other two habitats.We have thus provided evidence that litterfall and plant vegetation increase the abundance of alkane-degrading bacteria communities in natural habitats.

Exploring the relationship between plant vegetation and microbial communities using heat map analysis
Heat map analysis of plant species and read number of OTU versus habitat were performed to investigate the relationship between plant vegetation patterns and microbial communities in different habitats.The heat map presented in Fig. 7 was based on the microbial community abundance and plant productivity of litter-leaf from different habitats, with the color intensity indicating the number of reads or weight of litter-leaf (in logarithmic form).In Figure 7b, the color patterns of litter-leaf yields show that microbial abundance was similar in the windward and leeward habitats, demonstrating the similarity between the microbial community compositions in these two habitats.Furthermore, the litterfall data were also similar between the two habitats, which suggests a connection between plant vegetation species and microbial composition.The pattern similarities in the plant and bacterial species can be easily explained by substrate-specificinduced growth responses in the natural habitats.The data provided supporting evidence that diversity of bacteria in all OTUs was affected by plant vegetation, not only in quantity but also in species.In Figs. 2, 4, and 6, an increase in nalkane productivity is demonstrated to increase the number of alkB-gene-containing bacteria.In Fig. 7, we showed that a substrate-specific relationship may exist between plant vegetation and microbial communities.Thus, a leaf in a natural habitat may upregulate some alkB-gene-containing bacteria while downregulating other types of alkB-gene-containing bacteria.A future study could include a litter bag study that would enable better understanding of the relationship between individual factors within plant and bacterial species; for instance, one could study the interaction of a specific leaf species with a particular genus of n-alkane-degrading bacteria or microbial community in different habitats.

Conclusions
We confirmed that the annual litterfall productivity of the ravine habitat was approximately 30 % higher than that of the windward and leeward habitats, whereas the long-chain n-alkane productivity in the ravine habitat was twice that in the other habitats in a lowland subtropical rainforest of Taiwan.Plant vegetation not only affects n-alkane productivity; it also changes litter layer microbial communities.In the ravine habitat, high litter-leaf n-alkane productivity resulted in an increase in the n-alkane-degrading bacteria in the litterfall.Increases in alkB-gene-containing bacteria were identified primarily in the Proteobacteria and Actinobacteria phyla, with NGS-based analysis and semi-quantitative PCR employed for determining alkB gene concentrations.Biodiversity analyses showed that the diversity of bacteria containing the alkB gene was similar among the different habitats, suggesting that natural habitats harbor diverse alkBlineage-like bacteria with some species existing in very low copy numbers.Practical testing of alkB gene concentration by semi-quantitative PCR with primer coverage over Proteobacteria and Actinobacteria shows results consistent with the bioinformatics data, suggesting this method is helpful for future study on alkane bacteria activity assays from the environmental samples.Cross analyses regarding the quantity and quality of plant vegetation and structures of the microbial communities revealed that n-alkane-induced increases in alkB-gene-containing bacteria might have been substrate specific.To our knowledge, this is the first report on the regulation of n-alkane-degrading bacteria in litterfall in natural habitats.

Figure 1 .
Figure 1.Annual litterfall in three habitats of Nanjenshan Reserve.Annual productivity of litterfall in ravine habitat was higher than windward and leeward habitats (* p < 0.05).

Figure 3 .
Figure 2. (a) n-Alkane concentration in litter-leaf, litter layer, and surface soil.n-Alkane concentration in litter-leaf was higher in the ravine habitat than in the other two habitats (* p < 0.05).(b) Estimated annual n-alkane productivity generated by litter-leaf of litterfall in three habitats (* p < 0.05).

Figure 4 .
Figure 4.The relative abundance in microbial community of alkB gene lineages.White block indicated the relative abundance of alkB gene in phylum Actinobacteria while the striped block was in Proteobacteria.

Figure 5 .
Figure 5.The PCoA plot of OTUs data in phylum Actinobacteria and Proteobacteria in three habitats.The circle areas in A, B, C, and D are classes of α-Proteobacteria, β-Proteobacteria, γ -Proteobacteria, and Actinobacteria, respectively.

Figure 6 .
Figure 6.Semi-quantitative PCR of alkB gene results from samples of three habitats in Nanjenshan Reserve.(a) Agarose electrophoresis of alkB genes in different amplified cycles.Standard markers were loaded at both sides of samples.The samples from different habitats were separated by dash line.(b) The statistic results of semi-quantitative PCR after amplifying 30 cycles.The read numbers of alkB gene were significantly higher in the ravine habitat than in the other two habitats (* p < 0.05).

Figure 7 .
Figure 7. Heat map analysis of read number of OTUs and annual litterfall productions of each plant species in three habitats.(a) Heat map of read number of OTUs in three habitats.(b) Heat map of litterfall productions of plants in three habitats.