The algal strains of C. vulgaris and S. obliquus were identified by our results of sequencing of a partial region of 18S rDNA. The amplified 135 bp of 18S rDNA partial regions from C. vulgaris and S. obliquus DNA had 100% identity with Scenedesmus and Chlorella accessions in the NCBI’s GenBank (Fig. S1).
We also amplified a partial 700 bp (233 aa) region corresponding to the C-terminal domain of the iron hydrogenase large subunit. The identification of the hydA gene of C. vulgaris and S. obliquus provided us information about the genetic differences between phylogenetically distant groups of the algal hydrogenases, about non-synonymous changes in algal hydrogenases that could affect the protein’s functionality, and about the possibility of differential values in the production of H2 associated with the regulation in the transcription of the hydA gene in C. vulgaris and S. obliquus.
The hydrogenases of Scenedesmus obliquus and Chlorella vulgaris and their intra and interspecific evolutionary relationships with other algae
In order to establish genetic differences in the hydrogenase enzyme among algae, and to know if these non-synonymous changes affecting the 3D structure of the protein were related to the evolution of each species in the genus, we made a bioinformatic analysis including dendrogram construction and 3D structure modeling of the enzyme in 24 accessions, including 11 algal genera.
The results from our phylogenetic analysis indicated the formation of a main group including the genera Chlorella, Tetraspora, Scenedesmus, Tetradesmus, Monoraphidium, Raphidocells, Chlamydomonas, Volvox, and Coccomyxa. A second subgroup included the genus Tetraselmis. Some accessions of Chlorella, Raphidocelis, and Nannochloropsis formed a third group (Fig. 1), indicating that perhaps there are variants of hydrogenases in Chlorella and Rhapidocelis.
The dendrograms resulting from grouping algal accessions according to the characteristics of their hydrogenases did not in all cases correspond to the organisms’ phylogenetic classification. On one side, the genera Chlorella and Coccomyxa belong to the class Trebouxiophyceae but their accessions formed two subgroups within the main branch, similarly to the grouping pattern followed by accessions of Volvox and Tetraspora, both in the order Chlamydomonadales. On the other side, accessions of the genera Scenedesmus and Tetradesmus –belonging to the family Scenedesmaceae– grouped together reflecting a close relationship, and this same grouping pattern occurred between accessions in the genera Monoraphidium and Rhapidocelis both in the family Selenastraceae.
Our results from the 3D modeling of the HYD enzymes from C. vulgaris and S. obliquus showed found a 60% and 55% sequence identity with the [FeFe]-hydrogenase described for Chlamydomonas reinhardtii, respectively (Swanson et al. 2015) (Table S2). Almost all 3D models demonstrated similar structures, except that for Coccomyxa subellipsoidea XP_005643907.1 (Fig. 2, Fig. S2, Fig. S3). In almost all the 3D [FeFe]-hydrogenase structures modeled for the algal genera included in this study, we identified the aa residues responsible for the binding of the Iron/sulfur cluster (CPCACGCG; Fig. 2, Fig. S2, Fig. S3), but this sequence of residues was absent in the accessions Coccomyxa subellipsoidea XP_005643907.1 and Tetradesmus obliquus AA65921.1. Interestingly, our results also showed that accession Tetradesmus obliquus CAC34419.1 presented the binding aa sequence of the Iron/sulfur cluster, supporting the possible presence in algae of hydrogenase isoforms.
Isoforms of ferredoxins Fdx2–Fdx6 have been identified in Chlamydomonas reinhardtii that express differently depending on environmental conditions of oxygen, copper, iron, and ammonium concentrations (Meuser et al. 2011; Winkler et al. 2010) Furthermore, different isoforms of the genes coding for [FeFe] hydrogenases have been reported for C. variabilis (HydA1 and HydA2) and C. reinhardtii (HydE, HydF, and HydG). Also, studies of Bayesian inferred phylogeny reveal that the algal HydA are monophyletic and are nested between HydA genes of bacteria, fungi, and heterotrophic flagellates, which suggests that HydA emerged once early in chlorophyte evolution (Meuser et al. 2011) and also implies the possibility of finding genes coding for other isoforms of Hyd.
With the aim of knowing the variants of hyd genes, we built a network of haplotypes. For the construction of the minimal network of haplotypes of [FeFe]-hydrogenase we evaluated two structural domains in different microalgae species: the small subunit, and the terminal carbonyl domain (-C). The results of our evolutionary relationship analysis of both domains of the [FeFe]-hydrogenase in accessions of the genera Chlamydomonas, Nannochloropsis, Chlorella, Scenedesmus, Tetraselmis, Raphidocelis, Monoraphidium, Tetraspora, and Volvox indicated the presence of 17 genetic variants supported by 235 (small subunit domain) and 232 (C-terminal) total mutations. The small subunit variants of the [FeFe]-hydrogenase showed 192 discrete mutations among the analyzed genera and 43 divergent mutations in Chlorella ADK77883.1 and Tetraspora AMY63159.1, designated as V7CHLO and V14CHLO. The behavior of the C-terminal domain was similar, with 190 discrete and 42 divergent mutations. Accessions Chlorella fusca CAC83290.1, Tetradesmus obliquus CAC34419.1, Monoraphidium negletum XP_013906846.1, and Volvox carteri f. nagariensis XP_0029484897.1 presented the fixed variant (Fig. 3a).
Furthermore, our mutational analysis of [FeFe]-hydrogenase indicated conservation of 120 aa residues among accessions from species of Chlorella (CAC83290.1, ADK77883.1, PRW60372.1, AEA34989.1) and Scenedesmus (AXU2407.1) –including the V3CHLO and V4SCEN sequences we identified in this study– and distinguished an interspecific relationship between the V15CHL (small subunit domain and C-terminal) and V9CHLO (small subunit domain) variants of Chlorella (Fig. 3a and b). The variants showing a close mutational relationship were V8TETRA (Chlorella sorokiniana PRW60372.1) and V13TET (Tetradesmus obliquus AAG59621.1).
On the contrary, we observed a marked divergence among the V1CHLA (Chlamydomonas reinhardtii XP_00693376.1), V2NACH (Nannochloropsis gaditana XP_005854541.1), V5TETR (Tetraselmis sp. AHH85809.1), V11RAP (Raphidocelis subcapitata GBT94161.1), V12MON (Monoraphidium neglectum XP_013906846.1), V14TET (Tetraspora sp. AMY63159.1), V16NAN (Nannochloropsis salina TFJ80951.1), and V17VOL (Volvox carteri XP_002948487.1) variants (Peters et al. 2015; Sawyer and Winkler 2017)
Overall, this analysis demonstrates interspecific diversity of the [FeFe]-hydrogenase C-terminal and small subunit domains among microalgae taxa (Fig. 3). Interestingly, the evolutionary path of [FeFe]-hydrogenase indicated a definite relationship between species of Chlorella, Scenedesmus, and Tetradesmus in the minimal haplotype network, and the identification of eight missing variants and the reflected evolutionary pattern were supported by the analysis of the enzyme’s aa similarity in different species of microalgae (Fig. 1). Furthermore, the interspecific diversity we observed between species of microalgae suggests a divergent evolutionary strategy in the [FeFe]-hydrogenase enzyme, deriving in complementary functions –isofunctionalization and subfunctionalization– related to artificial selection pressure in some of the analyzed species (Peters et al. 2015).
Molecular hydrogen production and its relationship with the expression of hyd genes in Chlorella vulgaris and Scenedesmus obliquus
With the objective of knowing how the expression of the hyd genes correlates with the production of H2 in C. vulgaris and S. obliquus, we analyzed their H2 production and expression of hyd genes in conditions inducing the activity of the hydrogenase enzyme.
The determination of H2 in Chlorella vulgaris and Scenedesmus obliquus indicated some points to be highlighted (Table 1 and Fig. 4). When the microalgae were exposed to 1h of white light (140 µE m− 2 s− 1), a production of 2 ± 0.2 mL H2/L and 1 ± 0.30 of mL H2/L was determined, in Chlorella vulgaris and Scenedesmus obliquus, respectively. The H2 evolution kinetics of both microalgae were similar up to 5h of continuous exposure to white light, subsequently Scenedesmus obliquus showed better H2 values compared to Chlorella vulgaris, during the following 19h. During the evaluated time course (24h), it was observed that Scenedesmus obliquus (16 ± 0.50 mL H2/L) is a better H2 producer compared to Chlorella vulgaris (9.0 ± 0.40 mL H2/L) (Table 1 and Fig. 4a).
Table 1
Maximum productivity of H2 in cultures of C. vulgaris and S. obliquus under conditions of light (140 µE m− 2 s− 1)(data are shown as mean ± SD, n = 3).
Microalgae
|
Time (h)
|
Maximum production
H2 (mL L− 1)
|
Time
(min)
|
Maximum production
H2 (mL L− 1)
|
C. vulgaris
|
1
|
2.0a ± 0.20
|
30
|
1.0a ± 0.20
|
|
5
|
4.0b ± 0.36
|
60
|
2.0b ± 0.30
|
|
9
|
6.0c ± 0.24
|
90
|
2.0b ± 0.25
|
|
12
|
7.5c ± 0.30
|
120
|
2.0b ± 0.15
|
|
16
|
8.5d ± 0.28
|
150
|
2.0b ± 0.20
|
|
20
|
9.0d ± 0.30
|
180
|
2.0b ± 0.15
|
|
24
|
9.0d ± 0.40
|
210
|
3.0c ± 0.17
|
|
|
|
240
|
3.0c ± 0.23
|
|
|
|
270
|
3.0c ± 0.33
|
|
|
|
300
|
4.0c ± 0.32
|
S. obliquus
|
1
|
2.5A ± 0.40
|
30
|
2.0A ± 0.25
|
|
5
|
4.0B ± 0.23
|
60
|
2.0A ± 0.23
|
|
9
|
8.0C ± 0.33
|
90
|
2.0A ± 0.27
|
|
12
|
10D ± 0.23
|
120
|
2.0A ± 0.33
|
|
16
|
13E ± 0.32
|
150
|
2.0A ± 0.24
|
|
20
|
15E ± 0.15
|
180
|
2.0A ± 0.40
|
|
24
|
16E ± 0.50
|
210
|
2.0A ± 0.32
|
|
|
|
240
|
2.0A ± 0.40
|
|
|
|
270
|
4.0B ± 0.52
|
|
|
|
300
|
6.0C ± 0.32
|
On the other hand, monitoring the evolution of H2 during the first 5 h evaluated, with short intervals of 30 min (Fig. 4b), indicated that it is possible that H2 is detected in both microalgae from the first 30 min under white light exposure (Table 1 and Fig. 4b). In Chlorella vulgaris, significantly different H2 levels were detected after 180min (3h), while in Scenedesmus obliquus the H2 values were significantly different up to 270 min (4.5h), indicating differences between the kinetics of both microalgae. Ruiz Marin et al. (2020) had previously argued that both algae presented variations in the duration of the lag phase before initiating H2 production after being exposed to different light intensities, and suggested that the differences in the production of H2 is because they have different ability to adapt quickly to the new culture conditions. The adaptation of other microalgae to new culture conditions could depend in part on regulation of the enzymes involved in H2 production like HYD and FDX (Sun et al. 2013; Happe and Kaminski 2002). These variations in regulation will be dictated by the genome of each species analyzed.
Regulation of hyd genes in microalgae Scenedesmus obliquus and Chlorella vulgaris
The results of our analysis through qPCR of the relative expression of Hyd in the algae cultured in hydrogenase inducing media –made to know if the hyd gene is differentially regulated in C. vulgaris and S. obliquus– indicated that both strains showed higher levels of hyd expression at 1 h of culture (Fig. 5a). C. vulgaris showed 2.3 times more expression of the hyd gene than S. obliquus. Subsequently, in both strains, a decrease in the relative expression values was observed at 5 h of culture, until reaching a basal expression during the following 19 h.
It is known that the expression of the hyd genes is transient, and that their regulation can be inhibited by the products of photosynthesis (Weiner et al. 2018). Studies have shown that hydrogenase accumulates under dark anoxic adaptation conditions. Following such induction, exposure of algae to light supports high rates of H2 production, but H2 production ceases within a few minutes of illumination (Ghirardi 2015; Noone et al. 2017). This suggests that the observed increase in the expression of the hyd gene in C. vulgaris and S. obliquus within the first few hours after the transition from dark anoxia to light is the result of transcriptional regulation of the algal enzymes involved in this sudden change from dark to light, which takes place until the acclimatization of the algae and the activity of the hydrogenase becomes constant and basal (Fig. 5a).
With the aim of deepening the knowledge about regulation of hyd genes in the studied algae, we evaluated the relative expression of hyd genes in a 5 h period at 30 min intervals by qPCR. The results indicated higher expression levels at 30 min after the sudden change from darkness to light, reaching stable levels at 90 min (Fig. 5b).
This behavior pattern of hyd genes transcription in C. vulgaris and S. obliquus suggests that Hyd is involved in the production of H2 in the first 90 min, suggesting that possibly after the first 90 min hydrogenase is regulated by the levels of O2 as a product of the bio-photolysis of H2O that occurred during PSII in microalgae (Antal et al. 2003; Burlacot and Peltier 2018; Ghirardi et al. 2010). If we compare the levels of H2 production during the first 30 min in both microalgae, we observe that it is not differential between algae (2 mL H2/L) (Table 1 and Fig. 4b). Suggesting that the elevated expression levels in the first 30 minutes of light exposure (after a period of darkness) deregulates hydrogenase, enhancing its expression, later as the microalgae acclimatize to the new condition submitted, hydrogenase tends to regulate until to express sufficient levels to carry out its activity and contribute to the production of H2 in both microalgae (Fig. 5b)