Relationship between morphospecies and microcystinproducing genotypes of Microcystis species in Chinese freshwaters*

Yang LIU ,Youxin CHEN ,Haiyan FANG ,Hanyang LU ,Xingqiang WU ,Gongliang YU ,Shin-ichi NAKANO,Renhui LI

1 College of Life Science, Henan Normal University, Xinxiang 453007, China

2 Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China

3 Center for Ecological Research, Kyoto University, Shiga 520-2113, Japan

4 College of Life and Environmental Sciences, Wenzhou University, Wenzhou 325000, China

5 Journal of Henan Normal University, Xinxiang 453007, China

Abstract Twenty water bodies in China were sampled,and 186 strains of different Microcystis species were isolated,from which eight morphospecies were identified and 43 stains containing the mcyB gene were detected.Phylogenetic analysis based on the mcyB gene indicated that the microcystin (MC)-producing Microcystis in China could be divided into two groups (I and II) and showed significant differences between the two groups.The maximum sequence similarity was 69.1%.Microcystins (MCs) were measured by high-performance liquid chromatography (HPLC) analysis,and no microcystin-RR (MC-RR) was detected in some strains belonging to Group II.Compared to other regions of the world,the proportion of Chinese MC-producing Microcystis was different,and the regional differences were more obvious.A whole-cell polymerase chain reactio (PCR) assay was conducted to analyze the proportion of the mcyB gene in the laboratory cultured and field cultured Microcystis.The proportion of four morphospecies (M.vividis,M.ichthyoblabe,M.novacekii,and M.aeruginosa) that contained the mcyB gene exceeded 50% in the field cultured samples.Compared with former studies,M.aeruginosa was the most likely morphotype that can produce MCs in the world.This study provided new insight of Microcystis hazard assessment and field monitoring.

Keyword:Microcystis;morphospecies; mcyB;16S rDNA;microcystin;high-performance liquid chromatography (HPLC)

1 INTRODUCTION

Cyanobacterial blooms are frequently observed in eutrophic or estuarine waters worldwide.Several genera of cyanobacteria,includingMicrocystis,Anabaena,Planktothrix,Aphanizomenon,andNostoc,can produce a variety of secondary metabolites such as hepatotoxins,neurotoxins,and cytotoxins(Liu et al.,2016;Shan et al.,2019;Wang et al.,2019,2021).Among cyanobacteria,Microcystisis considered the most prominent contributor to the production of microcystin (MC),which are potent hepatotoxins.More than 90 MC variants have been reported so far (Hotto et al.,2008;Attard et al.,2018).These toxins could potentially inhibit the activities of eukaryotic serine/threonine protein phosphatase 1 and 2A,and have been implicated in several livestock and human poisonings (Mikalsen et al.,2003;Tanabe et al.,2004).Thus,the World Health Organization has published a preliminary guideline value for microcystin-LR (MC-LR) in drinking water,which is 1 μg/L(WHO,1998;Wu et al.,2007;Otten and Paerl,2011).

MCs are produced by a mixed polyketide synthase(PKS)/nonribosomal peptide synthetase (NRPS)complex via a thiotemplate mechanism (Mikalsen et al.,2003).MC synthetase (mcy) gene clusters encoding these biosynthetic enzymes have been identified and sequenced from the unicellular organismMicrocystisaeruginosa(Rantala et al.,2004).Themcygene cluster is composed of 10 genes arranged bidirectionally in two operons (mcyA-CandmcyD-J) (Rantala et al.,2004;Cirés et al.,2013).ThemcyA-Jgene cluster is essential for the production of MCs byMicrocystis.Therefore,molecular probes have been developed to detect and study the community composition and dynamics of MCproducing and non-MC-producingMicrocystisstrains in the field (Kurmayer et al.,2002;Via-Ordorika et al.,2004).ThemcyBgene,which encodes a peptide synthetase of 242 kDa,has been widely used as a marker for toxin biosynthesis and genetic diversity studies (Mikalsen et al.,2003;Dyble,2008).In contrast,a small number ofMicrocystisstrains tested positive for themcygene but lacked detectable MCs.Insertional inactivation,homology of themcyBregion or even the presence of themcygene cluster in some nontoxicMicrocystisstrains may explain these discrepancies (Via-Ordorika et al.,2004).Therefore,molecular methods combining chemical detection are necessary.These methods have been widely used to detect MC-producingMicrocystispopulations.

