摘 要：目的 運用MRI技術觀察人類和獼猴顳上回負責聽覺信息加工區域的結構連接的相似和差異，探討兩個物種結構連接模式的同源性，為跨物種比較和替代性實驗提供理論依據。材料與方法 選擇公開的腦影像數據集，人類被試22例和獼猴被試9例作為研究對象。對於預處理後的人類數據提取6個感興趣區域與18個目標區域，獼猴數據提取1個感興趣區域與18個目標區域，計算連接強度並構建連接指紋，利用曼哈頓距離和餘弦相似度分析連接模式的相似性與差異性，最後用排列置換測試對結果進行檢驗。結果 人類顳上回(area 38m)和獼猴顳上回前側的連接指紋的相似度值最高(cosine=0.966)，並且area 38m右側的結構連接模式與獼猴顳上回前側的連接模式更為相似。兩個區域均與情景記憶相關的鼻周皮層，獎賞、情緒和運動相關的腹側紋狀體以及記憶和學習相關的海馬區有較強的連接，與認知和運動相關的頂葉島蓋連接較弱。結論 人類顳上回(area38m)具有和獼猴顳上回前側相似的結構連接模式，可能具有共同的進化起源，這一發現為構建人類和獼猴跨物種比較框架提供了理論支持。
關鍵詞：顳上回 同源 連接模式 磁共振成像 人類 獼猴
Homology study of auditory region of superior temporal gyrus between human and macaque based on MRI
CHAI Jingwen WANG Qianshan YAO Rong WANG Yue LI Binqiang LI Haifang
College of Information and Computer science, Taiyuan University of Technology;
Abstract：Objective: To explore the homology of structural connectivity patterns between the two species and to provide a theoretical basis for cross-species comparisons and alternative experiments, the similarities and differences in the structural connections responsible for auditory information processing in the superior temporal gyrus of humans and macaques were observed using magnetic resonance imaging. Materials and Methods: A publicly available brain imaging dataset with 22 human subjects and 9 macaque subjects was selected as the study subjects. For the pre-processed human data, 6 regions of interest with 18 target regions were extracted, and for the macaque data, the STGr and 18 target regions was extracted. Connection strengths were calculated and connection fingerprints were constructed. Similarities and differences in connection patterns were analysed using Manhattan distance and cosine similarity, and the results were finally examined using a permutation test. Results: The result shows that the highest similarity values(cosine=0.966) were found in the connectivity fingerprints of the human superior temporal gyrus(area 38 m) and the macaque anterior superior temporal gyrus, and the structural connectivity patterns in the right side of the human superior temporal gyrus(area 38 m) were more similar to those in the macaque anterior superior temporal gyrus. Both regions had strong connections to the perinasal cortex associated with situational memory, the ventral striatum associated with reward, emotion and movement, and the hippocampus associated with memory and learning, and weaker connections to the parietal insula associated with cognition and movement. Conclusions: Therefore, the human area 38 m has a similar structural connectivity profile to the rostral superior temporal gyrus of macaque. In other words, macaque and human audition processing areas might have a similar precursor. This discovery provides theoretical support for constructing the cross-species comparison framework between human and macaques.
Keyword：superior temporal gyrus; homology; connectivity; magnetic resonance imaging; human; macaque;
感覺系統的功能是對外界環境進行快速有效地感知並確定感覺信息的來源。靈長類動物顳上回(superior temporal gyrus,STG)是位於顳葉上部的一道溝壑，是處理聽覺訊息的中樞。其中人類的聽覺皮層位於顳葉皮層，獼猴聽覺皮層的範圍是從核心聽覺區到顳上回前側(rostral superior temporal gyrus,STGr)和顳極背側[1]。人類和獼猴利用聲音信號來參與社會互動，當施加聲音刺激後發現獼猴在顳上回區域對聲音表現出與人類相似的反應[2]，並且獼猴顳上回的belt區域關注認知功能[3]，這與人類的顳上回可能有對應關係[4-5]。顳上回一直被視為大腦中聽覺和言語的重要樞紐，受損時可能會導致聽覺性失語症。因此，研究顳上回區域及其包含的聽覺功能對聽覺模式的理解是有重要意義的。
雖然獼猴與人類的諸多特徵具有相似性，但在形態、功能等方面也存在著較大的進化差異[6-7]。因此，動物實驗結論在物種間推廣的可行性需要經過嚴謹的跨物種比較分析。連接性是對一個區域如何連接到其他大腦區域的描述[8]。目前已有文獻基於連接的聚類方法推斷人類與獼猴組織原則的相似性[9]；連接指紋概念提供了一個比較不同大腦皮層區域的連接性的新方法，用連接模式表征大腦皮層區域[10]。已有研究基於連接對人類和獼猴進行跨物種比較推測出人類和獼猴的同源區[11-12]，但是關於顳上回區域人類和獼猴的同源性仍有待研究。
本研究利用彌散張量成像(diffusion tensor imaging,DTI)數據，並輔以纖維束追蹤的處理技術，基於大腦感興趣區域(region of interest,ROI)與目標區域建立腦區之間的連接。利用連接指紋比較方法對人類和獼猴顳上回區域的結構連接模式進行深入研究，有助於從連接模式的角度探究兩個物種之間的同源性，為探索人類特有高級認知功能和獼猴動物模型相關研究的精準展開提供支持，為替代性實驗提供理論依據。
1 材料和方法
1.1 感興趣區域提取
人類腦影像數據來源於人腦連接組項目(Human Connectome Project)數據中心，原始數據可以從網站//tu.ta7.cn/365img/Login.vm;jsessionid=5945C8EBAE8568584F5F9A78154E6686 magnetic resonance imaging,s MRI)數據。其中DTI數據包括多梯度的1.25 mm各向同性數據，s MRI數據包括0.7 mm高解析度的各向同性T1加權成像數據和T2加權成像數據。具體採集參數為：T1W:TE=2.14 ms,TR=2400 ms,TI=1000 ms，翻轉角=7°，體素解析度為0.7 mm×0.7 mm×0.7 mm;DTI:TE=89.5 ms,TR=5520 ms，彌散加權b值1000 s/mm2,2000 s/mm2,3000 s/mm2。
獼猴腦影像數據源於加利福尼亞大學戴維斯分校神經科學中心UC Davis，原始數據可以從網站http://fcon＿1000.projects.nitrc.org/indi/PRIMEdownloads.html下載[13]，飼養的獼猴符合UC-Davis IACUC道德認證，且該數據的採集實驗得到了當地倫理委員會的批准，並遵守歐盟關於保護用於科學目的的動物的指令(2010/63/EU)。採取同樣的入組和排除標準，研究選取了9例被試，被試年齡分布在18.5～22.5歲，體質量分布在7.28～14.95 kg的獼猴。選擇的腦影像數據包括彌散張量成像數據和結構磁共振成像數據。具體採集參數為：T1W:TE=3.65 ms,TR=2500 ms,TI=1100 ms，翻轉角=7°，體素解析度為0.3 mm×0.3 mm×0.3 mm;DTI:TE=115 ms,TR=6400 ms，彌散加權b值1600 s/mm2,800 s/mm2。
DTI數據主要基於牛津大學的FSL工具(https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/)完成。DTI數據處理包括：頭動渦流矯正，梯度方向校正，獲取大腦掩膜圖像，將圖譜配准到個體以及使用FMRIB提供的雙張量模型BEDPOSTX[14]估計彌散參數等操作。
研究顯示，中國科學院自動化研究所的人類腦網絡組圖譜對腦區精細劃分，有客觀精準的邊界定位，是第一次建立宏觀尺度上的活體全腦連接圖譜[15]。於是，本研究基於中國科學院自動化研究所腦網絡組圖譜和D99獼猴腦圖譜[16]，使用FSL根據圖譜腦區對應編號基於個體提取人類和獼猴顳上回的感興趣區域與目標區域。其中，參考早期跨物種研究Neubert等人提供的同源位點坐標，結合所選腦圖譜配准到Neubert提供坐標的空間，根據坐標點對應出腦區編號並提取出來作為目標區域。所用圖譜確定了18個同源腦區作為目標區域。目標區域包括腦島、輔助運動區、鼻周皮層、海馬體、海馬旁回等聽覺相關腦區以及其他與軀體感覺、運動、執行等功能相關的腦區，如表1所示。本文所選擇人類感興趣區域見圖1A所示，分別為：STG＿6＿1 (Area 38m，紅色區域)，STG＿6＿2 (Area 41/42，青色區域)，STG＿6＿3 (Area TE，黃色區域)，STG＿6＿4 (Area 22c，藍色區域)，STG＿6＿5(Area 38l，綠色區域)，STG＿6＿6 (Area 22r，粉色區域)。選擇獼猴感興趣區域為：STGr (見圖1B紅色區域)。
1.2 數據處理
使用牛津大學FMRIB軟體包提供的PROBTRACKX從感興趣區域到目標腦區進行機率纖維束追蹤，通過設置感興趣區域為種子點，目標區域為必經節點和終止節點來指導機率纖維束走向，估計感興趣區域和目標腦區間的機率連接情況，將追蹤結果值作為腦區之間的連接值。
連接指紋表示ROI與已確定獼猴和人類之間同源性的目標大腦區域的連接性。使用基於個體空間提取的ROI和18個同源腦區，進行機率纖維束追蹤。為了避免個體差異導致的實驗誤差，研究利用所有被試與目標區域連接的平均值作為構建連接指紋的連接值。標準化連接值後，建立ROI與目標區域的連接指紋圖，確定獼猴ROI連接指紋和每個人類ROI連接指紋之間結構連接模式的相似與差異。
1.3 同源性分析
在計算相似性之前，所有的連接指紋都需要標準化。連接指紋比較的是連接強度，即ROI與目標區域的連接模式。
常用的樣本相似性度量方法包括歐氏距離、曼哈頓距離和餘弦相似度等。在基於連接指紋圖的跨物種比較過程中，通常選擇曼哈頓距離和餘弦相似度兩個指標進行度量。曼哈頓距離又稱城市街區距離，通過距離的測量決定連接指紋的差異，其定義為
式中，p和q是連通指紋的兩個向量，i是指紋的n個目標區域。d的值越小表示兩個向量的相似性越高。而餘弦相似度是通過計算兩個向量的夾角餘弦值來評估他們的相似性，計算過程中關注的是指紋圖的角度，忽略特徵向量的長度，其定義為
式中，向量之間的夾角越小，cosine的值越接近1，表示兩個向量的相似性越高。
1.4 統計檢驗
由於所使用的被試數量較少，兩種被試的年齡、性別的分布不太均勻，因此這些被試不能完全反應物種整體水平在腦區連接結構方面的一些差異。針對以上問題，採用python以及Num Py庫對數據進行統計學分析。對人類與獼猴構建的連接指紋採用置換檢驗，通過對樣本進行順序修改置換5000次，重新計算檢驗統計量，構造經驗分布，對其可信度進行檢驗，並重複檢驗1000次，得到一組P值，使用Matlab基於linear step-up方法進行校正，在此基礎上求出P-value進行推斷。以P<0.05為差異有統計學意義。
2 結果
本文對人類和獼猴顳上回的感興趣區域和已知的人與獼猴同源區域建立連接指紋，通過計算曼哈頓距離和餘弦相似度，得到結構連接模式的相似性，如表2所示。
人類顳上回(area 38m)和獼猴STGr的連接指紋如圖2所示。人類area 38m和獼猴STGr連接指紋的曼哈頓距離結果如圖3所示，其中人類area 38m與獼猴感興趣區域的連接指紋距離值(d=0.733)是最小的，其對應的餘弦相似度的值是最大的(cosine=0.966)。
當考慮左右腦的差異，人腦左、右半球area 38m腦區連接指紋和排列置換測試結果如圖4所示，置換檢驗P<0.05，其中置換檢驗結果紅色代表真實值，藍色代表閾值，灰色代表測試的數據。area 38m左、右半球連接指紋對應的曼哈頓距離值分別為0.983和0.554，餘弦相似度的值分別為0.948和0.975。
3 討論
通過分析磁共振成像數據，本文對兩個物種顳上回聽覺區的結構連接模式進行同源性研究，發現人類顳上回(area 38m)區域與獼猴顳上回前側的結構連接相似度值最高，因此人類顳上回(area 38m)具有和獼猴顳上回前側相似的結構連接模式，可能具有共同的進化起源。
3.1 跨物種同源性分析
現有研究通過對兩個物種施加聽覺刺激，發現了兩個物種的顳葉前側均擁有語音區，並指出獼猴和人類的高級聽覺皮層的功能組織比目前已知的更為相似[17]。另一研究使用相同刺激對比獼猴和人類的反應得到猴腦的敏感性低於人腦，其原因可能是由於人腦功能的專門化[18]。本文在研究人類與獼猴顳上回參與聽覺處理區域的過程中，關注的是大腦解剖連接模式，而不是執行任務時大腦區域的功能激活。本文在人類的顳上回，搜索與獼猴顳上回前側具有相同結構連接模式的區域。實驗中用於識別人類和獼猴顳上回區域連接模式相似性的目標區域是已知同源的。分析發現人類的area 38m與獼猴顳上回前側具有最相似的結構連接模式。
利用結構連接指紋分析兩個物種的不同區域，可以解釋這些區域在物種進化過程中腦連接的差異。2018年劉義鵬等人就腦纖維可視化總結得到量化腦結構連接可以揭示在腦發育過程中的生理變化，有助於研究神經病理條件下腦網絡改變的結果。本研究結果顯示，在人類顳上回的6個感興趣區域中，area 38m與獼猴的連接指紋曼哈頓距離值最小(d=0.733)，對應餘弦相似度的值最大(cosine=0.966)。對兩個物種的連接模式分析得到：area 38m和獼猴STGr對鼻周皮層、運動皮層等區域的連接強度是一致的，其差異在於獼猴STGr與44v區域的連接比人類的更強，而人類area 38m與腹側紋狀體、海馬區的連接比獼猴更強。
雖然大腦左右半球在形態上大體對稱，但是在結構連接和功能上卻不完全對稱，即大腦偏側化。顯然，研究中顳上回area 38m＿L區域和area 38m＿R區域與獼猴顳上回前側都有著相似的結構連接模式：兩個區域均與情緒記憶相關的鼻周皮層[19]，獎賞、情緒和運動相關的腹側紋狀體以及記憶和學習相關的海馬區有較強的連接，與認知和運動相關的頂葉島蓋連接較弱。但也存在部分差異，兩個區域分別完成不同任務的處理：人類顳上回area 38m＿L區域相比顳上回area 38m＿R區域，顯示出與海馬旁回更強的連接。海馬旁回位於枕葉和顳葉下方的內側，與情感、學習及記憶等高級神經功能活動有關。這也驗證了先前對人類area 38m＿L區域的描述，它是更關注情感信息的區域。以上結果驗證了顳葉皮質在人腦中的擴張是優先於獼猴大腦的，得到人腦具有不同功能特徵的分區這一結論。
根據人類和獼猴連接指紋統計量的排列置換測試結果，證明人類area 38m與獼猴STGr區域有相似的結構連接模式這一結果是可信的。分別為左、右半球構建連接指紋後，曼哈頓距離和餘弦相似度結果均顯示人類和獼猴腦區之間的結構連接模式有較高相似性，但是左、右半球對應連接指紋的排列置換測試的結果卻不一致。area 38m＿R檢驗統計量的值是落在置信區間內的，證明該結果可信。而area38m＿L的檢驗統計量的值均位於置信區間外，這表示小機率事件發生，樣本所反應結果的可信度較低。觀察area 38m＿L區域所對應的閾值區間，發現該區域的餘弦相似度值較高(cosine=0.948)，且位於左閾值附近，即小機率事件呈現出了較低的跨物種一致性。
導致獼猴顳上回前側與人類area 38m＿R區域連接更為相似的原因可能是大腦的偏側化。在2018年梅丹丹等人根據先驗知識探究發現，人腦加工任務的半球偏側化導致左半球對言語功能、語言性聲音的處理、對聽覺語音語義信息加工更有優勢，其行為領域主要是情緒的加工；而右半球對不需要言語參加的空間知覺和形象思維、環境聲音與音樂的處理更有優勢。人類腦區結構和功能的變化造成與其他腦區連接關係的增強或減弱，因此導致了左右腦對於目標腦區連接的差異[20]。
3.2 研究意義
聽覺作為主要的感知方式之一，它是個體感知外界環境、維持生存和繁殖活動的重要途徑。大腦皮質顳上回主要是聽覺性語言中樞。顳上回受損傷可能出現聽覺性失語症或稱接受性失語症等併發症，也有引起癲癇發作的可能。因此，本文利用DTI數據，通過比較影像學量化物種間皮層區域在連接上的相似與差異，分析出人類與獼猴在顳上回聽覺處理相關區域的同源性，可以幫助了解活體的神經纖維走向，有助於深入了解人腦結構，為神經科學研究提供了理論支撐。由於許多精神病理綜合徵與大腦區域和連接之間的改變有關，構建連接指紋也可以作為精神病理學的診斷工具。對於有興趣研究聽覺處理的研究者而言，可以對顳上回前側區域進一步探究。
3.3 限制與不足
在研究人類和獼猴顳上回聽覺區域的同源性中，發現的兩個物種結構連接模式的相似性，能夠為探索人類特有高級認知功能和獼猴動物模型相關研究提供支持，也為替代性實驗提供了理論依據。但本文對兩個物種顳上回同源性的研究主要依據是大腦結構連接模式，其功能的連接差異需要進一步討論。目前研究發現人與猴子均優先使用左半球聽覺系統分析聲音，且存在左半球比右半球變化更強的功能不對稱現象[21]。現有研究發現人類顳上回在噪音中對語音感知有著關鍵作用，並且對低頻純音刺激的響應區域面積較大。因此，獼猴的左右半球是否存在與人類同樣的連接變化趨勢也有待研究。其次，在本研究中基於連接指紋對顳上回連接模式進行跨物種比較的方法的比較特徵精度和可擴展性較低，每個待比較腦區的連接指紋應有獨特的目標腦區來構建，後續研究可以探索更完善的高精度跨物種比較方法，進一步驗證和精細化現有結果。
綜上所述，本研究發現人類area 38m與獼猴顳上回前側STGr是同源的，即具有相似的連接模式。其中，人類右半球與獼猴的結構連接模式相似性比左半球更高，這些現象可能與大腦的偏側化相關。
作者利益衝突聲明：全體作者均聲明無利益衝突。
參考文獻
[1] Scott BH, Mishkin M, Yin P. Neural correlates of auditory short-term memory in rostral superior temporal cortex[J]. Curr Biol, 2014, 24(23):2767-2775. DOI:10.1016/j.cub.2014.10.004.
[2] Joly O, Pallier C, Ramus F, et al. Processing of vocalizations in humans and monkeys:A comparative fMRI study[J]. NeuroImage,2012, 62(3):1376-1389. DOI:10.1016/j.neuroimage.2012.05.070.
[3] Rauschecker JP. Where did language come from? Precursor mechanisms in nonhuman primates[J]. Curr Opin Behav Sci, 2018,21(1):195-204. DOI:10.1016/j.cobeha.2018.06.003.
[4] Yamamoto AK, Parker JO, Hope TMH, et al. A special role for the right posterior superior temporal sulcus during speech production[J].Neuro Image, 2019, 203:116184. DOI:10.1016/j.neuroimage.2019.116184.
[5] Nakai T, Koide-Majima N, Nishimoto S. Correspondence of categorical and feature-based representations of music in the human brain[J]. Brain Behav, 2021, 11(1):1936. DOI:10.1002/brb3.1936.
[6] Passingham R. How good is the macaque monkey model of the human brain?[J]. Curr Opin Neurobiol, 2009, 19(1):6-11. DOI:10.1016/j.conb.2009.01.002.
[7] Rilling JK. Comparative primate neurobiology and the evolution of brain language systems[J]. Curr Opin Neurobiol, 2014, 28(5):10-14.DOI:10.1016/j.conb.2014.04.002.
[8] Passingham RE, Stephan KE, Kötter R. The anatomical basis of functional localization in the cortex[J]. Nat Rev Neurosci, 2002, 3(8):606-616. DOI:10.1038/nrn893.
[9] Beckmann M, Johansen-Berg H, Rushworth MF. Connectivity-based parcellation of human cingulate cortex and its relation to functional specialization[J]. J Neurosci, 2009, 29(4):1175-1190. DOI:10.1523/JNEUROSCI.3328-08.2009.
[10] Mars RB, Verhagen L, Gladwin TE, et al. Comparing brains by matching connectivity profiles[J]. Neurosci Biobehav Rev, 2016,60(5):90-97. DOI:10.1016/j.neubiorev.2015.10.008.
[11] Xia XL, Fan LZ, Cheng C, et al. Interspecies Differences in the Connectivity of Ventral Striatal Components Between Humans and Macaques[J]. Front Neurosci, 2019, 13:623. DOI:10.1016/j.neubiorev.2015.10.008.
[12] Wang JJ, Becker B, Wang LJ, et al. Corresponding anatomical and coactivation architecture of the human precuneus showing similar connectivity patterns with macaques[J]. NeuroImage, 2019, 200:562-574. DOI:10.1016/j.neuroimage.2019.07.001.
[13] Milham MP, Ai L, Koo B, et al. An Open Resource for Non-human Primate Imaging[J]. Neuron, 2018, 100(1):61-74. DOI:10.1016/j.neuron.2018.08.039.
[14] Croxson PL, Forkel SJ, Cerliani L, et al. Structural Variability Across the Primate Brain:A Cross-Species Comparison[J]. Cereb Cortex,2018, 28(11):3829-3841. DOI:10.1093/cercor/bhx244.
[15] Fan LZ, Li H, Zhuo JJ, et al. The Human Brainnetome Atlas:A New Brain Atlas Based on Connectional Architecture[J]. Cereb Cortex,2016, 26(8):3508-3526. DOI:10.1093/cercor/bhw157.
[16] Reveley C, Gruslys A, Ye FQ, et al. Three-Dimensional Digital Template Atlas of the Macaque Brain[J]. Cereb Cortex, 2017, 27(9):4463-4477. DOI:10.1093/cercor/bhw248.
[17] Bodin C, Trapeau R, Nazarian B, et al. Functionally homologous representation of vocalizations in the auditory cortex of humans and macaques. Curr Biol, 2021, 31(21):4839-4844. DOI:10.1016/j.cub.2021.08.043.
[18] Dheerendra P, Baumann S, Joly O, et al. The Representation of Time Windows in Primate Auditory Cortex[J]. Cereb Cortex, 2021, 7:bhab434. DOI:10.1093/cercor/bhab434.
[19] Li XS, Li XH, Chen SJ, et al. Effect of emotional enhancement of memory on recollection process in young adults:the influence factors and neural mechanisms[J]. Brain Imaging Behav, 2020, 14(1):119-129. DOI:10.1007/s11682-018-9975-0.
[20] Xia XL, Gao F, Yuan Z. Species and individual differences and connectional asymmetry of Broca's area in humans and macaques[J].Neuroimage, 2021, 244:118583. DOI:10.1016/j.neuroimage.2021.118583.
[21] Ren JX, Xu T, Wang DH, et al. Individual Variability in Functional Organization of the Human and Monkey Auditory Cortex[J]. Cereb Cortex, 2021, 31(5):2450-2465. DOI:10.1093/cercor/bhaa366.
關鍵詞：顳上回 同源 連接模式 磁共振成像 人類 獼猴
Homology study of auditory region of superior temporal gyrus between human and macaque based on MRI
CHAI Jingwen WANG Qianshan YAO Rong WANG Yue LI Binqiang LI Haifang
College of Information and Computer science, Taiyuan University of Technology;
Abstract：Objective: To explore the homology of structural connectivity patterns between the two species and to provide a theoretical basis for cross-species comparisons and alternative experiments, the similarities and differences in the structural connections responsible for auditory information processing in the superior temporal gyrus of humans and macaques were observed using magnetic resonance imaging. Materials and Methods: A publicly available brain imaging dataset with 22 human subjects and 9 macaque subjects was selected as the study subjects. For the pre-processed human data, 6 regions of interest with 18 target regions were extracted, and for the macaque data, the STGr and 18 target regions was extracted. Connection strengths were calculated and connection fingerprints were constructed. Similarities and differences in connection patterns were analysed using Manhattan distance and cosine similarity, and the results were finally examined using a permutation test. Results: The result shows that the highest similarity values(cosine=0.966) were found in the connectivity fingerprints of the human superior temporal gyrus(area 38 m) and the macaque anterior superior temporal gyrus, and the structural connectivity patterns in the right side of the human superior temporal gyrus(area 38 m) were more similar to those in the macaque anterior superior temporal gyrus. Both regions had strong connections to the perinasal cortex associated with situational memory, the ventral striatum associated with reward, emotion and movement, and the hippocampus associated with memory and learning, and weaker connections to the parietal insula associated with cognition and movement. Conclusions: Therefore, the human area 38 m has a similar structural connectivity profile to the rostral superior temporal gyrus of macaque. In other words, macaque and human audition processing areas might have a similar precursor. This discovery provides theoretical support for constructing the cross-species comparison framework between human and macaques.
Keyword：superior temporal gyrus; homology; connectivity; magnetic resonance imaging; human; macaque;
感覺系統的功能是對外界環境進行快速有效地感知並確定感覺信息的來源。靈長類動物顳上回(superior temporal gyrus,STG)是位於顳葉上部的一道溝壑，是處理聽覺訊息的中樞。其中人類的聽覺皮層位於顳葉皮層，獼猴聽覺皮層的範圍是從核心聽覺區到顳上回前側(rostral superior temporal gyrus,STGr)和顳極背側[1]。人類和獼猴利用聲音信號來參與社會互動，當施加聲音刺激後發現獼猴在顳上回區域對聲音表現出與人類相似的反應[2]，並且獼猴顳上回的belt區域關注認知功能[3]，這與人類的顳上回可能有對應關係[4-5]。顳上回一直被視為大腦中聽覺和言語的重要樞紐，受損時可能會導致聽覺性失語症。因此，研究顳上回區域及其包含的聽覺功能對聽覺模式的理解是有重要意義的。
雖然獼猴與人類的諸多特徵具有相似性，但在形態、功能等方面也存在著較大的進化差異[6-7]。因此，動物實驗結論在物種間推廣的可行性需要經過嚴謹的跨物種比較分析。連接性是對一個區域如何連接到其他大腦區域的描述[8]。目前已有文獻基於連接的聚類方法推斷人類與獼猴組織原則的相似性[9]；連接指紋概念提供了一個比較不同大腦皮層區域的連接性的新方法，用連接模式表征大腦皮層區域[10]。已有研究基於連接對人類和獼猴進行跨物種比較推測出人類和獼猴的同源區[11-12]，但是關於顳上回區域人類和獼猴的同源性仍有待研究。
本研究利用彌散張量成像(diffusion tensor imaging,DTI)數據，並輔以纖維束追蹤的處理技術，基於大腦感興趣區域(region of interest,ROI)與目標區域建立腦區之間的連接。利用連接指紋比較方法對人類和獼猴顳上回區域的結構連接模式進行深入研究，有助於從連接模式的角度探究兩個物種之間的同源性，為探索人類特有高級認知功能和獼猴動物模型相關研究的精準展開提供支持，為替代性實驗提供理論依據。
1 材料和方法
1.1 感興趣區域提取
人類腦影像數據來源於人腦連接組項目(Human Connectome Project)數據中心，原始數據可以從網站//tu.ta7.cn/365img/Login.vm;jsessionid=5945C8EBAE8568584F5F9A78154E6686 magnetic resonance imaging,s MRI)數據。其中DTI數據包括多梯度的1.25 mm各向同性數據，s MRI數據包括0.7 mm高解析度的各向同性T1加權成像數據和T2加權成像數據。具體採集參數為：T1W:TE=2.14 ms,TR=2400 ms,TI=1000 ms，翻轉角=7°，體素解析度為0.7 mm×0.7 mm×0.7 mm;DTI:TE=89.5 ms,TR=5520 ms，彌散加權b值1000 s/mm2,2000 s/mm2,3000 s/mm2。
獼猴腦影像數據源於加利福尼亞大學戴維斯分校神經科學中心UC Davis，原始數據可以從網站http://fcon＿1000.projects.nitrc.org/indi/PRIMEdownloads.html下載[13]，飼養的獼猴符合UC-Davis IACUC道德認證，且該數據的採集實驗得到了當地倫理委員會的批准，並遵守歐盟關於保護用於科學目的的動物的指令(2010/63/EU)。採取同樣的入組和排除標準，研究選取了9例被試，被試年齡分布在18.5～22.5歲，體質量分布在7.28～14.95 kg的獼猴。選擇的腦影像數據包括彌散張量成像數據和結構磁共振成像數據。具體採集參數為：T1W:TE=3.65 ms,TR=2500 ms,TI=1100 ms，翻轉角=7°，體素解析度為0.3 mm×0.3 mm×0.3 mm;DTI:TE=115 ms,TR=6400 ms，彌散加權b值1600 s/mm2,800 s/mm2。
DTI數據主要基於牛津大學的FSL工具(https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/)完成。DTI數據處理包括：頭動渦流矯正，梯度方向校正，獲取大腦掩膜圖像，將圖譜配准到個體以及使用FMRIB提供的雙張量模型BEDPOSTX[14]估計彌散參數等操作。
研究顯示，中國科學院自動化研究所的人類腦網絡組圖譜對腦區精細劃分，有客觀精準的邊界定位，是第一次建立宏觀尺度上的活體全腦連接圖譜[15]。於是，本研究基於中國科學院自動化研究所腦網絡組圖譜和D99獼猴腦圖譜[16]，使用FSL根據圖譜腦區對應編號基於個體提取人類和獼猴顳上回的感興趣區域與目標區域。其中，參考早期跨物種研究Neubert等人提供的同源位點坐標，結合所選腦圖譜配准到Neubert提供坐標的空間，根據坐標點對應出腦區編號並提取出來作為目標區域。所用圖譜確定了18個同源腦區作為目標區域。目標區域包括腦島、輔助運動區、鼻周皮層、海馬體、海馬旁回等聽覺相關腦區以及其他與軀體感覺、運動、執行等功能相關的腦區，如表1所示。本文所選擇人類感興趣區域見圖1A所示，分別為：STG＿6＿1 (Area 38m，紅色區域)，STG＿6＿2 (Area 41/42，青色區域)，STG＿6＿3 (Area TE，黃色區域)，STG＿6＿4 (Area 22c，藍色區域)，STG＿6＿5(Area 38l，綠色區域)，STG＿6＿6 (Area 22r，粉色區域)。選擇獼猴感興趣區域為：STGr (見圖1B紅色區域)。
1.2 數據處理
使用牛津大學FMRIB軟體包提供的PROBTRACKX從感興趣區域到目標腦區進行機率纖維束追蹤，通過設置感興趣區域為種子點，目標區域為必經節點和終止節點來指導機率纖維束走向，估計感興趣區域和目標腦區間的機率連接情況，將追蹤結果值作為腦區之間的連接值。
連接指紋表示ROI與已確定獼猴和人類之間同源性的目標大腦區域的連接性。使用基於個體空間提取的ROI和18個同源腦區，進行機率纖維束追蹤。為了避免個體差異導致的實驗誤差，研究利用所有被試與目標區域連接的平均值作為構建連接指紋的連接值。標準化連接值後，建立ROI與目標區域的連接指紋圖，確定獼猴ROI連接指紋和每個人類ROI連接指紋之間結構連接模式的相似與差異。
1.3 同源性分析
在計算相似性之前，所有的連接指紋都需要標準化。連接指紋比較的是連接強度，即ROI與目標區域的連接模式。
常用的樣本相似性度量方法包括歐氏距離、曼哈頓距離和餘弦相似度等。在基於連接指紋圖的跨物種比較過程中，通常選擇曼哈頓距離和餘弦相似度兩個指標進行度量。曼哈頓距離又稱城市街區距離，通過距離的測量決定連接指紋的差異，其定義為
式中，p和q是連通指紋的兩個向量，i是指紋的n個目標區域。d的值越小表示兩個向量的相似性越高。而餘弦相似度是通過計算兩個向量的夾角餘弦值來評估他們的相似性，計算過程中關注的是指紋圖的角度，忽略特徵向量的長度，其定義為
式中，向量之間的夾角越小，cosine的值越接近1，表示兩個向量的相似性越高。
1.4 統計檢驗
由於所使用的被試數量較少，兩種被試的年齡、性別的分布不太均勻，因此這些被試不能完全反應物種整體水平在腦區連接結構方面的一些差異。針對以上問題，採用python以及Num Py庫對數據進行統計學分析。對人類與獼猴構建的連接指紋採用置換檢驗，通過對樣本進行順序修改置換5000次，重新計算檢驗統計量，構造經驗分布，對其可信度進行檢驗，並重複檢驗1000次，得到一組P值，使用Matlab基於linear step-up方法進行校正，在此基礎上求出P-value進行推斷。以P<0.05為差異有統計學意義。
2 結果
本文對人類和獼猴顳上回的感興趣區域和已知的人與獼猴同源區域建立連接指紋，通過計算曼哈頓距離和餘弦相似度，得到結構連接模式的相似性，如表2所示。
人類顳上回(area 38m)和獼猴STGr的連接指紋如圖2所示。人類area 38m和獼猴STGr連接指紋的曼哈頓距離結果如圖3所示，其中人類area 38m與獼猴感興趣區域的連接指紋距離值(d=0.733)是最小的，其對應的餘弦相似度的值是最大的(cosine=0.966)。
當考慮左右腦的差異，人腦左、右半球area 38m腦區連接指紋和排列置換測試結果如圖4所示，置換檢驗P<0.05，其中置換檢驗結果紅色代表真實值，藍色代表閾值，灰色代表測試的數據。area 38m左、右半球連接指紋對應的曼哈頓距離值分別為0.983和0.554，餘弦相似度的值分別為0.948和0.975。
3 討論
通過分析磁共振成像數據，本文對兩個物種顳上回聽覺區的結構連接模式進行同源性研究，發現人類顳上回(area 38m)區域與獼猴顳上回前側的結構連接相似度值最高，因此人類顳上回(area 38m)具有和獼猴顳上回前側相似的結構連接模式，可能具有共同的進化起源。
3.1 跨物種同源性分析
現有研究通過對兩個物種施加聽覺刺激，發現了兩個物種的顳葉前側均擁有語音區，並指出獼猴和人類的高級聽覺皮層的功能組織比目前已知的更為相似[17]。另一研究使用相同刺激對比獼猴和人類的反應得到猴腦的敏感性低於人腦，其原因可能是由於人腦功能的專門化[18]。本文在研究人類與獼猴顳上回參與聽覺處理區域的過程中，關注的是大腦解剖連接模式，而不是執行任務時大腦區域的功能激活。本文在人類的顳上回，搜索與獼猴顳上回前側具有相同結構連接模式的區域。實驗中用於識別人類和獼猴顳上回區域連接模式相似性的目標區域是已知同源的。分析發現人類的area 38m與獼猴顳上回前側具有最相似的結構連接模式。
利用結構連接指紋分析兩個物種的不同區域，可以解釋這些區域在物種進化過程中腦連接的差異。2018年劉義鵬等人就腦纖維可視化總結得到量化腦結構連接可以揭示在腦發育過程中的生理變化，有助於研究神經病理條件下腦網絡改變的結果。本研究結果顯示，在人類顳上回的6個感興趣區域中，area 38m與獼猴的連接指紋曼哈頓距離值最小(d=0.733)，對應餘弦相似度的值最大(cosine=0.966)。對兩個物種的連接模式分析得到：area 38m和獼猴STGr對鼻周皮層、運動皮層等區域的連接強度是一致的，其差異在於獼猴STGr與44v區域的連接比人類的更強，而人類area 38m與腹側紋狀體、海馬區的連接比獼猴更強。
雖然大腦左右半球在形態上大體對稱，但是在結構連接和功能上卻不完全對稱，即大腦偏側化。顯然，研究中顳上回area 38m＿L區域和area 38m＿R區域與獼猴顳上回前側都有著相似的結構連接模式：兩個區域均與情緒記憶相關的鼻周皮層[19]，獎賞、情緒和運動相關的腹側紋狀體以及記憶和學習相關的海馬區有較強的連接，與認知和運動相關的頂葉島蓋連接較弱。但也存在部分差異，兩個區域分別完成不同任務的處理：人類顳上回area 38m＿L區域相比顳上回area 38m＿R區域，顯示出與海馬旁回更強的連接。海馬旁回位於枕葉和顳葉下方的內側，與情感、學習及記憶等高級神經功能活動有關。這也驗證了先前對人類area 38m＿L區域的描述，它是更關注情感信息的區域。以上結果驗證了顳葉皮質在人腦中的擴張是優先於獼猴大腦的，得到人腦具有不同功能特徵的分區這一結論。
根據人類和獼猴連接指紋統計量的排列置換測試結果，證明人類area 38m與獼猴STGr區域有相似的結構連接模式這一結果是可信的。分別為左、右半球構建連接指紋後，曼哈頓距離和餘弦相似度結果均顯示人類和獼猴腦區之間的結構連接模式有較高相似性，但是左、右半球對應連接指紋的排列置換測試的結果卻不一致。area 38m＿R檢驗統計量的值是落在置信區間內的，證明該結果可信。而area38m＿L的檢驗統計量的值均位於置信區間外，這表示小機率事件發生，樣本所反應結果的可信度較低。觀察area 38m＿L區域所對應的閾值區間，發現該區域的餘弦相似度值較高(cosine=0.948)，且位於左閾值附近，即小機率事件呈現出了較低的跨物種一致性。
導致獼猴顳上回前側與人類area 38m＿R區域連接更為相似的原因可能是大腦的偏側化。在2018年梅丹丹等人根據先驗知識探究發現，人腦加工任務的半球偏側化導致左半球對言語功能、語言性聲音的處理、對聽覺語音語義信息加工更有優勢，其行為領域主要是情緒的加工；而右半球對不需要言語參加的空間知覺和形象思維、環境聲音與音樂的處理更有優勢。人類腦區結構和功能的變化造成與其他腦區連接關係的增強或減弱，因此導致了左右腦對於目標腦區連接的差異[20]。
3.2 研究意義
聽覺作為主要的感知方式之一，它是個體感知外界環境、維持生存和繁殖活動的重要途徑。大腦皮質顳上回主要是聽覺性語言中樞。顳上回受損傷可能出現聽覺性失語症或稱接受性失語症等併發症，也有引起癲癇發作的可能。因此，本文利用DTI數據，通過比較影像學量化物種間皮層區域在連接上的相似與差異，分析出人類與獼猴在顳上回聽覺處理相關區域的同源性，可以幫助了解活體的神經纖維走向，有助於深入了解人腦結構，為神經科學研究提供了理論支撐。由於許多精神病理綜合徵與大腦區域和連接之間的改變有關，構建連接指紋也可以作為精神病理學的診斷工具。對於有興趣研究聽覺處理的研究者而言，可以對顳上回前側區域進一步探究。
3.3 限制與不足
在研究人類和獼猴顳上回聽覺區域的同源性中，發現的兩個物種結構連接模式的相似性，能夠為探索人類特有高級認知功能和獼猴動物模型相關研究提供支持，也為替代性實驗提供了理論依據。但本文對兩個物種顳上回同源性的研究主要依據是大腦結構連接模式，其功能的連接差異需要進一步討論。目前研究發現人與猴子均優先使用左半球聽覺系統分析聲音，且存在左半球比右半球變化更強的功能不對稱現象[21]。現有研究發現人類顳上回在噪音中對語音感知有著關鍵作用，並且對低頻純音刺激的響應區域面積較大。因此，獼猴的左右半球是否存在與人類同樣的連接變化趨勢也有待研究。其次，在本研究中基於連接指紋對顳上回連接模式進行跨物種比較的方法的比較特徵精度和可擴展性較低，每個待比較腦區的連接指紋應有獨特的目標腦區來構建，後續研究可以探索更完善的高精度跨物種比較方法，進一步驗證和精細化現有結果。
綜上所述，本研究發現人類area 38m與獼猴顳上回前側STGr是同源的，即具有相似的連接模式。其中，人類右半球與獼猴的結構連接模式相似性比左半球更高，這些現象可能與大腦的偏側化相關。
作者利益衝突聲明：全體作者均聲明無利益衝突。
參考文獻
[1] Scott BH, Mishkin M, Yin P. Neural correlates of auditory short-term memory in rostral superior temporal cortex[J]. Curr Biol, 2014, 24(23):2767-2775. DOI:10.1016/j.cub.2014.10.004.
[2] Joly O, Pallier C, Ramus F, et al. Processing of vocalizations in humans and monkeys:A comparative fMRI study[J]. NeuroImage,2012, 62(3):1376-1389. DOI:10.1016/j.neuroimage.2012.05.070.
[3] Rauschecker JP. Where did language come from? Precursor mechanisms in nonhuman primates[J]. Curr Opin Behav Sci, 2018,21(1):195-204. DOI:10.1016/j.cobeha.2018.06.003.
[4] Yamamoto AK, Parker JO, Hope TMH, et al. A special role for the right posterior superior temporal sulcus during speech production[J].Neuro Image, 2019, 203:116184. DOI:10.1016/j.neuroimage.2019.116184.
[5] Nakai T, Koide-Majima N, Nishimoto S. Correspondence of categorical and feature-based representations of music in the human brain[J]. Brain Behav, 2021, 11(1):1936. DOI:10.1002/brb3.1936.
[6] Passingham R. How good is the macaque monkey model of the human brain?[J]. Curr Opin Neurobiol, 2009, 19(1):6-11. DOI:10.1016/j.conb.2009.01.002.
[7] Rilling JK. Comparative primate neurobiology and the evolution of brain language systems[J]. Curr Opin Neurobiol, 2014, 28(5):10-14.DOI:10.1016/j.conb.2014.04.002.
[8] Passingham RE, Stephan KE, Kötter R. The anatomical basis of functional localization in the cortex[J]. Nat Rev Neurosci, 2002, 3(8):606-616. DOI:10.1038/nrn893.
[9] Beckmann M, Johansen-Berg H, Rushworth MF. Connectivity-based parcellation of human cingulate cortex and its relation to functional specialization[J]. J Neurosci, 2009, 29(4):1175-1190. DOI:10.1523/JNEUROSCI.3328-08.2009.
[10] Mars RB, Verhagen L, Gladwin TE, et al. Comparing brains by matching connectivity profiles[J]. Neurosci Biobehav Rev, 2016,60(5):90-97. DOI:10.1016/j.neubiorev.2015.10.008.
[11] Xia XL, Fan LZ, Cheng C, et al. Interspecies Differences in the Connectivity of Ventral Striatal Components Between Humans and Macaques[J]. Front Neurosci, 2019, 13:623. DOI:10.1016/j.neubiorev.2015.10.008.
[12] Wang JJ, Becker B, Wang LJ, et al. Corresponding anatomical and coactivation architecture of the human precuneus showing similar connectivity patterns with macaques[J]. NeuroImage, 2019, 200:562-574. DOI:10.1016/j.neuroimage.2019.07.001.
[13] Milham MP, Ai L, Koo B, et al. An Open Resource for Non-human Primate Imaging[J]. Neuron, 2018, 100(1):61-74. DOI:10.1016/j.neuron.2018.08.039.
[14] Croxson PL, Forkel SJ, Cerliani L, et al. Structural Variability Across the Primate Brain:A Cross-Species Comparison[J]. Cereb Cortex,2018, 28(11):3829-3841. DOI:10.1093/cercor/bhx244.
[15] Fan LZ, Li H, Zhuo JJ, et al. The Human Brainnetome Atlas:A New Brain Atlas Based on Connectional Architecture[J]. Cereb Cortex,2016, 26(8):3508-3526. DOI:10.1093/cercor/bhw157.
[16] Reveley C, Gruslys A, Ye FQ, et al. Three-Dimensional Digital Template Atlas of the Macaque Brain[J]. Cereb Cortex, 2017, 27(9):4463-4477. DOI:10.1093/cercor/bhw248.
[17] Bodin C, Trapeau R, Nazarian B, et al. Functionally homologous representation of vocalizations in the auditory cortex of humans and macaques. Curr Biol, 2021, 31(21):4839-4844. DOI:10.1016/j.cub.2021.08.043.
[18] Dheerendra P, Baumann S, Joly O, et al. The Representation of Time Windows in Primate Auditory Cortex[J]. Cereb Cortex, 2021, 7:bhab434. DOI:10.1093/cercor/bhab434.
[19] Li XS, Li XH, Chen SJ, et al. Effect of emotional enhancement of memory on recollection process in young adults:the influence factors and neural mechanisms[J]. Brain Imaging Behav, 2020, 14(1):119-129. DOI:10.1007/s11682-018-9975-0.
[20] Xia XL, Gao F, Yuan Z. Species and individual differences and connectional asymmetry of Broca's area in humans and macaques[J].Neuroimage, 2021, 244:118583. DOI:10.1016/j.neuroimage.2021.118583.
[21] Ren JX, Xu T, Wang DH, et al. Individual Variability in Functional Organization of the Human and Monkey Auditory Cortex[J]. Cereb Cortex, 2021, 31(5):2450-2465. DOI:10.1093/cercor/bhaa366.
收藏