تسلل الأورام الدائري ووظيفة خلايا CD8+ T تحدد فعالية العلاج المناعي Circadian tumor infiltration and function of CD8+ T cells dictate immunotherapy efficacy

المجلة: Cell، المجلد: 187، العدد: 11
DOI: https://doi.org/10.1016/j.cell.2024.04.015
PMID: https://pubmed.ncbi.nlm.nih.gov/38723627
تاريخ النشر: 2024-05-01

خلية

التسلل الورمي اليومي ووظيفة CD8 تحدد خلايا T فعالية العلاج المناعي

ملخص رسومي

أهم النقاط

-كمية ووظيفة خلايا T المتسللة إلى الورم تعتمد على وقت اليوم
تغلغل الأورام الإيقاعي لخلايا T يعتمد على الساعة البيولوجية الوعائية
تعديل وقت الحقن خلال اليوم يحسن فعالية علاج خلايا CAR T
تحسين السيطرة على الورم بواسطة العلاج المضاد لـ PD-1 في الوقت المناسب يعتمد على نشاط خلايا T الدائرية

المؤلفون

تشين وانغ (王宸)، كون زينغ (曾群)، زينب ميليس غول، …، فيديريكو سيمونيتا، ميكائيل ج. بيتيت، كريستوف شيرمان

المراسلات

chen.wang.1@unige.ch (C.W.)، christoph.scheiermann@unige.ch (C.S.)

باختصار

يمكن تحسين فعالية علاج خلايا CAR T وحصار نقاط التفتيش المناعية في السرطان من خلال تعديل وقت العلاج خلال اليوم بسبب التسلل الإيقاعي ووظيفة خلايا T التي تنظمها الساعة البيولوجية في خلايا T والخلايا البطانية في بيئة الورم.

تسلل الأورام الدائري ووظيفة خلايا CD8+ T تحدد فعالية العلاج المناعي

تشن وانغ (وانغ تشن) كون زينغ (تسنج تشون) زينب ميلس غول، سي شي وانغ، روبرت بيك، فيل تشينغ، روبين بيل، يان وو، ستيفان ناولارتس، كولين بارنو بي تشون هسيوه، صوفي هيدلوند مولر، مارا سينيرينتي، منغزو صن، زيانغ سو، ستيفان جيميلان، فولوديمير بيترينكو، شارنا ديبنر، ستيفاني هيوز، كاميلا جاندوس، تشونغ وو لي، أوليفييه ميشيليان بينغ-تشي هو أبهيشيك دي. غارغ، فيديريكو سيمونيتا، ميكائيل ج. بيتيت، وكريستوف شيرمان قسم علم الأمراض والمناعة، كلية الطب، جامعة جنيف، 1211 جنيف، سويسرا مركز الأبحاث الانتقالية في الأورام – أمراض الدم (CRTOH)، جنيف 1211، سويسرا قسم الأورام وخدمة الأورام الدقيقة، مستشفيات جامعة جنيف، جامعة جنيف، جنيف 1211، سويسرا مركز أبحاث السرطان أغورا، لوزان 1011، سويسرا مركز علم الأحياء النظامية، معهد أبحاث مستشفى ماساتشوستس العام وكلية هارفارد الطبية، بوسطن، MA 02114، الولايات المتحدة الأمريكية المختبر الرئيسي لعلم الأورام والبحث الانتقالي، قسم علم الأمراض، مستشفى جامعة بكين للسرطان \ ومعهد، بكين 100142، الصين مختبر إجهاد الخلايا والمناعة، قسم الطب الخلوي والجزيئي، جامعة KU Leuven، لوفين 3000، بلجيكا قسم الأورام الأساسية، جامعة لوزان، لوزان 1066، سويسرا معهد لودفيغ لأبحاث السرطان، لوزان 1005، سويسرا قسم جراحة الصدر والغدد الصماء، قسم الجراحة، مستشفيات جامعة جنيف، جنيف 1205، سويسرا قسم فسيولوجيا الخلايا والتمثيل الغذائي، كلية الطب، جامعة جنيف، جنيف 1211، سويسرا مركز السكري، كلية الطب، جامعة جنيف، جنيف 1211، سويسرا معهد الوراثة والجينوم في جنيف (iGE3)، جنيف 1211، سويسرا مركز جنيف لأبحاث الالتهابات (GCIR)، جنيف 1211، سويسرا قسم أمراض الدم، قسم الأورام، مستشفيات جامعة جنيف، جامعة جنيف، جنيف، سويسرا المركز الطبي الحيوي (BMC)، معهد فسيولوجيا القلب والأوعية الدموية وعلم الأمراض، مركز والتر بريندل للطب التجريبي (WBex)، كلية الطب، جامعة لودفيغ ماكسيميليان (LMU) ميونيخ، بلانيغ-مارتينسريد 82152، ألمانيا ساهم هؤلاء المؤلفون بالتساوي جهة الاتصال الرئيسية*المراسلة: chen.wang.1@unige.ch (C.W.)، christoph.scheiermann@unige.ch (C.S.)https://doi.org/10.1016/j.cell.2024.04.015

الملخص

الملخص جودة وكمية اللمفاويات المتسللة إلى الورم، وخاصة تعتبر الخلايا معلمات مهمة للتحكم في نمو الورم واستجابته للعلاج المناعي. هنا، نوضح في سرطانات الفئران والبشر أن هذه المعلمات تظهر تذبذبات يومية، مدفوعة بكل من الساعة البيولوجية الداخلية لكريات الدم البيضاء وتسلل كريات الدم البيضاء الإيقاعي، الذي يعتمد على الساعة البيولوجية للخلايا البطانية في بيئة الورم الدقيقة. لاستغلال هذه الإيقاعات علاجياً، نوضح أن فعالية علاج الخلايا التائية المستقبلة لمستضدات الكيميرية وحجب نقاط التفتيش المناعية يمكن تحسينها من خلال ضبط وقت العلاج خلال اليوم. علاوة على ذلك، فإن توقيعات الخلايا التائية المعتمدة على وقت اليوم في نماذج أورام الفئران تتنبأ بالنجاة العامة لدى المرضى المصابين بالميلانوما وترتبط بالاستجابة للعلاج المضاد لـ PD-1. توضح بياناتنا الأهمية الوظيفية للديناميات اليومية في بيئة الورم الدقيقة وتقترح أهمية الاستفادة من هذه الميزات لتحسين تصميم التجارب السريرية المستقبلية ورعاية المرضى.

مقدمة

العلاج المناعي هو حالياً واحد من أكثر الوسائل فعالية في السيطرة على الأورام، وفي بعض الحالات، شفاء المرضى. تتضمن التدخلات الأكثر وعدًا حقن خلايا T المستقبلة للمستضدات الهجينة (CAR) أو مثبطات نقاط التفتيش المناعية (المعروفة أيضًا باسم blockers نقاط التفتيش المناعية [ICBs]). تتكون علاج خلايا CAR T من إعطاء خلايا T للمريض
التي تم تعديلها وراثيًا للتعرف على خلايا الورم. تهدف علاج ICB إلى تنشيط وظائف الخلايا التائية المضادة للورم؛ على سبيل المثال، فإن إعطاء جسم مضاد للبروتين 1 المبرمج لموت الخلايا (PD-1) يمنع الجزيئات المثبطة على سطح الخلايا التائية، مما يؤدي إلى تنشيطها.
لقد تم إثبات أن المناعة المضادة للأورام تعتمد مؤخرًا على وقت اليوم. على وجه التحديد، فإن الخلايا الشجرية أكثر فعالية في السيطرة على الأورام المزروعة في الفئران في أوقات محددة.
من اليوم. تخضع الوظائف العامة لجهاز المناعة لإيقاع يومي، تحكم. تخرج الكريات البيضاء من الدم إلى الأنسجة مثل العقد اللمفاوية بطريقة إيقاعية، الذي ارتبط بتحسين التعرف على المستضدات في أوقات محددة من اليوم. بالإضافة إلى خصائص الاتجار، فإن نمط الخلايا المناعية يتأثر أيضًا بوقت اليوم، بحيث يؤدي التحفيز في وقت محدد من اليوم إلى نتائج مختلفة في الاستجابات المناعية المرتبطة. لقد تم ربط ذلك بتفاعلات تعتمد على وقت اليوم تجاه اللقاحات في الفئران والبشر. وقد تم توسيعه مؤخرًا ليشمل أيضًا لقاحات مضادة للأورام. على الرغم من أن الجهاز المناعي يتمتع بإيقاع عالٍ في الحالة الثابتة وبعد التحدي الالتهابي، إلا أنه غير معروف ما إذا كانت الإيقاعات اليومية تنظم توزيع ونمط خلايا المناعة داخل الأورام.
نظرًا لأن الجهاز المناعي يخضع للتحكم اليومي، يبدو من المهم تحديد الأهمية المحتملة للإيقاع اليومي في سياق العلاج المناعي. هنا، نسعى لتحديد ما إذا كان بيئة الورم الدقيقة (TME) نفسها تظهر تذبذبات يومية. نحن نثبت أن الكريات البيضاء المتسللة إلى الورم (TILs) تظهر تذبذبات يومية في كل من العدد والنمط الظاهري، والتي يمكن استغلالها علاجيًا من خلال الإدارة الموقوتة لخلايا CAR T أو مثبطات المناعة (ICBs).

النتائج

تظهر تسرب الكريات البيضاء إلى الأورام تذبذبات يومية

استكشفنا ما إذا كان عدد الكريات البيضاء في الأورام يظهر اختلافات حسب الوقت من اليوم، من خلال حقن خلايا الميلانوما B16-F10 التي تعبر عن الألبومين البيضاوي (B16-F10-OVA) تحت الجلد (s.c.) في مجموعات من الفئران في وقت واحد من اليوم (أي في وقت زايجبير 9 [ZT9؛ 9 ساعات بعد بدء الضوء في بيئة 12 ساعة ضوء / 12 ساعة ظلام]) وجمع الأورام بعد 12 يومًا في 4 أوقات مختلفة من اليوم (ZT1 [“الصباح”], ZT7 [“النهار”], ZT13 [“المساء”], و ZT19 [“الليل”]). كان العدد الإجمالي لخلايا المناعة داخل الورم (TILs) يعتمد بشكل كبير على وقت جمع الأورام، حيث بلغ عدد الكريات البيضاء ذروته في المساء (ZT13) (الأشكال 1A و S1A). على وجه التحديد، اكتشفنا اختلافات حسب الوقت من اليوم في أعداد الأورام. و خلايا، خلايا خلايا بالإضافة إلى CD19 الخلايا (المتغيرة من حيث الطور عن عدد كريات الدم البيضاء)، في حين أن CD11b الخلايا و CD11b لم تظهر الخلايا تغييرات ملحوظة (الأشكال 1A و S1B و S1C). قمنا بتأكيد هذه البيانات من خلال صبغة المناعية للشرائح الورمية، حيث أظهرت أعدادًا أعلى من خلايا T في المساء (ZT13) مقارنةً بالصباح (ZT1) (الشكل 1B). لم نلاحظ اختلافات ملحوظة في حجم الورم خلال اليوم (الشكل S1D)، وعلى الرغم من أن الأورام استغرقت وقتًا أطول للنمو عند حصادها في الليل (ZT19) مقارنةً بالمساء (ZT13)، لاحظنا أعدادًا أعلى من الكريات البيضاء في المساء (الشكل 1A). وهذا يدل على أن الاختلافات في أعداد الكريات البيضاء كانت بسبب وقت حصاد الورم. كما لم نلاحظ أي تغييرات في كثافة الأوعية الدموية في الورم (الشكل S1E)، مما يشير إلى أن هذه الظاهرة خاصة بخلايا المناعة. قمنا بتأكيد الظاهرة في نموذج ورم الميلانوما العفوي (Tyr::Cre ; بتين ; براف ) (الأشكال 1C و S1A و S1F)، التي أظهرت معًا أن أعداد الخلايا اللمفاوية المتسللة (TILs) تظهر اختلافات حسب الوقت من اليوم، حيث تصل إلى ذروتها في المساء.
نقلنا الحيوانات بعد ذلك إلى ظروف الظلام التام، حيث تستمر الإيقاعات اليومية حتى في ظل ظروف بيئية ثابتة، مثل غياب أنماط الإضاءة الإيقاعية. لم يؤثر الظلام المستمر على الفروق الزمنية التي لوحظت خلال اليوم، مما يثبت أن هذه التذبذبات هي إيقاعات يومية حقيقية وليست مجرد استجابة لبيئة إضاءة:ظلام إيقاعية (الشكل 1D). ومع ذلك، أدى تحويل الفئران من دورة إضاءة:ظلام طبيعية إلى دورة مظلمة:مضيئة معكوسة مدتها 12 ساعة إلى عكس أعداد TIL (الأشكال 1E و1F). وهذا يدل على أن التذبذبات لم تكن تعتمد على الضوء بحد ذاته، ولكن يمكن تنسيقها مع نظام إضاءة جديد. هذه القدرة على التزامن مع إشارة خارجية هي ميزة إضافية للإيقاعات اليومية (الشكل 1F). كما أن تعريض الفئران لبروتوكول تأخر الطيران الحاد، الذي يتكون من تحويل الفئران إلى تأخير زمني مدته 6 ساعات وجمع الأنسجة في نقاط الوقت الصباحية والمسائية التي تم تقييمها سابقًا (الشكل 1E)، ألغى أيضًا الفروق الزمنية في TILs (الشكل 1F). معًا، توفر هذه البيانات دليلًا غير متوقع على أن TILs تظهر تذبذبات يومية في الخلوية، وهي عملية يتم تنسيقها بواسطة الضوء.

تتحكم خلايا البطانة في تسرب الكريات البيضاء الدائرية

لتقييم أسباب هذه التذبذبات، قمنا بتحديد ما إذا كانت تكاثر الكريات البيضاء، أو موتها، أو خروجها، و/أو تسللها تظهر اختلافات حسب الوقت من اليوم. لم نلاحظ اختلافات كبيرة في إجمالي تكاثر الكريات البيضاء، أو الاستماتة، أو الخروج (الأشكال S2A-S2C). بالمقابل، تم حظر حاد (24 ساعة) لتسلل الكريات البيضاء باستخدام مضاد لمولد الضد المرتبط بوظيفة اللمفاويات 1 (LFA-1)-إنتغرين (مضاد- ) الأجسام المضادة قللت بشكل كبير من عدد الكريات البيضاء في الورم في المساء (ZT13) وألغت الفروق الزمنية خلال اليوم (الشكل 2A). وهذا يشير إلى أن دخول الكريات البيضاء من الدم كان مسؤولاً عن الفروق الزمنية خلال اليوم وكان مطلوبًا للحفاظ على كثافة الخلايا المناعية في الورم. كما أشار إلى أن إمداد الكريات البيضاء إلى الورم من الدم كان ديناميكيًا بشكل مدهش على مدى فترة زمنية قصيرة. كانت تسلل الكريات البيضاء إلى الورم من الدم ذات صلة وظيفية حيث أن العلاج المزمن بمضاد LFA-1-إنتغرين زاد بشكل كبير من حجم الورم (الشكل S2D). للتحقق مما إذا كانت الكريات البيضاء تتسلل بالفعل إلى الورم بطريقة تعتمد على الوقت من اليوم، قمنا بإجراء اختبارات توطين قصيرة الأمد (ساعتين) مع كريات بيضاء تم نقلها بالتبني. تم جمع الخلايا من فئران المتبرعين في وقت معين من اليوم وتم حقنها عن طريق الوريد (i.v.) في متلقين تم تغيير توقيتهم، والذين تم وضعهم في خزائن مع أنظمة إضاءة مقلوبة، مما يحد من المتغير إلى التوقيت في المتلقي فقط. كان انتقال الكريات البيضاء إلى الورم يعتمد بشكل كبير على الوقت من اليوم، حيث وصلت المزيد من الخلايا إلى الورم في المساء، بما يتماشى مع بيانات حظر التسلل (الشكل 2B). وهذا أظهر أنه، في الفئران، كانت البيئة المجهرية للورم تظهر خصائص تعزز الهجرة في المساء، والتي لم تكن مدفوعة بفروق في كثافة الأوعية الدموية (الشكل S1E). نظرًا لأن -تمت معالجة تسرب الكريات البيضاء بواسطة مستقبلات الإنتغرين (الشكل 2A)، وقمنا بت quantifying مستويات التعبير عن الروابط الخاصة بها على خلايا بطانة الأوعية الدموية للورم، كحراس عند واجهة الدم-الورم. كانت جزيئات الالتصاق بين الخلايا 1 (ICAM-1) معبرة بشكل أعلى بكثير من قبل خلايا بطانة الأوعية في المساء (ZT13) مقارنة بالصباح (ZT1) (الشكل 2C)، في حين أن جزيئات الالتصاق الخلوية الوعائية 1 (VCAM-1) بالإضافة إلى E-selectin، وهي جزيئات التصاق إضافية لخلايا البطانة، لم تظهر اختلافات في وقت اليوم (الشكل S2E). كانت هذه الظاهرة لا تزال موجودة
الشكل 1. وقت الحصاد يحدد أعداد الكريات البيضاء المتسللة إلى الورم
(أ) الأعداد الكلية المعيارية لكريات الدم البيضاء المتسللة إلى الورم في ورم B16-F10-OVA، التي تم جمعها في 4 أوقات مختلفة من اليوم (وقت الزايتgeber [ZT])؛ 12 فأرًا من 4 تجارب مستقلة، تحليل كوزينور. المناطق المظللة تشير إلى مراحل الظلام.
(ب) التصوير (يسار) والتquantification (يمين) لـ CD4 و CD8 خلايا T في أورام B16-F10-OVA التي تم جمعها في ZT1 أو ZT13؛ فئران، اختبار ستودنت غير المقترن اختبار، قضبان القياس: .
(C) أعداد الكريات البيضاء المتسللة إلى الورم في Tyr::CreERT2، BRaf بتن الفئران التي تم جمعها في ZT1 أو ZT13؛ فئران، اختبار t لطلبة غير المقترنين.
(د) الأعداد الكلية المعنوية لكريات الدم البيضاء المتسللة إلى الورم في ورم B16-F10-OVA، التي تم جمعها في أوقات مختلفة من اليوم تحت ظروف ظلام دائم (الوقت اليومي [CT])؛ الفئران، من 3 تجارب مستقلة، تحليل التباين الأحادي.
جدول زمني خفيف من الضوء:الظلام (LD)، والظلام:الضوء المعكوس (DL)، وظروف اضطراب الساعة البيولوجية (JL). تشير المناطق المظللة إلى مراحل الظلام، وتدل الأرقام على أوقات الحصاد.
(F) أعداد الكريات البيضاء المتسللة إلى الورم في ورم B16-F10-OVA، التي تم جمعها في النقاط الزمنية المحددة (1 أو 13 ساعة في E) بعد بدء الدورة تحت الضوء: الظلام (LD، فئران)، مقلوب داكن:فاتح (DL، الفئران)، أو ما يتCorresponding jet lag (JL، ظروف (فئران) من 3 تجارب مستقلة، اختبار t لطلاب غير المتزاوج. المناطق المظللة تشير إلى المراحل المظلمة. جميع البيانات ممثلة كمتوسط SEM، غير دال إحصائيًا، . انظر أيضًا الشكل S1.
الشكل 2. خلايا البطانة تتحكم في تسرب الكريات البيضاء الدائرية وعلاج الخلايا التائية
(أ) الأعداد الكلية المعيارية لكريات الدم البيضاء المتسللة إلى الورم في نموذج B16-F10-OVA، التي تم جمعها في ZT1 أو ZT13، بعد 24 ساعة من علاج الأجسام المضادة للتحكم في النمط أو الأجسام المضادة المضادة لـ LFA-1؛ الفئران، من 5 تجارب مستقلة، اختبار t غير المتزاوج لستودنت.
(ب) استراتيجية البوابة وقياس الخلايا المنقولة بالتبني في أورام B16-F10-OVA بعد ساعتين من الحقن الوريدي عند ZT1 أو ZT13؛ الفئران، من تجربتين مستقلتين، اختبار ستودنت غير المقترن اختبار.
(ج) التصوير (يسار) والتquantification (يمين) للتعبير عن ICAM-1 على CD31 الخلايا البطانية في أورام B16-F10-OVA، التي تم جمعها في ZT1 أو ZT13 من المجموعة الضابطة (مجموعة من فئران B6 و Bmal1 الفئران)، NSG أو Bmal1 فئران؛ الأقسام، من الفئران، على التوالي، من 4 تجارب مستقلة، اختبار t لعينات غير مرتبطة وتحليل كوزينور، شريط القياس: .
(د) أعداد الكريات البيضاء في أورام B16-F10-OVA، التي تم جمعها في ZT1 أو ZT13 من مجموعة التحكم أو Bmal1 فئران؛ الفئران، من 4 تجارب مستقلة، اختبار ستودنت غير المقترن اختبار.
الأعداد المعيارية لكريات الدم البيضاء المنقولة بالتبني التي تم جمعها من أورام B16-F10-OVA بعد ساعتين من الحقن الوريدي في ZT1 أو ZT13 إلى مجموعة التحكم من نفس السلالة أو Bmal1 فئران؛ الفئران، من تجربتين مستقلتين، اختبار t غير المتزاوج لستودنت.
الأعداد المعيارية للخلايا المنقولة بالتبني OT-I في أورام B16-F10-OVA بعد ساعتين من الحقن في ZT1 أو ZT13؛ الفئران، من تجربتين مستقلتين، اختبار t غير المتزاوج لطلاب.
حجم الورم بعد الحقن الوريدي بـ خلايا OT-I المفعلة عند ZT1 ( الفئران) أو ZT13 ( فئران)، التحكم ( الفئران)، من 3 تجارب مستقلة، اختبار ت الطالب.
(H) حجم الورم بعد الحقن الوريدي لـ خلايا OT-I المفعلة عند ZT1 (Bmal1 الفئران، بمال1 الفئران) أو ZT13 ( فئران)، التحكم ( الفئران)، من تجربتين مستقلتين، اختبار t لستودنت، *Bmal1 ZT1 مقابل ZT13.
حجم الورم بعد الحقن الوريدي بـ CAR-CD28-CD3 البشري CD19 الخلايا في فئران NSG الحاملة لورم DoHH2 عند ZT1 (UT، خلايا T غير المعدلة، CAR T الفئران) أو ZT13 (UT، فئران، CAR T، فئران)، التحكم ( الفئران)، من تجربتين مستقلتين، تحليل التباين ثنائي الاتجاه.
(ج) عدد خلايا CAR T في أورام DoHH2 بعد 24 ساعة من الحقن الوريدي في ZT1 ( الفئران) أو ZT13 ( الفئران)، من تجربتين مستقلتين، اختبار t لطلاب غير المتزاوج. جميع البيانات ممثلة كمتوسط SEM، غير دال إحصائيًا، ; ; ; . انظر أيضًا الشكل S2.
في الفئران NOD scid gamma (NSG) التي تعاني من نقص المناعة الشديد، مما يشير إلى أن هذه الاختلافات لم تكن ناتجة عن التذبذبات الملحوظة في تسلل كريات الدم البيضاء في الورم (الشكل 2C). نحن
وبذلك تم التركيز على خلايا البطانة كوسائط محتملة لاختراق الكريات البيضاء الدائرية وزرع الأورام في حيوانات تظهر نقصًا قابلًا للتحفيز، محدد النسب، في الجين الدائري.
BMAL1 في الخلايا البطانية (Cdh5 :Bmal1 بمال1 ). في Bmal1 الفئران، تم إلغاء الفروق في تعبير ICAM-1 في خلايا بطانة الورم تمامًا (الشكل 2C)، ومع ذلك كانت كثافة الأوعية الدموية مشابهة للفئران من النوع البري (الشكل S1E). علاوة على ذلك، كانت الأعداد الإجمالية لكريات الدم البيضاء في الورم في Bmal1 توقفت عملية تسلل كريات الدم البيضاء المنقولة بالتبني إلى الورم في الفئران (الشكل 2D) عن الاعتماد على وقت اليوم (الشكل 2E). معًا، تُظهر هذه البيانات أن الخلايا البطانية تتحكم في عملية تسلل كريات الدم البيضاء الإيقاعية إلى الأورام، والتي تصل إلى ذروتها في المساء وتعتمد على آلية الساعة البيولوجية الذاتية للخلايا البطانية.

فعالية علاج خلايا CAR T تعتمد على وقت اليوم

نظرًا لأن خلايا T تمثل حاليًا الهدف الأكثر وعدًا للعلاج المناعي ضد السرطان، ركزنا في تجاربنا الوظيفية اللاحقة على هذه الخلايا. لتقييم جانب محتمل من الترجمة السريرية لاختراق خلايا T الإيقاعي في الأورام، قمنا في البداية بجمع خلايا CD8 المحددة لـ OVA. خلايا من فئران OT-I وأجرينا علاج الخلايا التائية على فئران مزروعة بأورام تعبر عن B16-F10-OVA. بعد النقل التبني الوريدي لخلايا OT-I المنشطة إلى الفئران الحاملة للأورام بعد 10 أيام من زراعة الورم، اكتشفنا زيادة ملحوظة في أعداد الخلايا المنقولة، عندما تم إعطاء حقن الخلايا التائية في ZT13 مقارنة بـ ZT1 (الشكل 2F). ومن المهم، أن هذا أدى إلى تحكم أفضل بشكل ملحوظ في الورم وتقليل عبء الورم (الأشكال 2G و S2F). بالمقابل، في Bmal1 في الفئران، تم إلغاء تأثيرات وقت اليوم، مما يدل على أن الفائدة العلاجية كانت محكومة بالساعة البيولوجية الذاتية للخلايا البطانية (الأشكال 2H و S2G). لتقييم جانب محتمل من الترجمة السريرية لاختراق خلايا T الإيقاعية في الأورام، استخدمنا خلايا CAR T البشرية، التي تحمل CAR مضاد للبشر CD19 (CD19 CAR-CD28-CD3). تم زراعة فئران NSG بخلايا DoHH2، وهي خط خلايا ورمية لسرطان الغدد الليمفاوية الكبيرة B المنتشر (DLBCL)، وهو المؤشر الرئيسي لعلاج خلايا CAR T في العيادة. تم علاج الفئران الحاملة للأورام بخلايا CAR T بطريقة تعتمد على وقت اليوم، حيث أظهرت حقن خلايا CAR T في ZT13 تقليلاً في حجم الورم، بينما لم يكن لحقنها في ZT1 أي تأثير على السيطرة على الورم (الأشكال 21 و S 2 H). تم ربط هذه الملاحظة بزيادة في توطين خلايا CAR T بعد 24 ساعة في ZT13 مقارنة بـ ZT1 في هذه الحيوانات (الشكل 2J)، بما يتماشى مع ملاحظتنا السابقة (الشكل 2B). توضح هذه البيانات بوضوح فعالية علاج خلايا T المعززة، ببساطة عن طريق تعديل وقت الحقن ليتناسب مع الوقت الأمثل لاختراق الكريات البيضاء خلال اليوم.

اختلافات يومية في نمط TIL

مع العلم أن عدد TILs يعتمد على وقت اليوم وأن هذه المتغير يمكن استغلاله لأغراض علاجية، قمنا بعد ذلك بتقييم ما إذا كان نمط TILs يعتمد أيضًا على وقت اليوم. لهذا الغرض، قمنا بإجراء تحليلات تسلسل RNA أحادي الخلية (scRNA-seq) لـ CD45 السكان، مرتبين حسب الأورام بعد 12 يومًا من الزرع وتم جمعهم في أربع أوقات مختلفة من اليوم، عند ZT1 و ZT7 و ZT13 و ZT19، مماثلة للتجارب السابقة. كشفت تحليلات الاقتراب والتقريب الموحد (UMAP) عن 18 مجموعة متميزة من كريات الدم البيضاء، مما يوفر رؤى حول المجموعات الفرعية وحالات كريات الدم البيضاء (الأشكال S3A-S3C و S4A و S4B). كانت المجموعة الأكثر بروزًا تتكون من البلعميات، التي شكلت
أعلى من التسلل، تليه خلايا القاتل الطبيعي (NK) ( ) وخلايا B ( (الأشكال S4A و S4B). تم التعبير عن جينات الساعة البيولوجية بشكل إيقاعي في جميع مجموعات الكريات البيضاء، حيث أظهر Per1 أكبر تذبذب، حيث بلغ ذروته عند ~ZT13 (الأشكال 3A و S4C). أظهرت البلعميات أعلى عدد من الجينات المتذبذبة، تليها وحيدات النواة وخلايا NK (الشكل S4D). ومع ذلك، أظهرت كل مجموعة من خلايا المناعة توقيعات جينية متذبذبة مميزة (الأشكال S4D-S4F؛ الجدول S1). أظهرت تحليلات علم وظائف الجينات أن الجينات المتذبذبة كانت غنية في الغالب في الأيض ولكنها أظهرت أيضًا اختلافات في التصاق الخلايا والتمايز بالإضافة إلى تنشيط خلايا T (الشكل S4E؛ الجدول S1). تشير هذه البيانات إلى أن عدد ونمط خلايا المناعة في البيئة المجهرية للورم (TME) متذبذب بشكل كبير ويعتمد بشدة على وقت اليوم.

اختلافات وقت اليوم في الورم نمط الخلية

مع التركيز على الخلايا اللمفاوية التائية، قمنا بإجراء تحليلات إضافية واكتشفنا 15 تحت مجموعة، مع الخلايا التي تمثل أكبر عدد من السكان (الأشكال 3B و3C وS5A-S5C). لاحظنا تغييرات ملحوظة تعتمد على وقت اليوم في نسبة تجمعات الخلايا المحددة على أنها “مضادة-” مقابل “مؤيدة للورم” اللمفاويات (انظر طرق STAR لتعريف التجمعات)، حيث بلغت ذروتها عند ZT13 وانخفضت عند ZT1 (الشكل 3D). علاوة على ذلك، كانت النسبة في التجمعات المحددة على أنها “غير مرهقة” مقابل “مرهقة” أظهرت الخلايا اختلافات، وهو نمط ظاهري لاحظناه أيضًا من خلال تحليل تدفق الخلايا (الأشكال 3E و S5D)، والذي يمكن تحسينه أكثر إلى نسبة الخلايا المحددة على أنها “ذاكرة فعالة” مقابل المستنفدة (الشكل 3D). وهذا يشير إلى أن أظهرت الخلايا نمطًا أكثر كبتًا في الصباح وتوقيعًا أكثر مضادًا للورم في المساء. للتحقيق في ذلك، قمنا بتعريف توقيعات الجينات لخلايا T المعبر عنها في المساء (ZT13) مقابل الصباح (ZT1) (الجدول S2). كان توقيع المساء غنيًا بمسارات التصاق الكريات البيضاء وتنشيط خلايا T (الشكل 3F) وكان مرتبطًا بشكل قوي بتحسين البقاء على قيد الحياة لدى مرضى الميلانوما البشرية من مجموعة بيانات أطلس جينوم السرطان (TCGA)، مما يشير إلى زيادة المناعية المضادة للورم لهذا التوقيع التعبيري للجينات (الشكل 3G). تظهر هذه البيانات أنه بالإضافة إلى التذبذبات في عدد خلايا T في البيئة المجهرية للورم، فإن وظيفتها تتغير أيضًا بطريقة تعتمد على وقت اليوم.

CD8 تذبذب نمط الخلية

تشير بيانات التعبير هذه إلى أن وظيفة خلايا T داخل بيئة الورم (TME) أظهرت اختلافات حسب الوقت من اليوم. أظهرت الجينات المرتبطة بتنشيط خلايا T، وكبحها، وهجرتها، ووظائفها الفعالة تذبذبًا في خلايا CD8 T الفعالة على مدار اليوم (الشكل 4A). ومن الجدير بالذكر أن تعبير جين نقطة التفتيش المناعية Pdcd1 (الذي يشفر PD-1) كان إيقاعيًا في مجموعة خلايا CD8 T هذه، حيث بلغ ذروته عند ZT1 وأظهر انخفاضًا عند ZT13 (الأشكال 4A و4B). وكان هذا متماشيًا مع تعبير Pdcd1 في إجمالي مجموعة خلايا CD8 T، كما تم تقييمه من خلال تحليلات qPCR لخلايا CD8 T المفروزة من الأورام وكذلك من خلال تحليلات pseudobulk لبيانات scRNA-seq (الأشكال 4C وS5E). كما تم التعبير عن PD-1 بشكل مختلف على مستوى بروتين سطح الخلية، بما يتماشى مع التحليلات النسخية، ولكنه كان متغير الطور، مشابهًا لما تم وصفه لزوجات mRNA-بروتين الأخرى المعبر عنها بشكل دوري. (الشكل 4D). يمكن إلغاء هذا الاختلاف في تعبير PD-1 عندما يتم حذف جين الساعة Bmal1 في
الشكل 3. الفروق اليومية في معايير خلايا T
(أ) تعبير جين الساعة البيولوجية Per1 في مجموعات كريات الدم البيضاء المتسللة إلى الورم مع منحنى كوسينور الملائم (خط) فترة الثقة (المنطقة المظللة).
(ب) تمثيل UMAP لتحليلات scRNA-seq لخلايا T التي تم جمعها في 4 أوقات مختلفة من اليوم.
(ج) الوفرة النسبية لمجموعات خلايا T عبر أوقات مختلفة من اليوم.
(د) نسبة وفرة مجموعات معينة من خلايا T في أوقات مختلفة من اليوم بواسطة تسلسل RNA أحادي الخلية. تشمل “الخلايا اللمفاوية المضادة للورم” خلايا NK، تشمل مجموعات “T”، و”الخلايا القاتلة الطبيعية T (NKT)”، و”CD8 الفعالة”، و”CD4 الساذجة”، و”CD8 ذاكرة الفعالة”، و”CD8 الانقسامية”، و”CD4 الفعالة”، بينما تشمل “الخلايا اللمفاوية المساعدة للورم” مجموعات “Treg”، و”Treg الانقسامية”، و”CD8 المستنفدة”.
نسبة طبيعية غير مستنفدة إلى مستنفدة الخلايا في أوقات مختلفة من اليوم بواسطة تحليل تدفق الخلايا؛ الفئران، من 5 تجارب مستقلة، تحليل كوزينور.
(F) تحليل علم الأحياء الجيني للجينات المتذبذبة في الخلايا.
(G) تحليل البقاء في مرضى الميلانوما من مجموعة بيانات TCGA باستخدام تعبير مرتفع أو منخفض لتوقيع جين ZT13 الفأري في الجدول S2، اختبار لوغرانك. جميع البيانات ممثلة كمتوسط SEM، * . انظر أيضًا الأشكال S3-S5 والجدول S1.
خلايا T (Cd4cre:Bmal1 بمال1 ) (الشكل 4E)، مما يشير إلى التعبير الإيقاعي لـ PD-1 في CD8 كانت الخلايا دورية وداخلية. لتأكيد ذلك، قمنا بتنشيط الخلايا في المختبر وتم تنسيق الخلايا (الأشكال 4F و S6A). تم تنشيط CD8 أظهرت خلايا T تعبيرًا متذبذبًا قويًا عن Pdcd1 بعد التزامن، والذي تم تخفيضه بشكل ملحوظ في خلايا T المجمعة من Bmal1. الفئران و Per1 الحيوانات (الشكل 4G). كما أظهر PDCD1 تذبذبات يومية في المختبر في خلايا بشرية متزامنة الخلايا، مما يشير إلى أن هذه التذبذبات لها أهمية في السياق البشري أيضًا (الشكل 4H). أظهرت هذه البيانات تعبير Pdcd1 في CD8 تكون خلايا T حقيقية من حيث إيقاع الساعة البيولوجية ومستقلة عن الخلايا، على مستوى كل من mRNA وبروتين السطح، مما قد يجعلها قابلة للاستهداف بواسطة مثبطات فحص المناعة المضادة لـ PD1 بطريقة تعتمد على وقت اليوم.

علاج مضاد PD-1 يعتمد على وقت اليوم

لذلك قمنا بتقييم أهمية التعبير المعتمد على وقت اليوم لبروتين PD-1 من خلال إجراء إعطاء مضادات PD-1 في أوقات محددة. كان لوقت إعطاء الأجسام المضادة تأثير قوي على نمو الورم في الميلانوما B16-F10-OVA.
نموذج الورم، حيث كانت الأورام تنمو بشكل أقل عندما تم إجراء العلاج في المساء مقارنةً بالصباح (الأشكال 5A و5B وS6B). في الواقع، كان للعلاج في الصباح تأثير ضئيل بشكل مدهش على عبء الورم، حتى عندما سبق إعطاء العلاج في الصباح إعطاء العلاج في المساء بـ 12 ساعة (الأشكال 5A و5B وS6B). يمكننا تأكيد هذه البيانات في نموذج سرطان القولون MC-38، مما يشير إلى أن هذه الظاهرة امتدت إلى نماذج ورم أخرى (الأشكال S6C وS6D). كان تأثير وقت اليوم معتمدًا بشكل حاسم على الخلايا، حيث لم تحقق علاج مضاد PD-1 أي فائدة في سيناريو حيث تم استنفاد الخلايا سابقًا (الأشكال 5C و S6E). علاوة على ذلك، كانت آلية ساعة الساعة البيولوجية لخلايا T مطلوبة حيث تم تغيير التذبذبات اليومية في Bmal1. الفئران، مما يدل على أن ساعة خلايا T تلعب دورًا في الوساطة الظاهرة (الأشكال 5D و S6F و S6G). باستخدام عدة دورات من العلاج المضاد لـ PD-1 في أوقات محددة وتركيبات من الإدارات المسائية والصباحية، لاحظنا أن وقت اليوم للعلاج الأول كان مسؤولًا عن الاختلاف، على عكس الجرعة اللاحقة (الأشكال 5E و S6H). وبالتالي، ركزنا على دورة الإدارة الأولى لتوضيح الآلية الأساسية وأجرينا تدفقًا
الشكل 4. تذبذب نمط الخلية
خريطة حرارية لتحليلات تسلسل RNA أحادي الخلية خلايا T الفعالة، تم جمعها في أربع أوقات مختلفة من اليوم.
(ب) تعبير Pdcd1 بواسطة تسلسل RNA أحادي الخلية لـ CD8 خلايا T الفعالة في أوقات مختلفة من اليوم.
(C) تعبير Pdcd1 بواسطة qPCR لفرز الخلايا المعتمد على الفلورية (FACS) لخلية CD8 الخلايا من أورام B16-F10-OVA عند ZT1 ( الفئران) أو ZT13 ( الفئران)، من 3 تجارب مستقلة، اختبار الطالب غير المتزاوج اختبار.
(د) تعبير PD-1 بواسطة تحليل تدفق الخلايا للورم الخلايا التي تم حصادها في أوقات مختلفة الفئران)، من 5 تجارب مستقلة، تحليل التباين الأحادي وتحليل الكوسينور.
تعبير PD-1 بواسطة تحليل التدفق الخلوي ل CD8 الورمي الخلايا، التي تم جمعها في ZT1 أو ZT13 من التحكم في الأشقاء الفئران) أو Bmal1 فئران ( فئران)، اختبار الطالب غير المتزاوج اختبار.
الضوء الحيوي النسبي المنشط الخلايا الناتجة من فئران PER2::LUC بعد التزامن، بيانات تمثيلية من فئران.
(G) تعبير Pdcd1 في CD8 المنشط الخلايا بعد التزامن من التحكم ( الفئران) أو بمال فئران ( الفئران) أو Per1 بير2 فئران ( الفئران) من 4 تجارب مستقلة، ***تحليل كوزينور. #تحليل التباين ثنائي الاتجاه.
(H) تعبير PDCD1 في خلايا CD8 المنشطة لدى البشر خلايا T بعد التزامن، ، من 4 تجارب مستقلة، تحليل كوزينور. جميع البيانات ممثلة كمتوسط SEM . انظر أيضًا الشكل S6.
تحليلات السيتومترية لـ الخلايا في الورم، بعد 24 ساعة من العلاج بمضاد PD-1 الذي تم في الصباح مقارنةً بالمساء. أظهرت النتائج أن إعطاء مضاد PD-1 في ZT13 وليس في ZT1 تسبب في زيادة كبيرة من تنشيط خلايا T وإفراز الحبيبات في الورم، كما يتضح من زيادة تعبير مستويات CD137 وCD107a، على التوالي (الشكل 5F). بالإضافة إلى ذلك، زادت معالجة مضاد PD-1 عند ZT13 أيضًا من إنتاج الإنترفيرون (IFN). وGranzyme B في الأورام المتسللة تمت مقارنة الخلايا مع ZT1 (الشكل 5G). لتعظيم التأثيرات المضادة للورم العلاجية الزمنية الملحوظة، قمنا بدمج علاج خلايا OT-I T في أوقات محددة مع حقن مضاد PD-1 في أوقات محددة. تم حقن خلايا OT-I T المنشطة مع جسم مضاد PD-1 في الفئران التي تم تغيير توقيتها في ZT13 أو في ZT1. وقد تم تعزيز التأثيرات المضادة للورم الآن، مع علاج ZT13.
تقليل عبء الورم بشكل أكبر مقارنةً بظروف ZT1 (الأشكال 5 H و S6I). وهذا يشير إلى أن كلا تحدد أعداد الخلايا والنمط الظاهري الفروق اليومية في الاستجابة للعلاج المناعي للأورام. معًا، تُظهر بياناتنا نمطًا ظاهريًا قويًا يعتمد على وقت اليوم لت infiltrate خلايا المناعة الورمية، والذي يمكن استغلاله مع الإدارة الموقوتة للعلاجات المناعية لتعظيم التأثيرات المضادة للأورام.

