Guanzhong Chen1,2†
Liwei Liu1,3†
Huanqiang Li1†Zhubin Lun1,4†Ziling Mai1,5Wenguang Lai1,5Enzhao Chen1Chunyun Zhou1Sijia Yu1,3Junqing Yang1Shiqun Chen1*Jiyan Chen1,2,3*
Yong Liu1,2,3*- 1Guangdong Provincial Key Laboratory of Coronary Heart Disease Prevention, Department of Cardiology, Guangdong Provincial People’s Hospital, Guangdong Cardiovascular Institute, Guangdong Academy of Medical Sciences, Guangzhou, China
- 2Guangdong Provincial People’s Hospital, School of Medicine, South China University of Technology, Guangzhou, China
- 3The Second School of Clinical Medicine, Southern Medical University, Guangzhou, China
- 4The First School of Clinical Medicine, Guangdong Medical University, Zhanjiang, China
- 5Guangdong Provincial People’s Hospital, School of Biology and Biological Engineering, South China University of Technology, Guangzhou, China
Background: Acute myocardial infarction (AMI), characterized by an event of myocardial necrosis, is a common cardiac emergency worldwide. However, the genetic mechanisms of AMI remain largely elusive.
Methods: A genome-wide association study dataset of AMI was obtained from the CARDIoGRAMplusC4D project. A transcriptome-wide association study (TWAS) was conducted using the FUSION tool with gene expression references of the left ventricle and whole blood. Significant genes detected by TWAS were subjected to Gene Ontology (GO) enrichment analysis. Then the TWAS results of AMI were integrated with mRNA expression profiling to identify common genes and biological processes. Finally, the identified common genes were validated by RT-qPCR analysis.
Results: TWAS identified 1,050 genes for the left ventricle and 1,079 genes for whole blood. Upon comparison with the mRNA expression profile, 4 common genes were detected, including HP (PTWAS = 1.22 × 10–3, PGEO = 4.98 × 10–2); CAMP (PTWAS = 2.48 × 10–2, PGEO = 2.36 × 10–5); TNFAIP6 (PTWAS = 1.90 × 10–2, PGEO = 3.46 × 10–2); and ARG1 (PTWAS = 8.35 × 10–3, PGEO = 4.93 × 10–2). Functional enrichment analysis of the genes identified by TWAS detected multiple AMI-associated biological processes, including autophagy of mitochondrion (GO: 0000422) and mitochondrion disassembly (GO: 0061726).
Conclusion: This integrative study of TWAS and mRNA expression profiling identified multiple candidate genes and biological processes for AMI. Our results may provide a fundamental clue for understanding the genetic mechanisms of AMI.
Introduction
Acute myocardial infarction (AMI), defined by an event of myocardial necrosis caused by the rupture of an atherosclerotic plaque, has become the leading cause of death and disability worldwide (Murray et al., 2015). Previous studies have identified several factors associated with the onset of AMI, including age, gender, hypertension, diabetes mellitus, smoking, alcohol consumption and physical exertion (Mittleman et al., 1993, 1995, 1999; Lipovetzky et al., 2007). Moreover, increasing evidence has demonstrated that genetic factors also contribute to the development of AMI. Genome-wide association studies (GWAS) have identified several genetic loci that confer susceptibility to AMI, such as CDKN2A/2B, CELSR2-PSRC1-SORT1, CXCL12, ABO, LDLR, and APOA5 (Qi et al., 2011; Do et al., 2015; Nikpay et al., 2015).
Although GWAS have successfully identified great numbers of loci, only a small proportion (approximately 10.6%) of phenotypic variance explains coronary artery disease or AMI (Deloukas et al., 2013). This is probably because significant single nucleotide polymorphisms (SNPs) identified by GWAS are mainly enriched in non-coding regulatory regions. Due to limited genomic annotation, it is difficult to explain the relative risk of diseases. Transcriptome-wide association study (TWAS) has identified significant gene expression-trait associations by integrating gene expression data with GWAS data (Gusev et al., 2016). Multiple studies have indicated that TWAS exhibited good performance in identifying novel genes of complex cardiovascular diseases, such as atrial fibrillation (Zhang et al., 2019) and calcific aortic valve stenosis (Thériault et al., 2018).
In this study, we conducted an integrative analysis of TWAS and mRNA expression profiles of AMI. TWAS was first performed on a large-scale GWAS from the CARDIoGRAMplusC4D Consortium (Nikpay et al., 2015) to identify novel associated genes. The AMI-associated genes identified by TWAS were further subjected to Gene Ontology (GO) analysis to explore the core biological processes or signaling pathways. TWAS results were finally compared to those of mRNA expression profiling to identify common genes.
Materials and Methods
GWAS Summary Datasets of AMI
We collected a large-scale GWAS summary data of coronary artery diseases from CARDIoGRAMplusC4D. Briefly, the study is a meta-analysis of 48 coronary artery disease studies. We chose one of the subgroup results (GWAS summary data for AMI, AMI patients, n = 43,676; controls, n = 1,28,199) to perform further TWAS analysis (Nikpay et al., 2015; Luo et al., 2019). SNP genotyping was conducted using a combination of array platforms, including Affymetrix GeneChip Human Mapping 500K, Affymetrix Genome-Wide Human SNP Array 6.0 and Illumina HumanHap300 BeadChip. In this dataset, an additive model where the log (genotype risk ratio) for a genotype was proportional to the number of risk alleles was used. For detailed information on the cohorts, genotyping, and quality control approaches, please refer to the published studies.
mRNA Expression Profiles of AMI
The mRNA expression profiles of AMI were obtained from two series of GEO datasets (GEO accession ID: GSE29532 and GSE61144). These two datasets were generated using Affymetrix Human Exon 1.0 ST Array (GPL5175) and Sentrix Human-6 v2 Expression BeadChip (GPL6106), respectively. For detailed information, please refer to original publications (Silbiger et al., 2013; Park et al., 2015). We used the online tool GEO2R to conduct differential gene expression analysis comparing the expression of whole blood mRNA between AMI cases and controls. Genes with fold change (FC) ≥2 or ≤0.5 and P < 0.05 were defined as differentially expressed genes (DEGs).
TWAS Analysis
TWAS analysis of AMI was performed using FUSION software. Briefly, FUSION used a relatively small set of individuals with both genotype and gene expression data to compute gene expression weights for different tissues. TWAS analysis uses pre-computed gene expression weights together with GWAS summary statistics to calculate the association of every gene to a given disease. The genetic values of expression are computed one probe set at a time using SNP genotyping data located 500 kb on either side of the gene boundary. In this study, gene expression weights of whole blood and left ventricle were derived from the FUSION website1, which is attached to the GTEx database. Genes with significant and suggestive association signals were identified by P < 0.05.
Functional Characterization of GO and KEGG Analysis
Associated genes identified by TWAS and mRNA expression profiles were subjected to GO and KEGG analysis. Analyses were evaluated using the Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics tool2. A GO term or KEGG pathway with adjusted P-value < 0.05 was considered statistically significant (Huang da et al., 2009). Visualization of GO analysis was conducted with R package “ggplot2.”
Animals
Adult male Sprague-Dawley (SD) rats weighing 200–220 g were purchased from the Experimental Animal Center at Guangzhou University of Chinese Medicine (Guangzhou, China). The animal facilities and protocols were approved by the Guangdong provincial people’s hospital Ethics Committee. Before the beginning of the experiments, all animals were allowed to acclimate to the facility in a standard SPF-class animal room for 1 week. Six rats were randomly divided into sham and cardiac ischemia and reperfusion groups, each group consisting of three rats. Before experiments, rats were fasted overnight and deprived of water.
Acute Myocardial Infarction Rat Model Establishment
Rats were weighed and anesthetized with 10% chloral hydrate (0.4 ml/100 g) by intraperitoneal injection. Tracheostomy was performed and rats were ventilated with room air from a positive pressure rodent ventilator. Lead II electrocardiography was recorded throughout the experiment. The left anterior descending coronary artery was then identified and ligated at approximately 2 mm distal to its origin after a left-side thoracotomy at the fourth intercostal space and heart exposure. A small vinyl tube was used to occlude the LAD by pulling the thread. A ligature for 30 min was then followed by a 24 h reperfusion in the AMI group. The same experimental performance was conduct in the sham group except occluding the LAD. ST-elevation from electrocardiography was used to confirm the establishment of the AMI rat model. At the end of the study, the animals were sacrificed and the hearts were collected for subsequent experiments.
Hematoxylin-Eosin (HE) Staining
Rat heart tissues were fixed in 10% neutral buffered formalin for more than 24 h and then embedded in paraffin. Tissue sections were cut as 4 mm-thick using a microtome. Histopathological evaluation of the tissue sections was conducted by hematoxylin-eosin staining. Finally, tissue morphology was observed.
Reverse Transcription Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) Analysis
RT-qPCR was used to verify the expression of the identified genes in heart tissues of both AMI and sham rats. Total RNA was extracted from the collected heart from the AMI rats and sham rats by Trizol Universal reagent (TIANGEN, DP424). Reverse transcription was then carried out using total RNA. cDNA was synthesized via reverse transcription using a PrimeScriptTM RT Reagent Kit (RR037) according to the manufacturer’s protocols (Takara, Japan). Then, quantitative real-time PCR was carried out using TB Green® Premix Ex TaqTM II (Tli RNaseH Plus, RR820A) according to the manufacturer’s protocols (Takara, Japan). The primers used to amplify genes in the reactions were synthesized by Generay (Shanghai, China). All RT-qPCR primer sequences are shown in Supplementary Table S4. qPCR was performed in a CFX Connect Real-Time System (Bio-Rad, United States). The relative expression levels of the miRNAs were normalized to those of the reference gene GAPDH and calculated using the 2–▲▲CT method. All reactions were repeated in triplicate.
Results
TWAS Analysis Results
A total of 1,862 genes were identified by TWAS across the two tissues (Supplementary Tables S1, S2). For the left ventricle, TWAS identified 1,050 genes with TWAS P < 0.05, such as SURF2 (PTWAS = 1.25 × 10–12), RP11-378J18.8 (PTWAS = 8.33 × 10–11), ABO (PTWAS = 1.70 × 10–7), ICA1L (PTWAS = 1.78 × 10–7), and MIA3 (PTWAS = 2.33 × 10–7). For whole blood, TWAS identified 1,079 genes with TWAS P < 0.05, such as PSRC1 (PTWAS = 1.18 × 10–10), CARF (PTWAS = 2.41 × 10–10), RP11-378J18.8 (PTWAS = 4.83 × 10–9), GGCX (PTWAS = 1.01 × 10–8), and POC1B (PTWAS = 1.07 × 10–8). Among all these genes, 266 were detected in both the left ventricle and whole blood. The top 20 AMI associated genes identified by TWAS are shown in Table 1.
Table 1. List of top 20 candidate genes identified by TWAS for AMI.
Integrative Analysis of TWAS and mRNA Expression Profiles
After further comparison of differentially expressed genes identified by TWAS and mRNA expression profiles, we identified four overlapping genes, including HP (PTWAS = 1.22 × 10–3, PGEO = 4.98 × 10–2); CAMP (PTWAS = 2.48 × 10–2, PGEO = 2.36 × 10–5); TNFAIP6 (PTWAS = 1.90 × 10–2, PGEO = 3.46 × 10–2); and ARG1 (PTWAS = 8.35 × 10–3, PGEO = 4.93 × 10–2). Detailed information on the common genes is shown in Table 2.
Table 2. List of overlapping candidate genes identified by TWAS and GEO mRNA profiling for AMI.
Functional Analyses
GO enrichment analysis of genes identified by TWAS detected 10 significant biological process terms, such as organelle disassembly (GO: 1903008), modification of synaptic structure (GO: 0099563), autophagy of mitochondrion (GO: 0000422) and mitochondrion disassembly (GO: 0061726). The associated GO terms identified by TWAS are listed in Figure 1 and Table 3. For KEGG pathway enrichment analysis (Supplementary Table S3), 5 significant pathways were detected, including bacterial invasion of epithelial cells (hsa05100, adjust P-value = 1.30 × 10–2), insulin resistance (hsa04931, adjust P-value = 1.30 × 10–2), proteoglycans in cancer (hsa05205, adjust P-value = 1.73 × 10–2), focal adhesion (hsa04510, adjust P-value = 4.80 × 10–2) and regulation of actin cytoskeleton (hsa04810, adjust P-value = 4.87 × 10–2).
Figure 1. GO enrichment analysis of the genes identified by TWAS.
Table 3. Gene ontology enrichment analysis results of AMI associated genes identified by TWAS.
qPCR Validation
To confirm myocardial infarction, HE staining was performed to detect the histopathological changes between the AMI group and sham group. Comparing with the normal cardiac histology in the sham group, the AMI group showed significant inflammatory cells infiltration, cardiomyocyte necrosis and fiber disruption (Figure 2A). RT-qPCR was used to confirm the expression of the identified genes in sham and AMI rat hearts. As shown in Figure 2B, the HP was significantly down-regulated, while TNFAIP6 and ARG1 were significantly up-regulated. However, the expression of CAMP was not significant. In general, RT-qPCR results confirmed the integrative analysis result of TWAS and mRNA profiles.
Figure 2. Validation of identified gene expression. (A) Pathological changes in rat myocardium tissues (HE staining × 200). (B) RT-qPCR analysis of four identified genes. Rat GAPDH was used as internal controls. The measurement data are expressed as the means ± SD. *P < 0.05.
Discussion
To demonstrate the implication and potential effects of genetic factors in the development of AMI, we performed TWAS by using large-scale GWAS data. Then, we compared the TWAS results with the mRNA expression profiles to explore common genes, in which 4 genes overlapped, including HP, CAMP, TNFAIP6, and ARG1. GO enrichment suggested that 10 biological processes, including autophagy of mitochondrion (GO: 0000422) and mitochondrion disassembly (GO: 0061726) were associated with AMI.
In our study, one of the notable genes identified was CAMP. CAMP encodes a member of the antimicrobial peptide family, characterized by a highly conserved N-terminal signal peptide containing a cathelin domain and a structurally variable cationic antimicrobial peptide (Zasloff, 2002; Hilchie et al., 2013). In addition to its antibacterial activities, the CAMP protein functions in cellular chemotaxis and inflammatory response regulation (Lande et al., 2007). Previous studies indicated that CAMP, acting as an activator of Akt, ERK and FoxO3a, protects against cardiomyocyte apoptosis in AMI (Bei et al., 2019). Moreover, Bei et al. (2019) showed that serum levels of LL-37 (human analog of CAMP) were significantly reduced in AMI patients. Low-level LL-37 might also be associated with a worse prognosis in AMI. Furthermore, other studies showed that CAMP was related to dyslipidemia (Chamorro et al., 2009; Dashty et al., 2014), indicating that CAMP might induce lipid abnormalities to increase the risk of AMI.
Another susceptibility gene was TNFIP6, also called TSG-6. TNFIP6 encodes a secretory protein of the hyaluronan-binding protein family, which is important in the protease network associated with inflammation (Milner et al., 2006; Heng et al., 2008). Accumulating evidence supports that inflammation plays an important role in the development of AMI (Wisniewski and Vilcek, 2004; Milner et al., 2006). A previous study showed that human multipotent stromal cells (hMSCs) enhance the repair of myocardial tissue during AMI. After knockdown of TSG-6 in hMSCs, positive improvements in inflammatory responses, infarct size, and cardiac function induced by hMSC injection were significantly attenuated (Lee et al., 2009). Therefore, TNFIP6 may represent an important gene involved in the development of AMI by suppressing the inflammatory response.
ARG1, also called arginase 1, was another gene identified by TWAS analysis. It was first studied in the urea cycle (Iyer et al., 1998), which can catalyze the hydrolysis of arginine and then convert arginine to ornithine and urea (Sin et al., 2015). However, ARG1 also could induce endothelial dysfunction and promote the progress of atherosclerosis by competing with nitric oxide synthase (NOS) for the common substrate-L-arginine and inhibited of biosynthesis of nitric oxide (NO) (Santhanam et al., 2007). It is these kinds of functions making ARG1 possible to participate in the development of AMI. Several studies have revealed this relationship between ARG1 and AMI. Zhu et al. (2015) found out that ARG1 can downregulate NO biosynthesis and promote the enrichment of inflammatory cells into infarcted heart areas in mice. In the latest research (Zhang et al., 2020), ARG1 was determined as a potential biomarker for AMI with the area under the receiver operating characteristic curve (AUC) up to 0.776. All together the application value of ARG1 in AMI is worth expecting.
The haptoglobin (Hp) gene locus at chromosome 16q22 is polymorphic with two alleles (Bowman and Kurosky, 1982). Interestingly, it is reported that Haptoglobin phenotypes differ in their ability in different diseases (Melamed-Frank et al., 2001; Bamm et al., 2004; Asleh et al., 2005). This gene has been revealed to play an important role in AMI. Hp 1-1 phenotype was determined to have an association with increased risk of AMI incidence in Chinese patients with T2D in the latest research (Gurung et al., 2019). Suleiman et al. (2005) have reported that DM individuals with the Hp 2 allele have significantly larger myocardial infarctions than DM individuals homozygous for the 1 allele, which was confirmed in mice by Blum et al. (2007). The underlying mechanism of this relationship may be due to increased oxidative stress in DM Hp 2 mice (Blum et al., 2007). Therefore, Hp is also an important susceptibility gene in the development and prognosis of AMI.
Functional enrichment analysis identified mitochondrial function related GO terms (GO: 0000422 and GO: 0061726) were enriched. It is well known that mitochondria are a source of energy for cardiomyocyte cells since they produce most of the adenosine triphosphate. Mitochondrial dysfunction at the onset of acute myocardial ischemia is a critical determinant of cell death in AMI patients (Hernandez-Resendiz et al., 2020). Also, efficient mitophagy response may allow cardiomyocytes to alleviate the hypoxia or nutritional stress (Lopez-Crisosto et al., 2017). Kubli et al. (2013) found that Parkin-deficient mice had significantly reduced mitophagy and accumulated dysfunctional mitochondria after myocardial infarction, leading to exacerbated cardiac injury and reduced survival. Along similar lines, pre- or post-conditioning could also mediate cardio-protective effects by activating the mitophagy during AMI (Huang et al., 2011; Wu et al., 2013; Yu et al., 2015). Additionally, a previous study showed that inflammatory conditions or arginine catabolism could affect mitochondrial function (Lee and Huttemann, 2014; Tang et al., 2020). Sun et al. indicated that the anti-inflammatory effects of CAMP were dependent its ability to stimulate mitochondrial biogenesis and maintain mitochondrial homeostasis (Sun et al., 2014). Tang et al. (2020) also found that the loss of ARG1 could cause arginine accumulation and mitochondrial dynamics disruption. Taken all together, the common genes shared by TWAS analysis and mRNA profiles may play key roles in the development of AMI by mediating the mitochondrial dysfunction.
To the best of our knowledge, this is the first large-scale comprehensive study in AMI combining TWAS and mRNA gene expression profiles. By integrating the GWAS data with gene expression profiles, common genes and potential pathways might be more convincing. However, some limitations of this study need to be noted. First, although TWAS was conducted to identify several casual genes associated with AMI, other loci that are independent of cis-gene expression effects may have been missed. Second, the novel candidate genes and pathways related to AMI resulted from bioinformatics analysis. Further functional studies are needed to confirm our findings and to clarify the potential biological processes underlying AMI.
Conclusion
In conclusion, we integrated a tissue-related TWAS and gene expression profile to identify 4 common genes and biological pathways related to AMI. These results may provide novel insights into investigating the pathogenesis of AMI and serve as a fundamental clue for understanding the genetic mechanisms of AMI.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by the Guangdong provincial people’s hospital Ethics Committee.
Author Contributions
GC, LL, YL, and JC: conception and design. HL, ZL, and SC: administrative support. GC, HL, and ZM: collection and assembly of data. GC and SC: data analysis and interpretation. All authors: wrote and approved the final manuscript.
Funding
This study was supported by the National Science Foundation of China (Grant Nos. 81670339 and 81970311), Cardiovascular Research Foundation Project of the Chinese Medical Doctor Association (SCRFCMDA201216), the Progress of Science and Technology Project in Guangdong Province (2017B020247060), Beijing Lisheng Cardiovascular Health Foundation (LHJJ20141751 and LHJJ201612127), and Dengfeng Project in Guangdong Province (DFJH201919 and DFJH2020026). The funding body plays no direct role in the design of the study, and collection, analysis, and interpretation of data, and in writing the 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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.616492/full#supplementary-material
Footnotes
References
Asleh, R., Guetta, J., Kalet-Litman, S., Miller-Lotan, R., and Levy, A. P. (2005). Haptoglobin genotype and diabetes-dependent differences in iron-mediated oxidative stress in vitro and in vivo. Circ. Res. 96, 435–441. doi: 10.1161/01.RES.0000156653.05853.b9
Bamm, V. V., Tsemakhovich, V. A., Shaklai, M., and Shaklai, N. (2004). Haptoglobin phenotypes differ in their ability to inhibit heme transfer from hemoglobin to LDL. Biochemistry 43, 3899–3906. doi: 10.1021/bi0362626
Bei, Y., Pan, L. L., Zhou, Q., Zhao, C., Xie, Y., Wu, C., et al. (2019). Cathelicidin-related antimicrobial peptide protects against myocardial ischemia/reperfusion injury. BMC Med. 17:42. doi: 10.1186/s12916-019-1268-y
Blum, S., Asaf, R., Guetta, J., Miller-Lotan, R., Asleh, R., Kremer, R., et al. (2007). Haptoglobin genotype determines myocardial infarct size in diabetic mice. J. Am. Coll. Cardiol. 49, 82–87. doi: 10.1016/j.jacc.2006.08.044
Bowman, B. H., and Kurosky, A. (1982). Haptoglobin: the evolutionary product of duplication, unequal crossing over, and point mutation. Adv. Hum. Genet. 12, 189-261, 453–454. doi: 10.1007/978-1-4615-8315-8_3
Chamorro, C. I., Weber, G., Grönberg, A., Pivarcsi, A., and Ståhle, M. (2009). The human antimicrobial peptide LL-37 suppresses apoptosis in keratinocytes. J. Invest. Dermatol. 129, 937–944. doi: 10.1038/jid.2008.321
Dashty, M., Motazacker, M. M., Levels, J., de Vries, M., Mahmoudi, M., Peppelenbosch, M. P., et al. (2014). Proteome of human plasma very low-density lipoprotein and low-density lipoprotein exhibits a link with coagulation and lipid metabolism. Thromb. Haemost. 111, 518–530. doi: 10.1160/TH13-02-0178
Deloukas, P., Kanoni, S., Willenborg, C., Farrall, M., Assimes, T. L., Thompson, J. R., et al. (2013). Large-scale association analysis identifies new risk loci for coronary artery disease. Nat. Genet. 45, 25–33. doi: 10.1038/ng.2480
Do, R., Stitziel, N. O., Won, H. H., Jørgensen, A. B., Duga, S., Angelica Merlini, P., et al. (2015). Exome sequencing identifies rare LDLR and APOA5 alleles conferring risk for myocardial infarction. Nature 518, 102–106. doi: 10.1038/nature13917
Gurung, R. L., Yiamunaa, M., Liu, S., Liu, J. J., Chan, C., Choo, R. W. M., et al. (2019). Association of haptoglobin phenotype with incident acute myocardial infarction in Chinese patients with type 2 diabetes. Cardiovasc. Diabetol. 18:65. doi: 10.1186/s12933-019-0867-4
Gusev, A., Ko, A., Shi, H., Bhatia, G., Chung, W., Penninx, B. W., et al. (2016). Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245–252. doi: 10.1038/ng.3506
Heng, B. C., Gribbon, P. M., Day, A. J., and Hardingham, T. E. (2008). Hyaluronan binding to link module of TSG-6 and to G1 domain of aggrecan is differently regulated by pH. J. Biol. Chem. 283, 32294–32301. doi: 10.1074/jbc.M804155200
Hernandez-Resendiz, S., Prunier, F., Girao, H., Dorn, G., and Hausenloy, D. J. (2020). Targeting mitochondrial fusion and fission proteins for cardioprotection. J. Cell. Mol. Med. 24, 6571–6585. doi: 10.1111/jcmm.15384
Hilchie, A. L., Wuerth, K., and Hancock, R. E. (2013). Immune modulation by multifaceted cationic host defense (antimicrobial) peptides. Nat. Chem. Biol. 9, 761–768. doi: 10.1038/nchembio.1393
Huang, C., Andres, A. M., Ratliff, E. P., Hernandez, G., Lee, P., and Gottlieb, R. A. (2011). Preconditioning involves selective mitophagy mediated by Parkin and p62/SQSTM1. PLoS One 6:e20975. doi: 10.1371/journal.pone.0020975
Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Iyer, R., Jenkinson, C. P., Vockley, J. G., Kern, R. M., Grody, W. W., and Cederbaum, S. (1998). The human arginases and arginase deficiency. J. Inherit. Metab. Dis. 21(Suppl. 1), 86–100. doi: 10.1023/a:1005313809037
Kubli, D. A., Zhang, X., Lee, Y., Hanna, R. A., Quinsay, M. N., Nguyen, C. K., et al. (2013). Parkin protein deficiency exacerbates cardiac injury and reduces survival following myocardial infarction. J. Biol. Chem. 288, 915–926. doi: 10.1074/jbc.M112.411363
Lande, R., Gregorio, J., Facchinetti, V., Chatterjee, B., Wang, Y. H., Homey, B., et al. (2007). Plasmacytoid dendritic cells sense self-DNA coupled with antimicrobial peptide. Nature 449, 564–569. doi: 10.1038/nature06116
Lee, I., and Huttemann, M. (2014). Energy crisis: the role of oxidative phosphorylation in acute inflammation and sepsis. Biochim. Biophys. Acta 1842, 1579–1586. doi: 10.1016/j.bbadis.2014.05.031
Lee, R. H., Pulin, A. A., Seo, M. J., Kota, D. J., Ylostalo, J., Larson, B. L., et al. (2009). Intravenous hMSCs improve myocardial infarction in mice because cells embolized in lung are activated to secrete the anti-inflammatory protein TSG-6. Cell Stem Cell 5, 54–63. doi: 10.1016/j.stem.2009.05.003
Lipovetzky, N., Hod, H., Roth, A., Kishon, Y., Sclarovsky, S., and Green, M. S. (2007). Emotional events and anger at the workplace as triggers for a first event of the acute coronary syndrome: a case-crossover study. Isr. Med. Assoc. J. 9, 310–315.
Lopez-Crisosto, C., Pennanen, C., Vasquez-Trincado, C., Morales, P. E., Bravo-Sagua, R., Quest, A. F. G., et al. (2017). Sarcoplasmic reticulum-mitochondria communication in cardiovascular pathophysiology. Nat. Rev. Cardiol. 14, 342–360. doi: 10.1038/nrcardio.2017.23
Luo, S., Au Yeung, S. L., Zhao, J. V., Burgess, S., and Schooling, C. M. (2019). Association of genetically predicted testosterone with thromboembolism, heart failure, and myocardial infarction: mendelian randomisation study in UK Biobank. BMJ 364:l476. doi: 10.1136/bmj.l476
Melamed-Frank, M., Lache, O., Enav, B. I., Szafranek, T., Levy, N. S., Ricklis, R. M., et al. (2001). Structure-function analysis of the antioxidant properties of haptoglobin. Blood 98, 3693–3698. doi: 10.1182/blood.v98.13.3693
Milner, C. M., Higman, V. A., and Day, A. J. (2006). TSG-6: a pluripotent inflammatory mediator? Biochem. Soc. Transac. 34(Pt 3), 446–450. doi: 10.1042/bst0340446
Mittleman, M. A., Maclure, M., Sherwood, J. B., Mulry, R. P., Tofler, G. H., Jacobs, S. C., et al. (1995). Triggering of acute myocardial infarction onset by episodes of anger. determinants of myocardial infarction onset study investigators. Circulation 92, 1720–1725. doi: 10.1161/01.cir.92.7.1720
Mittleman, M. A., Maclure, M., Tofler, G. H., Sherwood, J. B., Goldberg, R. J., and Muller, J. E. (1993). Triggering of acute myocardial infarction by heavy physical exertion. protection against triggering by regular exertion. determinants of myocardial infarction onset study investigators. N. Engl. J. Med. 329, 1677–1683. doi: 10.1056/nejm199312023292301
Mittleman, M. A., Mintzer, D., Maclure, M., Tofler, G. H., Sherwood, J. B., and Muller, J. E. (1999). Triggering of myocardial infarction by cocaine. Circulation 99, 2737–2741. doi: 10.1161/01.cir.99.21.2737
Murray, C. J., Barber, R. M., Foreman, K. J., Abbasoglu Ozgoren, A., Abd-Allah, F., Abera, S. F., et al. (2015). Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990-2013: quantifying the epidemiological transition. Lancet 386, 2145–2191. doi: 10.1016/s0140-6736(15)61340-x
Nikpay, M., Goel, A., Won, H. H., Hall, L. M., Willenborg, C., Kanoni, S., et al. (2015). A comprehensive 1,000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat. Genet. 47, 1121–1130. doi: 10.1038/ng.3396
Park, H. J., Noh, J. H., Eun, J. W., Koh, Y. S., Seo, S. M., Park, W. S., et al. (2015). Assessment and diagnostic relevance of novel serum biomarkers for early decision of ST-elevation myocardial infarction. Oncotarget 6, 12970–12983. doi: 10.18632/oncotarget.4001
Qi, L., Ma, J., Qi, Q., Hartiala, J., Allayee, H., and Campos, H. (2011). Genetic risk score and risk of myocardial infarction in Hispanics. Circulation 123, 374–380. doi: 10.1161/CIRCULATIONAHA.110.976613
Santhanam, L., Lim, H. K., Lim, H. K., Miriel, V., Brown, T., Patel, M., et al. (2007). Inducible NO synthase dependent S-nitrosylation and activation of arginase1 contribute to age-related endothelial dysfunction. Circ. Res. 101, 692–702. doi: 10.1161/circresaha.107.157727
Silbiger, V. N., Luchessi, A. D., Hirata, R. D., Lima-Neto, L. G., Cavichioli, D., Carracedo, A., et al. (2013). Novel genes detected by transcriptional profiling from whole-blood cells in patients with early onset of acute coronary syndrome. Clin. Chim. Acta 421, 184–190. doi: 10.1016/j.cca.2013.03.011
Sin, Y. Y., Baron, G., Schulze, A., and Funk, C. D. (2015). Arginase-1 deficiency. J. Mol. Med. 93, 1287–1296. doi: 10.1007/s00109-015-1354-3
Suleiman, M., Aronson, D., Asleh, R., Kapeliovich, M. R., Roguin, A., Meisel, S. R., et al. (2005). Haptoglobin polymorphism predicts 30-day mortality and heart failure in patients with diabetes and acute myocardial infarction. Diabetes 54, 2802–2806. doi: 10.2337/diabetes.54.9.2802
Sun, W., Zheng, Y., Lu, Z., Wang, H., Feng, Z., Wang, J., et al. (2014). LL-37 attenuates inflammatory impairment via mTOR signaling-dependent mitochondrial protection. Int. J. Biochem. Cell Biol. 54, 26–35. doi: 10.1016/j.biocel.2014.06.015
Tang, R., Wang, X., Zhou, J., Zhang, F., Zhao, S., Gan, Q., et al. (2020). Defective arginine metabolism impairs mitochondrial homeostasis in Caenorhabditis elegans. J. Genet. Genom. 47, 145–156. doi: 10.1016/j.jgg.2020.02.007
Thériault, S., Gaudreault, N., Lamontagne, M., Rosa, M., Boulanger, M. C., Messika-Zeitoun, D., et al. (2018). A transcriptome-wide association study identifies PALMD as a susceptibility gene for calcific aortic valve stenosis. Nat. Commun. 9:988. doi: 10.1038/s41467-018-03260-6
Wisniewski, H. G., and Vilcek, J. (2004). Cytokine-induced gene expression at the crossroads of innate immunity, inflammation and fertility: TSG-6 and PTX3/TSG-14. Cytokine Growth Factor Rev. 15, 129–146. doi: 10.1016/j.cytogfr.2004.01.005
Wu, X., He, L., Cai, Y., Zhang, G., He, Y., Zhang, Z., et al. (2013). Induction of autophagy contributes to the myocardial protection of valsartan against ischemia-reperfusion injury. Mol. Med. Rep. 8, 1824–1830. doi: 10.3892/mmr.2013.1708
Yu, P., Zhang, J., Yu, S., Luo, Z., Hua, F., Yuan, L., et al. (2015). Protective effect of sevoflurane postconditioning against cardiac ischemia/reperfusion injury via ameliorating mitochondrial impairment, oxidative stress and rescuing autophagic clearance. PLoS One 10:e0134666. doi: 10.1371/journal.pone.0134666
Zasloff, M. (2002). Antimicrobial peptides of multicellular organisms. Nature 415, 389–395. doi: 10.1038/415389a
Zhang, L., Liu, L., Ma, M., Cheng, S., Cheng, B., Li, P., et al. (2019). Integrative analysis of transcriptome-wide association study data and mRNA expression profiles identified candidate genes and pathways associated with atrial fibrillation. Heart Vessels 34, 1882–1888. doi: 10.1007/s00380-019-01418-w
Zhang, R., Ji, Z., Qu, Y., Yang, M., Su, Y., Zuo, W., et al. (2020). Clinical value of ARG1 in acute myocardial infarction patients: bioinformatics-based approach. Biomed. Pharmacother. 121:109590. doi: 10.1016/j.biopha.2019.109590
Keywords: acute myocardial infarction, transcriptome-wide association study, genome-wide association study, mRNA expression profiles, biological process
Citation: Chen G, Liu L, Li H, Lun Z, Mai Z, Lai W, Chen E, Zhou C, Yu S, Yang J, Chen S, Chen J and Liu Y (2021) Integrative Analysis of Transcriptome-Wide Association Study and mRNA Expression Profiles Identified Candidate Genes and Pathways Associated With Acute Myocardial Infarction. Front. Genet. 12:616492. doi: 10.3389/fgene.2021.616492
Received: 12 October 2020; Accepted: 08 January 2021;
Published: 02 February 2021.
Edited by:
Jordi Pérez-Tur, Superior Council of Scientific Investigations (CSIC), SpainReviewed by:
Rita Selvatici, University of Ferrara, ItalyBruna De Felice, University of Campania Luigi Vanvitelli, Italy
Wei Vivian Li, Rutgers, The State University of New Jersey, United States
Copyright © 2021 Chen, Liu, Li, Lun, Mai, Lai, Chen, Zhou, Yu, Yang, Chen, Chen and Liu. 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: Shiqun Chen, shiqunchen@126.com; Jiyan Chen, chenjiyandr@126.com; Yong Liu, liuyong@gdph.org.cn
†These authors have contributed equally to this work