Способ оценки риска возникновения заболеваний сердечно-сосудистой системы на основании средней скорректированной фазы сна по Мюнхенскому тесту
METHOD FOR ASSESSING THE RISK OF DEVELOPING DISEASES OF THE CARDIOVASCULAR SYSTEM BASED ON THE AVERAGE ADJUSTED SLEEP PHASE ACCORDING TO THE MUNICH TEST

Изобретение относится к области медицины, а именно к функциональной диагностике, общественному здоровью и здравоохранению, хронобиологии, может быть использовано при оценке риска возникновения некоторых заболеваний сердечно-сосудистой системы у лиц в возрасте от 18 до 65 лет включительно.
Скорость хода циркадных биологических часов оценивается величиной эндогенного циркадного периода. Величина эндогенного циркадного периода является видоспецифическим, а также индивидуальным генетически обусловленным признаком, колеблясь в диапазоне 24±4 часа при постоянных световых условиях, имея особенности у видов с дневным и ночным типом активности (Aschoff, 1960; Yan et al., 2020). Величина эндогенного циркадного периода остается постоянной в течение жизни (Duffy et al., 2001, 2011). У человека средняя скорость биологических часов несколько медленнее астрономического 24-часового цикла, средний эндогенный циркадный период составляет 24.1-24.2 часа, с половыми особенностями, у женщин достоверно ближе к 24 часам (Czeisler et al, 1999; Duffy et al., 2001; 2011; Wright et al., 2001, 2005; Woelders et al., 2017).
По данным литературы, эндогенный циркадный период взаимосвязан с устойчивостью к факторам десинхронизации биологических часов, а также с индивидуальной продолжительностью жизни (Wyse et al., 2010; Gutman et al., 2011; Libert et al., 2012; Hozer et al., 2020; Губин и Коломейчук, 2019). Несовпадение эндогенного циркадного периода с естественным суточным циклом увеличивает температуру тела, неэффективные энергетические затраты и снижает когнитивные показатели (Hozer & Pifferi, 2020) и коррелирует с подверженностью хроническим неинфекционным заболеваниям и смертностью от данных заболеваний (Губин и Коломейчук, 2019; Knutson & von Shantz, 2018).
Условия самоизоляции по причине пандемии COVID-19 предоставили уникальную возможность оценить скорректированную среднюю фазу сна в условиях нивелирования прессинга социального джетлага, оказывающего неравноценное воздействие на лиц с исходно разным хронотипом. С одной стороны, продолжительность сна большинства респондентов увеличилась, и скорректированная средняя фаза сна сместилась на более позднее время, что позволило частично компенсировать дефицит сна лиц с вечерним хронотипом (Korman et al., 2020). Однако результаты исследований также показали, что избыточный сдвиг скорректированной средней фазы сна на поздние часы коррелировал с ухудшением сна и самооценки здоровья.
УРОВЕНЬ ТЕХНИКИ
Для оценки суммарного риска возникновения заболеваний сердечно-сосудистой системы было разработано множество различных моделей.
Известна Шкала SCORE (Systematic COronary Risk Evaluation) разработаная для оценки риска смертельного сердечно-сосудистого заболевания в течение 10 лет. Для оценки суммарного риска смертельного сердечно-сосудистого заболевания с помощью калькулятора или специальных таблиц необходимо указать пол, возраст, уровень систолического АД, курит ли пациент и уровень общего холестерина. Полученная цифра представляет собой вероятность наступления смертельного исхода от сердечно-сосудистого заболевания в течение ближайших 10 лет, выраженную в процентах. В зависимости от полученного значения риска (в процентах) пациента следует отнести в одну из следующих категорий: низкий риск - менее 5%, высокий риск - 5% и более. (Conroy RM, Pyorala K, Fitzgerald AP, et al. Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. Eur Heart J 2003; 24:987-1003.)
К ограничениям метода можно отнести то, что он не применяется при установленном диагнозе сахарного диабета 1 и 2 типа, уровне холестерина выше 8,0 ммоль/л или АД больше 180/110 мм рт.ст. Не учитывает образ жизни пациента, требует проведения лабораторных и других исследований. Не учитывает хронобиологических показателей.
2
Применяется для оценки риска возникновения сердечно-сосудистых заболеваний индекс Кетле, так же называемый индексом массы тела. Индекс Кетле = массы тела (кг) : рост (м), по полученному значению определяется есть ли у человека ожирение и насколько высок риск возникновения у него сердечно-сосудистых заболеваний (National Institutes of Health (NIH), National Heart, Lung, and Blood Institute (NHLBI). The practical guide: identification, evaluation, and treatment of overweight and obesity in adults. Bethesda: National Institutes of Health. 2000, NIH publication 00-4084).
Данный метод не учитывает курение, уровень холестерина и глюкозы крови, а также образ жизни и другие факторы риска. Не учитывает хронобиологических показателей.
Известен способ выявления риска подверженности коронарной болезни сердца и инсульту (Патент RU 2356052 C2). Для осуществления способа определяют следующие параметры: концентрацию липопротеина низкой плотности в сыворотке крови или концентрацию С-реакционного протеина в плазме крови, систолическое или диастолическое давление крови, частоту биений сердца, радиус внутреннего канала артерии, температуру плазмы крови, угловой параметр артерии, параметр осевого положения диффузного потока в артерии. На основании полученных данных рассчитывают риски возникновения заболевания и определяют суммарный риск заболевания.
К недостаткам данного способа можно отнести необходимость инструментального, биохимического исследования и низкие возможности его применения в широкой медицинской практике. Не учитывает хронобиологических показателей.
Известен способ генетической диагностики подверженности к сердечно-сосудистым заболеваниям (Патент RU 2376372 C2). Способ основан на определении генотипов по полиморфным вариантам Al166C гена AGTR1, А-240Т и A2350G гена АСЕ, С677Т гена MTHFR, T174M гена AGT, C825T гена GNB3, VNTR4a/b и G894T гена NOS3, G1691A гена F5, PLA1/A2 гена ITGB3, G20210A гена F2 и оценке риска путем суммирования количества баллов, присвоенных каждому генотипу. При этом генотипам «низкого» риска сердечно-сосудистой патологии присваивается 0 баллов, к ним относятся генотипы 1166АА гена AGTR1, -240АА, 2350АА гена АСЕ, 677СС гена MTHFR, 174TT гена AGT, 825CC гена GNB3, VNTR4bb и 894GG гена NOS3, 1691GG гена F5, PLA1/A1 гена ITGB3, 20210GG гена F2. Генотипам «среднего» риска присваивается 0,5 баллов, к ним относятся генотипы 1166АС гена AGTR1, -240АТ и 2350AG гена АСЕ, 677СТ гена MTHFR, 174TM гена AGT, 825CT гена GNB3, VNTR4ab и 894GT гена NOS3. Генотипам «высокого» риска присваивается 1 балл, к ним относятся генотипы 1166СС гена AGTR1, -240ТТ и 2350GG гена АСЕ, 677ТТ гена MTHFR, 174MM гена AGT, 825TT гена GNB3, VNTR4aa и 894GT гена NOS3, 1691GA и 1691АА гена F5, PLA1/A2 и PLA2/A2 гена ITGB3, 20210GA и 20210АА гена F2. Риск сердечно-сосудистых болезней считается «низким» при сумме баллов от 0 до 3, средним - от 3,5 до 6, высоким - от 6,5 до 11 баллов.
К недостаткам метода можно отнести то, что комплексное генетическое исследование является дорогостоящим и в настоящее время не применимо при оценке рисков в широкой медицинской практике. Не учитывает хронобиологических показателей.
Известен способ прогнозирования гипертонической болезни и ишемической болезни сердца (Патент RU 2008803 C1). При котором диагностику осуществляют с использованием инструментального метода исследования, в частности эзофагогастродуоденоскопии и при наличии хронических эрозий слизистой оболочки желудка и луковицы двенадцатиперстной кишки прогнозируют развитие гипертонической болезни и ишемической болезни сердца.
К недостаткам данного способа можно отнести необходимость инструментального исследования и низкие возможности его применения в широкой медицинской практике. Не учитывает хронобиологических показателей.
Известен способ оценки суммарного риска развития сердечно-сосудистых заболеваний, специфический для российского населения (Патент RU 2352258 C2). При его применении строят регрессионную модель пропорционального риска , в которой в качестве факторов риска используют возраст индивидуума, систолическое артериальное давление, частоту сердечных сокращений, уровни общего холестерина и холестерина липопротеидов высокой плотности, индекс относительной массы тела, статус стенокардии напряжения и типичной для инфаркта миокарда боли, оценку электрокардиограммы в покое, закодированную по Миннесотскому коду, статусы курения и употребления алкоголя, определяемые в порядковой шкале. Рассчитывают абсолютный риск развития сердечно-сосудистых заболеваний по оригинальной математической формуле.
Несмотря на более высокую точность по сравнению с аналогами, данный метод требует обширного обследования пациента и довольно трудоемок, что так же затрудняет его использование в широкой медицинской практике. Не учитывает хронобиологических показателей.
2
Известен способ определения индивидуального профиля факторов суммарного сердечно-сосудистого риска у пациента трудоспособного возраста (Патент RU0002757752). Согласно которому проводят анамнестическое, клинико-инструментальное и лабораторные обследования пациента на предмет выявления у него значений следующих причинных факторов риска (ФР) ССЗ: систолическое артериальное давление (САД) в покое, мм рт.ст., частота сердечных сокращений (ЧСС) в покое, уд./мин, индекс массы тела (ИМТ), г/см, статус стенокардии напряжения, острый инфаркт миокарда (ОИМ) в анамнезе, острый инсульт (ОНМК) в анамнезе, сомнительные ЭКГ-признаки ишемической болезни сердца (ИБС), достоверно установленные ЭКГ-признаки ИБС, статус курения отдельно для мужчин и женщин, статус употребления алкоголя, уровень липопротеинов низкой плотности (ЛПНП), ммоль/л. Рассчитывают индивидуальный профиль факторов суммарного риска ССЗ пациента (ИСПР) на момент обследования. Для этого полученное при ранжировании значение каждого ФР умножают на соответствующий ему весовой коэффициент. Полученные произведения суммируют, получая значение ИСПР пациента. Дополнительно выявляют наличие у пациента подтвержденного диагноза сахарного диабета, либо, в его отсутствие, проводят клинико-инструментальные и лабораторные исследования на предмет его выявления у пациента. Принимая значение ИСПР пациента за 100%, определяют в нем долю - вклад в % каждого из полученных произведений индивидуальных значений ФР на их весовые коэффициенты, определяя таким образом персональный профиль причинных факторов сердечно-сосудистого риска данного пациента на момент обследования. В нем ФР с наибольшими значениями в % относят к главным индивидуальным причинным ФР ССЗ. Причем в случае наличия у пациента сахарного диабета, его также относят к главным индивидуальным причинным ФР данного пациента - требующим первоочередных профилактических или лечебных мер по снижению уровня индивидуального суммарного сердечно-сосудистого риска, в соответствии со значимостью выявленных причинных ФР в индивидуальном профиле пациента.
Хотя способ обеспечивает индивидуализацию главных причинных факторов риска заболеваний сердечно-сосудистой системы у конкретного пациента, в отличие от среднепопуляционного уровня риска, он требует обширного обследования пациента и довольно трудоемок, что так же затрудняет его использование в широкой медицинской практике. Не учитывает хронобиологических показателей.
Известен ряд способов прогноза риска сердечно-сосудистых заболеваний с применением технологии искусственного интеллекта: Метод и устройство прогнозирования риска на основе моделей, компьютерное оборудование и носитель информации (Патент CN113362952), Метод, модель и система прогнозирования риска на основе слияния гетерогенных признаков для ишемической болезни сердца (Патент CN109117864), Метод диагностики риска сердечно- сосудистых заболеваний на основе модели глубокой сверточной нейронной сети (Патент CN111000551) и ряд других.
К ограничениям данных методов следует отнести то, что обучение данных систем проводится на определенной когорте пациентов и не всегда представляется возможным экстраполировать эти данные на другие когорты. Не учитывает хронобиологических показателей.
За прототип взят способ оценки риска возникновения распространенных сопутствующих заболеваний на основании данных о хронотипе (крайние утренние, умеренные утренние, умеренные вечерние, крайние вечерние варианты хронотипа), определяемых с помощью шкалы хронотипов Morningness-Eveningness Questionnaire - MEQ (Horne & Ostberg, 1976) и данных полученных в ходе анализа более 430 тысяч анкет биобанка Великобритании (Knutson & von Shantz, 2018 DOI: 10.1080/07420528.2018.1454458). К недостаткам метода можно отнести отсутствие информации об особенностях СФС в рабочие и выходные дни в тесте MEQ и, соответствующей коррективы СФС.
Существующие способы оценки хронотипа трудоемки. Традиционный способ требует нахождения испытуемого в контролируемых лабораторных условиях в течение нескольких дней (Czeisler et al., 1999; Duffy et al., 2001; Wright et al., 2001, 2005; Chang et al., 2019).
В настоящее время с целью оценки хронотипа наиболее активно используется Мюнхенский тест (Munich Chrono-Type Questionnaire, MCTQ) (Roenneberg et al., 2003, 2019).
Технической  проблемой является отсутствие простого, не требующего дополнительных инструментальных и лабораторных исследований, способа оценки риска возникновения сердечно-сосудистых заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, только на основании данных о средней скорректированной фазе сна по Мюнхенскому тесту.
РАСКРЫТИЕ СУЩНОСТИ
Техническим результатом является оценка риска возникновения заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, на основании средней скорректированной фазы сна по Мюнхенскому тесту у лиц в возрасте от 18 до 65 лет включительно.
С использованием Мюнхенского теста MCTQ была получена уникальная наиболее масштабная глобальная база данных (Korman et al., 2020) проекта GCCS (Global Chrono Corona Survey) относительно особенностей хронотипа в условиях нивелирования социального прессинга раннего пробуждения, который в не равной мере действует на представителей разных хронотипов в повседневных условиях. Что позволило с более высокой точность определять эндогенный циркадный период человека по значению средней скорректированной фазы сна (далее СФС) по Мюнхенскому тесту MCTQ. Для максимально точной оценки исследование MCTQ следует проводить в конце отпускного (каникулярного) периода.
Расчет баллов риска проведен следующим образом. В основу взяты оригинальные данные (Knutson & von Shantz, 2018), после проведения дополнительно анализа данных (Gubin & Kolomeichuk, 2021) на основании выборки более 430 тыс. человек. Была проанализирована ассоциативная связь между хронотипами и природой заболеваемости. Данные по изученным патологиям были выражены в процентах, затем для каждого хронотипа была рассчитана средняя заболеваемость и расчет баллов риска возникновения заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, для лиц с различным хронотипом (Таблица 1).
Способ осуществляется следующим образом:
определяют пол, возраст, затем рассчитывают скорректированную среднюю фазу сна по Мюнхенскому тесту (далее СФС);
далее определяют риск возникновения заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, исходя из данных, представленных в таблице 1.
Таблица 1 Баллы риска возникновения заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, в зависимости от хронотипа.
<b>Эндогенный циркадный период (хронотипы)</b>
23,5 - 23,6
23,7 - 23,8
23,9 - 24,1
24,2 - 24,3
24,4 - 24,5
24,6 - 24,7
Крайние утренние
Умеренные утренние
Точные биоло-гические часы
Умеренные вечерние
Поздние вечерние
Крайние вечерние
<b>Процент риска относительно среднего</b>
106,1
102,2
98,0
97,0
98,5
99.7
<b>Оценка риска</b>
<b>2</b>
<b>1</b>
<b>0</b>
<b>0</b>
<b>0</b>
<b>1</b>
<b>Значение СФС (MSFsc)</b>
Мужчины от 18 до 22 включительно
3 ч 31 мин и ранее
3 ч 32 мин - 4 ч 05 мин
4 ч 06 мин - 4 ч 40 мин
4 ч 41 мин - 5 ч 14 мин
5 ч 15 мин - 5 ч 49 мин
5 ч 50 мин и позднее
Мужчины от 23 до 65 включительно
3 ч 01 мин и ранее
3 ч 02 мин - 3 ч 40 мин
3 ч 41 мин - 4 ч 19 мин
4 ч 20 мин - 4 ч 58 мин
4 ч 59 мин - 5 ч 38 мин
5 ч 39 мин и позднее
Женщины от 18 до 22 включительно
3 ч 20 мин и ранее
3 ч 21 мин - 3 ч 52 мин
3 ч 53 мин - 4 ч 24 мин
4 ч 25 мин - 4 ч 56 мин
4 ч 57 мин - 5 ч 28 мин
5 ч 29 мин и позднее
Женщины от 23 до 65 включительно
3 ч 00 мин и ранее
3 ч 01 мин - 3 ч 30 мин
3 ч 31мин - 4 ч 01 мин
4 ч 02 мин - 4 ч 31мин
4 ч 32 мин - 5 ч 02 мин
5 ч 03 мин и позднее
Интерпретация баллов повышения риска по данной нозологии: 0 баллов - риск отсутствует; 1 балл - минимальный риск (до 5%); 2 балла - небольшой риск (5-10 %).
ОСУЩЕСТВЛЕНИЕ ИЗОБРЕТЕНИЯ
Данный способ оценки риска возникновения заболеваний: ишемической болезни сердца, гипертонической болезни, атеросклероза, может применяться в направлениях общественное здоровье и здравоохранение, и хронобиология.
Пример осуществления изобретения:
Пример №1. Мужчина, возраст 19 лет, обратился в отделение медицинской профилактики для получения рекомендаций по коррекции образа жизни, рассчитанное значение скорректированной средней фазы сна по Мюнхенскому тесту MCTQ составило 4 ч 25 мин. Риск рассчитан - 0 баллов, оценен как отсутствие риска, рекомендации по коррекции образа жизни не требуются.
Пример №2. Мужчина, возраст 35 лет, обратился в отделение медицинской профилактики для получения рекомендаций по коррекции образа жизни, рассчитанное значение скорректированной средней фазы сна по Мюнхенскому тесту MCTQ составило 2 ч 55 мин. Риск рассчитан - 2 балла, оценен как небольшой. Рекомендовано: для снижения риска развития заболеваний дыхательной системы использование в первой половине дня средств светотерапии, биодинамического освещения, приближенного к естественному солнечному свету, пребывание до 17 часов под действием прямых солнечных лучей (на улице) не менее 60 минут в день, снижение экранного времени во второй половине дня.
Пример №3. Женщина, возраст 18 лет, обратилась в отделение медицинской профилактики для получения рекомендаций по коррекции образа жизни, рассчитанное значение скорректированной средней фазы сна по Мюнхенскому тесту MCTQ составило 3 ч 40 мин. Риск рассчитан - 1 балл, оценен как минимальный, рекомендации по коррекции образа жизни не требуются.
Пример №4. Женщина, возраст 38 лет, обратилась в отделение медицинской профилактики для получения рекомендаций по коррекции образа жизни, рассчитанное значение скорректированной средней фазы сна по Мюнхенскому тесту MCTQ составило 3 ч 50 мин. Риск рассчитан - 0 баллов, оценен как отсутствие риска, рекомендации по коррекции образа жизни не требуются.
ЛИТЕРАТУРА
1. Губин Д.Г., Коломейчук С.Н. ТОЧНОСТЬ БИОЛОГИЧЕСКИХ ЧАСОВ, ХРОНОТИП, ЗДОРОВЬЕ И ДОЛГОЛЕТИЕ. Хрономедицинский журнал (Тюменский медицинский журнал). 2019; 21(2): DOI: 10.36361/2307-4698-2019-21-2-14-27.
2. Aschoff J (1960) Exogenous and endogenous components in circadian rhythms. Cold Spring Harb Symp Quant Biol, 25, 11-28.
3. Brown S. A. et al. The period length of fibroblast circadian gene expression varies widely among human individuals. PLoS Biol 3, e338 (2005).
4. Chang AM, Duffy JF, Buxton OM, et al. Chronotype Genetic Variant in PER2 is Associated with Intrinsic Circadian Period in Humans. Sci Rep. 2019; 9(1):5350. doi:10.1038/s41598-019-41712-1.
5. Czeisler CA, Duffy JF, Shanahan TL, et al. Stability, precision, and near-24-hour period of the human circadian pacemaker. Science. 1999; 284(5423):2177-2181. doi:10.1126/science.284.5423.2177.
6. Duffy JF, Cain SW, Chang AM, et al. Sex difference in the near-24-hour intrinsic period of the human circadian timing system. Proc Natl Acad Sci U S A. 2011;  108 Suppl 3:15602-15608.
7. Duffy JF, Rimmer DW, Czeisler CA. Association of intrinsic circadian period with morningness-eveningness, usual wake time, and circadian phase. Behav Neurosci. 2001 Aug; 115(4):895-9. doi: 10.1037//0735-7044.115.4.895.
8. Gubin D.G., Kolomeichuk S.N., Weinert D. Circadian Clock Precision, Health, And Longevity. J. Chronomed. 2021; 23(1): 3-15. DOI: 10.36361/2307-4698-2020-23-1-3-15.
9. Gutman R, Genzer Y, Chapnik N, et al. Long-lived mice exhibit 24 h locomotor circadian rhythms at young and old age. Exp Gerontol. 2011; 46(7):606-9. doi: 10.1016/j.exger.2011.02.015.
10. Hida A, Kitamura S, Ohsawa Y, et al. In vitro circadian period is associated with circadian/sleep preference. Sci Rep. 2013; 3:2074. doi:10.1038/srep02074.
11. Horne JA, Ostberg OA. self-assessment questionnaire to determine morningness-eveningness in human circadian rhythms. Int J Chronobiol. 1976; 4(2):97-110.
12. Hozer C, Perret M, Pavard S, Pifferi F. Survival is reduced when endogenous period deviates from 24 h in a non-human primate, supporting the circadian resonance theory. Sci Rep. 2020; 10(1):18002. Published 2020 Oct 22. doi:10.1038/s41598-020-75068-8.
13. Hozer C, Pifferi F. Physiological and cognitive consequences of a daily 26 h photoperiod in a primate: exploring the underlying mechanisms of the circadian resonance theory. Proc Biol Sci. 2020 Jul 29;287(1931):20201079. doi: 10.1098/rspb.2020.1079.
14. Knutson KL, von Schantz M. Associations between chronotype, morbidity and mortality in the UK Biobank cohort. Chronobiol Int. 2018; 35(8):1045-1053. doi:10.1080/07420528.2018.1454458.
15. Korman M, Tkachev V, Reis C, et al. COVID-19-mandated social restrictions unveil the impact of social time pressure on sleep and body clock. Sci Rep. 2020; 10(1):22225. Published 2020 Dec 17. doi:10.1038/s41598-020-79299-7.
16. Libert S, Bonkowski MS, Pointer K, Pletcher SD, Guarente L. Deviation of innate circadian period from 24 h reduces longevity in mice. Aging Cell. 2012; 11(5):794-800. doi:10.1111/j.1474-9726.2012.00846.x.
17. Pagani L, Semenova EA, Moriggi E, et al. The physiological period length of the human circadian clock in vivo is directly proportional to period in human fibroblasts. PLoS One. 2010; 5(10):e13376. Published 2010 Oct 15. doi:10.1371/journal.pone.0013376.
18. Roenneberg T, Daan S, Merrow M. The art of entrainment. J Biol Rhythms. 2003 Jun; 18(3):183-94. doi: 10.1177/0748730403018003001.
19. Roenneberg T, Pilz LK, Zerbini G, Winnebeck EC. Chronotype and Social Jetlag: A (Self) Critical Review. Biology (Basel). 2019; 8(3):54. Published 2019 Jul 12. doi:10.3390/biology8030054.
20. Taillard J, Sagaspe P, Philip P, Bioulac S. Sleep timing, chronotype and social jetlag: Impact on cognitive abilities and psychiatric disorders. Biochem Pharmacol. 2021 Feb 2:114438. doi: 10.10 50.
21. Woelders T, Beersma DGM, Gordijn MCM, Hut RA, Wams EJ. Daily Light Exposure Patterns Reveal Phase and Period of the Human Circadian Clock. J Biol Rhythms. 2017; 32(3):274-286. doi:10.1177/0748730417696787.
22. Wright KP Jr, Gronfier C, Duffy JF, Czeisler CA. Intrinsic period and light intensity determine the phase relationship between melatonin and sleep in humans. J Biol Rhythms. 2005; 20(2):168-177. doi:10.1177/0748730404274265.
23. Wright KP Jr, Hughes RJ, Kronauer RE, Dijk DJ, Czeisler CA. Intrinsic near-24-h pacemaker period determines limits of circadian entrainment to a weak synchronizer in humans. Proc Natl Acad Sci USA. 2001; 98(24):14027-14032. doi:10.1073/pnas.201530198.
24. Wyse CA, Coogan AN, Selman C, Hazlerigg DG, Speakman JR. Association between mammalian lifespan and circadian free-running period: the circadian resonance hypothesis revisited. Biol Lett. 2010; 6(5):696-698. doi:10.1098/rsbl.2010.0152.
25. Yan L, Smale L, Nunez AA. Circadian and photic modulation of daily rhythms in diurnal mammals. Eur J Neurosci. 2020; 51(1):551-566. doi:10.1111/ejn.14172.