• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    In situ Assemblies of Bacteria and Nutrient Dynamics in Response to an Ecosystem Engineer, Marine Clam Scapharca subcrenata, in the Sediment of an Aquaculture Bioremediation System

    2020-11-30 04:43:14NICHOLAUSReganLUKWAMBEBetinaYANGWenZHUJinyongandZHENGZhongming
    Journal of Ocean University of China 2020年6期

    NICHOLAUS Regan, LUKWAMBE Betina, YANG Wen, ZHU Jinyong, and ZHENG Zhongming, *

    Assemblies of Bacteria and Nutrient Dynamics in Response to an Ecosystem Engineer, Marine Clam, in the Sediment of an Aquaculture Bioremediation System

    NICHOLAUS Regan1), 2), LUKWAMBE Betina1), 3), YANG Wen1), ZHU Jinyong1), and ZHENG Zhongming1), *

    1)School of Marine Sciences, Ningbo University, Ningbo 315832, China 2) Department of Natural Sciences, Mbeya University of Science and Technology, Mbeya 53000, Tanzania 3)Department of Food Science and Technology, University of Dar es Salaam, Dar es Salaam 35091, Tanzania

    The relationships between nutrient dynamics, microbial community, and macrofauna are important in bioremediation systems. In this study, we examined the effects of marine clamsedimentary activities on the microbial assemblages, benthic nutrient chemistry, and its subsequent remediation impacts on organic effluent in the sediment of an aquaculture wastewater bioremediation system. The results showed that microbial community composition differed significantly in the clam area (ANOSIM,=0.707,=0.037). Pyrosequencing of bacterial 16S rRNA gene revealed a total of 48 unique phyla, 79 classes, 107 orders, 197 families, and 321 genera amongst all samples. The most dominant bacterial assemblages were Proteobacteria, Bacteroidetes, Acidobacteria, Firmicutes, Verrucomicrobia, and Actinobacteria, with Bacteroidetes and Firmicutes significantly higher in all treatment samples than control (<0.001). All dominant phyla in the list were shared across all samples and accounted for 89% (control) and 97% (treatment) of the total 16S rRNA. The nutrient flux rates from the sediments into the water (treatment group) were 51% (ammonium), 88% (nitrate), 77% (nitrite) and 45% (phosphate) higher, relative to the control implying increased mineralization, degradability, and mobility of the benthic nutrients. Similarly, significantly increased oxygen consumption rates were evident in the clam area signifying improved oxygen distribution within the sediment. The organic effluent contents associated with total- organic matter, carbon, nitrogen, and, phosphate were lower among the clam treatments relative to the control. Our results describe the potential roles and mechanisms contributed by marine bivalveon benthic-bacterial-community assembly, nutrient balance, and effluent reduction in the sediments of aquaculture wastewaters bioremediation system.

    ; effluents; bioturbation: sediment chemistry; bioremediation; nutrient dynamics

    1 Introduction

    Intensive aquaculture has a significant adverse environmental impact due to the release of great amounts of organic wastes (Karakassis., 2000). For instance, wastewater generated from aquaculture productions is normally characterized by high chemical effluents associated with nitrogen wastes (ammonia, nitrate, and nitrite; Dauda., 2014) as well as uneaten feeds and excretion products. These products usually pollute the ecosystem of the receiving water bodies causing severe alterations of the physical-chemical properties and eventually accumulates in the sediments (Diaz, 2001). This understanding has introduced many apprehensions about the environmental impacts of aquaculture wastes. One of the key challenges for the sustainable development of aquaculture is the reduction of its environmental impacts. To maintain sustainable aquaculture production and its subsequent environmental pollution, various effluent remediation methods have been established including integrated aquaculture effluent bioremediation systems (Brito., 2018; Lukwambe., 2018).

    Integrated aquaculture effluent bioremediation systems involving macrofauna such as bivalves is one of the most cost-effective and rising biological wastewater treatment technology (Lukwambe., 2018; Nicholaus., 2019b). The effectiveness of this system largely depends on the interaction between the sediment microbial community, and the benthic macrofauna activities (., clam- burrowing) (Lukwambe., 2018; Nicholaus., 2019b). The practice of incorporating clam in aquaculture restoration systems to mitigate for the excess nutrients is common and predicated on the idea that clam will optimize and/or enhance the microbial growth and activities, and nutrient cycling in the system. Clams are the most important component in the system. They could not only reduce the nutrients by feeding on organic detritus but also remediate the sediment during the bioturbation process (., Nicholaus., 2019a, b). For example, the organic carbon produced by fish farming in the form of uneaten feed and excretions eventually deposit in the sediment and represent a great source of food for filter-feeding organisms, such as bivalves (Navarrete-Mier., 2010). Similarly, clams among other macrofauna inhabiting marine sediments markedly alter the physical structure and chemical composition and microbial biomass turnover (., Sinkko., 2011; Shen., 2017). In this case, clams played concurrent activities (feeding, bioirrigation, digging, mineralization, and burrowing) which are commonly referred to as ‘bioturbation’, a process of sediment reworkings by benthic animals (Aller, 2001). This process offers a significant ecological role in bioremediation practices (Nicholaus., 2019b).

    In the integrated model, sediments are responsible for transporting and cycling of a significant proportion of many nutrients and contaminants (Zhao., 2019). It is estimated that up to 80% of the nitrogen needed by the primary producers in the aquatic system is provided by sediment remineralization reactions (Dale and Prego, 2002). Conversely, marine sediments represent an important dynamic of inorganic nutrientsanaerobic decompositions regulated by microorganism under influence of bioturbating organisms (Lukwambe., 2018). Similarly, microorganisms are very useful and with major importance in aquaculture industrial wastewater treatment systems. They dictate several positive effects on the fate of aquaculture operations (Moriarty, 1997). The greatly microbial roles include the elimination of toxic environmental materials (ammonia, nitrite, and hydrogen sulfide), degradation of uneaten feed, and management of nutrient cycling and water quality (Moriarty, 1997; Yang., 2018). These among other functions make the interaction between clams and microorganisms key players in aquaculture bioremediation systems.

    Various studies including the commonly known ecosystem engineers like filter-feeder bivalves (Mussels: Zhang., 2011; Razor clam: Lukwambe., 2018; Zhao., 2019; blood cockles: Nicholaus., 2019b) have been previously reported to significantly influence nutrient dynamics and improve bioremediation environment. Moreover, to improve the previous studies, and collect sufficient reports that describe the unique roles and mechanisms played by intertidal shall burrower marine bivalves such as(, Zhang., 2020) in marine ecological sediment-based bioremediation system is necessary. Therefore, the present study aimed to 1) investigate and characterize the distribution and assembly of the microbial community in a field ecological bivalve- sediment based bioremediation system, 2) examine the ecosystem bioengineering mechanisms ofin boosting sediment microorganism, and 3) examine its potential effects on benthic biogeochemical nutrient fluxes and organic wastes mitigation in a field ecological aquaculture wastewater treatment system. This study will enlighten the role played by marine clam.in regulating the benthic microbial community assembly, nutrient dynamics, and organic effluent reduction in search of an optimized clam-sediment-bioengineering based bioremediation system, a recently rising and viable approach for aquaculture-wastes and contaminants management approach.

    2 Material and Methods

    2.1 Study Area and Experimental Design

    The ecological aquaculture wastewater treatment system studied is situated in Xiangshan Bay, in Ningbo, China. This project was built in 2017 purposely to improving and optimizing the aquaculture wastewater treatments. The studied system consists of 4 effluent treatment sections arranged serially: 1) Physical sedimentation; 2) Biofilms; 3) Bivalves; 4) Artificial wetlands (Lukwambe., 2018). The biofilters species (., bivalves and wetlands) were changed and the performance of the system was compared in different production seasons to get the optimized system. In this study, the clam cultured zone: about 800m2(2/5 of the total area) was chosen as the treatment area. And an area between biofilm and bivalve (clam) was chosen as the ‘control’. Eight partitions/enclosures, separated by PVC fabric mesh (100% polyester) were formed, 4 for control, and the other 4 for the treatment. Before the experiment, the bottoms of each enclosure were unified in order that they contain a similar initial condition. High vitality clams () with an average length of 1.8±0.3cm were stocked at an average stocking density of 360indm?2(Zhang., 2020). The average depth of each treatment enclosure was approximately 1.2m. and the size of the shrimp farm in which the effluents were produced is approximately 3300m2. The wastewater from the intensive aquaculture shrimp () ponds to the treatment system was discharged daily in the morning between 7:00–9:00 am. More than 165m3of effluent water was discharged daily. The average nutrients removal efficiency of the system is, 81.64% (total nitrogen (TN)), 82.34% (total phosphate (TP)), 75.68% (chemical oxygen demand), and 72.18% (ammonium (NH4+-N)) (Lukwambe., 2019b).

    2.2 Sample Collections

    2.2.1 Sediment microbial and biogeochemical samples

    During sampling, a total of 8 undisturbed sediment cores (approximately 8–9cm each), were collected using a handheld coring device (1.6m) with a removable round plexiglass tube (11cm i.d., 45cm height). Four for the control and four for the treatment to investigate 1) benthic biogeochemical nutrient fluxes, 2) sediment contents, and 3) the 16S rRNA microbial community assemblages. After sampling, the sediment cores (in the polypropylene tubes) and about 40L of thewater were brought back to the laboratory within 1h of collection. The plexiglass tubes were adapted with an injection port at 2cm above the sediment grind for drawing nutrient-flux water samples.

    2.2.2 Water column characteristics

    The dissolved oxygen (DO), temperature, salinity, and pH profiles of the water column were determinedusing a portable field multiparameter probe (YSI-550A oxygen meter, USA). Water column samples for phosphate (PO43?-P), nitrate (NO3?-N), nitrite (NO2?-N), and NH4+-N were also sampled for the determination of nutrient levels in the system.

    2.3 Biogeochemical Nutrient Incubation, Sampling and Flux Analytics

    The biogeochemical nutrient flux sampling and measurement procedures were carried as previously described (Zheng., 2009; Nicholaus and Zheng, 2014). Briefly, the sediment core tubes were slowly filled with fresh- water (from the site), and sealed with plastic stopper then incubated in water baths for 3h. Then rates of nutrient exchange were measured by extracting 45mL water samples with a plastic syringe at an interval of 1h through a hole bored across the tube at 1cm above the sediment surface. The concentration of collected nutrient samples was corrected for the composition of the water used for the replacement of extracted sample at each sampling time (Zheng., 2009). Extracted water samples for NH4+-N, NO2?-N, NO3?-N, and PO43?-P were collected in a screw cap polypropylene bottles, filtered (0.45μm cellulose acetate filter, GF) and immediately stored frozen (?20℃) until further analysis. To estimate the sediment oxygen flux, the DO level of the overlying water in each sediment core was measured at the beginning (0h) and the end (3h) of the incubation experiment using YSI- 550A oxygen meter (USA). Then all water samples for nutrient analyses were determined calorimetrically using a WESTCO SmartChem discrete analyzer (USA). The NO2?-N and NO3?-N were measured with the cadmium- copper reduction method. NH4+-N was measured by the indophenol blue method and the PO43?-P flux concentration was determined using the ammonium molybdate/ ascorbic acid method. The net biogeochemical nutrients and DO fluxes rates were calculated from the slopes of a linear regression concentration against time equation:

    where:(μmolm?2h?1) is the sediment O2or nutrients flux;(mgL?1) is the change in concentration of oxygen or nutrients (prior and after incubation);(m3) is the volume of the overlying water;(m2), is the cross-sectional area of the incubation chamber;(h) is the duration of incubation.

    2.4 Sediment TP, TOC, TN and OM

    The sediment samples were oven-dried at 80℃ to constant mass. The dried samples were ground to pass through a 1-mm mesh size before analyses. The analyses included the proportion of organic matter (OM) which was measured as a percentage loss upon ignition, LOI% for 4h at 550℃ after conversion of carbon. The TP contents were measured as previously described (Mudrich, 1996) while TN and total organic carbon (TOC) concentrations were determined by an Elemental NC Analyzer (Germany) after HCl acidification of the sediment samples to remove carbonate.

    2.5 Extraction and Sequencing of the Sediment Microbial DNA

    Post the nutrient flux incubation experiment, the intact sediment cores were sectioned at a different interval to make three cores samples: 0–2cm (surface), 2–4cm (middle), and 4–8cm (bottom). Before extraction, the sediment samples were freeze-dried by freeze dryer (Germany). Later, the genomic sediment DNA was extracted in triplicate from approximately 0.5g sediment of each replicate sample by using the manufacturer’s instructions given in the Power Soil DNA Isolation Kit (MoBio, USA). This process yielded a purified genomic sediment DNA for use in polymerase chain reaction (PCR) amplification. The extraction efficiency was improved by a disruptor (Bio101) FastPrep cell. DNA from the triplicate samples was pooled into a final volume of 50μL in Tris-EDTA buffer. DNA concentrations (quality and quantity) were quantified spectrophotometrically using an ND-2000 Nanodrop spectrometer (Thermo Fisher Scientific, Wilmington, DE) at 260nm absorption and 260/280 and 260/230nm absorbance ratios.

    High-throughput 16S rRNA gene sequencing was performed on the Illumina MiSeq platform (Illumina, San Diego, CA, USA). The V3–V4 hypervariable regions of the 16S rRNA gene were amplified by pyrosequencing to analyze the microbial community (composition, structure, and diversity) using fusion primers consisting of sequencer specific nucleotides, multiplex identifier tag, and template-specific nucleotide with the template specific sequences within the forward 341F (5’-CCTAYGGGRBG CASCAG-3’) and reverse 806R (5’-GGACTACNNG GGTATCTAAT-3’) barcoded PCR primers in 10ng DNA template. The V3-V4 hypervariable regions of the 16S ribosomal ribonucleic acid (rRNA) gene was amplified in triplicate PCR reactions (5μL volume in 50-μL reaction system). Briefly, according to Nicholaus. (2019b) the PCR cycle conditions for the 16S rRNA segment gene for each sample was carried out as follows: denaturation at 95℃ for 3min, and 27 cycles at 95℃ for 30s, 55℃ for 30s, 72℃ for 45s, plus a final extension of 10min at 72℃. Triplicate PCR amplicons were mixed after purification using a Takara purification kit (Takara, Japan). Libraries were then generated using TruSeq? DNA Sample Prep Kit for Illumina (Illumina, San Diego, CA, USA), following the manufacturer’s instructions. Library construction and sequencing were performed commercially (MajorBio Co. Ltd., Beijing).

    2.6 Bioinformatics Analysis

    The sequencing of paired reads was first joined with FLASH using default settings (Mago? and Salzberg, 2011). Then all analyses of the V3-V4 hypervariable regions of the microbial 16S rRNA amplicons (the Raw FASTQ files) were conducted by using the QIIME software package according to Caporaso. (2010). We demultiplexed and quality filtered the Raw 454 reads by removing low quality reads. Sequences shorter than 150 bp were discarded as well as the removal of the Roche adapters, linkers, primers, and sample barcodes (Caporaso., 2010). Microbial phylotypes were denoised to reduce possible sequencing errors and clustered at 97% identity in the operational taxonomic unit (OTU) using UCLUST (Edgar, 2010). One representative sequence of each cluster was aligned to the Greengenes database (release 13.8) with PyNAST 1.1 (DeSantis., 2006). To avoid artifacts resulting from the sample size, the chimeric sequences were detected using the Chimera Slayer algorithm in the QIIME 1.8.0 software package. OTUs observed in only one sample or represented by only one sequence (singleton OTUs) were removed. Only sequences flagged as non-chimeras by de novo and reference-based (Greengenes database) methods were retained for further analyses. Finally, the number of sequences assigned to each OTU was summarized in a table generated by QIIME. The OTU table produced by the quality filtering containing the taxonomic assignment of the OTUs were used for all downstream analyses.

    2.7 Statistical Analysis

    2.7.1 Microbial community composition and diversity

    Microbial community composition and diversity: Statistical approaches involving alpha diversity measures were calculated to provide information about the diversity’ within individual samples in all treatments (depths/sites). Sequences were normalized to the minimum number of reads per sample and alpha diversity parameters were computed using QIEME 1.8.0 (Caporaso., 2010). Community richness estimators; Chaol and Shannon diversity were estimated. Non-metric multidimensional scaling-NMDS, (metaMDS function, based on Bray-Curtis distance matrix) were adopted to visualize the dynamics of bacteria assemblage structure (Hallett., 2016). Before NMDS, Hellinger transformations were applied to microbial data. By using the function ‘a(chǎn)donis’ in the vegan package, the differences in bacterial communities among the sampling points were tested with permutational multivariate analysis of variance (PERMANOVA; Anderson, 2001).

    2.7.2 Microbial biomarker discovery and visualization

    The differentially significant taxonomic biomarkers/ features of the control and clam areas were identified by submitting the OTU lists of samples to the LEfSe (Linear discriminant analysis Effect Size) pipeline (LDA Effect size, http://huttenhower.sph.harvard.edu/galaxy/, Segata., 2011). This method applies robust statistical tests (Kruskal-Wallis) which involve biological consistency (Wilcoxon rank-sum tests, Segata., 2011) to predict taxonomic biomarkers based on effect-size estimation. Differentially abundant and biologically predominant taxa (phylum to genus level) were determined.

    2.7.3 Sediment nutrient fluxes

    A one way and 2-way analysis of variance (ANOVA) were used to test the variations in nutrient flux rates. Normal distribution and homogeneity of variance were tested for each data set, and data were log-transformed when necessary to fulfill the assumptions of ANOVA. Pearson’s correlation analysis was conducted to determine the links between the bacterial community and the environmental factors. Canonical Correspondence Analysis (CCA) was carried out to correlate the environmental variables with the sediment microbial community. Generally, all statistical analyses were performed using SPSS- software (SPSS, v16.0) and the open-source software R v.3.5.1 (R Core Team, 2017). The results were generally expressed as mean±standard error and differences were considered significant at alpha <0.05.

    3 Results

    3.1 Physicochemical Characteristics of the Sediment and Water

    The sediment cores from the control and treatment exhibited strong physical differences. The sediment from the clam zone was loose, soft, and fluffy whereas sediment texture in the control area appeared more compact, less fluffy, and coarser (physical touch and visual observations). Moreover, thecultured zone had increased brownish coloration suggesting an oxidized se- diment surface. The sediment TN, TOC, TP, and OM contents were markedly different between the treatment and the control area (ANOVA,<0.05; % dry weight-DW, Table 1). The TN, TOC, and TP contents in the clam bioturbated area were profoundly lower (<0.05). At 0–2 cm sediment dep- th (clam zone), the organic content was most altered accounting for 46% (TN), 17% (TOC), and 13% (TP) decrease compared to the control. Similarly, organic matter (OM) was reduced by 24% (DW, 0–2cm) relative to the control group.

    Lower sediment organic contents may signify an enhanced OM breakdown and transformation at the clam site. Nutrient levels in the water column between the control and the treatment area were differently distributed (Table 2).Other environmental variables such as DO, temperature, and salinity were measured. Temperature were normally 24.15±0.4℃ and salinity were 26.74ppt and the two showed similar values in both the control and the clam area. The concentration of DO fol- lowed similar trends with certain degrees of variations which were not significant among the treatments (>0.05).

    Table 1 The vertical sediment total nitrogen-TN, total organic carbon-TOC, total phosphorus-TP, and organic matter-OM distributions (%mean±SD) in each treatment during the present study. The different letter denotes significant differren-ces (P<0.05, ANOVA) among the same depth of different treatment groups and between depths of the same group

    Table 2 In situ nitrate (NO3?-N), nitrite (NO2?-N), ammonium (NH4+-N), phosphate (PO43?-P), Biological oxygen demand (BOD), Total nitrogen (TN), and Total phosphate (TP) distributions in the overlying water in the studied area before the incubation experiment

    3.2 Sediment-Water Fluxes of Nutrients and Sediment Oxygen Consumption (SOC)

    The net nutrient fluxes between control and treatment varied significantly (Fig.1,<0.05) during the benthic nutrient flux incubation experiment. The percentage flux rate change between the treatment and the control was 51%, 88%, 77%, and 45% for NH4+-N, NO3?-N, NO2?-N, and PO43?-P respectively (Fig.1). All the nutrients showed an out-flux release trend. Concerning SOC, the clam area showed the highest consumption with mean sediment oxygen uptake of 5.07±0.52μmol O2m?2h?1(treatment) compared to 1.94±0.21μmol O2m?2h?1(control) (Fig.1). The mean rate of SOC in the sediment cores (clam area) increased significantly compared to the control group (ANOVA,<0.05). A linear regression of the oxygen concentrations showed a significant change with time (<0.05). The SOC rate (mean±SD) was roughly 3.04-fold higher than that of the control.

    Fig.1 Sediment physicochemical nutrient and SOC fluxes during the incubation experiment between the control and the treatment. Data presented in μmol means±standard deviation (SD).

    3.3 Microbial Community Assemblage in the S. subcrenata Sediments

    Following a successful MiSeq sequencing, and after quality filtering a total of 18633 rRNA sequences representing 41584 OTUs () (clustered at 97% sequence similarity) with 48 phyla, 79 classes, 107 orders, 197 families, and 321 genera were detected amongst all samples. The number of sequences per sample ranged from 7935 to 82784, with a mean and standard deviation of 80436±26538. The highest average OTU counts (1610.91± 65.97) were obtained from the clam samples, compared to 544.39± 42.16 from the control (Table 3).

    Table 3 Mean±SD values for microbial abundance and alpha diversity variation present in the sediment samples collected from the clam and control during this study

    Notes: OTUs, Operational taxonomic units, Different superscript letters within the same row indicate significant differences (ANOVA;< 0.05)

    Further taxonomic analysis indicated 19 unique phyla (39.6%) were significantly abundant in the clam area, 10 unique phyla (20.8%) from the control, and 19 phyla (39.6%) shared in both groups. The prominent phylain the treatments were Proteobacteria, Firmicutes, Bacteroidetes, Verrucomicrobia, Nitrospirae, Chloroflexi, Plan- ctomycete (Fig.2a, b). The most abundant bacterial assemblies (phylum level, relative abundance > 0.5% in at least one sample) were observed in clam area (surface, 0–2cm) including Proteobacteria (41.3±5.24%), Chloroflexi (33.3±2.4%), Bacteroidetes (15.7±5.26%), Firmicute (9.05±1.62%), and WS3, (4.05±1.68 %) whereas same taxa existed in the control but with slightly lower reads (Proteobacteria; 32.5±15.26%, Chloroflexi; 26.3±11.13%, Bacteroidetes; 9.4±4.22%, Firmicute; 4.37±1.68%, and WS3; 1.12±1.21%) (Fig.2a) of the total 16S rRNA seque- nces.Deeper layers, 2–4 and 4–8cm were mostly dominated by Bacteroidetes, Nitrospira, Verrucomicrobia, Ch- loroflexi, and Firmicutes in both treatments but clam samples jointly shared 76% of the total gene sequence compared to control group (22%). At class level samples were enriched with Alphaproteobacteria (29.37%), Deltaproteobacteria (24.79%), and Bacteroides (18.29%) (Fig.2b).Analysis with LEfSe further indicated the significant distribution of microbial assemblies between the clam and the control (LEfSe, LD=3.5, Fig.3). The results of the LEfSe confirmed the tendency shown in Fig.3 (phylum to genus levels) and revealed that Acidobacteria (genus of Halophageles, Candidatus, and Solibacter), Bacteroidetes (genus Flavobacterium), Pseudomonas (genus of Pseudomonadales), and Proteobacteria (genus Novosphingobium) were significantly abundant in clam area.

    Fig.2 The relative abundance of the prominent bacterial taxa classi?ed on phylum level (a) and (b) at class level (including Proteobacteria sub-class) from the treatments and across the core vertical depths (0–8cm). Minor phyla accounting for less than 0.5% of total sequences are summarized in the group ‘Others’.

    Fig.3 The phylogenetic distribution of microbial lineages from the two study areas. Lineages with linear discriminating analysis values of 3.5 or higher determined by effect size measurements are displayed. The six rings of the cladogram stand for the domain (innermost), phylum, class, order, family, and genus. Enlarged circles in dark green and red are differentially abundant taxa identified as taxonomic biomarkers in the treatments (red=treatment, green=control).

    Fig.4 Non-metric multidimensional scaling (NMDS) ordination analysis plot of 16S rRNA amplicons of sediment-micro- bial-community composition of samples collected from the clam and control treatments. Community dissimilarity is based on the Bray-Curtis index. Shapes in a blue triangle, green circles, and red squares represent samples from the surface, middle and bottom respectively. Letters C and CL indicate control and clam samples respectively.

    3.4 Vertical Alpha Diversity of the Bacterial Community

    Bacterial diversity based on Chao1 (species richness), and Shannon (species richness and evenness), indexes of the observed OTUs varied significantly (ANOVA:<0.05, Table 3). The Simpson diversity index was non- significant (>0.05) throughout all samples, but with slightly higher reads in the clam than control samples. The Shannon index and Chao1 estimator confirmed that sediment samples at 0–2 and 2–4cm depth in the clam area experienced high microbial abundance and diversity compared to deeper depth (4–8cm) of the same treatment contrary to control group whereas no clear distribution patterns were observed. Generally, in all visualized samples, surface cores had the highest diversity indexes followed by middle and bottom cores in the clam treatment.

    Analysis with NMDS ordination (Fig.4) similarly showed a considerably high diverse structural and distribution of microbial community compositions in the clam area. Taxonomically (at the order level), the surface samples had greater abundances of Methylococcales, Cytophagales, SJA-15, and Rhodobacterales, while the bottom samples were enriched with Clostridiales, Desulfobacterales, and Myxococcales. At the genera level, taxa such as Nitrospina, Pseudomonas, Planctomyces, Desulfococcus, and Verrucomicrobium were evident after the heatmap plot of the 50 most prominent genera (Fig.5). To further describe the microbial assemblages, distribution, and stru- cture we carried out the analysis of similarity (ANOSIM) and PERMANOVA test that similarly indicated marked dissimilarities among the treatments, (ANOSIM,=0.707,=0.037 and PERMANOVA,<0.05, Table 4).

    Table 4 Bacterial community structure and composition determined by a permutational multivariate analysis of variance(PERMANOVA, ADONIS) for the treatment groups from the 16S rRNA sediment microbial community

    Fig.5 Heatmap of differentially abundant genus taxa identified in the sediment cores, at the surface, middle, and bottom among the treatments. Rows represent centered log ratio-transformed sequence counts for each of the samples.

    3.5 The Relationship Among Environmental Variables, Nutrient Fluxes, and the Bacterial Assemblies

    To evaluate the relationship between the bacterial communities, nutrient fluxes, and environmental variables, CCA and PCA were conducted in the overlying-water and the sediment depth (0–2, 2–4, and 4–8cm). CCA analysis uncovered nine environmental factors (EFs) which significantly correlated with the bacterial communities and jointly accounted for 89.7% of the total variation. OM, SOC, and NO2?-N contributed the most, followed by TOC, NH4+-N, NO3?-N, TP, and PO43?-P (Fig.6a). To further explore the interlinkages between bacterial community and the EFs, PCA was drawn and the two axes jointly explain 96.1% of the microbial variance at the genus level (Fig.6b). The correlations between bacterial assemblages and EFs in both ordinations were high, indicating a strong relationship between the genus taxa and EFs. Factors such as SOC, OM, and NO2?-N were amongst the strongest correlations between the bacterial communities and the EFs. Several potential genera like,,,andwere among the taxa that highly correlated with nutrient fluxes (Figs.6a, b).

    Fig.6 (a) Canonical correlation analysis (CCA) ordination plot of the microbial community composition and physicochemical variables within the treatments. (b) Principle Component Analysis (PCA) plots showing the relationship between differentially abundant genera and environmental parameters among the treatment samples during the study.

    4 Discussion

    Benthic macrofauna affects inorganic nutrients and microbial (processes, assemblages, and structure) in sediments due to their bioengineering activities (., Papaspyrou., 2006; Nicholaus., 2019a, b; Zhang., 2020). They are directly involved in sedimentary effluent bioremediation processes through benthic OM mineralization and aerobic metabolization (Lukwambe., 2020, Nicholaus., 2019b). These ecosystem-engi- neered roles tend to create a sustainable and ecologically balanced aquatic community. Besides, the disturbances and sediment reworking activities exerted by benthic macrofauna like bivalve filter-feeders can explain the potential mechanistic ecological process involved in microbial community assemblage, waste degradation, and the nutrient dynamics in aquatic systems (Shen., 2017; Nicholaus., 2019b).

    4.1 The Impacts of Clams in the Sediment During Bioremediation Activities

    The marine clamsubstantially influenced the sediment condition physiochemically and biologically. Compared to the control, the sediment in the clam area was composed of more fine particles, muddy, and silt. The organic contents associated with TN, TP, and TOC were significantly low in the clam area relative to the control (Table 1), which suggests a reduction of organic contents among the clam treatment. This finding corresponds with previous studies involving clam’s bioremediation in the sediments of aquaculture effluents (Nicho- laus., 2019b). Normallyactivities are associated with several processes such as burrowing, feeding, bioirrigation, biodeposition, resuspension, and mineralization. All these activities are believed to causing biogeochemical nutrient cycling which modifies the distribution of microorganisms responsible for bioremediation (Nicholaus., 2019b) and other microbial metabolic pathways such as sulfate reduction, sulfide oxidation, denitrification and nitrification (Kristensen and Kostka, 2005). The marine clamas a potential bioturbator induces a suitable microenvironment for OM decomposing bacteria such as(Lukwambe., 2020). Similarly, Satoh. (2007) reported a significant denitrification enhancement with ammonium oxidizing bacteria in burrows constructed by macrofauna. Microorganisms., bacteria play major roles in decomposing waste OM and removing carbonaceous compounds in the environment (Bender and Phillips, 2004). Because of the improved distribution of oxygenated water in the sediment, this may lead to increased microbial sediment OM degradations and reduction of organic effluents. These results are consistent with the previous findings which showed that the OM contents in sediments inhabited by razor clam (Lukwambe., 2018) and blood clam (Nicholaus., 2019b) were generally lower than those without

    Furthermore,this study revealed substantial impacts on the benthic biogeochemical nutrient fluxes and/or cycling in the bioturbated area. Increased release of ammonium, nitrate, nitrite, and phosphate from the sediment to the overlying water was observed accounting for 51%, 88%, 77%, and 45% respectively, higher in the treatment relative to the control (Fig.1). High rates of sediment- NH4+-N and PO43?-P release due to bivalve/clam activities have been previously reported (Nicholaus., 2019a; Zhao., 2019; Zhang., 2020). The level of NH4+- N, NO2?-N, and NO3?-N concentration is linked to the changes in oxygen concentration and abundance of microorganisms observed in this study. Similarly, its inferred that the relatively higher biogeochemical nutrient fluxes in the treatment area are attributable to the nutrient mineralization activities enhanced by the clams during bioturbation. Enhanced oxygen penetration mediated via macrofaunal activities (., burrowing) is reported to increase the bioavailability of benthic elements including phosphorus fractions (Nicholaus and Zheng, 2014; Nicho- laus., 2019a). Bioturbation can widely accelerate benthic process including nutrient fluxes, nitrification rates, shift of OM to deeper sediment layers and enhanced material cycling (Nicholaus., 2020).

    The SOC rates in the clam area were higher compared to that of the control (Fig.1, Table 2). The increase in oxygen influxes in the clam areas is attributable toburrowing activities in the sediments. Sediment rework by the bioturbating organism creates oxygen demand due to microbial activities such as waste degradation and mineralization. Similar results have been previously reported on the effects of clams on oxygen supply in the sediments (., Bertics and Ziebis, 2009; Ni- cholaus., 2019b). In the presence of the anaerobic condition, the process of nitrification is increased through microbial nitrate transformation which may result in a reduction of excess nitrogen and ammonium in the ecosystem.

    4.2 Role of S. subcrenata on the Sedimentary Bacterial Assemblages

    Macrofauna bioturbations create a fine spatial scale selective pressure that modulates significant rearrangements of sediment bacteria community assemblages and interactions (Booth., 2019). Similarly, burrows created by macrofauna can highly contain numerous potential bacteria assemblages with high abundance (Papaspyrou., 2006; Stauffert., 2014). In the present study, 48 unique bacterial phyla were widely differentiated and distributed across all samples. The clam inhabited area demonstrated a high diversity of the microbial community with significantly higher abundance of several potential bacterial phyla, such as Firmicutes, Bacteroidetes, Proteobacteria, Chloroflexi, Cyanobacteria, Lentisphaerae, Nitrospirae, Planctomycetes, Actinobacteria, Gemmatimonadetes, KSB3, and Verrucomicrobia as referred to control (Figs.2a, b). Among these taxa, the Verrucomicrobia, Chloroflexi, Firmicutes, and Planctomycetes, are one of the common taxa present in marine sediment and capable of degrading organic wastes (Lukwambe., 2018; Nicholaus., 2019b). This suggests thatimproved the sediment quality and exerted a measurable change in microbial diversity and abundance.is an intertidal clam that can potentially burrow up to 4cm deep (Nicholaus., unpublished). These burrow formations may be enriched with high oxygen-rich overlying water drawn from the surface (Bertics and Ziebis, 2009). Dissolved oxygen plus sediment median grain size, have been reported as the vital factor shaping benthic bacterial communities and OM mineralization.

    Further, the distribution of some functional groups/ strains associated with sulfur-reducing bacteria in the families Desulfobacteraceae and Desulfobulbaceae (phylum: Proteobacteria, class Deltaproteobacteria) were highly represented in different OTUs, in the clam area. Again, Novosphingobium of the class Alphaproteobacterial was expressed in all samples but more explicitly in S. subcrenata area signifying that denitrification was accelerated (Feng., 2012). The resultant bivalve bioengineering activities (burrow-walls, irrigation, ventilated-holes) can alter how the microbial assemblages are structured in the sediment-ecosystem (Papaspyrou., 2006; Bertics and Ziebis, 2009). The distinct formation of microbial assemblage structures and increased diversity in the clamarea is in correspondence with other previous clam bioturbation studies (Lukwambe., 2019a, b; Nicholaus., 2019a, b). According to Ieno. (2006) and Solan. (2008), bioturbation functional intensity increased benthic species richness a process that may result in improved sediment restoration and nutrient fluxes. Mermillod-Blondin. (2004) demonstrated that, due to its strong ecosystem engineering activities, the occurrence ofincreased total active bacterial abundance in benthic sediments. Also, Chen. (2017) showed that bioturbations by macrofauna can oxygenate the sediment, affect bacterial composition, and activity.

    Furthermore, the ecological process of thecan be described indirectly by the intermittent bio-irrigation of burrows that creates an oscillating oxic- anoxic ecosystem and enriching coexistence of oxic-anoxic microbial taxa (Volkenborn., 2012). The biodeposits from the clams (Fig.7) may result in high organic wastes enriched sediment which saves as food for the bacteria. Learman. (2016) demonstrated OM to be a driver of benthic microbial community structure across the surface sediment. More recently, Nicholaus. (2019b), suggested that nutrient cycling can be an important factor in altering the microorganism assemblies.

    Fig.7 A hypothetical sketch of sediment-bacteria-macrofauna ecological processes involved in the nutrient dynamics in the bivalve bioturbated sediment ecosystem.

    4.3 Relationship Between the Microbial Community and Nutrient Dynamics

    Microorganisms contribute significantly to the balancing and cycling of marine biogenic elements (Azam and Malfatti, 2007; Dyksma., 2016). In our study, the first two ordination axes of CCA, and PCA plot explained 42.3%, and 22.9%, 73.2% and 22.9% of the total variation in the bacterial community structure (Fig.6a). The results of CCA showed, the environmental variables such as nutrient fluxes, TOC, and SOC parameters, were correlated with bacterial communities (genus level, relative abundance >3%). The correlations between bacterial communities and nutrient fluxes were high, indicating a strong relationship between the microbial community and environmental variables. Further, an association was revealed between several bacterial taxa, associated with nitrification and denitrification such as Proteobacteria, Nitrospirae, Firmicutes, Bacteroidetes, Pseudomonas, and the benthic nutrient fluxes (Figs.5, 6b). The presence of these taxa signifies that probably the dinitrogen nitrate reduction to ammonium activities in the clam sediments was promoted (Moncada., 2019). Potential sulfate and phosphate-reducing taxa such as Deltaproteobacteria and Desulfomicrobium) were promoted in this study (Fig.6a, b) Desulfomicrobium can contribute indirectly to inorganic nutrient release (.., PO43?-P) via organic matter mineralization (Sinkko., 2011; Ki., 2018). The CCA results indicated that TOC, OM, NH4+-N, and NO3?-N were among factors correlated with several types of sediment bacteria. The OM in the sediment is expected to decompose into phosphate (PO43?-P) by the sulfate-reducing bacteria (SRB), such as,,, and, which are found in both the upper and deeper sediment layers. In most samples, bacteria relating to key functions such as Desulfatibacillum, Burkholderiales, Desulfobacterium, Desulfomicrobium, and Desulfosalsimonas, were dominant, and positively expressed in the treatments (Figs.5, 6b).,,play a key role in decomposing organic phosphorus indirectly to phosphate via OM mineralization (Zhao., 2012; Ki., 2018). Therefore, the close correlation among the different microbial functional genes and strong correlations between PCA and CCA in the clam treatment may suggest a strong link between benthic nutrient cycling processes and the microbial community which is key aspect in bioremediation systems The diverse microbial community and increased nutrient dynamics in the clam treatment is due to the fact that ecosystem engineering processes by benthic macrofauna bioturbation can cause steep chemical gradients along burrows, resulting in accelerated mineralization, sulfide oxidation, sulfate reduction, and denitrification coupled nitrification activities (Meysman., 2006).has been reported to highly alter the microbial community structure and increased nutrients flux and cycling in benthic sediments (Shen., 2017; Lukwambe., 2019b; Nicholaus., 2019b). Other clams with similar activities such as(Liu., 2007),(Nicholaus and Zheng, 2014), and(Zhao., 2019) has been shown to contribute on benthic nutrient dynamics which impact bioremediation. In the current study, due to close similarity in morphology and burrowing behavior,has equally expressed more or less similar roles. It’s well-known that the transformation and cycling of biogeochemical nutrients such as nitrogen and phosphorus are driven by sediment microorganisms (Hou., 2013). Decomposition processes involvingandcan have a strong effect on ammonification pathways (Dauda., 2014). Fig.7 describes the existing relationship between sediment bacterial community degradation processes, sediment contents, andbioturbation activities. Since the presence ofpositively accelerated the microbial biomass turnover may suggest that mineralization processes of the nitrogen associated compounds (ammonium, nitrite, and nitrate) and other organic components were similarly accelerated thus prompting high flux rates, nutrient cycling, reduced organic wastes, and improved response of microbial activities to nutrient load.

    5 Conclusions

    The current study demonstrates the contribution of clamon the microbial community assemblages and nutrient dynamics in the sediments of the aquaculture bioremediation system. The clam treated area showed a higher diverse range of microbial community assembly. Specifically, the,,, andwere significantly expressed bacterial taxa among the treatments. Substantially theaccelerated the flux rate of NH4+-N, NO3?-N and PO43?-P and SOC. Moreover, organic contents (TP, TOC, TP and TN) especially from the surface cores in the clam treated sediments were significantly lower. Our findings suggest thatbioturbation activities predominantly increase microbial community diversity and nutrient dynamics, promotes organic content transformations, and improve bioremediation activities in the sediments by indirectly affecting the nitrification-denitri- fication activities, and modified sediment micro-habitats.

    Acknowledgements

    This study was supported by the Zhejiang Public Welfare Technology Research Program of China (ZPWTP) (No. LGN18C190008),the Zhejiang Provincial Natural Science Foundation of China (No. LQ20C190003), the Ningbo Municipal Natural Science Foundation (Nos. 2019A610 421, 2019A610443), and the K. C. Wong Magna Fund in Ningbo University.

    Aller, R. C., 2001. Transport and reactions in the bioirrigated zone. In:. Boudreau, B., and J?rgensen, B. B., eds., Oxford University Press, New York, 269-301.

    Anderson, M. J., 2001. A new method for non-parametric multivariate analysis of variance., 26: 32-46.

    Azam, F., and Malfatti, F., 2007. Microbial structuring of marine ecosystems., 5: 782-791.

    Bender, J., and Phillips, P., 2004. Microbial mats for multiple applications in aquaculture and bioremediation., 94: 229-238.

    Bertics, V. J., and Ziebis, W., 2009. Biodiversity of benthic microbial communities in bioturbated coastal sediments is controlled by geochemical microniches., 3: 1269-1285.

    Booth, J. M., Fusi, M., Marasco, R., Mbobo, T., and Daffonchio, D., 2019. Fiddler crab bioturbation determines consistent changes in bacterial communities across contrasting environmental conditions., 9 (1): 3749.

    Brito, L. O., Cardoso Junior, L. de O., Lavander, H. D., Abreu, J. L., de Severi, W., and Gálvez, A. O., 2018. Bioremediation of shrimp biofloc wastewater using clam, seaweed, and fish., 34 (10): 901-913.

    Caporaso, J. G., Bittinger, K., Bushman, F. D., DeSantis, T. Z., Andersen, G. L., and Knight, R., 2010. PyNAST: A flexible tool for aligning sequences to a template alignment., 26: 266-267.

    Chen, X., Andersen, T. J., Morono, Y., Inagaki, F., J?rgensen, B. B., and Lever, M. A., 2017. Bioturbation as a key driver behind the dominance of bacteria over archaea in near-surface sediment., 7 (1): 2400.

    Dale, A. W., and Prego, R., 2002. Physico-biogeochemical controls on the benthic-pelagic coupling of nutrient fluxes and recycling in a coastal upwelling system., 235: 15-28.

    Dauda, A. B., and Akinwole, A. O., 2014. Interrelationships among water quality parameters in recirculating aquaculture system., 8: 20-25.

    DeSantis, T. Z., Hugenholtz, P., Keller, K., Brodie, E. L., Larsen, N., Piceno, Y. M., Phan, R., and Andersen, G. L., 2006. NAST: A multiple sequence alignment server for comparative analysis of 16S rRNA genes., 34: 394-399.

    Diaz, R. J., 2001. Overview of hypoxia around the world., 30: 275-281.

    Dyksma, S., Bischof, K., Fuchs, B. M., Hoffmann, K., Meier, D., Meyerdierks, A., Pjevac, P., Probandt, D., Richter, M., and Stepanauskas, R., 2016. Ubiquitousdominate dark carbon fixation in coastal sediments., 10: 1939-1953.

    Edgar, R. C., 2010. Search and clustering orders of magnitude faster than BLAST., 26 (19): 2460-2461.

    Feng, S., Xie, S., Zhang, X., Yang, Z., Ding, W., Liao, X., Liu, Y., and Chen, C., 2012. Ammonium removal pathways and microbial community in GAC-sand dual media filter in drinking water treatment., 24: 1587-1593.

    Hallett, L. M., Jones, S. K., MacDonald, A. A. M., Jones, M. B., Flynn, D. F. B., Ripplinger, J., Slaughter, P., Gries, C., and Collins, S. L., 2016. An r package of community dynamics metrics., 7: 1146-1151.

    Hou, J., Song, C., Cao, X., and Zhou, Y., 2013. Shifts between ammonia-oxidizing bacteria and archaea in relation to nitrification potential across trophic gradients in two large Chinese lakes (Lake Taihu and Lake Chaohu)., 47: 2285-2296.

    Ieno, E. N., Solan, M., Batty, P., and Pierce, G. J., 2006. How biodiversity affects ecosystem functioning: Roles of infaunal species richness, identity and density in the marine benthos., 311: 263-271.

    Karakassis, I., Tsapakis, M., Hatziyanni, E., Papadopoulou, K. N., and Plaiti, W., 2000. Impact of cage farming of fish on the seabed in three Mediterranean coastal areas., 57: 1462-1471.

    Kautsky, N., and Evans, S., 1987. Role of biodeposition byin the circulation of matter and nutrients in a Baltic coastal ecosystem., 38: 201-212.

    Ki, B. M., Huh, I. A., Choi, J. H., and Cho, K. S., 2018. Relationship of nutrient dynamics and bacterial community structure at the water-sediment interface using a benthic chamber experiment., 53 (5): 482-491.

    Kristensen, E., and Kostka, J. E., 2005. Macrofaunal burrows and irrigation in marine sediment: Microbiological and biogeochemical interactions. In:. Kristensen, E.,., eds., American Geophysical Union, Washington, D. C., 125-157.

    Learman, D. R., Henson, M. W., Thrash, J. C., Temperton, B., Brannock, P. M., Santos, S. R., Mahon, A. R., and Halanych, K. M., 2016. Biogeochemical and microbial variation across 5500km of Antarctic surface sediment implicates organic matter as a driver of benthic community structure., 7: 284.

    Liu, J., Chen, Z., Xu, S., and Zheng, X., 2007. Experimental research on the impact ofon DIN exchange at a tidal flat sediment-water interface., 25 (4): 434-443.

    Lukwambe, B., Nicholaus, R., Yang, W., and Zheng, Z., 2019a. Blood clams()bioturbation alter succession of bacterioplankton community and nutrients removal performance in an aquaculture wastewater bioremediation system.,734520.

    Lukwambe, B., Nicholaus, R., Zhao, L., Yang W., Zhu, J., and Zheng, Z., 2020. Microbial community and interspecies interaction during grazing of ark shell bivalve () in a full-scale bioremediation system of mariculture effluents., 158: 104956.

    Lukwambe, B., Yang, W., Zheng, Y., Nicholaus, R., Zhu, J., and Zheng, Z., 2018. Bioturbation by the razor clam () on the microbial community and enzymatic activities in the sediment of an ecological aquaculture wastewater treatment system., 643: 1098-1107.

    Lukwambe, B., Zhao, L., Nicholaus, R., Yang, W., Zhu, J., and Zheng, Z., 2019b. Bacterioplankton community in response to biological filters (clam, biofilm, and macrophytes) in an integrated aquaculture wastewater bioremediation system., 254: 113035.

    Mago?, T., and Salzberg, S. L., 2011. FLASH: Fast length adjustment of short reads to improve genome assemblies., 27: 2957-2963.

    Mermillod-Blondin, F., Rosenberg, R., Fran?ois-Carcaillet, F., Norling, K., and Mauclaire, L., 2004. Influence of bioturbation by three benthic infaunal species on microbial communities and biogeochemical processes in marine sediment., 36: 271-284.

    Meysman, F. J. R., Middelburg, J. J., and Heip, C. H. R., 2006. Bioturbation: A fresh look at Darwin’s last idea., 21: 688-695.

    Moncada, C., Hassenrück, C., G?rdes, A., and Conaco, C., 2019. Microbial community composition of sediments influenced by intensive mariculture activity., 19 (12): 1-12.

    Moriarty, D. J. W., 1997. The role of microorganisms in aquaculture ponds., 151 (1-4): 333-349.

    Mudrich, A., Azcue, J. M., and Mudroch, P., 1996.. CRC Press Inc., Boca Raton, Florida, USA, 1-12.

    Navarrete-Mier, F., Sanz-Lázaro, C., and Marín, A., 2010. Does bivalve mollusc polyculture reduce marine finfish farming environmental impact?, 306 (1-4): 101-107.

    Nicholaus, R., and Zheng, Z. M., 2014. The effects of bioturbation by the Venus clamon the fluxes of nutrients across the sediment-water interface in aquaculture ponds., 22 (2): 913-924.

    Nicholaus, R., Lukwambe, B., Lai, H., Yang, W., and Zheng, Z., 2019a. Nutrients cycling in ecological aquaculture wastewater treatment systems: Vertical distribution of benthic phosphorus fractions due to bioturbation activity by, 11: 469-480.

    Nicholaus, R., Lukwambe, B., Zhao, L., Yang, W., Zhu, J., and Zheng, Z., 2019b. Bioturbation of blood clamon benthic nutrient fluxes and microbial community in an aquaculture wastewater treatment system., 142:73-82.

    Papaspyrou, S., Gregersen, T., Kristensen, E., Christensen, B., and Cox, R., 2006. Microbial reaction rates and bacterial communities in sediment surrounding burrows of two nereidid polychaetes (and)., 148: 541-550.

    R Core Team, 2017. R: A Language and Environment for Statistical Computing. https://www.R-project.org/.

    Satoh, H., Nakamura, Y., and Okabe, S., 2007. Influences of in- faunal burrows on the community structure and activity of ammonia-oxidizing bacteria in intertidal sediments., 73: 1341-1348.

    Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., and Huttenhower, C., 2011. Metagenomic biomarker discovery and explanation., 12: R60.

    Shen, H., Jiang, G., Wan, X., Li, H., Qiao, Y., Thrush, S., and He, P., 2017. Response of the microbial community to bioturbation by benthic macrofauna on intertidal flats., 488: 44-51.

    Sinkko, H., Lukkari, K., Jama, A. S., Sihvonen, L. M., and Si- vonen, K., 2011. Phosphorus chemistry and bacterial community composition interact in brackish sediments receiving agricultural discharges., 6: e21555.

    Solan, M., Batty, P., Bulling, M. T., and Godbold, J. A., 2008. How biodiversity affects ecosystem processes: Implications for ecological revolutions and benthic ecosystem function., 2: 289-301.

    Stauffert, M., Duran, R., Gassie, C., and Cravo-Laureau, C., 2014. Response of archaeal communities to oil spill in bioturbated mudflat sediments., 67: 108-119.

    Volkenborn, N., Polerecky, L., Wethey, D. S., DeWitt, T. H., and Woodin, S. A., 2012. Hydraulic activities by ghost shrimpinduce oxic-anoxic oscillations in sediments., 455: 141-156.

    Yang, J., Gou, Y., Fang, F., Guo, J., Lu, L., Zhou, Y., and Ma, H., 2018. Potential of wastewater treatment using a concentrated and suspended algal-bacterial consortium in a photo mem- brane bioreactor., 335: 154- 160.

    Zhang, L., Shen, Q. S., Hu, H. Y., Shao, S. G., and Fan, C. X., 2011. Impacts ofon oxygen uptake and nutrient fluxes across the sediment-water Interface., 220: 399-411.

    Zhang, S., Fang, X., Zhang, J., Fang, Y., Hu, Z., Lizhen, W., and Daisuke, K., 2020. The effect of bioturbation activity of the ark clamon the fluxes of nutrient exchange at the sediment-water interface., 19: 232-240.

    Zhao, G., Du, J., Jia, Y., Lv, Y., Han, G., and Tian, X., 2012. The importance of bacteria in promoting algal growth in eutrophic lakes with limited available phosphorus., 42: 107-111.

    Zhao, L., Zheng, Y., Nicholaus, R., Lukwambe, B., Zhu, J., Yang, W., and Zheng, Z., 2019. Bioturbation by the razor clamaffects benthic nutrient fluxes in aquaculture wastewater treatment ecosystems., 11: 87-96.

    Zheng, Z. M., Dong, S. L., and Bai, P. F., 2009. Study on nutrient fluxes across sediment-water interface in different sea cucumber experimental enclosures., 39: 209-214 (in Chinese with English abstract).

    . E-mail:zhengzhongming@nbu.edu.cn

    January 21, 2020;

    May 10, 2020;

    September 24, 2020

    (Edited by Ji Dechun)

    嫩草影视91久久| 男女视频在线观看网站免费| 在线观看免费午夜福利视频| 国产 一区 欧美 日韩| 亚洲最大成人中文| 亚洲国产精品sss在线观看| 九九热线精品视视频播放| 一边摸一边抽搐一进一小说| 欧美日韩亚洲国产一区二区在线观看| 国产午夜精品论理片| 欧美一区二区精品小视频在线| 九九久久精品国产亚洲av麻豆 | 亚洲国产色片| 国产人伦9x9x在线观看| 手机成人av网站| 亚洲欧美日韩无卡精品| 亚洲九九香蕉| 在线看三级毛片| 国产精品野战在线观看| 人妻夜夜爽99麻豆av| 亚洲,欧美精品.| 成人无遮挡网站| 欧美中文综合在线视频| 国语自产精品视频在线第100页| 一级a爱片免费观看的视频| 成在线人永久免费视频| 很黄的视频免费| 国产激情偷乱视频一区二区| 999久久久国产精品视频| 99精品欧美一区二区三区四区| 国模一区二区三区四区视频 | www.999成人在线观看| 在线观看午夜福利视频| 国产 一区 欧美 日韩| 欧美日韩综合久久久久久 | 桃色一区二区三区在线观看| 亚洲成av人片在线播放无| 91字幕亚洲| 99热只有精品国产| 后天国语完整版免费观看| 最近最新中文字幕大全电影3| 黄色日韩在线| 国产精品久久久久久精品电影| 国产精品久久久久久人妻精品电影| 国产成人欧美在线观看| 久久久久国产精品人妻aⅴ院| 亚洲国产精品sss在线观看| 欧美性猛交黑人性爽| 51午夜福利影视在线观看| 免费人成视频x8x8入口观看| 嫁个100分男人电影在线观看| 后天国语完整版免费观看| 久久久久性生活片| 国产激情久久老熟女| 淫妇啪啪啪对白视频| 成人精品一区二区免费| 精品福利观看| 国产一级毛片七仙女欲春2| 国产精品影院久久| 男女做爰动态图高潮gif福利片| 久久天堂一区二区三区四区| 亚洲中文av在线| 日日干狠狠操夜夜爽| а√天堂www在线а√下载| 高清毛片免费观看视频网站| 国产69精品久久久久777片 | 国模一区二区三区四区视频 | 久久久国产精品麻豆| 国产麻豆成人av免费视频| 免费一级毛片在线播放高清视频| 亚洲精品在线观看二区| av在线天堂中文字幕| 夜夜夜夜夜久久久久| 亚洲无线观看免费| 欧美在线黄色| 国模一区二区三区四区视频 | 精品久久久久久久毛片微露脸| 很黄的视频免费| 国产午夜精品论理片| 成人三级黄色视频| 最近最新免费中文字幕在线| 国产在线精品亚洲第一网站| 午夜视频精品福利| av欧美777| 成人国产综合亚洲| 亚洲中文字幕日韩| 亚洲精品一区av在线观看| 国产主播在线观看一区二区| 日韩三级视频一区二区三区| 亚洲七黄色美女视频| 好男人在线观看高清免费视频| aaaaa片日本免费| 99久久成人亚洲精品观看| 中文字幕久久专区| 老司机在亚洲福利影院| 久久精品人妻少妇| 亚洲国产欧美人成| av片东京热男人的天堂| 国产精品99久久久久久久久| 男女做爰动态图高潮gif福利片| 亚洲精品久久国产高清桃花| 欧美最黄视频在线播放免费| 国产野战对白在线观看| 国产精品久久久久久人妻精品电影| 黄色丝袜av网址大全| 日韩中文字幕欧美一区二区| 欧美日韩中文字幕国产精品一区二区三区| 免费在线观看成人毛片| 亚洲欧美日韩东京热| 男人舔女人下体高潮全视频| 香蕉国产在线看| 国产激情欧美一区二区| 久久久国产欧美日韩av| 午夜福利在线观看吧| xxx96com| 午夜视频精品福利| 亚洲欧美一区二区三区黑人| 亚洲第一电影网av| 精品国产亚洲在线| 国产三级黄色录像| 中出人妻视频一区二区| 亚洲人成网站在线播放欧美日韩| АⅤ资源中文在线天堂| 欧美丝袜亚洲另类 | 69av精品久久久久久| www.熟女人妻精品国产| 国内少妇人妻偷人精品xxx网站 | 国产精品久久久久久久电影 | 欧美日韩一级在线毛片| 老司机深夜福利视频在线观看| 在线观看一区二区三区| 俺也久久电影网| 脱女人内裤的视频| 国产成人系列免费观看| a在线观看视频网站| 欧美最黄视频在线播放免费| 亚洲国产精品999在线| 亚洲五月天丁香| 亚洲在线自拍视频| 亚洲乱码一区二区免费版| 亚洲av熟女| 最新美女视频免费是黄的| 麻豆久久精品国产亚洲av| 九九热线精品视视频播放| 国产视频内射| 亚洲国产日韩欧美精品在线观看 | 18禁国产床啪视频网站| 69av精品久久久久久| 法律面前人人平等表现在哪些方面| 国产精品久久久久久久电影 | 老司机午夜十八禁免费视频| 精品福利观看| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久成人av| 久久久久久人人人人人| 欧美在线一区亚洲| 91老司机精品| 亚洲欧美日韩无卡精品| 欧美乱妇无乱码| 久久久久久久久中文| 啦啦啦免费观看视频1| 99riav亚洲国产免费| 真人做人爱边吃奶动态| 九九久久精品国产亚洲av麻豆 | 国产视频一区二区在线看| 亚洲熟妇熟女久久| 欧美一级毛片孕妇| 午夜福利在线观看免费完整高清在 | 日韩国内少妇激情av| 午夜免费观看网址| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 1024香蕉在线观看| 国产日本99.免费观看| 曰老女人黄片| 亚洲熟妇熟女久久| 欧美乱妇无乱码| 精品无人区乱码1区二区| 麻豆国产97在线/欧美| 欧美激情在线99| 国产成人一区二区三区免费视频网站| 国产成人啪精品午夜网站| 欧美绝顶高潮抽搐喷水| 亚洲av成人不卡在线观看播放网| 麻豆国产97在线/欧美| 性色avwww在线观看| 成人鲁丝片一二三区免费| 国产激情偷乱视频一区二区| 国产亚洲av高清不卡| 精品国产三级普通话版| 日本免费a在线| 久久久国产成人免费| 亚洲人成网站高清观看| 午夜两性在线视频| 国产伦精品一区二区三区视频9 | 三级毛片av免费| 国产精品国产高清国产av| 国产高清视频在线观看网站| 成熟少妇高潮喷水视频| 日韩欧美 国产精品| 一夜夜www| 国产成人av激情在线播放| 国产高清激情床上av| 一本久久中文字幕| 免费在线观看视频国产中文字幕亚洲| 88av欧美| 法律面前人人平等表现在哪些方面| 午夜福利高清视频| 欧美xxxx黑人xx丫x性爽| 小蜜桃在线观看免费完整版高清| 亚洲,欧美精品.| 久久久国产成人精品二区| 亚洲电影在线观看av| 日韩欧美免费精品| 黑人欧美特级aaaaaa片| 18禁黄网站禁片免费观看直播| 久久九九热精品免费| 免费在线观看视频国产中文字幕亚洲| 琪琪午夜伦伦电影理论片6080| 欧美另类亚洲清纯唯美| 亚洲精品一区av在线观看| 男插女下体视频免费在线播放| 久久久久免费精品人妻一区二区| 久久精品影院6| 欧美一级毛片孕妇| av天堂在线播放| 午夜免费激情av| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| 亚洲欧美日韩无卡精品| 1024手机看黄色片| 又紧又爽又黄一区二区| av在线天堂中文字幕| xxx96com| 99re在线观看精品视频| x7x7x7水蜜桃| 悠悠久久av| 免费人成视频x8x8入口观看| 亚洲av片天天在线观看| 十八禁网站免费在线| 精品日产1卡2卡| 国产精品综合久久久久久久免费| 国产成人影院久久av| 国产精品亚洲美女久久久| 老汉色av国产亚洲站长工具| 天堂影院成人在线观看| 男人舔女人下体高潮全视频| 成人鲁丝片一二三区免费| 免费在线观看影片大全网站| 男人的好看免费观看在线视频| 国产伦在线观看视频一区| 亚洲av美国av| 色尼玛亚洲综合影院| 欧美一区二区精品小视频在线| 色精品久久人妻99蜜桃| 欧美日韩瑟瑟在线播放| 叶爱在线成人免费视频播放| 国产1区2区3区精品| 色综合亚洲欧美另类图片| 欧美中文日本在线观看视频| 欧美+亚洲+日韩+国产| 十八禁人妻一区二区| 日韩国内少妇激情av| 久久精品综合一区二区三区| 淫妇啪啪啪对白视频| 一区福利在线观看| 亚洲第一欧美日韩一区二区三区| 精品人妻1区二区| 1024手机看黄色片| 午夜福利在线观看免费完整高清在 | 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 九九热线精品视视频播放| 国产精品一区二区免费欧美| 国产免费av片在线观看野外av| 国产亚洲欧美98| 一区二区三区高清视频在线| 青草久久国产| 精品不卡国产一区二区三区| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 19禁男女啪啪无遮挡网站| 一本精品99久久精品77| 黄色女人牲交| 久久伊人香网站| 2021天堂中文幕一二区在线观| 久久久国产成人免费| 国产亚洲av高清不卡| 国产精品自产拍在线观看55亚洲| 成人三级做爰电影| 不卡一级毛片| 亚洲精品粉嫩美女一区| 亚洲精品456在线播放app | 日本与韩国留学比较| 91在线精品国自产拍蜜月 | 精品久久久久久久人妻蜜臀av| 午夜福利在线在线| 欧美成人性av电影在线观看| 欧美一级毛片孕妇| 老司机在亚洲福利影院| 丰满人妻一区二区三区视频av | 国产私拍福利视频在线观看| 一本久久中文字幕| 99精品在免费线老司机午夜| 国产三级在线视频| 日本 欧美在线| 两个人看的免费小视频| 香蕉av资源在线| 久久久国产成人精品二区| 婷婷精品国产亚洲av| 欧美黑人欧美精品刺激| 俺也久久电影网| 欧美一级a爱片免费观看看| 一级黄色大片毛片| 国产欧美日韩一区二区三| 成人无遮挡网站| 男女那种视频在线观看| 国内少妇人妻偷人精品xxx网站 | 午夜久久久久精精品| 黄色成人免费大全| 巨乳人妻的诱惑在线观看| 啪啪无遮挡十八禁网站| 一本精品99久久精品77| ponron亚洲| 亚洲狠狠婷婷综合久久图片| 99久久成人亚洲精品观看| 日韩欧美 国产精品| 91av网站免费观看| 美女高潮的动态| 国产97色在线日韩免费| 老熟妇乱子伦视频在线观看| 18禁国产床啪视频网站| 桃色一区二区三区在线观看| 757午夜福利合集在线观看| 一进一出好大好爽视频| 精品欧美国产一区二区三| 中文字幕高清在线视频| 18禁黄网站禁片免费观看直播| 精品久久久久久久末码| 美女被艹到高潮喷水动态| 日韩欧美精品v在线| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 午夜福利高清视频| 欧美激情在线99| 夜夜躁狠狠躁天天躁| 亚洲专区中文字幕在线| 狠狠狠狠99中文字幕| 亚洲在线观看片| 女人高潮潮喷娇喘18禁视频| 国产av不卡久久| 亚洲成a人片在线一区二区| 又紧又爽又黄一区二区| 90打野战视频偷拍视频| 性欧美人与动物交配| 全区人妻精品视频| 日本免费a在线| 国产伦一二天堂av在线观看| 亚洲欧美一区二区三区黑人| 婷婷亚洲欧美| 国模一区二区三区四区视频 | cao死你这个sao货| 久久久久国产精品人妻aⅴ院| 18禁美女被吸乳视频| 国产精品99久久久久久久久| av视频在线观看入口| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 成人高潮视频无遮挡免费网站| 日本a在线网址| 欧美精品啪啪一区二区三区| 99国产精品一区二区蜜桃av| 三级男女做爰猛烈吃奶摸视频| 亚洲五月天丁香| 亚洲精品一卡2卡三卡4卡5卡| 最近最新免费中文字幕在线| 日本黄色视频三级网站网址| 麻豆国产av国片精品| 我的老师免费观看完整版| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 日本a在线网址| 亚洲国产精品合色在线| 欧美xxxx黑人xx丫x性爽| 欧美日韩综合久久久久久 | 日韩精品中文字幕看吧| 亚洲五月天丁香| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久国产高清桃花| 999精品在线视频| 国产真实乱freesex| 国内少妇人妻偷人精品xxx网站 | 精品国产美女av久久久久小说| 国内久久婷婷六月综合欲色啪| 他把我摸到了高潮在线观看| 国产精品一区二区免费欧美| 国产三级黄色录像| 欧美日韩国产亚洲二区| 99国产精品一区二区蜜桃av| 九九在线视频观看精品| 麻豆成人午夜福利视频| 欧美又色又爽又黄视频| 久久草成人影院| 少妇丰满av| 亚洲精品乱码久久久v下载方式 | 精品欧美国产一区二区三| 99久久久亚洲精品蜜臀av| 熟女人妻精品中文字幕| 国产野战对白在线观看| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 一级毛片高清免费大全| 一个人免费在线观看电影 | 后天国语完整版免费观看| 国产精品综合久久久久久久免费| 两人在一起打扑克的视频| 91av网一区二区| 欧洲精品卡2卡3卡4卡5卡区| 88av欧美| 日韩欧美在线乱码| av视频在线观看入口| 在线观看66精品国产| 两性夫妻黄色片| 18禁观看日本| 亚洲一区高清亚洲精品| 久久精品aⅴ一区二区三区四区| 国产av一区在线观看免费| 在线免费观看的www视频| 99热精品在线国产| 国内毛片毛片毛片毛片毛片| 网址你懂的国产日韩在线| 99riav亚洲国产免费| 亚洲一区二区三区色噜噜| 99久久国产精品久久久| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 亚洲,欧美精品.| 一级毛片精品| 日本免费a在线| 一个人看的www免费观看视频| 美女cb高潮喷水在线观看 | 成熟少妇高潮喷水视频| 亚洲,欧美精品.| 18美女黄网站色大片免费观看| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 在线视频色国产色| 天堂影院成人在线观看| 久久久久久国产a免费观看| 在线观看舔阴道视频| 国产成人精品无人区| 九九在线视频观看精品| 亚洲人成电影免费在线| 日韩成人在线观看一区二区三区| 久久久久精品国产欧美久久久| 亚洲av成人av| 成人永久免费在线观看视频| 免费看十八禁软件| 日韩人妻高清精品专区| 国产97色在线日韩免费| 国产又色又爽无遮挡免费看| 黑人欧美特级aaaaaa片| 嫩草影院精品99| 色老头精品视频在线观看| www.www免费av| 国产爱豆传媒在线观看| 国产成人欧美在线观看| 宅男免费午夜| 少妇的逼水好多| 亚洲欧美日韩卡通动漫| 欧美日韩黄片免| 亚洲人成伊人成综合网2020| 成人三级黄色视频| 欧美中文综合在线视频| 国产精品av视频在线免费观看| 老司机福利观看| 此物有八面人人有两片| 精品国产三级普通话版| 欧美色欧美亚洲另类二区| 1000部很黄的大片| 午夜免费激情av| 欧美不卡视频在线免费观看| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 一区二区三区国产精品乱码| 床上黄色一级片| 午夜精品一区二区三区免费看| 热99在线观看视频| 免费看日本二区| 男女午夜视频在线观看| 午夜激情福利司机影院| 日韩人妻高清精品专区| 午夜激情欧美在线| 亚洲 国产 在线| x7x7x7水蜜桃| 99国产极品粉嫩在线观看| 很黄的视频免费| 日本与韩国留学比较| 成年免费大片在线观看| 亚洲真实伦在线观看| 欧美色欧美亚洲另类二区| 亚洲真实伦在线观看| 亚洲国产精品久久男人天堂| 亚洲中文av在线| 在线十欧美十亚洲十日本专区| 一本久久中文字幕| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 免费观看人在逋| 夜夜看夜夜爽夜夜摸| av片东京热男人的天堂| 国产精品,欧美在线| 亚洲av成人一区二区三| 18美女黄网站色大片免费观看| 亚洲 国产 在线| 欧美中文综合在线视频| 亚洲av日韩精品久久久久久密| 手机成人av网站| 欧美日韩福利视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 国产探花在线观看一区二区| 最近最新中文字幕大全免费视频| 老汉色∧v一级毛片| 99久国产av精品| 老汉色av国产亚洲站长工具| 一级毛片精品| 午夜福利免费观看在线| 亚洲国产中文字幕在线视频| 亚洲专区字幕在线| 最近视频中文字幕2019在线8| 免费一级毛片在线播放高清视频| 日韩三级视频一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美精品v在线| 国产69精品久久久久777片 | 麻豆av在线久日| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 叶爱在线成人免费视频播放| 精品一区二区三区视频在线观看免费| 性欧美人与动物交配| 国产高清视频在线观看网站| 欧美日韩瑟瑟在线播放| 熟女人妻精品中文字幕| 男女做爰动态图高潮gif福利片| 天堂av国产一区二区熟女人妻| 狠狠狠狠99中文字幕| 国产99白浆流出| 国产69精品久久久久777片 | 亚洲自偷自拍图片 自拍| 青草久久国产| 国产真人三级小视频在线观看| 日日夜夜操网爽| 日日干狠狠操夜夜爽| 国产精品综合久久久久久久免费| 国产久久久一区二区三区| 露出奶头的视频| 黑人巨大精品欧美一区二区mp4| 日本黄色片子视频| 国产欧美日韩精品亚洲av| 国产毛片a区久久久久| 熟女人妻精品中文字幕| 他把我摸到了高潮在线观看| 国产伦在线观看视频一区| 久久久色成人| 午夜福利欧美成人| 波多野结衣巨乳人妻| 国模一区二区三区四区视频 | 国产aⅴ精品一区二区三区波| 国产一区二区三区在线臀色熟女| 亚洲国产高清在线一区二区三| 国产v大片淫在线免费观看| 亚洲av成人av| 国产精品一区二区精品视频观看| xxxwww97欧美| 亚洲精品美女久久久久99蜜臀| 欧美日韩瑟瑟在线播放| 狂野欧美白嫩少妇大欣赏| 久久久成人免费电影| 国产伦一二天堂av在线观看| 亚洲国产色片| 国产精品久久久人人做人人爽| 欧美日韩国产亚洲二区| 欧美成狂野欧美在线观看| 日韩 欧美 亚洲 中文字幕| 真人做人爱边吃奶动态| 俺也久久电影网| www日本黄色视频网| 亚洲av中文字字幕乱码综合| 国产成人影院久久av| 巨乳人妻的诱惑在线观看| 小蜜桃在线观看免费完整版高清| 日本一二三区视频观看| 午夜福利欧美成人| netflix在线观看网站| tocl精华| 九九久久精品国产亚洲av麻豆 | 长腿黑丝高跟| 国产精品一区二区精品视频观看| 国产亚洲欧美98| 变态另类丝袜制服| 黑人巨大精品欧美一区二区mp4| 久久草成人影院| 久久久久久大精品| 国产99白浆流出| 欧美国产日韩亚洲一区| 香蕉丝袜av| 1000部很黄的大片| 精品久久久久久久久久久久久| 亚洲欧美一区二区三区黑人| 此物有八面人人有两片| 久久久精品大字幕| 麻豆av在线久日|