不同分期矽肺患者气道细菌群落微生态研究_荆茹.pdf

实用医学杂志2023年第39卷第9期The Journal of Practical Medicine2023 Vol.39 No.9不同分期矽肺患者气道细菌群落微生态研究荆茹1孟凡利2李旭东1，21广州医科大学公共卫生学院（广州 511436）；2广东省职业病防治院，广东省职业健康工程技术研究中心，广东省职业病防治重点实验室（广州 510300）【摘要】目的探讨不同分期矽肺患者气道细菌微生物构成的差异。方法采用立意抽样方法，选择29例职业性矽肺病患者，13例健康受试者为研究对象，其中矽肺壹期患者6例（矽肺壹期组），矽肺贰期患者15例（矽肺贰期组），矽肺叁期患者8例（矽肺叁期组）。采集42例受试者痰液样本，采用16S rDNA测序技术检测痰液样品中细菌微生物，分析4组受试者细菌微生物的相对丰度差异，对健康受试者与矽肺患者组间存在差异的细菌微生物进行物种多样性指数（alpha多样性指数、beta多样性指数）差异分析；绘制受试者工作特征曲线（ROC）评估细菌微生物作为生物标记物的可行性。结果健康受试者与矽肺患者痰液样品中奈瑟菌属（Neisseria）、莫拉氏菌属（Moraxella）、葡萄球菌属（Staphylococcus）、棒状杆菌属（Corynebacterium）、狡诈球菌属（Dolosigranulum）、密螺旋体菌属（Treponema）相对丰度分别比较，差异均有统计学意义（P 0.05）。LEfSe分析结果显示，矽肺壹期组增加的物种包括Neisseria、奈瑟氏菌科（Neisseriaceae）、假单胞菌目（Pseudomonadales）和莫拉菌科（Moraxellaceae）；矽肺贰期组增加的物种有芽孢杆菌科（Bacillacea）、芽孢杆菌属（Bacillus）、Staphylococcus、葡萄球菌科（Staphyloccaceae）和金黄色葡萄球菌（Staphylococcus_aureus）。ROC曲线分析结果显示，采用Neisseria相对丰度对矽肺壹期与矽肺贰期患者进行分期诊断的ROC曲线下面积为0.87（95%置信区间0.710.99），灵敏度100.0%，特异度73.3%，约登指数0.73，P=0.01。结论Neisseria可能与矽肺患者的疾病进程有关，并可能具备对矽肺分期进行区别的生物标记物的潜力。【关键词】矽肺；细菌；微生物；高通量测序【中图分类号】R135Microecological study of airway bacterial community in patients with silicosis in different disease stagesJING Ru*，MENG Fanli，LI Xudong.School of Public Health，Guangzhou Medical University，Guangzhou 511436，ChinaCorresponding author：LI XudongEmail：【Abstract】ObjectiveTo investigate the composition in bacterial microbial composition in silicosis patients with different stages.Methods13 healthy subjects，29 patients with occupational silicosis，including 6 patients with silicosis I，15 patients with silicosis and 8 patients with silicosis.The sputum of four groups ofsubjects was collected，16S rDNA sequencing technology was used to detect bacterial microorganisms in the sputum samples，and the relative abundance of bacterial microorganisms and the difference of microbial diversity index（alpha diversity index，beta diversity index）were analyzed；the subject working characteristic curve（ROC）was drawn to assess whether bacterial microorganisms have the potential of biomarkers.ResultsThe relative abundance of Neisseria，Moraxella，Staphylococcus，Corynebacterium，Dolosigranulum and Treponema in sputum samples of subjects（P values 18岁。排除标准：罹患其他肺部疾病（如肺部肿瘤、肺结节、囊性纤维化等）或既往进行过肺部手术者；罹患其他炎性疾病者；近 4 周内因其他疾病使用了抗生素或接受过输血治疗者；正在接受全身激素治疗者；正在参加双盲药物临床试验者；长期食用益生菌类食品者。本研究通过广东省职业病防治院医学伦理委员会审核批准；研究对象均签署知情同意书。1.2痰液样品采集采集研究对象入院24 h内的自主咳痰的痰液样品2 mL/人。样本采集时，嘱患者咳深部痰于无菌痰杯中，于1 h内送至-20 的冰箱内保存，1 周内转运至-80 的冰箱内保存待检。1.316S rDNA 测序技术检测痰液中细菌微生物采用十六烷基三甲基溴化铵法（cetyltrimethylammonium bromide，CTAB）从 痰 液 样 品 中 提 取DNA，取2 mL样品DNA用无菌水稀释至质量浓度为1 g/L，以此为模板，采用带有barcode序列的特异性引物（美国Bioo Scientific公司）进行目标区域的聚合酶链式反应（polymerase chain reaction，PCR）（中国诺禾致源公司，Novo G100）扩增。扩增区域为V4区；扩增引物名称：515F，806R；引物序列（方向 5端-3端）为：GTGCCAGCMGCCGCGGTAA，GGACTACHVGGGTWTCTAAT。采用 TruSeqDNA PCRFree 试剂盒（美国Illumina公司）按照说明书步骤生成测序文库22。扩增产物由北京诺禾致源科技股份有限公司基于Illumina NovaSeq6000PE250平台（美国Illumina公司）进行文库测序。1.4生物信息学分析（1）采用 R 语言软件 R.V2.5.3对物种丰度矩阵文本库绘制物种累积箱型图，以衡量测序深度和样品数量是否能有效反映样品中的群落结构情况23。（2）采用 Qiime 软件（V1.9.1）计算 3 组矽肺患者痰液样品中细菌微生物群落Alpha多样性的相关指数，并进行Beta多样性 的 主 坐 标 分 析（principal coordinate analysis，PCoA）。PCoA图展示样品间细菌微生物的总体相似性，为细菌微生物的可视化聚类情况，由第一主坐标和第二主坐标分别解释总变异的程度24-25。Alpha 多样性表示在生物样品中微生物物种在空间和时间上的相对丰度，其描述指标主要包括Shannon指数、Simpon指数、Chao1指数和ACE指数等，可用来描述生态失调26；其中，Simpson指数与1106实用医学杂志2023年第39卷第9期The Journal of Practical Medicine2023 Vol.39 No.9Shannon 指数可衡量微生物群落的多样性；Chao1指数与ACE指数用于衡量微生物群落丰度，指数越大说明群落丰富度越高2729。Beta多样性分析用以衡量细菌群落组成如何随时间或空间而变化，可分析不同群落之间物种多样性的变化趋势30。（3）采用非参数检验的分子方差分析（analysis of molecular variance，AMOVA）进行基于非加权Unifrac距离时3组组间差异检验31。（4）采用线性判别分析（linear discriminant analysis，LDA）效应大小（LDA effect size，LEfSe）方法分析 3 组患者痰液样品中细菌菌属的相对丰度，筛选关键菌属，默认的LDA分数的筛选值为432。1.5统计学分析采用 SPSS 25.0 软件进行统计分析。计量资料经正态性检验不符合正态分布者，采用中位数（M）描述。两组组间 M 比较采用MannWhitney U 检验，多组组间 M 比较采用 KruskalWalis H 检验。采用受试者工作特征（receiveroperator characteristic curve，ROC）曲线评估关键菌属作为生物标记物对矽肺分期进行区别的可行性；当 ROC 曲线下面积（the area under the ROCcurve，AUC）0.50 时，证明该诊断试验具有一定的诊断价值；AUC 值越接近 1.000，证明诊断试验的真实性越好3335。检验水准=0.05（双侧）。2结果2.1各组细菌菌群相对丰度比较健康受试者与3组矽肺患者痰液样品中奈瑟菌属（Neisseria）、莫拉菌属（Moraxella）、葡萄球菌属（Staphylococcus）、棒状杆菌属（Corynebacterium）、狡诈球菌属（Dolosigranulum）、密螺旋体菌属（Treponema）相对丰度分别比较，差异均有统计学意义（P 0.05）。见表1。表14组受试者痰液样品中细菌微生物相对丰度比较结果Tab.1Comparative results of relative abundance of bacterial microorganisms in sputum samples%组别健康对照组矽肺壹期组矽肺贰期组矽肺叁期组H值P值患者数（例）136158Neisseria3.40（0.74，11.69）abc1.66（0.92，1.20）0.41（0.024，0.78）a0.39（0.22，2.62）10.88 0.01Moraxella0.00（0.00，0.01）bc0.79（0.00，6.59）0.47（0.02，0.87）0.01（0.00，0.02）b19.52 0.01Staphylococcus0.22（0.20，3.93）0.68（0.46，2.10）0.94（0.62，1.88）0.21（0.14，0.56）b9.65 0.01Corynebacterium0.15（0.09，0.33）b0.23（0.19，0.34）0.73（0.34，1.64）a0.25（0.16，0.46）15.81 0.01Dolosigranulum0.04（0.02，0.13）b1.04（0.41，2.33）0.42（0.30，0.45）0.29（0.01，0.30）a18.28 0.01Treponema0.59（0.36，1.52）bc0.51（0.13，2.08）0.21（0.11，0.31）0.03（0.01，0.10）ab15.59 0.01注：与矽肺壹期组比较，AP 0.05；与矽肺贰期组比较，BP 0.05；与矽肺叁期组比较，CP 矽肺叁期 矽肺贰期。见表2。2.4Beta 多样性分析基于非加权 Unifrac 距离进行 PCoA，结果显示，第一主坐标差异贡献率为30.27%，第二主坐标差异贡献率为10.81%，该两个主坐标共解释了 41.08%的差异，即 3 组矽肺患者组间群落多样性有差异。见图2。AMOVA分析结1375913111715211923272529样本编号5001 000 1 5002 000物种注释种类图1矽肺患者痰液样品中细菌微生物群落累积箱型图Fig.1Species accumulation boxplot of bacterial microbialcommunities in sputum samples from silicosis patients表23组患者痰液样品细菌微生物组Alpha多样性指数比较结果Tab.2Alpha diversity index of bacterial microbiome ofsputum samples in 3 groups组别矽肺壹期组矽肺贰期组矽肺叁期组患者数（例）6158Shannon指数18.0013.8714.88Simpon指数16.7514.2015.19Chao1指数19.8312.0716.88ACE指数18.6712.2017.501107实用医学杂志2023年第39卷第9期The Journal of Practical Medicine2023 Vol.39 No.9果显示，3组矽肺患者细菌微生物群落物种丰富度比较，差异有统计学意义（离均差平方和=0.55，均方=0.28，F=2.20，P 0.01）；其中，矽肺壹期组与矽肺贰期组患者细菌微生物群落物种丰富度比较，差异有统计学意义（离均差平方和=0.40，均方=0.40，F=3.48，P 0.01），矽肺叁期组患者细菌微生物群落物种丰富度分别与矽肺壹期组、矽肺贰期组比较，差异均无统计学意义（离均差平方和分别为0.29、0.16，均方分别为0.29、0.16，F值分别为1.73、1.40，P值分别为0.08、0.09）。2.5分析矽肺壹期组和矽肺贰期组组间存在的差异菌属基于AMOVA分析，对细菌微生物物种丰富度比较差异有统计学意义的矽肺壹期组和矽肺贰期组进一步进行LEfSe分析。结果显示，矽肺壹期组相较矽肺贰期组增加的物种包括Neisseria、奈瑟氏菌科（Neisseriaceae）、假单胞菌目（Pseudomonadales）和莫拉菌科（Moraxellaceae）；矽肺贰期组相较矽肺壹期组增加的物种有芽孢杆菌科（Bacillaceae）、芽孢杆菌属（Bacillus）、Staphylococcus、葡萄球菌科（Staphylococaceae）、金黄色葡萄球菌（Staphylococcus_aureus）。见图3。2.6绘制差异菌属的ROC曲线结合LEfSe分析结果及 2.2 小节中部分矽肺壹期与贰期细菌微生物相对丰度差异有统计学意义的结果，对Neisseria进行两组矽肺患者的ROC曲线绘制分析。结果显示，Neisseria 相对丰度对矽肺壹期与矽肺贰期进行分期的区别性分析 AUC 值为 0.87（95%置信区间0.710.99，标准误0.08，灵敏度100.0%，特异度73.3%，约登指数 0.73，P=0.01），即 Neisseria 在矽肺壹期组和矽肺贰期组患者痰液样品中具备生物标记物的潜力。见图4。3讨论矽肺发病机制尚未阐明，且由于尚无法通过保守治疗进行根治，可导致过早死亡，对患者及其家庭造成的巨大疾病负担3638。近年来，针对矽肺患者分期鉴别的有关生物标志物研究成为热点，对包括矽肺在内的职业性尘肺病的诊断与鉴别诊断具有一定的价值；其中，矽肺患者痰液样本中的生物标志物也是研究方向之一3940。痰液采集是一种从气道收集细胞的方法，可直接评估气道炎症41。传统的痰培养技术只能描述肺部疾病已有现象，缺乏揭示疾病潜在生物学信息的能力。本研究对矽肺患者的痰液样品进行高通量测序，并对测序文库进行统计学分析后，进一步检验微生物物种多样性的相关指数，如表1、2所示，矽肺壹期、贰期、叁期组患者与健康受试者组间微生物群落存在分布差异。通过检验微生物物种多样性的多个指数，结果表明：与矽肺壹期组比较，矽肺贰期组患者痰液样本中奈瑟菌属（Neisseria）、棒状杆菌属（Corynebacterium）相对丰度存在差异，矽肺叁期组患者痰液样本中狡诈球菌属（Dolosigranulum）、密螺旋体菌属（Treponema）相对丰度存在差异；与矽肺贰期组患者比较，矽肺叁期组患者痰液-1.0-0.50.00.5第一主坐标0.20.00.20.4第二主坐标矽肺壹期矽肺贰期矽肺叁期图2细菌微生物基于非加权Unifrac距离PCoA图Fig.2Bacterial microbes based on the unweighted Unifracdistance PCoA plot-4-2024LDAg_Neisseriaf_Neisseriaceaeo_Pseudomonadalesf_Moraxellaceaef_Bacillaceaeg_Bacillusg_Staphylococcusf_Staphylococcaceaes_Staphylococcus_aureus矽肺贰期矽肺壹期图3两组矽肺患者LEfSe分析图Fig.3LEfSe analysis of 2 groups图4矽肺壹期6例与矽肺贰期15例样本中的Neisseria物种相对丰度ROC曲线Fig.4ROC curves of relative abundance of Neisseria in 6silicosis and 15 silicosis 0.00.20.40.60.81.01-特异度AUC值：0.867（0.712-1.000）0.00.20.40.60.81.0敏感度ROC曲线Neisseria参考线1108实用医学杂志2023年第39卷第9期The Journal of Practical Medicine2023 Vol.39 No.9样 本 中 莫 拉 氏 菌 属（Moraxella）、葡 萄 球 菌 属（Staphylococcus）、密螺旋体菌属（Treponema）相对丰度存在差异。对存在分布差异的细菌微生物采用LEfSe分析后，根据差异菌绘制ROC曲线，结果显示，矽肺壹期组患者痰液样品中Neisseria相对丰度高于矽肺贰期组（P 0.05）。Neisseria是一群革兰染色阴性双球菌，菌属包括脑膜炎奈瑟菌（N.meningitidis）、淋病奈瑟菌（N.gonorrhoeae）、干燥奈瑟菌（N.sicca）、微黄奈瑟菌（N.subflava）、浅黄奈瑟菌（N.flavescens）和黏液奈瑟菌（N.mucosa）等。其中，可导致人体发病的只有N.meningitidis 和N.gonorrhoeae；非致病性的N.sicca、N.subflava、N.flavescens 和 N.mucosa 在实际临床工作中很少被考虑，尽管有证据表明它们可引起侵袭性机会性感染42。STANCIL 等43-44研究结果显示，特发性肺纤维化患者和对照受试者中Neisseria、链球菌属和嗜血杆菌属分布存在差异。KITSIOS等45研究结果显示，在预后不良的的急性期慢性阻塞性肺气肿患者痰液样品中Neisseria 相对丰度较阴性对照组降低。非致病菌的细菌相对丰度增高，导致气道微生态发生变化，人类宿主与共生菌之间的平衡随之改变，而患者和非患者间的差异菌属，如非致病性Neisseria 在免疫疗法方面可能能够显示出潜力46。目前有研究尝试通过微生物对矽肺患者进行治疗，可能成为辨别矽肺患者疾病进展及预后治疗的一个未来研究方向。痰液样本简便易得，其特征化有助于辅助临床医生对矽肺患者疾病发生发展程度的辨别，建议将痰液微生物分析加入诊断过程中，可以帮助医师认识到患者可能面临纤维化进一步加重的风险，进而加强用药。本研究首次对不同分期矽肺患者气道细菌微生物存在的差异进行了研究，但受限于矽肺患者产生的痰液量较少，导致组内异质性较大；所采集痰液样本为自主咳痰样本，存在杂菌污染的可能，从而导致缺乏对于结论的验证完整性，未来可能需要通过进一步补充样本量及采用其他能够减少杂菌的痰液采集方式进行后续研究。综上所述，矽肺壹期、贰期、叁期患者痰液中的细菌菌属构成存在差异，特征细菌菌属 Neisseria可能与矽肺患者的疾病进程有关，尚有待进一步加大基础研究和临床验证，以促进其临床应用。（志谢：感谢对本项研究提供帮助的樊春月医生、刘海月老师。）【Author contributions】JING Ru performed the analysis of dataand wrote the article.MENG Fanli performed the experiments.LIXudong designed the study and reviewed the article.All authors readand approved the final manuscript as submitted.参考文献1HOY R F，CHAMBERS D C.Silicarelated diseases in the modern world J.Allergy，2020，75（11）：28052817.2TAN S，CHEN S.The mechanism and effect of autophagy，apoptosis，and pyroptosis on the progression of silicosis J.Int J MolSci，2021，22（15）：8110.3KREFFT S，WOLFF J，ROSE C.Silicosis：An update andguide for clinicians J.Clin Chest Med，2020，41（4）：709722.4COLLABORATORS G C O D.Global，regional，and national agesexspecific mortality for 282 causes of death in 195 countriesand territories，19802017：A systematic analysis for the globalburden of disease study 2017 J.Lancet，2018，392（10159）：17361788.5TAN S，CHEN S.Macrophage autophagy and silicosis：Currentperspective and latest insights J.Int J Mol Sci，2021，22（1）：453.6中华人民共和国国家卫生健康委员会.2021年我国卫生健康事业发展统计公报 EB/OL 20220712.http：/ J.中国职业医学，2022，49（3）：280286.8中华人民共和国国家卫生和计划生育委员会.职业性尘肺病的诊断：GBZ 702015S.北京：中国标准出版社，2015.9HOY R，CHAMBERS D C.Silicosis：An ancient disease inneed of a dose of modern medicineJ.Respirology，2020，25（5）：464465.10 PINKERTON J W，KIM R Y，ROBERTSON A A B，et al.Inflammasomes in the lung J.Mol Immunol，2017，86：4455.11 ADAMCAKOVA J，MOKRA D.New insights into pathomechanisms and treatment possibilities for lung silicosis J.Int J MolSci，2021，22（8）：4162.12 乔洁洁，张璐，韩科，等.职业性矽肺患者肺泡巨噬细胞中铁蛋白轻链蛋白表达分析 J.中国职业医学，2021，48（5）：519522.13 张玉，耿晓，赛林霖，等.自噬在矽肺模型大鼠肺部炎症反应中的作用 J.中国职业医学，2022，49（2）：153157，164.14 INSERRA A，ROGERS G B，LICINIO J，et al.The microbiotainflammasome hypothesis of major depressionJ.Bioessays，2018，40（9）：e1800027.15 FLUME P A，CHALMERS J D，OLIVIER K N.Advances inbronchiectasis：Endotyping，genetics，microbiome，and diseaseheterogeneity J.Lancet，2018，392（10150）：880890.16 GALIANA A，AGUIRRE E，RODRIGUEZ J C，et al.Sputummicrobiota in moderate versus severe patients with copd J.EurRespir J，2014，43（6）：17871790.17 GOLEVA E，JACKSON L P，HARRIS J K，et al.The effects ofairway microbiome on corticosteroid responsiveness in asthmaJ.Am J Respir Crit Care Med，2013，188（10）：11931201.18 CARMODY L A，ZHAO J，KALIKIN L M，et al.The daily dynamics of cystic fibrosis airway microbiota during clinical stability and at exacerbation J.Microbiome，2015，3：12.19 MAUGHAN H，CUNNINGHAM K S，WANG P W，et al.Pul1109实用医学杂志2023年第39卷第9期The Journal of Practical Medicine2023 Vol.39 No.9monary bacterial communities in surgically resected noncystic fibrosis bronchiectasis lungs are similar to those in cystic fibrosisJ.Pulm Med，2012，2012：746358.20 LEVY S E，MYERS R M.Advancements in nextgeneration sequencing J.Annu Rev Genomics Hum Genet，2016，17：95115.21 YATERA K，NOGUCHI S，MUKAE H.The microbiome in thelower respiratory tract J.Respir Investig，2018，56（6）：432439.22 BOWMAN S K，SIMON M D，DEATON A M，et al.Multiplexed illumina sequencing libraries from picogram quantities ofDNA J.BMC Genomics，2013，14：466.23 GREMBI J A，NGUYEN L H，HAGGERTY T D，et al.Gut microbiota plasticity is correlated with sustained weight loss on alowcarb or lowfat dietary intervention J.Sci Rep，2020，10（1）：1405.24 CAPORASO J G，KUCZYNSKI J，STOMBAUGH J，et al.Qiime allows analysis of highthroughput community sequencing data J.Nat Methods，2010，7（5）：335336.25 BECKERS B，OP DE BEECK M，WEYENS N，et al.Structuralvariability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of fieldgrown poplar trees J.Microbiome，2017，5（1）：25.26 AVALOSFERNANDEZ M，ALIN T，MTAYER C，et al.Therespiratory microbiota alphadiversity in chronic lung diseases：First systematic review and meta analysisJ.Respir Res，2022，23（1）：214.27 EINARSSON G G，COMER D M，MCILREAVEY L，et al.Community dynamics and the lower airway microbiota in stable chronic obstructive pulmonary disease，smokers and healthy non smokers J.Thorax，2016，71（9）：795803.28 SVERRILD A，KIILERICH P，BREJNROD A，et al.Eosinophilic airway inflammation in asthmatic patients is associatedwith an altered airway microbiome J.J Allergy Clin Immunol，2017，140（2）：407417.29 FEIGELMAN R，KAHLERT C R，BATY F，et al.Sputum DNAsequencing in cystic fibrosis：Noninvasive access to the lung microbiome and to pathogen details J.Microbiome，2017，5（1）：20.30 MA Y，YANG X，HUA G，et al.Contribution of gut microbiomes and their metabolomes to the performance of dorper and tansheep J.Front Microbiol，2022，13：1047744.31 HUANG K，WANG T，DUNN D W，et al.A generalized framework for amova with multiple hierarchies and ploidies J.IntegrZool，2021，16（1）：3352.32 CHIANG T Y，YANG Y R，ZHUO M Y，et al.Microbiome profiling of nasal extracellular vesicles in patients with allergic rhinitis J.World Allergy Organ J，2022，15（8）：100674.33 方积乾.生物医学研究的统计方法 M.2版.北京：高等教育出版社，2020：193196.34 NAHM F S.Receiver operating characteristic curve：Overviewand practical use for clinicians J.Korean J Anesthesiol，2022，75（1）：2536.35 MUSCHELLI J.Roc and auc with a binary predictor：A potentially misleading metric J.J Classif，2020，37（3）：696708.36 ALTUNDA HATMAN E，ACAR KARAGL D，KUMAN OYMAN E，et al.Premature deaths due to silicosis in turkey，20062017：A twelve year longitudinal studyJ.Balkan Med J，2021，38（6）：374381.37 WANG F X，GAO F Y，LIU X，et al.Long noncoding rna expression in silicosis and mrak050699 function in epithelialmesenchymal transition J.Hum Exp Toxicol，2021，40（12）：763774.38 DESCATHA A，SEMBAJWE G，BAER M，et al.Who/ilo workrelated burden of disease and injury：Protocol for systematic reviews of exposure to long working hours and of the effect of exposure to long working hours on strokeJ.Environ Int，2018，119：366378.39 董文雅，魏文静，谭卫国，等.常规抗结核化疗对肺结核患者痰液菌群的影响 J.实用医学杂志，2021，37（3）：336342.40 蔡文臣，张诗汇，魏中秋，等.Acsdkp抑制spp1/tgf1信号拮抗实验性矽肺的机制 J.实用医学杂志，2021，37（2）：159163.41 DICKER A J，HUANG J T J，LONERGAN M，et al.The sputum microbiome，airway inflammation，and mortality in chronicobstructive pulmonary diseaseJ.J Allergy Clin Immunol，2021，147（1）：158167.42 HUMBERT M V，CHRISTODOULIDES M.Atypical，yet not infrequent，infections with neisseria speciesJ.Pathogens，2019，9（1）：10.43 STANCIL I T，MICHALSKI J E，DAVISHALL D，et al.Pulmonary fibrosis distal airway epithelia are dynamically and structurally dysfunctional J.Nat Commun，2021，12（1）：4566.44 NATALINI J G，SINGH S，SEGAL L N.The dynamic lung microbiome in health and disease J.Nat Rev Microbiol，2022，10（16）：114.45 KITSIOS G D，FITCH A，MANATAKIS D V，et al.Respiratorymicrobiome profiling for etiologic diagnosis of pneumonia in mechanically ventilated patientsJ.Front Microbiol，2018，9：1413.46 DOREY R B，THEODOSIOU A A，READ R C，et al.The nonpathogenic commensal neisseria：Friends and foes in infectiousdise