Target biomarkers for H2 at both the protein and genome levels are still unclear. In this study, quantitative proteomics acquired from a mouse model were first analyzed. At the same time, functional pathway analysis helped identify functional pathways at the protein level. Then, bioinformatics on mRNA sequencing data were conducted between sepsis and normal mouse models. Differential expressional genes with the closest relationship to disease status and development were identified through module correlation analysis. Then, common biomarkers in proteomics and transcriptomics were extracted as target biomarkers. Through analyzing expression quantitative trait locus (eQTL) and genome-wide association studies (GWAS), colocalization analysis on Apoa2 and sepsis phenotype was conducted by summary-data-based Mendelian randomization (SMR). Then, two-sample and drug-target, syndrome Mendelian randomization (MR) analyses were all conducted using the Twosample R package. For protein level, protein quantitative trait loci (pQTLs) of the target biomarker were also included in MR. Animal experiments helped validate these results. As a result, Apoa2 protein or mRNA was identified as a target biomarker for H2 with a protective, causal relationship with sepsis. HDL and type 2 diabetes were proven to possess causal relationships with sepsis. The agitation and inhibition of Apoa2 were indicated to influence sepsis and related syndromes. In conclusion, we first proposed Apoa2 as a target for H2 treatment.
Keywords: GWAS; Mendelian randomization; PWAS; genome; molecular hydrogen; proteomics; septic lung injury.