Community Characteristics Analysis of Eukaryotic Microplankton Via ITS Gene Metabarcoding Based on Environmental DNA in Lower Reaches of Qiantang River, China
Total Page:16
File Type:pdf, Size:1020Kb
Open Journal of Animal Sciences, 2021, 11, 105-124 https://www.scirp.org/journal/ojas ISSN Online: 2161-7627 ISSN Print: 2161-7597 Community Characteristics Analysis of Eukaryotic Microplankton via ITS Gene Metabarcoding Based on Environmental DNA in Lower Reaches of Qiantang River, China Aiju Zhang1, Jun Wang1, Yabin Hao1, Shanshi Xiao1, Wei Luo1, Ganxiang Wang2, Zhiming Zhou1 1Agriculture Ministry Key Laboratory of Healthy Freshwater Aquaculture, Key Laboratory of Freshwater Aquaculture Genetic and Breeding of Zhejiang Province, Zhejiang Research Center of East China Sea Fishery Research Institute, Zhejiang Institute of Freshwater Fisheries, Huzhou, China 2Pinghu Fisheries Technology Promotion Center, Pinghu, China How to cite this paper: Zhang, A.J., Wang, Abstract J., Hao, Y.B., Xiao, S.S., Luo, W., Wang, G.X. and Zhou, Z.M. (2021) Community Eukaryotic microplankton plays an important role in water biotic community Characteristics Analysis of Eukaryotic and in maintaining the stability of water ecosystems. Environmental DNA Microplankton via ITS Gene Metabarcod- metabarcoding provides the opportunity to integrate traditional and emerg- ing Based on Environmental DNA in Low- er Reaches of Qiantang River, China. Open ing approaches to discover more new species, and develop molecular biotic Journal of Animal Sciences, 11, 105-124. indices that can be more rapidly, frequently, and robustly used in water qual- https://doi.org/10.4236/ojas.2021.112009 ity assessments. In order to examine assemblages of eukaryotic microplank- ton in lower reaches of Qiantang River, ITS gene metabarcoding technology Received: January 13, 2021 based on environmental DNA was carried out. As a result, various species of Accepted: March 30, 2021 Published: April 2, 2021 phytoplankton, fungi and zooplankton were annotated on. More phylum, classes and specieses of eukaryotic phytoplankton and zooplankton were Copyright © 2021 by author(s) and found after compared communities taxa based on metabarcoding with that Scientific Research Publishing Inc. obtained from morphological examination. Nevertheless, Chlorophyceae was This work is licensed under the Creative Commons Attribution International the most common assemblage both identified by using these two methods, License (CC BY 4.0). also Mesocyclops leuckarti and Acanthocyclops bicuspidatus were both found http://creativecommons.org/licenses/by/4.0/ to be the dominant species of Cyclopoida in the river. Additionally, the reads Open Access proportions of phytoplankton and zooplankton at the three freshwater sam- pling sites (Tonglu, Fuyang and Wenyan) decreased as temperature drop. Meanwhile, twenty classes of fungi were annotated on, of which the commu- nity characteristic was first researched in the river. There were significant spatial differences in values of Chao1 index for eukaryotic microplankton. Cluster analysis and Non-metric multidimensional scaling ordination further DOI: 10.4236/ojas.2021.112009 Apr. 2, 2021 105 Open Journal of Animal Sciences A. J. Zhang et al. confirmed that the community composition of eukaryotic microplankton at class level for Jiashao-September sample had the most dissimilarity with the others. Keywords Microplankton, eDNA, Metabarcoding, ITS, Diversity, Qiantang River 1. Introduction Eukaryotic microplankton is a group of plankton with particle sizes less than 20-µm [1]. It plays an important role in water biotic community and in main- taining the stability of water ecosystems due to its important link in energy flow and material circulation [1]. As the basis of the food chain, analyzing its meta- community structure is very important to assess the status and development tendency of an ecosystem [2]. Traditional practices for biological surveys of inland waters usually center on a common set of ecological indicators or indices/measurements of biodiversity [3], requiring many preparations, such as morphological taxonomic expertise, intact specimens and adequate time [4] [5]. Environmental DNA (eDNA) meta- barcoding provides the opportunity to integrate traditional and emerging ap- proaches to discover more new species [6] [7], and develop molecular biotic in- dices that can be more rapidly, frequently, and robustly used in water quality as- sessments [6] [8] [9]. So far, metabarcoding technology has been widely used for biodiversity monitoring in biological and environmental samples [10] [11] [12] [13]. Despite some limitations (each marker region might import some biases and the blast sequence database is incomplete), the practices of metabarcod- ing-based analysis for estimating diversity and relative abundance of taxonomic groups in aquatic systems will likely increase as technology improved [14]. Now, marker genes used for the analysis of plankton communities via eDNA metabarcoding often focuses on ribosomal DNA (rDNA), ribulose bisphosphate carboxylase large subunit gene (rbc L) and cytochrome C oxidase subunit gene 1 (cox1) [8] [15] [16] [17] [18]. Among these, the internal transcribed spacer (ITS), as one gene fragment of rDNA, has been widely used due to its fast evolu- tion and high specificity [17], and its relevant metabarcoding practices for the purpose of monitoring diversity of eukaryon communities mainly focused on soil, plant and marine systems [19] [20] [21] [22]. Here, an ITS metabarcoding assessing assemblages of eukaryotic microplank- ton were conducted. Specifically, ITS gene sequence analysis was performed on surface water samples collected from four sites in Qiantang River lower reaches (Zhejiang, China), an important freshwater fishing water used for drinking, electricity generation, flood control and recreation, in order to analyze the community diversity of eukaryotic microplankton in the section, to assess the DOI: 10.4236/ojas.2021.112009 106 Open Journal of Animal Sciences A. J. Zhang et al. utility of this approach for monitoring diversity of freshwater eukaryotic micro- plankton community. As a whole, this study provided a theoretical basis for fur- ther study on the function of eukaryotic microplankton in freshwater ecosystems. 2. Materials and Methods 2.1. Water Sample Collection and Physic Chemical Analysis A total of eight water samples were collected in September and November 2019, at four separate sampling sites in Qiantang River lower reaches, including Ton- glu, Fuyang, Wenyan and Jiashao sites, hereafter referred to as TL, FY, WY and JS separately (Figure 1). All sampling, filtering, and other equipments were ste- rilized before use. 2 L of surface water was collected at a set time by boat from each site. Simul- taneously, five environmental variables, such as transparency, salinity, water temperature (WT), pH and dissolved oxygen (DO), were measured in situ by using a secchi disk and a portable water quality detector (Hach, USA). For each (a) (b) Figure 1. Map shows the study site. (a) The location of Qiantang River in Zhejiang prov- ince of China; (b) Sampling sites in Qiantang River lower reaches. The red dot represents the water sampling site, the black dot represents the local government seat. DOI: 10.4236/ojas.2021.112009 107 Open Journal of Animal Sciences A. J. Zhang et al. water sample, a 1.5 L water subsample was used to test six environmental va- riables immediately after transported to the laboratory within 24 hours at 4˚C, including Chemical oxygen demand (COD), Chlorophyll a (Chla), total nitrogen (TN), total phosphorus (TP), ammonium (NH4-N) and nitrate-nitrite (NO3-N), the recording of which were following the scheme of [23]. All environmental va- riables were measured in triplicate. 2.2. DNA Sample Processing and High-Throughput Sequencing (HTS) For each water sample, a residual 500-mL water subsample was then filtered through a 0.22-μm cellulose acetate filter paper using a peristaltic pump in the field. Then, each paper was placed inside a commercial sterile centrifuge tube and stored in a container filled with liquid nitrogen until subsequent manipula- tions were performed. DNA was extracted from filters using EZNA water DNA kit (Omega, USA) following the manufacturer’s protocol. The concentration and purity of DNA were determined using NanoDrop 2000c spectrophotometer (Thermo, USA), then followed by multiplex PCR using the universal primers for ITS belonging to eukaryotic mitochondrial DNA fragments, ITS-F (5'-GTGA ATCATCGARTC-3), ITS-R (5'-TCCTCCGCTTATTGAT-3') [22]. Each eDNA sample was subjected to PCRs in triplicate [17.35 μL dH2O, 2.5 μL 10 × PCR buffer, 0.15 μL 5 U/μL Thermo scientific Taq DNA polymerase, 2 μL dNTPs mix (2.5 mM), 1 μL of each primer (10 mM), and 1 μL sample eDNA in a total vo- lume of 25 μL]. PCR cycling parameters were as follows: 96˚C for 3 min, fol- lowed by 35 cycles of 96˚C for 30 s, 50˚C for 30 s, and 72˚C for 90 s, with a final elongation step at 72˚C continued for 7 min. Negative control was conducted simultaneously. After subjecting the PCR products to 1% polyacrylamide gels (see Figure S1), the quantified, size-selected libraries were constructed and con- tinuously diluted to a concentration suitable for sequencing. The libraries were finally sequenced on the Illumina MiSeq 2000 platform by following the manu- facturer’s protocols step by step. 2.3. Phytoplankton Samples Collection and Treatment Phytoplankton samples were also simultaneously collected at the four sites. For phytoplankton counts, 1.0 L of water samples were sampled each time and pre- served with 1% Lugol’s iodine solution. Phytoplankton samples were concen- trated to a final volume of 30 ml after sedimentation for 48 h. Thereafter, the taxa were verified and counted under 200× and 400× magnifications for at least 500 specimens [23]. The data were made to compare with that collected from ITS gene metabarcoding method. 2.4. Bioinformatics and Sequencing Data Upload The raw sequencing FASTQ file was transformed to a FASTA file by the Fastx toolkit V0.0.1 [24]. Clean reads were gained after trimming the low quality se- DOI: 10.4236/ojas.2021.112009 108 Open Journal of Animal Sciences A. J. Zhang et al. quences and PCR chimeras by using Fqtrim V0.9.4 (http://ccb.jhu.edu/software/fqtrim/) and Vsearch 2.3.4 [25].