الظاهرة الدورية لخلايا المناعة في السرطانات البشرية

للاستفسار عما إذا كانت تقلبات وقت اليوم في تسلل خلايا المناعة موجودة أيضًا في الأورام البشرية، قمنا بتحليل عينات من الميلانوما البشرية من مرضى خضعوا لإزالة جراحية للورم في أوقات مختلفة من اليوم. باستخدام التحليلات المناعية النسيجية، وجدنا أن أعداد و
الشكل 5. الفعالية اليومية لعلاج مضاد PD-1
حجم الورم بعد علاج مضاد PD-1 (الأسهم) المطبق في ZT1 أو ZT13، مضاد PD-1 ( فئران)، التحكم ( الفئران)، من 3 تجارب مستقلة، تحليل التباين ثنائي الاتجاه.
(ب) حجم الورم بعد 14 يومًا من زراعة الورم، مضاد PD-1 ( فئران)، التحكم ( الفئران)، من 3 تجارب مستقلة، تحليل كوزينور.
حجم الورم بعد علاج anti-PD-1 (الأسهم) المطبق في ZT1 أو ZT13 مع أو بدون نقص CD8 في الفئران.
حجم الورم بعد علاج مضاد PD-1 (الأسهم) المطبق في ZT1 أو ZT13 في المجموعة الضابطة (ZT1، فئران، ZT13، الفئران) أو في Bmal1 (ZT1، فئران، ZT13، الفئران)، من 3 تجارب مستقلة، تحليل التباين ثنائي الاتجاه.
حجم الورم بعد جرعتين من علاج مضاد PD-1 (الأسهم) المقدم عند ZT1 أو ZT13؛ الفئران، من تجربتين مستقلتين، تحليل التباين ثنائي الاتجاه. (F) النمط الظاهري للورم خلايا تم جمعها في ZT1 أو ZT13، 24 ساعة بعد التحكم المتساوي الفئران) أو علاج مضاد PD-1 (ZT1، فئران، ZT13، الفئران)، من 4 تجارب مستقلة، تحليل التباين ثنائي الاتجاه، واختبار سيداكس بعد ذلك.
إنتاج السيتوكينات بواسطة خلايا CD8 الورمية خلايا T التي تم جمعها في ZT1 ( الفئران) أو ZT13 ( الفئران) 24 ساعة بعد علاج anti-PD-1، من 4 تجارب مستقلة، اختبار t غير المتزاوج.
حجم الورم بعد السيطرة ( الفئران) أو العلاج المركب من مضاد PD-1 وإدارة خلايا OT-I (السهم)، تم إعطاؤه في ZT1 ( الفئران) أو ZT13 ( 5 فئران)، من تجربتين مستقلتين، تحليل التباين الثنائي الاتجاه. الخطوط المنقطة تشير إلى نمو الورم من الشكل 2F. جميع البيانات ممثلة كمتوسط SEM، غير دال إحصائيًا، . انظر أيضًا الشكل S6.
CD8 أظهرت الخلايا نمطًا ملحوظًا يعتمد على وقت اليوم مع ذروة في وقت الظهيرة المبكرة (الأشكال 6A و6B). وهذا يدل على أن وقت اليوم كان له تأثير أيضًا على أعداد خلايا المناعة المتسللة في البشر. للحصول على معلومات أوسع فيما يتعلق بتوقيت تسلل خلايا المناعة في الأورام البشرية، استخدمنا مجموعات بيانات تسلسل RNA أحادي الخلية من الفئران مع مجموعات بيانات تسلسل RNA أحادي الخلية من جلد الفئران المتاحة علنًا والموسومة زمنياً، بالإضافة إلى مجموعات بيانات تسلسل RNA الكلي، لتدريب نموذج تعلم آلي قادر على توقيع مجموعات بيانات التسلسل حيث لم تكن المعلومات عن وقت اليوم متاحة (الشكل 6C). قمنا بتدريب خوارزمية ZeitZeiger، باستخدام إما جميع الجينات أو المعاملات المحفوظة (57 جينًا)، التي تم التنبؤ بها بواسطة ZeitZeiger لت oscillate (الشكل 6C؛ الجدول S3). كلا الطريقتين توقعتا وقت بيانات الاختبار بشكل جيد بما في ذلك الجلد الصحي وكذلك عينات الميلانوما – مع متوسط خطأ قدره (الأشكال S6J-S6L؛ الجدول S3). قمنا أيضًا بالتحقق من صحة نموذج ZeitZeiger على مجموعات بيانات RNA-seq وميكروأري جلد الفئران العامة وحققنا مرة أخرى توقيتًا جيدًا (الأشكال 6D و S6M). وهذا سمح لنا بتحديد الوقت من اليوم بدقة للأنسجة الصحية في الفئران وكذلك عينات الأورام لدينا بطريقة عمياء مع خطأ قدره باستخدام نموذج ZeitZeiger المدرب، هدفنا بعد ذلك إلى توقع وقت الحصاد أيضًا في خزعات بشرية. في الواقع، باستخدام البيانات العامة
البيانات المتاحة من المصفوفات الدقيقة من جلد الإنسان السليم حيث كانت معلومات وقت اليوم معروفة، تمكنا من التنبؤ بشكل كافٍ بوقت اليوم في هذه العينات مع خطأ قدره (الأشكال 6E و S6N). مع العلم أن الخوارزمية قامت بتعيين وقت اليوم بشكل مناسب لعينات البشر، قمنا بعد ذلك بتعيين معلومات وقت اليوم لـ 103 عينة من الميلانوما الأولية من مجموعة بيانات TCGA. ومن المثير للاهتمام أننا لاحظنا اختلافات في وقت اليوم في نسبة الخلايا المستنفدة الخلايا مقابل غير المستنفدة خلايا T، مماثلة لبيانات الفئران، ولكن – كما هو متوقع – متغيرة الطور بسبب التذبذبات المناعية المعكوسة بين البشر النهاريين والفئران الليلية، مع نسبة أعلى لمكافحة الورم في الصباح مقارنة بعد الظهر (الشكل 6F). تشير هذه البيانات إلى أن كل من أعداد الكريات البيضاء والنمط الظاهري في الخلايا اللمفاوية المتسللة إلى الورم البشرية تعتمد على وقت اليوم.
أخيرًا، تحققنا مما إذا كانت بصمة RNA المسائية التي تم تحديدها في خلايا الميلانوما في الفئران كافية للت correlate مع استجابة المرضى للعلاج المناعي. في الواقع، كشفت البصمة عن فرق كبير في معدل استجابة المرضى، مع المرضى الذين أظهروا تعبيرًا أعلى عن هذه العلامة في خلايا T أيضًا يظهرون استجابة إيجابية لعلاج ICB (الشكل 6G؛ الجدول S2). معًا، توفر بياناتنا دليلًا قويًا على أن حجرة خلايا المناعة من
الشكل 6. نمط الخلايا المناعية اليومية في السرطانات البشرية (A و B) التصوير (يسار) والتquantification (يمين) لـ CD4 المرضى) و CD8 خلايا T في الميلانوما البشرية، تم جمعها في أوقات مختلفة من اليوم، تحليل كوزينور، شريط القياس: تلوين CD8 أو CD4 باللون الأحمر.
(C) مخطط لتدريب الخوارزمية والتحقق منها للتنبؤ بوقت اليوم في العينات غير الموقعة زمنياً. تم استخدام خوارزمية نهائية تحتوي على 57 جينًا للتحقق الخارجي والتنبؤ.
(D و E) خطأ (ساعات) خوارزمية توقع الوقت في التحقق الخارجي في الفأر ( الفئران) (D) والبشر ( عينات) (E) الأنسجة.
(F) مرهق/غير مرهق نسبة خلايا T في أنسجة الميلانوما الأولية من مجموعة بيانات TCGA للميلانوما، تم رسمها باستخدام الخوارزمية للتنبؤ بوقت الخزعات، صباحاً، المرضى؛ بعد الظهر، المرضى، طالب اختبار.
(G) إثراء تعبير جين ZT13 في الفئران لدى المرضى الذين يستجيبون أو لا يستجيبون لعلاج ICB (غير المستجيبين، الخلايا؛ المستجيبين، الخلايا). يتم تمثيل جميع البيانات كوسيط وربع. انظر أيضًا الشكل S6 والجدولين S2 وS3.
سرطانات الميلانوما تعتمد بشكل مفاجئ على وقت اليوم، ويمكن استغلال ذلك علاجياً من خلال تدخلات مؤقتة مثل حقن خلايا CAR T أو حجب نقاط التفتيش المناعية.

نقاش

هنا، نقدم أدلة على الديناميات اليومية في الميكروبيئة المناعية للورم، مع تداعيات وظيفية للنقطة الزمنية المثلى لإدارة العلاجات المناعية. في المرضى الذين يعانون من الميلانوما، مرضى سرطان الرئة غير صغير الخلايا (NSCLC) والمرضى الذين يعانون من سرطان الخلايا الحرشفية النقيلي في المريء، أظهرت البيانات الاستعادية لمجموعات متطابقة تم علاجها بالعلاج المناعي القائم على الأجسام المضادة مؤخرًا أنها حساسة لوقت اليوم، حيث لوحظت معدلات بقاء أفضل بعد حقن العلاج المناعي في الصباح مقارنةً بعد الظهر. هذه الملاحظات السريرية تتماشى مع ما نحدده هنا بتفاصيل آلية. ستكون هناك حاجة إلى تجارب سريرية عشوائية ومستقبلية لدعم هذا الادعاء.
سؤال رئيسي في هذا المجال هو آلية تأثيرات وقت اليوم هذه. تظل مستويات الأجسام المضادة مرتفعة على مدى أسابيع في المرضى، مما قد يشير إلى أن التأثيرات اليومية الطويلة الأمد الملحوظة يتم توجيهها من خلال توقيت التدخلات الأولى. تشير بياناتنا إلى أن الدورة الأولى من العلاج تحقق نتائج يومية مختلفة وأن توقيت الجرعة (الجرعات) اللاحقة قد لا يكون ذا أهمية كبيرة. تدعم بيانات من تجارب سريرية استعادية الأهمية المحتملة لتوقيت الحقن الأولى لتحقيق الفوائد العلاجية. تشير بياناتنا إلى أن الأجسام المضادة التي تُحقن في أوقات مختلفة من اليوم تواجه اختلافات في كل من عدد الخلايا المناعية الموجودة في الورم وكذلك مستويات PD-1 لديها. وهذا يعني أن ناتج التغيرات اليومية في الأعداد ونمط الخلايا المناعية المتسللة، وبالتالي الهدف من الجسم المضاد، يظهر تذبذبات يومية، والتي قد تكون مسؤولة عن هذا التأثير الزمني العلاجي. ومع ذلك، لا يُعرف ما إذا كانت مستويات الأجسام المضادة المحلية تتذبذب، كما هو موجود داخل الفضاء بين الخلايا في الورم. لقد أظهرت الدراسات أن البلعميات الورمية تقوم بسرعة بإزالة الأجسام المضادة من خلايا T المجاورة. الذي قد يكون آلية إضافية لـ
حساسية وقت اليوم، نظرًا لأننا نرى أيضًا تذبذبات في نمط الخلايا البلعمية في بيئة الورم. لقد ركزنا في هذه الدراسة على خلايا CD8 T، حيث تمثل حاليًا أهدافًا علاجية رئيسية في العيادة. تظهر بياناتنا أن إعطاء الأجسام المضادة يغير نمط خلايا CD8 T داخل الورم بطريقة تعتمد على وقت اليوم، بعد 24 ساعة من الحقن، مما يشير إلى أن هذه الخلايا تشارك في الاستجابات اليومية الأولية. ومع ذلك، قد تكون التذبذبات في سلالات خلايا المناعة الأخرى التي نصفها، ولا سيما خلايا النخاع، ذات أهمية وظيفية إضافية للتأثير اليومي. لا يزال يتعين التحقيق في مساهمة مجموعات خلايا المناعة المختلفة في هذه الاستجابة المناعية المضادة للورم اليومية في الدراسات المستقبلية.
في النهاية، يبدو أن تحدي الجهاز المناعي بأجسام مضادة في وقت محدد من اليوم لا يغير فقط الكمية ولكن أيضًا جودة الاستجابة، بحيث أن الجهاز المناعي، بمجرد تحفيزه في الوقت “الخاطئ”، قد لا يكون قادرًا على الاستجابة بنفس المستوى والجودة كما هو الحال مع جهاز مناعي تم تحديه في الوقت “الصحيح” – بفارق 12 ساعة فقط. بطريقة مماثلة للاستجابة المضادة للأورام الموصوفة هنا، يوجد تأثير مخزن في اللقاحات، حيث يكون المستضد موجودًا على مدى عدة أيام. ومع ذلك، يتم ملاحظة تأثيرات وقت اليوم أيضًا في دراسات التطعيم، في الإعدادات ما قبل السريرية. وفي المرضى، مع وجود حساسية مناعية أعلى في فترات زمنية متشابهة جداً. لا تزال الآليات الدقيقة التي تكمن وراء هذه الظواهر بحاجة إلى تحديد.
تم الإبلاغ عن اختلافات بين الجنسين في علاج السرطان المناعي، حيث تشير التحليلات الاستعادية إلى أن الرجال غالبًا ما يستجيبون بشكل أفضل مقارنة بالنساء. من الجدير بالذكر، فيما يتعلق بالعلاج المناعي الزمني، أن تحليلات المجموعات الفرعية أظهرت أن النساء حققن فائدة عالية بشكل خاص في البقاء العام عند استخدام مثبطات نقاط التفتيش المناعية في وقت مبكر من اليوم مقارنةً بعد الظهر. إلى أي مدى وما هي الآليات التي تحتاج إلى التحقيق فيها بمزيد من التفصيل في الدراسات المستقبلية. يجب أن تثبت الديناميات اليومية في البيئة المجهرية للورم التي نكشف عنها هنا أنها مفيدة في العيادة، حيث تمثل خزعات الورم حاليًا الوسيلة الرئيسية لتقييم مستوى تسلل خلايا المناعة وتحديد النظام العلاجي الناتج. نظرًا لأن تسلل خلايا المناعة يعتمد بشكل كبير على وقت اليوم، قد يختلف العلاج اعتمادًا على متى تم إجراء الخزعات. معًا، نوضح بتفصيل آلي أن خلايا المناعة المتسللة إلى الورم تظهر تذبذبات يومية في كل من العدد والنمط الظاهري، والتي يمكن استغلالها علاجيًا من خلال الإدارة الموقوتة لخلايا CAR T أو مثبطات المناعة.

قيود الدراسة

ركزنا في هذه الدراسة على الميلانوما، ولا يزال يتعين إثبات ما إذا كان هذا ينطبق أيضًا على أورام أخرى في الدراسات المستقبلية. تشير ملاحظاتنا في الفئران إلى أن على الأقل الفروق الزمنية خلال اليوم في العلاج المناعي القائم على الأجسام المضادة تمتد أيضًا إلى سرطان القولون. علاوة على ذلك، تشير البيانات السريرية الاستعادية إلى أن تحسين توقيت الحقن قد يوفر فوائد لعدة أنواع من السرطان، بما في ذلك الميلانوما وسرطان الخلايا الحرشفية الميتاستاتي في المريء. على الرغم من أن بياناتنا تشير إلى أن النقطة الزمنية المثلى للحقن هي المساء في الفئران، قد تكون أفضل نقطة زمنية للبشر هي ساعات الصباح، نظرًا للاختلافات في فترات الراحة والسلوك النشط بين النوعين. ومع ذلك، لا يزال الوقت الدقيق المطلوب
محدد، حيث أن البيانات الحالية تقيم فقط تأثيرات الصباح مقابل بعد الظهر.
على الرغم من أن بياناتنا تشير إلى أن الحذف المحدد للسلالة والقابل للتحفيز لجين Bmal1 في الخلايا البطانية يؤثر على سلوك هذه الخلايا المتذبذب، وبالتالي على تجنيد الكريات البيضاء، قد تكون هناك عوامل إضافية يمكن أن تسهم في الإلغاء الملحوظ لهجرة الكريات البيضاء المعتمدة على وقت اليوم.

طرق النجوم *

تُقدم طرق مفصلة في النسخة الإلكترونية من هذه الورقة وتشمل ما يلي:
  • جدول الموارد الرئيسية
  • توافر الموارد
  • جهة الاتصال الرئيسية
  • توفر المواد
  • توفر البيانات والشيفرة
  • نموذج تجريبي وتفاصيل المشاركين في الدراسة – الحيوانات
  • تفاصيل الطريقة
  • خطوط خلايا الورم والتطعيم
  • تنشيط خلايا T في الفئران
  • تزامن خلايا T في الفئران
  • إنسان تنشيط خلايا T والتزامن
  • تعبير لوكفيراز خلايا
  • هضم الأنسجة وتحضير الخلايا المفردة
  • تدفق الخلايا
  • نقل كريات الدم البيضاء من الفئران
  • اختبار EdU
  • فرز خلايا CD8 T واستخراج RNA
  • استخراج RNA، النسخ العكسي و qPCR
  • علاجات الأجسام المضادة في الجسم الحي
  • توليد خلايا CAR T البشرية
  • نقل خلايا CAR T البشرية
  • تصوير المناعة الفلورية
  • التشخيص المناعي النسيجي البشري
  • تسلسل الخلايا المفردة والتحليل
  • تحليل زايتيزايجر
  • تقدير تركيبة الخلايا في مرضى الميلانوما في TCGA
  • التكميم والتحليل الإحصائي

