ヨーロッパ中央部新石器時代最初期における農耕民と狩猟採集民との関係
取り上げるのが遅れてしまいましたが、ヨーロッパ中央部新石器時代最初期における農耕民と狩猟採集民との関係に関する研究(Nikitin et al., 2019)が公表されました。線形陶器文化(Linear Pottery、Linearbandkeramik、略してLBK)は、ヨーロッパ中央部の新石器時代の始まりにおいて重要な役割を果たしました。文化的・経済的・遺伝的に、LBKは究極的にはアナトリア半島西部に起源がありますが、在来のヨーロッパ中石器時代狩猟採集民社会の明確な特徴も示します。LBKの起源に関しては、いくつかのモデルが長年にわたって提案されてきました。在来モデルでは、LBKはアジア西部の新石器時代一括要素への適応を通じて、在来の中石器時代狩猟採集民集団により、境界での接触および文化的拡散を通じて確立された、と示唆されます。
統合モデルでは、LBKの形成は植民化・境界移動および接触のようなメカニズムを通じての、中石器時代狩猟採集民の農耕牧畜生活様式への統合として説明されます。このモデルでは、スタルチェヴォ・ケレス・クリシュ(Starčevo-Körös-Criş)文化(SKC)と関連する小集団、おそらくはヨーロッパにおけるLBKの先行者が、祖先の大半がアナトリア半島から早期に到来したバルカンの故地を離れ、北西部へと新たな領域に定着した、と想定されます。在来の中石器時代集団との接触および生産物の交換は、狩猟採集民の農耕民共同体への同化をもたらし、そこでは農耕慣行が採用されました。そうした相互作用の証拠は、明確なSKC関連墓地のあるハンガリー北東部のティスザスゼレス・ドマハザ(Tiszaszőlős-Domaháza)遺跡に存在し、ほぼ狩猟採集民の遺伝的系統の個体群の埋葬を含みます。
移民モデルでは、中石器時代ヨーロッパ中央部の低人口密度地域が、SKC文化と関連する開拓者の農耕牧畜集団に取って代わられ、しだいに在来の狩猟採集民集団は撤退させられていき、その狩猟採集民は到来してきたスタルチェヴォ移住民に顕著な影響を与えなかった、と想定されます。このモデルでは、在来集団の物質文化特徴の組み込みなしに、新たな到来者たちはその祖先的物質文化を新たに定着した領域で複製した、と想定されます。新たな環境と資源への技術革新と適応に起因して、いくつかの変異的文化が生まれ、石器や土器や建築物質技術において変化が見られます。同時に、装飾的デザインのような象徴的な体系は変化ないままでした。このモデルは1950年代末に出現し、20世紀後半に広く支持されました。
現在まで、古代DNA研究により、新石器時代ヨーロッパ農耕民集団は、おもにアナトリア半島中央部および西部の新石器時代農耕民(ANF)の遺伝的子孫と示されてきました。それらの遺伝的痕跡は、在来のヨーロッパ中央部の中石器時代狩猟採集民(WHG)とは、ミトコンドリアDNA(mtDNA)やY染色体のような単系統でもゲノム規模でも異なります。それにも関わらず、新たな到来者たちが文化的および遺伝的に在来の狩猟採集民とどの程度相互作用したのか、不明確なままです。つまり、統合主義もしくは移行主義のモデルがどの程度正確なのか、不明です。
遺伝的に、新石器時代ヨーロッパ中央部農耕民はWHG集団の遺伝的祖先特徴をわずかしか有していませんが、アナトリア半島農耕民のヨーロッパ新石器時代の子孫の遺伝子プールにおけるWHG混合の程度と時期は、ヨーロッパ中央部全域で様々です。ヨーロッパ新石器時代農耕民におけるWHG系統の量は、現在のハンガリーとドイツと他のヨーロッパ地域において、新石器時代を通じて増加する傾向にありますが(関連記事)、その混合の最初の程度は未解決のままです。それは部分的には、新石器時代農耕民移住民の最初の段階と同時代の人類遺骸が少ないためです。
オーストリアのウィーン南方に位置する、ブルン複合遺跡(the Brunn am Gebirge, Wolfholz archaeological complex)の一部であるブルン2(Brunn 2)遺跡は、オーストリアで最古の新石器時代遺跡で、ヨーロッパ中央部でも最古級となります。ブルン2遺跡はLBKの最初の段階に区分され、形成期と呼ばれています。放射性炭素年代測定法では、ブルン2遺跡の較正年代は紀元前5670~紀元前5350年前です。ブルン2遺跡形成期のおもな特徴は、洗練された土器の欠如と明確なスタルチェヴォ文化の特徴を有する粗放な土器の使用です。ヨーロッパ最初の農耕民の文化的属性の形成におけるアナトリア半島からの移民の主導的役割は、ブルン2遺跡の人工物の比較類型学的分析で明らかです。
ブルン2遺跡では、土器や石器といった豊富な人工物とともに人形や笛などいくつかの象徴的遺物が発見されており、儀式活動が行なわれた大規模なLBK共同体の「中央集落」の一部だった、と示唆されています。ブルン2遺跡では4人の埋葬が確認されています。この4人の被葬者は、個体1(I6912)・2(I6913)・3(I6914)・4(I6915)です。4人全員の歯はひじょうに摩耗しており、おもに植物性食料を摂取していた、と示唆されます。放射性炭素年代測定法により、この4人の年代はブルン複合遺跡の最初の段階と確認され、最初のヨーロッパ中央部新石器時代農耕民となります。本論文は、これら4人の遺伝子と同位体を分析し、その起源と食性と移動を調査します。
ブルン2遺跡の4人のうち、3人(I6912・I6913・I6914)で有用な遺伝的データが得られました。この3人は男性で、mtDNAハプログループ(mtHg)は、I6912がJ1、I6913がU5a1、I6914がK1b1aです。Y染色体ハプログループ(YHg)は、I6912がBT、I6913がCT、I6914がG2a2a1aです。I6912とI6913は、網羅率が低いため、さらに詳細に区分できませんでした。全ゲノム配列からは一塩基多型データが得られました。核DNAの網羅率は、I6912が0.035倍、I6913が0.006倍、I6914が0.497倍です。
主成分分析では、I6912とI6914がANFおよびそれと密接に関連するヨーロッパ新石器時代農耕民(ENF)と集団化する一方で、I6913はWHGに最も近いものの、ENFとANFにより近づいています。ただ、I6912とI6913、とくに後者の網羅率は低いので、主成分分析における位置づけには注意する必要がある、と本論文は指摘します。また本論文はf統計を使用し、WHG関連系統の割合を、I6912は12±3%、I6913は57±8%、I6914は1%未満と推定しています。さらに、I6912のWHG系統は、ヨーロッパ南東部の狩猟採集民よりもヨーロッパ西部および中央部の狩猟採集民の方に近い、と推定されました。I6913では、この区別が正確にはできませんでした。I6914はANFおよびヨーロッパ南東部のスタルチェヴォ関連個体群とほぼ対称的に関連している一方で、ヨーロッパ中央部の他のLBK集団と過剰にアレル(対立遺伝子)を共有しています。I6914と他のLBK個体群とのI6912およびI6913と比較しての高い遺伝的類似性は、これまでに研究されているアナトリア半島およびヨーロッパ南東部の農耕民とは共有されないわずかな遺伝的浮動を経験した集団出身である、と示唆します。
安定同位体(炭素13および窒素15)の有用なデータが得られたのはI6914と I6915です。安定同位体分析では、ブルン2遺跡個体群はアナトリア半島およびヨーロッパの新石器時代農耕民の範囲内に収まり、C3植物もしくはそれを食べた草食動物におもに食資源を依存していた、と推定されます。歯のエナメル質のストロンチウム同位体分析では、幼児期の場所が推定されます。ブルン2遺跡個体群では、I6912はブルン2遺跡一帯の出身で、I6913は外部出身と推定されます。同位体分析では、窒素15の値から、より新しい個体で高くなっている、と示されます。これは、ヨーロッパ中央部の初期農耕民ではアナトリア半島の栽培植物が気候の違いから不作となり、農作物よりも動物性タンパク質に依存するようになったことを示しているかもしれません。あるいは、家畜の増加によるものである可能性もあります。
ヨーロッパ中央部の早期新石器時代の遺骸は比較的豊富ですが、同時期の狩猟採集民の遺骸はほとんど知られていないので、とくにヨーロッパ新石器化の最初期段階における、狩猟採集民の生活様式や拡散してきたアナトリア半島起源の農耕民との統合に関する理解は難しくなっています。上述のように、この時期のヨーロッパにおける拡散してきた農耕民と在来の狩猟採集民との混合は限定的と推測されていますが、農耕民から狩猟採集民への生産物の流通が報告されてきており、両者の交流は少なくとも一定以上存在した、と考えられます。
ブルン2遺跡個体群の生物考古学的分析は、ヨーロッパ中央部における早期ENFの生活史の推測を可能とし、アナトリア半島起源の移住してきた農耕民と在来の狩猟採集民との間の最初期の相互作用の証拠を提示します。ブルン2遺跡の3人のmtHgのうち2系統(J1とK1b1a)は、近東の新石器時代個体群とヨーロッパにおけるその子孫たちで一般的に見られるものです。一方、I6913のmtHg-U5a1などU5は、ヨーロッパの狩猟採集民に特徴的と考えられてきましたが、最近、アナトリア半島中央部のチャタルヒュユク(Çatalhöyük)遺跡の個体でU5b2の個体が確認されています(関連記事)。I6914はYHg-G2a2a1aで、ANFおよびENF集団に特徴的なG2a系統です。カルパチア盆地とヨーロッパ南東部の早期新石器時代遺跡のLBK関連遺骸の以前の研究では、YHg-G2aは早期ENFで優勢とされています。同時期のmtHgの高い多様性は、早期LBK共同体におけるYHgの減少を示唆します。しかし、LBK集団における性比偏りの証拠は見つかっていません(関連記事)。
I6914は遺伝的にはほぼANF関連系統となり、LBKやSKCを含むENF関連個体のほぼ全員と一致します。WHG関連系統は、I6914にはほとんど見られませんが、I6912とI6913では確認されます。アナトリア半島起源のANF関連系統の個体群がヨーロッパ南東部を経てヨーロッパ中央部に拡散する過程で、ヨーロッパ在来のWHG関連系統個体群と混合したと考えられますが、その年代は特定できなかったので、それがブルン2遺跡一帯とヨーロッパ中央部到来前のどちらで起きたのか、不明です。あるいは、ブルン2遺跡の移民が、ブルン2遺跡よりも600年ほど早くアナトリア半島からバルカン半島に移住し、WHG関連系統個体群を取り込んだ集団と交流したか、その子孫だった可能性も考えられます。しかし、I6913ではWHG関連系統の割合が高く、それよりは低いものの有意にWHG関連系統を有するI6912では、ヨーロッパ南東部よりもヨーロッパ西部および中央部の狩猟採集民の方に近いWHG関連系統と推定されているので、ヨーロッパ中央部に到達してからの最近の混合である可能性が高そうです。
ブルン2遺跡の石器群からは、早期ENFと在来の狩猟採集民との間の活発な相互作用が示唆されます。ブルン複合遺跡の15000個におよぶ石器群は、新石器時代農耕民が在来の狩猟採集民との交易のために狩猟用として製作した、と本論文は推測しています。早期LBK地域では暴力の証拠が欠如しており、LBK農耕民と在来の狩猟採集民との間の関係が悪くなかったことを示唆します。上述のようにI6913は外部出身と考えられますが、I6913の墓の石器は形成期LBKの集落が発見されたハンガリーのバラトン湖(Lake Balaton)の石で製作されたので、I6913はバラトン湖周辺地域の出身かもしれません。また本論文は、ブルン2遺跡の石器の製作には狩猟採集民社会出身者が関わっている可能性を指摘しており、そうだとすると、I6913における高い割合のWHG関連系統を説明できるかもしれません。
ブルン2遺跡は、ヨーロッパ中央部における最初期新石器時代農耕民が、在来の狩猟採集民と文化的・遺伝的にどのような関係を築いたのか、検証するための格好の事例を提供します。ブルン2遺跡の個体群は、ANFとWHGの混合の第一世代だった可能性もあり、ヨーロッパ中央部の最初期新石器時代において、外来の農耕民集団と在来の狩猟採集民との間で、後期新石器時代程ではないとしても、一定以上混合が進んでいたことを示唆します。本論文は、ヨーロッパにおける新石器化に関して、アナトリア半島からの農耕民の移住が関わっており、さまざまな程度で在来集団と混合していった、と指摘します。この在来集団には、狩猟採集民だけではなく、より早期にバルカン半島に進出していた農耕民も含まれるかもしれません。今後の課題は、ヨーロッパにおけるANFと在来の狩猟採集民との関係をより広範に調査することです。
参考文献:
Nikitin AG. et al.(2019): Interactions between earliest Linearbandkeramik farmers and central European hunter gatherers at the dawn of European Neolithization. Scientific Reports, 9, 19544.
https://doi.org/10.1038/s41598-019-56029-2
統合モデルでは、LBKの形成は植民化・境界移動および接触のようなメカニズムを通じての、中石器時代狩猟採集民の農耕牧畜生活様式への統合として説明されます。このモデルでは、スタルチェヴォ・ケレス・クリシュ(Starčevo-Körös-Criş)文化(SKC)と関連する小集団、おそらくはヨーロッパにおけるLBKの先行者が、祖先の大半がアナトリア半島から早期に到来したバルカンの故地を離れ、北西部へと新たな領域に定着した、と想定されます。在来の中石器時代集団との接触および生産物の交換は、狩猟採集民の農耕民共同体への同化をもたらし、そこでは農耕慣行が採用されました。そうした相互作用の証拠は、明確なSKC関連墓地のあるハンガリー北東部のティスザスゼレス・ドマハザ(Tiszaszőlős-Domaháza)遺跡に存在し、ほぼ狩猟採集民の遺伝的系統の個体群の埋葬を含みます。
移民モデルでは、中石器時代ヨーロッパ中央部の低人口密度地域が、SKC文化と関連する開拓者の農耕牧畜集団に取って代わられ、しだいに在来の狩猟採集民集団は撤退させられていき、その狩猟採集民は到来してきたスタルチェヴォ移住民に顕著な影響を与えなかった、と想定されます。このモデルでは、在来集団の物質文化特徴の組み込みなしに、新たな到来者たちはその祖先的物質文化を新たに定着した領域で複製した、と想定されます。新たな環境と資源への技術革新と適応に起因して、いくつかの変異的文化が生まれ、石器や土器や建築物質技術において変化が見られます。同時に、装飾的デザインのような象徴的な体系は変化ないままでした。このモデルは1950年代末に出現し、20世紀後半に広く支持されました。
現在まで、古代DNA研究により、新石器時代ヨーロッパ農耕民集団は、おもにアナトリア半島中央部および西部の新石器時代農耕民(ANF)の遺伝的子孫と示されてきました。それらの遺伝的痕跡は、在来のヨーロッパ中央部の中石器時代狩猟採集民(WHG)とは、ミトコンドリアDNA(mtDNA)やY染色体のような単系統でもゲノム規模でも異なります。それにも関わらず、新たな到来者たちが文化的および遺伝的に在来の狩猟採集民とどの程度相互作用したのか、不明確なままです。つまり、統合主義もしくは移行主義のモデルがどの程度正確なのか、不明です。
遺伝的に、新石器時代ヨーロッパ中央部農耕民はWHG集団の遺伝的祖先特徴をわずかしか有していませんが、アナトリア半島農耕民のヨーロッパ新石器時代の子孫の遺伝子プールにおけるWHG混合の程度と時期は、ヨーロッパ中央部全域で様々です。ヨーロッパ新石器時代農耕民におけるWHG系統の量は、現在のハンガリーとドイツと他のヨーロッパ地域において、新石器時代を通じて増加する傾向にありますが(関連記事)、その混合の最初の程度は未解決のままです。それは部分的には、新石器時代農耕民移住民の最初の段階と同時代の人類遺骸が少ないためです。
オーストリアのウィーン南方に位置する、ブルン複合遺跡(the Brunn am Gebirge, Wolfholz archaeological complex)の一部であるブルン2(Brunn 2)遺跡は、オーストリアで最古の新石器時代遺跡で、ヨーロッパ中央部でも最古級となります。ブルン2遺跡はLBKの最初の段階に区分され、形成期と呼ばれています。放射性炭素年代測定法では、ブルン2遺跡の較正年代は紀元前5670~紀元前5350年前です。ブルン2遺跡形成期のおもな特徴は、洗練された土器の欠如と明確なスタルチェヴォ文化の特徴を有する粗放な土器の使用です。ヨーロッパ最初の農耕民の文化的属性の形成におけるアナトリア半島からの移民の主導的役割は、ブルン2遺跡の人工物の比較類型学的分析で明らかです。
ブルン2遺跡では、土器や石器といった豊富な人工物とともに人形や笛などいくつかの象徴的遺物が発見されており、儀式活動が行なわれた大規模なLBK共同体の「中央集落」の一部だった、と示唆されています。ブルン2遺跡では4人の埋葬が確認されています。この4人の被葬者は、個体1(I6912)・2(I6913)・3(I6914)・4(I6915)です。4人全員の歯はひじょうに摩耗しており、おもに植物性食料を摂取していた、と示唆されます。放射性炭素年代測定法により、この4人の年代はブルン複合遺跡の最初の段階と確認され、最初のヨーロッパ中央部新石器時代農耕民となります。本論文は、これら4人の遺伝子と同位体を分析し、その起源と食性と移動を調査します。
ブルン2遺跡の4人のうち、3人(I6912・I6913・I6914)で有用な遺伝的データが得られました。この3人は男性で、mtDNAハプログループ(mtHg)は、I6912がJ1、I6913がU5a1、I6914がK1b1aです。Y染色体ハプログループ(YHg)は、I6912がBT、I6913がCT、I6914がG2a2a1aです。I6912とI6913は、網羅率が低いため、さらに詳細に区分できませんでした。全ゲノム配列からは一塩基多型データが得られました。核DNAの網羅率は、I6912が0.035倍、I6913が0.006倍、I6914が0.497倍です。
主成分分析では、I6912とI6914がANFおよびそれと密接に関連するヨーロッパ新石器時代農耕民(ENF)と集団化する一方で、I6913はWHGに最も近いものの、ENFとANFにより近づいています。ただ、I6912とI6913、とくに後者の網羅率は低いので、主成分分析における位置づけには注意する必要がある、と本論文は指摘します。また本論文はf統計を使用し、WHG関連系統の割合を、I6912は12±3%、I6913は57±8%、I6914は1%未満と推定しています。さらに、I6912のWHG系統は、ヨーロッパ南東部の狩猟採集民よりもヨーロッパ西部および中央部の狩猟採集民の方に近い、と推定されました。I6913では、この区別が正確にはできませんでした。I6914はANFおよびヨーロッパ南東部のスタルチェヴォ関連個体群とほぼ対称的に関連している一方で、ヨーロッパ中央部の他のLBK集団と過剰にアレル(対立遺伝子)を共有しています。I6914と他のLBK個体群とのI6912およびI6913と比較しての高い遺伝的類似性は、これまでに研究されているアナトリア半島およびヨーロッパ南東部の農耕民とは共有されないわずかな遺伝的浮動を経験した集団出身である、と示唆します。
安定同位体(炭素13および窒素15)の有用なデータが得られたのはI6914と I6915です。安定同位体分析では、ブルン2遺跡個体群はアナトリア半島およびヨーロッパの新石器時代農耕民の範囲内に収まり、C3植物もしくはそれを食べた草食動物におもに食資源を依存していた、と推定されます。歯のエナメル質のストロンチウム同位体分析では、幼児期の場所が推定されます。ブルン2遺跡個体群では、I6912はブルン2遺跡一帯の出身で、I6913は外部出身と推定されます。同位体分析では、窒素15の値から、より新しい個体で高くなっている、と示されます。これは、ヨーロッパ中央部の初期農耕民ではアナトリア半島の栽培植物が気候の違いから不作となり、農作物よりも動物性タンパク質に依存するようになったことを示しているかもしれません。あるいは、家畜の増加によるものである可能性もあります。
ヨーロッパ中央部の早期新石器時代の遺骸は比較的豊富ですが、同時期の狩猟採集民の遺骸はほとんど知られていないので、とくにヨーロッパ新石器化の最初期段階における、狩猟採集民の生活様式や拡散してきたアナトリア半島起源の農耕民との統合に関する理解は難しくなっています。上述のように、この時期のヨーロッパにおける拡散してきた農耕民と在来の狩猟採集民との混合は限定的と推測されていますが、農耕民から狩猟採集民への生産物の流通が報告されてきており、両者の交流は少なくとも一定以上存在した、と考えられます。
ブルン2遺跡個体群の生物考古学的分析は、ヨーロッパ中央部における早期ENFの生活史の推測を可能とし、アナトリア半島起源の移住してきた農耕民と在来の狩猟採集民との間の最初期の相互作用の証拠を提示します。ブルン2遺跡の3人のmtHgのうち2系統(J1とK1b1a)は、近東の新石器時代個体群とヨーロッパにおけるその子孫たちで一般的に見られるものです。一方、I6913のmtHg-U5a1などU5は、ヨーロッパの狩猟採集民に特徴的と考えられてきましたが、最近、アナトリア半島中央部のチャタルヒュユク(Çatalhöyük)遺跡の個体でU5b2の個体が確認されています(関連記事)。I6914はYHg-G2a2a1aで、ANFおよびENF集団に特徴的なG2a系統です。カルパチア盆地とヨーロッパ南東部の早期新石器時代遺跡のLBK関連遺骸の以前の研究では、YHg-G2aは早期ENFで優勢とされています。同時期のmtHgの高い多様性は、早期LBK共同体におけるYHgの減少を示唆します。しかし、LBK集団における性比偏りの証拠は見つかっていません(関連記事)。
I6914は遺伝的にはほぼANF関連系統となり、LBKやSKCを含むENF関連個体のほぼ全員と一致します。WHG関連系統は、I6914にはほとんど見られませんが、I6912とI6913では確認されます。アナトリア半島起源のANF関連系統の個体群がヨーロッパ南東部を経てヨーロッパ中央部に拡散する過程で、ヨーロッパ在来のWHG関連系統個体群と混合したと考えられますが、その年代は特定できなかったので、それがブルン2遺跡一帯とヨーロッパ中央部到来前のどちらで起きたのか、不明です。あるいは、ブルン2遺跡の移民が、ブルン2遺跡よりも600年ほど早くアナトリア半島からバルカン半島に移住し、WHG関連系統個体群を取り込んだ集団と交流したか、その子孫だった可能性も考えられます。しかし、I6913ではWHG関連系統の割合が高く、それよりは低いものの有意にWHG関連系統を有するI6912では、ヨーロッパ南東部よりもヨーロッパ西部および中央部の狩猟採集民の方に近いWHG関連系統と推定されているので、ヨーロッパ中央部に到達してからの最近の混合である可能性が高そうです。
ブルン2遺跡の石器群からは、早期ENFと在来の狩猟採集民との間の活発な相互作用が示唆されます。ブルン複合遺跡の15000個におよぶ石器群は、新石器時代農耕民が在来の狩猟採集民との交易のために狩猟用として製作した、と本論文は推測しています。早期LBK地域では暴力の証拠が欠如しており、LBK農耕民と在来の狩猟採集民との間の関係が悪くなかったことを示唆します。上述のようにI6913は外部出身と考えられますが、I6913の墓の石器は形成期LBKの集落が発見されたハンガリーのバラトン湖(Lake Balaton)の石で製作されたので、I6913はバラトン湖周辺地域の出身かもしれません。また本論文は、ブルン2遺跡の石器の製作には狩猟採集民社会出身者が関わっている可能性を指摘しており、そうだとすると、I6913における高い割合のWHG関連系統を説明できるかもしれません。
ブルン2遺跡は、ヨーロッパ中央部における最初期新石器時代農耕民が、在来の狩猟採集民と文化的・遺伝的にどのような関係を築いたのか、検証するための格好の事例を提供します。ブルン2遺跡の個体群は、ANFとWHGの混合の第一世代だった可能性もあり、ヨーロッパ中央部の最初期新石器時代において、外来の農耕民集団と在来の狩猟採集民との間で、後期新石器時代程ではないとしても、一定以上混合が進んでいたことを示唆します。本論文は、ヨーロッパにおける新石器化に関して、アナトリア半島からの農耕民の移住が関わっており、さまざまな程度で在来集団と混合していった、と指摘します。この在来集団には、狩猟採集民だけではなく、より早期にバルカン半島に進出していた農耕民も含まれるかもしれません。今後の課題は、ヨーロッパにおけるANFと在来の狩猟採集民との関係をより広範に調査することです。
参考文献:
Nikitin AG. et al.(2019): Interactions between earliest Linearbandkeramik farmers and central European hunter gatherers at the dawn of European Neolithization. Scientific Reports, 9, 19544.
https://doi.org/10.1038/s41598-019-56029-2
この記事へのコメント