Common cis-eQTLs do not explain chromosome 21 expression deviations from the ploidy expectation in Down syndrome

技術的要約：一般的なシス型eQTLは、21トリソミーにおける染色体21の遺伝子発現偏差を説明しない

問題提起
ダウン症候群（DS）は、ヒト21番染色体（HSA21）の三倍体化である21トリソミー（T21）によって引き起こされる。ほとんどのHSA21遺伝子は、DS患者において、正常二倍体（D21）と比較して予想される1.5倍の発現増加を示すが、一部の遺伝子集合はこの期待値から逸脱し、二倍体に近い発現レベルを示す。Hunterらによる家系ベースの研究を含む先行研究では、一般的な遺伝的調節変異（具体的にはシス型発現量的形質遺伝子座、すなわちシス型eQTL）が、特定の個人におけるこれらの偏差を説明する可能性が示唆されていた。しかし、一般的なシス型調節変異が、集団レベルにおいてこれらの偏差を説明できるかどうかは不明なままであった。本研究では、T21集団とD21集団の間で、一般的なシス型eQTLの頻度や効果の違いが、HSA21の遺伝子量効果（dosage effects）における遺伝子特異的なバッファリング（緩衝作用）や増幅を説明するかどうかを調査する。

方法論
著者らは、Human Trisome Project (HTP) に由来する、274名のT21参加者と80名の二倍体（D21）参加者の全血バルクRNA-seqデータからなる、現在までに存在する最大規模のDS特異的トランスクリプトームデータセットを利用した。三倍体の文脈における正確性を確保するため、著者らは厳格な倍数性考慮型の解析パイプラインを採用した：

1. 倍数性を考慮したジェノタイピング： 標準的な生殖細胞系列のジェノタイピング・ワークフローは二倍体を前提としており、HSA21上のアレル・ドーズ（アレル数）を誤判定する。著者らは、Octopusを用いてWGSデータから直接バリアントを再コールし、その際、T21サンプルのHSA21の倍数性を3に、D21サンプルの倍数性を2に明示的に設定した。これにより、バリアント呼び出しのための正確なアレルカウント（0–3コピー）を生成した。
2. 倍数性を考慮した発現変動解析： DESeq2を用い、HSA21遺伝子の期待されるフォールドチェンジをデフォルトの1.0ではなく1.5（log2 1.5 ≈ 0.585）に設定したモデルで、差分的発現解析を行った。共変量には、組成の違いを制御するために、年齢、性別、BMI、サンプルソース、および質量分析法から導出された19種類の細胞型比率を含めた。
3. 偏差の特定： 遺伝子が「偏差している」と分類されるのは、倍数性補正後のフォールドチェンジが1.5から少なくとも33%以上異なる場合とした（調整済み p < 0.01）。
4. シス型eQTL関連解析： 各偏差遺伝子について、GTEx v10の全血データセットからの一般的なシス型eQTL（±1 Mbのウィンドウ）をクロスリファレンスした。T21コホート内において、共変量調整済みの発現をジェノタイプ・ドーズ（0–3コピー）に対して回帰した。有意性は、多重検定を調整した上での遺伝子レベルの置換テスト（1,000回の置換）によって評価された。
5. 比較解析： 偏差している遺伝子におけるシス型eQTLの検出率と効果サイズを、期待される1.5倍の変化を示すHSA21遺伝子と比較した。また、T21参加者とGTExドナー間のアレル頻度の違いが、数学的に観察された発現偏差を説明できるかどうかをモデル化した。

主な結果

• 発現の偏差： 低リードカウントおよび高リピート含有量に対するフィルタリングを行い、倍数性補正を適用した後、標的とした318個のHSA21遺伝子のうち、13個（8.5%）が期待される1.5倍の変化から有意に偏差した。テスト対象となった10個の偏差遺伝子のうち、5個は期待より高く、5個は期待より低かった。
• シス型eQTLの検出： テスト可能であった（GTExにデータとバリアントが存在した）10個の偏差遺伝子のうち、8個でT21コホートにおけるシス型eQTLが検出された。しかし、期待される1.5倍の変化を示した99個のHSA21遺伝子のうち、59個も有意なシス型eQTLを持っていた。
• 効果サイズと頻度： T21参加者におけるシス型eQTLの大きさは、GTExで観察されたものと同様であった（偏差遺伝子のメディアン・アレリック・フォールドチェンジは約0.73、期待される遺伝子量を示す遺伝子では約0.76）。決定的なことに、T21参加者におけるバリアントのアレル頻度は、GTExドナーのそれと系統的な差はなかった（メディアン差 －0.0005）。
• 説明力： 数学的モデリングにより、たとえT21参加者が代替アレルを3コピーすべて、あるいは全く持っていなかったとしても、観察されたアレル頻度の違いは、観察された発現偏差の規模を生み出すことはできないことが示された。最良のバリアントは、T21参加者間の発現分散のわずか2%から33%しか説明できなかった。
• 特定の事例： 期待よりも低い発現を示した2つの遺伝子（BACE2およびOLIG2）は、GTExでは強力なeQTLを有していたにもかかわらず、T21コホートではシス型eQTLが検出されなかった。しかし、検出の欠如は検出力の低さによるものではなく、彼らの偏差パターンは共変量モデル全体で堅牢であった。

本研究は、一般的なシス型調節変異は、T21におけるHSA21遺伝子に対して、二倍体におけるものと同様に作用するものの、どの特定の遺伝子が偏差するかを集団レベルで説明することはない、と結論付けている。

• 個人 vs 集団： 著者らは、個人のジェノタイプがその特定の偏差を説明できること（Hunterらによる家系ベースの研究で示されている通り）を明確にしているが、一般的なシス型eQTLは、生存出生個体としてのT21集団における系統的な差異を生み出すことはない。つまり、偏差は生存出生のT21集団における調節アレルの頻度のシフトによって引き起こされるものではない。
• 遺伝的仮説の棄却： これらの知見は、一般的なシス型調節変変異が、ダウン症における遺伝子特異的なバッファリングまたは増幅の主要なメカニズムであるという仮説を支持しない。
• 代替メカニズム： 著者らは、これらの偏差を駆動するメカニズムは他の場所に存在する可能性が高いと考えており、それはオートレギュラトリー（自己調節的）なフィードバックループ、三倍体化されたmiRNAによる転写後サイレンシング、選択的スプライシング、クロマチン・リモデリング、あるいは一般的なeQTLマッピングでは捉えきれない稀な調節バリアントといった、非遺伝的な要因である可能性を示唆している。

本研究は、適切な倍数性考慮型の手法を用いて解析すれば、大部分のHSA21遺伝子は期待される1.5倍の遺伝量増加に従うことを強調しており、少数の遺伝子で見られるこれらの偏差は、一般的なシス型調節変異を超えた調査を必要とする未解決の課題として残されている。