معلومات إضافية

يمكن العثور على معلومات إضافية على الإنترنت فيhttps://doi.org/10.1016/j/cell. 2024.04.015.

شكر وتقدير

نشكر أورس ألبريخت ويورغن ريبيرجر على توفير فئران Per1-/-Per2-/ . نشكر منصة IGE3 للجينوم ومرافق تحليل الخلايا في جامعة جنيف. نشكر الدكتورة ليديا لوتس، والدكتورة لور غارنييه، والدكتورة فلور سينتوريل، والدكتور ليكينغ تشينغ على مدخلاتهم القيمة والنقاشات. تلقت هذه الدراسة دعمًا من المجلس الأوروبي للبحث (ERC CoG 101001233، CIRCADYN [C.S.])، ومن مؤسسة العلوم الوطنية السويسرية (SNSF، 310030_182417/1 [C.S.] و310030_185255 [S.H.])، ومن رابطة السرطان السويسرية (KLS-4836-08-2019 [C.S.])، ومن رابطة جنيف للسرطان (2106 [C.S.])، ومن برنامج الاتحاد الأوروبي ITN (813284، INTEGRATA [C.S.])، ومن مؤسسة HUG الخاصة (RC05-03 [C.S.]). يتلقى C.J. دعمًا من زمالة SNSF PRIMA (PR00P3_179727)، ومن رابطة السرطان السويسرية (KFS 5250-02-2021)، ومن رابطة جنيف للسرطان (GCL، 2007). M.C. هو الحاصل على جائزة راتب دكتوراه iGE3. حصل P.-C. Ho على دعم من معهد أبحاث السرطان، ومن تحالف أبحاث الميلانوما، ومن أبحاث السرطان لودفيغ، ومن مؤسسة العلوم الوطنية السويسرية (310030L_208130). حصل M.J.P. على دعم
من مؤسسة ISREC، أبحاث لودفيغ للسرطان، ومن منح المعهد الوطني للصحة P01CA240239 وR01-CA218579. تم تمويل R.B. من خلال زمالة تنقل ما بعد الدكتوراه ومنحة العودة من SNSF (P400PM_183852 وP5R5PM_203164). يتم دعم الأبحاث في مختبر A.D.G. من قبل مؤسسة أبحاث فلاندرز (FWO) (منحة البحث الأساسي، G0B4620N؛ منحة FWO SBO لمجموعة “ANTIBODY”)، جامعة KU Leuven (منحة C1، C14/19/098، ومنح C3، C3/21/037 وC3/22/022)، Kom op Tegen Kanker (KOTK/2018/11509/1 وKOTK/2019/11955/1)، VLIR-UOS (منحة iBOF، iBOF/21/048، لمجموعة “MIMICRY”)، وصندوق أبحاث أوليفيا هندريكس (OHRF Immuunbiomarkers). تم تمويل C.D. من خلال منحة مؤسسة العلوم الوطنية السويسرية 310030_184708/1، مؤسسة فوندوبل، مؤسسة نوفارتس لصحة المستهلك، برنامج EFSD/Novo Nordisk لأبحاث السكري في أوروبا، مؤسسة Swiss Life، مؤسسة أولغا ماينفيش، مؤسسة الابتكار في السرطان وعلم الأحياء، الرابطة الرئوية الجنيفية، رابطة السرطان السويسرية KFS-5266-02-2021-R، مؤسسة فيلوكس، مؤسسة ليناارد، مؤسسة ISREC، ومؤسسة جيرترود فون ميسنر (V.P. وC.D.). حصل Y.W. على دعم من مؤسسة العلوم لمستشفى سرطان بكين (2022-21). تم إنشاء البيانات الموسعة الشكل 6A والملخص الرسومي باستخدام BioRender.

مساهمات المؤلفين

تصور: C.W.، Z.L.، S.H.، C.J.، C.D.، A.D.G.، P.-C. Ho، O.M.، F.S.، M.J.P.، و C.S.؛ المنهجية: C.W.، Q.Z.، P.C.، S.N.، Z.S.، و A.D.G.؛ التحقيق: C.W.، Q.Z.، S.W.، R.B.، Y.W.، C.B.، R.P.، M.C.، M.S.، P.-C. Hsueh، S.H.M.، S.J.، V.P.، F.S.، Z.L.، و Z.M.G.؛ التصور: C.W.، Q.Z.، P.C.، S.N.، A.D.G.، و C.S.؛ الإشراف: C.S.؛ الكتابة: C.W.، Q.Z.، M.J.P.، و C.S.

إعلان المصالح

كان C.S. مستشارًا لشركة باير وتلقى أتعاب متحدث من Abbvie. كان M.J.P. مستشارًا لشركات أسترازينيكا، ديبيوفارم، إيلستار ثيرابيوتيكس، إيميونونشيا، KSQ ثيرابيوتيكس، ماكسي فاك، ميرك، مولكيولار بارتنرز، ثيرد روك فنتشرز، وتايدال. تلقى F.S. أتعاب استشارية من BMS/Celgene، إنسايت، كايت/جيلياد، أتعاب متحدث من كايت/جيلياد، إنسايت، دعم سفر من كايت/جيلياد، نوفارتس، أسترازينيكا، نيفي، جانسن، وتمويل بحثي من كايت/جيلياد، نوفارتس، وBMS/Celgene. تلقى A.D.G. أتعاب استشارية/استشارية/محاضرات أو تمويل بحثي من بويرنجر إنجلهايم، SOTIO، ميلتينيي بيوتيك، نوفيجينيكس، وإيزوبليكسيس. تلقى R.B. أتعاب متحدث من جانسن وهو متدرب في برنامج ENDEAVOUR-Breast لشركة دايشي سانكيو. جميع هذه العلاقات غير مرتبطة بهذه الدراسة.
تاريخ الاستلام: 12 أكتوبر 2023
تمت المراجعة: 2 فبراير 2024
تم القبول: 16 أبريل 2024
نُشر: 8 مايو 2024

