$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Bu analiz hattının potansiyelini göstermek için, fonksiyonel manyetik rezonans görüntüleme (fMRI) taramaları sırasında farelere işitsel uyaranlar sunulurken, kafaya sabitlenmiş dişi farelerin - özellikle anneler ve bakirelerin - davranışsal videoları elde edildi. Uyaranlar, yavru çağrıları ve pasif olarak sunulan saf tonlardan oluşuyordu, bu nedenle herhangi bir talimat verilmiş davranışsal çıktı olmadan. Yavru çağrıları, yuvalarından geçici olarak izole edilen 6 günlük fare yavruları tarafından yayılan ultrasonik seslendirmelerin kayıtlarıydı. Bu yavru çağrıları tipik olarak, annenin yerini bulduğu, doğru hareket ettiği, araştırdığı, daha sonra yavruyu aldığı ve yuvanın güvenliğine geri götürdüğü anne yavru alma eylemini ortaya çıkarır - bakire dişilerin tipik olarak sergilemediği bir davranış11. Önceki çalışmalar, yavru çağrılarının annelerin birincil işitsel korteksinde güçlü uyarılmış aktivite ortaya çıkardığını, ancak bakirelerin ortaya çıkmadığını, saf ton tepkilerinin ise grup farkı göstermediğini göstermiştir12. Bu nedenle, yavru çağrılarına verilen davranışsal tepkilerdeki grup farklılıklarının, ancak saf tonların değil, tanımlanacağı varsayılmıştır. Bununla birlikte, kafası sabit bir ortamda yavru ipuçlarına yanıt olarak dişi farelerin spesifik davranışları hakkında çok az literatür vardır. Bu nedenle, belirli bir davranışsal çıktı beklenmemiştir, bu da bu deneyi önerilen veriye dayalı davranışsal analizler için mükemmel bir test haline getirmiştir.
Sekiz gün boyunca, tüm hayvanlar yavaş yavaş deneyci kullanımına, kafa fiksasyonuna ve MRI ortamına alıştırıldı. Baş fiksasyonuna ve deneysel ortama alışma, sunulan uyaranlara verilen davranışsal tepkileri değerlendirmek için çok önemlidir. Çevreye uygun şekilde alıştırılmazlarsa, hayvanlar yalnızca stres tepkileri gösterebilir ve aksi takdirde videografi yoluyla ayrıştırılabilecek uyaran güdümlü etkileri yıkayabilir.
Davranışsal cihaz kurulumu genellikle her veri toplama oturumu için aynı olsa da, bir hayvanın fMRI taramaları için kafası her sabitlendiğinde kamera görüş alanının (FOV) hafifçe kayması mümkündür (Şekil 1A'daki örnek FOV'lara bakınız). Bu muhtemelen kamera yuvasının konumundaki küçük değişikliklerin yanı sıra her hayvan için baş direği kafatası bağlantısındaki bireysel farklılıklardan kaynaklanıyordu. Bu nedenle, taramalardan elde edilen videoların, taramalar ve hayvanlar arasında içerdikleri uzamsal bilgilerin karşılaştırılmasına izin vermek için doğrusal ortak kayıt yoluyla birbiriyle uyumlu hale getirilmesi gerekiyordu. Ortak kayıt, temsili verilere göre uyarlandı. Örneğin, standart bir tarama gününde, hayvan başına 3-4 tarama yapıldı ve bu da hayvan başına günde 3-4 video ile sonuçlandı. Her tarama arasında, kamera ve kafa sabitleme bileşenleri de dahil olmak üzere beşik bileşenlerinde herhangi bir hareket yoktu. Böylece, hayvan başına günde bir video için ortak kayıt dönüşümü belirlendikten sonra, aynı gün o hayvanın diğer 2-3 videosuna uygulanabilir. Bu pipeline'ın ortak kayıt adımında zaman tasarrufu sağlanırken, gerekirse mekansal dönüşümler de her video için ayrı ayrı hesaplanabilir. Her videoda bulunan üç özellik, mekansal dönüşümü etiketlemek ve hesaplamak için nokta olarak seçildi. Gösterilen temsili verilerde, bu üç nokta baş direğinin merkezi, hayvanın sağ gözünün önden görünümü ve hayvanın sağ gözünün profil görünümüydü (45 derecelik bir açıyla yerleştirilmiş bir aynada görülebilir). Bu adımın etkisini göstermek için bu deney sırasında hem ortak kayıttan önce hem de sonra çekilen tüm videolardan hesaplanan ortalama kareyi gösteren Şekil 1A,B'de video ortak kayıt örnekleri gösterilmektedir.
Bu boru hattı, büyük ölçüde, ardışık kareler13 arasındaki görünür hızlarını tahmin ederek bir videodaki nesnelerin hareketini tahmin etmek için kullanılan bir bilgisayarla görme yöntemi olan optik akışa dayanır. Optik akış, önceden belirlenmiş herhangi bir vücut parçası veya ilgilenilen eylemler olmadan hareketin - davranışsal tepkilerin bir vekili - nicelleştirilmesine izin verdiği için seçildi. Ek olarak, bu yöntem, deney sırasında erişilebilen tek MR uyumlu kamera ile elde edilen temsili videoların sınırlı görüntü kalitesine uygundu. Bu boru hattında, küresel, yoğun optik akış tahmininin Horn-Schunck algoritması kullanılır; bununla birlikte, Lucas-Kanade algoritması gibi diğer algoritmalar, sağlanan MATLAB komut dosyalarında küçük değişikliklerle kolayca kullanılabilir 14,15,16. Şekil 2A, her piksel için üst üste bindirilmiş bağıl optik akış hızı vektörlerine sahip örnek bir video çerçevesini göstermektedir. Daha büyük vektörlerin, hayvanın burnu ve pençeleri gibi hareketin beklendiği alanlarda göründüğünü unutmayın.
Boru hattı, tam FOV'daki tüm videolar için optik akışı tahmin ederken, kalan analiz adımları, FOV içinde bir ilgi alanı (ROI) seçilerek hesaplama yükü ve veri depolama ihtiyaçları açısından çok daha yönetilebilir hale getirildi. Yatırım getirisi, deneye ve önceden belirlenmiş ilgilenilen davranışlara göre veya daha veriye dayalı bir yaklaşımla seçilebilir. Temsili verilerin a priori izlemek için belirli bir davranışsal çıktı içermediği göz önüne alındığında, veriye dayalı bir yaklaşım izlenmiştir. Optik akış büyüklüğünün standart sapması, Şekil 2B'de görselleştirildiği gibi, her pikseldeki tüm videolarda hesaplandı. En yüksek standart sapma alanları, göz çevresi ve burun gibi hayvanın kıvrımlarını içeriyordu ve bu da gözlemlenen optik akış dalgalanmalarının videolardaki gürültü tarafından değil, hayvan hareketi tarafından yönlendirildiğine dair güven veriyordu. Daha düşük standart sapma alanları, tarama sırasında kamera ve yuvada meydana gelen küçük titreşimlerin sonucu olabilecek beşiğin dış hatlarını içeriyordu. MRG ortamı, beşiğin yansıtıcı bileşenlerinde gözlemlenen optik akış dalgalanmalarında ortaya çıkabilen gradyanların değişmesi nedeniyle veri toplama sırasında kaçınılmaz olarak titreşimlerle doludur. Neyse ki, bu titreşimler veri toplama boyunca sabittir, bu nedenle uyaran durumundan bağımsızdır ve davranışsal analiz sonuçlarını etkilemesi beklenmez. Aynanın piksellerinde, hayvanın yüzünün profil görünümüne karşılık gelen ve temsili veriler için ROI seçimine rehberlik eden özellikle yüksek optik akış standart sapmaları gözlendi. Bu aynı zamanda, çırpma ve koklama gibi davranışsal tepkilerin, yavru çağrıları yayan izole bir yavru arayan dişi farelerin tipik davranış repertuarının bir parçası olarak beklenebileceği bir bölgeydi17.
Bir ROI seçtikten ve tüm videolar için optik akışın ortalama kare büyüklüğünü çıkardıktan sonra, ilgilenilen dönemlerdeki optik akış, gruplar ve koşullar arasında karşılaştırılabilir. Her videonun optik akış vektörünün ilk 20 saniyesinin (s) parlaklık sabitleme için maskelendiğini unutmayın; ancak bu, deneysel ihtiyaçlara uyacak şekilde ayarlanabilir. Her video için, seçilen ROI için kare bazında optik akış vektörü Z-skorlandı, ardından yüz ifadesi sınıflandırmasındaki önceki çalışmaların ardından, kameranın 30 Hz kare hızının beklenen herhangi bir davranışsal dalgalanmadan daha hızlı olmasını hesaba katmak için 5 hertz'de (Hz) düşük geçiren filtrelendi2. Son olarak, Z skorlu ve filtrelenmiş vektörler, uyaran sunumunun optik akış üzerindeki etkisini değerlendirmek için uyaran sunum dönemlerine ayrıldı. Her dönem için, optik akışı uyaran öncesi döneme normalleştirmek için ortalama uyaran öncesi temel sinyal çıkarıldı. Şekil 3A, Şekil 3B-E'deki özet grup verileriyle birlikte, bir anne ve bir bakire olmak üzere iki video için örnek optik akış zaman serilerini göstermektedir. Şekil 3B ve Şekil 3D, taban çizgisine göre uyaran sunumu sırasında zaman içindeki kümülatif optik akışı gösterirken, Şekil 3C ve Şekil 3E, uyaran başlangıcından 2.5 s sonraki kümülatif optik akışı özetlemektedir. Kümülatif optik akış, spontan davranışsal tepkilerin zaman kilitli bir şekilde gerçekleşeceğini varsaymadan, zaman içindeki genel hareketi yakalamak için hesaplandı. Genel olarak, bu temsili sonuçlar, tahmin edildiği gibi, yavru çağrılarının, ancak saf tonların değil, anne ve bakire dişi farelerden önemli ölçüde farklı davranışsal tepkiler uyandırdığını göstermektedir (grup arasında Mann-Whitney U testi: yavru çağrıları: p = 0.026; saf tonlar: p = 0.093). Bununla birlikte, bu uyaran etkisi 2 yönlü bir ANOVA'dan sağ çıkamazken, grup etkisi hayatta kaldı (uyaran: F(1,10) = 0.19, p = 0.67; grup: F(1,10) = 8.61, p = 0.015). Genel olarak, bu sonuçlar annelerin bakirelere kıyasla daha yüksek uyaranla uyarılmış hareket gösterdiğini ve annenin yavru çağrılarına verdiği tepkinin saf tonlardan daha tutarlı olduğunu göstermektedir. Bu, annelerde artan dikkat veya stresin yanı sıra, saf tonların aksine, doğal ortamlarda annelerde yavru geri alma tepkisini uyandıran yavru çağrılarının davranışsal alaka düzeyini yansıtabilir. Birlikte ele alındığında, bu temsili sonuçlar, dış uyaran sunumu sırasında optik akış tahmininin, nüanslı, spontan davranışsal tepkiler hakkında bilgi çıkarabileceğini göstermektedir.
Son olarak, temsili verilerde yakalanan davranışın mekansal özelliklerini ortaya çıkarmak için daha keşifsel bir analiz yapıldı. Uyaran sunumu sırasında ROI'deki hangi piksellerin koordineli bir şekilde dalgalandığını araştırmak için, zaman içinde uzamsal optik akış bilgileri üzerinde temel bileşen analizi (PCA) yapılmıştır. Bu analiz, Şekil 4'te gösterildiği gibi, ilk PC'ye en çok katkıda bulunan piksellerin yanı sıra her bir uyaran koşulu için grup farklılıkları gösteren pikselleri ortaya çıkardı. Şekil 4A,B'nin en sağdaki panelleri, her iki uyaran tipinin sunumu sırasında annelerin bakirelere göre burun içinde daha fazla hareket gösterdiğini göstermektedir. İki grup ve iki uyaran koşulu arasında, ilk PC, optik akış analizindeki toplam varyansın %5.41 ± %0.59'unu açıkladı. Analizimizin bu kısmı için ayna yatırım getirisi korunurken, gelecekteki analizler, uyaranlara yanıt olarak yüzün ötesindeki koordineli hareketleri karakterize etmek için FOV'un daha geniş bir bölümüne genişleyebilir. Örneğin, pençe hareketini karşılaştırmak, yavru çağrılarının tipik olarak annelerde yavru alımını başlattığı, ancak bakirelerde başlatmadığı ve pençelerin hayvanın kafasından daha özgürce hareket edebildiği göz önüne alındığında, daha önemli grup farklılıklarını ortaya çıkarabilir.
Şimdiye kadarki temsili sonuçlar, bu boru hattının kısıtlı videografi ortamlarında dış uyaranlara talimatsız davranışsal tepkileri değerlendirebileceğini öne sürse de, optik akışta gözlemlenen dalgalanmaların gerçekten anlamlı hayvan davranışını yansıtıp yansıtmadığına dair bir soru kaldı. Bu soruyu yanıtlamak için, aynı işlem hattı kullanılarak ayrı bir doğrulama veri kümesi analiz edildi. Ayrı bir deneyde, su kısıtlı erkek fareler, hafif bir işareti 6 μL ("yüksek ödül") veya 1 μL ("düşük ödül") su ödülünün verilmesiyle ilişkilendirmek üzere eğitildi. Özellikle, işitsel stimülasyon deneyinden farklı olarak, bu deneyin a priori davranışsal bir okuması vardı: yalama hızı. Su musluğuna yakın piksel parlaklığının analizi ile kolaylaştırılan, yalamaların video tabanlı tespiti yoluyla bir lickometre oluşturuldu. Bu nedenle, lickometer tarafından sağlanan davranışsal okuma, bu boru hattının optik akış tahmini tarafından sağlanan davranışsal okumasını karşılaştırmak için kullanılabilir ve spontan davranışları tespit etmedeki geçerliliğini destekler. Video ortak kaydı, optik akış tahmini, ROI seçimi (bir kez daha ayna FOV'u içeren) ve optik akış nicelemesinin ardından, optik akış büyüklüğünün karşılaştırılması, yüksek ve düşük ödüllere verilen davranışsal tepkiler arasında önemli bir fark olduğunu ortaya koydu. Bu analizin sonuçları Şekil 5A,B'de gösterilmiştir, burada Şekil 5A grup ortalaması alınmış zaman serisini gösterir ve Şekil 5B, uyaran başlangıcından 2.5 s sonraki kümülatif optik akışı özetlemektedir (koşul eşleştirilmiş Wilcoxon işaretli sıra testi: p = 0.031). Optik akış değerlerinin, temsili verilerde kaydedilen spontan davranışa kıyasla daha büyük olduğunu ve talimatsız, nüanslı davranışsal tepkileri değerlendirmenin zorluğunu daha da vurguladığını unutmayın. Şekil 5C, iki ödül koşulunun her biri için lickometer tarafından kaydedilen gerçek yalama oranını gösterirken, Şekil 5D, ödül teslimini takip eden 2,5 saniye içinde düşük ödül koşuluna kıyasla yüksek ödül koşulu sırasında daha fazla yalamanın kaydedildiğini göstermektedir (koşul eşleştirilmiş Wilcoxon işaretli sıra testi: p = 0.031). Toplamda, her iki analiz de yüksek ve düşük ödül tepkilerinin karşılaştırılmasında benzer bir eğilim ortaya çıkardı ve koşullar arasında hayvan davranışındaki anlamlı farklılıkları yakalamak için sunulan video analiz hattını doğruladı.
Birlikte ele alındığında, burada gösterilen temsili sonuçlar, yavru çağrılarının sunumunun, anne ve bakire dişi farelerde önemli ölçüde farklı tepkiler uyandırdığını, saf tonların ise uyandırmadığını göstermektedir. Bir doğrulama veri kümesinin analizi, optik akışta gözlemlenen farklılıkların, dış uyaranlara verilen davranışsal tepkilerdeki anlamlı farklılıkları yansıttığına dair güven verir.

Şekil 1: Video ortak kaydı. (A) Her biri farklı hayvanların farklı videolarından alınmış, üç kare halinde etiketlenmiş önceden belirlenmiş özellikleri (baş direğinin merkezi, sağ gözün önden görünümü, sağ gözün profil görünümü) gösteren video ortak kayıt örnekleri. Katmanda, üç noktanın nasıl hizalanmadığına dikkat edin, bu da davranışsal kurulumun veri toplama oturumları arasında nasıl biraz değiştiğini gösterir. (B) Video ortak kaydı gerçekleştirmeden önce (solda) ve sonra (sağda) her videonun son 10 saniyesindeki ortalama kare. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 2: Optik akış tahmini. (A) Mavi ile kaplanmış bağıl optik akış vektörlerine sahip örnek çerçeve. (B) Tüm videolarda optik akışın piksel bazında ortalama standart sapması. Ayna aracılığıyla hayvanın yüzünün profil görünümü etrafında seçilen yatırım getirisi macenta ile özetlenmiştir. Renk çubuğu standart sapmaya karşılık gelir. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 3: Gruplar ve koşullar arasında optik akışın görselleştirilmesi ve karşılaştırılması. (A) Bir annenin bir videosu ve bir bakirenin videosu için örnek Z-skorlu ve 5 Hz alçak geçiren filtreli optik akış zaman serisi. (B) ve (D) Yavru çağrılarının (B) ve saf tonların (D) uyaran periyodu sırasında ölçülen kümülatif optik akış. Gölgelendirme, ortalamanın (SEM) standart hatasını temsil eder. (C) ve (E) Yavru çağrıları (C) ve saf tonlar (E) için uyaran sunumunun ilk 2,5 saniyesi sırasında kümülatif optik akış. * P < 0.05, grup arası Mann-Whitney U testi (yavru çağrıları: p = 0.026; saf tonlar: p = 0.093), N = 6 grup başına gösterir. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 4: Gruplar arası PC1 karşılaştırması. (A,B) Gruplar arasında yavru çağrıları (A) ve saf tonlar (B) sırasında optik akışın ilk ana bileşeninin (PC1) yüklerini gösteren haritalar, çoklu karşılaştırmalar için düzeltme olmaksızın istatistiksel olarak p < 0.05'te eşiklenmiştir. Isı haritası, her pikselin optik akış dalgalanmasının diğer piksellere göre PC1'e ne kadar katkıda bulunduğunu gösterir. Küçük siyah metin, PC1 tarafından açıklanan varyans yüzdesini gösterir (annelerin yavru çağrıları: %4,95 ± %0,72; bakirelerin yavru çağrıları: %4,44 ± %0,59; annelerin saf tonları: %4,95 ± %0,39; bakirelerin saf tonları: %4,72 ± %1,24 (ortalama ± standart sapma)). Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 5: Davranışsal analiz hattının doğrulanması. (A) Yüksek ve düşük ödül dağıtımı sırasında ortalama Z-puanlı, uyaran öncesi taban çizgisi çıkarılmış ve 5 Hz düşük geçiren filtrelenmiş optik akış zaman serisi. Gölgelendirme SEM'i temsil eder, beyaz dikey çizgi ödül dağıtımını gösterir ve kırmızı yıldız işaretleri, ışık işaretinin başlangıcı ve ofseti nedeniyle geçici optik akış artefaktlarını gösterir. (B) Ödülün teslim edilmesinden sonraki ilk 2,5 saniye boyunca kümülatif optik akışın özeti. * P < 0.05'i, koşul eşleştirilmiş Wilcoxon işaretli sıra testi arasında gösterir (p = 0.031). (C) Yüksek ve düşük ödül dağıtımı sırasında lickometre ile elde edilen ortalama yalama oranı. Gölgelendirme SEM'i temsil eder ve beyaz dikey çizgi ödül dağıtımını gösterir. (D) Ödülün teslim edilmesinden sonraki ilk 2,5 saniye içinde kaydedilen yalama sayısı. * p < 0.05, koşul arası eşleştirilmiş Wilcoxon işaretli sıra testi (p = 0.031), N = 6'yı gösterir. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Video 1: Örnek video 1. PCR_Br011_20231015_1842_output.avi. Bu videoyu indirmek için lütfen buraya tıklayın.
Video 2: Örnek video 2. PCR_Br014_20231015_1722_output.avi. Bu videoyu indirmek için lütfen buraya tıklayın.
Ek Dosya 1: Video 1 için olay zamanlama dosyası. PCR_Br011_20231015_1842_output_videoTimestamps. Mat. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Dosya 2: Video 2 için olay zamanlama dosyası. PCR_Br014_20231015_1722_output_videoTimestamps. Mat. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Kodlama Dosyası 1: script1_videocoreg.m. Bu komut dosyası, ortak kayıt yoluyla tüm videoları uzamsal olarak birbirine hizalar. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Kodlama Dosyası 2: script2_optflow_roiselect.m. Bu komut dosyası, dönüştürülen her videonun tam FOV'u için optik akışı tahmin eder ve işlem hattının geri kalanı için bir ROI seçimine olanak tanır. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Kodlama Dosyası 3: script3_optflow_analysis.m. Bu komut dosyası, seçilen yatırım getirisi için farklı gruplar/koşullar arasında optik akış büyüklüğünü karşılaştırır. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Kodlama Dosyası 4: script4_optflow_pca.m. Bu komut dosyası, seçilen ROI'nin tahmini optik akışı üzerinde PCA gerçekleştirir. Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Kodlama Dosyası 5: script5_optflow_pca_analysis.m. Bu komut dosyası, seçilen ROI için PCA sonuçlarını gruplar/koşullar arasında karşılaştırır. Bu dosyayı indirmek için lütfen buraya tıklayın.