Likewise,these methods are helpful for the rapid detection of the toxicMicrocystisin freshwater bodies,and also provide a support for the studies on the relationships among the morphology,colony size,and the occurrence ofmcygene-containing colonies.Relationship betweenMicrocystiscolony size and the occurrence ofmcygene-containing colonies was studied earlier (Via-Ordorika et al.,2004;Dyble et al.,2008).By far,a consensus is that the largest size colony has a maximum proportion ofmcy(Jungmann et al.,1996;Via-Ordorika et al.,2004).A new study of the occurrence of MCs biosynthesis genes in natural populations of differentMicrocystisspecies in Russian reservoirs (Sidelev et al.,2020) show that no correlation between colony size and the frequency ofmcygenes for individual morphospeciesM.aeruginosaandM.flos-aquae.However,they identified only five morphotypes ofMicrocystisand two of them were chosen in 10 reservoirs.At present,at least 51 species ofMicrocystishave been reported in the world,and approximately 13 species have been described in China (Otsuka et al.,1999a;Wu et al.,2007;Xu et al.,2008;Xiao et al.,2017;Pérez-Carrascal et al.,2019).Therefore,more morphospecies and water bodies are needed to be used to investigate and analyze the morphology and toxigenic capacity ofMicrocystisspecies.At present,Microcystisare dominant species in most freshwater bodies of China,such as Taihu Lake,Caohu Lake,and Dianchi Lake (Xu et al.,2008;Shen et al.,2018).Most species of the genusMicrocystishave the ability to produce MCs,and a quick judgment to identify potential toxicMicrocystisblooms by morphology is also needed for water quality management.In this study,20 water bodies (in 13 cities) across several climatic zones of China were surveyed to examine the relationship between morphospecies and MC-producing genotypes of the genusMicrocystis.Furthermore,themcygenes and toxin content were also analyzed by Polymerase Chain Reactio (PCR) and high-performance liquid chromatography (HPLC).

2 MATERIAL AND METHOD

2.1 Sampling and isolation

