Transcriptome analysis
| Gene | Log2 fold change | Log CPM | p-value | False discovery rate |
| HSPA7 | −6.6745634 | 2.94357766 | 2.8625E−56 | 4.325E−52 |
| CCL26 | −4.0022854 | 1.53369684 | 4.5982E−36 | 3.4737E−32 |
| RN7SKP48 | −8.4063178 | 0.00054165 | 4.2481E−34 | 2.1395E−30 |
| ARC | −6.0359691 | 2.80438403 | 1.3802E−29 | 5.2134E−26 |
| ANKRD1 | −3.5612232 | 2.00265271 | 1.2001E−17 | 3.6264E−14 |
| CRYAB | −1.4560917 | 7.69968071 | 1.2704E−15 | 3.199E−12 |
| HSPA6 | −6.1119842 | 5.60205649 | 1.0489E−14 | 2.264E−11 |
| VTRNA1-3 | −5.0654191 | 2.21568985 | 8.9014E−13 | 1.6811E−09 |
| SPP1 | −2.8175604 | 0.80469766 | 3.4552E−12 | 5.8005E−09 |
| HSPA1B | −3.291132 | 8.25091346 | 1.2859E−10 | 1.9428E−07 |
| XIRP1 | −4.3644087 | 0.39829287 | 4.8336E−10 | 6.6391E−07 |
| NEFM | −2.8494359 | 2.86373374 | 9.6932E−10 | 1.2205E−06 |
| FOSB | −2.3239735 | 2.72522439 | 2.473E−09 | 2.8742E−06 |
| RRAD | −2.2480685 | 4.46413378 | 5.0727E−09 | 5.4745E−06 |
| HSPA1A | −2.5310616 | 8.45620337 | 2.8287E−08 | 2.8492E−05 |
| IL1RL1 | −2.2992064 | 1.98026678 | 6.8739E−08 | 6.4911E−05 |
| RASD1 | −2.2417413 | 1.63816909 | 8.5819E−08 | 7.6273E−05 |
| VTRNA1-2 | −4.2152004 | 2.37699757 | 2.3762E−07 | 0.00019946 |
| IL16 | 0.9270627 | 3.3688119 | 3.1214E−07 | 0.00024822 |
| SPON1 | 0.76087863 | 5.63037004 | 1.4456E−06 | 0.00109207 |
| COL15A1 | 1.0153401 | 6.86806305 | 1.5531E−06 | 0.00111743 |
| RP11-96K19.4 | −2.3607353 | −0.5695676 | 1.8955E−06 | 0.0012725 |
| ABCA8 | 0.6821761 | 4.63139662 | 1.9371E−06 | 0.0012725 |
| HIST1H2BK | −0.7249494 | 4.90831468 | 2.3208E−06 | 0.00146106 |
| HBEGF | −1.0262606 | 2.6512004 | 2.5521E−06 | 0.00150211 |
| UPP1 | −0.6218694 | 4.65882924 | 2.5849E−06 | 0.00150211 |
| MASP1 | 0.81818236 | 7.05487073 | 9.1183E−06 | 0.00510251 |
| C1RL | 0.44446462 | 5.41931433 | 1.1084E−05 | 0.00598114 |
| DUSP5 | −0.4901746 | 4.66863971 | 1.2827E−05 | 0.00668276 |
| TNXB | 0.57515764 | 5.54831679 | 1.3871E−05 | 0.00698593 |
| THBD | −0.5948758 | 5.0121393 | 2.0903E−05 | 0.01018796 |
| ADAMTS8 | 0.79724554 | 5.51175195 | 2.3182E−05 | 0.01065668 |
| MXRA5 | 0.62849027 | 7.08606778 | 2.3276E−05 | 0.01065668 |
| DUSP1 | −0.7788436 | 7.04579431 | 2.7687E−05 | 0.01230367 |
| MMP12 | −2.2223219 | 3.62749486 | 2.9065E−05 | 0.01254702 |
| SLITRK6 | 0.892711 | 5.72816647 | 3.1127E−05 | 0.01306376 |
| LINC00324 | −2.2447718 | 0.41686191 | 3.659E−05 | 0.01494156 |
| PLXDC2 | 0.52042951 | 5.75171529 | 4.2403E−05 | 0.01685985 |
| ABCA6 | 0.7034962 | 3.44538556 | 4.5537E−05 | 0.01764146 |
| RP11-64B16.2 | −2.7001276 | −0.3550856 | 4.9958E−05 | 0.01887036 |
| LAMA4 | 0.62652521 | 8.74391096 | 5.2261E−05 | 0.01925895 |
| VTRNA2-1 | −5.2769676 | 0.22835068 | 5.4686E−05 | 0.01928778 |
| METTL7A | 0.47147354 | 5.2757973 | 5.4893E−05 | 0.01928778 |
| LCE3D | −1.2948825 | 1.28528292 | 5.8792E−05 | 0.02018826 |
| WISP1 | 0.823069 | 3.56387389 | 6.5606E−05 | 0.02183693 |
| LCE3E | −1.2785048 | 1.1649956 | 6.6483E−05 | 0.02183693 |
| GADD45B | −1.2909076 | 5.87248869 | 7.0692E−05 | 0.02272509 |
| RGS16 | −1.3795921 | 0.16021976 | 8.2909E−05 | 0.02609742 |
| RP11-367F23.2 | −2.5075731 | 0.32741435 | 8.6844E−05 | 0.026778 |
| AF131217.1 | 0.51616864 | 4.24716742 | 9.5197E−05 | 0.02876658 |
CPM: counts per million. Nine primary epithelial cell isolates obtained before and after experimental bronchial thermoplasty were exposed to 37°C versus 65°C and next-generation sequencing was performed as described in the Methods section. Transcripts with the highest fold change between groups are shown at the top of the table. p-values have been adjusted using the Benjamini–Hochberg false discovery rate to correct for multiple testing.