AUTHOR=Gong Weiming , Guo Ping , Liu Lu , Guan Qingbo , Yuan Zhongshang 

TITLE=Integrative Analysis of Transcriptome-Wide Association Study and mRNA Expression Profiles Identifies Candidate Genes Associated With Idiopathic Pulmonary Fibrosis

JOURNAL=Frontiers in Genetics

VOLUME=11

YEAR=2020

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2020.604324

DOI=10.3389/fgene.2020.604324

ISSN=1664-8021

ABSTRACT=<p>Idiopathic pulmonary fibrosis (IPF) is a type of scarring lung disease characterized by a chronic, progressive, and irreversible decline in lung function. The genetic basis of IPF remains elusive. A transcriptome-wide association study (TWAS) of IPF was performed by FUSION using gene expression weights of three tissues combined with a large-scale genome-wide association study (GWAS) dataset, totally involving 2,668 IPF cases and 8,591 controls. Significant genes identified by TWAS were then subjected to gene ontology (GO) and pathway enrichment analysis. The overlapped GO terms and pathways between enrichment analysis of TWAS significant genes and differentially expressed genes (DEGs) from the genome-wide mRNA expression profiling of IPF were also identified. For TWAS significant genes, protein–protein interaction (PPI) network and clustering modules analyses were further conducted using STRING and Cytoscape. Overall, TWAS identified a group of candidate genes for IPF under the Bonferroni corrected <italic>P</italic> value threshold (0.05/14929 = 3.35 × 10<sup>–6</sup>), such as <italic>DSP</italic> (<italic>P</italic><sub><italic>TWAS</italic></sub> = 1.35 × 10<sup>–29</sup> for lung tissue), <italic>MUC5B</italic> (<italic>P</italic><sub><italic>TWAS</italic></sub> = 1.09 × 10<sup>–28</sup> for lung tissue), and <italic>TOLLIP</italic> (<italic>P</italic><sub><italic>TWAS</italic></sub> = 1.41 × 10<sup>–15</sup> for whole blood). Pathway enrichment analysis identified multiple candidate pathways, such as herpes simplex infection (<italic>P</italic> value = 7.93 × 10<sup>–5</sup>) and antigen processing and presentation (<italic>P</italic> value = 6.55 × 10<sup>–5</sup>). 38 common GO terms and 8 KEGG pathways shared by enrichment analysis of TWAS significant genes and DEGs were identified. In the PPI network, 14 genes (<italic>DYNLL1</italic>, <italic>DYNC1LI1</italic>, <italic>DYNLL2</italic>, <italic>HLA-DRB5</italic>, <italic>HLA-DPB1</italic>, <italic>HLA-DQB2</italic>, <italic>HLA-DQA2</italic>, <italic>HLA-DQB1</italic>, <italic>HLA-DRB1</italic>, <italic>POLR2L</italic>, <italic>CENPP</italic>, <italic>CENPK</italic>, <italic>NUP133</italic>, and <italic>NUP107</italic>) were simultaneously detected by hub gene and module analysis. In conclusion, through integrative analysis of TWAS and mRNA expression profiles, we identified multiple novel candidate genes, GO terms and pathways for IPF, which contributes to the understanding of the genetic mechanism of IPF.</p>