a, IRRs (95% CIs) количества осложнений для госпитализированных (
n = 114.864) и негоспитализированных (
n = 20.297) групп COVID-19 по годам. IRRs были оценены в сравнении с контрольной группой без инфекции (
n = 5.206.835). Точки в центре полос погрешностей представляют скорректированные IRR, оцененные с использованием числа постострых осложнений, а полосы погрешностей соответствуют 95% CIs.
b, Количество осложнений на 1.000 человек от инфекции SARS-CoV-2 в госпитализированных и негоспитализированных группах COVID-19 по годам.
c, DALYs на 1.000 человек вследствие инфекции SARS-CoV-2 в госпитализированных и негоспитализированных группах COVID-19 по годам.
d, Кумулятивное число осложнений на 1.000 человек в связи с инфекцией SARS-CoV-2 в группах госпитализированных и негоспитализированных COVID-19 в течение 3 лет.
e, Кумулятивные DALYs на 1.000 человек вследствие инфекции SARS-CoV-2 для госпитализированных и негоспитализированных групп COVID-19 за 3 года.
a заказывается IRRs среди госпитализированных групп COVID-19 каждый год. Красными и синими точками показаны IRR в госпитализированных и негоспитализированных группах COVID-19 значительно больше 1, а полосы погрешностей связаны с ними CIs.
B
b,
c,
d и
e, красными точками обозначены абсолютные показатели в группах госпитализированных COVID-19, значительно выше, чем в контрольной группе без инфекции; синие точки представляют абсолютные показатели в негоспитализированных группах COVID-19, значительно выше, чем в контрольной группе без инфекции; а серые точки представляют абсолютные показатели в госпитализированных или негоспитализированных группах COVID-19, которые не были значительно выше, чем в контрольной группе без инфекции. Более толстые горизонтальные столбцы представляют избыточные показатели в группе госпитализированных с COVID-19 по сравнению с группой без госпитализации с COVID-19, где красные столбцы указывают на значительно отличающиеся показатели, а серые столбцы указывают на отсутствие статистической разницы. Более тонкие горизонтальные полосы, расположенные ближе к оси y, представляют абсолютный показатель в группе без госпитализации COVID-19 по сравнению с контрольной группой без инфекции. Системы органов в
группах b,
c,
d и
e были отсортированы на основе статистической значимости и величины различий между госпитализированными и негоспитализированными группами COVID-19 (горизонтальные красно-серые полосы между двумя точками).
Анализ чувствительностиМы провели несколько анализов чувствительности.
(1) Мы построили двойную устойчивую модель корректировки, в которой ковариаты использовались как в моделях воздействия, так и в моделях результатов, вместо основного подхода, где ковариаты применялись только в модели воздействия.
(2) В первичном подходе мы построили модели Пуассона с нулевым раздуванием вместо моделей Poisson (Пуассона).
(3) Мы не подвергали цензуре участников групп с COVID-19 при повторном заражении, в то время как в первичном подходе участники групп с COVID-19 подвергались цензуре при повторном заражении.
(4) Мы дополнительно скорректировали 100 алгоритмически выбранных многомерных ковариат, вместо того чтобы использовать только набор предварительно определенных ковариат в нашем первичном подходе.
(5) Вместо того, чтобы определять госпитализацию во время острой фазы как дату госпитализации в течение 7 дней до или в течение 30 дней после положительного теста в основных анализах, мы использовали альтернативное определение госпитализации как дату госпитализации в день положительного теста или в течение 30 дней после положительного теста.
(6) Мы усекли веса оценки склонности на уровне 99,5% процентилей, а не 99,9% процентилей в основных анализах.
(7) Мы оценили IRRs среди выборки с полными данными по всем ковариатам (
n = 4.432.414, 83,0% от полной выборки) для проверки согласованности результатов с результатами, полученными с помощью множественного импутации отсутствующих данных.
(8) Мы оценили риски на основе моделей Файн-Грея, в которых смерть и инфекция SARS-CoV-2 во время наблюдения рассматривались как конкурирующие риски.
(9) Мы применили обратную вероятность цензурирования для учета неслучайной цензуры в связи со смертью или инфекцией SARS-CoV-2 во время наблюдения в трех группах (контрольная группа без инфекции, негоспитализированная группа COVID-19 и группа госпитализированных COVID-19).
(10) В качестве альтернативы мы использовали более узкое определение PASC, которое включало 73 исхода вместо 80 исходов, включенных в первичный анализ. Результаты этих анализов чувствительности согласуются с результатами основных анализов (Дополнительная таблица
7).
Анализ контроля негативных исходовМы изучили связь между COVID-19 и заболевшими новообразованиями в качестве контроля негативного исхода. Результаты показали нейтральную связь между COVID-19 и контролем негативного исхода у негоспитализированных участников (IRRs 1,03 (0,98–1,08) в первый год, 0,94 (0,87–1,01) во второй год и 0,95 (0,88–1,03) в третий год) и госпитализированных пациентов (IRRs 0,93 (0,82–1,05) в первый год, 0,92 (0,72–1,12) во второй год и 0,93 (0,68–1,08) в третий год).
ОбсуждениеВ этом исследовании 135.161 человека с инфекцией SARS-CoV-2 и 5.206.835 человек в контрольной группе, за которыми наблюдали в течение 3 лет, мы показываем, что среди негоспитализированных лиц повышенный риск смерти больше не наблюдался после первого года инфекции, а DALYs и риск постострых осложнений значительно снизились за 3 года: 41,2 осложнения на 1.000 человек и 9,6 DALYs на 1.000 человек, начиная с третьего года. Среди госпитализированных пациентов риск смерти снизился в течение 3 лет, но оставался значительно повышенным на третий год после инфицирования (на 29% повышенный риск и избыточное бремя смертности – 8,16 на 1000 человек). Риск развития постострых осложнений с годами снижался, но существенный остаточный риск сохранялся и на третий год, что привело к 252,8 случаям на 1000 человек и 90,0 DALYs на 1000 человек. В целом, наши результаты показывают снижение рисков в течение 3 лет наблюдения, но стойко повышенный риск серьезных неблагоприятных исходов среди госпитализированных лиц.
Это 3-летнее последующее исследование дополняет результаты нашего предыдущего 2-летнего последующего анализа [20]. За 2 года PASC обеспечил около 81 и 676 DALYs на 1000 человек у негоспитализированных и госпитализированных лиц соответственно. На третий год наблюдения PASC дополнительно показал 9,6 и 90,0 DALYs на 1000 человек у негоспитализированных и госпитализированных лиц соответственно. Эти результаты согласуются с единственным 3-летним исследованием PASC, которое было опубликовано на сегодняшний день, показывающим стойкую симптоматику в этих системах органов [21] и вкладом Taquet et al. [15], демонстрируя изменение горизонта риска нейропсихиатрических осложнений через 2 года. Представленные здесь эпидемиологические данные о стойкости риска также согласуются с исследованиями Peluso et al. [22], в которых они показывают, что инфекция SARS-CoV-2 может приводить к стойкой активации Т-клеток в различных тканях организма, которая может все еще проявляться через 2–3 года после первоначальной легкой инфекции и может быть связана с PASC, предполагая, что даже отдаленная и клинически легкая инфекция SARS-CoV-2 может иметь долгосрочные последствия для тканевого иммунного гомеостаза.
Уроки, извлеченные из постострых инфекционных синдромов, позволяют сделать два ключевых наблюдения:
(1) риски для многих состояний в постострой фазе со временем снижаются (т.е. риск становится равным в инфицированной и контрольной группах через несколько месяцев или несколько лет после первоначального заражения) и
(2) Могут существовать латентные эффекты, примером которых являются вирус Эпштейна-Барр, приводящий к рассеянному склерозу, или вирус полиомиелита, приводящий к постполиомиелитному синдрому через много лет или десятилетий после первоначальной инфекции [23, 24, 25].
Поскольку SARS-CoV-2 является новым вирусом (который был обнаружен только в конце 2019 года), последующее наблюдение за инфицированными людьми доступно только в течение нескольких лет. При 3-летнем наблюдении в нашей когорте не было доказательств в отношении новых неблагоприятных исходов, которые ранее не проявлялись. Долгосрочные исследования с тщательным наблюдением за инфицированными лицами для оценки траекторий риска имеют решающее значение для лучшего понимания горизонта риска последствий для здоровья и выявления предполагаемых скрытых эффектов [15].
Предыдущая эпидемиологическая характеристика PASC сообщала о повышенном риске в зависимости от тяжести инфекции
; наш анализ также добавляет, что риск не только выше среди госпитализированных по сравнению с негоспитализированными группами в каждый момент времени, но и что риск сохраняется дольше среди госпитализированных пациентов, усугубляя их кумулятивное бремя потери здоровья. В то время как 3-летнее кумулятивное бремя DALYs из-за PASC было высоким среди негоспитализированной группы (91,2 DALYs), оно было почти в 8,4 раза выше среди госпитализированной группы (766,2 DALYs), что отражает значительные потери здоровья в этой группе. Механизмы более длительного сохранения риска у людей с более тяжелым острым заболеванием не совсем ясны. Объяснение может быть связано, в частности, с уязвимостью людей, у которых развивается тяжелая форма COVID-19, в отношении более сопутствующих заболеваний, дисфункции иммунной системы или генетической предрасположенности [27, 28, 29, 30, 31].
Исследования показали, что при тяжелом течении COVID-19 SARS-CoV-2 приводит к системной инфекции с репликацией вируса во внелегочных участках и персистенцией RNA в различных тканях, включая системы органов, которые демонстрируют стойкие риски в наших анализах [32, 33] . Неясно, играет ли и в какой степени повреждение во время острой фазы и/или персистенция вируса механистическую роль в пролонгированном горизонте риска в этих органах. Как более высокий, так и более длительный горизонт риска постострых осложнений, а также устойчиво повышенный риск смерти через 3 года среди госпитализированных с COVID-19 по сравнению с группой без госпитализации COVID-19 позволяют предположить, что тяжесть острой инфекции является ключевым фактором выражения долгосрочных неблагоприятных исходов для здоровья. Снижение риска госпитализации (за счет более широкого охвата вакцинацией и противовирусными препаратами), профилактика развития осложнений у госпитализированных лиц, а также раннее распознавание и лечение PASC у пострадавших лиц являются ключевыми стратегиями снижения риска хронической потери здоровья у людей с инфекцией SARS-CoV-2.
Несмотря на то, что мы подчеркиваем высокую потерю здоровья из-за PASC у людей с тяжелым течением COVID-19, и хотя риск PASC (в относительной шкале) меньше у людей с легкой формой COVID-19, их абсолютное число намного больше (чем у госпитализированных) из-за гораздо более высокой доли людей с легким течением
заболевания [26]. Следовательно, большая часть бремени PASC в популяциях связана с легкой инфекцией.
Согласно анализу, проведенному сотрудниками Глобального бремени болезней (GBD – Global Burden of Disease), около 90% людей с PASC протекали в легкой форме COVID-19, что позволяет предположить, что, хотя профилактика тяжелого заболевания имеет важное значение, также необходимы стратегии по снижению риска постострой и долгосрочной потери здоровья у людей с легкой формой COVID-19 [1, 26].
Настоящее исследование имеет несколько сильных сторон. Он использовал широту и глубину общенациональной интегрированной системы здравоохранения в Соединенных Штатах и собрал большую когорту из 5.341.996 участников и наблюдал за ними в течение 3 лет. Мы использовали передовые статистические подходы для балансировки исходных характеристик; Мы оценили риск развития широкого спектра осложнений в 10 системах органов у людей, которые не были госпитализированы и были госпитализированы во время острой фазы заболевания; Мы также предоставили оценки риска как в относительной, так и в абсолютной шкале, оценив количество осложнений в каждом году и в совокупности в 3 года. Кроме того, мы оценили DALYs, которые обеспечивают более полную оценку нагрузки, учитывающую возникновение осложнений и их влияние на общее состояние здоровья. Мы провели несколько анализов чувствительности, которые дали стабильные результаты; и тестирование контроля негативных исходов дало результаты, соответствующие априорным ожиданиям.
Это исследование также имеет несколько ограничений. Популяция VA в основном пожилая, белая и мужская, что не может быть обобщено на другие популяции, не являющиеся ветеранами. Чтобы убедиться, что у участников есть данные наблюдения за 3 годами, мы включили людей, у которых были положительные тесты на SARS-CoV-2 в 2020 году – эпохе, которая предшествовала введению вакцин против COVID-19. Следовательно, полученные здесь результаты не отражают долгосрочных последствий для здоровья инфицированных SARS-CoV-2 среди вакцинированных лиц, которые, учитывая эффективность вакцин в снижении риска PASC, могут быть ниже. Риски, представленные в этих результатах, также могут не отражать риски людей с противовирусными препаратами или тех, кто перенес инфекцию SARS-CoV-2 более поздними вариантами (например, Omicron). Несмотря на то, что мы использовали большое количество предварительно определенных ковариат из нескольких источников данных, включая диагнозы, лекарства и результаты лабораторных исследований, и сбалансировали эти ковариаты между группами COVID-19 и контрольной группой без инфекции, мы не можем полностью исключить смещение классификации и остаточное смешение. Мы заранее определили всеобъемлющий набор исходов, о которых сообщалось как о последствиях инфекции SARS-CoV-2; последствия, которые еще предстоит охарактеризовать, не включены в этот анализ. Мы использовали всеобъемлющее определение PASC, охватывающее 80 последствий; Более узкие или более широкие определения могут привести к различным оценкам. Оценки бремени в абсолютном масштабе также отражают влияние исходных показателей. Поскольку количество участников варьируется между негоспитализированной и госпитализированной группами, 95% CIs следует интерпретировать в контексте размера выборки в каждой группе; Аналогичным образом, оценки рисков и бремени в 1, 2 и 3 годах и в совокупности через 3 года следует интерпретировать вместе с их интервалами неопределенности. Мы оценили риски и бремя в 1, 2 и 3 годы после заражения SARS-CoV-2, и наши оценки отражают средний риск за эти периоды времени. Несмотря на то, что ресурсы данных VA COVID-19 всесторонне собирают результаты тестов на COVID-19 из широкого спектра источников данных (включая данные VA и другие данные), мы не можем исключить возможность недиагностированного COVID-19 или положительных результатов на COVID-19, которые не были зарегистрированы системами данных VA. Если это произошло в контрольной группе, это могло снизить расчетное бремя PASC. Мы использовали данные и методологии исследования GBD для присвоения коэффициента бремени здоровья для каждого исхода здоровья [26]. Этот подход предполагает, что состояние здоровья имеет одинаковый коэффициент нагрузки на здоровье независимо от его причины (например, COVID-19 или другие факторы).
В целом, наши результаты свидетельствуют о значительном снижении рисков PASC в течение 3-летнего периода наблюдения. Среди негоспитализированных лиц сохраняется небольшой остаточный риск, а среди госпитализированных лиц наблюдается нетривиальное бремя смерти и потери здоровья.
МетодыЗаявление об этикеЭто исследование было одобрено институциональным наблюдательным советом системы здравоохранения VA St. Louis, который предоставил отказ от информированного согласия (протокол No 1606333).
Дизайн исследования и обстановкаДанное исследование проводилось с использованием электронных баз данных здравоохранения Министерства по делам ветеранов США. VA управляет крупнейшей общенациональной интегрированной системой здравоохранения в Соединенных Штатах, включающей 1321 медицинское учреждение (172 медицинских центра и 1138 амбулаторных пунктов), обслуживающих более 9 миллионов ветеранов США каждый год. Ветераны, зарегистрированные в системе здравоохранения VA, имеют доступ к широкому спектру медицинских услуг, включая амбулаторную помощь, стационарную помощь, рецепты, психиатрическую помощь, медицинское обслуживание на дому, первичную медико-санитарную помощь, специализированную помощь, гериатрический и расширенный уход, реабилитационные услуги, медицинское оборудование и протезирование.
Источники данныхЭлектронные медицинские базы данных Управления по делам ветеранов содержат исчерпывающую информацию об амбулаторных и стационарных приемах, результатах лабораторных исследований и лекарственных препаратах во время плановых медицинских осмотров и обновляются ежедневно
дней [25]. Источники данных также включали систему поиска идентификационных записей бенефициаров VA, файл о жизненно важном состоянии Medicare, мастер-файл Администрации социального обеспечения и Национальную администрацию кладбищ (VA Beneficiary Identification Record Locator System, the Medicare Vital Status File, the Social Security Administration Master File and the National Cemetery Administration.). Для получения информации о диагнозах и процедурах мы использовали стационарные и амбулаторные домены баз данных VA Corporate Data Warehouse [34, 35, 36]. Для получения данных об аптечных записях использовались домены амбулаторной аптеки и администрирования лекарств со штрих-кодом. Для получения данных по лабораторным измерениям использовалась область лабораторных результатов. Результаты тестов на инфекцию SARS-CoV-2 были получены из общего ресурса данных VA COVID-19, который состоял из результатов тестов полимеразной цепной реакции, тестов на антигены, проведенных в VA, или тестов, сообщенных в VA [37]. Данные стационарных и амбулаторных пациентов Medicare были получены от
Центров VA по услугам Medicare и Medicaid (CMS – Centers for Medicare and Medicaid Services). Индекс депривации территории (ADI – area deprivation index) использовался в качестве суммарного показателя контекстуального социально-экономического неблагополучия (доход, образование, занятость и качество жилья) [38].
КогортаБлок-схема, показывающая построение когорты, представлена на рис.
5. Впервые мы идентифицировали группу контакта с первой инфекцией SARS-CoV-2 в период с 1 марта по 31 декабря 2020 года (
n = 149.459). Затем мы отобрали тех людей, которые являются пользователями VA, определяемых как имевшие по крайней мере два обращения за медицинской помощью с интервалом не менее 180 дней в течение 2 лет до инфекции (
n = 143.034). Чтобы изучить риск PASC, мы отобрали людей, которые были живы через 30 дней после заражения, в результате чего была сформирована аналитическая когорта из 135.161 человека в группе COVID-19. Госпитализация в острой фазе определялась как дата госпитализации в течение 7 дней до или в течение 30 дней после положительного теста. Затем группа COVID-19 была дополнительно разделена по учреждениям медицинской помощи во время острой фазы инфекции на негоспитализированные (
n = 114.864) и госпитализированные (
n = 20 297) группы COVID-19. Дата положительного результата теста на SARS-CoV-2 была определена как
T0, а последующее наблюдение началось через 30 дней после
T0. За участниками наблюдали до смерти, повторной инфекции SARS-CoV-2, через 1.080 дней после первого заражения или до 31 декабря 2023 года.
Чтобы составить контрольную группу без инфекции SARS-CoV-2, мы сначала выявили 6.231.638 человек, которые были живы на 1 марта 2020 года и не имели положительного результата теста на SARS-CoV-2 в период с 1 марта 2020 года по 30 марта 2021 года. Затем мы случайным образом распределили
T0 для контрольной группы на основе распределения
T0 в общей группе COVID-19, чтобы убедиться, что доля участников, начавших лечение в определенную дату, была одинаковой между группой COVID-19 и контрольной группой; 6 194 973 участника были живы в случайно распределенной группе.
T0. Как и в группе с COVID-19, мы также потребовали, чтобы контрольная группа сталкивалась с системой здравоохранения VA не менее двух раз с интервалом не менее 180 дней в течение 2 лет до назначенного
T0, в результате чего в контрольной группе без инфекции была получена итоговая аналитическая когорта из 5.206.835 участников. Наблюдение началось через 30 дней после
T0, и участники наблюдались до смерти от инфекции SARS-CoV-2, через 1.080 дней после
T0 или 31 декабря 2023 года.
РезультатыМы предварительно определили список из 80 отдельных исходов, которые являются хорошо охарактеризованными последствиями инфекции SARS-CoV-2 на основе предыдущих доказательств [1, 4, 5, 6, 8, 9, 11, 12, 13, 20, 25, 27]. Эти исходы были определены с использованием кодов диагноза), медицинских назначений и лабораторных измерений [20, 25]. Исходы инцидентов были идентифицированы, когда они не происходили в течение 2 лет до
T0 и были первыми случаями с 30 дней после
T0 до конца наблюдения. Затем отдельные исходы были сгруппированы в 10 систем органов: сердечно-сосудистую, коагуляционную и гематологическую, утомляющую, желудочно-кишечную, почечную, психическую, метаболическую, скелетно-мышечную, неврологические и легочные. Для общего PASC и исходов на уровне системы органов мы оценили количество осложнений как сумму возникновения отдельных исходов, включенных в составной исход. Кроме того, мы использовали методологии GBD для оценки DALYs, которые представляют собой меру бремени болезни, учитывающую количество включенных осложнений и их влияние на общее состояние здоровья [20, 26, 39]. Для каждого отдельного исхода был присвоен коэффициент нагрузки на здоровье
[20, 25, 26]. Затем DALYs оценивались по взвешенной сумме всех отдельных исходов, включенных в составной исход, где вес – коэффициент нагрузки на здоровье для каждого отдельного исхода [20, 25].
КовариатыНабор предварительно определенных ковариат был выбран на основе предварительных знаний о потенциальных искажающих факторах, которые могут повлиять на взаимосвязь между инфекцией SARS-CoV-2 и PASC [1, 4, 6, 7, 8, 12, 13, 20, 25]. Демографические ковариаты включали возраст, пол, расовую принадлежность (белый, черный и другие), ADI по адресу проживания и статус курильщика (никогда, бывшие и нынешние курильщики). Дополнительные ковариаты включали
расчетную скорость клубочковой фильтрации (eGFR – estimated glomerular filtration rate), систолическое и диастолическое артериальное давление и индекс массы тела, измеренные до и ближе всего к
T0. Набор переменных, определяющих использование медицинских услуг, включал использование долгосрочного ухода в течение года до пандемии, получение вакцинации против сезонного гриппа каждый год в течение 5 лет до пандемии
T0, количество стационарных и амбулаторных посещений в рамках программы Medicare за год до пандемии, количество стационарных и амбулаторных уникальных медицинских рецептов, а также количество стационарных и амбулаторных лабораторных панелей в медицинской системе VA, разделенных интервалами в 180 дней. Сопутствующие заболевания включали тревогу, сердечно-сосудистые заболевания, цереброваскулярные заболевания, хронические заболевания почек, хронические заболевания легких, деменцию, депрессию, диабет, ослабленный иммунитет (трансплантация органов в анамнезе, терминальная стадия заболевания почек, рак, HIV (ВИЧ) или назначение кортикостероидов или иммунодепрессантов) и заболевания периферических артерий. Для учета пространственно-временных вариаций учитывалась календарная неделя заражения SARS-CoV-2 для групп COVID-19 или присвоенный
T0 для контрольной группы, а также географическое расположение медицинской службы. Пропущенные значения включают 9,2% для eGFR, 4,9% для систолического и диастолического артериального давления и 10,4% для индекса массы тела. Отсутствующие данные были условно исчислены с помощью многомерного условного расчета с помощью цепных уравнений и метода сопоставления с прогностическим средним при условии всех ковариат в группах COVID-19 и контрольной группе по отдельности [25]. Все ковариаты были измерены с использованием периода ретроспективного анализа в 2 года до
T0, если не указано иное.
Статистический анализГруппа COVID-19 была разделена по учреждениям медицинской помощи во время острой фазы на две взаимоисключающие группы: негоспитализированные и госпитализированные группы COVID-19. Сообщалось о исходных характеристиках групп COVID-19 и контрольной группы без инфекции. Непрерывные переменные были представлены как средние (стандартные отклонения), а категориальные переменные – кaк частоты (проценты). Стандартизированные различия средних были рассчитаны для оценки ковариатного баланса между группами COVID-19 и контрольной группой без инфекции, где значение менее 0,1 считалось свидетельством хорошего ковариатного баланса. Аналитическая блок-схема представлена на рис.
6.
Для балансировки исходных различий между двумя группами COVID-19 и контрольной группой без инфекции было использовано обратное вероятностное взвешивание [20]. Были построены логистические регрессионные модели для оценки вероятности отнесения к целевой группе с учетом заранее определенных ковариат (оценка склонности). Чтобы обеспечить репрезентативную оценку риска, мы выбрали общую популяцию (группы COVID-19 и контрольную группу) в качестве целевой популяции. Затем были рассчитаны обратные веса вероятностей для всех трех групп как оценка склонности, деленная на (1 − оценка склонности). Мы усекли веса оценки предрасположенности на уровне 99,9% процентилей в каждой группе (контрольная группа, группа без госпитализации с COVID-19 и группа с госпитализированным COVID-19), чтобы уменьшить влияние чрезмерно больших весов на аналитические результаты. Мы оценивали риск смерти и риск осложнений на уровнях общего PASC, систем органов и индивидуальных исходов в взвешенных когортах в течение трех временных периодов: 30–360 дней (первый год), 361–720 дней (второй год) и 721–1080 дней (третий год) после
T0. Для оценки риска исхода инцидента в каждом периоде участники рассматривались как подверженные риску, если исследуемый исход не имел места в предыдущем периоде. Мы оценивали веса шкалы предрасположенности независимо друг от друга в течение каждого периода и применяли веса разных периодов в одной модели результатов для оценки рисков и кумулятивной нагрузки. Участники подвергались цензуре на момент смерти или инфекции SARS-CoV-2 во время наблюдения как в группах COVID-19 (негоспитализированных и госпитализированных), так и в контрольной группе.
IRR, абсолютные коэффициенты, абсолютные различия коэффициентов, кумулятивные коэффициенты, кумулятивные различия коэффициентов смертности, количество осложнений и DALYsоценивались из взвешенных обобщенных оценочных уравнений с использованием логарифмической связи и распределения Пуассона. Различия в показателях смертности, числе осложнений и DALYs в целом и по системам органов между группами COVID-19 и контрольной группой без инфекции рассматривались как исходы, обусловленные COVID-19. Процентный вклад числа осложнений и DALYsв каждый год в течение периода наблюдения оценивался для общего PASC и по системам органов. 95% CIs были получены из 2,5-го и 97,5-го процентилей параметрического бутстрепа 1000 раз на основе точечных оценок и ковариационной матрицы обобщенных оценочных уравнений
20. Количество осложнений и DALYs сообщается в расчете на 1 000 человек.
Во всех анализах 95% CI IRR, исключающие единицу или количество последовательностей, и DALYs, исключающие ноль, считались свидетельством статистической значимости. Анализ проводился с использованием SAS Enterprise Guide версии 8.3 (SAS Institute), а результаты визуализировались с помощью R версии 4.3.2.
Анализ чувствительностиМы провели несколько анализов чувствительности.
(1) Мы оценили двойную устойчивость моделей корректировки, в которых ковариаты использовались как в моделях воздействия, так и в моделях результатов, вместо основного подхода, где ковариаты применялись только в модели воздействия.
(2) Вместо моделей Poisson (Пуассона) в первичном подходе мы построили модели Poisson(Пуассона) с нулевым раздуванием, чтобы оценить, как большое количество нулевых результатов влияет на соответствие модели.
(3) Вместо основного подхода, при котором участники подвергались цензуре при заражении SARS-CoV-2 во время наблюдения, мы не подвергали цензуре участников в группах COVID-19 при повторном заражении и рассматривали повторное заражение как естественный исход первого заражения.
(4) Вместо того, чтобы использовать только предопределенный набор ковариат в нашем первичном подходе, мы дополнительно скорректировали 100 алгоритмически выбранных высокомерных ковариат [40].
(5) Вместо того, чтобы определять госпитализацию во время острой фазы как дату госпитализации в течение 7 дней до или в течение 30 дней после положительного теста в основных анализах, мы использовали альтернативное определение госпитализации как дату госпитализации в день положительного теста или в течение 30 дней после положительного теста.
(6) Мы усекли вес оценки склонности на уровне 99,5% процентилей, а не 99,9% процентилей в основных анализах.
(7) Мы оценили результаты среди подвыборки с полными данными по всем ковариатам (
n = 4.432.414, 83,0% от полной выборки) для проверки согласованности результатов с результатами, полученными с помощью множественной импутации отсутствующих данных.
(8) Мы оценили риски на основе моделей Fine – Gray (Файн-Грея), в которых смерть и инфекция SARS-CoV-2 во время наблюдения рассматривались как конкурирующие риски [41].
(9) Мы дополнительно применили обратную вероятность цензурирования для учета неслучайной цензуры в связи со смертью или инфекцией SARS-CoV-2 во время наблюдения в трех группах (контрольная группа без инфекции и группы без госпитализации COVID-19 и госпитализированные группы COVID-19) [42].
(10) В качестве альтернативы мы использовали более узкое определение PASC, которое включало 73 исхода вместо 80 исходов, включенных в первичные анализы.
Контроль негативных исходовМы использовали тот же аналитический подход (описанный выше) для изучения связи между COVID-19 и инцидентами новообразования в качестве контроля негативного исхода в каждом году в течение 3 лет наблюдения [43]. Нет механистической поддержки или клинических доказательств причинно-следственной связи между инфекцией SARS-CoV-2 и риском возникновения новообразований. Воспроизведение априори ожидаемой нулевой связи между COVID-19 и контролем негативных исходов может снизить опасения по поводу возможных ложных предубеждений [43].
Сводка отчетностиДополнительная информация о дизайне исследования доступна в
Nature Portfolio Reporting Summary по ссылке на эту статью.
Доступность данныхДанные, подтверждающие выводы этого исследования, доступны в Министерстве по делам ветеранов США. Данные VA находятся в свободном доступе для исследователей за брандмауэром VA с утвержденным протоколом исследования VA. Для получения дополнительной информации посетите https://www.virec. research.va.gov или свяжитесь с Информационным ресурсным центром VA по адресу VIReC@va.gov.
Доступность кодаАналитический код доступен по адресу
https://github.com/caimiao 0714/ COVID3year.