REFERENCES

  1. Mellman, I., Coukos, G., and Dranoff, G. (2011). Cancer immunotherapy comes of age. Nature 480, 480-489. https://doi.org/10.1038/ nature10673.
  2. Tawbi, H.A., Schadendorf, D., Lipson, E.J., Ascierto, P.A., Matamala, L., Castillo Gutiérrez, E., Rutkowski, P., Gogas, H.J., Lao, C.D., De Menezes, J.J., et al. (2022). Relatlimab and Nivolumab versus Nivolumab in Untreated Advanced Melanoma. N. Engl. J. Med. 386, 24-34. https://doi. org/10.1056/NEJMoa2109970.
  3. Cappell, K.M., and Kochenderfer, J.N. (2023). Long-term outcomes following CAR T cell therapy: what we know so far. Nat. Rev. Clin. Oncol. 20, 359-371. https://doi.org/10.1038/s41571-023-00754-1.
  4. Flugel, C.L., Majzner, R.G., Krenciute, G., Dotti, G., Riddell, S.R., Wagner, D.L., and Abou-El-Enein, M. (2023). Overcoming on-target, off-tumour toxicity of CAR T cell therapy for solid tumours. Nat. Rev. Clin. Oncol. 20, 49-62. https://doi.org/10.1038/s41571-022-00704-3.
  5. Tumeh, P.C., Harview, C.L., Yearley, J.H., Shintaku, I.P., Taylor, E.J.M., Robert, L., Chmielowski, B., Spasic, M., Henry, G., Ciobanu, V., et al.
    (2014). PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature 515, 568-571. https://doi.org/10.1038/nature13954.
  6. Wang, C., Barnoud, C., Cenerenti, M., Sun, M., Caffa, I., Kizil, B., Bill, R., Liu, Y., Pick, R., Garnier, L., et al. (2023). Dendritic cells direct circadian anti-tumour immune responses. Nature 614, 136-143. https://doi.org/ 10.1038/s41586-022-05605-0.
  7. Scheiermann, C., Gibbs, J., Ince, L., and Loudon, A. (2018). Clocking in to immunity. Nat. Rev. Immunol. 18, 423-437. https://doi.org/10.1038/ s41577-018-0008-4.
  8. Arjona, A., Silver, A.C., Walker, W.E., and Fikrig, E. (2012). Immunity’s fourth dimension: approaching the circadian-immune connection. Trends Immunol. 33, 607-612. https://doi.org/10.1016/j.it.2012.08.007.
  9. Curtis, A.M., Bellet, M.M., Sassone-Corsi, P., and O’Neill, L.A.J. (2014). Circadian clock proteins and immunity. Immunity 40, 178-186. https:// doi.org/10.1016/j.immuni.2014.02.002.
  10. Cermakian, N., Stegeman, S.K., Tekade, K., and Labrecque, N. (2022). Circadian rhythms in adaptive immunity and vaccination. Semin. Immunopathol. 44, 193-207. https://doi.org/10.1007/s00281-021-00903-7.
  11. Man, K., Loudon, A., and Chawla, A. (2016). Immunity around the clock. Science 354, 999-1003. https://doi.org/10.1126/science.aah4966.
  12. Palomino-Segura, M., and Hidalgo, A. (2021). Circadian immune circuits. J. Exp. Med. 218, e20200798. https://doi.org/10.1084/jem. 20200798.
  13. Wang, C., Lutes, L.K., Barnoud, C., and Scheiermann, C. (2022). The circadian immune system. Sci. Immunol. 7, eabm2465. https://doi.org/10. 1126/sciimmunol.abm2465.
  14. He, W., Holtkamp, S., Hergenhan, S.M., Kraus, K., de Juan, A., Weber, J., Bradfield, P., Grenier, J.M.P., Pelletier, J., Druzd, D., et al. (2018). Circadian Expression of Migratory Factors Establishes Lineage-Specific Signatures that Guide the Homing of Leukocyte Subsets to Tissues. Immunity 49, 1175-1190.e7. https://doi.org/10.1016/j.immuni.2018.10.007.
  15. Druzd, D., Matveeva, O., Ince, L., Harrison, U., He, W., Schmal, C., Herzel, H., Tsang, A.H., Kawakami, N., Leliavski, A., et al. (2017). Lymphocyte Circadian Clocks Control Lymph Node Trafficking and Adaptive Immune Responses. Immunity 46, 120-132. https://doi.org/10.1016/j.immuni. 2016.12.011.
  16. Shimba, A., Cui, G., Tani-ichi, S., Ogawa, M., Abe, S., Okazaki, F., Kitano, S., Miyachi, H., Yamada, H., Hara, T., et al. (2018). Glucocorticoids Drive Diurnal Oscillations in T Cell Distribution and Responses by Inducing Interleukin-7 Receptor and CXCR4. Immunity 48, 286-298.e6. https:// doi.org/10.1016/j.immuni.2018.01.004.
  17. Casanova-Acebes, M., Nicolás-Ávila, J.A., Li, J.L., García-Silva, S., Balachander, A., Rubio-Ponce, A., Weiss, L.A., Adrover, J.M., Burrows, K., AGonzález, N., et al. (2018). Neutrophils instruct homeostatic and pathological states in naive tissues. J. Exp. Med. 215, 2778-2795. https://doi.org/ 10.1084/jem. 20181468.
  18. Nguyen, K.D., Fentress, S.J., Qiu, Y., Yun, K., Cox, J.S., and Chawla, A. (2013). Circadian gene Bmal1 regulates diurnal oscillations of Ly6C(hi) inflammatory monocytes. Science 341, 1483-1488. https://doi.org/10. 1126/science. 1240636.
  19. Cervantes-Silva, M.P., Carroll, R.G., Wilk, M.M., Moreira, D., Payet, C.A., O’Siorain, J.R., Cox, S.L., Fagan, L.E., Klavina, P.A., He, Y., et al. (2022). The circadian clock influences cell responses to vaccination by regulating dendritic cell antigen processing. Nat. Commun. 13, 7217. https:// doi.org/10.1038/s41467-022-34897-z.
  20. Ince, L.M., Barnoud, C., Lutes, L.K., Pick, R., Wang, C., Sinturel, F., Chen, C.S., de Juan, A., Weber, J., Holtkamp, S.J., et al. (2023). Influence of circadian clocks on adaptive immunity and vaccination responses. Nat. Commun. 14, 476. https://doi.org/10.1038/s41467-023-35979-2.
  21. Nobis, C.C., Dubeau Laramée, G., Kervezee, L., Maurice De Sousa, D., Labrecque, N., and Cermakian, N. (2019). The circadian clock of CD8 T cells modulates their early response to vaccination and the rhythmicity of related signaling pathways. Proc. Natl. Acad. Sci. USA 116, 2007720086. https://doi.org/10.1073/pnas.1905080116.
  22. Gibbs, J., Ince, L., Matthews, L., Mei, J., Bell, T., Yang, N., Saer, B., Begley, N., Poolman, T., Pariollaud, M., et al. (2014). An epithelial circadian clock controls pulmonary inflammation and glucocorticoid action. Nat. Med. 20, 919-926. https://doi.org/10.1038/nm.3599.
  23. Scheiermann, C., Kunisaki, Y., Lucas, D., Chow, A., Jang, J.E., Zhang, D., Hashimoto, D., Merad, M., and Frenette, P.S. (2012). Adrenergic nerves govern circadian leukocyte recruitment to tissues. Immunity 37, 290-301. https://doi.org/10.1016/j.immuni.2012.05.021.
  24. Hazan, G., Duek, O.A., Alapi, H., Mok, H., Ganninger, A., Ostendorf, E., Gierasch, C., Chodick, G., Greenberg, D., and Haspel, J.A. (2023). Biological rhythms in COVID-19 vaccine effectiveness in an observational cohort study of 1.5 million patients. J. Clin. Invest. 133, e167339. https://doi.org/ 10.1172/JCl167339.
  25. Long, J.E., Drayson, M.T., Taylor, A.E., Toellner, K.M., Lord, J.M., and Phillips, A.C. (2016). Morning vaccination enhances antibody response over afternoon vaccination: A cluster-randomised trial. Vaccine 34, 26792685. https://doi.org/10.1016/j.vaccine.2016.04.032.
  26. Ho, P.C., Bihuniak, J.D., Macintyre, A.N., Staron, M., Liu, X., Amezquita, R., Tsui, Y.C., Cui, G., Micevic, G., Perales, J.C., et al. (2015). Phosphoenolpyruvate Is a Metabolic Checkpoint of Anti-tumor T Cell Responses. Cell 162, 1217-1228. https://doi.org/10.1016/j.cell.2015.08.012.
  27. Raskov, H., Orhan, A., Christensen, J.P., and Gögenur, I. (2021). Cytotoxic CD8+ T cells in cancer and cancer immunotherapy. Br. J. Cancer 124, 359-367. https://doi.org/10.1038/s41416-020-01048-4.
  28. National Comprehensive Cancer Network (2024). B-Cell Lymphomas. https://www.nccn.org/guidelines/guidelines-detail?category=1&id=1480.
  29. Neelapu, S.S., Dickinson, M., Munoz, J., Ulrickson, M.L., Thieblemont, C., Oluwole, O.O., Herrera, A.F., Ujjani, C.S., Lin, Y., Riedell, P.A., et al. (2022). Axicabtagene ciloleucel as first-line therapy in high-risk large B-cell lymphoma: the phase 2 ZUMA-12 trial. Nat. Med. 28, 735-742. https://doi. org/10.1038/s41591-022-01731-4.
  30. Arjona, A., and Sarkar, D.K. (2005). Circadian oscillations of clock genes, cytolytic factors, and cytokines in rat NK cells. J. Immunol. 174, 76187624. https://doi.org/10.4049/jimmunol.174.12.7618.
  31. Hughey, J.J., Hastie, T., and Butte, A.J. (2016). ZeitZeiger: supervised learning for high-dimensional data from an oscillatory system. Nucleic Acids Res. 44, e80. https://doi.org/10.1093/nar/gkw030.
  32. Zhao, Y., Liu, M., Chan, X.Y., Tan, S.Y., Subramaniam, S., Fan, Y., Loh, E., Chang, K.T.E., Tan, T.C., and Chen, Q. (2017). Uncovering the mystery of opposite circadian rhythms between mouse and human leukocytes in humanized mice. Blood 130, 1995-2005. https://doi.org/10.1182/blood-2017-04-778779.
  33. Sade-Feldman, M., Yizhak, K., Bjorgaard, S.L., Ray, J.P., de Boer, C.G., Jenkins, R.W., Lieb, D.J., Chen, J.H., Frederick, D.T., Barzily-Rokni, M., et al. (2018). Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma. Cell 175, 998-1013.e20. https://doi. org/10.1016/j.cell.2018.10.038.
  34. Naulaerts, S., Datsi, A., Borras, D.M., Antoranz Martinez, A., Messiaen, J., Vanmeerbeek, I., Sprooten, J., Laureano, R.S., Govaerts, J., Panovska, D., et al. (2023). Multiomics and spatial mapping characterizes human CD8 T cell states in cancer. Sci. Transl. Med. 15, eadd1016. https://doi.org/ 10.1126/scitransImed.add1016.
  35. Qian, D.C., Kleber, T., Brammer, B., Xu, K.M., Switchenko, J.M., Jano-paul-Naylor, J.R., Zhong, J., Yushak, M.L., Harvey, R.D., Paulos, C.M., et al. (2021). Effect of immunotherapy time-of-day infusion on overall survival among patients with advanced melanoma in the USA (MEMOIR): a propensity score-matched analysis of a single-centre, longitudinal study. Lancet Oncol. 22, 1777-1786. https://doi.org/10.1016/S1470-2045(21) 00546-5.
  36. Yeung, C., Kartolo, A., Tong, J., Hopman, W., and Baetz, T. (2023). Association of circadian timing of initial infusions of immune checkpoint inhibitors with survival in advanced melanoma. Immunotherapy 15, 819-826. https://doi.org/10.2217/imt-2022-0139.
  37. Gonçalves, L., Gonçalves, D., Esteban-Casanelles, T., Barroso, T., Soares de Pinho, I., Lopes-Brás, R., Esperança-Martins, M., Patel, V., Torres, S., Teixeira de Sousa, R., et al. (2023). Immunotherapy around the Clock: Impact of Infusion Timing on Stage IV Melanoma Outcomes. Cells 12, 2068. https://doi.org/10.3390/cells12162068.
  38. Karaboué, A., Collon, T., Pavese, I., Bodiguel, V., Cucherousset, J., Zakine, E., Innominato, P.F., Bouchahda, M., Adam, R., and Lévi, F. (2022). Time-Dependent Efficacy of Checkpoint Inhibitor Nivolumab: Results from a Pilot Study in Patients with Metastatic Non-Small-Cell Lung Cancer. Cancers (Basel) 14, 896. https://doi.org/10.3390/cancers14040896.
  39. Rousseau, A., Tagliamento, M., Auclin, E., Aldea, M., Frelaut, M., Levy, A., Benitez, J.C., Naltet, C., Lavaud, P., Botticella, A., et al. (2023). Clinical outcomes by infusion timing of immune checkpoint inhibitors in patients with advanced non-small cell lung cancer. Eur. J. Cancer 182, 107-114. https://doi.org/10.1016/j.ejca.2023.01.007.
  40. Nomura, M., Hosokai, T., Tamaoki, M., Yokoyama, A., Matsumoto, S., and Muto, M. (2023). Timing of the infusion of nivolumab for patients with recurrent or metastatic squamous cell carcinoma of the esophagus influences its efficacy. Esophagus 20, 722-731. https://doi.org/10.1007/s10388-023-01006-y.
  41. Landré, T., Karaboué, A., Buchwald, Z.S., Innominato, P.F., Qian, D.C., Assié, J.B., Chouaïd, C., Lévi, F., and Duchemann, B. (2024). Effect of immunotherapy-infusion time of day on survival of patients with advanced cancers: a study-level meta-analysis. ESMO Open 9, 102220. https://doi. org/10.1016/j.esmoop.2023.102220.
  42. Martin, A.M., Cagney, D.N., Catalano, P.J., Alexander, B.M., Redig, A.J., Schoenfeld, J.D., and Aizer, A.A. (2018). Immunotherapy and Symptomatic Radiation Necrosis in Patients With Brain Metastases Treated With Stereotactic Radiation. JAMA Oncol. 4, 1123-1124. https://doi.org/10. 1001/jamaoncol.2017.3993.
  43. Postow, M.A., Goldman, D.A., Shoushtari, A.N., Betof Warner, A., Callahan, M.K., Momtaz, P., Smithy, J.W., Naito, E., Cugliari, M.K., Raber, V., et al. (2022). Adaptive Dosing of Nivolumab + Ipilimumab Immunotherapy Based Upon Early, Interim Radiographic Assessment in Advanced Melanoma (The ADAPT-IT Study). J. Clin. Oncol. 40, 1059-1067. https://doi. org/10.1200/JCO.21.01570.
  44. Arlauckas, S.P., Garris, C.S., Kohler, R.H., Kitaoka, M., Cuccarese, M.F., Yang, K.S., Miller, M.A., Carlson, J.C., Freeman, G.J., Anthony, R.M., et al. (2017). In vivo imaging reveals a tumor-associated macrophagemediated resistance pathway in anti-PD-1 therapy. Sci. Transl. Med. 9, eaal3604. https://doi.org/10.1126/scitransImed.aal3604.
  45. Conforti, F., Pala, L., Bagnardi, V., De Pas, T., Martinetti, M., Viale, G., Gelber, R.D., and Goldhirsch, A. (2018). Cancer immunotherapy efficacy and patients’ sex: a systematic review and meta-analysis. Lancet Oncol. 19, 737-746. https://doi.org/10.1016/S1470-2045(18)30261-4.
  46. Lévi, F.A., Okyar, A., Hadadi, E., Innominato, P.F., and Ballesta, A. (2024). Circadian Regulation of Drug Responses: Toward Sex-Specific and Personalized Chronotherapy. Annu. Rev. Pharmacol. Toxicol. 64, 89-114. https://doi.org/10.1146/annurev-pharmtox-051920-095416.
  47. Duan, J., Ngo, M.N., Karri, S.S., Tsoi, L.C., Gudjonsson, J.E., Shahbaba, B., Lowengrub, J., and Andersen, B. (2023). tauFisher accurately predicts circadian time from a single sample of bulk and single-cell transcriptomic data. Preprint at bioRxiv. https://doi.org/10.1101/2023.04.04.535473.
  48. Welz, P.-S., Zinna, V.M., Symeonidi, A., Koronowski, K.B., Kinouchi, K., Smith, J.G., Guillén, I.M., Castellanos, A., Furrow, S., Aragón, F., et al. (2019). BMAL1-Driven Tissue Clocks Respond Independently to Light to Maintain Homeostasis. Cell 177, 1436-1447.e12. https://doi.org/10. 1016/j.cell.2019.05.009.
  49. Wang, H., van Spyk, E., Liu, Q., Geyfman, M., Salmans, M.L., Kumar, V., Ihler, A., Li, N., Takahashi, J.S., and Andersen, B. (2017). Time-Restricted Feeding Shifts the Skin Circadian Clock and Alters UVB-Induced DNA Damage. Cell Rep. 20, 1061-1072. https://doi.org/10.1016/j.celrep. 2017.07.022.
  50. Tsujihana, K., Tanegashima, K., Santo, Y., Yamada, H., Akazawa, S., Nakao, R., Tominaga, K., Saito, R., Nishito, Y., Hata, R.-I., et al. (2022). Circadian protection against bacterial skin infection by epidermal CXCL14mediated innate immunity. Proc. Natl. Acad. Sci. USA 119, e2116027119. https://doi.org/10.1073/pnas.2116027119.
  51. Geyfman, M., Kumar, V., Liu, Q., Ruiz, R., Gordon, W., Espitia, F., Cam, E., Millar, S.E., Smyth, P., Ihler, A., et al. (2012). Brain and muscle Arnt-like protein-1 (BMAL1) controls circadian cell proliferation and susceptibility to UVB-induced DNA damage in the epidermis. Proc. Natl. Acad. Sci. USA 109, 11758-11763. https://doi.org/10.1073/pnas.1209592109.
  52. Spörl, F., Korge, S., Jürchott, K., Wunderskirchner, M., Schellenberg, K., Heins, S., Specht, A., Stoll, C., Klemz, R., Maier, B., et al. (2012). Krüp-pel-like factor 9 is a circadian transcription factor in human epidermis that controls proliferation of keratinocytes. Proc. Natl. Acad. Sci. USA 109, 10903-10908. https://doi.org/10.1073/pnas.1118641109.
  53. del Olmo, M., Spörl, F., Korge, S., Jürchott, K., Felten, M., Grudziecki, A., de Zeeuw, J., Nowozin, C., Reuter, H., Blatt, T., et al. (2022). Inter-layer and inter-subject variability of diurnal gene expression in human skin. NAR Genom. Bioinform. 4, Iqac097. https://doi.org/10.1093/nargab/Iqac097.
  54. Wu, G., Ruben, M.D., Schmidt, R.E., Francey, L.J., Smith, D.F., Anafi, R.C., Hughey, J.J., Tasseff, R., Sherrill, J.D., Oblong, J.E., et al. (2018). Popula-tion-level rhythms in human skin with implications for circadian medicine. Proc. Natl. Acad. Sci. USA 115, 12313-12318. https://doi.org/10.1073/ pnas. 1809442115.
  55. Sherrill, J.D., Finlay, D., Binder, R.L., Robinson, M.K., Wei, X., Tiesman, J.P., Flagler, M.J., Zhao, W., Miller, C., Loftus, J.M., et al. (2021). Transcriptomic analysis of human skin wound healing and rejuvenation following ablative fractional laser treatment. PLoS One 16, e0260095. https://doi.org/10.1371/journal.pone. 0260095.
  56. Hao, Y., Hao, S., Andersen-Nissen, E., Mauck, W.M., 3rd, Zheng, S., Butler, A., Lee, M.J., Wilk, A.J., Darby, C., Zager, M., et al. (2021). Integrated analysis of multimodal single-cell data. Cell 184, 3573-3587.e29. https:// doi.org/10.1016/j.cell.2021.04.048.
  57. Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., Feng, T., Zhou, L., Tang, W., Zhan, L., et al. (2021). clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb) 2, 100141. https://doi. org/10.1016/j.xinn.2021.100141.
  58. Carlucci, M., Kriščiūnas, A., Li, H., Gibas, P., Koncevičius, K., Petronis, A., and Oh, G. (2019). DiscoRhythm: an easy-to-use web application and package for discovering rhythmicity. Bioinformatics 36, 1952-1954. https://doi.org/10.1093/bioinformatics/btz834.
  59. Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14, 7. https://doi.org/10.1186/1471-2105-14-7.
  60. Colaprico, A., Silva, T.C., Olsen, C., Garofano, L., Cava, C., Garolini, D., Sabedot, T.S., Malta, T.M., Pagnotta, S.M., Castiglioni, I., et al. (2016). TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 44, e71. https://doi.org/10.1093/nar/ gkv1507.
  61. Leek, J.T., Johnson, W.E., Parker, H.S., Jaffe, A.E., and Storey, J.D. (2012). The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882-883. https://doi.org/10.1093/bioinformatics/bts034.
  62. Holtkamp, S.J., Ince, L.M., Barnoud, C., Schmitt, M.T., Sinturel, F., Pilorz, V., Pick, R., Jemelin, S., Mühlstädt, M., Boehncke, W.H., et al. (2021). Circadian clocks guide dendritic cells into skin lymphatics. Nat. Immunol. 22, 1375-1381. https://doi.org/10.1038/s41590-021-01040-x.
  63. Wittenbrink, N., Ananthasubramaniam, B., Münch, M., Koller, B., Maier, B., Weschke, C., Bes, F., de Zeeuw, J., Nowozin, C., Wahnschaffe, A., et al. (2018). High-accuracy determination of internal circadian time from a single blood sample. J. Clin. Invest. 128, 3826-3839. https://doi.org/ 10.1172/JCI120874.
  64. Newman, A.M., Steen, C.B., Liu, C.L., Gentles, A.J., Chaudhuri, A.A., Scherer, F., Khodadoust, M.S., Esfahani, M.S., Luca, B.A., Steiner, D., et al. (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 37, 773-782. https://doi. org/10.1038/s41587-019-0114-2.

طرق النجوم

جدول الموارد الرئيسية

كاشف أو مورد المصدر معرف
الأجسام المضادة
مضاد الفأر CD45 النسخة 30-F11 بي دي بيوساينس 564279/AB_2651134
مضاد الماوس CD3e النسخة KT3.1.1 بايو ليجند 155617/AB_2832541
النسخة GK1.5 المضادة للفأر CD4 بي دي بيوساينس 563232/AB_2738083
النسخة 53-6.7 من CD8a المضاد للفأر بي دي بيوساينس 563152/AB_2738030
النسخة N418 من مضاد CD11c للفأر بايو ليجند 117307/AB_313776
مضاد الفأر CD19 النسخة 1D3 بي دي بيوساينس 566412/AB_2744315
النسخة المضادة للفأر NK1.1 PK136 بايو ليجند 108715/AB_493591
النسخة المضادة لمستضدات الكريات البيضاء البشرية من نوع MHCII M5/114.15.2 بايو ليجند 107643/AB_2565976
مضاد الفأر Ly6G النسخة 1A8 بايو ليجند 127645/AB_2566317
مضاد الفأر Ly6C النسخة HK1.4 بايو ليجند 128023/AB_10640119
مضاد الفأر CD137 النسخة 17B5 بايو ليجند 106106/AB_2287565
مضاد الفأر CD107a النسخة 1D4B بايو ليجند 121610/AB_571991
مضاد PD-1 الفأري النسخة 29F.1A12 بايو ليجند 135213/AB_10689633
مضاد للجرانزيم B الفأري النسخة QA16A02 بايو ليجند 372207/AB_2687031
مضاد IFN- الفأر استنساخ XMG1.2 بايو ليجند 505837/AB_11219004
النسخة MEC13.3 من مضاد الفأر CD31 بايو ليجند 102520/AB_2563319
نسخة مضادة للفأر CD4، GK1.5 بايو ليجند 100405/AB_312690
نسخة مضادة للفأر CD54، YN1/1.7.4 بايو ليجند 116114/AB_493495
النسخة 53-6.7 من مضاد الفأر CD8 بايو ليجند 100724/AB_389326
النسخة 429 من الأجسام المضادة ضد VCAM-1 للفئران بايو ليجند 105710/AB_493427
مضاد الفأر مضاد CD11a النسخة M17/4 بايو إكس سيل BE0006/ AB_1107578
مضاد الفأر مضاد-CD18 النسخة M18/2 بايو إكس سيل BE0009/ AB_1107607
مضاد PD-1 للفأر، النسخة 29F.1A12 بايو إكس سيل BE0273/ AB_2687796
مضاد الفأر CD8a، النسخة YTS 169.4 بايو إكس سيل BE0117/ AB_10950145
النسخة UZ6 المضادة للبروتين E-selectin للفئران إنفيتروجين MA1-06506/ AB_2186701
نسخة مضادة للبشر CD3 BW264/56 ميلتيني بيولوجي 130-113-687/AB_2726228
نسخة مضادة للبشر CD8 RPA-T8 بي دي بيوساينس 563795/AB_2722501
نسخة CD4 المضادة للبشر SK3 بي دي بيوساينس 565995/AB_2739446
النسخة المضادة للبشر CD19 HIB19 بايو ليجند 302219/AB_389313
نسخة مضادة للبشر CD4 EP204 جسر تشونغشان الذهبي للتكنولوجيا الحيوية ZA-0519
نسخة CD8 المضادة للبشر SP16 جسر تشونغشان الذهبي للتكنولوجيا الحيوية ZA-0508
تحكم إيزوتيب IgG2a للفأر، النسخة 2A3 بايو إكس سيل BE0089/ AB_1107769
تحكم إيزوتيب IgG2b للفأر، النسخة LTF-2 بايو إكس سيل BE0090/ AB_1107780
دينابيدز ت Activator CD3/CD28 البشري جيبكو 11132D/AB_2916088
نسخة مضادة للبشر CD3 BW264/56 بايو ليجند 317348/ AB_2571995
دينابيدز نشاط T-Activator CD3/CD28 للفئران لتوسيع وتنشيط خلايا T جيبكو 11453D
أجسام مضادة ثانوية ماعز مضادة لجرذان IgG (H+L) تم امتصاصها بشكل متقاطع إنفيتروجين A21247/ AB_141778
المواد الكيميائية، الببتيدات، والبروتينات المؤتلفة
DRAQ7 بيوليجند 424001
عدّ الخرز ثيرمو فيشر C36950
(يتبع في الصفحة التالية)
مستمر
كاشف أو مورد المصدر معرف
إي بايوساينس صبغة القابلية للإصلاح eFluor إي بايوساينس 65-0865-18
مجموعة صبغة عامل النسخ Foxp3 eBioscience 00-5523-00
فوربول 12-ميريستات 13-أسيتات سيغما-ألدريتش P1585
أيونيوميسين سيغما-ألدريتش I3909
جولجي بلوق بي دي بيوساينس 555029
حدث الخلية كاسبيز-3/7 الأخضر جاهز للتدفق كاشف إنفيتروجين R37167
متتبع الخلايا صبغة حمراء داكنة إنفيتروجين C34565
محلول تحليل كريات الدم الحمراء بيوليجند 420302
إدو أليكس فلور مجموعة اختبار تحليل تدفق الخلايا إنفيتروجين C10646
رياجنت تريزول إنفيتروجين 15596018
إنترلوكين-2 البشري بيبروتيك 200-02
IL-2 الفأر بيوليجند 575406
ريترونيكتين تاكارا T100B
ستربتافيدين بيوليجند 405207
بروتين L المربوط بالبيوتين ثيرمو فيشر ٢٩٩٩٧
كولاجيناز IV شركة وورثينغتون للمواد الكيميائية الحيوية LS004189
كولاجيناز د روش ١١٠٨٨٨٦٦٠٠١
دي إنز I روش 04716728001
DAPI بيوليجند 422801
بلوك إيد حل الحجب إنفيتروجين B10710
مجموعة كشف البوليمر بوند ريد لايكا دي إس 9390
اختبارات تجارية حاسمة
مجموعة كواشف كروميم أحادي الخلية 3′ v3.1 مع مؤشرات مزدوجة 10x جينومكس 1000269
طقم عزل خلايا T CD8+، بشري ميلتيني بيولوجي ١٣٠-٠٩٦-٤٩٥
مجموعة عزل خلايا T CD8a+، فئران ميلتيني بيولوجي ١٣٠-١٠٤-٠٧٥
البيانات المودعة
تسلسل RNA أحادي الخلية في الوقت المحدد من الخلايا اللمفاوية التائية المستخرجة من الميلانوما الفأرية هذه الورقة جيو: GSE260641
تسلسل RNA أحادي الخلية من البشرة في الفأر دوان وآخرون جيو: GSE223109
تسلسل RNA لجلد الفأر ويلز وآخرون جيو: GSE115104، GSE114943
تسلسل RNA لجلد الفأر وانغ وآخرون جيو: GSE83855
تسلسل RNA في البشرة للفأر في ظلام دائم تسوجيهانا وآخرون جيو: GSE174155
مصفوفة ميكروية لجلد الفأر جيفمان وآخرون جيو: GSE38622، GSE38623
ميكروأري الفقاعات الشفطية للبشرة البشرية سبورل وآخرون جيو: GSE35635
ميكروأري جلد الإنسان دل أولمو وآخرون جيو: GSE205155
ميكروأري جلد الساعد البشري وو وآخرون جيو: GSE112660
ميكروأري جلد الإنسان (علاج بالليزر الجزئي) شيريل وآخرون جيو: GSE139305
تسلسل RNA أحادي الخلية لميلانوما الإنسان ساد-فيلدمن وآخرون جيو: GSE120575
بيانات داعمة أخرى هذه الورقة (https://doi.org/10.26037/yareta: 5rmtImd5sreo7bmcbe44bbijfm).
نماذج تجريبية: خطوط الخلايا
B16-F10-OVA ستيفاني هيوغ، جامعة جنيف، سويسرا غير متوفر
مستمر
كاشف أو مورد المصدر معرف
خط الخلايا السرطانية القولونية MC38 من الفئران مارك ج. سميث RRID: CVCL_B288
لمفوما بشرية DoHH2 فرانشيسكو بيرتوني، IOR، بيلينزونا، سويسرا غير متوفر
نماذج تجريبية: الكائنات/السلالات
فأر: WT C57BL/6J مختبر جاكسون السلالة # 000664 RRID:IMSR_JAX:000664
فأر: WT C57BL/6N نهر تشارلز السلالة # CRL027 RRID:IMSR_CRL:027
فأر: B6.129S4(Cg)-Bmal1tm1Weit/J مختبر جاكسون
السلالة #:007668
RRID:IMSR_JAX:007668
فأر: B6.Cg-Tg(Cd4-cre)1Cwi/BfluJ مختبر جاكسون
السلالة #:022071
RRID:IMSR_JAX:022071
فئران Cdh5-cre/ERT2 رالف آدامز، MPI مونستر، ألمانيا غير متوفر
NSG (NOD.Cg-PrkdcscidIL2rgtm1Wjl/SzJ) نهر تشارلز السلالة #:614
BRaf بتن تير::كريERT2 بينغ-تشي هو، جامعة الأمم المتحدة، سويسرا هو وآخرون
CD45.1 OT-I والتر رايث، جامعة جنيف، سويسرا غير متوفر
أوليغونوكليوتيدات
فأر Pdcd1 F
CGGTTTCAAGGCATGGTCATTGG
مايكروسينث غير متوفر
فأر Pdcd1 R
TCAGAGTGTCGTCCTTGCTTCC
مايكروسينث غير متوفر
فأر Rplp0
GAGACTGAGTACACCTTCCCAC
مايكروسينث غير متوفر
فأر Rplp0 R
ATGCAGATGGATCAGCCAGG
مايكروسينث غير متوفر
فأر Rpl32 F
ACAATGTCAAGGAGCTGGAG
مايكروسينث غير متوفر
فأر Rpl32 R
TTGGGATTGGTGACTCTGATG
مايكروسينث غير متوفر
PDCD1-F البشري
CCAGGATGGTTCTTAGACTCCC
مايكروسينث غير متوفر
بروتين PDCD1 البشري
TTTAGCACGAAGCTCTCCGAT
مايكروسينث غير متوفر
البروتين البشري Rplp0 F
TGGTCATCCAGCAGGTGTTCGA
مايكروسينث غير متوفر
البروتين البشري Rplp0
ACAGACACTGGCAACATTGCGG
مايكروسينث غير متوفر
البرمجيات والخوارزميات
فاسك ديفا 6 بي دي بيوساينس RRID:SCR_001456
فلو جو 10 بي دي بيوساينس RRID:SCR_008520
سلايدبوك 3 آي RRID:SCR_014300
إيميج جي فيجي RRID:SCR_002285
سيل رينجر 6.1.2 10X جينوميكس RRID:SCR_017344
سورات 4.3.0 هاو وآخرون https://satijalab.org/seurat/
R 4.2.2 مؤسسة R https://www.r-project.org/
كلاستر بروفايلر وو وآخرون https://bioconductor.org/packages/ release/bioc/html/clusterProfiler.html
كوزينور المصدر المفتوح https://github.com/sachsmc/cosinor
إيقاع الديسكو كارلوتشي وآخرون https://bioconductor.org/packages/ release/bioc/html/DiscoRhythm.html
GSVA هانزلمان وآخرون https://bioconductor.org/packages/ release/bioc/html/GSVA.html
(يتبع في الصفحة التالية)
مستمر
المُعَايِن أو المورد المصدر معرف
TCGAbiolinks كولابريكو وآخرون https://bioconductor.org/packages/ release/bioc/html/TCGAbiolinks.html
زيت زايجر هيوي وآخرون https://github.com/hugheylab/zeitzeiger
سفا ليك وآخرون https://bioconductor.org/packages/ release/bioc/html/sva.html
رمز التحليل هذه الدراسة https://github.com/zqun1/circadian_immune_mouse_melanoma

توافر الموارد

جهة الاتصال الرئيسية

يجب توجيه المزيد من المعلومات وطلبات الموارد أو الكواشف إلى جهة الاتصال الرئيسية، كريستوف شيرمان.christoph.scheiermann@unige.ch).

توفر المواد

لم تنتج هذه الدراسة مواد كيميائية فريدة جديدة.

توفر البيانات والشيفرة

  • تم إيداع بيانات تسلسل RNA أحادي الخلية في GEO: GSE260641 وهي متاحة للجمهور اعتبارًا من تاريخ النشر. أرقام الوصول للبيانات العامة المستخدمة في هذه الدراسة مدرجة في جدول الموارد الرئيسية. جميع البيانات التي تدعم استنتاجات هذه الورقة متاحة على الإنترنت.https://doi.org/10.26037/yareta:5rmtlmd5sreo7bmcbe44bbijfm).
  • تم إيداع جميع الشيفرات الأصلية في GitHubhttps://github.com/zqun1/circadian_immune_mouse_melanomaوهو متاح للجمهور اعتبارًا من تاريخ النشر.
  • أي معلومات إضافية مطلوبة لإعادة تحليل البيانات المبلغ عنها في هذه الورقة متاحة من جهة الاتصال الرئيسية عند الطلب.

تفاصيل النموذج التجريبي وشارك الدراسة

الحيوانات

تم شراء فئران C57BL/6 و NSG من تشارلز ريفر أو مختبر جاكسون. خطوط الفئران المتحولة Bmal1 (تم شراؤه من مختبرات جاكسون) و Cdh5 (هدية من الدكتور رالف آدامز، معهد ماكس بلانك للطب الحيوي الجزيئي في مونستر، ألمانيا) تم تهجينها وتربيتها في ENVIGO. تم الحفاظ على الفئران المعدلة وراثيًا كمتجانسة لجين Bmal وheterozygous للـ Cre المعني. تم تربية فئران CD45.1 OT-I (هدية من والتر رايث) في المنزل. كانت جميع الفئران المستخدمة تتراوح أعمارها بين 8-12 أسبوعًا. إعادة التركيب بواسطة Cre في Cdh5 تم تحفيز الفئران عن طريق الحقن داخل البطن بالتاموكسيفين (1 ملغ لكل حقنة) لمدة 5 أيام متتالية. ما لم يُذكر خلاف ذلك، كانت الفئران تعيش في بيئة تحت جدول ضوء:ظلام مع الطعام والماء حسب الرغبة. عندما تم التحقيق في نقاط زمنية متعددة في وقت واحد، تم استخدام خزائن محكمة الإغلاق (Techniplast) لنقل الحيوانات إلى المرحلة المعنية قبل التجارب. تتوافق أوقات العلاج مع وقت زايجبير (ZT) وتشير إلى التوقيت بالنسبة لتشغيل الأضواء في منشأة الحيوانات بحيث يكون ZT1 بعد ساعة من تشغيل الأضواء (صباحًا)، و ZT13 بعد ساعة من إطفاء الأضواء (مساءً). تم تحفيز نموذج ورم عفوي في BRaf. بتن فئران Tyr::CreERT2 كما تم وصفه سابقًا عن طريق العلاج الموضعي بـ -هيدروكسي تاموكسيفين ( على سطح الجلد. خلال المرحلة المظلمة، تم إجراء جميع المناورات تحت ضوء أحمر خافت. تم الموافقة على جميع إجراءات وتجارب الحيوانات وتم تنفيذها وفقًا لإرشادات لجنة أبحاث الحيوانات في جنيف، سويسرا، والسلطة البيطرية في كانتون فود، سويسرا، أو لجنة رعاية واستخدام الحيوانات في MGH (IACUC) وتم تنفيذها وفقًا للوائح IACUC في MGH.

تفاصيل الطريقة

خطوط خلايا الورم والتطعيم

تم الحفاظ على خطوط خلايا الميلانوما الفأرية B16-F10 (ATCC) وB16-F10-OVA (هدية من ستيفاني هيوجز، جامعة جنيف، سويسرا) في وسط RPMI (جيبكو) معززة بـ مصل فتيان البقر المعطل حرارياً (جيبكو) البنسلين-ستربتوميسين (جيبكو)، و من -ميركابتوإيثانول (جيبكو). تم الحفاظ على خلايا سرطان القولون الغدي MC38 (هدية من ميكائيل بيتيت، جامعة جنيف، سويسرا) في DMEM (جيبكو)، مصل فسيولوجي معطل حرارياً البنسلين-ستربتوميسين). تم استخدام خطوط الخلايا حتى المرور 10 وتم اختبارها سلبية لميكوبلازما. ما لم يُذكر خلاف ذلك، خلايا الورم في تم حقن PBS تحت الجلد في الجانب الأيمن من الفئران، تحت تخدير الإيزوفلوران. لتجربة علاج مضاد PD-1، خلايا الورم في PBS أو في تم حقن PBS تحت الجلد. خط خلايا اللمفوما البشرية DoHH2 (هدية من فرانشيسكو بيرتوني، المعهد)
تم زراعة (IOR) أبحاث الأورام، بيلينزونا، سويسرا في RPMI (جيبكو) معززة بـ FCS المعطل حرارياً (جيبكو). إجمالي تم حقن الخلايا تحت الجلد في الجانب الأيمن من فئران NSG. تم مراقبة حجم الورم كل 1 إلى 2 يوم باستخدام كاليبر وحُسب بواسطة الطول. عرض عرض .

تنشيط خلايا T في الفئران

خلايا الطحال من فئران C57BL/6 من النوع البري، Cd4cre:Bmal1 وتم جمع الفئران المتحكم بها من نفس السلالة، أو فئران OT-I، عن طريق سحق الطحال من خلال منخل خلايا شبكي. تم تحلل كريات الدم الحمراء باستخدام محلول تحلل كريات الدم الحمراء (Biolegend، 420302) لمدة 5 دقائق على الثلج. لتفعيل خلايا من C57BL/6 أو Cd4cre:Bmal1 فئران تم إثراء الخلايا باستخدام مجموعة عزل الخلايا (ميلتينيي) وفقًا لتعليمات الشركة المصنعة. تم عد خلايا T المعزولة وإعادة تعليقها بتركيز من في وسط خلايا T الكامل (RPMI، 10% مصل فتيان معطل حرارياً، 2 مللي مولار L-glutamine، البنسلين-ستربتوميسين -ميركابتوإيثانول، 1 مللي مول من صوديوم البيروفات) مدعوم بـ إنترلوكين-2 (IL-2) من الفأر (Biolegend) في طبق 24 بئر. داينابيدز تم استخدام كريات Mouse T-Activator CD3/CD28 (Gibco، 11452D) لتنشيط خلايا T وفقًا لتعليمات الشركة المصنعة. لتنشيط خلايا OT-I، تم إعادة تعليق الخلايا الطحالية الكاملة بتركيز من في وسط خلايا T الكامل (RPMI، مصل فسيولوجي معطل حرارياً، 2 مللي مول من L-glutamine، 1% من البنسيلين-ستربتوميسين، -ميركابتوإيثانول، 1 مللي مول من صوديوم البيروفات) مدعوم بـ الإنترلوكين-2 (IL-2) من الفأر (Biolegend) و1 نانومتر من ببتيد OVA257-264 (SIINFEKL، InvivoGen) في طبق 6 آبار. بعد 3 أيام من التنشيط، تم جمع الخلايا وتنقيتها عن طريق إزالة دينا بيدز (لخلايا متعددة النسائل) أو باستخدام طقم عزل الخلايا (ميلتيني) (لـ OT-I). للتجارب داخل الجسم، تم حقن خلايا OT-I عن طريق الوريد (i.v) في الفئران الحاملة للأورام.

تزامن خلايا T في الفئران

تم تنسيق الخلايا باستخدام مصل الخيل كما هو موصوف سابقًا. باختصار، تم جمع خلايا T المنشطة (كما هو موضح أعلاه) وزرعها في صفيحة 24 بئر بكثافة في وسط خلايا T الكامل، معزز بـ إنترلوكين-2 (IL-2) للفأر (بايو ليجند) و الأجسام المضادة المضادة لـ CD3 الفأري (Invitrogen، 16-0031-85) للحفاظ على إشارة TCR. تم تسخين حجم متساوي من مصل الخيل (Sigma، h1270) مسبقًا وأضيف مباشرة إلى الصفيحة (صدمة المصل). بعد 2 ساعة من الحضانة في مع تم غسل الخلايا في CO2 وإعادة تعليقها في وسط كامل لخلايا T مع IL-2 والأجسام المضادة anti-CD3.

إنسان تنشيط خلايا T والتزامن

إنسان تم عزل خلايا T باستخدام مجموعة عزل خلايا T (Miltenyi 130-096-495) من الخلايا الوحيدة النواة في الدم المحيطي من طبقة البوفي (البوفي كوت) للمتبرعين الأصحاء (المقدمة من مستشفيات جامعة جنيف (HUG)، سويسرا) باستخدام فيكول-باك بلس (Cytiva). تم إعادة تعليق الخلايا عند وتم تنشيطها لمدة 3 أيام باستخدام Dynabeads T-Activator البشري CD3/CD28 (Gibco، 11132D) في وجود IL-2 البشري المؤتلف ( بيبرتيك، 200-02) في وسط الخلايا التائية الكامل عند مع تم تنسيق الخلايا باستخدام مصل الخيل كما هو موضح أعلاه، مع إضافة إنترلوكين-2 البشري أجسام مضادة مضادة للبشر CD3 (إنفيتروجين، النسخة OKT3، 16-0037-81).

تعبير لوكفيراز خلايا

تم عزل الخلايا من الطحال لفئران Per2::Luc وتم تنشيطها كما هو موضح أعلاه. تم تنشيط CD8 تم تنسيق الخلايا وإعادة تعليقها في في 2 مل من وسط خلايا T الكاملة المحتوية على لوسيفرين (Abcam، ab45164). تم إغلاق الأطباق باستخدام بارافيلم وتم تسجيلها باستخدام كاشف لومي سايكل (أكتيميتريكس).

هضم الأنسجة وتحضير الخلايا المفردة

تم جمع نسيج الورم وقطعه إلى قطع صغيرة، وتم إضافة وسط الهضم (RPMI يحتوي على كولاجيناز IV (شركة وورثينغتون للكيماويات الحيوية) DNase I (روش 04716728001) و مصل فتي (FCS) المعطل حرارياً وتم حضنه لمدة 30 دقيقة عند تم شطف الخلايا من خلال منخل الخلايا للحصول على تعليقات خلوية مفردة. لإزالة الحطام، تم استخدام Lympho-lyte(R)-M (CEDARLANE) وفقًا لتعليمات الشركة المصنعة.

تدفق الخلايا

تم إعداد تعليقات خلوية مفردة وتم تحضينها مع حاصرات مستقبلات Fc للفئران أو البشر (مضاد الفئران CD16/32 من Biolegend، مادة حجب مستقبلات Fc البشرية من Miltenyi Biotec) لمدة 10 دقائق في درجة حرارة الغرفة. بعد التحضين، ما لم يُذكر خلاف ذلك، تم إضافة مزيج الأجسام المضادة مباشرة إلى تعليق الخلايا وتم تحضينه لمدة 15 دقيقة في .
تم غسل الخلايا وإعادة تعليقها في محلول عازل FACs مع صبغة حيوية (DRAQ7، بايوليجند، ) وتمت ميزتها باستخدام جهاز BD LSR Fortessa ذو 18 لونًا (BD Biosciences) أو Beckman Coulter Cytoflex. تم تحليل البيانات المكتسبة باستخدام FACSDiva 6 (BD Biosciences) وFlowJo 10 (BD). تم حساب عدد الخلايا باستخدام كرات العد (C36950، C36995، ThermoFisher).
لصبغ الخلايا داخل الخلايا، تم صبغ الخلايا بصبغة الحياة (eBioscience صبغة القابلية للإصلاح eFluor 780، 65-0865-18)، تليها تلون السطح كما تم وصفه سابقًا. تم تثبيت الخلايا وتمريرها باستخدام مجموعة محلول صبغ عوامل النسخ Foxp3 (eBioscience، 00-5523-00). بعد الغسل بمحلول التمرير، تمت إضافة الأجسام المضادة داخل الخلايا وتمت المعالجة لمدة 30 دقيقة في درجة حرارة الغرفة. تم استخدام الأجسام المضادة التالية لصبغ داخل الخلايا: مضاد غرانزيم الفأر ومضاد IFN- الفأر. “. لاكتشاف السيتوكينات (Granzyme B، IFN- )، تم تحفيز الخلايا قبل التلوين في وسط RPMI الكامل مع فوروبيول 12-ميريستات 13-أسيتات (PMA، ; سيغما-ألدريتش)، أيونوميسين ( ; سيغما-ألدريتش)، وبي دي غولجي بلوق ( ; BD
العلوم الحيوية، 555029) وتم حضنها لمدة 4 ساعات في ثاني أكسيد الكربون. كاسبيز- تم إجراء الكشف باستخدام CellEvent كاسبيز-3/7 الأخضر جاهز للتدفق المادة الكيميائية (Invitrogen R37167) وفقًا لتعليمات الشركة المصنعة.

نقل كريات الدم البيضاء من الفئران

تم جمع خلايا المتبرعين من نخاع العظم (تم شطفه بمحلول PBS) والطحال (عن طريق سحقه من خلال) مصفاة الخلايا) من الفئران المتبرعة وتم تحللها باستخدام محلول تحلل كريات الدم الحمراء (420302 بايو ليجند). تم وسم الخلايا بـ CellTracker صبغة حمراء داكنة في PBS، )، وتم حقنه عن طريق الوريد ( خلايا نخاع العظام و الخلايا الطحالية). ما لم يُذكر خلاف ذلك، تم جمع الأنسجة بعد ساعتين من حقن الخلايا. بالنسبة لاختبارات الخروج، كان المجموع الكلي لـ خلايا ( خلايا نخاع العظام و الخلايا الطحالية) في تم حقن PBS داخل الورم. تم جمع الأنسجة بعد 24 ساعة من الحقن.

اختبار EdU

EdU (1 ملغ لكل فأر، Click-iT بلاس إدو أليكسافلور تم حقن مجموعة اختبار تدفق الخلايا، إنفيتروجين، C10646) عن طريق البطن في الفئران قبل 4 ساعات من جمع الأنسجة. تم جمع أنسجة الورم وهضمها كما هو موضح أعلاه. تم صبغ الخلايا أولاً بصبغة الحياة (eBioscience صبغة القابلية للإصلاح eFluor 780، 65-0865-18)، تليها تلطيخ السطح كما تم وصفه سابقًا. ثم تم تثبيت الخلايا، وتم اختراقها، وتم الكشف عن EdU وفقًا لتعليمات الشركة المصنعة.

فرز خلايا T CD8 واستخراج RNA

