Ectothiorhodospira lacustris sp. nov., a New Purple Sulfur Bacterium from Low-Mineralized Soda Lakes That Contains a Unique Pathway for Nitric Oxide Reduction

Several strains of a Gram-negative, anaerobic photoautotrophic, motile, rod-shaped bacterium, designated as B14B, A-7R, and A-7Y were isolated from biofilms of low-mineralized soda lakes in central Mongolia and Russia (southeast Siberia). They had lamellar stacks as photosynthetic structures and bacteriochlorophyll a as the major photosynthetic pigment. The strains were found to grow at 25–35 °C, pH 7.5–10.2 (optimum, pH 9.0), and with 0–8% (w/v) NaCl (optimum, 0%). In the presence of sulfide and bicarbonate, acetate, butyrate, yeast extract, lactate, malate, pyruvate, succinate, and fumarate promoted growth. The DNA G + C content was 62.9–63.0 mol%. While the 16S rRNA gene sequences confirmed that the new strains belonged to the genus Ectothiorhodospira of the Ectothiorhodospiraceae, comparison of the genome nucleotide sequences of strains B14B, A-7R, and A-7Y revealed that the new isolates were remote from all described Ectothiorhodospira species both in dDDH (19.7–38.8%) and in ANI (75.0–89.4%). The new strains are also genetically differentiated by the presence of a nitric oxide reduction pathway that is lacking from all other Ectiothiorhodospiraceae. We propose to assign the isolates to the new species, Ectothiorhodospira lacustris sp. nov., with the type strain B14BT (=DSM 116064T = KCTC 25542T = UQM 41491T).

A large number of genome sequences of different species and strains of Ectothiorhodospiraceae isolated from different habitats at different times are now available [1,9]. The genome sequences of strains B14B, A-7R, and A-7Y, isolated in the late 20th century from low-salinity soda lakes in central Mongolia and southeastern Siberia (Russia), were recently analyzed to establish their taxonomic position [1,10]. Based on a comparison of genome length and G + C content (mol%) in DNA, they belong to a group of species related to Ect. shaposhnikovii and Ect. vacuolata. On the phylogenetic tree of total genomes, strains B14B, A-7R, and A-7Y form a separate branch and have ANI values <90% (75.0-89.4%) relative to any of the other strains of the genus Ectothiorhodospira, but high ANI values between themselves (99.0-99.9%). Therefore, they are considered to represent a new species [1].
Although these strains were isolated over two decades ago, their specific characteristics had not yet been described. Given their unique genomic placement, a more detailed characterization of these strains was warranted. This paper presents morphological and physiological characterization of strains B14B, A-7R, and A-7Y and describes the taxonomic position of the new species. Strain B14B is designated as the type of strain.
Pure cultures were obtained by serial dilutions after repeated transfers of well-isolated colonies grown in 0.5-0.7% (w/v) agar medium. In liquid media, the cultures were grown anaerobically in screw-capped glass bottles at 25-35 • C and illumination of 2000 lx.

Morphological, Physiological, and Biochemical Characteristics
Cell morphology was studied using an Olympus BX 41 phase contrast microscope (Olympus, Tokyo, Japan). Whole cells and ultrathin sections were examined under a JEOL JEM-100C electron microscope (JEOL, Tokyo, Japan) at 80 kV. Whole cells were contrasted with 1% phosphotungstic acid (PTA). For ultrathin sectioning, the material was treated according to Ryter and Kellenberger [12], dehydrated, and embedded in Epon. The sections were placed on formvar-coated copper grids and contrasted according to Reynolds [13].
Absorption spectra of the pigments were determined within the 350-1100-nm range using an SF-56A spectrophotometer (LOMO, St. Petersburg, Russia) for both whole cells (resuspended in 50% glycerol) and in acetone extracts.
Growth at different temperature, pH, and salinity values was studied using the basic medium with variations in the tested parameters and components. The required pH was adjusted by adding HCl or Na 2 CO 3 (18.25 or 53 g/L solutions, respectively). The concentration of organic substrates was 0.5 g/L. Biomass accumulation was determined at the stationary growth phase as OD 650 using a KFK-3 spectrophotometer (ZOMZ, Sergiyev Posad, Russia).
Response to oxygen was determined by growth in agar stabs (0.7%, w/v) in cottonplugged test tubes. The distance of noticeable bacterial growth from the surface of the medium was recorded after two weeks.

