A molecular single-cell lung atlas of lethal COVID-19.
影響因子： 42.778
PMID：33915568
期刊年卷：Nature 2021 Apr 29;
綜合性期刊一區 綜合性期刊 Q1 1/64
DOI：10.1038/s41586-021-03569-1
https://www.nature.com/articles/s41586-021-03569-1#Abs1
摘要
呼吸衰竭是嚴重 SARS-CoV-2 感染患者死亡的主要原因，但對肺組織水平的宿主反應知之甚少。在這里，我們對 19 名死于 COVID-19 并接受快速尸檢的個體和 7 名對照個體的肺中的約 116,000 個核進行了單核 RNA 測序。綜合分析確定了細胞組成、轉錄細胞狀態和細胞間相互作用的重大變化，從而提供了對致命 COVID-19 生物學的深入了解。COVID-19 患者的肺部高度發炎，異常活化的單核細胞衍生的巨噬細胞和肺泡巨噬細胞密集浸潤，但 T 細胞反應受損。與其他病毒和細菌引起的肺炎相比，單核細胞/巨噬細胞來源的白細胞介素-1β 和上皮細胞來源的白細胞介素-6 是 SARS-CoV-2 感染的獨特特征。肺泡 2 型細胞采用炎癥相關的瞬時祖細胞狀態，未能完全轉變為肺泡 1 型細胞，導致肺再生受損。此外，我們確定了最近描述的擴展CTHRC1+病理性成纖維細胞導致 COVID-19 中迅速發生的肺纖維化。蛋白質活性和配體-受體相互作用的推斷確定了破壞有害回路的假定藥物靶標。該圖譜可以解析致命的 COVID-19，可能有助于我們了解 COVID-19 幸存者的長期并發癥，并為治療開發提供重要資源。
引言
在全球范圍內，由 SARS-CoV-2 感染引起的 COVID-19 大流行已導致超過 1.45 億病例（美國為 3200 萬）和 310 萬例死亡（美國為 570,000；截至 4 月 26 日的數據，2021) 1 . 大約 15% 的感染者會發展為嚴重疾病，這可能表現為急性呼吸窘迫綜合征 (ARDS)，并與大量發病率和死亡率相關2 , 4。
以前，單細胞RNA測序（scRNA-SEQ）健康個體的分析已經揭示所必需的SARS-CoV的-2條目宿主受體的組織分布5，6，7，和從患者的支氣管肺泡灌洗液和血液檢查不同嚴重程度的 COVID-19 已經確定了 SARS-CoV-2 感染對免疫反應和細胞因子失調的影響8 , 9 , 10 , 11 , 12. 然而，由于獲取患者組織的實際限制，SARS-CoV-2在肺組織水平的影響尚不清楚。一系列尸檢研究檢查了死于 COVID-19 的個體的福爾馬林固定石蠟包埋 (FFPE) 組織切片，擴展了我們對病毒器官性的理解，但這些研究在發現潛力方面受到低重分析（例如例如，免疫組織化學）和/或延長的驗尸間隔（的PMI），該RNA的質量產生不利的影響13，14，15。
我們建立了快速尸檢計劃，并根據機構審查委員會批準的協議，在死亡數小時內從 COVID-19 患者身上收集速凍器官標本。我們對死于 COVID-19 的個體和對照個體的肺樣本進行了單核 RNA-seq (snRNA-seq)，以構建一個圖譜，可以深入了解 COVID-19 的病理生理學，并為進一步研究提供關鍵資源。
COVID-19 中的肺細胞景觀
COVID-19 隊列由 19 名患者（12 名男性和 7 名女性）組成，他們的中位年齡為 72 歲（范圍，58 歲至超過 89 歲）（補充表1，擴展數據圖1a），并接受了快速尸檢。中位數尸檢間隔 (PMI) 為 4 小時（范圍，2-9 小時）。所有患者都存在與嚴重 COVID-19 16風險增加相關的潛在合并癥（補充表1）。對照組包括 7 名個體（4 男 3 女），中位年齡為 70 歲（范圍，67 至 79 歲），他們在 COVID-19 之前的時代接受了肺切除術或活檢（補充表1）。
使用 snRNA-seq 17和集成的質量控制管道（參見 方法），我們生成了一個肺圖譜，其中包含 116,314 個細胞核，其中 79,636 個來自 COVID-19 感染的肺，36,678 個來自對照肺（圖1a）。我們使用三管齊下的方法進行細胞類型識別：聚類標記的無偏識別、使用來自報告圖譜的特征發現細胞類型，以及使用專家知識手動管理以對細胞群和細胞狀態進行細分（參見 方法）。我們報告了三個粒度級別的細胞類型分配：主要細胞類型、中間粒度和細粒度（補充表2）。我們使用UMAP（圖1b、c、擴展數據圖1b-d）來可視化降維數據。我們確定了九種主要細胞類型：上皮細胞 ( n = 30,070 個細胞)、髓細胞 ( n = 29,632)、成纖維細胞 ( n = 22,909)、內皮細胞 ( n = 5,386)、T 和自然殺傷 (NK) 淋巴細胞（n = 16,751)、B 淋巴細胞和漿細胞 ( n = 7,236)、神經元細胞 ( n = 2,017)、肥大細胞 ( n = 1,464) 和抗原呈遞細胞 (APC；主要是樹突細胞) ( n = 849）。在最細粒度的水平上，我們確定了 41 種不同的細胞類型（補充表2）。
圖 1：研究設計和細胞景觀
a，研究設計概述。b，UMAP 中的主要簇和各自的細胞類型分配。c，與b 中嵌入相同的細胞的起源。d，對照（n = 7）和 COVID-19 肺（n = 19）中主要細胞類型的比例。中線，中位數；框邊緣，第 25 個和第 75 個百分位數；胡須，不超過 ±1.5 × 四分位距 (IQR) 的最極端點。Wilcoxon 秩和檢驗。
我們發現 COVID-19 和對照肺之間的細胞分數在全局（圖1d）和免疫區室和非免疫區室（擴展數據圖2a-c）之間存在顯著差異。由于肺泡 II 型 (AT2) 和 I 型 (AT1) 細胞的丟失，以及單核細胞/巨噬細胞、成纖維細胞和神經元細胞的增加，上皮細胞區室減少；這些觀察結果與供體性別無關（擴展數據圖3a、b）。
我們發現COVID-19 和對照肺之間的ACE2、CD147（也稱為BSG）、NPR1、TMPRSS2、FURIN或CTSL的表達沒有重大差異（擴展數據圖3c-f）。這表明細胞類型比例的變化與對病毒進入很重要的受體或假定蛋白酶的表達無關，盡管我們不能排除病毒介導的細胞死亡選擇性消耗這些基因高表達的細胞的可能性。我們在兩名患者（補充表3）中檢測到 SARS-CoV-2 讀數，其中一名患有 HIV/AIDS（CD4 +入院時T 細胞計數 29/mm ；在 28 個細胞中檢測到 662 個獨特的分子標識符），這表明原則上可以捕獲病毒讀數。
髓樣細胞的異常激活
髓樣細胞代表了COVID-19肺中的主要細胞成分，而且比對照肺（圖更普遍有1D，擴展數據圖2的A，C，圖4a）。我們鑒定了單核細胞 ( n = 3,176)、單核細胞衍生的巨噬細胞 (MDMs; n = 9,534)、過渡 MDMs ( n = 4,203) 和常駐肺泡巨噬細胞 (AMs; n = 12,511)，它們作為擴散成分的不同軌跡被回收（DC）分析和COVID-19肺（圖更頻繁2a-c中，擴展數據圖4B-I，補充表2，4，5）。來自 COVID-19 患者的骨髓細胞被高度異常激活。例如，COVID-19 肺中的 MDM 差異表達激活基因（例如CTSB、CTSD、CTSZ、PSAP）和兩種長鏈非編碼 RNA，NEAT1和MALAT1，它們與巨噬細胞異常激活和 T 細胞免疫受損相關18（擴展數據圖5a，補充表5）。AMs 由胎兒單核細胞產生并可以自我更新19，在 COVID-19 肺中富集和高度激活（圖2c，擴展數據圖5a））。值得注意的是，COVID-19 AM 顯示出腫瘤相關巨噬細胞受體 AXL 的 mRNA 和蛋白質表達強烈降低（圖2d，擴展數據圖5b、c），這是一種受體酪氨酸激酶，對于協調清除凋亡細胞（胞吞作用）很重要) 和隨后的組織再生過程中的抗炎調節20。這些數據表明，髓樣細胞是 COVID-19 炎癥失調的主要來源。
圖 2：COVID-19 中的免疫反應
a，突出免疫細胞簇的UMAP投影。b，使用前三個 DC 可視化髓樣細胞。插圖表示組分配。c，對照（n = 7）和 COVID-19 肺（n = 19）的髓樣細胞分數。中線，中位數；框邊緣，第 25 個和第 75 個百分位數；晶須：不超過±1.5 × IQR 的最極端點。Wilcoxon 秩和檢驗。d。對照和 COVID-19 肺組織中 CD169、AXL 和 DAPI（大圖）的代表性免疫熒光染色；頂部，帶有疊加層的選定區域；底部，個人頻道。比例尺，20 微米。e , f , Top 20 經常檢測到的IGHV –COVID -19 ( e ) 中的IGLV組合和相應的組注釋 (f )。先前描述的抗 RBD 抗體的組合21。g，T/NK細胞的UMAP；插圖，小組作業。h，i，GZMB ( h ) 和MKI67 ( i ) 在與g相同的嵌入中的RNA 表達（對數歸一化）。
血漿和 T 細胞反應
為了深入了解肺中針對 SARS-CoV-2 感染的體液免疫，我們基于每個細胞的鏈和同種型, 通過確定可變重 ( IGHV ) 和輕 ( IGLV ) 的mRNA 共表達來鑒定漿細胞并重建免疫球蛋白（擴展數據圖6a-c）（參見 方法；擴展數據圖6d-k，補充表6）。IGHV1-18–IGLV3-20產生針對 SARS-CoV-2 刺突蛋白的受體結合域 (RBD) 的中和抗體 (S309)21，是常見的IGHV–IGLV 之一組合，這表明發生了協調的抗體反應（圖2e，f，擴展數據圖6l，m）。在 T/NK 細胞簇（圖2g）中，我們區分了 CD8 + T 細胞（n = 3,561）、T 調節（T reg）細胞（n = 649）、其他 CD4 + T 細胞（n = 7,586）和NK 細胞 ( n = 2,141)。我們發現 COVID-19 肺中 T 細胞豐度沒有顯著增加，只有與 T 細胞活化和組織駐留相關的細胞因子和程序適度上調（圖2g-i，擴展數據圖7a-i））。盡管免疫反應模式變化很大（擴展數據圖7j，k），但這些數據表明，在主要保留的體液免疫反應的背景下，受損的 T 細胞反應可能會導致 COVID-19 的致死結果。
肺泡上皮再生受損
在上皮區室中，我們確定了肺泡上皮細胞（AT1 和 AT2 細胞；n = 20,949）、氣道上皮細胞（基底細胞、纖毛細胞、棒狀細胞、杯狀細胞和粘液細胞；n = 7,223），這是一個以炎癥表達為特征的簇和細胞周期基因，包括IRF8、B2M、MKI67和TOP2A（循環上皮；n = 609），以及一個顯示細胞外基質 (ECM) 成分COL6A3、COL1A2和COL3A1（ECM高上皮）高表達的簇; n= 1,179) (圖3a, b，擴展數據圖8A-C，補充表2，7）。
圖 3：肺再生受損和炎癥來源
a、b、研究的肺泡和氣道上皮細胞的 UMAP ( a ) 和相應的組分配 ( b )。c，來自 COVID-19 和對照肺的 AT1 和 AT2 細胞的差異基因表達（對數歸一化，縮放；參見 方法）。列，單個單元格；行，最高調節基因的表達。左欄，AT1（紫色）和 AT2（粉色）細胞的譜系標記。顏色編碼的頂部泳道表示簽名的表達強度（對數歸一化；參見 方法）和組分配，如右側所示。exp.，表達。d , e , ETV5和CAV1 的小提琴圖分別在 AT2 和 AT1 細胞中的 mRNA 表達（對數歸一化）；帶有 Bonferroni 校正的 Wilcoxon 秩和檢驗。f，UMAP嵌入AT1和AT2細胞并鑒定DATP；插圖表示小組分配。g，AT1和AT2細胞中DATP特征表達（對數歸一化）的小提琴圖。Wilcoxon 秩和檢驗。h，前三個 DC 顯示了 AT2 和 AT1 細胞和 DATP 的主要軌跡、DATP 特征的表達和組分配（插圖）。i，對照 ( n = 7) 和 COVID-19 肺 ( n = 19)中 DATP 和 AT 細胞的比例。中線，中位數；框邊緣，第 25 個和第 75 個百分位數；晶須，不超過 ±1.5 × IQR 的最極端點。Wilcoxon 秩和檢驗。j，對照和 COVID-19 肺組織中 pro-SPC、KRT8 和 DAPI 的代表性免疫熒光染色；頂部，具有覆蓋的代表性區域；底部，帶有選定區域各個通道的小圖像。比例尺，50 微米。k , l，健康肺組織中 IL-1β ( k ) 和 IL-6 ( l ) 的組織質量細胞計數定量，以及來自具有不同感染病因的供體的樣本。每個點代表感興趣區域 (ROI) 中 IL-1β 和 IL-6 的量化；帶有 Benjamini-Hochberg 錯誤發現率 (FDR) 調整的雙邊 Mann-Whitney U檢驗。
AT2 細胞在肺再生22期間充當 AT1 細胞的祖細胞。對照肺中的 AT2 和 T1 細胞形成不同的簇（圖3a、b）并證明了差異基因表達 (DGE) 分析的預期變化，包括AT2 細胞中譜系標記SFTPC和SFTPB的表達，以及AT1 中的CLIC5和AGER細胞（圖3c，補充表7）。相比之下，COVID-19 肺中 AT2 和 AT1 細胞的聚類不那么離散，很大一部分細胞不與其對照對應物重疊（圖3b））。來自 COVID-19 肺的 AT2 和 AT1 細胞均顯示出定義標志物的整體表達降低（圖3c）。COVID-19 AT2 細胞顯示ETV5表達降低（圖3d），這是一種維持 AT2 細胞身份所需的轉錄因子。降低的ETV5表達與向 AT1 細胞的分化相關23，表明 AT2 細胞已啟動再生程序（圖3d，擴展數據圖8d）。CAV1是 AT1 晚期成熟的標志物24，在來自 COVID-19 肺的 AT1 細胞中以顯著較低的水平表達（圖3e））。總體而言，這些數據表明 COVID-19 肺中 AT2 到 AT1 細胞的轉化不完全。
最近的研究表明，炎癥可以誘導細胞狀態，其特征是不能完全轉變為 AT1 細胞；這已被稱為“損傷相關瞬態祖細胞”（DATPs），“肺泡分化中間”（ADI），或“預AT1移行細胞狀態”（PATS）25，26，27（以下簡稱為DATPs）。我們使用 DATP 標記基因（KRT8、CLDN4和CDKN1A）25 的表達來開發 DATP 特征（參見 方法；擴展數據圖8e-h，補充表8）并發現來自 COVID-19 肺的肺泡上皮細胞在表達該特征方面的得分明顯高于來自對照肺的細胞（圖3f、g，擴展數據圖8i）。DC 分析將主要軌跡從 AT2 分離到 AT1 細胞，而 DATP 主要位于 AT2 和 AT1 細胞之間（圖3h，擴展數據圖8j-n）。與分化的 AT2 或 AT1 細胞相比，DATP 的基因集富集分析 (GSEA) 顯示 TNFα 和 p53 信號的富集，以及通過 HIF-1α 的缺氧反應（擴展數據圖8o），與 DATP 中涉及的途徑一致在小鼠模型中27. 與 p53 信號的過度表達一致，大多數 DATP 沒有經歷細胞分裂（擴展數據圖8p），表明它們在 DATP 細胞狀態中停滯。
DATP 在 COVID-19 中比對照肺更頻繁（圖3i）。相應組織的免疫熒光染色顯示 KRT8 +和 CLDN4 + DATP 在 COVID-19 肺中的頻率較高（圖3j，擴展數據圖8r，s），我們觀察到隨著時間的推移，AT1 細胞豐度逐漸喪失。癥狀發作至死亡（擴展數據圖8t）。總體而言，這些數據表明，除了病毒感染直接破壞肺泡上皮外，COVID-19 患者的肺再生過程也會受損。
我們接下來確定了導致 DATP 細胞狀態，更一般地說，導致 COVID-19 肺部高炎癥環境的炎癥來源。在 mRNA 水平捕獲炎性細胞因子白細胞介素 (IL)-1β（和其他）可能是有限的，因為 IL-1β 的生物活性形式在觸發 DATP 中起主要作用25，是通過從 pro-IL 裂解產生的-1β 炎癥小體激活時；因此，蛋白質水平評估提供了補充信息。為此，我們利用了最近發布的高分辨率成像質量細胞計數數據集，該數據集對來自 23 個個體（包括健康對照）的 237 個組織區域進行了分析；流感肺炎、細菌性肺炎或ARDS患者；和 10 名死于 COVID-19 的患者28. IL-1β 在 COVID-19 患者的單核細胞和巨噬細胞中的表達比健康個體或其他疾病組的患者更強（圖3k，擴展數據圖9a-c）。IL-6 是 COVID-19 病理生理學中調用的另一種關鍵炎癥細胞因子，它在 COVID-19 患者的上皮細胞中更豐富，但與其他疾病組的患者相比，這些患者的巨噬細胞中沒有差異表達（圖 1）。3l，擴展數據圖9d-f）。最后，我們發現 COVID-19 患者的 I 型干擾素和干擾素反應基因在各種細胞類型（包括 AT2 細胞、單核細胞和巨噬細胞）中的表達強于對照供體（擴展數據圖9g、h） . 總之，這些數據表明，與其他病毒性或細菌性肺炎相比，髓源性 IL-1β 可能是 COVID-19 的一個顯著特征，并且可能有助于誘導和維持 DATP 細胞狀態。
COVID-19 中的異位簇狀細胞
中捕獲的氣道上皮細胞，我們回收四個不同的軌跡：KRT5+TP63+基底（? = 534），俱樂部（? = 1232），和杯狀細胞（? = 1757），和一個軌跡具有較少細胞（? = 110）這主要在 COVID-19 肺中發現，我們將其確定為推定的簇狀細胞（擴展數據圖10a-e）。簇狀細胞參與氣道炎癥和腸組織再生29，但它們在病毒性肺炎中的作用仍不清楚。簇絨單元的數量（CHAT +或 POU2F3 +) 在 COVID-19 患者的上呼吸道中增加了三倍，并且它們異位存在于 COVID-19 的肺實質中，但不存在于對照肺中（擴展數據圖10f-k）。為了開始闡明簇狀細胞在病毒性肺炎中的假定作用，我們用 PR8（一種實驗室適應的 H1N1 流感病毒株）感染了缺乏簇狀細胞的野生型和Pou2f3-/-小鼠（參見 方法）。與對照相比，Pou2f3-/-小鼠的肺顯示巨噬細胞浸潤減少和趨化基因（包括Ccl3和Ccl8) 也參與將髓細胞募集到死于 COVID-19 的個體的肺部（擴展數據圖9 g、h、11a-l）。盡管需要進一步研究它們的作用，但這些異位簇狀細胞可能有助于 COVID-19 的病理生理學（補充討論）。
病理性成纖維細胞和肺纖維化
COVID-19肺中的成纖維細胞明顯多于對照肺（圖1d）；α-平滑肌肌動蛋白 (α-SMA) 的免疫組織化學染色證實了這一發現（擴展數據圖12a-d）。纖維化程度（由 Sirius red 纖維化評分確定，參見 方法）與疾病持續時間相關（圖4a），表明 COVID-19 中肺纖維化隨時間增加。我們確定了五種成纖維細胞亞型：肺泡 ( n = 4,670)、外膜 ( n = 3,773)、病理性 ( n = 2,322)、中間病理性 ( n = 8,779) 和其他 ( n = 1,099)（圖4b）, 擴展數據圖12e )。成纖維細胞簇差異的主要驅動因素是與對照肺相比，COVID-19 肺中病理或中間病理成纖維細胞（以下統稱為 pFB）的頻率增加（圖4c，擴展數據圖12f）。pFB 強烈表達*****CTHRC1*****，這是最近描述的定義這些細胞的標志基因，以及病理性 ECM 3 基因，包括COL1A1和COL3A1（擴展數據圖12e，補充表9）。在小鼠模型和特發性肺纖維化 (IPF) 或硬皮病3患者中，pFB 是肺纖維化的關鍵驅動因素。它們增加的頻率表明 pFB 促進了 COVID-19 患者快速發展的肺纖維化。**
圖 4：COVID-19 中的病理性成纖維細胞和隨后的纖維化。
a，COVID-19 樣本（n = 16，參見 方法）中從癥狀出現到死亡和纖維化評分的天數測定系數 ( R 2 )。誤差帶，皮爾遜相關性的 95% se 間隔。b，成纖維細胞（FB）亞群的 UMAP；插圖表示小組分配。path., 病態的。c，對照（n = 7）和COVID-19肺（n = 19）的所有成纖維細胞中病理成纖維細胞的分數。中線，中位數；框邊緣，第 25 個和第 75 個百分位數；晶須，不超過 ±1.5 × IQR 的最極端點。Wilcoxon 秩和檢驗。
鑒于成纖維細胞在肺生態系統重塑中的重要性，我們接下來研究了所有主要細胞類型（包括成纖維細胞）的配體 - 受體相互作用（參見 方法）。在所有細胞中豐富的推斷配體 - 受體相互作用中有 TGFβ1-TGFβ 受體 2 和 BMP6-ACVR1（擴展數據圖12g-i，補充表 10），它們分別屬于 TGFβ 家族和超家族。TGFβ 信號傳導在促進肺纖維化方面具有重要作用，并且與成纖維細胞介導的 ADI 維持有關27，這與 DATP 細胞狀態密切相關。為了研究針對 pFB 的潛在治療策略，我們從單核轉錄組推斷蛋白質活性，然后將 pFB 與其他成纖維細胞進行比較。該分析預測 pFB 將顯示 JunB 和 JunD 的活性增加（擴展數據圖12j，補充表 11），它們通過增強的 TGFβ 和 STAT3 信號在小鼠模型中誘導肺纖維化，并與增加的 IL-1β 產生相關30。最后，我們推斷了 pFB 中的可成藥靶標（參見 方法），并將 MMP14 和 STAT3 確定為取消 pFB 中有害程序的潛在靶標（擴展數據圖12j，補充表 11 )。
討論
我們使用短 PMI 尸檢標本和對照肺樣本生成了 COVID-19 的單細胞轉錄組肺圖譜。我們的分析提供了對致命 COVID-19 的細胞景觀、細胞程序和細胞回路的廣泛普查。蛋白質活性和細胞間相互作用的額外推斷，以及使用成像質譜數據分析各種細胞類型的炎性細胞因子，為 SARS-CoV-2 感染肺部的有害后果提供了一個詳細的視角。
我們的分析表明，產生 IL-1β 的異常激活的單核細胞/巨噬細胞之間存在相互作用，炎癥誘導的肺泡上皮再生受損，以及促進纖維化并可能損害再生的病理性成纖維細胞的擴增（擴展數據圖12f，k，補充討論）。除了這些有害事件外，我們的數據表明，盡管可能存在足夠的體液免疫反應（補充討論），但死于 COVID-19 的個體肺部 T 細胞反應不足。最近的一項研究表明，感染 COVID-19 的癌癥患者的 B 細胞功能受損與死亡率增加無關31，但缺乏足夠的 CD8+ T 細胞反應（即使存在足夠的體液免疫）與較差的病毒控制和增加的死亡率有關31。盡管我們的 COVID-19 隊列不包括癌癥患者，但這些數據表明，雖然在針對 SARS-CoV-2 的 T 細胞免疫充足的情況下，體液免疫可能是可有可無的，但我們的患者可能缺乏適當的 T 細胞反應導致了致命的后果。
盡管我們的研究深入了解了宿主對致命 SARS-CoV-2 感染的反應，但它受到樣本量較小的限制。然而，通過協調努力，我們的工作將為一系列研究做出貢獻，包括 TM Delorey 等人的配套論文。32，具有簡化的協議和統一的元數據，以實現集成和組合分析，并將有助于解釋重要的協變量。此外，由于我們的分析側重于死于 COVID-19 的患者的肺組織，因此我們僅檢查了潛在疾病表型的一個子集。盡管如此，一些觀察結果，例如肺纖維化的快速發展（補充討論）)，可能與從嚴重 COVID-19 中幸存下來的患者相關，并可能讓我們了解這些人33的長期并發癥。
總之，我們從短 PMI 組織標本中生成了分子單細胞肺圖譜，并確定了致命的 COVID-19 的病理回路。該圖譜為研究宿主對 SARS-CoV-2 的反應和了解 COVID-19 引起的潛在長期肺部后遺癥建立了重要資源，并為嚴重疾病的治療開發提供了基礎。