للحصول على تم هضم أنسجة الورم ومعالجتها كما هو موضح أعلاه. ثم تم صبغ الكريات البيضاء المتسللة إلى الورم باستخدام مضاد CD45 الفأري (النسخة 30-F11 BUV 395)، CD3e (النسخة KT3.1.1، BV421)، CD4 (النسخة GK1.5 BV650)، CD8a (النسخة 53-6.7 BV605)، CD19 (النسخة 1D3، BB700)، وNK1.1 (النسخة PK136، PE/Cy5) مع DRAQ7 وCD45. تم فرز الخلايا باستخدام جهاز BD Aria II. تم فرز CD8 باستخدام تحليل تدفق الخلايا. تم جمع خلايا T في محلول RNAprotect Cell Reagent (رقم الكاتالوج #76526، كيوجين). تم عزل RNA باستخدام مجموعة RNeasy Plus Micro (رقم الكاتالوج #74034، كيوجين) وفقًا لتعليمات الشركة المصنعة. تم تقييم سلامة وكمية RNA باستخدام جهاز Bioanalyzer (أجيلة تكنولوجيز) أو جهاز Nanodrop 2000 (ثيرمو فيشر).

استخراج RNA، النسخ العكسي و qPCR

تم جمع الخلايا وتحليلها باستخدام مادة Trizol. تم تحليل كمية وجودة RNA باستخدام جهاز Nanodrop 2000 (ThermoFisher) أو Bioanalyzer. تم إجراء النسخ العكسي باستخدام PrimeScript. مجموعة كواشف RT (تاكارا) وفقًا للتعليمات المقدمة. تم إجراء تحليلات Q-PCR باستخدام PowerUp SYBR Green (أبلايد بيوسيستمز). تم إجراء تقدير النسخة باستخدام طريقة باستخدام Rplp0 أو Rpl32 كجينات مرجعية داخلية.

علاجات الأجسام المضادة في الجسم الحي

لمنع LFA-1، مضاد CD11a (نسخة M17/4، 100 ) و anti-CD18 (النسخة M18/2، تم إعطاء الأجسام المضادة عن طريق الوريد (i.v.). جسم مضاد مضاد لـ PD-1 (مضاد الفأر PD-1، النسخة تم إعطاؤه عن طريق الحقن داخل الصفاق (i.p) في اليوم الرابع أو السادس (B16-F10OVA) أو اليوم السادس أو السابع (MC-38) بعد زراعة الورم. لتقليل CD8 خلايا T، جسم مضاد للإزالة (مضاد للفأر CD8a، النسخة YTS 169.4، تم حقنها عن طريق البطن في نفس يوم إعطاء PD-1 وتكرارها كل يومين. تم استخدام الأجسام المضادة للتحكم من نوع IgG2a (النسخة 2A3) وIgG2b (النسخة LTF-2). تم شراء جميع الأجسام المضادة من BioXCell.

إنتاج خلايا CAR T البشرية

تم جمع خلايا الدم المحيطية الوحيدة النواة (PBMCs) من طبقة البوفي (البوفي كوت) للمتبرعين الأصحاء (المقدمة من مستشفيات جامعة جنيف (HUG)، سويسرا) باستخدام فيكول-باك بلس (Cytiva)، ثم تم تجميدها في جرعات وتخزينها في النيتروجين السائل. تم إذابة PBMCs المجمدة وإعادة تعليقها عند خلايا لكل مل، تم تنشيطها لمدة 3 أيام باستخدام Dynabeads T-Activator CD3/CD28 (Gibco، 11132D) في وجود IL-2 البشري المؤتلف ، بيبروتك، 200-02) في وسط RPMI1640 (جيبكو) معزّز بـ مصل فتي (FBS) المعطل حرارياً (جيبكو) و البنسلين-ستربتوميسين عند مع تم نقل خلايا T المنشطة عن طريق زراعتها لمدة 48 ساعة في أطباق 24 بئر مغطاة بالريترونيكتين (تاكارا، T100B) مع فيروس رجعي يحمل جين CAR من الجيل الثاني المضاد لـ CD19-CD28-CD3. تم توفير بناء CAR بلطف من قبل الدكتورة كريستال ماكال (جامعة ستانفورد). تمت إزالة دينابيدز، وبعد الغسل، تم زراعة الخلايا في وسط كامل مع IL-2. في اليوم الرابع عشر، تم جمع الخلايا وتجميدها حتى الاستخدام لاحقًا. تم تقييم فعالية نقل تعبير CAR على خلايا T من خلال ارتباط ستربتافيدين المرتبط بـ APC (Biolegend، 405207) بالبروتين البيوتيني. (ThermoFisher، 29997). قبل الاستخدام، CAR تم إذابة الخلايا، وغسلها، وتركها للراحة لمدة 4 ساعات في .

نقل خلايا CAR T البشرية

تم زراعة خلايا ليمفوما DoHH2 البشرية كما هو موضح أعلاه. من أجل تجارب العلاج على المدى الطويل، تم حقن خلايا CAR T أو عدد متساوٍ من الخلايا غير المعدلة وراثيًا عن طريق الوريد في الفئران الحاملة للأورام في ZT1 أو ZT13 في اليوم الثاني عشر بعد زراعة الورم. من أجل تجربة التوجه القصيرة الأمد، تم حقن خلايا CAR T عن طريق الوريد في الساعة ZT1 أو ZT13 في اليوم الخامس عشر بعد زراعة الورم. بعد أربع وعشرين ساعة من النقل، تم جمع الأورام وهضمها باستخدام وسط الهضم (RPMI يحتوي على كولاجيناز IV (شركة وورثينغتون للكيماويات الحيوية) كولاجيناز D (روش، 11088866001)، DNase I (روش 04716728001 ) و مصل فتي (FCS) المعطل حرارياً وتم حضنه لمدة 30 دقيقة في تم تحليل تعليق الخلايا المفردة كما هو موضح أعلاه.

تصوير المناعة الفلورية

تم جمع أنسجة الورم وتضمينها في كتل OCT (CellPath) وتم الاحتفاظ بها في تم قطع الأورام إلى شرائح باستخدام جهاز التبريد. تم تثبيت الشرائح بعد ذلك بـ يرجى العثور على 10 دقائق في درجة حرارة الغرفة. بعد ثلاث غسلات مع PBS، تم تحضينها مع BlockAid محلول الحجب (إنفيتروجن ، B10710) لمدة ساعتين عند درجة حرارة الغرفة. بعد PBS، تم صبغ الشرائح بمضاد أجسام مخفف في محلول الحجب في بين عشية وضحاها. تم استخدام الأجسام المضادة التالية: AF594 مضاد للفأر CD31 (نسخة MEC13.3، Biolegend 102520)، FITC مضاد للفأر CD4 (نسخة GK1.5، Biolegend 100405)، AF647 مضاد للفأر CD54 (نسخة YN1/1.7.4، Biolegend 116114)، AF647 مضاد للفأر CD8 (نسخة 53-6.7، Biolegend 100724)، AF488 مضاد للفأر VCAM-1 (نسخة 429، Biolegend 105710)، E-selectin (نسخة UZ6، Invitrogen MA1-06506)، مضاد ثانوي ماعز مضاد للفأر IgG (H+L) معاد امتصاصه، أليكسا فلور. (Invitrogen A21247). تم غسل الشرائح 3 مرات باستخدام PBS وتم صبغها بـ DAPI (300 نانومتر، Biolegend 422801) لمدة 5 دقائق في درجة حرارة الغرفة. تم الحصول على الصور كشرائح باستخدام ميكروسكوب زيس أكسيو إكزامينر.Z1 المجهز بمصادر ليزر 405 و488 و561 و640 نانومتر. تم إجراء جميع تحليلات الصور في ImageJ وSlidebook (ابتكارات التصوير الذكي، 3i).

التشخيص المناعي النسيجي البشري

تم صبغ مقاطع الميلانوما المثبتة بالفورمالين والمضمنة في البارافين (FFPE) بواسطة الكيمياء المناعية باستخدام الأجسام المضادة الأولية المضادة لـ CD4 (EP204) أو المضادة لـ CD8 (SP16) (جسر الذهب زونغشان للتكنولوجيا الحيوية، بكين، الصين). تم إجراء جميع الصبغات على منصة ليكا بوند III الآلية (ليكا بيوسيستمز، نيوكاسل، المملكة المتحدة) باستخدام مجموعة كشف بوند بوليمر ريفين ريد (DS9390). تم مسح الشرائح باستخدام ماسح الشرائح بانوراميك P250 FLASH (3DHistech، بودابست، المجر).

تسلسل الخلايا الفردية والتحليل

تم جمع الأورام في أوقات مختلفة من اليوم ومعالجتها كما هو موضح أعلاه. تم فرز الكريات البيضاء المتسللة إلى الورم (CD45 DRAQ7-) باستخدام جهاز فرز Biorad S3 في PBS مع تم تحميل تعليق الخلايا على أداة الجينوميات كروميم. تم إعداد مكتبات RNA-Seq أحادية الخلية باستخدام مجموعة كواشف كروميم أحادية الخلية 3′ v3.1 مع فهارس مزدوجة وفقًا لبروتوكول الشركة المصنعة. تم إجراء تقدير كمية المكتبة وتقييم الجودة باستخدام جهاز قياس الفلورية Qubit (ثيرمو فيشر ساينتيفيك) وجهاز Tapestation (شريحة حساسية عالية للحمض النووي – أجيلنت تكنولوجيز). تم تسلسل المكتبات على جهاز إيلومينا نوفا سيك 6000 باستخدام طريقة النهاية المزدوجة. تم استخدام وضع التسلسل bp. في المجموع، تم تسلسل 30693 خلية مفردة (ZT1: 9764؛ ZT7: 7933؛ ZT13: 6255؛ ZT19: 6651).
تم إنشاء مصفوفات العد الخام باستخدام Cell Ranger v6.1.2 مع بناء الجينوم الفأري mm 10 كمرجع. تم إزالة الجينات المعبر عنها في أقل من 20 خلية. تم استبعاد الخلايا التي تعبر عن أقل من 200 ميزة، أو أكثر من 7500 ميزة، أو أكثر من تمت إزالة الجينات الميتوكوندرية من التحليل الإضافي. تم إجراء تقليل الأبعاد وتحليل التجميع باستخدام Seurat v4. تم توضيح الخلايا بواسطة علامات الخلايا التقليدية: B (Cd79a، Igkc)، cDC1 (Clec9a، Xcr1)، cDC2 (H2-Aa، H2-Ab1، Itgax، Clec10a)، cDC المنشطة (Ccr7، Fscn1)، pDC (Siglech، Bst2)، البلعميات (Cd68، Itgam)، أحادية النواة (Ly6c2، Ccr2)، العدلات (G0s2، Csf3r)، NK (Nkg7، Klrb1c)، (Trdc, Tcrg-C1)، خلايا T الظهارية الشجرية (Trdv4، Tcrg-C1)، CD4 (Trac، Cd4)، CD8 (Trac، Cd8a)، Treg (Foxp3، Ctla4)، الخلايا الميلانية (Mlana، Dct)، والخلايا الانقسامية (Birc5، Top2a). تشمل ‘الخلايا اللمفاوية المضادة للورم’ NK، تشمل ‘الخلايا اللمفاوية المضادة للورم’ مجموعات NKT و CD8 الفعالة و CD4 الساذجة و CD8 ذاكرة فعالة و CD8 انقسامية و CD4 الفعالة، بينما تشمل ‘الخلايا اللمفاوية المؤيدة للورم’ مجموعات Treg و Treg انقسامية و CD8 المستنفدة. تم تصنيف خلايا تحت المجموعات CD4 و CD8 بناءً على علامات الاستنفاد (Tox و Lag3 و Pdcd1 و Havcr2) و الفعالية (Gzmb و Prf1 و Ifng و II2 و II4) و الذاكرة (Gzmk و Itga4) و الذاكرة/ السذاجة (Tcf7 و Sell و Ccr7 و Lef1). تم تحديد الجينات المتذبذبة بواسطة Discorhythm. مع اعتبار كل خلية كقراءة واحدة. تم إجراء تحليلات تمثيل مجموعة الجينات الزائدة باستخدام clusterProfiler. تم حساب درجات توقيع الجينات بواسطة GSVA. تم إجراء تحليلات البقاء باستخدام حزمة TCGAbiolinks. تم إجراء ملاءمة كوزينور للتعبير الجيني باستخدام حزمة كوزينور.

تحليل ZeitZeiger

ساعة الزمن تم تدريب النموذج باستخدام بيانات تسلسل RNA أحادي الخلية لورم الميلانوما في الفئران من الدراسة الحالية، وبيانات تسلسل RNA أحادي الخلية لجلد الفئران من داخل المؤسسة، وبيانات تسلسل RNA أحادي الخلية لجلد الفئران من GSE223109، وبيانات تسلسل RNA الكمي من جلد الفأر من GSE115104 و GSE83855 مجموعات البيانات. لتحقيق التوازن في تمثيل بيانات scRNAseq مقابل RNAseq وبيانات الميلانوما مقابل بيانات الجلد السليمة، تم أخذ عينة عشوائية من مجموعة بيانات scRNAseq الخاصة بالميلانوما بنسبة 60% أربع مرات، وتم أخذ عينة عشوائية من مجموعة بيانات scRNAseq الخاصة بالجلد في المختبر بنسبة مرتين لتشكيل مجموعة التدريب. بالنسبة لمجموعات بيانات scRNAseq، تم حساب التعبير الزائف الكثيف في كل فأر وضربه في 100 لتوليد النسخ لكل مليون كما في البيانات لتدريب ZeitZeiger. تم تقسيم كل مجموعة بيانات عشوائيًا إلى مجموعة تدريب ومجموعة اختبار بنسبة تم تصحيح تأثير الدفعة لبيانات التدريب وبيانات الاختبار بواسطة Combat من حزمة SVA. تم تلخيص المعاملات المحفوظة من 20 جولة من التحقق المتقاطع مع بذور عشوائية مختلفة ومعلمات sumabsv و nSpc المختلفة. تم إضافة الجينات KIf9 و Nr1d1 و Cry2 و Per1 و Gapdh و Ppia و Psmb2 يدويًا إلى المعاملات، كما تم وصفه سابقًا. تم إجراء مزيد من التحقق من صحة نموذج ZeitZeiger المدرب باستخدام بيانات RNAseq أو بيانات المصفوفة الدقيقة المعيرة بواسطة RMA من كل من الفأر (GSE114943، جي إس إي 174155 جي إس إي 38622 جي إس إي 38623 ) والإنسان (GSE35635، جي إس إي 205155 جي إس إي 112660 و GSE139305 لتصحيح الفرق السلوكي بين الإنسان والفأر خلال النهار والليل، تم توليد أوقات أخذ العينات المتوقعة للإنسان عن طريق إضافة 12 ساعة.

تقدير تركيبة الخلايا في مرضى الميلانوما في TCGA

لتقدير نسبة المستهلكين مقابل غير المستهلكين الخلايا، بالإضافة إلى المجموعات المناعية الرئيسية، خلايا من GSE120575 تمت الإشارة إليها كما هو موضح سابقًا وتم استخدام أفضل 20 جينًا معبرًا بشكل مختلف لتقدير تركيبة خلايا المناعة بواسطة cibersortx.

التكميم والتحليل الإحصائي

ما لم يُذكر خلاف ذلك، تم رسم جميع البيانات من تكرارات بيولوجية مستقلة. تم تحليل البيانات باستخدام برنامج بريزم 9 وبريزم 10 (غرافيك باد). .

الرسوم التوضيحية الإضافية

الشكل S1. وقت الحصاد يحدد أعداد الكريات البيضاء المتسللة إلى الورم، مرتبط بالشكل 1
استراتيجية البوابة لكريات الدم البيضاء في B16-F10-OVA و Tyr::Cre ; بتين ; براف أورام.
(ب) الأعداد الكلية الطبيعية لمجموعات الكريات البيضاء في أورام B16-F10-OVA، التي تم جمعها في 4 أوقات مختلفة من اليوم (وقت الزايتgeber [ZT])؛ فئران من 4 تجارب مستقلة، تحليل كوزينور.
(ج) أعداد كريات الدم البيضاء المأخوذة في 4 أوقات مختلفة من اليوم من فئران تحمل ورم B16-F10-OVA، فئران من 5 تجارب مستقلة، تحليل كوزينور.
حجم ورم B16-F10-OVA، الذي تم قياسه 3 مرات في اليوم، فئران.
(هـ) قياس CD31 منطقة في أورام B16-F10-OVA، التي تم جمعها في ZT1 أو ZT13، الفئران، تحليل التباين الأحادي، واختبار الطالب غير المتزاوج اختبار.
(F) الأعداد المعيارية لكريات الدم البيضاء المأخوذة عند ZT1 أو ZT13 من Tyr::Cre ; بتين ; BRaf فئران تحمل الأورام فئران من تجربتين مستقلتين، اختبار t لطلبة المزدوج. جميع البيانات ممثلة كمتوسط SEM، غير دال إحصائيًا. *** ; **** .
الشكل S2. تسلل الكريات البيضاء الدائري وعلاج الخلايا التائية، متعلق بالشكل 2
(أ) استراتيجية البوابة وقياس EdU الكريات البيضاء في أورام B16-F10-OVA، فئران، اختبار t لطلبة غير المقترنين.
(ب) استراتيجية البوابة وقياس الكاسبيز-3/-7 في كريات الدم البيضاء في أورام B16-F10-OVA، التي تم جمعها في ZT1 ( الفئران) أو ZT13 ( الفئران)، من تجربتين مستقلتين، اختبار ستودنت غير المقترن اختبار.
(ج) الأعداد المعيارية للخلايا المتبقية في الأورام بعد 24 ساعة من الحقن داخل الورم الكريات البيضاء عند ZT1 الفئران) أو ZT13 ( الفئران)، من تجربتين مستقلتين، اختبار ستودنت غير المتزاوج اختبار.
حجم الورم بعد علاج مضاد LFA-1 ( فئران)، التحكم في النوع المتماثل ( الفئران)، من تجربتين مستقلتين، تحليل التباين ثنائي الاتجاه.
(E) قياس VCAM-1 ( أقسام من 3 فئران) أو E-selectin ( تعبير الأقسام من 4 فئران على CD31 الخلايا في أورام B16-F10-OVA، التي تم جمعها في ZT1 أو ZT13، من تجربتين مستقلتين، اختبار t لطلاب غير المتزاوج. جميع البيانات ممثلة كمتوسط SEM، غير دال إحصائيًا.
حجم الورم بعد الحقن الوريدي بـ خلايا OT-I المفعلة عند ZT1 الفئران) أو ZT13 ( الفئران)، من 3 تجارب مستقلة.
الشكل S3. تحديد مجموعات TIL، المتعلقة بالشكل 3
(A-C) رسم ميزات وخريطة حرارية للجينات لتحديد كل مجموعة من كريات الدم البيضاء التي تم الحصول عليها من تسلسل RNA أحادي الخلية لأورام B16-F10-OVA، التي تم جمعها في 4 أوقات مختلفة من اليوم.
خلية
مقالة
الوصول المفتوح
الشكل S4. تذبذبات TILs اليومية، المتعلقة بالشكل 3
تمثيل UMAP لتحليلات scRNA-seq لـ TILs التي تم جمعها في 4 أوقات مختلفة من اليوم.
(ب) الوفرة النسبية لخلايا المناعة لكل مجموعة عبر أوقات مختلفة من اليوم.
(ج) تعبير الجين الدائري Per1 في كل مجموعة من كريات الدم البيضاء التي تم الحصول عليها من تسلسل RNA أحادي الخلية (scRNA-seq) لأورام B16-F10-OVA، التي تم جمعها في 4 أوقات مختلفة من اليوم.
(د) عدد الجينات المتذبذبة بشكل ملحوظ التي تم اكتشافها في كل مجموعة باستخدام تحليلات كوزينور (CS) ودورة JTK.
(E) المسارات الجينية المتذبذبة الشائعة في مجموعات الخلايا المناعية المختلفة.
خريطة حرارية للجينات المتذبذبة لكل مجموعة. تم تمييز الجينات الدائرية.
الشكل S5. الجينات المتذبذبة يوميًا في TILs، المتعلقة بالشكل 3
رسم ميزات (A و B) ورسم نقطي للجينات لتحديد مجموعات خلايا T المختلفة التي تم الحصول عليها من تسلسل RNA أحادي الخلية لأورام B16-F10-OVA، التي تم جمعها في 4 أوقات مختلفة من اليوم.
(ج) الوفرة النسبية لكل تجمع الخلايا في أوقات مختلفة من اليوم.
(د) استراتيجية البوابة للمستنفدين وغير المستنفدين الخلايا، محددة مسبقًا على خلايا CD8 T.
(E) تعبير Pdcd1 تم قياسه من خلال تحليل pseudobulk لخلايا CD8 العنقود الذي تم الحصول عليه من تسلسل RNA أحادي الخلية لأورام B16-F10-OVA، التي تم جمعها في أربع أوقات مختلفة من اليوم. لتقليل تشتت بيانات تسلسل RNA أحادي الخلية، تم أخذ عينة عشوائية من العنقود بواسطة لمدة 10 مرات باستخدام بذور مختلفة. تم حساب متوسط مستويات التعبير لجين Pdcd1 في 4 أوقات مختلفة من اليوم داخل كل خلية زائفة.
الشكل S6. تذبذب نمط الخلية في الفئران والبشر، مرتبط بالأشكال 4 و5 و6
نظرة عامة تخطيطية لـ تزامن الخلايا في المختبر.
(ب) حجم الورم بعد علاج مضاد PD-1 المطبق في ZT1 أو ZT13 في أورام B16-F10-OVA؛ فئران؛ تحكم الفئران، من 3 تجارب مستقلة.

Journal: Cell, Volume: 187, Issue: 11
DOI: https://doi.org/10.1016/j.cell.2024.04.015
PMID: https://pubmed.ncbi.nlm.nih.gov/38723627
Publication Date: 2024-05-01

Cell

Circadian tumor infiltration and function of CD8 T cells dictate immunotherapy efficacy

Graphical abstract

Highlights

-The quantity and function of tumor-infiltrating T cells are time-of-day dependent
-Rhythmic tumor infiltration of T cells depends on the endothelial circadian clock
-Adjusting the time of injection during the day improves CAR T cell therapy efficacy
-Enhanced tumor control by timed anti-PD-1 therapy depends on circadian T cell activity

Authors

Chen Wang(王宸),Qun Zeng(曾群), Zeynep Melis Gül,...,Federico Simonetta, Mikaël J.Pittet,Christoph Scheiermann

Correspondence

chen.wang.1@unige.ch(C.W.), christoph.scheiermann@unige.ch(C.S.)

In brief

Efficacy of CAR T cell therapy and immune checkpoint blockade in cancer can be improved by adjusting the time of treatment during the day due to rhythmic infiltration and function of T cells regulated by the circadian clock in T cells and endothelial cells in the tumor microenvironment.

Circadian tumor infiltration and function of cells dictate immunotherapy efficacy

Chen Wang(王宸), Qun Zeng(曾群), Zeynep Melis Gül, Sisi Wang, Robert Pick, Phil Cheng, Ruben Bill, Yan Wu, Stefan Naulaerts, Coline Barnoud, Pei-Chun Hsueh, Sofie Hedlund Moller, Mara Cenerenti, Mengzhu Sun, Ziyang Su, Stéphane Jemelin, Volodymyr Petrenko, Charna Dibner, Stéphanie Hugues, Camilla Jandus, Zhongwu Li, Olivier Michielin, Ping-Chih Ho, Abhishek D.Garg, Federico Simonetta, Mikaël J.Pittet, and Christoph Scheiermann Department of Pathology and Immunology,Faculty of Medicine,University of Geneva, 1211 Geneva,Switzerland Translational Research Centre in Onco-Hematology(CRTOH),Geneva 1211,Switzerland Department of Oncology and Precision Oncology Service,Geneva University Hospitals,University of Geneva,Geneva 1211,Switzerland AGORA Cancer Research Center,Lausanne 1011,Switzerland Center for Systems Biology,Massachusetts General Hospital Research Institute and Harvard Medical School,Boston,MA 02114,USA Key Laboratory of Carcinogenesis and Translational Research,Department of Pathology,Peking University Cancer Hospital &Institute, Beijing 100142,China Laboratory of Cell Stress &Immunity,Department of Cellular &Molecular Medicine,KU Leuven,Leuven 3000,Belgium Department of Fundamental Oncology,University of Lausanne,Lausanne 1066,Switzerland Ludwig Institute for Cancer Research,Lausanne 1005,Switzerland Division of Thoracic and Endocrine Surgery,Department of Surgery,Geneva University Hospitals,Geneva 1205,Switzerland Department of Cell Physiology and Metabolism,Faculty of Medicine,University of Geneva,Geneva 1211,Switzerland Diabetes Center,Faculty of Medicine,University of Geneva,Geneva 1211,Switzerland Institute of Genetics and Genomics of Geneva(iGE3),Geneva 1211,Switzerland Geneva Centre for Inflammation Research(GCIR),Geneva 1211,Switzerland Division of Hematology,Department of Oncology,Geneva University Hospitals,University of Geneva,Geneva,Switzerland Biomedical Center(BMC),Institute for Cardiovascular Physiology and Pathophysiology,Walter Brendel Center for Experimental Medicine (WBex),Faculty of Medicine,Ludwig-Maximilians-Universität(LMU)Munich,Planegg-Martinsried 82152,Germany These authors contributed equally Lead contact*Correspondence:chen.wang.1@unige.ch(C.W.),christoph.scheiermann@unige.ch(C.S.)https://doi.org/10.1016/j.cell.2024.04.015

Abstract

SUMMARY The quality and quantity of tumor-infiltrating lymphocytes,particularly cells,are important param- eters for the control of tumor growth and response to immunotherapy.Here,we show in murine and hu- man cancers that these parameters exhibit circadian oscillations,driven by both the endogenous circa- dian clock of leukocytes and rhythmic leukocyte infiltration,which depends on the circadian clock of endothelial cells in the tumor microenvironment.To harness these rhythms therapeutically,we demon- strate that efficacy of chimeric antigen receptor T cell therapy and immune checkpoint blockade can be improved by adjusting the time of treatment during the day.Furthermore,time-of-day-dependent T cell signatures in murine tumor models predict overall survival in patients with melanoma and correlate with response to anti-PD-1 therapy.Our data demonstrate the functional significance of circadian dy- namics in the tumor microenvironment and suggest the importance of leveraging these features for improving future clinical trial design and patient care.

INTRODUCTION

Immunotherapy is currently one of the most effective means of controlling tumors and,in some cases,curing patients. The most promising interventions involve the injection of chimeric an- tigen receptor(CAR)T cells or immune checkpoint inhibitors (also known as immune checkpoint blockers[ICBs]). CAR T cell therapy consists of administering to the patient T cells
that are genetically engineered to recognize tumor cells. ICB therapy aims to activate the antitumor functions of endogenous T cells;for example,the administration of an anti-programmed cell death protein 1 (PD-1)antibody inhibits suppressive mole- cules on the surface of T cells,thereby activating them.
Anti-tumor immunity has recently been demonstrated to be time-of-day dependent. Specifically,dendritic cells are more effective in controlling tumors grafted into mice at specific times
of the day. Overall functions of the immune system are under a circadian, control. Leukocytes emigrate from blood to tissues such as lymph nodes in a rhythmic fashion, which has been associated with better antigen recognition at specific times of the day. In addition to trafficking properties, the phenotype of immune cells is also time-of-day gated, so that stimulation at a specific time of the day yields different outcomes in the associated immune responses. This has been linked to time-of-day dependent reactions to vaccines in mice and humans and was recently extended also to anti-tumor vaccinations. Although the immune system is highly rhythmic in the steady state and after inflammatory challenge, it is unknown whether circadian rhythms regulate the distribution and phenotype of immune cells within tumors.
Considering that the immune system is under circadian control, it seems important to define the possible relevance of circadian rhythm in an immunotherapy context. Here, we set out to define whether the tumor microenvironment (TME) itself exhibits circadian oscillations. We demonstrate that tumor-infiltrating leukocytes (TILs) exhibit circadian oscillations in both number and phenotype, which can be exploited therapeutically with timed administration of CAR T cells or ICBs.

RESULTS

Leukocyte infiltration of tumors exhibits circadian oscillations

We explored whether the number of leukocytes in tumors exhibits time-of-day differences, by injecting B16-F10 melanoma cells expressing ovalbumin (B16-F10-OVA) subcutaneously (s.c.) into cohorts of mice at one time of the day (that is, at zeitgeber time 9 [ZT9; 9 h after light onset in a 12 h light/ 12 h dark environment]) and harvested tumors 12 days later at 4 different times of the day (ZT1 [“morning”], ZT7 [“day”], ZT13 [“evening”], and ZT19 [“night”]). The overall number of TILs was highly dependent on the time of day of tumor harvest, with total leukocyte numbers peaking in the evening (ZT13) (Figures 1A and S1A). Specifically, we detected time-of-day differences in the numbers of tumor and cells, cells, cells, as well as CD19 cells (phase-shifted from blood leukocyte counts), whereas CD11b cells and CD11b cells did not exhibit significant changes (Figures 1A, S1B, and S1C). We confirmed these data by immunofluorescence staining of tumor sections, showing higher T cell numbers in the evening (ZT13) compared with the morning (ZT1) (Figure 1B). We did not observe significant differences in tumor volume within a day (Figure S1D), and although tumors had longer times to grow when harvested at night (ZT19) compared with the evening (ZT13), we observed higher leukocyte numbers in the evening (Figure 1 A ). This demonstrated that differences in leukocyte numbers were due to the time of day of tumor harvest. We further observed no changes in vascular density in the tumor (Figure S1E), indicating this phenotype to be immune cell specific. We confirmed the phenotype in a spontaneous melanoma tumor model (Tyr::Cre ; Pten ; BRaf ) (Figures 1C, S1A, and S1F), which together demonstrated that numbers of TILs exhibit time-of-day differences, peaking in the evening.
We next transferred animals to complete darkness conditions, as circadian rhythms persist even under constant environmental conditions, such as in the absence of rhythmic lighting patterns. Constant darkness did not alter the observed time-of-day differences, demonstrating these oscillations to be bona fide circadian in nature and not simply a response to a rhythmic light:dark environment (Figure 1D). However, switching mice from a normal light-dark to a 12 h inverted dark-light cycle inversed TIL numbers (Figures 1E and 1F). This demonstrated that the oscillations were not dependent on light per se but that they could be synchronized to a new lighting regime. This ability to be synchronized by an external cue is an additional feature of circadian rhythms (Figure 1F). Subjecting mice to an acute jet-lag protocol, which consisted of shifting mice to a 6 h phase-delay and harvesting tissues at the previously assessed respective morning and evening time points thereafter (Figure 1E), furthermore abrogated time-of-day differences in TILs (Figure 1F). Together, these data provide unexpected evidence that TILs exhibit circadian oscillations in cellularity, a process that is entrained by light.

Endothelial cells gate circadian leukocyte infiltration

To assess the etiology of these oscillations, we quantified whether leukocyte proliferation, death, egress, and/or infiltration exhibited time-of-day differences. We did not observe significant differences in total leukocyte proliferation, apoptosis, or egress (Figures S2A-S2C). By contrast, an acute (24 h) block of leukocyte infiltration with an anti-lymphocyte function-associated antigen 1 (LFA-1)-integrin (anti- ) antibody greatly reduced leukocyte counts in the tumor in the evening (ZT13) and abrogated the time-of-day differences (Figure 2A). This suggested that leukocyte input from the blood accounted for the time-ofday differences and was required to maintain overall TIL cellularity. It further indicated that leukocyte supply to the tumor from the blood was surprisingly dynamic over such a short time frame. Leukocyte infiltration into the tumor from blood was of functional relevance since chronic anti-LFA-1-integrin treatment strongly increased tumor volume (Figure S2D). To investigate whether leukocytes indeed infiltrated the tumor in a time-of-day-dependent manner, we performed short-term (2 h) homing assays with adoptively transferred leukocytes. Cells were harvested from donor mice at one time of the day and intravenously (i.v.) injected into phase-shifted recipients, housed in cabinets with inverted lighting regimes, thus limiting the variable to timing in the recipient only. Leukocyte trafficking to the tumor was strongly time-of-day dependent, with more cells reaching the tumor in the evening, in line with the infiltration block data (Figure 2B). This demonstrated that, in mice, the TME exhibited pro-migratory properties in the evening, which were not driven by differences in vascular density (Figure S1E). Since the -integrin mediated leukocyte infiltration (Figure 2A), we quantified expression levels of its ligands on endothelial cells of the tumor vasculature, as gatekeepers located at the blood-tumor interphase. Intercellular adhesion molecule 1 (ICAM-1) was much higher expressed by endothelial cells in the evening (ZT13) compared with the morning (ZT1) (Figure 2C), whereas vascular cell adhesion molecule 1 (VCAM-1) as well as E-selectin, additional endothelial cell adhesion molecules, did not exhibit time-of-day differences (Figure S2E). This phenotype was still present
Figure 1. Time of day of harvest dictates numbers of tumor-infiltrating leukocytes
(A) Normalized total numbers of tumor-infiltrating leukocytes in a B16-F10-OVA tumor, harvested at 4 different times of the day (zeitgeber time [ZT]); , 12 mice from 4 independent experiments, cosinor analysis. Shaded areas indicate dark phases.
(B) Imaging (left) and quantification (right) of CD4 and CD8 T cells in B16-F10-OVA tumors harvested at ZT1 or ZT13; mice, unpaired Student’s test, scale bars: .
(C) Numbers of tumor-infiltrating leukocytes in Tyr::CreERT2, BRaf , Pten mice harvested at ZT1 or ZT13; mice, unpaired Student’s t test.
(D) Normalized total numbers of tumor-infiltrating leukocytes in a B16-F10-OVA tumor, harvested at different times of the day under constant darkness conditions (circadian time [CT]); mice, from 3 independent experiments, one-way ANOVA.
(E) Light schedule of light:dark (LD), inverted dark:light (DL), and jet lag (JL) conditions. Shaded areas indicate dark phases, and numbers indicate the harvest times.
(F) Numbers of tumor-infiltrating leukocytes in a B16-F10-OVA tumor, harvested as the indicated time points (1 or 13 h in E) after the onset of the cycle under light:dark (LD, mice), inverted dark:light (DL, mice), or corresponding jet lag (JL, mice) conditions, from 3 independent experiments, unpaired Student’s t test. Shaded areas indicate dark phases. All data are represented as mean SEM, ns, not significant, . See also Figure S1.
Figure 2. Endothelial cells gate circadian leukocyte infiltration and T cell therapy
(A) Normalized total numbers of tumor-infiltrating leukocytes in a B16-F10-OVA model, harvested at ZT1 or ZT13, 24 h after the treatment of isotype control or anti-LFA-1 antibodies; mice, from 5 independent experiments, unpaired Student’s t test.
(B) Gating strategy and quantification of adoptively transferred cells in B16-F10-OVA tumors 2 h post i.v. injection at ZT1 or ZT13; mice, from 2 independent experiments, unpaired Student’s test.
(C) Imaging (left) and quantification (right) of ICAM-1 expression on CD31 endothelial cells in B16-F10-OVA tumors, harvested at ZT1 or ZT13 from control (pooled of B6 mice and Bmal1 mice), NSG or Bmal1 mice; sections, from mice, respectively, from 4 independent experiments, unpaired Student’s t test and cosinor analysis, scale bars: .
(D) Numbers of leukocytes in B16-F10-OVA tumors, harvested at ZT1 or ZT13 from littermate control or Bmal1 mice; mice, from 4 independent experiments, unpaired Student’s test.
(E) Normalized numbers of adoptively transferred leukocytes harvested from B16-F10-OVA tumors 2 h post i.v. injection at ZT1 or ZT13 into littermate control or Bmal1 mice; mice, from 2 independent experiments, unpaired Student’s t test.
(F) Normalized numbers of adoptively transferred OT-I cells in B16-F10-OVA tumors 2 h post injection at ZT1 or ZT13; mice, from 2 independent experiments, unpaired Student’s t test.
(G) Tumor volume after i.v. injection of activated OT-I cells at ZT1 ( mice) or ZT13 ( mice), control ( mice), from 3 independent experiments, Student’s t test.
(H) Tumor volume after i.v. injection of activated OT-I cells at ZT1 (Bmal1 mice, Bmal1 mice) or ZT13 ( mice), control ( mice), from 2 independent experiments, Student’s t test, *Bmal1 ZT1 vs. ZT13.
(I) Tumor volume after i.v. injection of human CD19 CAR-CD28-CD3 cells in DoHH2 tumor-bearing NSG mice at ZT1 (UT, untransduced T cells, , CAR T, mice) or ZT13 (UT, mice, CAR T, mice), control ( mice), from 2 independent experiments, two-way ANOVA.
(J) Number of CAR T cells in DoHH2 tumors 24 h post i.v. injection at ZT1 ( mice) or ZT13 ( mice), from 2 independent experiments, unpaired Student’s t test. All data are represented as mean SEM, ns, not significant, ; ; ; . See also Figure S2.
in highly immune-deficient NOD scid gamma (NSG) mice, indicating that these differences were not caused by the observed oscillations in the tumor leukocyte infiltrate (Figure 2C). We
thus focused on endothelial cells as potential mediators of circadian leukocyte infiltration and grafted tumors into animals exhibiting an inducible, lineage-specific lack of the circadian gene
BMAL1 in endothelial cells (Cdh5 :Bmal1 , Bmal1 ). In Bmal1 mice, time-of-day differences in ICAM-1 expression in tumor endothelial cells were completely abrogated (Figure 2C), yet vascular density was similar to wild-type mice (Figure S1E). Furthermore, total tumor leukocyte numbers in Bmal1 mice (Figure 2D) as well as adoptively transferred leukocyte infiltration to the tumor ceased to be time-of-day dependent (Figure 2E). Together, these data demonstrate that endothelial cells control a rhythmic leukocyte infiltration process into tumors, which peaks in the evening and is dependent on the endothelial cellautonomous circadian clock machinery.

Efficacy of CAR T cell therapy is time-of-day dependent

Given that T cells currently represent the most promising target of anti-cancer immunotherapies, we focused our subsequent functional assays on these cells. To evaluate a potential clinical translational aspect of rhythmic T cell infiltration into tumors, we initially harvested OVA-specific CD8 cells from OT-I mice and performed T cell therapy on mice engrafted with B16-F10-OVA-expressing tumors. After i.v. adoptive transfer of activated OT-I cells into tumor-bearing mice 10 days after tumor engraftment, we detected significantly increased numbers of transferred cells, when T cell infusion had been administered at ZT13 compared with ZT1 (Figure 2F). Of relevance, this yielded significantly better tumor control and reduced tumor burden (Figures 2G and S2F). By contrast, in Bmal1 mice, the time-of-day effects were abrogated, demonstrating that the therapeutic benefit was governed by the endothelial cell-autonomous circadian clock (Figures 2H and S2G). To further evaluate a potential clinical translational aspect of rhythmic T cell infiltration into tumors, we employed human CAR T cells, bearing an antihuman CD19 CAR (CD19 CAR-CD28-CD3). NSG mice were engrafted with DoHH2 cells, a tumor cell line for human diffuse large B cell lymphoma (DLBCL), the major indication for CAR T cell therapy in the clinic. Tumor-bearing mice were treated with CAR T cells in a time-of-day-dependent manner, with infusion of CAR T cells at ZT13 showing a reduction in tumor volume, while infusion at ZT1 had no effect on tumor control (Figures 21 and S 2 H ). This observation was linked to increased CAR T cell homing after 24 h at ZT13 compared with ZT1 in these animals (Figure 2J), in line with our previous observation (Figure 2B). These data demonstrate clearly enhanced T cell therapy efficacy, simply by adjusting the time of injection to the optimal leukocyte infiltration time of the day.

Diurnal differences in TIL phenotype

Knowing that the number of TILs depends on the time of day and that this variable could be exploited for therapeutic purposes, we next assessed whether the phenotype of TILs was also time-ofday dependent. To this end, we performed single-cell RNA sequencing (scRNA-seq) analyses of the CD45 population, sorted from tumors 12 days after engraftment and harvested at four different times of the day, at ZT1, ZT7, ZT13, and ZT19, analogous to the previous experiments. Uniform manifold approximation and projection (UMAP) analyses revealed 18 distinct leukocyte populations, providing insights into subsets and states of leukocytes (Figures S3A-S3C, S4A, and S4B). The most prominent population consisted of macrophages, which made
up of the infiltrate, followed by natural killer (NK) cells ( ) and B cells ( ) (Figures S4A and S4B). Circadian clock genes were expressed rhythmically in all leukocyte subsets with Per1 exhibiting the most wide-spread oscillation, peaking at ~ZT13 (Figures 3A and S4C). Macrophages displayed the highest number of oscillating genes, followed by monocytes and NK cells (Figure S4D). However, each immune cell population exhibited distinct oscillatory gene signatures (Figures S4D-S4F; Table S1). Gene ontology analyses showed that oscillating genes were mostly enriched in metabolism but also showed differences in cell adhesion and differentiation as well as T cell activation (Figure S4E; Table S1). These data indicate that the number and phenotype of immune cells in the TME are highly oscillatory and strongly depend on the time of day.

Time-of-day differences in the tumor cell phenotype

Focusing on T lymphocytes, we performed further analyses and detected 15 subclusters, with cells representing the largest population (Figures 3B, 3C, and S5A-S5C). We observed significant time-of-day-dependent changes in the ratio of cell clusters defined as “anti-” versus “pro-tumorigenic” lymphocytes (see STAR Methods for cluster definition), peaking at ZT13 and troughing at ZT1 (Figure 3D). Furthermore, the ratio in clusters defined as “non-exhausted” versus “exhausted” cells exhibited differences, a phenotype we also observed by flow cytometry (Figures 3E and S5D), which could be further refined to the ratio of cells defined as “effector memory” versus exhausted (Figure 3D). This suggested that cells exhibited a more suppressed phenotype in the morning and a more anti-tumorigenic signature in the evening. To investigate this, we defined the gene signatures of T cells expressed in the evening (ZT13) versus the morning (ZT1) (Table S2). The evening signature was enriched in leukocyte adhesion and T cell activation pathways (Figure 3F) and was indeed strongly correlated with better survival in human melanoma patients from the The Cancer Genome Atlas (TCGA) dataset, indicating higher anti-tumorigenic immunogenicity of this gene expression signature (Figure 3G). These data show that in addition to oscillations in the number of T cells in the TME, their function also changes in a time-of-day-dependent manner.

CD8 cell phenotype oscillates

These expression data indicated that T cell function within the TME exhibited time-of-day differences. Genes associated with T cell activation, suppression, migration, and effector functions showed an oscillation in effector CD8 T cells across the day (Figure 4A). Notably, expression of the immune checkpoint gene Pdcd1 (encoding for PD-1) was rhythmic in this CD8 T cell cluster, peaking at ZT1 and exhibiting a trough at ZT13 (Figures 4A and 4B). This was in line with expression of Pdcd1 in the whole CD8 T cell population, as assessed by qPCR analyses of sorted CD8 T cells from tumors as well as by pseudobulk analyses of the scRNA-seq data (Figures 4C and S5E). PD-1 was also differently expressed at the cell surface protein level, in line with the transcriptomic analyses, yet phase-shifted, similar to what has been described for other circadian-expressed mRNA-protein pairs (Figure 4D). This difference in PD-1 expression could be abrogated when the clock gene Bmal1 was deleted in
Figure 3. Diurnal differences in T cell parameters
(A) Circadian gene Per1 expression in total tumor infiltrating leukocyte clusters with fitted cosinor curve (line) and confidence interval (shaded area).
(B) UMAP representation of scRNA-seq analyses of T cells harvested at 4 different times of the day.
(C) Relative abundance of T cell subsets across different times of the day.
(D) Abundance ratio of specific T cell subsets across different times of the day by scRNA-seq. “Anti-tumorigenic lymphocytes” include NK, T, natural killer T (NKT), CD8 effector, CD4 naive, CD8 effector memory, CD8 mitotic, and CD4 effector clusters, whereas “pro-tumorigenic lymphocytes” include Treg, Treg mitotic, and CD8 exhausted clusters.
(E) Normalized ratio of non-exhausted over exhausted cells across different times of the day by flow cytometry; mice, from 5 independent experiments, cosinor analysis.
(F) Gene ontology analysis of genes oscillating in cells.
(G) Survival analysis in melanoma patients from the TCGA dataset using a high or low expression of the murine ZT13 gene signature in Table S2, logrank test. All data are represented as mean SEM, * . See also Figures S3-S5 and Table S1.
T cells (Cd4cre:Bmal1 , Bmal1 ) (Figure 4E), indicating that rhythmic expression of PD-1 in CD8 cells was circadian and cell-intrinsic. To confirm this, we activated cells in vitro and synchronized the cells (Figures 4F and S6A). Activated CD8 T cells showed a strong oscillatory expression of Pdcd1 after synchronization, which was significantly dampened in T cells harvested from Bmal1 mice and Per1 animals (Figure 4G). PDCD1 also exhibited circadian oscillations in vitro in synchronized human cells, indicating these oscillations to be of relevance in the human setting as well (Figure 4H). These data demonstrated Pdcd1 expression in CD8 T cells to be bona fide circadian and cell autonomous in nature, at both the mRNA and surface protein levels, which might be targetable by anti-PD1 immune check inhibitors in a time-of-day-dependent manner.

Anti-PD-1 therapy is time-of-day dependent

We therefore assessed the relevance of time-of-day-dependent expression of PD-1 by performing timed administration of anti-PD-1 antibodies. Time of day of antibody administration had a strong impact on tumor growth in the B16-F10-OVA melanoma
tumor model, with tumors growing less when therapy was performed in the evening compared with the morning (Figures 5A, 5B, and S6B). In fact, therapy in the morning had surprisingly little effect on tumor burden, even when morning administration preceded that of evening administration by 12 h (Figures 5A, 5B, and S6B). We could confirm these data in an MC-38 colon carcinoma model, indicating that this phenotype extended to other tumor models (Figures S6C and S6D). The time-of-day effect was critically dependent on cells, as anti-PD-1 therapy yielded no benefit in a scenario where cells had previously been depleted (Figures 5C and S6E). Furthermore, the T cell circadian clock machinery was required since circadian oscillations were altered in Bmal1 mice, demonstrating the T cell clock to be mediating the phenotype (Figures 5D, S6F, and S6G). Using several cycles of timed anti-PD-1 therapy and combinations of evening and morning administrations, we observed that time of day of the first treatment was responsible for the difference, in contrast to the subsequent dose (Figures 5E and S6H). We thus focused on the first administration cycle to further elucidate the underlying mechanism and performed flow
Figure 4. cell phenotype oscillates
(A) Heatmap of scRNA-seq analyses of effector T cells, harvested at 4 different times of the day.
(B) Pdcd1 expression by scRNA-seq of CD8 effector T cells at different times of the day.
(C) Pdcd1 expression by qPCR of fluorescence-activated cell sorting (FACS)-sorted CD8 cells from B16-F10-OVA tumors at ZT1 ( mice) or ZT13 ( mice), from 3 independent experiments, unpaired Student’s test.
(D) PD-1 expression by flow cytometry of tumor cells harvested at different time ( mice), from 5 independent experiments, one-way ANOVA and cosinor analysis.
(E) PD-1 expression by flow cytometry of tumor CD8 cells, harvested at ZT1 or ZT13 from littermate control ( mice) or Bmal1 mice ( mice), unpaired Student’s test.
(F) Relative bioluminescence of activated cells generated from PER2::LUC mice post synchronization, representative data from mice.
(G) Pdcd1 expression of activated CD8 cells post synchronization from control ( mice) or Bmal mice ( mice) or Per1 Per2 mice ( mice) from 4 independent experiments, ***cosinor analysis. #Two-way ANOVA.
(H) PDCD1 expression in human activated CD8 T cells post synchronization, , from 4 independent experiments, cosinor analysis. All data are represented as mean SEM, . See also Figure S6.
cytometry analyses of cells in the tumor, 24 h after anti-PD-1 treatment performed in the morning compared with the evening. The results indicated that the administration of anti-PD-1 at ZT13 but not at ZT1 caused a significant increase of T cell activation and degranulation in the tumor, as shown by higher expression CD137 and CD107a levels, respectively (Figure 5F). Additionally, anti-PD-1 treatment at ZT13 also increased the production of interferon (IFN) and Granzyme B in tumor-infiltrating cells compared with ZT1 (Figure 5G). To maximize the observed chrono-therapeutic anti-tumor effects, we combined timed OT-I T cell therapy with timed anti-PD-1 infusion. Activated OT-I T cells were co-injected with anti-PD-1 antibody into phase-shifted mice at ZT13 or at ZT1. Anti-tumor effects were now enhanced, with ZT13 treatment
further reducing tumor burden compared with ZT1 conditions (Figures 5 H and S6I). This indicates that both cell numbers and phenotype dictate circadian differences in the response to tumor immunotherapy. Together, our data demonstrate a strong time-of-day-dependent phenotype of the tumor immune cell infiltrate, which can be harnessed with timed administration of immunotherapies to maximize anti-tumor effects.

Circadian immune cell phenotype in human cancers

To interrogate whether time-of-day oscillations in the immune cell infiltrate also exist in human tumors, we analyzed human melanoma samples from patients that had undergone surgical removal of the tumor at different times of the day. Using immunohistochemical analyses, we found that numbers of and
Figure 5. Circadian efficacy of anti-PD-1 treatment
(A) Tumor volume after anti-PD-1 treatment (arrows) administered at ZT1 or ZT13, anti-PD-1 ( mice), control ( mice), from 3 independent experiments, two-way ANOVA.
(B) Tumor volume 14 days post tumor engraftment, anti-PD-1 ( mice), control ( mice), from 3 independent experiments, cosinor analysis.
(C) Tumor volume after anti-PD-1 treatment (arrows) administered at ZT1 or ZT13 with or without ( mice) CD8 depletion.
(D) Tumor volume after anti-PD-1 treatment (arrows) administered at ZT1 or ZT13 in control (ZT1, mice, ZT13, mice) or in Bmal1 (ZT1, mice, ZT13, mice), from 3 independent experiments, two-way ANOVA.
(E) Tumor volume after two doses of anti-PD-1 treatment (arrows) administered at ZT1 or ZT13; mice, from 2 independent experiments, two-way ANOVA. (F) Phenotype of tumor cell harvested at ZT1 or ZT13, 24 h post isotype control ( mice) or anti-PD-1 treatment (ZT1, mice, ZT13, mice), from 4 independent experiments, two-way ANOVA, and Sidak’s post test.
(G) Cytokine production by tumor CD8 T cells harvested at ZT1 ( mice) or ZT13 ( mice) 24 h post anti-PD-1 treatment, from 4 independent experiments, unpaired Student’s t test.
(H) Tumor volume after control ( mice) or combination therapy of anti-PD-1 and OT-I cell administration (arrow), administered at ZT1 ( mice) or ZT13 ( 5 mice), from 2 independent experiments, two-way ANOVA. Dotted lines indicate tumor growth from Figure 2F. All data are represented as mean SEM, ns, not significant, . See also Figure S6.
CD8 cells exhibited a striking time-of-day-dependent profile with a peak in the early afternoon (Figures 6A and 6B). This indicated that the time of day also dictated numbers of TILs in humans. To obtain broader information with respect to the timed immune cell infiltrate in human tumors, we utilized our mouse scRNA-seq datasets together with publicly available, timestamped mouse skin scRNA-seq and bulk RNA-seq datasets to train a machine learning model to be able to time stamp sequencing datasets where information on the time of day was not available (Figure 6C). We trained a ZeitZeiger algorithm, using either all genes or conserved coefficients ( 57 genes), which were predicted by ZeitZeiger to oscillate (Figure 6C; Table S3). Both approaches predicted the time of testing data wellincluding healthy skin as well as melanoma samples-with an average error of (Figures S6J-S6L; Table S3). We further validated the ZeitZeiger model on public mouse skin RNA-seq and microarray datasets and again achieved good time stamping (Figures 6D and S6M). This allowed us to correctly identify the time of day of healthy tissues in mice as well as our mouse tumor samples in a blinded manner with an error of . Using the trained ZeitZeiger model, we next aimed to predict the time of day of harvest also in human biopsies. Indeed, using publicly
available microarray data from healthy human skin where time of day information was known, we could adequately predict the time of day in these samples with an error of (Figures 6E and S6N). Knowing that the algorithm adequately assigned time of day to human samples, we next assigned time-of-day information to 103 primary melanoma samples from the TCGA dataset. Interestingly, we observed time-of-day differences in the ratio of exhausted cells versus non-exhausted T cells, analogous to the mouse data, but-as expected-phase-shifted due to inverted immune oscillations between diurnal humans and nocturnal mice, with a higher anti-tumor ratio in the morning compared with the afternoon (Figure 6F). These data indicate that both leukocyte numbers and phenotype in human TILs are time-of-day dependent.
Lastly, we investigated whether the evening RNA signature identified in mouse melanoma T cells was sufficient to correlate with the response of patients to ICBs. Indeed, the signature unveiled a significant difference in the response rate of patients, with patients that exhibited a higher expression of this signature in T cells also showing favorable response to ICB treatment (Figure 6G; Table S2). Together, our data provide compelling evidence that the immune cell compartment of
Figure 6. Circadian immune cell phenotype in human cancers (A and B) Imaging (left) and quantification (right) of CD4 patients) and CD8 patients) T cells in human melanoma, harvested at different times of the day, cosinor analysis, scale bar: , CD8 or CD4 staining in red.
(C) Schematic for the algorithm training and validation to predict time of day in non-time stamped samples. A final algorithm containing 57 genes was used for external validation and prediction.
(D and E) Error (hours) of the time prediction algorithm in external validation in mouse ( mice) (D) and human ( samples) (E) tissues.
(F) Exhausted/non-exhausted T cell ratio in primary melanoma tissues from the TCGA melanoma dataset, plotted by using the algorithm to predict time of biopsies, morning, patients; afternoon, patients, Student’s test.
(G) Enrichment of the murine ZT13-gene expression signature in patients responding or not responding to ICB treatment (non-responders, cells; responders, cells). All data are represented as median and quartile. See also Figure S6 and Tables S2 and S3.
melanomas is surprisingly time-of-day dependent and that this can be harnessed therapeutically with timed interventions of CAR T cell infusions or immune checkpoint blockade.

DISCUSSION

Here, we provide evidence for circadian dynamics in the tumor immune microenvironment, with functional implications for the optimal time point for the administration of immunotherapies. In patients with melanoma, non-small cell lung cancer (NSCLC) patients, and patients with metastatic squamous cell carcinoma of the esophagus, retrospective data of matched cohorts treated with ICB was recently indicated to be time-of-day sensitive, with better survival observed after ICB infusion in the morning compared with the afternoon. These clinical observations are in line with what we define here in mechanistic detail. Future randomized and prospective clinical trials will be required to support that claim.
A key question in the field is the mechanism of these time-ofday effects. Antibody titers remain high over weeks in patients, which may indicate that the observed long-term circadian effects are mediated via the timing of the first interventions. Our data indicate that the first cycle of therapy yields different circadian outcomes and that the timing of the subsequent dose(s) may not be as relevant. Data from retrospective clinical trials support the potential importance of the timing of the first infusion(s) for therapeutic benefits. Our data suggest that an antibody infused at different times of the day encounters differences in both the number of immune cells present in the tumor as well as their levels of PD-1. This would mean that the product of circadian changes in numbers and the phenotype of TILsand thus the target of the antibody-exhibits circadian oscillations, which could be responsible for this chronotherapeutic impact. However, whether local levels of antibody oscillate, such as found within the tumor interstitium, is not known. Tumor macrophages have been shown to quickly strip antibodies off adjacent T cells, which may be an additional mechanism for
time-of-day sensitivity, given that we also see oscillations in macrophage phenotype in the TME. We have focused on CD8 T cells in this study, as they currently represent major therapeutic targets in the clinic. Our data show that antibody administration alters the phenotype of CD8 T cells within the tumor in a time-of-day-dependent manner, already 24 h after the infusion, indicating that these cells are involved in the initial circadian responses. Nevertheless, oscillations in other immune cell lineages that we describe, notably of myeloid cells, may be of additional functional relevance for the circadian effect. The contribution of different immune cell populations in this circadian anti-tumor response remains to be investigated in future studies.
Ultimately, it appears that challenging the immune system with an antibody at a specific time of day not only changes the quantity but also the quality of the response so that the immune system, once stimulated at the “wrong” time, may not be able to respond anymore to the same level and quality as an immune system challenged at the “right” time-just 12 h apart. In an analogous manner to the anti-tumor response described here, a depot effect exists in vaccines, with antigen being present over several days. Yet, time-of-day effects are observed also in vaccination studies, in the pre-clinical setting, and in patients, with higher immune sensitivity present at very similar time windows. The precise mechanisms at play underlying these phenomena remain to be identified.
Sex differences have been reported in cancer immunotherapy, with retrospective analyses indicating that men often respond better compared with women. Of note, with respect to chronoimmunotherapy, subgroup analyses did reveal that women had a particularly high benefit in overall survival with immune checkpoint inhibition infused earlier in the day compared with in the afternoon. To what extent and what the mechanisms are will need to be investigated in more detail in future studies. The circadian dynamics in the TME that we uncover here should prove to be of use in the clinic, as tumor biopsies represent currently the main means of assessing the level of immune cell infiltration and defining the ensuing therapeutic regimen. Given that this immune cell infiltration is strongly time-of-day dependent, therapy may differ, depending on when biopsies were performed. Together, we demonstrate in mechanistic detail that TILs exhibit circadian oscillations in both number and phenotype, which can be exploited therapeutically with timed administration of CAR T cells or ICBs.

Limitations of the study

We focused on melanoma in this study, and it remains to be demonstrated in future studies whether this also applies to other tumors. Our observations in mice indicate that at least time-ofday differences in ICB also extend to colon carcinoma. Furthermore, retrospective clinical data indicate that optimization of infusion timing may provide benefits for several types of cancer, including melanoma, NSLSC, and metastatic esophageal squamous cell carcinoma patients. Although our data indicate the optimal time point of infusion to be the evening in mice, the best time point for humans may be the morning hours, given the different periods of behavioral rest and activity in the two species. However, the precise time window remains to be
defined, as current data only assessed morning versus afternoon effects.
Although our data indicate that the lineage-specific, inducible deletion of Bmal1 in endothelial cells affects the oscillatory behavior of these cells and, by extension, the recruitment of leukocytes, there may be additional factors that could contribute to the observed abrogation of time-of-day-dependent leukocyte migration.

STAR * METHODS

Detailed methods are provided in the online version of this paper and include the following:
  • KEY RESOURCES TABLE
  • RESOURCE AVAILABILITY
  • Lead contact
  • Materials availability
  • Data and code availability
  • EXPERIMENTAL MODEL AND STUDY PARTICIPANT DETAILS – Animals
  • METHOD DETAILS
  • Tumor cell lines and inoculation
  • Mouse T cell activation
  • Mouse T cell synchronization
  • Human T cell activation and synchronization
  • Luciferase-expressing cells
  • Tissue digestion and single-cell preparation
  • Flow Cytometry
  • Mouse leukocyte adoptive transfer
  • EdU assay
  • CD8 T cell sorting and RNA extraction
  • RNA extraction, reverse transcription and qPCR
  • In vivo antibody treatments
  • Human CAR T cell generation
  • Human CAR T cell adoptive transfer
  • Immunofluorescence imaging
  • Human immunohistochemistry
  • Single cell sequencing and analysis
  • ZeitZeiger analysis
  • Estimation of cell composition in TCGA melanoma patients
  • QUANTIFICATION AND STATISTICAL ANALYSIS

SUPPLEMENTAL INFORMATION

Supplemental information can be found online at https://doi.org/10.1016/j.cell. 2024.04.015.

ACKNOWLEDGMENTS

We thank Urs Albrecht and Jürgen Ripperger for providing Per1-/-Per2-/mice. We thank the IGE3 Genomics Platform and the flow cytometry core facility of the University of Geneva. We thank Dr. Lydia Lutes, Dr. Laure Garnier, Dr. Flore Sinturel, and Dr. Liqing Cheng for invaluable input and discussions. This study received support from the European Research Council (ERC CoG 101001233, CIRCADYN [C.S.]), the Swiss National Science Foundation (SNSF, 310030_182417/1 [C.S.] and 310030_185255 [S.H.]), the Swiss Cancer League (KLS-4836-08-2019 [C.S.]), the Geneva Cancer League (2106 [C.S.]), the EU ITN (813284, INTEGRATA [C.S.]), and the Fondation privé des HUG (RC05-03 [C.S.]). C.J. is supported by a SNSF PRIMA fellowship (PR00P3_179727), the Swiss Cancer League (KFS 5250-02-2021), and the Geneva Cancer League (GCL, 2007). M.C. is the recipient of an iGE3 PhD salary award. P.-C. Ho received support from the Cancer Research Institute, the Melanoma Research Alliance, the Ludwig Cancer Research, and the Swiss National Science Foundation (310030L_208130). M.J.P. received support
from the ISREC Foundation, Ludwig Cancer Research, and NIH grants P01CA240239 and R01-CA218579. R.B. was funded by a Postdoc Mobility fellowship and a Return Grant of the SNSF (P400PM_183852 and P5R5PM_203164). Research in the A.D.G. lab is supported by Research Foundation Flanders (FWO) (Fundamental Research Grant, G0B4620N; FWO SBO grant for “ANTIBODY” consortium), KU Leuven (C1 grant, C14/19/098, and C3 grants, C3/21/ 037 and C3/22/022), Kom op Tegen Kanker (KOTK/2018/11509/1 and KOTK/ 2019/11955/1), VLIR-UOS (iBOF grant, iBOF/21/048, for “MIMICRY” consortium), and Olivia Hendrickx Research Fund (OHRF Immuunbiomarkers). C.D. was funded by Swiss National Science Foundation grant 310030_184708/1, the Vontobel Foundation, the Novartis Consumer Health Foundation, EFSD/Novo Nordisk Programme for Diabetes Research in Europe, Swiss Life Foundation, the Olga Mayenfisch Foundation, the Fondation pour l’innovation sur le cancer et la biologie, the Ligue Pulmonaire Genevoise, Swiss Cancer League KFS-5266-02-2021-R, the Velux Foundation, the Leenaards Foundation, the ISREC Foundation, and the Gertrude von Meissner Foundation (V.P. and C.D.). Y.W. received support from the Science Foundation of Peking University Cancer Hospital (2022-21). Extended data Figure 6A and the graphical abstract were created with BioRender.

AUTHOR CONTRIBUTIONS

Conceptualization: C.W., Z.L., S.H., C.J., C.D., A.D.G., P.-C. Ho, O.M., F.S., M.J.P., and C.S.; methodology: C.W., Q.Z., P.C., S.N., Z.S., and A.D.G.; investigation: C.W., Q.Z., S.W., R.B., Y.W., C.B., R.P., M.C., M.S., P.-C. Hsueh, S.H.M., S.J., V.P., F.S., Z.L., and Z.M.G.; visualization: C.W., Q.Z., P.C., S.N., A.D.G., and C.S.; supervision: C.S.; writing: C.W., Q.Z., M.J.P., and C.S.

DECLARATION OF INTERESTS

C.S. has been a consultant for Bayer and received speaker fees from Abbvie. M.J.P. has been a consultant for AstraZeneca, Debiopharm, Elstar Therapeutics, ImmuneOncia, KSQ Therapeutics, MaxiVax, Merck, Molecular Partners, Third Rock Ventures, and Tidal. F.S. received consulting fees from BMS/ Celgene, Incyte, Kite/Gilead, speaker fees from Kite/Gilead, Incyte, travel support from Kite/Gilead, Novartis, AstraZeneca, Neovii, Janssen, and research funding from Kite/Gilead, Novartis, and BMS/Celgene. A.D.G. received consulting/advisory/lecture honoraria or research funding from Boehringer Ingelheim, SOTIO, Miltenyi Biotec, Novigenix, and IsoPlexis. R.B. received speaker fees from Janssen and is a mentee of the ENDEAVOUR-Breast program of Daichii Sankyo. All these relationships are unrelated to this study.
Received: October 12, 2023
Revised: February 2, 2024
Accepted: April 16, 2024
Published: May 8, 2024

REFERENCES

  1. Mellman, I., Coukos, G., and Dranoff, G. (2011). Cancer immunotherapy comes of age. Nature 480, 480-489. https://doi.org/10.1038/ nature10673.
  2. Tawbi, H.A., Schadendorf, D., Lipson, E.J., Ascierto, P.A., Matamala, L., Castillo Gutiérrez, E., Rutkowski, P., Gogas, H.J., Lao, C.D., De Menezes, J.J., et al. (2022). Relatlimab and Nivolumab versus Nivolumab in Untreated Advanced Melanoma. N. Engl. J. Med. 386, 24-34. https://doi. org/10.1056/NEJMoa2109970.
  3. Cappell, K.M., and Kochenderfer, J.N. (2023). Long-term outcomes following CAR T cell therapy: what we know so far. Nat. Rev. Clin. Oncol. 20, 359-371. https://doi.org/10.1038/s41571-023-00754-1.
  4. Flugel, C.L., Majzner, R.G., Krenciute, G., Dotti, G., Riddell, S.R., Wagner, D.L., and Abou-El-Enein, M. (2023). Overcoming on-target, off-tumour toxicity of CAR T cell therapy for solid tumours. Nat. Rev. Clin. Oncol. 20, 49-62. https://doi.org/10.1038/s41571-022-00704-3.
  5. Tumeh, P.C., Harview, C.L., Yearley, J.H., Shintaku, I.P., Taylor, E.J.M., Robert, L., Chmielowski, B., Spasic, M., Henry, G., Ciobanu, V., et al.
    (2014). PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature 515, 568-571. https://doi.org/10.1038/nature13954.
  6. Wang, C., Barnoud, C., Cenerenti, M., Sun, M., Caffa, I., Kizil, B., Bill, R., Liu, Y., Pick, R., Garnier, L., et al. (2023). Dendritic cells direct circadian anti-tumour immune responses. Nature 614, 136-143. https://doi.org/ 10.1038/s41586-022-05605-0.
  7. Scheiermann, C., Gibbs, J., Ince, L., and Loudon, A. (2018). Clocking in to immunity. Nat. Rev. Immunol. 18, 423-437. https://doi.org/10.1038/ s41577-018-0008-4.
  8. Arjona, A., Silver, A.C., Walker, W.E., and Fikrig, E. (2012). Immunity’s fourth dimension: approaching the circadian-immune connection. Trends Immunol. 33, 607-612. https://doi.org/10.1016/j.it.2012.08.007.
  9. Curtis, A.M., Bellet, M.M., Sassone-Corsi, P., and O’Neill, L.A.J. (2014). Circadian clock proteins and immunity. Immunity 40, 178-186. https:// doi.org/10.1016/j.immuni.2014.02.002.
  10. Cermakian, N., Stegeman, S.K., Tekade, K., and Labrecque, N. (2022). Circadian rhythms in adaptive immunity and vaccination. Semin. Immunopathol. 44, 193-207. https://doi.org/10.1007/s00281-021-00903-7.
  11. Man, K., Loudon, A., and Chawla, A. (2016). Immunity around the clock. Science 354, 999-1003. https://doi.org/10.1126/science.aah4966.
  12. Palomino-Segura, M., and Hidalgo, A. (2021). Circadian immune circuits. J. Exp. Med. 218, e20200798. https://doi.org/10.1084/jem. 20200798.
  13. Wang, C., Lutes, L.K., Barnoud, C., and Scheiermann, C. (2022). The circadian immune system. Sci. Immunol. 7, eabm2465. https://doi.org/10. 1126/sciimmunol.abm2465.
  14. He, W., Holtkamp, S., Hergenhan, S.M., Kraus, K., de Juan, A., Weber, J., Bradfield, P., Grenier, J.M.P., Pelletier, J., Druzd, D., et al. (2018). Circadian Expression of Migratory Factors Establishes Lineage-Specific Signatures that Guide the Homing of Leukocyte Subsets to Tissues. Immunity 49, 1175-1190.e7. https://doi.org/10.1016/j.immuni.2018.10.007.
  15. Druzd, D., Matveeva, O., Ince, L., Harrison, U., He, W., Schmal, C., Herzel, H., Tsang, A.H., Kawakami, N., Leliavski, A., et al. (2017). Lymphocyte Circadian Clocks Control Lymph Node Trafficking and Adaptive Immune Responses. Immunity 46, 120-132. https://doi.org/10.1016/j.immuni. 2016.12.011.
  16. Shimba, A., Cui, G., Tani-ichi, S., Ogawa, M., Abe, S., Okazaki, F., Kitano, S., Miyachi, H., Yamada, H., Hara, T., et al. (2018). Glucocorticoids Drive Diurnal Oscillations in T Cell Distribution and Responses by Inducing Interleukin-7 Receptor and CXCR4. Immunity 48, 286-298.e6. https:// doi.org/10.1016/j.immuni.2018.01.004.
  17. Casanova-Acebes, M., Nicolás-Ávila, J.A., Li, J.L., García-Silva, S., Balachander, A., Rubio-Ponce, A., Weiss, L.A., Adrover, J.M., Burrows, K., AGonzález, N., et al. (2018). Neutrophils instruct homeostatic and pathological states in naive tissues. J. Exp. Med. 215, 2778-2795. https://doi.org/ 10.1084/jem. 20181468.
  18. Nguyen, K.D., Fentress, S.J., Qiu, Y., Yun, K., Cox, J.S., and Chawla, A. (2013). Circadian gene Bmal1 regulates diurnal oscillations of Ly6C(hi) inflammatory monocytes. Science 341, 1483-1488. https://doi.org/10. 1126/science. 1240636.
  19. Cervantes-Silva, M.P., Carroll, R.G., Wilk, M.M., Moreira, D., Payet, C.A., O’Siorain, J.R., Cox, S.L., Fagan, L.E., Klavina, P.A., He, Y., et al. (2022). The circadian clock influences cell responses to vaccination by regulating dendritic cell antigen processing. Nat. Commun. 13, 7217. https:// doi.org/10.1038/s41467-022-34897-z.
  20. Ince, L.M., Barnoud, C., Lutes, L.K., Pick, R., Wang, C., Sinturel, F., Chen, C.S., de Juan, A., Weber, J., Holtkamp, S.J., et al. (2023). Influence of circadian clocks on adaptive immunity and vaccination responses. Nat. Commun. 14, 476. https://doi.org/10.1038/s41467-023-35979-2.
  21. Nobis, C.C., Dubeau Laramée, G., Kervezee, L., Maurice De Sousa, D., Labrecque, N., and Cermakian, N. (2019). The circadian clock of CD8 T cells modulates their early response to vaccination and the rhythmicity of related signaling pathways. Proc. Natl. Acad. Sci. USA 116, 2007720086. https://doi.org/10.1073/pnas.1905080116.
  22. Gibbs, J., Ince, L., Matthews, L., Mei, J., Bell, T., Yang, N., Saer, B., Begley, N., Poolman, T., Pariollaud, M., et al. (2014). An epithelial circadian clock controls pulmonary inflammation and glucocorticoid action. Nat. Med. 20, 919-926. https://doi.org/10.1038/nm.3599.
  23. Scheiermann, C., Kunisaki, Y., Lucas, D., Chow, A., Jang, J.E., Zhang, D., Hashimoto, D., Merad, M., and Frenette, P.S. (2012). Adrenergic nerves govern circadian leukocyte recruitment to tissues. Immunity 37, 290-301. https://doi.org/10.1016/j.immuni.2012.05.021.
  24. Hazan, G., Duek, O.A., Alapi, H., Mok, H., Ganninger, A., Ostendorf, E., Gierasch, C., Chodick, G., Greenberg, D., and Haspel, J.A. (2023). Biological rhythms in COVID-19 vaccine effectiveness in an observational cohort study of 1.5 million patients. J. Clin. Invest. 133, e167339. https://doi.org/ 10.1172/JCl167339.
  25. Long, J.E., Drayson, M.T., Taylor, A.E., Toellner, K.M., Lord, J.M., and Phillips, A.C. (2016). Morning vaccination enhances antibody response over afternoon vaccination: A cluster-randomised trial. Vaccine 34, 26792685. https://doi.org/10.1016/j.vaccine.2016.04.032.
  26. Ho, P.C., Bihuniak, J.D., Macintyre, A.N., Staron, M., Liu, X., Amezquita, R., Tsui, Y.C., Cui, G., Micevic, G., Perales, J.C., et al. (2015). Phosphoenolpyruvate Is a Metabolic Checkpoint of Anti-tumor T Cell Responses. Cell 162, 1217-1228. https://doi.org/10.1016/j.cell.2015.08.012.
  27. Raskov, H., Orhan, A., Christensen, J.P., and Gögenur, I. (2021). Cytotoxic CD8+ T cells in cancer and cancer immunotherapy. Br. J. Cancer 124, 359-367. https://doi.org/10.1038/s41416-020-01048-4.
  28. National Comprehensive Cancer Network (2024). B-Cell Lymphomas. https://www.nccn.org/guidelines/guidelines-detail?category=1&id=1480.
  29. Neelapu, S.S., Dickinson, M., Munoz, J., Ulrickson, M.L., Thieblemont, C., Oluwole, O.O., Herrera, A.F., Ujjani, C.S., Lin, Y., Riedell, P.A., et al. (2022). Axicabtagene ciloleucel as first-line therapy in high-risk large B-cell lymphoma: the phase 2 ZUMA-12 trial. Nat. Med. 28, 735-742. https://doi. org/10.1038/s41591-022-01731-4.
  30. Arjona, A., and Sarkar, D.K. (2005). Circadian oscillations of clock genes, cytolytic factors, and cytokines in rat NK cells. J. Immunol. 174, 76187624. https://doi.org/10.4049/jimmunol.174.12.7618.
  31. Hughey, J.J., Hastie, T., and Butte, A.J. (2016). ZeitZeiger: supervised learning for high-dimensional data from an oscillatory system. Nucleic Acids Res. 44, e80. https://doi.org/10.1093/nar/gkw030.
  32. Zhao, Y., Liu, M., Chan, X.Y., Tan, S.Y., Subramaniam, S., Fan, Y., Loh, E., Chang, K.T.E., Tan, T.C., and Chen, Q. (2017). Uncovering the mystery of opposite circadian rhythms between mouse and human leukocytes in humanized mice. Blood 130, 1995-2005. https://doi.org/10.1182/blood-2017-04-778779.
  33. Sade-Feldman, M., Yizhak, K., Bjorgaard, S.L., Ray, J.P., de Boer, C.G., Jenkins, R.W., Lieb, D.J., Chen, J.H., Frederick, D.T., Barzily-Rokni, M., et al. (2018). Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma. Cell 175, 998-1013.e20. https://doi. org/10.1016/j.cell.2018.10.038.
  34. Naulaerts, S., Datsi, A., Borras, D.M., Antoranz Martinez, A., Messiaen, J., Vanmeerbeek, I., Sprooten, J., Laureano, R.S., Govaerts, J., Panovska, D., et al. (2023). Multiomics and spatial mapping characterizes human CD8 T cell states in cancer. Sci. Transl. Med. 15, eadd1016. https://doi.org/ 10.1126/scitransImed.add1016.
  35. Qian, D.C., Kleber, T., Brammer, B., Xu, K.M., Switchenko, J.M., Jano-paul-Naylor, J.R., Zhong, J., Yushak, M.L., Harvey, R.D., Paulos, C.M., et al. (2021). Effect of immunotherapy time-of-day infusion on overall survival among patients with advanced melanoma in the USA (MEMOIR): a propensity score-matched analysis of a single-centre, longitudinal study. Lancet Oncol. 22, 1777-1786. https://doi.org/10.1016/S1470-2045(21) 00546-5.
  36. Yeung, C., Kartolo, A., Tong, J., Hopman, W., and Baetz, T. (2023). Association of circadian timing of initial infusions of immune checkpoint inhibitors with survival in advanced melanoma. Immunotherapy 15, 819-826. https://doi.org/10.2217/imt-2022-0139.
  37. Gonçalves, L., Gonçalves, D., Esteban-Casanelles, T., Barroso, T., Soares de Pinho, I., Lopes-Brás, R., Esperança-Martins, M., Patel, V., Torres, S., Teixeira de Sousa, R., et al. (2023). Immunotherapy around the Clock: Impact of Infusion Timing on Stage IV Melanoma Outcomes. Cells 12, 2068. https://doi.org/10.3390/cells12162068.
  38. Karaboué, A., Collon, T., Pavese, I., Bodiguel, V., Cucherousset, J., Zakine, E., Innominato, P.F., Bouchahda, M., Adam, R., and Lévi, F. (2022). Time-Dependent Efficacy of Checkpoint Inhibitor Nivolumab: Results from a Pilot Study in Patients with Metastatic Non-Small-Cell Lung Cancer. Cancers (Basel) 14, 896. https://doi.org/10.3390/cancers14040896.
  39. Rousseau, A., Tagliamento, M., Auclin, E., Aldea, M., Frelaut, M., Levy, A., Benitez, J.C., Naltet, C., Lavaud, P., Botticella, A., et al. (2023). Clinical outcomes by infusion timing of immune checkpoint inhibitors in patients with advanced non-small cell lung cancer. Eur. J. Cancer 182, 107-114. https://doi.org/10.1016/j.ejca.2023.01.007.
  40. Nomura, M., Hosokai, T., Tamaoki, M., Yokoyama, A., Matsumoto, S., and Muto, M. (2023). Timing of the infusion of nivolumab for patients with recurrent or metastatic squamous cell carcinoma of the esophagus influences its efficacy. Esophagus 20, 722-731. https://doi.org/10.1007/s10388-023-01006-y.
  41. Landré, T., Karaboué, A., Buchwald, Z.S., Innominato, P.F., Qian, D.C., Assié, J.B., Chouaïd, C., Lévi, F., and Duchemann, B. (2024). Effect of immunotherapy-infusion time of day on survival of patients with advanced cancers: a study-level meta-analysis. ESMO Open 9, 102220. https://doi. org/10.1016/j.esmoop.2023.102220.
  42. Martin, A.M., Cagney, D.N., Catalano, P.J., Alexander, B.M., Redig, A.J., Schoenfeld, J.D., and Aizer, A.A. (2018). Immunotherapy and Symptomatic Radiation Necrosis in Patients With Brain Metastases Treated With Stereotactic Radiation. JAMA Oncol. 4, 1123-1124. https://doi.org/10. 1001/jamaoncol.2017.3993.
  43. Postow, M.A., Goldman, D.A., Shoushtari, A.N., Betof Warner, A., Callahan, M.K., Momtaz, P., Smithy, J.W., Naito, E., Cugliari, M.K., Raber, V., et al. (2022). Adaptive Dosing of Nivolumab + Ipilimumab Immunotherapy Based Upon Early, Interim Radiographic Assessment in Advanced Melanoma (The ADAPT-IT Study). J. Clin. Oncol. 40, 1059-1067. https://doi. org/10.1200/JCO.21.01570.
  44. Arlauckas, S.P., Garris, C.S., Kohler, R.H., Kitaoka, M., Cuccarese, M.F., Yang, K.S., Miller, M.A., Carlson, J.C., Freeman, G.J., Anthony, R.M., et al. (2017). In vivo imaging reveals a tumor-associated macrophagemediated resistance pathway in anti-PD-1 therapy. Sci. Transl. Med. 9, eaal3604. https://doi.org/10.1126/scitransImed.aal3604.
  45. Conforti, F., Pala, L., Bagnardi, V., De Pas, T., Martinetti, M., Viale, G., Gelber, R.D., and Goldhirsch, A. (2018). Cancer immunotherapy efficacy and patients’ sex: a systematic review and meta-analysis. Lancet Oncol. 19, 737-746. https://doi.org/10.1016/S1470-2045(18)30261-4.
  46. Lévi, F.A., Okyar, A., Hadadi, E., Innominato, P.F., and Ballesta, A. (2024). Circadian Regulation of Drug Responses: Toward Sex-Specific and Personalized Chronotherapy. Annu. Rev. Pharmacol. Toxicol. 64, 89-114. https://doi.org/10.1146/annurev-pharmtox-051920-095416.
  47. Duan, J., Ngo, M.N., Karri, S.S., Tsoi, L.C., Gudjonsson, J.E., Shahbaba, B., Lowengrub, J., and Andersen, B. (2023). tauFisher accurately predicts circadian time from a single sample of bulk and single-cell transcriptomic data. Preprint at bioRxiv. https://doi.org/10.1101/2023.04.04.535473.
  48. Welz, P.-S., Zinna, V.M., Symeonidi, A., Koronowski, K.B., Kinouchi, K., Smith, J.G., Guillén, I.M., Castellanos, A., Furrow, S., Aragón, F., et al. (2019). BMAL1-Driven Tissue Clocks Respond Independently to Light to Maintain Homeostasis. Cell 177, 1436-1447.e12. https://doi.org/10. 1016/j.cell.2019.05.009.
  49. Wang, H., van Spyk, E., Liu, Q., Geyfman, M., Salmans, M.L., Kumar, V., Ihler, A., Li, N., Takahashi, J.S., and Andersen, B. (2017). Time-Restricted Feeding Shifts the Skin Circadian Clock and Alters UVB-Induced DNA Damage. Cell Rep. 20, 1061-1072. https://doi.org/10.1016/j.celrep. 2017.07.022.
  50. Tsujihana, K., Tanegashima, K., Santo, Y., Yamada, H., Akazawa, S., Nakao, R., Tominaga, K., Saito, R., Nishito, Y., Hata, R.-I., et al. (2022). Circadian protection against bacterial skin infection by epidermal CXCL14mediated innate immunity. Proc. Natl. Acad. Sci. USA 119, e2116027119. https://doi.org/10.1073/pnas.2116027119.
  51. Geyfman, M., Kumar, V., Liu, Q., Ruiz, R., Gordon, W., Espitia, F., Cam, E., Millar, S.E., Smyth, P., Ihler, A., et al. (2012). Brain and muscle Arnt-like protein-1 (BMAL1) controls circadian cell proliferation and susceptibility to UVB-induced DNA damage in the epidermis. Proc. Natl. Acad. Sci. USA 109, 11758-11763. https://doi.org/10.1073/pnas.1209592109.
  52. Spörl, F., Korge, S., Jürchott, K., Wunderskirchner, M., Schellenberg, K., Heins, S., Specht, A., Stoll, C., Klemz, R., Maier, B., et al. (2012). Krüp-pel-like factor 9 is a circadian transcription factor in human epidermis that controls proliferation of keratinocytes. Proc. Natl. Acad. Sci. USA 109, 10903-10908. https://doi.org/10.1073/pnas.1118641109.
  53. del Olmo, M., Spörl, F., Korge, S., Jürchott, K., Felten, M., Grudziecki, A., de Zeeuw, J., Nowozin, C., Reuter, H., Blatt, T., et al. (2022). Inter-layer and inter-subject variability of diurnal gene expression in human skin. NAR Genom. Bioinform. 4, Iqac097. https://doi.org/10.1093/nargab/Iqac097.
  54. Wu, G., Ruben, M.D., Schmidt, R.E., Francey, L.J., Smith, D.F., Anafi, R.C., Hughey, J.J., Tasseff, R., Sherrill, J.D., Oblong, J.E., et al. (2018). Popula-tion-level rhythms in human skin with implications for circadian medicine. Proc. Natl. Acad. Sci. USA 115, 12313-12318. https://doi.org/10.1073/ pnas. 1809442115.
  55. Sherrill, J.D., Finlay, D., Binder, R.L., Robinson, M.K., Wei, X., Tiesman, J.P., Flagler, M.J., Zhao, W., Miller, C., Loftus, J.M., et al. (2021). Transcriptomic analysis of human skin wound healing and rejuvenation following ablative fractional laser treatment. PLoS One 16, e0260095. https://doi.org/10.1371/journal.pone. 0260095.
  56. Hao, Y., Hao, S., Andersen-Nissen, E., Mauck, W.M., 3rd, Zheng, S., Butler, A., Lee, M.J., Wilk, A.J., Darby, C., Zager, M., et al. (2021). Integrated analysis of multimodal single-cell data. Cell 184, 3573-3587.e29. https:// doi.org/10.1016/j.cell.2021.04.048.
  57. Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., Feng, T., Zhou, L., Tang, W., Zhan, L., et al. (2021). clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb) 2, 100141. https://doi. org/10.1016/j.xinn.2021.100141.
  58. Carlucci, M., Kriščiūnas, A., Li, H., Gibas, P., Koncevičius, K., Petronis, A., and Oh, G. (2019). DiscoRhythm: an easy-to-use web application and package for discovering rhythmicity. Bioinformatics 36, 1952-1954. https://doi.org/10.1093/bioinformatics/btz834.
  59. Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14, 7. https://doi.org/10.1186/1471-2105-14-7.
  60. Colaprico, A., Silva, T.C., Olsen, C., Garofano, L., Cava, C., Garolini, D., Sabedot, T.S., Malta, T.M., Pagnotta, S.M., Castiglioni, I., et al. (2016). TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 44, e71. https://doi.org/10.1093/nar/ gkv1507.
  61. Leek, J.T., Johnson, W.E., Parker, H.S., Jaffe, A.E., and Storey, J.D. (2012). The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882-883. https://doi.org/10.1093/bioinformatics/bts034.
  62. Holtkamp, S.J., Ince, L.M., Barnoud, C., Schmitt, M.T., Sinturel, F., Pilorz, V., Pick, R., Jemelin, S., Mühlstädt, M., Boehncke, W.H., et al. (2021). Circadian clocks guide dendritic cells into skin lymphatics. Nat. Immunol. 22, 1375-1381. https://doi.org/10.1038/s41590-021-01040-x.
  63. Wittenbrink, N., Ananthasubramaniam, B., Münch, M., Koller, B., Maier, B., Weschke, C., Bes, F., de Zeeuw, J., Nowozin, C., Wahnschaffe, A., et al. (2018). High-accuracy determination of internal circadian time from a single blood sample. J. Clin. Invest. 128, 3826-3839. https://doi.org/ 10.1172/JCI120874.
  64. Newman, A.M., Steen, C.B., Liu, C.L., Gentles, A.J., Chaudhuri, A.A., Scherer, F., Khodadoust, M.S., Esfahani, M.S., Luca, B.A., Steiner, D., et al. (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 37, 773-782. https://doi. org/10.1038/s41587-019-0114-2.

STAR*METHODS

KEY RESOURCES TABLE

REAGENT or RESOURCE SOURCE IDENTIFIER
Antibodies
Anti-mouse CD45 clone 30-F11 BD Biosciences 564279/AB_2651134
Anti-mouse CD3e clone KT3.1.1 BioLegend 155617/AB_2832541
Anti-mouse CD4 clone GK1.5 BD Biosciences 563232/AB_2738083
Anti-mouse CD8a clone 53-6.7 BD Biosciences 563152/AB_2738030
Anti-mouse CD11c clone N418 BioLegend 117307/AB_313776
Anti-mouse CD19 clone 1D3 BD Biosciences 566412/AB_2744315
Anti-mouse NK1.1 clone PK136 BioLegend 108715/AB_493591
Anti-mouse MHCII clone M5/114.15.2 BioLegend 107643/AB_2565976
Anti-mouse Ly6G clone 1A8 BioLegend 127645/AB_2566317
Anti-mouse Ly6C clone HK1.4 BioLegend 128023/AB_10640119
Anti-mouse CD137 clone 17B5 BioLegend 106106/AB_2287565
Anti-mouse CD107a clone 1D4B BioLegend 121610/AB_571991
Anti-mouse PD-1clone 29F.1A12 BioLegend 135213/AB_10689633
Anti-mouse Granzyme B clone QA16A02 BioLegend 372207/AB_2687031
Anti-mouse IFN- clone XMG1.2 BioLegend 505837/AB_11219004
Anti-mouse CD31 clone MEC13.3, BioLegend 102520/AB_2563319
Anti-mouse CD4 clone, GK1.5 BioLegend 100405/AB_312690
Anti-mouse CD54 clone, YN1/1.7.4 BioLegend 116114/AB_493495
Anti-mouse CD8 clone 53-6.7 BioLegend 100724/AB_389326
Anti-mouse VCAM-1 clone 429 BioLegend 105710/AB_493427
Anti-mouse anti-CD11a clone M17/4 BioXCell BE0006/ AB_1107578
Anti-mouse anti-CD18 clone M18/2 BioXCell BE0009/ AB_1107607
Anti-mouse PD-1, clone 29F.1A12 BioXCell BE0273/ AB_2687796
Anti-mouse CD8a, clone YTS 169.4 BioXCell BE0117/ AB_10950145
Anti-mouse E-selectin clone UZ6 Invitrogen MA1-06506/ AB_2186701
Anti-human CD3 clone BW264/56 Miltenyi Biotec 130-113-687/AB_2726228
Anti-human CD8 clone RPA-T8 BD Biosciences 563795/AB_2722501
Anti-human CD4 clone SK3 BD Biosciences 565995/AB_2739446
Anti-human CD19 clone HIB19 BioLegend 302219/AB_389313
Anti-human CD4 clone EP204 Zhongshan Golden Bridge Biotechnology ZA-0519
Anti-human CD8 clone SP16 Zhongshan Golden Bridge Biotechnology ZA-0508
rat IgG2a isotype control, clone 2A3 BioXCell BE0089/ AB_1107769
rat IgG2b isotype control, clone LTF-2 BioXCell BE0090/ AB_1107780
Dynabeads human T-Activator CD3/CD28 Gibco 11132D/AB_2916088
Anti-human CD3 clone BW264/56 BioLegend 317348/ AB_2571995
Dynabeads Mouse T-Activator CD3/ CD28 for T-Cell Expansion and Activation Gibco 11453D
Goat anti-Rat IgG (H+L) Cross-Adsorbed Secondary Antibody Invitrogen A21247/ AB_141778
Chemicals, peptides, and recombinant proteins
DRAQ7 Biolegend 424001
Counting Beads ThermoFisher C36950
(Continued on next page)
Continued
REAGENT or RESOURCE SOURCE IDENTIFIER
eBioscience Fixable Viability Dye eFluor eBioscience 65-0865-18
Foxp3 / Transcription Factor Staining Buffer Set eBioscience 00-5523-00
Phorbol 12-myristate 13-acetate Sigma-Aldrich P1585
Ionomycin Sigma-Aldrich I3909
GolgiPlug BD Biosciences 555029
CellEvent Caspase-3/7 Green Ready Flow Reagent Invitrogen R37167
CellTracker Deep Red Dye Invitrogen C34565
RBC lysis buffer Biolegend 420302
EdU Alexa Fluor Flow Cytometry Assay Kit Invitrogen C10646
Trizol Reagent Invitrogen 15596018
Human IL-2 PeproTech 200-02
Mouse IL-2 Biolegend 575406
retronectin Takara T100B
streptavidin Biolegend 405207
biotinylated protein L ThermoFisher 29997
Collagenase IV Worthington Biochemical Corporation LS004189
Collagenase D Roche 11088866001
DNase I Roche 04716728001
DAPI Biolegend 422801
BlockAid Blocking Solution Invitrogen B10710
Bond Polymer Refine Red Detection kit Leica DS9390
Critical commercial assays
Chromium Single Cell 3′ v3.1 Reagent Kit with dual indexes 10x Genomics 1000269
CD8+ T Cell Isolation Kit, human Miltenyi Biotec 130-096-495
CD8a+ T Cell Isolation Kit, mouse Miltenyi Biotec 130-104-075
Deposited data
Timed scRNA-seq of TILs from murine melanoma This paper GEO: GSE260641
Mouse epidermal scRNA-seq Duan et al. GEO: GSE223109
Mouse skin RNA-seq Welz et al. GEO: GSE115104, GSE114943
Mouse skin RNA-seq Wang et al. GEO: GSE83855
Mouse epidermal RNA-seq in constant darkness Tsujihana et al. GEO: GSE174155
Mouse skin microarray Geyfman et al. GEO: GSE38622, GSE38623
Human epidermal suction blister microarray Spörl et al. GEO: GSE35635
Human skin microarray del Olmo et al. GEO: GSE205155
Human forearm skin microarray Wu et al. GEO: GSE112660
Human skin microarray (fractional laser treatment) Sherrill et al. GEO: GSE139305
Human melanoma scRNA-seq Sade-Feldman et al. GEO: GSE120575
Other supporting data This paper (https://doi.org/10.26037/yareta: 5rmtImd5sreo7bmcbe44bbijfm).
Experimental models: Cell lines
B16-F10-OVA Stéphanie Hugues, UNIGE, Switzerland N/A
Continued
REAGENT or RESOURCE SOURCE IDENTIFIER
Murine MC38 colorectal carcinoma cell line Mark J. Smyth RRID: CVCL_B288
Human lymphoma DoHH2 Francesco Bertoni, IOR,Bellinzona, Switzerland N/A
Experimental models: Organisms/strains
Mouse: WT C57BL/6J The Jackson Laboratory Strain # 000664 RRID:IMSR_JAX:000664
Mouse: WT C57BL/6N Charles River Strain # CRL027 RRID:IMSR_CRL:027
Mouse: B6.129S4(Cg)-Bmal1tm1Weit/J The Jackson Laboratory
Strain #:007668
RRID:IMSR_JAX:007668
Mouse: B6.Cg-Tg(Cd4-cre)1Cwi/BfluJ The Jackson Laboratory
Strain #:022071
RRID:IMSR_JAX:022071
Cdh5-cre/ERT2 mice Ralf Adams, MPI Munster, Germany N/A
NSG (NOD.Cg-PrkdcscidIL2rgtm1Wjl/SzJ) Charles River Strain #:614
BRaf , Pten , Tyr::CreERT2 Ping-Chih Ho, UNIL, Switzerland He et al.
CD45.1 OT-I Walter Reith, UNIGE, Switzerland N/A
Oligonucleotides
Mouse Pdcd1 F
CGGTTTCAAGGCATGGTCATTGG
Microsynth N/A
Mouse Pdcd1 R
TCAGAGTGTCGTCCTTGCTTCC
Microsynth N/A
Mouse Rplp0 F
GAGACTGAGTACACCTTCCCAC
Microsynth N/A
Mouse Rplp0 R
ATGCAGATGGATCAGCCAGG
Microsynth N/A
Mouse Rpl32 F
ACAATGTCAAGGAGCTGGAG
Microsynth N/A
Mouse Rpl32 R
TTGGGATTGGTGACTCTGATG
Microsynth N/A
Human PDCD1-F
CCAGGATGGTTCTTAGACTCCC
Microsynth N/A
Human PDCD1-R
TTTAGCACGAAGCTCTCCGAT
Microsynth N/A
Human Rplp0 F
TGGTCATCCAGCAGGTGTTCGA
Microsynth N/A
Human Rplp0
ACAGACACTGGCAACATTGCGG
Microsynth N/A
Software and algorithms
FACSDiva 6 BD Biosciences RRID:SCR_001456
FlowJo 10 BD Biosciences RRID:SCR_008520
Slidebook 3 i RRID:SCR_014300
ImageJ FIJI RRID:SCR_002285
Cell Ranger 6.1.2 10X Genomics RRID:SCR_017344
Seurat 4.3.0 Hao et al. https://satijalab.org/seurat/
R 4.2.2 The R Foundation https://www.r-project.org/
clusterProfiler Wu et al. https://bioconductor.org/packages/ release/bioc/html/clusterProfiler.html
Cosinor Open source https://github.com/sachsmc/cosinor
DiscoRhythm Carlucci et al. https://bioconductor.org/packages/ release/bioc/html/DiscoRhythm.html
GSVA Hanzelmann et al. https://bioconductor.org/packages/ release/bioc/html/GSVA.html
(Continued on next page)
Continued
REAGENT or RESOURCE SOURCE IDENTIFIER
TCGAbiolinks Colaprico et al. https://bioconductor.org/packages/ release/bioc/html/TCGAbiolinks.html
zeitZeiger Hughey et al. https://github.com/hugheylab/zeitzeiger
sva Leek et al. https://bioconductor.org/packages/ release/bioc/html/sva.html
Analysis code This study https://github.com/zqun1/circadian_ immune_mouse_melanoma

RESOURCE AVAILABILITY

Lead contact

Further information and requests for resources or reagents should be directed to and will be fulfilled by the lead contact, Christoph Scheiermann (christoph.scheiermann@unige.ch).

Materials availability

This study did not generate new unique reagents.

Data and code availability

  • Single-cell RNA-seq data have been deposited at GEO: GSE260641 and are publicly available as of the date of publication. Accession numbers of public data used in this study are listed in the key resources table. All data supporting the conclusions of this paper are available online (https://doi.org/10.26037/yareta:5rmtlmd5sreo7bmcbe44bbijfm).
  • All original code has been deposited at GitHub https://github.com/zqun1/circadian_immune_mouse_melanoma and is publicly available as of the date of publication.
  • Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

EXPERIMENTAL MODEL AND STUDY PARTICIPANT DETAILS

Animals

C57BL/6 and NSG mice were purchased from Charles River or the Jackson Laboratory. Transgenic mouse lines Bmal1 (purchased from Jackson Labs) and Cdh5 (gift from Dr. Ralf Adams, Max-Planck-Institute for Molecular Biomedicine Münster, Germany) were crossed and bred at ENVIGO. Transgenic mice were maintained as homozygous for Bmal and heterozygous for the relevant Cre. CD45.1 OT-I (gift from Walter Reith) mice were bred in house. All mice used were at 8-12 weeks of age. Cre recombination in Cdh5 mice was induced by intraperitoneal injection of tamoxifen ( 1 mg per injection) for 5 consecutive days. Unless specified, mice were housed under a light:dark schedule with food and water ad libitum. When multiple time points were investigated simultaneously, light-tight cabinets (Techniplast) were used to shift animals to the respective phase prior to the experiments. Treatment times correspond to Zeitgeber time (ZT) and indicate timing relative to lights on in the animal facility such that ZT1 is 1 h after lights on (morning), and ZT13 is 1 h after lights off (evening). A spontaneous tumor model was induced in BRaf , Pten , Tyr::CreERT2 mice as previously described by topically treatment with -hydroxytamoxifen ( in ethanol) on the skin surface. During the dark phase, all manipulations were performed under dim red light. All animal procedures and experiments were approved and performed in accordance with the guidelines of the animal research committee of Geneva, Switzerland, the veterinary authority of Canton de Vaud, Switzerland, or the MGH Institutional Animal Care and Use Committee (IACUC) and were performed in accordance with MGH IACUC regulations.

METHOD DETAILS

Tumor cell lines and inoculation

Mouse melanoma cell lines B16-F10 (ATCC) and B16-F10-OVA (gift from Stéphanie Hugues, University of Geneva, Switzerland) were maintained in RPMI (Gibco) supplemented with heat-inactivated FCS (Gibco), penicillin-streptomycin (Gibco), and of -mercaptoethanol (Gibco). MC38 murine colon adenocarcinoma cells (gift from Mikaël Pittet, University of Geneva, Switzerland) were maintained in DMEM (Gibco), heat-inactivated FCS and penicillin-streptomycin). Cell lines were used by passage 10 and tested negative for Mycoplasma. Unless otherwise specified, tumor cells in PBS were injected s.c. into the right flank of mice, under isoflurane anesthesia. For anti-PD-1 therapy experiment, tumor cells in PBS or in PBS were injected s.c. Human lymphoma DoHH2 cell line (gift from Francesco Bertoni, Institute
of Oncology Research (IOR), Bellinzona, Switzerland) were cultured in RPMI (Gibco) supplemented with heat-inactivated FCS (Gibco). A total of cells were injected s.c. into the right flank of NSG mice. Tumor volume was monitored every 1 to 2 days using a caliper and calculated by length width width .

Mouse T cell activation

Splenocytes from C57BL/6 wild type mice, Cd4cre:Bmal1 and their littermate controls, or OT-I mice were collected by smashing the spleen through a mesh cell strainer. Red blood cells were lysed using RBC lysis buffer (Biolegend, 420302) for 5 min on ice. To activate cells from C57BL/6 or Cd4cre:Bmal1 mice, cells were enriched using a cell isolation kit (Miltenyi) according to the manufacturer’s instructions. The isolated T cells were counted and resuspended in a concentration of in T cell complete medium (RPMI, 10% heat-inactivated FCS, 2 mM L-glutamine, penicillin-streptomycin, -mercaptoethanol, 1 mM pyruvate sodium) supplemented with mouse interleukin-2 (IL-2) (Biolegend) in 24-well plate. Dynabeads Mouse T-Activator CD3/CD28 beads (Gibco, 11452D) were used to activate T cells according to the manufacturer’s instructions. For OT-I cell activation, whole splenocytes were resuspended in a concentration of in T cell complete medium (RPMI, heatinactivated FCS, 2 mM L-glutamine, 1% penicillin-streptomycin, -mercaptoethanol, 1 mM pyruvate sodium) supplemented with mouse interleukin-2 (IL-2) (Biolegend) and 1 nM of OVA257-264 peptide (SIINFEKL, InvivoGen) in a 6-well plate. After 3 days of activation, cells were collected and purified by removing Dynabeads (for polyclonal cells) or using a cell isolation kit (Miltenyi) (for OT-I). For in vivo experiments, OT-I cells were injected intravenously (i.v) into tumor-bearing mice.

Mouse T cell synchronization

Cells were synchronized using horse serum as previously described. In brief, activated T cells (as described above) were collected and seeded into 24 -well plate with a density of in T cell complete medium, supplemented with mouse interleukin-2 (IL-2) (Biolegend) and anti-mouse CD3 antibody (Invitrogen, 16-0031-85) to maintain the TCR signal. An equal volume of horse serum (Sigma, h1270) was pre-warmed and added directly to the plate (serum shock). After 2 h incubation at with CO2, cells were washed and resuspended in T cell complete medium with IL-2 and anti-CD3 antibody.

Human T cell activation and synchronization

Human T cells were isolated using T Cell Isolation Kit (Miltenyi 130-096-495) from the peripheral blood mononuclear cells from healthy donors’ buffy coat (provided by University Hospitals of Geneva (HUG), Switzerland) using Ficoll-Paque Plus (Cytiva). Cells were resuspended at and activated for 3 days with Dynabeads human T-Activator CD3/CD28 (Gibco, 11132D) in the presence of recombinant human IL-2 ( , PeproTech, 200-02) in T cell complete medium at with . Cells were synchronized using horse serum as described above, supplemented with human IL- 2 and anti-human CD3 antibody (Invitrogen, Clone OKT3, 16-0037-81).

Luciferase-expressing cells

cells were isolated from the spleen of Per2::Luc mice and activated as described above. Activated CD8 cells were synchronized and resuspended at in 2 ml T cell complete medium containing luciferin (Abcam, ab45164). Dishes were sealed with parafilm and recorded using a LumiCycle detector (Actimetrics).

Tissue digestion and single-cell preparation

Tumor tissue was collected and chopped into small pieces, digestion medium was added (RPMI containing collagenase IV (Worthington Biochemical Corporation), DNase I (Roche 04716728001) and heat-inactivated FCS) and incubated for 30 minutes at . Cells were rinsed through a cell strainer to obtain single-cell suspensions. To remove debris, Lympho-lyte(R)-M (CEDARLANE) was used according to manufacturer’s instructions.

Flow Cytometry

Single-cell suspensions were prepared and incubated with mouse or human Fc receptor block (anti-mouse CD16/32 Biolegend, human FcR blocking reagent, Miltenyi Biotec) for 10 minutes at room temperature (RT). After incubation, unless specified otherwise, the antibody mix was added directly into the cell suspension and incubated for 15 min at .
Cells were washed and resuspended in FACs buffer with viability dye (DRAQ7, Biolegend, ) and characterized using an 18-colour BD LSR Fortessa (BD Biosciences) or Beckman Coulter Cytoflex. Acquired data were analyzed using FACSDiva 6 (BD Biosciences) and FlowJo 10 (BD). Cell counts were calculated using Counting Beads (C36950, C36995, ThermoFisher).
For intracellular staining, cells were stained with viability dye (eBioscience Fixable Viability Dye eFluor 780, 65-0865-18), followed by surface staining as previously described. Cells were fixed and permeabilized using Foxp3 / Transcription Factor Staining Buffer Set (eBioscience, 00-5523-00). Upon wash with permeabilization buffer, the intracellular antibody was added and incubated for 30 min at room temperature. The following antibodies were used for intracellular staining: anti-mouse Granzyme and anti-mouse IFN- . For cytokine detection (Granzyme B, IFN- ), cells were stimulated before staining in complete RPMI medium with phorbol 12-myristate 13-acetate (PMA, ; Sigma-Aldrich), ionomycin ( ; Sigma-Aldrich), and BD GolgiPlug ( ; BD
Biosciences, 555029 ) and incubated 4 hours at CO2. Caspase- detection was performed using CellEvent Caspase-3/ 7 Green Ready Flow Reagent (Invitrogen R37167) according to the manufacturer’s instruction.

Mouse leukocyte adoptive transfer

Donor cells were collected from bone marrow (flushed with PBS) and spleen (smashing through a cell strainer) from donor mice and lysed with RBC lysis buffer ( 420302 Biolegend). Cells were labeled with CellTracker Deep Red Dye ( in PBS, ), and injected i.v. ( bone marrow cells and splenocytes). Unless specified, tissues were harvested after 2 h of cell injection. For egress assays, a total of cells ( of bone marrow cells and of splenocytes) in PBS were injected intratumorally. Tissues were harvested 24 h after the injection.

EdU assay

EdU ( 1 mg per mouse, Click-iT Plus EdU Alexa Fluor Flow Cytometry Assay Kit, Invitrogen, C10646) was injected intraperitoneally into mice 4 h before tissue collection. Tumor tissue was collected and digested as described above. Cells were first stained with viability dye (eBioscience Fixable Viability Dye eFluor 780, 65-0865-18), followed by surface staining as previous described. Then cells were fixed, permeabilized, and EdU detected according to the manufacturer’s instruction.

CD8 T cell sorting and RNA extraction

To obtain T cells, tumor tissues were digested and processed as described above. Tumor-infiltrating leukocytes were then stained using anti-mouse CD45 (clone 30-F11 BUV 395), CD3e (clone KT3.1.1, BV421), CD4 (clone GK1.5 BV650), CD8a (clone 53-6.7 BV605), CD19 (clone 1D3, BB700), and NK1.1 (clone PK136, PE/Cy5) together with DRAQ7, and CD45 cells were sorted using a BD Aria II. Flow cytometry sorted CD8 T cells were collected in RNAprotect Cell Reagent (cat. #76526, Qiagen). RNA was isolated using an RNeasy Plus Micro Kit (cat. #74034, Qiagen) according to the manufacturer’s instructions. RNA integrity and quantity were assessed with a Bioanalyzer (Agilent Technologies) or a Nanodrop 2000 (ThermoFisher).

RNA extraction, reverse transcription and qPCR

Cells were collected and lysed using Trizol Reagent. RNA quantity and quality was analyzed using a Nanodrop 2000 (ThermoFisher) or Bioanalyzer. Reverse transcription was performed using PrimeScript RT Reagent Kit (Takara) according to the provided instructions. Q-PCR analyses were performed using PowerUp SYBR Green (Applied Biosystems). Quantification of the transcript was performed using the method using Rplp0 or Rpl32 as internal reference genes.

In vivo antibody treatments

To block LFA-1, anti-CD11a (clone M17/4, 100 ) and anti-CD18 (clone M18/2, ) antibodies were administered intravenously (i.v.). Anti-PD-1 antibody (anti-mouse PD-1, clone ) was administrated intraperitoneally (i.p) at day 4 or 6 (B16-F10OVA) or day 6 or 7 (MC-38) after tumor engraftment. To deplete CD8 T cells, depletion antibody (anti-mouse CD8a, clone YTS 169.4, ) was injected i.p. at the same day of PD-1 administration and repeated every 2 days. Isotype control rat IgG2a (clone 2A3) and rat IgG2b (clone LTF-2) antibodies were used. All antibodies were purchased from BioXCell.

Human CAR T cell generation

Human peripheral blood mononuclear cells (PBMCs) were collected from healthy donors’ buffy coat (provided by University Hospitals of Geneva (HUG), Switzerland) using Ficoll-Paque Plus (Cytiva), then frozen in aliquots and stored in liquid nitrogen. Cryopreserved PBMCs were thawed and resuspended at cells per ml , activated for 3 days with Dynabeads human T-Activator CD3/CD28 (Gibco, 11132D) in the presence of recombinant human IL-2 ( , PeproTech, 200-02) in RPMI1640 (Gibco) media supplemented with heat-inactivated FBS (Gibco) and penicillin-streptomycin at with . Activated T cells were transduced by culturing them for 48 hours in retronectin (Takara, T100B)-coated 24-well plates with retrovirus encoding second generation anti-CD19-CD28-CD3ఢ CAR. The CAR construct was kindly provided by Dr. Crystal Mackall (Stanford University). Dynabeads were removed, and after washing, cells were cultured in complete media with IL-2. At day 14, cells were harvested and frozen until further use. Transduction efficacy of CAR expression on T cells was evaluated by APC-conjugated streptavidin (Biolegend, 405207) binding to biotinylated protein (ThermoFisher, 29997). Before use, CAR cells were thawed, washed and rested for 4 hours at .

Human CAR T cell adoptive transfer

Human DoHH2 lymphoma cells were engrafted as described above. For long-term treatment experiments, CAR T cells or equal number of untransduced cells were injected i.v. into tumor-bearing mice at ZT1 or ZT13 on day 12 after tumor engraftment. For short-term homing experiment, CAR T cells were injected i.v. at ZT1 or ZT13 on day 15 after tumor engraftment. Twenty-four hours post transfer, tumors were harvested and digested using digestion medium (RPMI containing collagenase IV (Worthington Biochemical Corporation), collagenase D (Roche, 11088866001), DNase I (Roche 04716728001 ) and heat-inactivated FCS) and incubated for 30 minutes at . Single cell suspension were analyzed as described above.

Immunofluorescence imaging

Tumor tissues were collected and embedded into OCT blocks (CellPath) and kept at . Tumors were sectioned into sections using a Cryostat. Sections were postfixed with PFA for 10 min at room temperature. Following three washes with PBS, they were incubated with BlockAid Blocking Solution (Invitrogen , B10710) for 2 h at RT. After PBS, sections were stained with an antibody diluted in Blocking Solution at overnight. The following antibodies were used, AF594 anti-mouse CD31 (clone MEC13.3, Biolegend 102520), FITC anti-mouse CD4 (clone, GK1.5, Biolegend 100405), AF647 anti-mouse CD54 (clone, YN1/1.7.4, Biolegend 116114), AF647 anti-mouse CD8 (clone 53-6.7, Biolegend 100724), AF488 anti mouse VCAM-1 (clone 429, Biolegend 105710), E-selectin (clone UZ6, Invitrogen MA1-06506), Goat anti-Rat IgG (H+L) Cross-Adsorbed Secondary Antibody, Alexa Fluor (Invitrogen A21247). Slides were washed 3 times using PBS and stained with DAPI ( 300 nM , Biolegend 422801 ) for 5 min at room temperature. Images were obtained as sections using a Zeiss Axio Examiner.Z1 confocal spinning disk microscope equipped with 405-, 488-, 561- and 640-nm laser sources. All image analyses were performed in ImageJ and Slidebook (Intelligent Imaging Innovations, 3i).

Human immunohistochemistry

Formalin-fixed, paraffin-embedded (FFPE) sections of melanoma were stained by immunohistochemistry using anti-CD4 (EP204) or anti-CD8 (SP16) (Zhongshan Golden Bridge Biotechnology, Beijing, China) primary antibodies. All staining was performed on the Leica Bond III automated platform (Leica Biosystems, Newcastle, UK) using Bond Polymer Refine Red Detection kit (DS9390). Slides were scanned with a Pannoramic P250 FLASH slide scanner (3DHistech, Budapest, Hungary).

Single cell sequencing and analysis

Tumors were harvested at different times of the day and processed as described above. Tumor infiltrating leukocytes were sorted (CD45 DRAQ7-) using a Biorad S3 sorter in PBS with BSA. Cellular suspension was loaded on a Genomics Chromium instrument. Single-cell RNA-Seq libraries were prepared using Chromium Single Cell 3′ v3.1 Reagent Kit with dual indexes according to manufacturer’s protocol. Library quantification and quality assessment was performed using a Qubit fluorometer (ThermoFisher Scientific) and a Tapestation (DNA High sensitivity chip – Agilent Technologies). Libraries were sequenced on an Illumina NovaSeq 6000 using paired-end bp as sequencing mode. In total, 30693 single cells were sequenced (ZT1: 9764; ZT7: 7933; ZT13: 6255; ZT19: 6651).
Raw count matrices were generated with Cell Ranger v6.1.2 with mouse genome build mm 10 as reference. Genes expressed in less than 20 cells were removed. Cells expressing less than 200 features, more than 7500 features, or more than of mitochondrial genes were removed from further analysis. Dimensionality reduction and clustering analysis were done using Seurat v4. Cells were annotated by conventional cell markers: B (Cd79a, Igkc), cDC1 (Clec9a, Xcr1), cDC2 (H2-Aa, H2-Ab1, Itgax, Clec10a), activated cDC (Ccr7, Fscn1), pDC (Siglech, Bst2), macrophage (Cd68, Itgam), monocyte (Ly6c2, Ccr2), neutrophil (G0s2, Csf3r), NK (Nkg7, Klrb1c), (Trdc, Tcrg-C1), dendritic epidermal T cells (Trdv4, Tcrg-C1), CD4 (Trac, Cd4), CD8 (Trac, Cd8a), Treg (Foxp3, Ctla4), melanocyte (Mlana, Dct), and mitotic cells (Birc5, Top2a). ‘Anti-tumorigenic lymphocytes’ include NK, , NKT, CD8 effector, CD4 naïve, CD8 effector memory, CD8 mitotic, and CD4 effector clusters, whereas ‘Pro-tumorigenic lymphocytes’ include Treg, Treg mitotic and CD8 exhausted clusters. Further annotation of CD4 and CD8 T cell subclusters cells were based on exhaustion (Tox, Lag3, Pdcd1, Havcr2), effector (Gzmb, Prf1, Ifng, II2, II4), memory (Gzmk, Itga4), and memory/ naïve (Tcf7, Sell, Ccr7, Lef1) markers. Oscillatory genes were identified by Discorhythm with each cell as one read. Gene set over-representation analyses were performed with clusterProfiler and gene signature scores were calculated by GSVA. Survival analyses were performed by using the TCGAbiolinks package. Cosinor fitting of gene expression was performed by using the cosinor package.

ZeitZeiger analysis

The ZeitZeiger model was trained using mouse melanoma scRNAseq data from the present study, mouse skin scRNAseq data from in-house, mouse skin scRNAseq data from GSE223109, and mouse skin bulk RNAseq data from GSE115104 and GSE83855 data sets. To balance the representation of scRNAseq vs RNAseq and melanoma vs healthy skin data, the melanoma scRNAseq dataset was randomly sampled by 60% four times and the in-house skin scRNAseq dataset was randomly sampled by two times to form the training set. For scRNAseq datasets, pseudobulk expression in each mouse was calculated and multiplied by 100 to generate the transcript per million as in data for training ZeitZeiger. Each dataset was randomly divided into a training set and a testing set at a ratio of . Training data and testing data were batch effect corrected by combat from the sva package. Conserved coefficients were summarized from 20 runs of cross validation with different random seeds and different sumabsv and nSpc parameters. The genes KIf9, Nr1d1, Cry2, Per1, Gapdh, Ppia, and Psmb2 were manually added to the coefficients, as previously described. Further validation of the trained ZeitZeiger model was performed using RNAseq or RMA-normalized microarray data from both mouse (GSE114943, GSE174155, GSE38622, GSE38623 ) and human (GSE35635, GSE205155, GSE112660 and GSE139305 ). To correct for the diurnal and nocturnal behavioral difference between human and mouse, predicted human sampling time were generated by adding 12 h .

Estimation of cell composition in TCGA melanoma patients

To estimate the ratio of exhausted vs non-exhausted cells, as well as major immune subsets, cells from GSE120575 were annotated as described earlier and the top 20 differentially expressed genes were used to estimate the immune cell composition by cibersortx

QUANTIFICATION AND STATISTICAL ANALYSIS

Unless specified, all data were plotted from independent biological replicates. Data was analyzed using Prism 9 and Prism 10 (GraphPad). .

Supplemental figures

Figure S1. Time of day of harvest dictates numbers of tumor-infiltrating leukocytes, related to Figure 1
(A) Gating strategy of leukocytes in B16-F10-OVA and Tyr::Cre ; Pten ; BRaf tumors.
(B) Normalized total cell numbers of leukocyte subsets in B16-F10-OVA tumors, harvested at 4 different times of the day (zeitgeber time [ZT]); mice from 4 independent experiments, cosinor analysis.
(C) Blood leukocyte numbers sampled at 4 different times of the day from B16-F10-OVA tumor-bearing mice, mice from 5 independent experiments, cosinor analysis.
(D) Tumor volume of B16-F10-OVA tumors, measured 3 times a day, mice.
(E) Quantification of the CD31 area in B16-F10-OVA tumors, harvested at ZT1 or ZT13, mice, one-way ANOVA, and unpaired Student’s test.
(F) Normalized numbers of blood leukocyte sampled at ZT1 or ZT13 from Tyr::Cre ; Pten ; BRaf tumor-bearing mice, mice from 2 independent experiments, paired Student’s t test. All data are represented as mean SEM, ns, not significant. *** ; **** .
Figure S2. Circadian leukocyte infiltration and T cell therapy, related to Figure 2
(A) Gating strategy and quantification of EdU leukocytes in B16-F10-OVA tumors, mice, unpaired Student’s t test.
(B) Gating strategy and quantification of caspase-3/-7 leukocytes in B16-F10-OVA tumors, harvested at ZT1 ( mice) or ZT13 ( mice), from 2 independent experiments, unpaired Student’s test.
(C) Normalized numbers of remaining cells in tumors 24 h post intratumoral injection of leukocytes at ZT1 ( mice) or ZT13 ( mice), from 2 independent experiments, unpaired Student’s test.
(D) Tumor volume after anti-LFA-1 treatment ( mice), isotype control ( mice), from 2 independent experiments, two-way ANOVA.
(E) Quantification of VCAM-1 ( sections from 3 mice) or E-selectin ( sections from 4 mice) expression on CD31 cells in B16-F10-OVA tumors, harvested at ZT1 or ZT13, from 2 independent experiments, unpaired Student’s t test. All data are represented as mean SEM, ns, not significant.
(F) Tumor volume after i.v. injection of activated OT-I cells at ZT1 ( mice) or ZT13 ( mice), from 3 independent experiments.
Figure S3. Identification of TIL subsets, related to Figure 3
(A-C) Feature plot and heatmap of genes to identify each leukocyte cluster obtained by scRNA-seq of B16-F10-OVA tumors, harvested at 4 different times of the day.
Cell
Article
OPEN ACCESS
Figure S4. Diurnal oscillating TILs, related to Figure 3
(A) UMAP representation of scRNA-seq analyses of TILs harvested at 4 different times of the day.
(B) Relative abundance of immune cells per cluster across different times of the day.
(C) Expression of the circadian gene Per1 in each leukocyte cluster obtained by scRNA-seq of B16-F10-OVA tumors, harvested at 4 different times of the day.
(D) Number of significantly oscillating genes detected in each cluster using cosinor (CS) and JTK cycle analyses.
(E) Top oscillatory Gene Ontology pathways common in the different immune cell clusters.
(F) Heatmap of oscillating genes of each cluster. Circadian genes are highlighted.
Figure S5. Diurnal oscillating genes in TILs, related to Figure 3
(A and B) Feature plot and dot plot of genes to identify different T cell clusters obtained by scRNA-seq of B16-F10-OVA tumors, harvested at 4 different times of the day.
(C) Relative abundance of each cell cluster across different times of the day.
(D) Gating strategy of exhausted and non-exhausted cells, pre-gated on CD8 T cells.
(E) Expression of Pdcd1 quantified by pseudobulk analysis of the CD8 cluster obtained by scRNA-seq of B16-F10-OVA tumors, harvested at 4 different times of the day. To reduce the sparseness of scRNA-seq data, cluster was randomly sampled by for 10 times with different seeds. The expression levels of Pdcd1 at 4 different times of the day were averaged within each pseudobulk cell.
Figure S6. cell phenotype oscillates in mice and humans, related to Figures 4, 5, and 6
(A) Schematic overview of cell synchronization in vitro.
(B) Tumor volume after anti-PD-1 treatment administered at ZT1 or ZT13 in B16-F10-OVA tumors; mice; control mice, from 3 independent experiments.