Genome Sequencing, 16S rRNA Gene Sequencing, Phylogenetic Analyses
The genome sequencing of Ect. strains B14B, A-7R and A-7Y was described in Imhoff et al., 2022 [1]. Average percentage nucleotide identity (ANIb) between the whole genomes was calculated using JSpecies [14]. Digital DNA-DNA hybridization (dDDH) data were obtained using the Type (Strain) Genome Server (TYGS) web server (https: //tygs.dsmz.de accessed on 5 May 2023) [15]. The program used the distance formula d4 to calculate a similarity based on sequence identity.
For synteny analysis, comparative genome regions were generated in BV-BRC using global PATRIC PGFam families to determine a set of genes that matched a focus gene [16]. All Ectothiorhodospira genomes were used in the search and compared to the B14B genome. The gene set was compared to the focus gene using BLAST and sorted by BLAST scores within BV-BRC [16]. The 'periplasmic sensor signal transduction histidine kinase PGF_03758479' was used as the focus gene to analyze synteny of the gene cluster.
The multiple sequence alignment for the 16S rRNA comparisons was performed using MUSCLE [20]. The evolutionary history was inferred by using the maximum likelihood method and Tamura-Nei model [21]. The bootstrap consensus tree inferred from 1000 replicates was taken to represent the evolutionary history of the taxa analyzed. Initial tree(s) for the heuristic search were obtained automatically by applying the maximum parsimony method. A discrete gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 0.6304)). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 57.63% sites). The branch lengths were measured in the number of substitutions per site. This analysis used 16S rRNA sequences that were obtained from the genome sequences or downloaded from Genbank when genomes were not available, with the following accession numbers: AM902494 ("Ect. imhoffii") and X93480 (Ect. sp. DSM 239). All positions with less than 95% site coverage were eliminated, i.e., fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position (partial deletion option). There were a total of 1464 positions in the final dataset. Evolutionary analyses were conducted in MEGA X [22], and iTOL was used for tree visualization [19]. Pairwise 16S rRNA comparisons were performed using LALIGN (https://www.ebi.ac.uk/Tools/psa/lalign/ accessed on 6 May 2023) using the genome-derived 16S rRNA sequences, and standard parameters were used.
The nucleotides of the complete DNA sequences of the genomes of the strains B14B, A-7R, and A-7Y were deposited in DDBJ/ENA/GenBank under accession no. JAJNQM000000000, JAJNQO000000000, and JAJNQN000000000, respectively [1].
The 16S rRNA gene sequences of the strains B14B, A-7R, and A-7Y were deposited with GenBank under accession numbers OQ618219, OQ618217, and OQ618218, respectively.

Morphology and Ultrastructure
The cell morphology of all new strains was similar. Cells were oval-to rod-shaped and sometimes slightly bent ( Figure 1A,B). The cells varied from 0.5 to 1.3 µm in width and 0.7 to 2.7 µm in length. Young cells of the strains B14B and A-7R were motile by means of a tuft of polar flagella ( Figure 1B). Motility was not found in strain A-7Y. Cells multiplied by binary fission. In the later stages of growth, the cells formed a mucous capsule ( Figure 1C). Ultrathin sections of cells of strain B14B showed a Gram-negative type of cell wall ( Figure 1D). The internal photosynthetic membranes were lamellae and located in ordered stacks, as in other species of the Ectothiorhodospiraceae ( Figure 1C,D). Dense oval inclusions, characteristic of polyphosphates, were seen regularly in the cells ( Figure 1C).
The 16S rRNA gene sequences of the strains B14B, A-7R, and A-7Y were deposited with GenBank under accession numbers OQ618219, OQ618217, and OQ618218, respectively.

Morphology and Ultrastructure
The cell morphology of all new strains was similar. Cells were oval-to rod-shaped and sometimes slightly bent ( Figure 1A,B). The cells varied from 0.5 to 1.3 µm in width and 0.7 to 2.7 µm in length. Young cells of the strains B14B and A-7R were motile by means of a tuft of polar flagella ( Figure 1B). Motility was not found in strain A-7Y. Cells multiplied by binary fission. In the later stages of growth, the cells formed a mucous capsule ( Figure 1C). Ultrathin sections of cells of strain B14B showed a Gram-negative type of cell wall ( Figure 1D). The internal photosynthetic membranes were lamellae and located in ordered stacks, as in other species of the Ectothiorhodospiraceae ( Figure 1C,D). Dense oval inclusions, characteristic of polyphosphates, were seen regularly in the cells ( Figure 1C).
The main pigment was bacteriochlorophyll a. Its presence on the absorption spectrum of living cells was indicated by absorption maxima at 378, 594, 799, 862 nm and a shoulder at (882-892) nm ( Figure 2). On the absorption spectrum of the acetone extract, bacteriochlorophyll a had a major maximum at 772 nm ( Figure 2). Bacteria contained carotenoids of the spirilloxanthin series, the presence of which was confirmed by the maximum at 491 and the shoulder (510-525) nm on the absorption spectrum of whole-cell pigments ( Figure 2). All known representatives of the genus Ectothiorhodospira are characterized by the presence of bacteriochlorophyll a and carotenoids of the spirilloxanthin series. The main pigment was bacteriochlorophyll a. Its presence on the absorptio trum of living cells was indicated by absorption maxima at 378, 594, 799, 862 n shoulder at (882-892) nm ( Figure 2). On the absorption spectrum of the acetone bacteriochlorophyll a had a major maximum at 772 nm ( Figure 2). Bacteria co carotenoids of the spirilloxanthin series, the presence of which was confirmed maximum at 491 and the shoulder (510-525) nm on the absorption spec whole-cell pigments ( Figure 2). All known representatives of the genus Ectothiorh are characterized by the presence of bacteriochlorophyll a and carotenoids of th loxanthin series.

Physiological Characteristics
Bacteria of the studied strains were anaerobic phototrophs. Photoaut growth with sulfide was weak (growth yields were roughly half that of phot trophically grown cultures). Photoautotrophically grown cells could be subc several times. Sulfide was oxidized to sulfate with the formation of S° as an o intermediate, which was deposited outside the cells. Thiosulfate and sulfite w oxidized under photoautotrophic conditions. Growth on thiosulfate and sulfite w sible only in the presence of organic substrates. Bacteria grew photoheterotrop with sulfide, thiosulfate and sulfite as a sulfur source for biosynthesis. Assimilat fate reduction was not possible. Under microaerobic conditions in the dark, grow only possible on organic substrates with sulfide as a source of reduced sulfur (Ta Under phototrophic conditions on a medium with sulfide as a sulfur sourc B14B used acetate, butyrate, yeast extract, lactate, malate, pyruvate, succina fumarate. In the presence of arginine, aspartate, glycerol, glucose, glutamate, propionate, peptone, propanol, sucrose, sorbitol, formate, fructose, citrate, growth was not observed. Nitrogen sources were ammonium, nitrate, and urea ( Growth factors were not required, but vitamin complex and yeast extract sti growth.
The best growth occurred at pH 9.0; the growth range was between pH 7.5

Physiological Characteristics
Bacteria of the studied strains were anaerobic phototrophs. Photoautotrophic growth with sulfide was weak (growth yields were roughly half that of photoheterotrophically grown cultures). Photoautotrophically grown cells could be subcultured several times. Sulfide was oxidized to sulfate with the formation of S • as an oxidation intermediate, which was deposited outside the cells. Thiosulfate and sulfite were not oxidized under photoautotrophic conditions. Growth on thiosulfate and sulfite was possible only in the presence of organic substrates. Bacteria grew photoheterotrophically, with sulfide, thiosulfate and sulfite as a sulfur source for biosynthesis. Assimilatory sulfate reduction was not possible. Under microaerobic conditions in the dark, growth was only possible on organic substrates with sulfide as a source of reduced sulfur (Table 1).
Under phototrophic conditions on a medium with sulfide as a sulfur source, strain B14B used acetate, butyrate, yeast extract, lactate, malate, pyruvate, succinate, and fumarate. In the presence of arginine, aspartate, glycerol, glucose, glutamate, lactose, propionate, peptone, propanol, sucrose, sorbitol, formate, fructose, citrate, ethanol growth was not observed. Nitrogen sources were ammonium, nitrate, and urea (Table 1). Growth factors were not required, but vitamin complex and yeast extract stimulated growth.
The best growth occurred at pH 9.0; the growth range was between pH 7.5 and 10.2 ( Figure 3). Growth was possible at NaCl concentrations ranging from 0% (maximum growth) to 8%; the best growth occurred at 0% NaCl ( Figure 3, Table 1). The studied strains were mesophiles. Optimum growth was at 25-35 • C; the temperature limits of growth were not determined.

Taxonomic Position
According to the results of whole-genome sequencing, the G + C content in the DNA of the studied strains B14B, A-7R, and A-7Y was 62.9-63.0 mol% [1]. Phylogenetic analysis based on the comparison of 16S rRNA gene sequences showed that these strains belonged to the genus Ectothiorhodospira and revealed the greatest affinity with species Ect. shaposhnikovii and Ect. vacuolata but places the new strains on a separate clade on the 16S rRNA tree ( Figure  4) [1]. , the dDDH and ANI results revealed low similarity between the studied strains and the related species that did not exceed accepted thresholds for bacterial species [14,24]. It is generally well known that whole genome-based comparisons of dDDH and ANI provide a higher resolution for species differentiation than single gene-based comparisons. dDDH values between strain B14B and the known Ectothiorhodospiraceae species varied from 19.7% with Ets. mongolicus DSM 15479 T to 38.8% with Ectothiorhodospira vacuolata DSM 2111 T ( Table 2).

Taxonomic Position
According to the results of whole-genome sequencing, the G + C content in the DNA of the studied strains B14B, A-7R, and A-7Y was 62.9-63.0 mol% [1]. Phylogenetic analysis based on the comparison of 16S rRNA gene sequences showed that these strains belonged to the genus Ectothiorhodospira and revealed the greatest affinity with species Ect. shaposhnikovii and Ect. vacuolata but places the new strains on a separate clade on the 16S rRNA tree (Figure 4)  , the dDDH and ANI results revealed low similarity between the studied strains and the related species that did not exceed accepted thresholds for bacterial species [14,24]. It is generally well known that whole genome-based comparisons of dDDH and ANI provide a higher resolution for species differentiation than single gene-based comparisons. dDDH values between strain B14B and the known Ectothiorhodospiraceae species varied from 19.7% with Ets. mongolicus DSM 15479 T to 38.8% with Ectothiorhodospira vacuolata DSM 2111 T ( Table 2). On the total genome phylogenetic tree, strains B14B, A-7R, and A-7Y formed a separate branch ( Figure 5) and had ANI values <90% (75.0-89.4%) relative to any of the other strains of the genus Ectothiorhodospira, but high ANI values between themselves (99.0-99.9%) [1]. ANI percentages > 95% have been suggested to recognize bacteria as belonging to the same species, while bacteria with ANI <90% would be recognized in most cases as sep-arate species [14]. The ANI with the most closely related species Ect. shaposhnikovii and Ect. vacuolata was 87.8-89.4%, which are all well below the arbitrary species cutoff of 95%. rRNA tree ( Figure  4) [1]. Although an LALIGN analysis (https://www.ebi.ac.uk/Tools/psa/lalign/ accessed on 6 May 2023) showed the highest 16S rRNA identity with Ect. shaposhnikovii DSM 243 T (99.4%; 1530/1540 bp); Ect. vacuolata DSM 2111 T (99.0%; 1511/1526 bp), Ect. haloalkaliphila ATCC 51935 T (99.0%; 1445/1459 bp), and Ect. magna DSM 22250 T (98.8%; 1454/1471 bp), the dDDH and ANI results revealed low similarity between the studied strains and the related species that did not exceed accepted thresholds for bacterial species [14,24]. It is generally well known that whole genome-based comparisons of dDDH and ANI provide a higher resolution for species differentiation than single gene-based comparisons. dDDH values between strain B14B and the known Ectothiorhodospiraceae species varied from 19.7% with Ets. mongolicus DSM 15479 T to 38.8% with Ectothiorhodospira vacuolata DSM 2111 T (Table 2).  the tree as outgroups. The phylogenetic tree was calculated in MegaX, and iTOL was used to draw the phylogenetic tree. Bootstrap values were generated from 1000 bootstrapping rounds. The new strains are indicated in red. On the total genome phylogenetic tree, strains B14B, A-7R, and A-7Y formed a separate branch ( Figure 5) and had ANI values <90% (75.0-89.4%) relative to any of the other strains of the genus Ectothiorhodospira, but high ANI values between themselves (99.0-99.9%) [1]. ANI percentages > 95% have been suggested to recognize bacteria as belonging to the same species, while bacteria with ANI <90% would be recognized in most cases as separate species [14]. The ANI with the most closely related species Ect. shaposhnikovii and Ect. vacuolata was 87.8-89.4%, which are all well below the arbitrary species cutoff of 95%.  A comparative systems analysis in BV-BRC between genomes of all known species of Ectothiorhodospiraceae and Halorhodospiraceae revealed that strains B14B, A-7R, and A-7Y had 15 unique PGFams that are not found in any other Ectothiorhodospiraceae genome, and they all contained a gene cluster for nitric oxide reductase (NOR), which is absent in all other species. The norCBQD cluster contains two NO reductase activation proteins as well as small and large NOR subunits (C and B) that show similarity to cytochrome oxidases. Nitric oxide reductase catalyzes bacterial denitrification by reducing nitric oxide (NO) to nitrous oxide (N 2 O), an obligatory step involved in the sequential reduction of nitrate to dinitrogen. The presence of norCB subunits is often used as a functional marker gene to identify denitrifying bacteria [25]. In addition to the PGFam analysis, a BLAST gene search using each of the 4 nor genes confirmed that they were only present in strains B14B, A-7R, and A-7Y and absent in all other Ectothiorhodospiraceae strains. There are no known studies of nitric oxide reduction in Ectothiorhodospiraceae, but strains B14B, A-7R, and A-7Y at least have the genetic capacity for nitric oxide reduction. An NCBI protein BLAST search of the NorQ protein showed Wenzhouxiangella sediminis (81.44% identity), Halomonas korlensis (76.03% identity), and Halomonas nitroreducans (74.53% identity) as the closest named relatives. All three of these are moderately halophilic, denitrifying Gammaproteobacteria isolated from saline, alkalophilic sediments, or soils [26][27][28].
The norCBQD gene cluster is found upstream from a 'periplasmic sensor signal transduction histidine kinase' (PGF_03758479) and a 'two-component transcriptional response regulator' (PGF_09666309). A synteny analysis in BV-BRC showed that these two genes are only conserved in six other Ectothiorhodospiraceae strains ( Figure 6), which, interestingly, are not the closest related to the B14B clade on the whole genome and 16S rRNA trees (Figures 4 and 5). This indicates that the norCBQD gene cluster in the three strains, B14B, A-7R, and A-7Y may have been obtained in a more recent evolutionary event. In addition to the general genomic difference, these genetic differences are another reason to consider these strains as a separate species [1]. A comparative systems analysis in BV-BRC between genomes of all known species of Ectothiorhodospiraceae and Halorhodospiraceae revealed that strains B14B, A-7R, and A-7Y had 15 unique PGFams that are not found in any other Ectothiorhodospiraceae genome, and they all contained a gene cluster for nitric oxide reductase (NOR), which is absent in all other species. The norCBQD cluster contains two NO reductase activation proteins as well as small and large NOR subunits (C and B) that show similarity to cytochrome oxidases. Nitric oxide reductase catalyzes bacterial denitrification by reducing nitric oxide (NO) to nitrous oxide (N2O), an obligatory step involved in the sequential reduction of nitrate to dinitrogen. The presence of norCB subunits is often used as a functional marker gene to identify denitrifying bacteria [25]. In addition to the PGFam analysis, a BLAST gene search using each of the 4 nor genes confirmed that they were only present in strains B14B, A-7R, and A-7Y and absent in all other Ectothiorhodospiraceae strains. There are no known studies of nitric oxide reduction in Ectothiorhodospiraceae, but strains B14B, A-7R, and A-7Y at least have the genetic capacity for nitric oxide reduction. An NCBI protein BLAST search of the NorQ protein showed Wenzhouxiangella sediminis (81.44% identity), Halomonas korlensis (76.03% identity), and Halomonas nitroreducans (74.53% identity) as the closest named relatives. All three of these are moderately halophilic, denitrifying Gammaproteobacteria isolated from saline, alkalophilic sediments, or soils [26][27][28].
The norCBQD gene cluster is found upstream from a periplasmic sensor signal transduction histidine kinase (PGF_03758479) and a two-component transcriptional response regulator (PGF_09666309). A synteny analysis in BV-BRC showed that these two genes are only conserved in six other Ectothiorhodospiraceae strains ( Figure 6), which, interestingly, are not the closest related to the B14B clade on the whole genome and 16S rRNA trees (Figures 4 and 5). This indicates that the norCBQD gene cluster in the three strains, B14B, A-7R, and A-7Y may have been obtained in a more recent evolutionary event. In addition to the general genomic difference, these genetic differences are another reason to consider these strains as a separate species [1]. Based on the combination of phenotypic and phylogenetic properties of strains B14B, A-7R, and A-7Y, it was proposed to classify them as a new species for which the name Ectothiorhodospira lacustris sp. nov. is proposed with strain B14B T (=DSM 116064 T =KCTC 25542 T = UQM 41491 T ) as type strain. Based on the combination of phenotypic and phylogenetic properties of strains B14B, A-7R, and A-7Y, it was proposed to classify them as a new species for which the name Ectothiorhodospira lacustris sp. nov. is proposed with strain B14B T (=DSM 116064 T =KCTC 25542 T = UQM 41491 T ) as type strain.
Cells are oval-to rod-shaped, sometimes slightly bent rods, 0.5-1.3 × 0.7-2.7 µm in size. Young cells are motile by means of polar flagella. Some strains are not mobile. Reproduction occurs by binary fission. Internal photosynthetic membranes are lamellar stacks. The color of cell suspensions is deep red. Photosynthetic pigments are bacteriochlorophyll a and carotenoids of the spirilloxanthin series. Absorption spectra of living cells have maxima at 378, 462, 491, 594, 799, and 862 nm, with shoulders at (510-525), (710-730), and (882-892) nm.
Growth occurs anaerobic in the light. Photoautotrophic growth with sulfide is weak (growth yield is roughly half that of photoheterotrophically grown cultures). Photoautotrophically grown cells can be subcultured several times. Sulfide is oxidized to sulfate, with S • , which is deposited outside the cells, as an intermediary product. Thiosulfate and sulfite are not oxidized under photoautotrophic conditions. Photoheterotrophic growth is possible with sulfide, thiosulfate, and sulfite as a sulfur source for biosynthesis under anoxic conditions. Assimilatory sulfate reduction is not possible. Chemolithotrophic growth is not possible under microoxic conditions in the dark in the presence of sulfide. Chemoheterotrophic growth is possible under microoxic conditions in the dark in the presence of sulfide.
Under phototrophic conditions and in the presence of sulfide, acetate, butyrate, yeast extract, lactate, malate, pyruvate, succinate, and fumarate are used. Arginine, aspartate, glycerol, glucose, glutamate, lactose, propionate, peptone, propanol, sucrose, sorbitol, formate, fructose, citrate, and ethanol are not used. Ammonia, nitrate, and urea are used as nitrogen sources and contain genes encoding a nitrogenase, so should be able to grow on N 2 . Growth factors are not required, but vitamin complex and yeast extract stimulate growth.
The bacterium is a mesophile with optimum growth at 25-35 • C. The temperature limits of growth have not been determined.
Habitat: low-mineralized soda lakes. The type strain B14B was isolated from biofilms of the shallow-water soda lake Bezymyannoye in central Mongolia (pH 9.0, mineralization 1-2 g/L, of which 80% sodium bicarbonate and 20% NaCl from the total mineralization of water).
The DNA G + C content of the species is 62.9-63.0 mol% according to genomic DNA data. The DNA G + C content of the type strain B14B is 62.9 mol% according to genomic DNA data.
The type strain B14B is deposited in the German Collection of Microorganisms and Cell Cultures GmbH (DSMZ) under accession number DSM 116064, in the Korean Collection for Type Cultures under accession number KCTC 25542, and in the Collection of Unique and Extremophilic Microorganisms under accession number UQM 41491.
The full genome DNA sequence of the type strain B14B T is registered under accession no. JAJNQM000000000.
The 16S rRNA gene sequence of the type strain B14B T is deposited at GenBank under accession no. OQ618219. Data Availability Statement: All sequence data are available in a publicly accessible repository of NCBI Genbank.