Samples were collected with a plankton net (25#,200 mesh) from 20 freshwater bodies,located in 13 cities of China (Fig.1) from tropical to temperate zones during 2006 and 2007.There are two types of samples.One is from field experiments and the other is from laboratory experiments.In the field experiment,a 500-mL sample was collected from each water body,and all the types ofMicrocystisfrom each water body were distinguished by classical taxonomic identification.Fifty samples of eachMicrocystistype were selected for whole-cell detection.In the laboratory experiment,single colonies ofMicrocystisfrom each water body were isolated using the Pasteur capillary pipette method under a dissecting microscope,and each colony that was chosen was washed three times in sterile water,and each colony was kept in one well of a 24-well plate in 2-mL MA culture medium (Ichimura,1979).To successfully isolate single colonies ofMicrocysits,we picked up 24 colonies in each freshwater body based on our experience.All cultures were grown at(25±1) °C with a 12-h∶12-h light∶dark cycle under a photon irradiance of 25 μmol photons/(m2·s) provided by daylight fluorescent lamps.Taxonomic identifications were assigned following the descriptions of Komárek and Anagnostidis (1998) by a light microscope (Olympus BX51,Japan).

Fig.1 Map of the sampling sites

2.2 DNA extraction,cloning and sequencing of 16S rRNA,mcyB

Genomic DNAs (the pure culture grew to a certain biomass) were extracted following the modified cetyltrimethylammonium bromide (CTAB) method as described by Meißner et al.(1996).The primers for amplification of theMicrocystis16S rDNA fragment were MicroF (5ʹ-GCCGCRAGGT-GAAAMCTAA-3ʹ)and MicroR (5ʹ-AATCCAAAG-ACCTTCCTCCC-3ʹ)(Neilan et al.,1997).To identify the micorcystinproducing ability ofMicrocystisstrains,we detected a region within themcyBgene required for microcystin biosynthesis by PCR.The primer sequences were tox4f (5ʹ-GGATATCCTC-TCAGATTCGG-3ʹ) and tox4r (5ʹ-CACTAACCC-CTATTTTGGATACC-3ʹ)(Kurmayer et al.,2002),which have been suggested to be specific for the genusMicrocystis.The PCRs were carried out in 50-μL mixtures containing 5–10 ng of genomic DNA,25 μL of 2×PCR Master Mix,10 pmol of primer and sterile water to 50 μL.The reactions were run in a BioRad Thermal Cycler(BioRad,USA) with one cycle at 94 °C for 5 min,35 cycles at 94 °C for 30 s,annealing at 50 °C for 30 s(MicroF/MicroR) or at 55 °C for 30 s (tox4f/tox4r),strand extension at 72 °C for 25 s,and final extension at 72 °C for 5 min.PCR products were detected by 1% agarose gel electrophoresis,and amplified bands were detected under 300 nm ultraviolet light.Then,PCR products were purified with the kits (Qiagen,Germany) according to the manufacturer’s protocol.Then,the purified products were cloned into the PMD18-T vector (TaKaRa,Japan) and transformed intoE.coliDH5α competent cells (TaKaRa,Japan)on an ice bath.Finally,a single positive clone was selected for sequencing (Invitrogen,Shanghai,China)for each isolate.The accession numbers for the obtainedmcyBsequences are KJ818165–KJ818201,and for 16S rDNA gene sequences are KJ818122-KJ818163.

2.3 Sequence alignment and phylogenetic analysis

Sequences were initially corrected and aligned by using BioEdit software (http://www.mbio.ncsu.edu/bioedit/bioedit.html).Twenty-one strains ofMicrocystiscovering 6 different species were downloaded from GenBank (www.blast.ncbi.nlm.nih.gov),and 37 sequences obtained from this study were used for the 16S rDNA phylogenetic analysis.The 16S rDNA sequence ofSynechocystisPCC6803 was used as the out-group.Phylogenetic distances in the trees were estimated using the algorithm of Kimura (1980),and 1 000 bootstrap trials were performed.The phylogenetic tree was constructed using the neighbor-joining (NJ) method in MEGA X(Kumar et al.,2018).Phylogenetic analysis formcyBwas performed the same as above.A total of 67 sequences were used to construct the phylogenetic tree,including 24 sequences from GenBank,and 42 sequences from the present study,andNostocsp.was used as the outgroup.Finally,the intra-group and between group similarity were performed by using“sequence identity matrix”option in BioEdit software.

2.4 HPLC analysis

Strains ofMicrocystisthat have been shown to have themcyBfragment were large-scale cultured and collected,as previously described by Liu et al.(2011).TheMicrocystiscells were lyophilized,and approximately 50 mg of these freeze-dried cells were extracted for the analysis of MCs.Chromatographic separation was performed using a Waters Alliance 2695 HPLC system,Waters 2996 photo diode array detector and Waters empower chromatography software (Waters,USA).The wavelength was 238 nm,and the chromatographic column was a Synergi Hydro-RP C18 column (4 μm,250 mm×4.6 mm).The mobile phase consisted of pure methanol (solvent A)and 0.05% aqueous trifluoroacetic acid (solvent B).The linear gradient was 50% solvent A at 0 min,80%solvent A at 15 min,50% solvent A at 16 min,and 50% solvent A at 20 min.The detector resolution was set at 1.2 nm and the flow rate was set at 1 mL/min with a column temperature of 30 °C.

2.5 Whole cell PCR in field samples

A whole cell PCR method was performed to increase the effi ciency of detecting MC from field samples (Howitt,1996;Pan et al.,2001).Microcysitscolony for the whole-cell PCR was isolated from each sample by using the Pasteur capillary pipette method.The isolates were washed one to three times with distilled water and rechecked under an inverted microscope.Finally,1–2-μL isolates as the DNA template were directly used in PCR amplification to detect themcyBgene.The method used was the same as that mentioned in Section 2.2.Finally,50 samples were collected and detected to evaluate the ratio of MC producing strains in the field.

3 RESULT

3.1 Morphospecies of 186 Microcystis strains

A total of 186 strains ofMicrocystiswere isolated from 20 water bodies in China.The strains were distinguished by classical taxonomic identification,since it is still an irreplaceable method.It is necessary to identify the strains as soon as they were collected due to their high phenotypic plasticity.Finally,eight morphospecies were identified according to the descriptions of Komárek and Komárková (2002)(Fig.2):M.aeruginosa(Kuëtzing) Kuëtzing,M.botrysTeiling,M.flos-aquae(Wittrock) Kirchner,M.novacekii(Komaèrek) Compère,M.ichthyoblabeKuëtzing,M.wesenbergii(Komárek) Komárek,M.viridis(A.Brown) Lemmermann,andM.smithiiKomárek &Anagnostidis. The number of each morphospecies was shown in Table 1.

Table 1 Information of Microcystis strains used in this study and the results of HPLC analysis

Fig.2 Colonies of eight Microcystis morphospecies in China

3.2 Phylogenetic analysis

From the 186 isolates,only 43 showed the presence ofmcyBby PCR.The 16S rRNA gene sequences of the 43Microcystisstrains withmcyBgenes were amplified.Each sequence showed a product of 1 444 bp,which corresponded to the length of the 16S rDNA,and an alignment of 1 402 nucleotides after excluding positions with gaps was used.The NJ phylogenetic tree based on the 16S rRNA gene (Fig.3),includes all 43 strains ofMicrocystisin this study and 21 strains ofMicrocystisobtained from GenBank.Amplification of themcyBgene with primers tox4f and tox4r revealed that 43 strains ofMicrocystisgave PCR products in size of 1 312 bp and an alignment of 1 116 nucleotides after excluding positions with gaps was used.The NJ phylogenetic tree based on themcyBgene (Fig.4)includes all 43 strains ofMicrocystisin this study and 24 strains ofMicrocystisobtained from GenBank.Finally,the sequence similarities in Goup I and Group II (Fig.4) were calculated separately.The high similarity within each group was 99%,and the maximum similarity ofmcyBsequences between the two groups was 69.1%.

Fig.3 Phylogenetic tree (NJ) based on 16S rRNA sequences

Fig.4 Phylogenetic tree (NJ) based on mcyB sequences

3.3 HPLC analysis of MC occurrence in different morphospecies

Analysis of themcyBgene from 186 strains ofMicrocystisby PCR showed that only 43 strains produced a PCR product formcyB.Then,strains of the same region and morphology were removed (strains from the same place and sampling site were grouped,and 1–2 isolates were randomly selected),and 20 strains were finally enlarged in large-scale culture to analyze the type and content of microcystin in different morphospecies.The results (Table 1) show that 7 morphospecies sampled from 10 different waterbodies contained different microcystin variants and contents.All testedM.viridiscontained two types of microcystin(MC-LR and MC-RR) except the strain from Fuxian Lake.The contents of MC-RR were much higher inM.viridisthan MC-LR.AllM.flos-aquaestrains contained MC-LR but no MC-RR was observed.

3.4 The proportion of Microcystis containing the mcyB gene

IndividualMicrocystiscolonies,consisting of 8 different morphospecies from indoor culture and field samples,were tested for themcyBgene using a whole-cell PCR assay.Thirty individualMicrocystiscolonies of each morphospecies were used to detect themcyBgene.The results are shown in Table 2.All of theM.wesenbergiimorphospecies produced no MC,which is consistent with previous studies (Xu et al.,2008).In our survey,the proportion of four morphospecies (M.viridis,M.ichthyoblabe,M.novacekii,andM.aeruginosa) that contained themcyBgene exceeded 50% in the field samples.Accordingly,these species are most likely to produce MCs in Chinese water bodies.

Table 2 The percentage of microcystin-producing Microcystis spp.in the samples

4 DISCUSSION

4.1 Morphospecies of MCs producing strains

Although cyanobacterial blooms have broken out all over the world,the compositions ofMicrocystisin different regions were varied.Six morphospecies ofMicrocystiswere distinguished in Japanese waters(Otsuka et al.,1999b),five morphotypes ofMicrocystiswere identified in Russia (Sidelev et al.,2020),and seven morphospecies were found in European waters(Kurmayer et al.,2002;Via-Ordorika et al.,2004).The eight species identified in this study were strains common in most water bodies in China and could produce MCs,except forM.wesenbergii.This result is consistent with previous studies in that nomcyBfragment was detected (Xu et al.,2008).Although the seven species were able to produce MCs,some of them may lack detectable MCs (Tanabe et al.,2004;Via-Ordorika et al.,2004).It was still unclear to what extend the morphological characteristics were linked to the ability of producing MCs.However,the colonies assigned toM.aeruginosawere identified as having the highest frequency of toxin production(Kurmayer et al.,2002;Via-Ordorika et al.,2004).A previous research found that toxicity ofMicrocystisunder natural conditions could be associated mainly with large-diameter colonies (Jungmann et al.,1996),which may be consistent with our finding thatM.aeruginosaalways had a large size in the field.

4.2 Phylogenetic analysis

The phylogenetic tree (Fig.3) shows low bootstrap values as a result of the extremely high similarity(99.2%–100%) of sequences determined.The species ofMicrocystisisolated from different places were mixed together according to the topology tree.Moreover,different morphotype strains were also mixed;for example,M.botrys,M.viridis,andM.wesenbergiiwere mixed with the“M.aeruginosacomplex”(Otsuka et al.,1998,1999a;Tan et al.,2010).Although the 16S rDNA sequence is insuffi ciently variable for phylogenetic analysis ofMicrocystisspp.at the species level,it is still valid at the genus level.Therefore,this may appeal to the findings of population genomics that predicatedMicrocystisrepresents a single,globally distributed,and homogeneous gene pool (Pérez-Carrascal et al.,2019).However,our study agreed with the view that more“ecophenic”and/or phenotypic forms have been described than genospecies exist (Palinska et al.,1996;Le Ai Nguyen et al.,2012).

The phylogenetic tree (Fig.4) could be divided into two groups (Groups I &II),and high similarity (99%)of sequences within each group was found.The maximum similarity ofmcyBsequences between the two groups was 69.1%.Group I covered more of the sequences in GenBank than Group II.A large proportion of species in Group II were uncultured cyanobacteria.The three strains ofM.flos-aquaein Group II (CHAB543,CHAB546,and CHAB547)contained one type of MC-LR but no MC-RR.ThemcyBregion was more variable in the MC biosynthesis gene cluster and showed high mutation rates(Kurmayer et al.,2002;Mikalsen et al.,2003).This was related to the production of various MC isoforms.We use the same sequences from GenBank (AY034602 and AB092806) that covered the MC synthetase module encoding themcyB1domain to perform phylogenetic analysis (Dyble et al.2008).The results of this study confirm thatmcyBsequences inMicrocystiscould be grouped into two major clusters(mcyB1(B) andmcyB1(C)) (Mikalsen et al.,2003).The (B) variants ofmcyB1produced various MC-LR isoforms.The (C) variant of themcyB1module synthesized MC-RR,and the members of some subgroups produced MC-LR in addition to MC-RR(Mikalsen et al.,2003;Dyble et al.,2008;Otten and Paerl,2011).These findings were very interesting as Group I and Group II havemcyB1(C) andmcyB1(B),respectively,each group with the same properties to produce MC isoforms.Such findings suggested that the genotypic variations within themcyBmodule were strongly correlated with their isoforms of MC.

4.3 Analysis of the mcyB gene and MCs

ThemcyBsequences of theseM.flos-aquaewere also unique and had a low similarity with others in this study.Previous studies have revealed several phylogenetic subgroups of the adenylation domain in themcyBgene.Mikalsen et al.(2003) named the two major clusters ofmcyBsequences inMicrocystismcyB1(B) andmcyB1(C).Strains that synthesized MC-RR generally belonged to the C-type subgroup,and the B-type subgroup generally produced MC-LR.This was the reason why these strains were divided into two groups in the phylogenetic tree.Some strains from GenBank could group with theseM.flos-aquaestrains,and whether they contain no MC-RR was unknown.The contents of MC-RR were much higher than that of MC-LR in Group I.These findings also revealed that geographical segregation existed on the basis of the homology ofmcyBsequences.These conclusions generally supported the results of previous studies showing thatMicrocystisstrains may have geographical distribution in the world (Wu et al.,2007).In our study,all strains ofM.aeruginosahad both MC-LR and MC-RR.TheM.botrysfrom different sites contained distinct MCs.It was also found that evenMicrocystisfrom the same place contained distinct MCs.Again,this study supported the results of previous studies showing that genetic variation in themcyB1module generated different MC isoforms (Mikalsen et al.,2003).This is the first to report two strains ofM.novacekiiproducing both MC-LR and MC-RR that were isolated in China.Previous studies found that noM.novacekiistrains isolated in Japan produced MC,but two toxic strains ofM.novacekiiwere isolated in Thailand (Otsuka et al.,1999a,2000;Tsujimura et al.,2000).Likewise,we found thatM.smithiiproduced MC-RR,which we covered in previous article (Liu et al.,2011).Although theM.ichthyoblabemorphotype in Taihu Lake consistently lacked themcyBgene,all four strains possessed themcyBgene and produced both MC-LR and MC-RR except one strain.This strain (CHAB181)contained only MC-RR,but its MC content was the highest,up to 2 825.44 μg/g.Previously studies reported little about the characteristics of toxin production byM.ichthyoblabe,and the main reason may be that it was diffi cult to distinguish these four close morphospecies (M.aeruginosa,M.flos-aquae,M.ichthyoblabe,andM.novacekii) (Otsuka et al.,1999b;Tan et al.,2010).These four species were presumed to be“M.aeruginosacomplex”.M.aeruginosahas been extensively studied.All the strains in Table 1 are consistent with B-type variants ofmcyB1producing various MC-LR isoforms and most C-type variants of themcyB1module synthesizing MC-RR comparing with previous studies (Mikalsen et al.,2003;Dyble et al.,2008).

4.4 The proportion of the mcyB gene in cultured and field Microcystis in Chinese freshwaters

Kurmayer et al.(2002) reported that in Lake Wannsee of Germany,the proportion ofMicrocystiscolonies with themcyBgene were as follows:73%M.aeruginosa,16%M.ichthyoblabe,37%M.flosaquaeand 45%M.botrys.It was interesting that the proportion of toxicM.botrysandM.aeruginosawas almost the same.Another study in European freshwater bodies also found that more than 75% of theM.botrysandM.aeruginosacontained themcygenes,and approximately 20% of theM.ichthyoblabeandM.vividisgave a PCR product of themcygenes(Via-Ordorika et al.,2004).Although the proportion of toxicM.ichthyoblabeandM.vividiswas different between China and Europe,M.aeruginosawas the most likely morphotype to produce MC (Vasconcelos and Pereira,2001;Cirés et al.,2013).In this study,the proportions ofM.flos-aquaeandM.smithiiin the field were 20% and 11%,respectively.We also studied the proportion of cultured species,and the results from the indoor cultures are slightly different from those of the field samples.We inferred that it might be associated with the culture conditions because temperature,light,and nutrients had some effects on the production of MCs (Song et al.,1998;Wiedner et al.,2003).Moreover,mostMicrocystismorphotypes could maintain their colony characteristics for approximately one month (Otsuka et al.,2000).Accordingly,it seemed to be obscure to find distinction among these morphospecies.This may support the idea of population genomics thatMicrocystispopulations are not a coherent species (Pérez-Carrascal et al.,2019).Inconsistencies between traditional and population genomic classifications exist and need a further study.The strains ofMicrocystisused in this study were not sterile.Therefore,this may be another reason why the two results are different.However,our findings still provide new insight intoMicrocystishazard assessment and field monitoring.

5 CONCLUSION

The results show a relationship between morphospecies and MC-producing genotypes of commonMicrocystisspecies in Chinese freshwaters.The data also show thatmcyBgene-containingMicrocystiscolonies of China could be divided into two groups (Group I and Group II).MostMicrocystisstrains of Group I contained two types of MC isoforms,while Group II contained only MC-LR.The proportion of four morphospecies (M.vividis,M.ichthyoblabe,M.novacekii,andM.aeruginosa)that contained themcyBgene exceeded 50% in the field samples.It was consistent with previous studies thatM.aeruginosawas the most likely morphotype to produce MC in the world.It could help water quality managers for a quick judgment on potential toxicMicrocystisblooms.

6 DATA AVAILABILITY STATEMENT

The data generated in this study are available from the corresponding author on reasonable request.

7 ACKNOWLEDGMENT

We are grateful to the lab members from Harmful Alga Group for assistances with samples collection and strains isolation.We are also thankful to Dr.Shen LIN,Li YANG,Zhongjie WANG,Mengling ZHU,and Prof.Zhongxing WU for your help in this research.