- 1Department of Cardiac Surgery, Xiangya Hospital, Central South University, Changsha, China
- 2Department of Endocrinology, Xiangya Hospital, Central South University, Changsha, China
- 3The Institute of Medical Science Research, Xiangya Hospital, Central South University, Changsha, China
- 4Department of Cardiology, Xiangya Hospital, Central South University, Changsha, China
- 5National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
Background: Atrial fibrillation (AF) increases the risk of stroke and heart failure. Postoperative AF (POAF) increases the risk of mortality after cardiac surgery. This study aims to explore mechanisms underlying AF, analyze infiltration of immune cells in left atrium (LA) from patients with AF, and identify potential circular RNA (circRNA) biomarkers for POAF.
Methods: Raw data of GSE797689, GSE115574, and GSE97455 were downloaded and processed. AF-related gene co-expression network was constructed using weighted gene correlation network analysis and enrichment analysis of genes in relevant module was conducted. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were applied to investigate pathways significantly enriched in AF group. Infiltration of immune cells was analyzed using single-sample GSEA. Differentially expressed genes (DEGs) between patients with or without AF were identified and competing endogenous RNA (ceRNA) networks of DEGs were constructed. To screen biomarkers for POAF, differentially expressed circRNAs (DEcircRNAs) between patients with or without POAF were identified. Intersection between DEcircRNAs and circRNAs in ceRNA networks of DEGs were extracted and circRNAs in the intersection were further screened using support vector machine, random forest, and neural network to identify biomarkers for POAF.
Results: Three modules were found to be relevant with AF and enrichment analysis indicated that genes in these modules were enriched in synthesis of extracellular matrix and inflammatory response. The results of GSEA and GSVA suggested that inflammatory response-related pathways were significantly enriched in AF group. Immune cells like macrophages, mast cells, and neutrophils were significantly infiltrated in LA tissues from patients with AF. The expression levels of immune genes such as CHGB, HLA-DRA, LYZ, IGKV1-17 and TYROBP were significantly upregulated in patients with AF, which were correlated with infiltration of immune cells. ceRNA networks of DEGs were constructed and has_circ_0006314 and hsa_circ_0055387 were found to have potential predictive values for POAF.
Conclusion: Synthesis of extracellular matrix and inflammatory response were main processes involved in development and progression of AF. Infiltration of immune cells was significantly different between patients with or without AF. Has_circ_0006314 and hsa_circ_0055387 were found to have potential predictive values for POAF.
Introduction
Atrial fibrillation (AF) is one of the most common cardiac arrhythmia diseases, which affects 2%–3% of the Western population (Heijman et al., 2018). Traditional risk factors of AF include older age, coronary artery disease, male sex, hypertension, obesity, diabetes mellitus, and obstructive sleep apnea (Staerk et al., 2017). And there are studies revealed that AF exhibited a genetic predisposition (Palatinus and Das, 2015; Tucker et al., 2016). AF could increase the risk of stroke, heart failure, and dementia (Michaud and Stevenson, 2021). These severer complications of AF significantly increased the mortality by a factor of 2.4 among men and by a factor of 3.5 among women (Benjamin et al., 1998).
Structural remodeling and electrical remodeling are the key pathogenesis of AF. Development of atrial fibrosis and enlargement of left atrium are the hallmark of structural remodeling and atrial fibrosis is considered to be the substrate for AF initiation and maintenance (Dzeshka et al., 2015). Atrial fibrosis is characterized by progressive accumulation of extracellular matrix (ECM) in myocardium, which is predominantly comprised of collagen (Ueha et al., 2012). Cardiac fibroblasts (CFs) account for up to 60% of cells in cardiac muscle and play pivotal roles in synthesis and formation of ECM (Souders et al., 2009). When activated by stimulators of collagen synthesis like transforming growth factor beta-1 (TGF-β1), CFs differentiate into myofibroblasts which have approximate 2-fold higher capacity to synthesize collagen than CFs (Baum and Duffy, 2011). The differentiation of CFs and the synthesis of collagen eventually lead to atrial fibrosis. Electrical remodeling of AF is associated with the dysregulation of cardiac ion channels such as Ca2+ and K+ channels (Christ et al., 2004). The altered expression or distribution of connexin-40 and connexin-43 could result in heterogeneous conduction and hampering cell to cell communication (Yan et al., 2018). In addition, inflammation is closely related to AF initiation and maintenance. Macrophages and T cells were found infiltrating in the atrium of patients with AF (Yamashita et al., 2010). Researches indicated that the increasements of inflammatory cytokines like C-reactive protein (CRP), interleukin (IL)-1β, IL-6, and tumor necrosis factor (TNF)-α were correlated with the progression of AF (Li and Brundel, 2020). Furthermore, inflammation-related transcript factors like nuclear factor kappa B (NF-κB) could upregulate the expression level of TGF-β1 and slow conduction by decreasing the expression of sodium channel, both of which accelerate the progression of AF (Ziolo and Mohler, 2015; Faria and Persaud, 2017). But the exact mechanism underlying AF initiation and maintenance is largely unknown and needs further exploration.
Postoperative AF (POAF) is characterized by new onset AF in the immediate period after cardiac or non-cardiac surgery, which increases the intensive care unit time and in hospital cost (Greenberg et al., 2017). In this study, we focused on POAF after cardiac surgery. POAF is an important complication of coronary artery bypass grafting and valvular surgery. And the prevalence of POAF ranges between 20% and 40% in different studies (Dobrev et al., 2019). According to recent researches, older age, male sex, congestive heart failure, hypertension, obesity, and cross-clamp time are independently associated with POAF (Dobrev et al., 2019). POAF significantly increases the risk of stroke, subsequent AF recurrence, and cardiovascular death (Ahlsson et al., 2010). Animal experiments suggested that dysregulation of autonomic nervous system, inflammation, and oxidative stress were related with POAF (Dobrev et al., 2019). However, the mechanism of POAF remains unclear. Therefore, it is hard to predict and prevent POAF.
This study aimed to construct AF-related gene co-expression networks, analyze the infiltration of immune cells in AF, construct competing endogenous RNA (ceRNA) networks of differentially expressed genes (DEGs), and identify potential biomarkers for prediction of POAF using bioinformatics.
Materials and methods
Data acquisition and processing
The analysis process of this study is shown in Figure 1. The raw data of GSE79768 (including six left atrium samples from patients with sinus rhythm (SR) and seven left atrium samples from patients with AF) and GSE115574 (including 15 left atrium samples from patients with SR and 14 left atrium samples from patients with AF) were downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds). These datasets were based on the same platform, GPL570 (Affymetrix Human Genome U133 Plus 2.0 Array). The expression matrix was extracted using “affy” package in R software (version 4.1.2; R Foundation for Statistical Computing, Vienna, Austria) and the data were annotated using the “dplyr” and “limma” packages in R software (version 4.1.2). After annotation, the data in expression matrix were normalized using “limma” package in R software (version 4.1.2). Then, the expression matrix of two datasets was merged and the batch effect between two datasets was removed using the “sva” package in R software (version 4.1.2). The merged expression matrix was used for further analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g001.jpg)
FIGURE 1. Flow chart diagram for the analysis process in this study. SR: sinus rhythm, AF: atrial fibrillation, POAF: postoperative atrial fibrillation, WGCNA: weighted gene correlation network analysis, GSEA: gene set enrichment analysis, GSVA: gene set variation enrichment analysis, ceRNA: competing endogenous RNA.
The raw data of GSE97455 (including 15 plasma samples from patients without POAF and 15 plasma samples from patients with POAF) were downloaded from GEO database (https://www.ncbi.nlm.nih.gov/gds). In GSE97455, plasma samples were collected from patients who were going to undergo coronary artery bypass grafting surgery. The time of collection was before the start of the surgery and the occurrence of AF was recorded after surgery. The data were annotated using the “dplyr” and “limma” packages in R software (version 4.1.2) and then normalized using “limma” package in R software (version 4.1.2). Finally, the names of circular RNA (circRNA) in the dataset were transformed into the format of circRNA ID in circBase (http://www.circbase.org/) for further analysis.
Gene co-expression network construction by weighted gene correlation network analysis (WGCNA)
Gene co-expression networks were constructed using the WGCNA package in R software (version 4.1.2) (Langfelder and Horvath, 2008). Soft-thresholding power was used to construct a weighted adjacency matrix. Relationships between a single gene and others in the analysis were incorporated, and the adjacency matrix was transformed into the topological matrix (TOM). Then, a hierarchical clustering analysis of genes was performed using 1—TOM as the distance measure. Thereafter, modules were detected using a dynamic tree cut algorithm with a minimum module size of 50 and a minimum cut height of 0.99. The correlation between each module and AF was calculated and shown in a heatmap. Finally, the most relevant gene modules were selected for further analysis.
Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis of genes in the relevant gene modules
The Database for Annotation, Visualization and Integrated Discovery (DAVID, 2021 Update, https://david.ncifcrf.gov/home.jsp) were used to conduct GO and KEGG pathways enrichment analysis of genes in the relevant gene modules.
Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA)
The gene set files used in this study were downloaded from the Molecular Signatures Database version 7.5.1 (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp). The enrichment scores of GO and KEGG pathways terms in each group were calculated using the GSEA software (version 4.2.3), and terms enriched in the AF group were identified. A nominal p value of <0.05 and false-discovery rate q value of <0.25 were considered as significantly enriched in the AF group.
GSVA was applied to evaluate GO and KEGG pathway terms enriched in each sample by converting the gene expression matrix into a gene set expression matrix using the “GSVA” package in R (version 4.1.2). After that, the differentially enriched terms between two groups were identified using R (version 4.1.2) with the threshold of p < 0.05. The differentially enriched terms were visualized using the “pheatmap” package in R (version 4.1.2).
Single-sample gene set enrichment analysis (ssGSEA)
The gene set files used in this study were downloaded from the Molecular Signatures Database version 7.5.1 (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp). The scores of 16 immune cells and four cardiac cells were calculated using “gsva” package in R (version 4.1.2).
Identification of differentially expressed genes (DEGs)
DEGs and differentially expressed circRNA (DEcircRNA) between two groups were identified using R (version 4.1.2) with the threshold of p < 0.05 and absolute value of log fold change (FC) > 0.5. The immune and inflammation-related genes (IRGs) list was downloaded from Immport (https://www.immport.org/home). And differentially expressed IRGs (DEIRGs) were identified using R (version 4.1.2) with the threshold of p < 0.05 and absolute value of log fold change (FC) > 0.2.
Construction of competing endogenous RNA (ceRNA) networks of DEGs
The Encyclopedia of RNA Interactomes (ENCORI, https://starbase.sysu.edu.cn/) (Li et al., 2014) database was used to construct ceRNA networks of DEGs. The micro-RNA (miRNA) which could interact with the DEGs was predicted using ENCORI with the threshold of number of supporting crosslinking-immunoprecipitation and high-throughput sequencing (CLIP-seq) experiments≥2 and number of target-predicting programs≥2. Then, the interaction pairs of miRNA-circRNA were predicted using ENCORI with the threshold of number of supporting CLIP-seq experiments≥5 and supporting degradome-seq experiments≥1. Finally, the circRNA-miRNA-mRNA interaction network was visualized using Cytoscape software (version 3.9.1).
Identification of candidate circRNA biomarkers for POAF
The intersection between circRNA in ceRNA network and DEcircRNA was extracted. Upregulated circRNA which is associated with upregulated DEGs was selected as candidate biomarker for POAF. Similarly, downregulated circRNA which is associated with downregulated DEGs was selected as candidate biomarker for POAF.
Construction of predictive model for POAF and identification of biomarkers for POAF
The candidate circRNAs were selected as features to construct predictive models. Support vector machine (SVM), random forest (RF), and neural network were applied to construct predictive models for POAF using “e1071”, “randomForest”, and “neuralnet” packages in R (version 4.1.2). Feature importance of each predictive model was analyzed and ranked. Thereafter, the intersection of top five circRNA in each predictive model was identified as potential biomarkers for POAF.
Statistical analysis
The relative expression levels of mRNA were presented as mean ± standard deviation (SD). The difference in scores of immune cells and cardiac cells between different groups was analyzed by the Mann-Whitney U test. Pearson’s correlation analysis was used to analyze the relationship between inflammation genes and the scores of immune cells and cardiac cells. A value of p < 0.05 was considered to be statistically significant. Statistical analyses were performed using SPSS version 19 (IBM Corporation, Armonk, NY, United States) and R (version 4.1.2).
Results
Construction of AF-related gene co-expression networks and enrichment analysis of genes in the relevant modules
AF-related gene co-expression networks were constructed using WGCNA and six gene modules were obtained (Figure 2A). The correlation between each gene module and the occurrence of AF was analyzed, and coefficients and p value are shown in Figure 2B. The absolute values of coefficients of green, blue, and brown modules were close, therefore, enrichment analysis of genes in these modules were conducted, respectively. Genes in green module were significantly enriched in regulation of bone morphogenetic protein (BMP) signaling pathway, negative regulation of DNA recombination, and protein digestion and absorption (Figures 2C,D). The enrichment analysis of genes in blue module revealed that the top three terms in the GO biological process (BP) subdivision were inflammatory response, immune response, and detoxification of copper ion, while the top three KEGG pathways were tumor necrosis factor (TNF) signaling pathway, complement and coagulation cascades, and IL-17 signaling pathway (Figures 2E,F). The enrichment analysis of genes in brown module indicated that the top three terms in GO BP subdivision were collagen fibril organization, extracellular matrix organization, and cell adhesion, while the top three KEGG pathways were complement and coagulation cascades, ECM-receptor interaction, and focal adhesion (Figures 2G,H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g002.jpg)
FIGURE 2. Construction of AF-related gene co-expression network and the enrichment analysis of genes in the AF-related gene modules. (A) Dendrogram and clustering for identification of gene co-expression modules. (B) Correlation analysis of gene co-expression modules with AF. The numbers above brackets were correlation coefficients and the numbers in brackets were p values. (C,D) GO-BP and KEGG enrichment analysis of genes in green module. (E,F) GO-BP and KEGG enrichment analysis of genes in blue module. (G,H) GO-BP and KEGG enrichment analysis of genes in brown module. AF: atrial fibrillation, GO-BP: gene ontology biological process, KEGG: Kyoto Encyclopedia of Genes and Genomes.
GSEA and GSVA
GSEA was conducted to further explore the terms of GO BP and pathways enriched in AF group. As shown in Figures 3A–C, the terms of GO BO like B cell receptor signaling pathway, immune response regulating cell surface receptor signaling pathway, and cellular defense response were significantly enriched in AF group. In terms of KEGG pathway, FC gamma R mediated phagocytosis, B cell receptor signaling pathway, and FC epsilon RI signaling pathway were significantly enriched in AF group (Figures 3D–F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g003.jpg)
FIGURE 3. The results of GSEA and GSVA. (A–C) The results of GSEA revealed the GO-BP terms significantly enriched in AF group. (D–F) The results of GSEA revealed the KEGG pathways significantly enriched in AF group. (G) The variations of GO-BP terms between different samples were analyzed using GSVA and the top 20 significantly different GO-BP terms were shown in the heatmap. (H) The variations of KEGG pathways between different samples were analyzed using GSVA and the significantly different KEGG pathways were shown in the heatmap. GSEA: gene set enrichment analysis, GSVA: gene set variation analysis, GO-BP: gene ontology biological process, KEGG: Kyoto Encyclopedia of Genes and Genomes, AF: atrial fibrillation.
GSVA was utilized to analyze the variations of pathway between different samples and the differentially activated pathways were identified. The top 10 upregulated and downregulated terms of GO BP were displayed in Figure 3G. The results indicated that terms like Toll like receptor seven signaling pathway and FC epsilon receptor signaling pathway were significantly activated in AF group while detoxification of copper ion and stress response to metal ion were significantly downregulated in AF group (Figure 3G). The heatmap of variations of KEGG pathway between different samples was shown in Figure 3H suggesting the activation of FC gamma R mediated phagocytosis and FC epsilon RI signaling pathway in samples with AF, which is consistent with the results of GSEA.
The infiltration of immune cells in different samples and the correlation between immune-related genes and immune cells
The infiltration of immune cells in different samples was investigated using ssGSEA and the results showed that the scores of macrophages, mast cells, neutrophils, tumor infiltrating lymphocyte (TIL), and regulatory T cells (Treg) were significantly higher in AF group, while the score of immature dendritic cells (iDCs) was significantly decreased in AF group (Figure 4A). In addition, scores of four cardiac cells were calculated using ssGSEA, and the results revealed that the score of epicardial cell was significantly decreased in AF group, while the score of cardiac fibroblasts was significantly increased in AF group. Then, the differentially expressed immune-related genes (DEIRGs) between SR group and AF group were identified (Figure 4C) and the top five upregulated DEIRGs were CHGB, HLA-DRA, LYZ, IGKV1-17, and TYROBP (Figure 4D). The correlation between DEIRGs and different types of cells was determined and shown in Figure 4E.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g004.jpg)
FIGURE 4. The analysis of immune cells and cardiac cells in LA samples and the correlation between immune genes, immune cells, and cardiac cells. (A) The scores of immune cells in LA samples. (B) The scores of cardiac cells in LA samples. (C) The volcano plot of differentially expressed immune genes between patients with or without AF. Dots in green: downregulated gene, dots in red: upregulated gene. (D) The expression levels of top five upregulated immune genes in SR group and AF group. (E) The correlation between immune genes, immune cells, and cardiac cells. Dots in blue: SR samples, dots in red: AF samples. LA: left atrium, SR: sinus rhythm, AF: atrial fibrillation. *, p < 0.05, **, p < 0.01, ***, p < 0.001.
Construction of ceRNA networks of DEGs
The DEGs between SR group and AF group were analyzed and the top 10 upregulated and downregulated genes were shown in Table 1. Out of top 10 downregulated DEGs, ceRNA networks of BCHE, FAM110C, KRT18, and LRRC49 were constructed successfully (Figure 5A). The prediction results of other six genes showed that there was no interaction pair of miRNA-mRNA according to the threshold of prediction. Similarly, out of top 10 upregulated DEGs, ceRNA networks of FHL2, RELN, and LBH were constructed successfully and the ceRNA network of FHL2 was displayed in Figure 5B as representative network. miRNA which could interact with other seven upregulated genes was not found with the threshold of prediction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g005.jpg)
FIGURE 5. ceRNA networks of DEGs. (A) The ceRNA networks of downregulated DEGs. (B) The ceRNA networks of upregulated DEGs, FHL2. ceRNA: competing endogenous RNA, DEGs: differentially expressed genes.
Screening of candidate circRNA biomarkers for POAF
The differentially expressed circRNA (DEcircRNA) between no-POAF group and POAF group was identified with the threshold of p < 0.05 and absolute value of log fold change (FC) > 0.5 (Figure 6A). There were 210 upregulated circRNA and 293 downregulated circRNA. The intersection between upregulated circRNA and the circRNA in ceRNA networks of upregulated DEGs was extracted as candidate circRNA biomarkers for POAF, and the intersection between downregulated circRNA and the circRNA in ceRNA networks of downregulated DEGs was extracted as candidate circRNA biomarkers for POAF, too (Figure 6B). Finally, 10 candidate circRNA biomarkers for POAF were obtained and the expression levels of these circRNA in different groups of GSE97455 were shown in Figure 6C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g006.jpg)
FIGURE 6. Screening of candidate circRNA biomarkers for POAF. (A) The volcano plot of differentially expressed circRNAs between patients with or without POAF. Dots in green: downregulated circRNAs, dots in red: upregulated circRNAs. (B) The ceRNA networks of DEGs after screening. The label in red: upregulated DEGs or differentially expressed circRNAs, the label in blue: downregulated DEGs or differentially expressed circRNAs, the label in green: miRNAs. (C) The expression levels of candidate circRNA biomarkers for POAF in patients with or without POAF. POAF: postoperative atrial fibrillation, DEGs: differentially expressed genes. *, p < 0.05, **, p < 0.01, ***, p < 0.001.
Construction of predictive models for POAF and identification of biomarkers for POAF
10 candidate circRNA biomarkers for POAF were selected as features of predictive models. SVM, RF, and neural network were applied to construct predictive models for POAF based on the data of GSE97455. In order to achieve highest accuracy with least features in model, the algorithm of SVM chose hsa_circ_0006314, hsa_circ_0001349, hsa_circ_0055387, hsa_circ_0003069, and hsa_circ_0002211 out of 10 candidate circRNA biomarkers as features of predictive model (accuracy = 1). The algorithm of RF used 10 candidate circRNA biomarkers as features of predictive model (accuracy = 0.933) and the rank of feature importance was shown in Figure 7A. Top five features in the RF model were hsa_circ_0055387, hsa_circ_0006314, hsa_circ_0092283, hsa_circ_0001349, and hsa_circ_0006623. Similarly, the algorithm of neural network selected 10 candidate circRNA biomarkers to construct predictive model (accuracy = 1). The importance of features in neural network model was analyzed and the top five features were hsa_circ_0006314, hsa_circ_0041631, hsa_circ_0092283, hsa_circ_0055387, and hsa_circ_0004470 (Figure 7B). Finally, the intersection between top five features of three models was extracted. Hsa_circ_0006314 and hsa_circ_0055387 were identified as potential biomarkers for POAF.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1003366/fgene-13-1003366-HTML/image_m/fgene-13-1003366-g007.jpg)
FIGURE 7. Identification of biomarkers for POAF. (A) The rank of feature importance in RF predictive model for POAF. (B) The rank of feature importance in neural network predictive model for POAF. (C) The intersection between top five features of each model was shown in the Venn diagram. POAF: postoperative atrial fibrillation, RF: random forest.
Discussion
AF is the most common cardiac arrhythmia disease which is characterized by irregular and abnormally fast contraction of atrial cardiomyocytes (Hindricks et al., 2020). Symptoms of AF include irregular heart rate, palpitation, dizziness, and tiredness (Hindricks et al., 2020). Moreover, AF could lead to severer complications like stroke, heart failure, and cognitive impairment, both of which are associated with the increasement of morbidity and mortality (Zoni-Berisso et al., 2014). According to the epidemiological study, the prevalence of AF would keep increasing in the next 3 decades (Colilla et al., 2013). It was suggested that by the year 2050, approximately 5.2 million men and 3.1 million women with age older than 60 would suffer from AF in China (Zhang et al., 2021). AF is a serious threat to patients. Therefore, it is necessary to further explore the mechanisms underlying AF initiation and maintenance and find out potential therapeutic targets for AF.
As mentioned above, atrial structural remodeling is the key substrate of AF. Left atrium (LA) and right atrium (RA) exhibited different susceptibilities in AF initiation and maintenance (Tsai et al., 2016). Research revealed that the high frequency reentrant sources (or rotors) were more distributed in LA (Kalifa et al., 2006). Besides, Sahadevan et al. reported that the dominant frequency of AF was present in LA more frequently (Sahadevan et al., 2004). Moreover, Lin et al. demonstrated that complex fractionated atrial electrograms were more extensively present in LA (Lin et al., 2009). There were other researches suggested that LA is more important than RA in AF maintenance (Haissaguerre et al., 2008; Oral et al., 2008). In consideration of the importance of LA in AF initiation and maintenance, gene expression data of LA samples from patients with or without AF were extracted from GSE79768 and GSE115574 to investigate the potential mechanisms of AF. After construction of AF-related gene co-expression networks, the results of correlation analysis between each gene module and the presence of AF indicated that the green, blue, and brown modules were significantly relevant to AF. Genes in green module were significantly enriched in regulation of BMP signaling pathway. BMPs are members of TGF-β family and involve in the control of cell behavior, proliferation, and differentiation (Kraj, 2022). Chen et al. reported that BMP7 could inhibit atrial fibrosis via mediating TGF-β1/SMAD3 signaling pathway (Chen et al., 2016). In terms of blue and brown modules, enrichment analysis revealed that genes in these modules were significantly associated with synthesis and organization of ECM and inflammatory response. Accumulation of ECM in atrium is the key feature of atrial fibrosis which is mainly mediated by TGF-β1 signaling pathway, renin-angiotensin-aldosterone axis, matrix metalloproteinases (MMPs), and tissue inhibitors of MMPs (Sohns and Marrouche, 2020). CFs are the pivotal cell type in the synthesis and organization of ECM, which account for up to 60% of cells in the cardiac muscle (Souders et al., 2009). In this study, the score of CFs calculated by ssGSEA was significantly increased in AF group, further validating the important role of CFs in AF. In addition to affecting atrial structure, atrial fibrosis could lead to slow conduction, increase heterogeneity of conduction, and turn the atrium into the substrate of AF initiation (Li et al., 1999). Atrial fibrosis is also related to the occurrence and perpetuation of focal and re-entry arrhythmia, which contributes to the maintenance of AF (Burstein and Nattel, 2008). Thus, future research in therapy of AF should take atrial fibrosis into consideration to improve clinical outcome (Sohns and Marrouche, 2020).
There were researches demonstrated that patients with AF exhibited systemic inflammation with the increasement of inflammatory markers like CRP, IL-6, and IL-8 in peripheral blood (Hu et al., 2015). Inflammatory cytokines contribute to the development of electrical remodeling directly. TNF could induce dysregulated Ca2+ in pulmonary vein cardiomyocytes (Lee et al., 2007) and inhibit the expression level of sarcoplasmic/endoplasmic reticulum Ca2+ ATPase 2a (SERCA2a) (Kao et al., 2010). Moreover, it was reported that TNF could decrease the expression levels of connexin-40 and connexin-43 and alter the intracellular distribution of connexin-43 in mice model (Sawaya et al., 2007; Liew et al., 2013). Inflammation also involves in the development and progression of structural remodeling. Inflammatory cytokines like CRP, heat shock protein (HSP) 27, and IL6 are positively associated with the atrial size (Hu et al., 2015). TNF could mediate structural remodeling via activating TGF-β signaling pathway and upregulating the expression level of MMP2 and MMP9 (Liew et al., 2013). In this study, the infiltration of immune cells and DEIRGs were analyzed to further investigate the role of inflammatory response in AF. The results indicated that the scores of macrophages, mast cells, neutrophils, TIL, and Treg were significantly increased in AF group. S. Kugler et al. found that the infiltration of macrophage in LA was significantly increased in patients with AF (Kugler et al., 2022). Mast cell frequently infiltrates in the fibrotic tissue and is an important cell type in mediating tissue fibrosis (Legere et al., 2019). Inhibition of mast cell was able to rescue cardiac fibrosis of AF mouse model (Legere et al., 2019). Zhang et al. reported that the infiltration of neutrophils in atrium was associated with the increased vulnerability to AF in a canine sterile pericarditis model (Zhang et al., 2011). An increased ratio of neutrophil to lymphocyte in peripheral blood was related to the increased incidence of new onset AF in patients after surgery. But the roles of TIL and Treg in AF are unclear and need further investigation. After that, DEIRGs between SR group and AF group were identified and the top five upregulated DEIRGs were CHGB, HLA-DRA, LYZ, IGKV1-17, and TYROBP. The roles of these IRGs in AF need further exploration, and these IRGs might be potential therapeutic targets for AF.
MiRNAs are noncoding RNAs with length of 21–24 nucleotide and mediate gene expression by binding to the 3′-untranslated region of mRNA which inhibits mRNA translation or inducing mRNA degradation (Blazejowska et al., 2021). CircRNAs, new type of endogenous non-coding RNAs can form covalent closed-loop structures through reverse splicing with longer half-life than mRNA (Guo et al., 2022). Both of miRNAs and circRNAs involve in the development and progression of cardiac diseases and might be potential biomarkers for cardiac diseases (Lu et al., 2022). The theory of ceRNA network suggests that non-coding RNAs like circRNAs contain site for binding to miRNAs and circRNAs could sequester and inactive miRNAs by acting as miRNA sponge, which further mediates the expression level of target mRNA of miRNAs (Guo et al., 2022). The construction of ceRNA network would benefit to identification of potential therapeutic targets or biomarkers for diseases. In this study, ceRNA networks of DEGs were constructed. circRNAs in these networks might be associated with the formation of substrate for AF and have potential predictive values for AF, which need further validation.
POAF is a common complication after cardiac surgery which is associated with increased risk of subsequent AF and mortality after cardiac surgery (Ahlsson et al., 2010). POAF could also occur after non-cardiac surgery, but we only discussed POAF after cardiac surgery in this study. Anderson et al. reported that reactive oxygen species (ROS) were increased in the atriums from patients who were susceptible to POAF (Anderson et al., 2014). The expression level of NF-κB was significantly upregulated in atrial tissues from patients with POAF (Li et al., 2009). Both ROS and inflammatory factors are key mediators for development and progression of AF (Hu et al., 2015). Transient factors during cardiac surgery like activation of autonomic nervous system, inflammation, and cardiopulmonary bypass significantly contribute to the initiation of POAF (Maesen et al., 2012). Besides, pre-existing atrial substrate like atrial fibrosis, inflammation, and electrical remodeling in patients also involves in the presence of POAF (Dobrev et al., 2019). In summary, the occurrence of POAF is largely due to the acute transient factors during cardiac surgery and the pre-existing atrial substrates which facilitates the initiation of AF (Dobrev et al., 2019). Therefore, we speculated that circRNAs in the ceRNA networks of DEGs might have potential predictive values for POAF, but these circRNAs need to be further screened and validated. To screen circRNAs with potential predictive values for POAF, circRNAs expression data of GSE97455 which included 15 plasma samples from patients without POAF and 15 plasma samples from patients with POAF were used and DEcircRNAs between patients with or without AF were identified. Then, the intersection between DEcircRNAs and circRNAs in the ceRNA networks of DEGs was extracted. After that, Upregulated circRNAs which is associated with upregulated DEGs and downregulated circRNAs which is associated with downregulated DEGs were selected as candidate biomarkers for POAF and 10 circRNAs were obtained. In order to identify hub circRNAs out of 10 circRNAs, predictive models for POAF were constructed using 10 circRNAs as features by SVM, RF, and neural network, both of which are machine learning algorithms and had been used in diagnosis of cardiovascular disease (Ahsan and Siddique, 2022). Top five circRNAs of each predictive models were identified according to the feature importance ranked by algorithms. Finally, the intersection between top five circRNAs of each predictive model was extracted. Hsa_circ_0006314 and hsa_circ_0055387 were identified as hub circRNAs and might be potential biomarkers for POAF. According to the ceRNA networks, hsa_circ_0055387 might be able to affect the expression of LBH and RELN via binding has-miR-27a-3p, and hsa_circ_0006314 might mediate the expression of RELN by interacting with has-miR-545-3p. But the exact roles of hsa_circ_0006314 and hsa_circ_0055387 in POAF and predictive values of hsa_circ_0006314 and hsa_circ_0055387 for POAF need further validation.
In conclusion, synthesis of ECM and inflammatory response were main processes involved in the development and progression of AF. Immune cells like macrophages, mast cells, neutrophils, TIL, and Treg were significantly infiltrated in LA tissues from patients with AF. The expression levels of IRGs such as CHGB, HLA-DRA, LYZ, IGKV1-17, and TYROBP were significantly upregulated in patients with AF, which were also correlated with infiltration of immune cells in patients. After screening, hsa_circ_0006314 and hsa_circ_0055387 were found to have potential predictive values for POAF.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: GSE79768: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE79768 GSE115574: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE115574 GSE97455: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE97455.
Author contributions
YC analyzed the data, completed figures, and wrote manuscript. TO, YY, and CF processed the raw data and predicted circRNA-miRNA and miRNA-mRNA interaction pairs. C-ET, JL, and FL designed the research. JL and FL reviewed and edited the manuscript.
Funding
This study was funded by the National Natural Science Foundation of China (81873494 and 82070352) and the Hunan Key R&D Programs (2021SK2022).
Acknowledgments
We thank Yangjie Zhou from Xiangya Hospital, Central South University for advices on writing this manuscript.
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
Ahlsson, A., Fengsrud, E., Bodin, L., and Englund, A. (2010). Postoperative atrial fibrillation in patients undergoing aortocoronary bypass surgery carries an eightfold risk of future atrial fibrillation and a doubled cardiovascular mortality. Eur. J. Cardiothorac. Surg. 37 (6), 1353–1359. doi:10.1016/j.ejcts.2009.12.033
Ahsan, M. M., and Siddique, Z. (2022). Machine learning-based heart disease diagnosis: A systematic literature review. Artif. Intell. Med. 128, 102289. doi:10.1016/j.artmed.2022.102289
Anderson, E. J., Efird, J. T., Davies, S. W., O'Neal, W. T., Darden, T. M., Thayne, K. A., et al. (2014). Monoamine oxidase is a major determinant of redox balance in human atrial myocardium and is associated with postoperative atrial fibrillation. J. Am. Heart Assoc. 3 (1), e000713. doi:10.1161/JAHA.113.000713
Baum, J., and Duffy, H. S. (2011). Fibroblasts and myofibroblasts: What are we talking about? J. Cardiovasc. Pharmacol. 57 (4), 376–379. doi:10.1097/FJC.0b013e3182116e39
Benjamin, E. J., Wolf, P. A., D'Agostino, R. B., Silbershatz, H., Kannel, W. B., and Levy, D. (1998). Impact of atrial fibrillation on the risk of death: The framingham heart study. Circulation 98 (10), 946–952. doi:10.1161/01.cir.98.10.946
Blazejowska, E., Urbanowicz, T., Gasecka, A., Olasinska-Wisniewska, A., Jaguszewski, M. J., Targonski, R., et al. (2021). Diagnostic and prognostic value of miRNAs after coronary artery bypass grafting: A review. Biol. (Basel) 10 (12), 1350. doi:10.3390/biology10121350
Burstein, B., and Nattel, S. (2008). Atrial fibrosis: Mechanisms and clinical relevance in atrial fibrillation. J. Am. Coll. Cardiol. 51 (8), 802–809. doi:10.1016/j.jacc.2007.09.064
Chen, X., Xu, J., Jiang, B., and Liu, D. (2016). Bone morphogenetic protein-7 antagonizes myocardial fibrosis induced by atrial fibrillation by restraining transforming growth factor-β (TGF-β)/Smads signaling. Med. Sci. Monit. 22, 3457–3468. doi:10.12659/msm.897560
Christ, T., Boknik, P., Wohrl, S., Wettwer, E., Graf, E. M., Bosch, R. F., et al. (2004). L-type Ca2+ current downregulation in chronic human atrial fibrillation is associated with increased activity of protein phosphatases. Circulation 110 (17), 2651–2657. doi:10.1161/01.CIR.0000145659.80212.6A
Colilla, S., Crow, A., Petkun, W., Singer, D. E., Simon, T., and Liu, X. (2013). Estimates of current and future incidence and prevalence of atrial fibrillation in the U.S. adult population. Am. J. Cardiol. 112 (8), 1142–1147. doi:10.1016/j.amjcard.2013.05.063
Dobrev, D., Aguilar, M., Heijman, J., Guichard, J. B., and Nattel, S. (2019). Postoperative atrial fibrillation: Mechanisms, manifestations and management. Nat. Rev. Cardiol. 16 (7), 417–436. doi:10.1038/s41569-019-0166-5
Dzeshka, M. S., Lip, G. Y., Snezhitskiy, V., and Shantsila, E. (2015). Cardiac fibrosis in patients with atrial fibrillation: Mechanisms and clinical implications. J. Am. Coll. Cardiol. 66 (8), 943–959. doi:10.1016/j.jacc.2015.06.1313
Faria, A., and Persaud, S. J. (2017). Cardiac oxidative stress in diabetes: Mechanisms and therapeutic potential. Pharmacol. Ther. 172, 50–62. doi:10.1016/j.pharmthera.2016.11.013
Greenberg, J. W., Lancaster, T. S., Schuessler, R. B., and Melby, S. J. (2017). Postoperative atrial fibrillation following cardiac surgery: A persistent complication. Eur. J. Cardiothorac. Surg. 52 (4), 665–672. doi:10.1093/ejcts/ezx039
Guo, L., Jia, L., Luo, L., Xu, X., Xiang, Y., Ren, Y., et al. (2022). Critical roles of circular RNA in tumor metastasis via acting as a sponge of miRNA/isomiR. Int. J. Mol. Sci. 23 (13), 7024. doi:10.3390/ijms23137024
Haissaguerre, M., Wright, M., Hocini, M., and Jais, P. (2008). The substrate maintaining persistent atrial fibrillation. Circ. Arrhythm. Electrophysiol. 1 (1), 2–5. doi:10.1161/CIRCEP.108.764233
Heijman, J., Guichard, J. B., Dobrev, D., and Nattel, S. (2018). Translational challenges in atrial fibrillation. Circ. Res. 122 (5), 752–773. doi:10.1161/CIRCRESAHA.117.311081
Hindricks, G., Potpara, T., Dagres, N., Arbelo, E., Bax, J. J., Blomstrom-Lundqvist, C., et al. (2020). ESC Guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European Association for Cardio-Thoracic Surgery (EACTS): The Task Force for the diagnosis and management of atrial fibrillation of the European Society of Cardiology (ESC) Developed with the special contribution of the European Heart Rhythm Association (EHRA) of the ESC. Eur. Heart J. 42 (5), 373–498.
Hu, Y. F., Chen, Y. J., Lin, Y. J., and Chen, S. A. (2015). Inflammation and the pathogenesis of atrial fibrillation. Nat. Rev. Cardiol. 12 (4), 230–243. doi:10.1038/nrcardio.2015.2
Kalifa, J., Tanaka, K., Zaitsev, A. V., Warren, M., Vaidyanathan, R., Auerbach, D., et al. (2006). Mechanisms of wave fractionation at boundaries of high-frequency excitation in the posterior left atrium of the isolated sheep heart during atrial fibrillation. Circulation 113 (5), 626–633. doi:10.1161/CIRCULATIONAHA.105.575340
Kao, Y. H., Chen, Y. C., Cheng, C. C., Lee, T. I., Chen, Y. J., and Chen, S. A. (2010). Tumor necrosis factor-alpha decreases sarcoplasmic reticulum Ca2+-ATPase expressions via the promoter methylation in cardiomyocytes. Crit. Care Med. 38 (1), 217–222. doi:10.1097/CCM.0b013e3181b4a854
Kraj, P. (2022). Bone morphogenetic proteins shape Treg cells. Front. Immunol. 13, 865546. doi:10.3389/fimmu.2022.865546
Kugler, S., Onodi, Z., Ruppert, M., Sayour, A. A., Olah, A., Benke, K., et al. (2022). Inflammasome activation in end-stage heart failure-associated atrial fibrillation. Esc. Heart Fail. 9 (4), 2747–2752. doi:10.1002/ehf2.13972
Langfelder, P., and Horvath, S. (2008). Wgcna: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Lee, S. H., Chen, Y. C., Chen, Y. J., Chang, S. L., Tai, C. T., Wongcharoen, W., et al. (2007). Tumor necrosis factor-alpha alters calcium handling and increases arrhythmogenesis of pulmonary vein cardiomyocytes. Life Sci. 80 (19), 1806–1815. doi:10.1016/j.lfs.2007.02.029
Legere, S. A., Haidl, I. D., Legare, J. F., and Marshall, J. S. (2019). Mast cells in cardiac fibrosis: New insights suggest opportunities for intervention. Front. Immunol. 10, 580. doi:10.3389/fimmu.2019.00580
Li, D., Fareh, S., Leung, T. K., and Nattel, S. (1999). Promotion of atrial fibrillation by heart failure in dogs: Atrial remodeling of a different sort. Circulation 100 (1), 87–95. doi:10.1161/01.cir.100.1.87
Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 42, D92–D97. doi:10.1093/nar/gkt1248
Li, J. Y., Lai, Y. J., Yeh, H. I., Chen, C. L., Sun, S., Wu, S. J., et al. (2009). Atrial gap junctions, NF-kappaB and fibrosis in patients undergoing coronary artery bypass surgery: The relationship with postoperative atrial fibrillation. Cardiology 112 (2), 81–88. doi:10.1159/000141012
Li, N., and Brundel, B. (2020). Inflammasomes and proteostasis novel molecular mechanisms associated with atrial fibrillation. Circ. Res. 127 (1), 73–90. doi:10.1161/CIRCRESAHA.119.316364
Liew, R., Khairunnisa, K., Gu, Y., Tee, N., Yin, N. O., Naylynn, T. M., et al. (2013). Role of tumor necrosis factor-alpha in the pathogenesis of atrial fibrosis and development of an arrhythmogenic substrate. Circ. J. 77 (5), 1171–1179. doi:10.1253/circj.cj-12-1155
Lin, Y. J., Tai, C. T., Kao, T., Chang, S. L., Lo, L. W., Tuan, T. C., et al. (2009). Spatiotemporal organization of the left atrial substrate after circumferential pulmonary vein isolation of atrial fibrillation. Circ. Arrhythm. Electrophysiol. 2 (3), 233–241. doi:10.1161/CIRCEP.108.812024
Lu, P., Ding, F., Xiang, Y. K., Hao, L., and Zhao, M. (2022). Noncoding RNAs in cardiac hypertrophy and heart failure. Cells 11 (5), 777. doi:10.3390/cells11050777
Maesen, B., Nijs, J., Maessen, J., Allessie, M., and Schotten, U. (2012). Post-operative atrial fibrillation: A maze of mechanisms. Europace 14 (2), 159–174. doi:10.1093/europace/eur208
Michaud, G. F., and Stevenson, W. G. (2021). Atrial fibrillation. N. Engl. J. Med. 384 (4), 353–361. doi:10.1056/NEJMcp2023658
Oral, H., Chugh, A., Good, E., Crawford, T., Sarrazin, J. F., Kuhne, M., et al. (2008). Randomized evaluation of right atrial ablation after left atrial ablation of complex fractionated atrial electrograms for long-lasting persistent atrial fibrillation. Circ. Arrhythm. Electrophysiol. 1 (1), 6–13. doi:10.1161/CIRCEP.107.748780
Palatinus, J. A., and Das, S. (2015). Your father and grandfather's atrial fibrillation: A review of the genetics of the most common pathologic cardiac dysrhythmia. Curr. Genomics 16 (2), 75–81. doi:10.2174/1389202916666150108222031
Sahadevan, J., Ryu, K., Peltz, L., Khrestian, C. M., Stewart, R. W., Markowitz, A. H., et al. (2004). Epicardial mapping of chronic atrial fibrillation in patients: Preliminary observations. Circulation 110 (21), 3293–3299. doi:10.1161/01.CIR.0000147781.02738.13
Sawaya, S. E., Rajawat, Y. S., Rami, T. G., Szalai, G., Price, R. L., Sivasubramanian, N., et al. (2007). Downregulation of connexin40 and increased prevalence of atrial arrhythmias in transgenic mice with cardiac-restricted overexpression of tumor necrosis factor. Am. J. Physiol. Heart Circ. Physiol. 292 (3), H1561–H1567. doi:10.1152/ajpheart.00285.2006
Sohns, C., and Marrouche, N. F. (2020). Atrial fibrillation and cardiac fibrosis. Eur. Heart J. 41 (10), 1123–1131. doi:10.1093/eurheartj/ehz786
Souders, C. A., Bowers, S. L. K., and Baudino, T. A. (2009). Cardiac fibroblast: The renaissance cell. Circ. Res. 105, 1164–1176. doi:10.1161/CIRCRESAHA.109.209809
Staerk, L., Sherer, J. A., Ko, D., Benjamin, E. J., and Helm, R. H. (2017). Atrial fibrillation: Epidemiology, pathophysiology, and clinical outcomes. Circ. Res. 120, 1501–1517. doi:10.1161/CIRCRESAHA.117.309732
Tsai, F. C., Lin, Y. C., Chang, S. H., Chang, G. J., Hsu, Y. J., Lin, Y. M., et al. (2016). Differential left-to-right atria gene expression ratio in human sinus rhythm and atrial fibrillation: Implications for arrhythmogenesis and thrombogenesis. Int. J. Cardiol. 222, 104–112. doi:10.1016/j.ijcard.2016.07.103
Tucker, N. R., Clauss, S., and Ellinor, P. T. (2016). Common variation in atrial fibrillation: Navigating the path from genetic association to mechanism. Cardiovasc. Res. 109 (4), 493–501. doi:10.1093/cvr/cvv283
Ueha, S., Shand, F. H., and Matsushima, K. (2012). Cellular and molecular mechanisms of chronic inflammation-associated organ fibrosis. Front. Immunol. 3, 71. doi:10.3389/fimmu.2012.00071
Yamashita, T., Sekiguchi, A., Iwasaki, Y. K., Date, T., Sagara, K., Tanabe, H., et al. (2010). Recruitment of immune cells across atrial endocardium in human atrial fibrillation. Circ. J. 74 (2), 262–270. doi:10.1253/circj.cj-09-0644
Yan, J., Thomson, J. K., Zhao, W., Wu, X., Gao, X., DeMarco, D., et al. (2018). The stress kinase JNK regulates gap junction Cx43 gene expression and promotes atrial fibrillation in the aged heart. J. Mol. Cell. Cardiol. 114, 105–115. doi:10.1016/j.yjmcc.2017.11.006
Zhang, J., Johnsen, S. P., Guo, Y., and Lip, G. Y. H. (2021). Epidemiology of atrial fibrillation: Geographic/ecological risk factors, age, sex, genetics. Card. Electrophysiol. Clin. 13 (1), 1–23. doi:10.1016/j.ccep.2020.10.010
Zhang, Z., Zhang, C., Wang, H., Zhao, J., Liu, L., Lee, J., et al. (2011). n-3 polyunsaturated fatty acids prevents atrial fibrillation by inhibiting inflammation in a canine sterile pericarditis model. Int. J. Cardiol. 153 (1), 14–20. doi:10.1016/j.ijcard.2010.08.024
Ziolo, M. T., and Mohler, P. J. (2015). Defining the role of oxidative stress in atrial fibrillation and diabetes. J. Cardiovasc. Electrophysiol. 26 (2), 223–225. doi:10.1111/jce.12560
Keywords: atrial fibrillation, postoperative atrial fibrillation, immune cells, inflammation, circRNA, biomarker
Citation: Chen Y, Ouyang T, Yin Y, Fang C, Tang C-E, Luo J and Luo F (2022) Analysis of infiltrated immune cells in left atriums from patients with atrial fibrillation and identification of circRNA biomarkers for postoperative atrial fibrillation. Front. Genet. 13:1003366. doi: 10.3389/fgene.2022.1003366
Received: 26 July 2022; Accepted: 22 November 2022;
Published: 09 December 2022.
Edited by:
Laura Iop, University of Padua, ItalyReviewed by:
Changfa Guo, Fudan University, ChinaSaeid Hosseini, Iran University of Medical Sciences, Iran
Copyright © 2022 Chen, Ouyang, Yin, Fang, Tang, Luo and Luo. 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: Jingmin Luo, jmluo99@126.com; Fanyan Luo, drlfy@csu.edu.cn