Alia Ebrahim1
Manal A. Alfwuaires2
Mohammad H. Abukhalil3,4
Fawaz Alasmari5Fawad Ahmad1Rui Yao1Ying Luo1
Ying Huang1*†- 1Jiangsu Key Laboratory for Microbes and Genomics, School of Life Sciences, Nanjing Normal University, Nanjing, China
- 2Department of Biological Sciences, Faculty of Science, King Faisal University, Al-Ahsa, Saudi Arabia
- 3Department of Medical Analysis, Princess Aisha Bint Al-Hussein College of Nursing and Health Sciences, Al-Hussein Bin Talal University, Ma’an, Jordan
- 4Department of Biology, College of Science, Al-Hussein Bin Talal University, Ma’an, Jordan
- 5Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
Due to iron’s essential role in cellular metabolism, most organisms must maintain their homeostasis. In this regard, the fission yeast Schizosaccharomyces pombe (sp) uses two transcription factors to regulate intracellular iron levels: spFep1 under iron-rich conditions and spPhp4 under iron-deficient conditions, which are controlled by spGrx4. However, bioinformatics analysis to understand the role of the spGrx4/spFep1/spPhp4 axis in maintaining iron homeostasis in S. pombe is still lacking. Our study aimed to perform bioinformatics analysis on S. pombe proteins and their sequence homologs in Aspergillus flavus (af), Saccharomyces cerevisiae (sc), and Homo sapiens (hs) to understand the role of spGrx4, spFep1, and spPhp4 in maintaining iron homeostasis. The three genes’ expression patterns were also examined at various iron concentrations. A multiple sequence alignment analysis of spGrx4 and its sequence homologs revealed a conserved cysteine residue in each PF00085 domain. Blast results showed that hsGLRX3 is most similar to spGrx4. In addition, spFep1 is most closely related in sequence to scDal80, whereas scHap4 is most similar to spFep1. We also found two highly conserved motifs in spFep1 and its sequence homologs that are significant for iron transport systems because they contain residues involved in iron homeostasis. The scHap4 is most similar to spPhp4. Using STRING to analyze protein-protein interactions, we found that spGrx4 interacts strongly with spPhp4 and spFep1. Furthermore, spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that the three proteins play cooperative roles in iron homeostasis. At the highest level of Fe, spgrx4 had the highest expression, followed by spfep1, while spphp4 had the lowest expression; a contrast occurred at the lowest level of Fe, where spgrx4 expression remained constant. Our findings support the notion that organisms develop diverse strategies to maintain iron homeostasis.
1 Introduction
Due to iron’s critical role in cell metabolism, maintaining its balance is essential for most organisms (Maio et al., 2022). Low-intracellular iron concentration inhibits the activity of Fe-dependent enzymes, while high iron accumulation stimulates reactive oxygen species (ROS) formation and disrupts other metal-trafficking pathways (Vachon et al., 2012; Chen et al., 2022). Therefore, cells must tightly regulate iron uptake to maintain an equilibrium between adequate and excessive concentrations (Liu et al., 2022). In this regard, the fission yeast Schizosaccharomyces pombe (sp) uses two transcription factors, spFep1 and spPhp4, to regulate iron levels within cells (Mercier et al., 2006; Jbel et al., 2009). spFep1 is activated when iron is abundant (Jbel et al., 2009), whereas spPhp4 is activated when iron is scarce (Mercier et al., 2006; Attarian et al., 2018; Misslinger et al., 2021). spFep1 and spPhp4 collaborate to regulate each other’s transcriptomic expression levels in an iron-dependent manner (Vachon et al., 2012; Berndt et al., 2021).
When the amount of intracellular iron is high, spFep1 binds to chromatin, preventing the expression of several genes involved in iron acquisition (Pelletier et al., 2002; Labbé et al., 2007; Mao and Chen, 2019; Robinson et al., 2021). However, iron transporter gene transcription at low iron concentrations is activated due to spFep1’s inability to bind chromatin (Jbel et al., 2009; Berndt et al., 2021). Several fungi possess iron-dependent regulators with similar sequences and functions to spFep1, but not S. pombe (Chao et al., 2008; Haas et al., 2008; Hwang et al., 2012; Alkafeef et al., 2020). It was found that spFep1 interacts with monothiol glutaredoxin spGrx4 (Vachon et al., 2012; Miele et al., 2021). For instance, spGrx4 is essential for inhibiting spFep1 function when cells transition from iron-rich to iron-poor conditions (Jbel et al., 2011; Kim et al., 2011; Daniel et al., 2020). spGrx4’s glutaredoxin (GRX) domain forms a connection with spFep1 via its N-terminal (Gupta and Outten, 2020; Talib and Outten, 2021). The spFep1 DNA-binding domain is inactivated by the spFep1-GRX-binding domain, which in turn prevents the spFep1 DNA-binding domain from binding to chromatin, thereby inhibiting the expression of the gene of interest (Miele et al., 2021).
spphp4 encodes a subunit of the CCAAT-binding protein complex that plays a critical role in regulating iron homeostasis (Mühlenhoff et al., 2020). This complex also contains three additional subunits, spPhp2, spPhp3, and spPhp5, which are highly expressed (Jbel et al., 2009; Jbel et al., 2011; Ohtsuka et al., 2021; Protacio et al., 2022). However, spphp4 transcript levels are increased in iron-deficient conditions and decreased in iron-rich environments (Jbel et al., 2009). Although spPhp4 is not required for the complex’s DNA-binding activity, it is responsible for the Php complex’s transcriptional repression ability (Brault and Labbé, 2020). Furthermore, researchers used genome-wide microarray analysis and iron-starvation conditions and discovered that spPhp4 can coordinate the suppression of 86 genes (Jung et al., 2010; Martínez-Pastor and Puig, 2020). Many of these genes produce proteins involved in iron-dependent metabolic processes, such as heme biosynthesis, sulfur-iron cluster assembly, the tricarboxylic acid pathway, and mitochondrial respiration. Under iron-deficient conditions, DNA microarray analysis also revealed that the gene encoding the iron-responsive transcriptional inhibitor spFep1 exhibited spPhp4-dependentdown-regulation (Martínez-Pastor and Puig, 2020).
To our knowledge, bioinformatics analysis to understand the role of spGrx4, spFep1, and spPhp4 in maintaining iron homeostasis is still elusive. In this study, we conducted a bioinformatics analysis on S. pombe and their sequence homologs in Aspergillus flavus (af), Saccharomyces cerevisiae (sc), and Homo sapiens (hs) to understand the role of spGrx4, spFep1, and spPhp4 in maintaining iron homeostasis. Under various ferrous concentrations, expression analysis was also conducted to gain a deeper understanding of the critical roles of spGrx4, spPhp4, and spFep1 in iron homeostasis. In light of this, we proposed a model in which spFep1 and spPhp4 interact to regulate the expression of each other, thus controlling the strict regulation of intracellular iron concentrations (Jung et al., 2010). Our analyses of spGrx4 revealed a conserved cysteine residue in each PF00085 domain, and that human hsGLRX3 is most similar to spGrx4. spFep1 is highly correlated with scDal80. The scHap4 is most similar to spPhp4. Our analysis also showed that Grx4 interacts strongly with spPhp4 and spFep1, and that spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that our three proteins play cooperative roles in iron homeostasis. At the highest level of Fe, spgrx4 had the highest expression, followed by spfep1, while spphp4 had the lowest expression; a contrast occurred at the lowest level of Fe, where spgrx4 expression remained constant.
2 Materials and methods
2.1 Proteins data set and sequence analyses
The protein sequences of spGrx4, spPhp4, and spFep1 sequences were used as queries in BLAST searches against fungal and human genome databases, such as the NCBI, OrthoFinder and PHYLome, to identify sequence homologs. The E-value threshold of 0.05 was used to retrieve all protein sequences. We used TeUniProt Knowledgebase (UniProtKB) to determine whether the candidate proteins contained the conserved active site. Using SMART, Simple Modular Analysis Research Tool, the presence of conserved domains in Te candidate proteins was determined (Letunic and Bork, 2018). We utilized the Multiple EM for Motif Elicitation (MEME Suite) Version 5.4.1 (Peng et al., 2018) with the default settings to identify conserved sequence motifs (El-Sappah et al., 2021a; El-Sappah et al., 2021b).
2.2 Phylogenetic analysis
Multiple amino acid sequences were aligned to construct a phylogenetic tree using the MUSCLE algorithm incorporated into MEGA X software; version 10.0.5 (Kumar et al., 2018; Li et al., 2022). The phylogenetic tree was constructed using maximum likelihood (Jones-Taylor-Thornton model, 100 bootstrap replications) (Li et al., 2021).
2.3 Protein–protein interactions, prediction, and protein modeling
According to Kelley et al. (2015), the Phyre2 online webserver (sbg.bio.ic.ac.uk/phyre2/) was used for protein prediction, modeling, and analysis in the intensive mode. The STRING database (https://string-db.org) was used to analyze protein-protein interactions using the amino acid sequences.
2.4 Growth conditions and ferroptosis treatments
S. pombe cells were cultured overnight in YES medium at 30°C before diluting to an OD600 of 0.2 and re-culturing for the same period. The cells were then assembled and assigned an OD600value of 3.0. A tenfold dilution series was prepared, and each dilution was spotted on plates containing 3 μl YES media containing 0.25, 0.5, 3, or 6.0 mM Fe2 (SO4)3, or was left untreated.
2.5 RNA isolation and real-time quantitative PCR
Following the manufacturer’s instructions, total RNA was isolated from S. pombe cells using the High Pure RNA Isolation kit (Roche). Next, quantitative RT-PCR was conducted (Stratagene) using a FastStart SYBR Green Master kit (FSSGM; Roche) and the Mx3005P system. Four specific primers were employed for the following genes: act1, grx4, fep1, and php4, as follows (act1: Forward (F) 5′-ACTATGTATCCCGGTATTGCC-3′, Reverse (R) 5′-GACAGAGTATTTACGCTCAGG-3′; grx4: (F) 5′-GCGTCGCCTATTGTGCAAGG-3′, (R) 5′-GTTGGTTGTCGGCAGGTTCG-3′; fep1: (F) 5′-CCCGAAGAGCCACCCTCAAA-3′, (R) 5′-CCGGCGAGGTAGAGGATTGG-3′; php4: (F) 5′-CAAGACCTGGGCCGACTTCA-3′, (R) 5′-AGCCTGCAGCATCTCCAACT-3′). We constructed the first-strand cDNA template in accordance with the manufacturer’s guidelines. Next, 50 µl of ultra-pure nuclease-free H2O was added to a mixture of 25 ml of FSSGM, 0.2 µM of each F and R primer, and 0.1 µg of cDNA. PCR was configured as follows: Pre-incubation at 95°C for 10 min, then 40 cycles of 95°C for 10 s, 55°C for 10 s, and 72°C for 20 s. The final step consisted of gradually increasing the temperature from 55°C to 95°C at a rate of 1°C/10 s to collect melting curve data (Islam et al., 2022). The reference (actin gene, act1) and target genes were serially diluted (1:10, 1:100), along with a non-template control, and added to each assay. Using melting curve analysis, we validated the specificity of each amplification reaction. In addition, we adjusted the expression levels for the reference gene act1, and Pfaffle (2001) calculated the relative gene expression levels using Pfaffle’s (2001) method.
3 Results
3.1 spGrx4, spFep1, and spPhp4 alignment and domain analysis
3.1.1 spGrx4
As part of our study, we used the protein sequences of spGrx4, spFep1, and spPhp4 as query sequences, and performed BLAST searches against the human and fungal genome databases to determine the functional and structural diversity of these proteins. For further examination of spGrx4, eight spGrx4 homologs were chosen: two from S. cerevisiae (scGrx4 and scGrx3), three from A. flavus (afCA14001870, af F9C072156032, and afAFLA70126g002401), and three from humans. We found that hsGLRX5 was the most closely related by sequence homology to spGrx4 with a 50% identification rate (Table 1), while hsGLRX3 was predicted to be the functional homolog of spGrx4 based on the two major ortholog predictor tools (OrthoFinder and PHYLome). Using the maximum-likelihood method, we constructed a phylogenetic tree of Grx4 based on their similarity to estimate their evolutionary relationships (Figure 1). Our research demonstrates that spGrx4 is strongly associated with scGrx4. We also aligned the sequences of spGrx4 and its sequence homologs in humans, A. flavus, and S. cerevisiae to comprehend the functional and structural diversity of the proteins (Figure 2). Our multiple sequence alignment revealed that all selected Grx4 PF00085 domains contain a conserved cysteine (Cys) residue (Cys in spGrx4, afCA14001870, af F9C072156032, scGrx4, scGrx3, hsGLRX3 and hsGLRX5, and selenocysteine (Sec) in hsGPX1, hsGPX4, and hsGPX6). Moreover, our analysis revealed the presence of Cys172 in Grx4 PF00462 (All of selected proteins except in hsTXN).
TABLE 1. Sequence comparison of spGrx4 to its sequence homologs in A. flavus, H. sapiens, and S. cerevisiae.
FIGURE 1. Bayesian evolutionary tree of spGrx4, spFep1, and spPhp4 and their sequence homologs in A. flavus, S. cerevisiae, and humans. spGrx4, spFep1, and spPhp4 are colored red. In the analysis, the strict clock model was implemented. Using the MCMC algorithm, the trustworthiness of phylogeny nodes was assessed. At the nodes, Bayesian posterior probabilities are displayed. The scale bar indicates that there are 0.2 amino acid substitutions at each site.
FIGURE 2. Multiple sequence alignment of the PF00085 and PF00462 domain sequences in spGrxp4 and its sequence homologs in A. flavus, S. cerevisiae, and humans using the MEGAX10.1.8 software. The Grx4 domains were discovered through a search of the SMART database. The two characteristic residues (a Cys residue in the PF00085 domain and a Cys172 residue in the PF00462 domain) are enclosed in square brackets. They are highlighted according to the biochemical properties of the amino acid residues.
3.1.2 spFep1
Ten spFep1 homologs were chosen (Table 2): two from S. cerevisiae (scGat1 and scDal80), two from A. flavus (afSreA and afSreP), and six from humans (hsGATA1 to hsGATA6). According to the NCBI basic sequence blast results and the two major ortholog predictors (OrthoFinder and PHYLome) tools, scDal80 is most similar to spFep1, 60.9% identity (Table 2). To determine the evolutionary relationships between these Fep1 sequences, we constructed a phylogenetic tree of spFep1 and its sequence homologs using the maximum-likelihood method and a similarity score (Figure 1). Our study revealed that spFep1 is most similar to scDal80. Additionally, multiple sequence alignments of spFep1 and its sequence homologs in humans, A. flavus and S. cerevisiae were performed. This analysis identified four conserved Cys residues in all selected Fep1 PF00320 domains (Figure 3). According to the SMART database, our analysis revealed that spFep1 and its sequence homologs in human, A. flavus and S. cerevisiae all contain the GATA zinc finger (PF00320) domain, which is depicted in Figure 3. In addition, we discovered two highly conserved motifs in spFep1 (Figure 4) that are significant for iron transport systems because they contain residues involved in iron homeostasis (Figure 5).
TABLE 2. Sequence comparison of spFep1 to its sequence homologs in A. flavus, H. sapiens, and S. cerevisiae.
FIGURE 3. Multiple sequence alignment of sequences of the PF00320 domain in spFep1 and its sequence homologs in A. flavus, S. cerevisiae, and humans using the MEGA_X_10.1.8 software. The domain was identified by searching the SMART database. The four characteristic Cys residues are boxed while amino acid residues are highlighted according to the biochemical properties of the amino acid residues.
FIGURE 4. The conserved motif analysis of selected proteins in Figure 1 using the MEME Suite, where the motif composition of these proteins using 10 conserved motifs is represented by a unique color mentioned in the box on the top right.
FIGURE 5. The 10 conserved motifs of selected proteins in Figure 1. The sequence logos was created by the MEME Suite.
3.1.3 spPhp4
We chose 5 spPhp4 homologous: one from S. cerevisiae (scHap4), one from A. flavus (afHapX), and three from humans (hsHMGN3, hsNOLC1, and hsOPTN). The scHap4 is more similar to spPhp4 than others (Table 3). It shared the highest similarity with spFep1. The Php4 phylogenetic tree (Figure 1) demonstrates that spPhp4 is most related to scHap4. Furthermore, our multiple sequence alignment revealed the presence of two conserved Cys residues in spPhp4 and its homologs, despite the absence of these residues in its PF10297 domains (Figure 6). Based on the SMART database, our study demonstrated that A. flavus, S. cerevisiae, and spPhp4 all contain the Hap4_Hap_bind (PF10297) domain, as depicted in Figure 6. In addition, the MEME Suite identified a highly conserved motif in these Php4 (Figures 4, 5).
TABLE 3. Sequence comparison of spPhp4 to its sequence homologs in A. flavus, H. sapiens, and S. cerevisiae.
FIGURE 6. Multiple sequence alignment of sequences of the PF10297 domain in spPhp4 and its homologs in A. flavus, S. cerevisiae and humans. The domain was discovered through a search of the SMART database. They are highlighted according to the biochemical properties of the amino acid residues.
3.2 3D structure and protein-protein interaction
The spGrx4 protein contains 244 amino acids, with 28 negatively charged residues (Asp and Glu) and 19 positively charged residues (Arg and Lys). At its N-terminus, spGrx4 has an additional domain containing a WAAPCK motif, similar to the WCGPCK motif in thioredoxin’s active site (Figure 7). spFep1 contains 564 amino acids, 43 of which are positively charged (Arg and Lys) and 53 negatively charged (Asp and Glu). Moreover, two Cys2-Cys2 zinc fingers are present in the spFep1 protein (Figure 7). In addition, the Cys-X5-Cys-X8-Cys-X2-Cys motifs are positioned between the two zinc fingers (Figure 7). spPhp4 is composed of 295 amino acids, 38 (Asp and Glu) and 46 (Arg and Lys) of which are positively charged. Subsequently, using STRING to analyze protein-protein interactions, it was determined that spGrx4 interacts strongly with spPhp4 and spFep1. Furthermore, spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that our three proteins play cooperative roles in iron homeostasis (Figure 7).
FIGURE 7. Protein analysis; (A) The predicted three-dimensional structures of spGrx4, spFep1, and spPhp4. (B) The Protein-protein interaction network building among spGrx4, spFep1, and spPhp4 was performed using the Phyre2 server in an intensive mode. Rainbow colors visualized protein models were visualized by rainbow colors in the direction from N to C terminus. Unknown 3D structures are represented by empty nodes, known or predicted 3D structures by filled nodes, query proteins, and the first shell of interactors by colored nodes, and the second shell of interactors by white nodes.
3.3 spgrx4, spfep1, and spphp4 expression under different iron conditions
At the highest level of iron, spgrx4 expressed itself more than spfep1, while spphp4 expressed itself the least. By decreasing the iron concentration, spgrx4 continued to exhibit the highest expression level, while fep1 expression decreased and spphp4 expression increased (Figure 8). When iron levels are low, S. pombe cells inhibit the transcription of many genes that encode iron-utilizing proteins. spPhp4 represses iron utilization/storage genes by interacting with the CCAAT-binding Php complex composed of spPhp2, spPhp3, and spPhp5.
FIGURE 8. The qRT–PCR analysis of spgrx4, spfep1 and spphp4 at varying concentrations of ferrous. The β-actin gene served as the internal control. Error bars were used to depict the standard error of the three scientific repeats. The various letters represent statistically significant differences (p < 0.05) between the treatment and the control group.
4 Discussion
Iron is a key element required for several fundamental organismal and cellular processes, including nucleic acid replication, mitochondrial respiration, oxygen transport, intermediary and xenobiotic metabolism, cell signaling, and host defense (Wang and Babitt, 2018). The property of iron as a transition metal to easily donate or accept electrons to contribute to oxidation reduction reactions is the basis on which these biological functions depend (Wang and Babitt, 2018). Cells must tightly regulate iron uptake to balance adequate and excessive iron levels (Katsarou and Pantopoulos, 2020). S. pombe uses two transcription factors, spFep1 under iron-abundant conditions and spPhp4 under iron-abundant conditions iron-deficient conditions, which are controlled by spGrx4 (Alkafeef et al., 2020). spFep1 and spPhp4 reciprocally regulate the expression of each iron-dependent genes, thereby coordinating their activities at the transcriptional level. Our analysis of the domain, protein structure, and gene expression patterns of the spGrx4, spFep1, and spPhp4 under different iron concentrations will help in understanding their roles in maintaining iron homeostasis. Our analysis showed that the human hsGLRX3 is most similar to the spGrx4, and that spGrx4 is strongly similar to scGrx4. All selected TRX (PF00085) domains in spGrx4 and its homologs contain a conserved cysteine (Cys) residue. In addition, our analysis revealed the presence of Cys172 in the GRX domain (PF00462). The highly conserved residues C172GFS, which are required for Grx4’s cellular functions, are located in the C-terminal GRX-like domain of the protein (Herrero and De La Torre-Ruiz, 2007). Furthermore, the strength of the interaction between spGrx4 and spPhp4 bound to iron depends on certain conserved Cys residues, as spGrx4 depends on Cys172 and spPhp4 depends on both Cys221 and Cys227 (Dlouhy et al., 2017). spGrx4 and its sequence homologs in human, A. flavus and S. cerevisiae all contain the thioredoxin (PF00085) and GRX (PF00462) domains (Figure 2), which contribute electrons to ribonucleotide reductase and are involved in thiol-redox and ribonucleotide reductase (Holmgren, 1989). In addition, yeast two-hybrid mapping reveals that the C-terminal region of spPhp4 physically interacts with spGrx4 via the C-terminal GRX-like and N-terminal TRX-like domains of spGrx4 (Vachon et al., 2012). Our analysis revealed that, scDal80 is the most similar to spFep1. However, scDal80 and spFep1 are not orthologs because they start serving different functions. scDal80 participates in nitrogen metabolism, whereas spFep1 participates in iron metabolism. scDal80 inhibits a number of independent inducer genes associated with nitrogen metabolism (Cunningham and Cooper, 1993). scDal80 negatively regulates a wide range of nitrogen-catabolic genes, and it encodes a protein with a zinc finger motif that bears extensive homology to the metazoan GATA-binding family of transcriptional activators (Cunningham and Cooper, 1993). It is worth to noting that scGzf3, which is a paralog of scDal80, inhibits nitrogen catabolic gene expression by competing with scGat1 for GATA site binding (Palavecino-Ruiz et al., 2017). These results are consistent with the notion that S. cerevisiae and S. pombe use different mechanisms to maintain iron homeostasis.
Our analysis identified four conserved Cys residues in all selected GATA zinc-finger domains (PF00320) in spFep1 and its homologs. A region containing four conserved Cys residues that regulate iron homeostasis in S. pombe is flanked by the DNA-binding motifs of these GATA-type transcription factors, which typically contain one or two zinc finger motifs (Attarian et al., 2018). We also found that spFep1 and its sequence homologs in human, A. flavus, and S. cerevisiae all contain the GATA zinc finger (PF00320) domain (Figure 3). This domain regulates iron homeostasis (Blaiseau et al., 2001). When cells are iron-replete, spFep1 binds to GATA sequences in these genes and inhibits their transcription, thereby downregulating iron transport systems (Dlouhy et al., 2017). Moreover, we observed two highly conserved motifs in spFep1 that are significant for iron transport systems because they contain residues involved in iron homeostasis.
We also showed that the scHap4 is more similar to spPhp4 and shared the highest degree of similarity with spFep1, and that the Cys-rich protein plays a significant iron-binding role. Comparable to a similar Cys-rich region with an iron-sensing function in HapX, the Php4 protein contains a Cys residue close to its carboxyl terminus (Hortschansky et al., 2007). By using genome-wide microarray analysis, Mercier et al. (2008) determined that spPhp4 is involved in the inhibition of 86 genes under iron-starvation conditions. Our study also demonstrated that spPhp4 and its sequence homologs in A. flavus, S. cerevisiae, and humans all contain the Hap4_Hap_bind (PF10297) domain. This domain is essential for scHap4 to form the scHap complex by allowing it to bind to scHap2, scHap3, and scHap5 (McNabb and Pinto, 2005). spGrx4 interacts with spPhp4 through its C terminus, which includes amino acid (aa) residues 152–254 (Mercier and Labbé, 2009). spPhp4 and scHap4 share a small amount of sequence homology overall. scHap4 contains a conserved domain, which is required for to link scHap4 to scHap2/scHap3/scHap5 (McNabb and Pinto, 2005). Furthermore, the MEME Suite identified a highly conserved motif in these spPhp4 (Figure 4). This motif is important for iron transport systems because it contains iron homeostasis-involved residues (Figure 5). Unlike spPhp4, scHap4 plays a role in controlling the fermentation/respiration switch (Lascaris et al., 2003). These results also support the notion that S. cerevisiae and S. pombe use different strategies to control iron homeostasis.
Our study also revealed that spGrx4 interacts strongly with spPhp4 and spFep1, and that spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that our three proteins play cooperative roles in iron homeostasis (Gupta and Outten, 2020; Misslinger et al., 2021). spGrx3, spGrx4, and spGrx5 are three monothiol glutaredoxins found in S. pombe (Chung et al., 2005; Berndt et al., 2021). They are primarily localized at the nuclear rim, throughout the cell (cytosol and nucleus), and in mitochondria, respectively (Chung et al., 2005; Mercier and Labbé, 2009). At the active site of spGrx3, spGrx4, and spGrx5, which is part of the GRX-like domain, one Cys residue is highly conserved (Chung et al., 2005; Berndt et al., 2021). Earlier research discovered that the spGrx3 and spGrx4 monothiol glutaredoxins contain the N-terminal TRX-like domain (Herrero and de la Torre-Ruiz, 2007). This domain is required for the nuclear localization of monothiol glutaredoxins containing the TRX domain (Molina et al., 2004). The spGrx4 protein is associated with oxidative, osmotic, nitrosative, and iron-dependent stress responses (Chung et al., 2005; Kim et al., 2005; Mercier and Labbé, 2009).
At the highest level of Fe, spgrx4 expressed itself more than spfep1, while spphp4 expressed itself the least. By decreasing the Fe concentration, grx1 continued to exhibit the highest expression level, while spfep1 expression decreased and spPhp4 expression increased (Figure 8). When iron levels are low, S. pombe cells inhibit the transcription of many genes that encode iron-utilizing proteins (Labbé et al., 2007; Liu et al., 2022). Interacting with the CCAAT-binding complex composed of spPhp2, spPhp3, and spPhp5, spPhp4 modulates transcription (Mercier et al., 2006). spPhp4 is inactivated under iron-rich conditions, allowing iron-dependent proteins to be translated (Mercier et al., 2006). S. pombe uses spFep1 and spPhp4 in iron-rich and iron-deficient conditions, respectively (Brault et al., 2015) to regulate the amount of iron in cells. Furthermore, spfep1 and spphp4 mutually regulate one another’s expression using iron-dependent mechanisms (Labbé et al., 2013). Moreover, spFep1, a GATA-type transcriptional inhibitor, primarily regulates iron acquisition related genes (Pelletier et al., 2002). Under iron-replete conditions, spFep1 binds to the GATA sequences in iron transport genes and inhibits their transcription, thereby downregulating the iron transport systems (Jbel et al., 2009; Jung et al., 2021; Miele et al., 2021). In S. pombe, spPhp2, spPhp3, and spPhp5 operate as heterotrimers to bind CCAAT cis-acting elements and activate gene expression (Mercier et al., 2006). This complex regulates the expression of genes encoding proteins involved in iron-dependent metabolic processes such as amino acid biosynthesis, iron-sulfur cluster assembly, the TCA cycle, mitochondrial respiration, and heme biosynthesis (Mercier et al., 2008). spPhp4 binds to the CCAAT-binding complex and transforms it from an activator to a suppressor in the absence of iron (Mercier et al., 2006). Therefore, Php4’s primary function is iron conservation, as depleted cellular iron inhibits iron-dependent protein-coding gene expression (Brault et al., 2016).
The activities of spFep1 and spPhp4 are also controlled post-translationally by the cytosolic CGFS glutaredoxin spGrx4 (Jbel et al., 2009; Mercier and Labbé, 2009; Encinar del Dedo et al., 2015). spGrx4 inhibits binding of spFep1 to the promoters of iron transport genes under iron-limiting conditions (Jbel et al., 2011; Encinar del Dedo et al., 2015). spGrx4 regulates interactions between spPhp4 and the spPhp2/3/5 complex (Mercier and Labbé, 2009). Under iron-depleted conditions, spPhp4 binds to the CCAAT-binding HAP complex (Php2/3/5) to suppress iron utilization/storage genes (Mercier and Labbé, 2009). However, when iron is present in high concentrations, spPhp4 dissociates from spPhp2/3/5 and is then re-localized to cytoplasm by spCrm1 (Mercier and Labbé, 2009). A yeast two-hybrid study revealed that spGrx4 physically interacts with spPhp4’s C-terminal region via its N- and C-termini TRX and GRX domains (Vachon et al., 2012). In contrast, iron-dependent binding to the TRX domain requires Cys172 at the conserved CGFS active site, whereas binding to the GRX domain is weak and constitutive. Cys172 is essential for spGrx4-dependent inhibition of spPhp4 in strains expressing C172 S/A spGrx4 mutants (Kim et al., 2011; Encinar del Dedo et al., 2015; Dlouhy et al., 2017). Cys221 and Cys227 in spPhp4 may be involved in the iron-mediated inactivation of spPhp4 by interacting with the GRX domain of spGrx4 (Dlouhy et al., 2017).
5 Conclusion
Here, we performed a bioinformatics analysis on S. pombe and their sequence homologs in A. flavus, S. cerevisiae, and H. sapiens, and an expression analysis with different iron concentrations to gain a deeper understanding of the critical roles of spGrx4, spPhp4, and spFep1 in iron homeostasis. Our analyses of Grx4 revealed a conserved cysteine residue in each TRX domain, and that human hsGLRX3 is most similar to spGrx4. Moreover, spFep1 is highly correlated with scDal80, whereas scHap4 is most similar to spFep1 and to spPhp4. Our analysis also showed that spGrx4 interacts strongly with spPhp4 and spFep1, and that spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that our three proteins play cooperative roles in iron homeostasis. At the highest level of Fe, spgrx4 had the highest expression, followed by spfep1, while spphp4 had the lowest expression; a contrast occurred at the lowest level of Fe, where spgrx4 expression remained constant. Our analysis improves our understanding of the structure and function of spGrx4, spFep1, and spPhp4.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
AE and YH designed and performed the experiments. AE, MAA, MHA, FzA, FdA, RY, and YL conducted bioinformatics analyses. AE, MAA, MHA, FzA, RY, and YL analyzed the data. AE and YH prepared the manuscript. All authors read, reviewed and approved the final version of the manuscript.
Funding
This work was supported by the Researchers supporting project number (RSP 2022R235), King Saud University, Riyadh, Saudi Arabia.
Acknowledgments
We are grateful to Valerie Wood from PomBase for the fruitful discussion.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alkafeef, S. S., Lane, S., Yu, C., Zhou, T., Solis, N. V., and Filler, S. G. (2020). Proteomic profiling of the monothiol glutaredoxin Grx3 reveals its global role in the regulation of iron dependent processes. PLoS Genet 16, e1008881. doi:10.1371/journal.pgen.1008881
Attarian, R., Hu, G., Sánchez-León, E., Caza, M., Crol, I. D., Do, E., et al. (2018). The monothiol glutaredoxin Grx4 regulates iron homeostasis and virulence in Cryptococcus neoformans. MBio 9, e02377-18–e02318., doi:10.1128/mBio.02377-18
Berndt, C., Christ, L., and Rouhier, N. (2021). Glutaredoxins with iron-sulphur clusters in eukaryotes - structure, function and impact on disease. Biochim. Biophys. Acta. Bioenerg 1862 (1), 148317. doi:10.1016/j.bbabio.2020.148317
Blaiseau, P.-L., Lesuisse, E., and Camadro, J.-M. (2001). Aft2p, a novel iron-regulated transcription activator that modulates, with Aft1p, intracellular iron use and resistance to oxidative stress in yeast. J. Biol. Chem 276, 34221–34226. doi:10.1074/jbc.M104987200
Brault, A., and Labbé, S. (2020). Iron deficiency leads to repression of a non-canonical methionine salvage pathway in Schizosaccharomyces pombe. Mol. Microbiol 114, 46–65. doi:10.1111/mmi.14495
Brault, A., Mourer, T., and Labbé, S. (2015). Molecular basis of the regulation of iron homeostasis in fission and filamentous yeasts. IUBMB life 67, 801–815. doi:10.1002/iub.1441
Brault, A., Rallis, C., Normant, V., Garant, J.-M., Bahler, J., and Labbe, S. (2016). Php4 is a key player for iron economy in meiotic and sporulating cells. G3 (Bethesda) 6 (10), 3077–3095. doi:10.1534/g3.116.031898
Chao, L. Y., Marletta, M. A., and Rine, J. (2008). Sre1, an iron-modulated GATA DNA-binding protein of iron-uptake genes in the fungal pathogen Histoplasma capsulatum. Biochemistry 47, 7274–7283. doi:10.1021/bi800066s
Chen, G.-H., Song, C.-C., Pantopoulos, K., Wei, X.-L., Zheng, H., and Luo, Z. (2022). Mitochondrial oxidative stress mediated Fe-induced ferroptosis via the NRF2-ARE pathway. Free Radic. Biol. Med 180, 95–107. doi:10.1016/j.freeradbiomed.2022.01.012
Chung, W.-H., Kim, K.-D., and Roe, J.-H. (2005). Localization and function of three monothiol glutaredoxins in Schizosaccharomyces pombe. Biochem. Biophys. Res. Commun 330, 604–610. doi:10.1016/j.bbrc.2005.02.183
Cunningham, T. S., and Cooper, T. G. (1993). The Saccharomyces cerevisiae DAL80 repressor protein binds to multiple copies of GATAA-containing sequences (URSGATA). J. Bacteriol 175 (18), 5851–5861. doi:10.1128/jb.175.18.5851-5861.1993
Daniel, T., Faruq, H. M., Laura Magdalena, J., Manuela, G., and Christopher Horst, L. (2020). Role of GSH and iron-sulfur glutaredoxins in iron metabolism—review. Molecules 25, 3860. doi:10.3390/molecules25173860
Dlouhy, A. C., Beaudoin, J., Labbé, S., and Outten, C. E. (2017). Schizosaccharomyces pombe Grx4 regulates the transcriptional repressor Php4 via [2Fe–2S] cluster binding. Metallomics 9, 1096–1105. doi:10.1039/c7mt00144d
El-Sappah, A. H., Elbaiomy, R. G., Elrys, A. S., Wang, Y., Zhu, Y., and Huang, Q. (2021). Genome-wide identification and expression analysis of metal tolerance protein gene family in Medicago truncatula under a broad range of heavy metal stress. Front. Genet 12, 713224. doi:10.3389/fgene.2021.713224
El-Sappah, A. H., Elrys, A. S., Desoky, E.-S. M., Zhao, X., Bingwen, W., and El-Sappah, H. H. (2021). Comprehensive genome wide identifica-tion and expression analysis of MTP gene family in tomato (Solanum lycopersicum) under multiple heavy metal stress. Saudi J. Biol. Sci 28, 6946–6956. doi:10.1016/j.sjbs.2021.07.073
Encinar Del Dedo, J., Gabrielli, N., Carmona, M., Ayte, J., and Hidalgo, E. (2015). A cascade of iron-containing proteins governs the genetic iron starvation response to promote iron uptake and inhibit iron storage in fission yeast. PLoS Genet 11, e1005106. doi:10.1371/journal.pgen.1005106
Gupta, M., and Outten, C. E. (2020). Iron–sulfur cluster signaling: The common thread in fungal iron regulation. Curr. Opin. Chem. Biol 55, 189–201. doi:10.1016/j.cbpa.2020.02.008
Haas, H., Eisendle, M., and Turgeon, B. G. (2008). Siderophores in fungal physiology and virulence. Annu. Rev. Phytopathol 46, 149–187. doi:10.1146/annurev.phyto.45.062806.094338
Herrero, E., and De La Torre-Ruiz, M. A. (2007). Monothiol glutaredoxins: A common domain for multiple functions. Cell. Mol. Life Sci 64, 1518–1530. doi:10.1007/s00018-007-6554-8
Holmgren, A. (1989). Thioredoxin and glutaredoxin systems. J. Biol. Chem 264, 13963–13966. doi:10.1016/s0021-9258(18)71625-6
Hortschansky, P., Eisendle, M., Al-Abdallah, Q., Schmidt, A. D., Bergmann, S., and Thön, M. (2007). Interaction of HapX with the CCAAT-binding complex—A novel mechanism of gene regulation by iron. EMBO J 26, 3157–3168. doi:10.1038/sj.emboj.7601752
Hwang, L. H., Seth, E., Gilmore, S. A., and Sil, A. (2012). SRE1 regulates iron-dependent and-independent pathways in the fungal pathogen Histoplasma capsulatum. Eukaryot. Cell 11, 16–25. doi:10.1128/EC.05274-11
Islam, M., Qi, S., Zhang, S., Amin, B., Yadav, V., and El-Sappah, A. H. (2022). Genome-wide identification and functions against tomato spot-ted wilt tospovirus of PR-10 in Solanum lycopersicum. Int. J. Mol. Sci 23, 1502. doi:10.3390/ijms23031502
Jbel, M., Mercier, A., and Labbé, S. (2011). Grx4 monothiol glutaredoxin is required for iron limitation-dependent inhibition of Fep1. Eukaryot. Cell 10, 629–645. doi:10.1128/EC.00015-11
Jbel, M., Mercier, A., Pelletier, B., Beaudoin, J., and Labbé, S. (2009). Iron activates in vivo DNA binding of Schizosaccharomyces pombe transcription factor Fep1 through its amino-terminal region. Eukaryot. Cell 8, 649–664. doi:10.1128/EC.00001-09
Jung, W. H., Saikia, S., Hu, G., Wang, J., Fung, C. K.-Y., and D'souza, C. (2010). HapX positively and negatively regulates the transcriptional response to iron deprivation in Cryptococcus neoformans. PLoS Pathog 6, e1001209. doi:10.1371/journal.ppat.1001209
Jung, W. H., Sánchez-León, E., and Kronstad, J. W. (2021). Coordinated regulation of iron metabolism in cryptococcus neoformans by GATA and CCAAT transcription factors: Connections with virulence. Curr. Genet. 67, 583–593. doi:10.1007/s00294-021-01172-5
Katsarou, A., and Pantopoulos, K. (2020). Basics and principles of cellular and systemic iron homeostasis. Mol. Aspects Med. 75, 100866. doi:10.1016/j.mam.2020.100866
Kelley, L., Mezulis, S., Yates, C., Wass, M., and Sternberg, M. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc 10, 845–858. doi:10.1038/nprot.2015.053
Kim, H.-G., Kim, B.-C., Park, E.-H., and Lim, C.-J. (2005). Stress-dependent regulation of a monothiol glutaredoxin gene from Schizosaccharomyces pombe. Can. J. Microbiol 51, 613–620. doi:10.1139/w05-034
Kim, K.-D., Kim, H.-J., Lee, K.-C., and Roe, J.-H. (2011). Multi-domain CGFS-type glutaredoxin Grx4 regulates iron homeostasis via direct interaction with a repressor Fep1 in fission yeast. Biochem. Biophys. Res. Commun 408, 609–614. doi:10.1016/j.bbrc.2011.04.069
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). Mega X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol 35, 1547–1549. doi:10.1093/molbev/msy096
Labbé, S., Khan, M. G. M., and Jacques, J.-F. (2013). Iron uptake and regulation in Schizosaccharomyces pombe. Curr. Opin. Microbiol 16, 669–676. doi:10.1016/j.mib.2013.07.007
Labbé, S., Pelletier, B., and Mercier, A. (2007). Iron homeostasis in the fission yeast Schizosaccharomyces pombe. Biometals 20, 523–537. doi:10.1007/s10534-006-9056-5
Lascaris, R., Bussemaker, H. J., Boorsma, A., Piper, M., van der Spek, H., Grivell, L., et al. (2018). Hap4p overexpression in glucose-grown Saccharomyces cerevisiae induces cells to enter a novel metabolic state. Genome Biol. 4 (1), R3. doi:10.1186/gb-2002-4-1-r3
Letunic, I., and Bork, P. (2018). 20 years of the SMART protein domain annotation resource. Nucleic Acids Res 46, D493-D496–D496. doi:10.1093/nar/gkx922
Li, J., Zhang, L., Elbaiomy, R. G., Chen, L., Wang, Z., and Jiao, J. (2022). Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae). PeerJ 10, e12880. doi:10.7717/peerj.12880
Li, J., Zhang, L., Yuan, Y., Wang, Q., Elbaiomy, R. G., and Zhou, W. (2021). In silico functional prediction and expression analysis of C2H2 zinc-finger family transcription factor revealed regulatory role of ZmZFP126 in maize growth. Front. Genet 12, 770427. doi:10.3389/fgene.2021.770427
Liu, Z., Ebrahim, A., Wu, X., Li, M., and Huang, Y. (2022). Loss of PPR protein Ppr2 induces ferroptosis-like cell death in Schizosaccharomyces pombe. Arch. Microbiol 204, 360–369. doi:10.1007/s00203-022-02970-2
Maio, N., Zhang, D., Ghosh, M. C., Jain, A., SantaMaria, A. M., and Rouault, T. A. (2022). Mechanisms of cellular iron sensing, regulation of erythropoiesis and mitochondrial iron utilization. Semin. Hematol 58, 161–174. doi:10.1053/j.seminhematol.2021.06.001
Mao, Y., and Chen, C. (2019). The hap complex in yeasts: Structure, assembly mode, and gene regulation. Front. Microbiol 10, 1645. doi:10.3389/fmicb.2019.01645
Martínez-Pastor, M. T., and Puig, S. (2020). Adaptation to iron deficiency in human pathogenic fungi. Biochim. Biophys. Acta. Mol. Cell Res 1867, 118797. doi:10.1016/j.bbamcr.2020.118797
Mcnabb, D. S., and Pinto, I. (2005). Assembly of the hap2p/hap3p/hap4p/hap5p-DNA complex in Saccharomyces cerevisiae. Eukaryot. Cell 4, 1829–1839. doi:10.1128/EC.4.11.1829-1839.2005
Mercier, A., and Labbé, S. (2009). Both Php4 function and subcellular localization are regulated by iron via a multistep mechanism involving the glutaredoxin Grx4 and the exportin Crm1. J. Biol. Chem 284, 20249–20262. doi:10.1074/jbc.M109.009563
Mercier, A., Pelletier, B., and Labbé, S. (2006). A transcription factor cascade involving Fep1 and the CCAAT-binding factor Php4 regulates gene expression in response to iron deficiency in the fission yeast Schizosaccharomyces pombe. Eukaryot. Cell 5, 1866–1881. doi:10.1128/EC.00199-06
Mercier, A., Watt, S., BäHler, J. R., and Labbé, S. (2008). Key function for the CCAAT-binding factor Php4 to regulate gene expression in response to iron deficiency in fission yeast. Eukaryot. Cell 7, 493–508. doi:10.1128/EC.00446-07
Miele, A. E., Cervoni, L., Le Roy, A., Cutone, A., Musci, G., and Ebel, C. (2021). Biophysical characterization of the complex between the iron-responsive transcription factor Fep1 and DNA. Eur. Biophys. J 50 (3-4), 501–512. doi:10.1007/s00249-020-01489-y
Misslinger, M., Hortschansky, P., Brakhage, A. A., and Haas, H. (2021). Fungal iron homeostasis with a focus on Aspergillus fumigatus. Biochim. Biophys. Acta. Mol. Cell Res. 1868, 118885. doi:10.1016/j.bbamcr.2020.118885
Molina, M. M., Bellí, G., De La Torre, M. A., Rodríguez-Manzaneque, M. T., and Herrero, E. (2004). Nuclear monothiol glutaredoxins of Saccharomyces cerevisiae can function as mitochondrial glutaredoxins. J. Biol. Chem 279, 51923–51930. doi:10.1074/jbc.M410219200
Mühlenhoff, U., Braymer, J. J., Christ, S., Rietzschel, N., Uzarska, M. A., and Weiler, B. D. (2020). Glutaredoxins and iron-sulfur protein biogenesis at the interface of redox biology and iron metabolism. Biol. Chem 401 (12), 1407–1428. doi:10.1515/hsz-2020-0237
Ohtsuka, H., Shimasaki, T., and Aiba, H. (2021). Genes affecting the extension of chronological lifespan in Schizosaccharomyces pombe (fission yeast). Mol. Microbiol 115, 623–642. doi:10.1111/mmi.14627
Palavecino-Ruiz, M., Bermudez-Moretti, M., and Correa-Garcia, S. (2017). Unravelling the transcriptional regulation of Saccharomyces cerevisiae UGA genes: The dual role of transcription factor Leu3. Microbiology 163 (11), 1692–1701. doi:10.1099/mic.0.000560
Pelletier, B., Beaudoin, J., Mukai, Y., and Labbé, S. (2002). Fep1, an iron sensor regulating iron transporter gene expression in Schizosaccharomyces pombe. J. Biol. Chem 277, 22950–22958. doi:10.1074/jbc.M202682200
Peng, S., Cheng, M., Huang, K., Cui, Y., Zhang, Z., and Guo, R. (2018). Efficient computation of motif discovery on intel many integrated core (mic) architecture. BMC Bioinforma 19, 282–110. doi:10.1186/s12859-018-2276-1
Pfaffle, M. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 29, e45–e2007. doi:10.1093/nar/29.9.e45
Protacio, R. U., Mukiza, T. O., Davidson, M. K., and Wahls, W. P. (2022). DNA sequence-mediated, evolutionarily rapid redistribution of meiotic recombination hotspots. Genetics 220, 685–694. iyab212. doi:10.1534/genetics.111.134130
Robinson, J. R., Isikhuemhen, O. S., and Anike, F. N. (2021). Fungal–metal interactions: A review of toxicity and homeostasis. J. Fungi 7, 225. doi:10.3390/jof7030225
Talib, E. A., and Outten, C. E. (2021). Iron-sulfur cluster biogenesis, trafficking, and signaling: Roles for CGFS glutaredoxins and BolA proteins. Biochim. Biophys. Acta. Mol. Cell Res 1868, 118847. doi:10.1016/j.bbamcr.2020.118847
Vachon, P., Mercier, A., Jbel, M., and Labbé, S. (2012). The monothiol glutaredoxin Grx4 exerts an iron-dependent inhibitory effect on Php4 function. Eukaryot. Cell 11, 806–819. doi:10.1128/EC.00060-12
Keywords: iron homeostasis, Schizosaccharomyces pombe, Grx4, Fep1, Php4
Citation: Ebrahim A, Alfwuaires MA, Abukhalil MH, Alasmari F, Ahmad F, Yao R, Luo Y and Huang Y (2022) Schizosaccharomyces pombe Grx4, Fep1, and Php4: In silico analysis and expression response to different iron concentrations. Front. Genet. 13:1069068. doi: 10.3389/fgene.2022.1069068
Received: 13 October 2022; Accepted: 16 November 2022;
Published: 07 December 2022.
Edited by:
Ahmed H. El-Sappah, Zagazig University, EgyptReviewed by:
Huang Qiulan, Chengdu Institute of Biology (CAS), ChinaAli B. Roomi, Al-Ayen University, Iraq
Ahmad Ali, Jiangsu University, China
Copyright © 2022 Ebrahim, Alfwuaires, Abukhalil, Alasmari, Ahmad, Yao, Luo and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ying Huang, yhuang@njnu.edu.cn
†ORCID: Ying Huang, orcid.org/0000-0001-8469-5185