イベリコ交雑豚の成長、肥満、肉質形質に関するSNPの発見と関連研究

Scientific Reports volume 12、記事番号: 16361 (2022) この記事を引用

1192 アクセス

2 引用

1 オルトメトリック

メトリクスの詳細

イベリコ豚とその交配は、高品質の肉製品を得るために生産されています。この研究の目的は、筋肉の成長、脂肪、肉の品質、代謝に関連する形質との関連性について、生物学的および機能的基準によって選択された幅広い DNA マーカーのパネルを評価することでした。全ゲノム配列決定とSNP発見のために、出生後成長パターンが異なる18頭の交雑種イベリコ豚を使用し、1,300万以上の変異が検出されました。我々は、注釈付き遺伝子上に位置し、著しく異なる成長パターンを持つブタ間で異なる対立遺伝子頻度を示す 1,023 個のミスセンス SNP を選択しました。このパネルを、文献マイニングおよび筋肉 RNAseq データから得られた 192 個の候補 SNP で補完しました。選択されたマーカーは、商業集団からの 480 頭のイベリコ × デュロック豚で遺伝子型決定され、表現型が得られ、分離を示す遺伝子型決定に成功した 1005 頭の SNP について関連研究が実行されました。その結果、候補遺伝子（LEPR、ACACA、FTO、LIPE、SCDなど）におけるいくつかの既知のSNPが肥満、成長、脂肪酸組成に及ぼす影響が確認され、LRIG3、DENND1Bなどのあまり知られていない遺伝子における新しいSNPの興味深い影響も明らかになりました。、SOWAHB、EPHX1、またはNFE2L2は体重、さまざまな年齢における1日の平均増加量および肥満度に影響を与え、KRT10、NLE1、KCNH2、またはAHNAKは肥満度とFA組成に影響を与えます。この結果は、ブタにおけるマーカー支援選抜戦略の将来の実施のための貴重な基礎を提供し、関連形質の遺伝的構造のより良い理解に貢献する。

イベリコ豚は、脂肪生成の可能性が高く、食欲があり、優れた肉品質を特徴とする、自生の脂肪の多い品種です。これらの特徴は、その遺伝学と伝統的な大規模な生産システムの結果であり、両方ともオレイン酸と抗酸化物質を多く含む皮下脂肪および筋肉内脂肪の沈着に寄与しています1,2。イベリコ豚×デュロック交雑種豚は、集中的な生産システムでイベリコ肉製品を生産するために使用される主な遺伝子型ですが、感覚的な肉の品質は純粋なイベリコ豚よりも低いと考えられています3、4、5。さらに、交雑されたイベリコ豚は、非常に不均一な発育、生産性、および品質の形質を特徴とし、根底にある遺伝子を同定する大きな可能性を秘めています。関連する形質の分子遺伝学的基礎に関する知識は、マーカー支援選抜プログラムの設計に役立つであろう。

生産的形質の遺伝的構造を深化させるために、主に純血種または実験的交配を用いて、低スケールの構造的アプローチを用いて、イベリコ豚集団においてさまざまな候補遺伝子研究が行われてきました。関連する遺伝子と変異はこれまでこのアプローチで同定されており、主に LEPR7,8、ELOVL69、ACSL410,11 などの IBMAP 実験集団で検出されたものです。最近、Fernandez-Barroso et al.12 は、PRKAG3、CAPN1、および他のいくつかの遺伝子の変異が純粋なイベリコ動物の肉の品質形質に有意に関連していることを発見しました。また、機能的遺伝研究がイベリコ豚で行われ、表現型の変異に潜在的に関与する機能的候補遺伝子および経路に関する情報が提供されている5、13、14、15。一方、ゲノムワイド関連研究 (GWAS) 用の市販の高密度 SNP アレイの使用は、複雑な形質の遺伝的基盤をさらに深く調べるための系統的かつ強力なアプローチを提供します。このアプローチにより、イベリア人を含むさまざまな集団や品種において、多くの興味深いゲノム領域と候補遺伝子が同定されました16、17、18、19、20、21、22。それにもかかわらず、市販の SNP アレイは遺伝的変異の一部しか識別できないため、その適用性は限られています。これは、それらが少数の国際品種で検出された SNP に基づいて設計されているため、特に在来品種の分析では情報量と検出力が限られているためです 22,23。対照的に、全ゲノム配列決定 (WGS) は、原因となるものを含むすべての遺伝的変異を検出できる可能性があり 24、現地品種の関連形質の遺伝的基盤を明らかにするのに役立つ、カスタマイズされた有益な SNP パネルの設計が可能になる可能性があります。このような進歩により、非国際的豚品種の関連する生産的および品質的形質の遺伝的改善が達成可能になります。

この研究の目的は、イベリコ交雑種豚の筋肉の成長、脂肪の沈着と組成、代謝、肉質に関連する形質と関連する幅広い DNA マーカーのパネルを発見し、評価することでした。この目的のために、全ゲノム配列データ、トランスクリプトーム配列データ、候補遺伝子アプローチの組み合わせから取得および選択された SNP を使用して設計された中密度カスタムジェノタイピングプロトコルが使用されました。

この研究は、研究に使用される動物の保護に関する欧州連合指令 2010/63/UE に準拠したスペインの動物保護政策 RD53/2013 に従って実施されました。この実験は、INIA の施設内動物管理使用委員会 (IACUC) と呼ばれる INIA 動物研究倫理委員会によって特別に評価され、承認されました (レポート CEEA 2012/036)。この研究は、ARRIVE ガイドライン (https://arriveguidelines.org/arrive-guidelines) に従って実施されました。

イベリア人×デュロック人の交雑母集団を使用しました。動物は商業農場、イベリコスデアラウゾ 2004 SL (スペイン、サラマンカ、ゾリタデラフロンテーラ) で飼育されました。この研究には、47頭の3回目および4回目の出産のイベリコ豚から生まれた合計477頭の交雑種子豚(雌50%、雄50%)が参加した。イベリコ豚には、Duroc PIC 雄豚 (Genus plc、英国) からの冷却された精液が授精されました。雌豚は電子耳標で個別に識別され、妊娠101日目までグループで飼育され、その後哺乳期が終わるまで個別の囲いに移された。 De Blas et al.25 の推奨事項に基づいて、個々の 1 日の必要量を満たすように調整されたイベリコ品種の標準食を雌豚と子豚に与えました。子豚には出生時に個別にタグが付けられ、離乳するまで雌豚と一緒に過ごしました。生後 2 日以内に、交雑養育によって同腹子が均等化され、雄子豚は外科的に去勢されました。アザペロン（1 mg/kg LW; Stresnil®、Ecuphar Veterinaria、SLU、スペイン、バルセロナ）による鎮静および局所麻酔（塩酸リドカイン 2% 溶液; Lidocaina Normon®、20 mg/mL、Laboratorios Normon、マドリッド）後に外科的精巣切除術を実施しました。、スペイン）皮膚を通して精索に直接挿入された25G針を使用して適用されます。子豚は平均生後 24 日で離乳され、移行期には性別と体重 (BW) ごとに分配された 12 頭/囲いの子豚のグループで飼育されました。成長肥育期（生後 70 日から屠殺まで）では、動物を性別と体重別に 40 頭/檻のグループに分けて飼育した。雌は、それぞれ生後120日と148日に2回のワクチン接種によるVACSINCEL（Zoetis Inc.、米国ニュージャージー州）を使用した免疫去勢を受けた。生後 240 日目から 340 日目まで、豚は最低市場体重に達するたびに市場に出され、イベリコ交雑種豚の場合は枝肉重量 115 kg に設定されています。生後 340 日目に、残りの豚は体重に関係なくすべて市場に送られました。表現型の評価は、この研究に関与した動物の一部を使用して、以前の論文ですでに記載されている方法に従いました6。

さまざまな測定を以下の時点で実施した：出生、離乳（生後24日）、平均年齢の70、110、150、180、215および240日目、および屠殺時。体重はこれらすべての時点で個別に測定されました。形態学的測定値（後頭鼻の長さ、両頭頂骨の直径、体幹の長さ、腹部および胸部の周囲および胸部の最大直径）を、すべての子豚の出生時および離乳時にメジャーテープを使用して記録した。背脂肪の厚さとロース肉の直径は、離乳時と生後 110 日および 215 日に測定されました。屠殺場では枝肉の重量と長さ、背脂肪の厚さを記録した。

平均 1 日体重増加 (ADWG) は、哺乳期 (1 ～ 24 歳)、移行期 (25 ～ 70 歳)、および生後 240 日 (70 歳) までの成長期から肥育期の 5 つの選択された期間中に計算されました。生後 –110、111 ～ 150、151 ～ 180、181 ～ 215、および 216 ～ 240 日。各期間はその最終日によって名前が付けられます。飼料転換率 (FCR) は、成長期から肥育期の選択された 5 つの期間中にペンで計算されました。一日の平均体重増加と FCR も、それぞれ出生から屠殺までと離乳から屠殺までの全期間について計算されました。

血液サンプルは生体内 (110 日および 215 日齢) および屠殺時に採取され、DNA 抽出と代謝状態の評価に使用されました。背最長筋（LD）、肝右葉、皮下脂肪のサンプルを屠殺時に採取し、脂肪酸（FA）組成が分析されるまで-20℃でバイオバンクしました。筋肉のドリップロス分析も実施されました26。

代謝状態は、生後 110 日および 215 日の時点と屠殺時に評価されました。血液サンプルは、グルコース（グルコースおよびフルクトサミン）の代謝および脂質プロファイル（総コレステロール、高密度リポタンパク質コレステロール（HDL-c）、低密度リポタンパク質コレステロール（LDL-c）およびトリグリセリド）の代謝に関連するパラメータを決定するために分析されました。臨床化学分析装置（Saturno 300 plus、Crony Instruments srl、ローマ、イタリア）による。

肝脂肪 (LF) および筋肉内脂肪 (IMF) は、Segura および Lopez-Bote27 の記載に従って抽出され、乾物のパーセンテージとして表されました。これらの脂質は、脂肪組織の 2 つの主要な部分、つまり中性脂質 (NL) と極性脂質 (PL) に分別されました 1,27。皮下脂肪を抽出し、外層と内層に分けて個別に分析しました。個々の FA 値、飽和、一価不飽和、多価不飽和 FA (SFA、MUFA、PUFA) の割合、不飽和指数 (UI)、および n-3 FA (n-3) と n-6 FA の合計の合計 ( n-6) とその比 (n-6/n-3) を計算しました28。さらに、ステアロイル CoA デサチュラーゼ酵素 1 (SCD1) の活性は、2 つの不飽和指数、C18:1/C18:0 および MUFA/SFA29 の比を使用して推定されました。

交雑したイベリコ豚集団の中から、出生後成長パターンが異なる 18 頭の雄豚が全ゲノム配列決定のために選択されました。 240 日の出生体重と体重は、低出生体重および低最終体重 (LL; 6 匹)、低出生体重だが正常最終体重 (LN; 6 匹)、および正常値の 3 つの異なる成長カテゴリーに対応する動物を選択するために考慮されました。出生体重と最終体重の両方 (NN; 6 匹) (図 1)。選択された動物は、さまざまな同腹仔から生まれました。

WGS 用に選択された 18 匹の動物の成長パターン。各成長グループには 6 匹の動物が含まれます: 低出生体重児および低最終体重 (LL; ピンク)。低出生体重と正常最終体重（LN; 紫）、および出生体重と最終体重の両方の正常値（NN; 緑色）。母集団の平均は黒で表示されます。 BWの体重。

NucleoSpin® Blood Kit (Macherey-Nagel、デューレン、ドイツ) を使用して、すべての動物の血液サンプルからゲノム DNA を取得しました。 NanoDrop 2000 分光光度計 (Thermo Fisher Scientific、米国マサチューセッツ州ウォルサム) を使用して、DNA の濃度と品質を測定しました。選択した 18 匹の動物の全ゲノム再配列決定を、CNAG-CRG (スペイン、バルセロナ) の Illumina HiSeq 2000 プラットフォームで実行しました。約 1 μg のゲノム DNA がランダムに断片化されました。ゲノムライブラリーは、TruSeq DNA サンプル調製キットを使用して、製造業者のプロトコール (Illumina、True Seq DNA 調製ガイド) に従って調製されました。ペアエンドライブラリーは、v4 ケミストリーと 2 × 125 リードを使用して Illumina HiSeq 2000 でシーケンスされました。

FastQC ソフトウェアを使用して、生データの品質管理を実行しました (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/)。 Trimmomatic パッケージバージョン 0.3830 は、ペアエンド読み取りトリミングに使用されました。リードをトリミングしてアダプターを除去し、リードの開始または終了で低品質または N 塩基を除去し、4 塩基スライディングウィンドウで 3' 末端からリードをスキャンし、平均 Phred 品質が 20 を下回ったときにカットしました。トリミングにより、リードが 40 bp より短い場合は削除されました。フィルタリングされたペアリードは、「bwa-mem」アルゴリズムを採用した Burrows-Wheeler Aligner (BWA バージョン 0.7.10)32 によってブタ参照ゲノムビルド 1131 とアライメントされました。 SAMtools バージョン 1.633 は、マップされていないリードまたは不適切にマップされたリードのフィルタリング、並べ替え、インデックス付け、リードグループの追加に使用されました。 Picard ツールバージョン 2.18.9 (http://broadinstitute.github.io/picard/) を使用して重複をマークしました。このステップからバリアントコールフォーマット (VCF) ファイルまで、読み取りはゲノム解析ツールキット (GATK バージョン 4.0.6.0) を使用して処理されました 34,35。 BaseRecalibrator を使用して塩基品質が再調整され、HaplotypeCaller を使用してサンプルごとにバリアント呼び出しが実行され、その結果 Genomic Variant Call Format ファイル (GVCF) が生成され、デフォルトのパラメーター値を使用して CombineGVCF を使用して VCF ファイルに結合されました。最後に、デフォルトのパラメーターを使用して GenotypeGVCF を使用してバリアントを呼び出しました。合計 14,871,405 個と 3,613,138 個の生の SNP と INDELS が検出されました。 SNP は、GATK VariantFiltration ツールを使用して品質フィルター処理され、パラメーター「QUAL < 40 || QD < 5.0 || FS > 60.0 || MQ < 40.0 || DP < 30」を持つものは除外され、フィルター処理された SNP は 13,324,328 個になりました。 VCFtools36 を使用して、マイナー対立遺伝子頻度 (MAF) が 0.05 未満の SNP をフィルター処理し、フィルター処理された SNP の最終セット 11,756,179 を導き出しました。

Ensembl バリアント効果予測子 (VEP バージョン 99)37 を使用して、WGS データセットのバリアントに VEP データベース (マージされたキャッシュファイル Sscrofa11.1 バージョン 99) 内の情報を注釈付けしました。 VEP の -- distance 5 kb のデフォルトパラメータは、上流と下流のバリアントを定義するために適用されました。また、ミスセンスバリアントの SIFT スコア 38 の予測に --sift b オプションを適用しました。

私たちは、注釈付き遺伝子上に位置し、増殖パターンが異なるグループ間で対立遺伝子頻度の違いを示すミスセンス SNP に焦点を当てました。 VEPによると、1170万個のうち、合計3万2688個のマーカーがミスセンスとして分類された。 LD プルーニングは PLINK 1.939 で実装され、1000 kb ウィンドウで r2 > 0.6 の合計 20,049 個の SNP を削除しました。最後に、常染色体上にマッピングされた合計 8,419 個の SNP が下流解析に使用されました。成長グループ内の対立遺伝子頻度は、PLINK 1.9 を使用して計算されました。少なくとも 1 つの比較 (LL 対 LN、LL 対 NN、または LN 対 NN) において、グループ間の対立遺伝子頻度の差が 0.25 を超える SNP を維持しました。その結果、合計 2259 個の SNP が残りました。

これらの事前に選択された SNP のうち、1023 は Agriseq の遺伝子型解析に技術的に適切でした。これらはデザインに含まれており、イベリア × デュロック集団における Genotyping by Sequencing (GBS) によって遺伝子型が特定され、そこから表現型も入手できました。 GBS によって遺伝子型決定されたマーカーのリストは、補足表 S1 として含まれています。 GBS の場合、Agriseq HTS ライブラリーキットを使用してライブラリーを作成し、Ion GeneStudio S5 Sequencer 上の Ion 540 チップでシーケンスしました。

さらに、文献マイニングおよび純粋およびデュロック交配イベリア筋の RNAseq データから得られた 192 の候補 SNP でこのパネルを補完しました 13、14、40 (補足表 S2)。これらの SNP は、特注の TaqMan® OpenArray™ (OA) ジェノタイピングプラットフォーム (Thermo Fisher Scientific、米国マサチューセッツ州ウォルサム) を使用して、同じイベリア人 × デュロック人集団において遺伝子型決定されました。ジェノタイピングは、Veterinary Molecular Genetics Service (SVGM、UAB、スペイン、バルセロナ) の QuantStudio™ 12 K Flex リアルタイム PCR システム (Thermo Fisher Scientific、米国マサチューセッツ州ウォルサム) で実施されました。メーカーの指示に従って、DNA サンプルをアレイにロードし、増幅しました。対立遺伝子特異的シグナル強度の検出は OpenArray NT Imager を使用して実行され、遺伝子型は OpenArray SNP Genotyping 解析ソフトウェアを使用して呼び出されました。さらに、SNP 画像を視覚的に検査して、クラスタリングの問題を検出しました。

合計 SNP データセット (1215 マーカー) は、MAF 値 (0.1) による PLINK 1.9 を使用し、ジェノタイピングエラーのあるマーカーを除外して関連解析のためにフィルターにかけられ、フィルター後に 1005 個の SNP が残りました。これらの分析は、ジェノタイピングおよび表現型データを使用して 477 人の個人に対して実行されました。 SNP の影響は、ソフトウェア GCTA41 の Leave-One-chromosome-out オプションを使用して分析されました。このオプションは、次の一般モデルを使用して、混合線形モデルベースの関連分析を実装します。

ここで、y は各形質に対応する表現型値、a は平均項、b には分析された SNP の相加効果 (固定項) が含まれ、X には 0、1、または 2 としてコード化された遺伝子型指標変数が含まれ、f には分析された形質のグループに応じて異なる残りの固定効果と共変量 (補足表 S3)、および F その対応する発生率行列、g- はマッピングされたものを除くすべての SNP (ゲノム関係行列、GRM によって捕捉される) の累積効果です。は分析された SNP と同じ染色体上にあり、e は残差です。遺伝的分散 var(g-) はヌルモデル y = a + g- + e に基づいて推定され、各 SNP と形質の間の関連性がテストされながら固定されます。この分散は、GRM の計算時に染色体が除外されるたびに再推定されます。 Rstudio43 のライブラリ q-value42 を使用することで誤発見率 (FDR) を制御しました。 p 値と q 値がそれぞれ 0.05 および 0.10 未満の SNP は、テストされた形質と有意に関連しているとみなされました。

Illumina シーケンシングによりサンプルごとに平均 8,300 万のペアリードが生成され、その範囲は 6,700 万から 1 億の範囲で、平均カバレッジは 8.4 × でした。品質管理後、平均してリードの 98% が参照ゲノム (Sscrofa11.1) にマッピングされ、1,400 万を超えるバリアントが検出されました。いくつかの品質フィルターと機能的および技術的な選択基準を適用した後、最終的に 1,023 個のミスセンス多型が GBS によって遺伝子型特定されました。

GBS 品質指標により、ジェノタイピング手順の信頼性が確認されました。平均サンプルとマーカーコールは両方とも 93.9% で、マーカーの 90.2% のコール率は 90% を超えていました。平均カバレッジは 183 × で、ジェノタイピング目的で Agriseq が推奨するカバレッジ 100 × を大幅に上回りました。均一性（平均深さの少なくとも 0.2 倍を持つターゲット塩基の割合）は平均 92.9% でした。

GBS マーカーパネルに含まれる 1,023 個の SNP のうち、120 個の多型が分離の欠如 (MAF < 0.10) または遺伝子型解析エラーにより廃棄され、関連研究に有用なデータを含む 903 個の GBS マーカーが得られました。

並行して、候補遺伝子の 192 個のマーカーがオープンアレイによって遺伝子型決定されました。それらのうち、60 個は単型であり、合計 30 個は低い変動性 (0 > MAF < 0.10) または遺伝子型決定エラーにより廃棄されたため、102 個の OA マーカーが情報として残りました。

GBS と OA のジェノタイピングを合わせた結果、1005 個のマーカーが関連研究に含まれる結果となりました (GBS と OA からそれぞれ 903 個と 102 個)。イベリコ交雑種集団における遺伝子型別マーカーの最小対立遺伝子頻度 (MAF) 分布を図 2 に示します。これは、WGS データに由来する GBS マーカーパネル内の中間頻度を持つマーカーのかなりの割合を反映していますが、文献目録または RNAseq データによって提案された候補遺伝子に基づく OA パネル。

Genotyping by Sequencing (WGS データから取得) または OpenArray (文献目録および RNAseq データから取得) によって遺伝子型決定されたマーカーの MAF 分布。

関連研究は、遺伝子型特定に成功した1005のSNPに対して実行され、母集団内での分離と、利用可能なすべての成長、肥満、FA組成および代謝形質を示しました（補足表S4に含まれる記述統計）。有意な関連結果（q < 0.10）は補足表S5に含まれており、それらの結果のグラフ表示は図3に示されています。補足表S6には、少なくとも1つの有意な関連性を示す遺伝子について、有意な名目p値との関連結果が含まれています(q < 0.10)。

さまざまな形質について得られた有意な関連結果 (q < 0.10) の概要。形質の種類/段階 (x 軸) および多型の原因 (y 軸) ごとにグループ化されています。色付きのバーは、各遺伝子と有意に関連する形質の数を示します。 ADG の 1 日あたりの平均増体量、FCR 飼料変換率、WGS 全ゲノム配列。

成長形質と肥満形質に関しては、72 の有意な関連性が観察され、そのほとんどは、以前の生物学的または機能的仮説が利用可能な OA プラットフォーム (n = 59) で遺伝子型決定された候補遺伝子に関係していました。

発達の初期段階（誕生および離乳）における体重および形態学的測定に影響を与えるSNPを有する候補遺伝子には、ACACA、EPHX1、およびFREM2遺伝子が含まれます（補足表S5）。 ACACA 遺伝子内で分析された 2 つの SNP は、出生時に記録されたすべての身体測定値に体系的に影響を与えました (表 1)。離乳時には、EPHX1 と腹囲、および FREM2 と離乳時体重について時間通りの関連性が検出され、有意性閾値に達しました。両方の遺伝子の他の関連形質について、いくつかの名目レベルの関連性が観察されました（補足表S6）。 FREM2 は初期発達形質とのいくつかの示唆的な関連を示しましたが、EPHX1 は生涯の成長および発達形質との示唆的な関連を示しました。生後 150 日以降、体重、ADG、および FCR に影響を与える主な遺伝子は、NMNAT1、SOWAHB、EPHX1、TFRC、および SCD でした。それらのうち、SCD 遺伝子内で分析された 3 つの SNP は最も大きな影響を示し、たとえば、屠殺重量では 29 kg (3.2 SD)、枝肉重量では 10 kg (1.25 SD)、および 0.125 kg/d に達しました。 1.9 SD) グローバル ADG の場合 (表 2)。合計 13 の関連性には、WGS で検出された SNP が含まれ、コードされたタンパク質に対する潜在的な影響および成長グループ間の対立遺伝子頻度の違いに基づいて選択されました。 WGS データから得られた SNP に対応する結果は、LRIG3、DENND1B、TMCC1、MFRP など、出生時のみの体重と身体測定に関連する新しい遺伝子を提供します。出生形態的形質に関連する遺伝子の中で、ACACA (表 1) と TMCC1 が最も大きな影響を及ぼした遺伝子でした。

背脂肪は、LEPR、NFE2L2、FTO、およびLIPE候補遺伝子の影響を受けました（補足表S5）。 LEPR は、背脂肪に対して最も明白な影響を示し、成長段階におけるさまざまな脂肪層の深さおよびその他のいくつかの形質に対する有意または示唆的な影響 (表 3)、および屠殺時の最大規模の影響 (0.7 SD) を示しました。さらに、WGS データから得られた CFAP157/STXBP1 に位置する SNP は、屠殺時の背脂肪の厚さと関連しており、関連する大きさ (0.6 SD) の影響を及ぼしました。

脂肪酸組成は、背脂肪 (内層と外層)、および背最長筋と肝臓組織 (中性および極性脂質画分) で測定されました。これらの形質について合計 100 の有意な関連性が見つかりました (補足表 S5)。そのうち 45 の関連は背脂肪 FA 組成に対応し (13 は内層、32 は外層)、36 はロース FA に対応しました (25 は中性11は極性脂質画分）、19は肝臓FAプロファイルに対応しました（12は中性脂質画分、7は極性脂質画分）。成長形質と肥満形質に関しては、ほとんどの関連性が、WGS 由来のもの（n = 7）と比較して、OA によって遺伝子型決定された生物学的および機能的候補遺伝子にあるマーカーと一致しました（n = 93）。 100 の有意な関連性は、35 遺伝子に位置する 42 の SNP に関連しており、さまざまな形質または組織に対する 1 ～ 18 の有意な効果を示していました。少なくとも 2 つの有意な効果を持つ SNP の結果を表 4 に示します。FA 形質との有意な関連性が最も多かった遺伝子は KRT10 (表 4 および 5) で、18 件の有意な関連がありました。およびNLE1は、9つの有意な関連性を有し(表4および6)、両遺伝子は、主に短鎖飽和FAおよび一価不飽和FAに関して、背脂肪および筋肉中性脂質のFAプロファイルに影響を与える。 KCNH2 遺伝子の変異は、背脂肪外層においてのみいくつかの主要な PUFA 測定値 (C18:2n-6、C18:3n-3、PUFA、n-3、n-6) と有意に関連していた。一方、AHNAK、LIPE、EPHX1、TNF などの他の遺伝子は、同様の形質に影響を与える異なる場所に影響を及ぼしました (表 4)。

SCD 遺伝子に位置する 3 つの異なる SNP は、パルミトレイン酸および肝臓の中性脂質画分 (C16:1n-9、C18:1/C18:0) のみにおけるオレイン酸/ステアリン酸の不飽和比と関連していました (表 2、4)。ただし、内部背脂肪 FA 形質については名目レベルでの有意な効果が観察され、複数の検定で補正された有意性閾値には達しませんでした (補足表 S6)。また、脂肪生成の重要な遺伝子、FASN および ACACA は、背脂肪外層と筋肉の中性脂質の両方で、それぞれミリスチン酸 (C14:0) およびパルミトレイン酸 (C16:1n-7) と関連していました (表 4)。）。

さらに、IRX3、ADIPOQ、MSTN、MYOD1など、肥満関連形質に影響を与えることが知られている他の関連遺伝子は、1つの形質のみに重大な影響を及ぼしました（補足表S5）。

代謝形質に関連して、生後 110 日および 215 日の血漿中の HDL 濃度には 2 つの有意な関連性のみが見つかり、どちらも WGS マーカーに関連していました (補足表 S5)。 110 日齢で見つかった関連には TBC1D16 遺伝子に位置する SNP が関与しており、215 日齢で見つかった関連には遺伝子 CD4 が関与しています。

この研究に含まれる他の肉の品質特性（筋肉内脂肪含量、pH、およびドリップロス）については、複数のテスト補正後に有意な影響は観察されませんでした。

この研究では、関連する生産的形質に関連する DNA 多型を明らかにするために、最初の SNP 発見ステップの後に中スループットのジェノタイピングアプローチが実行されました。分析された SNP は主に WGS データから取得され、文献および RNAseq データから取得された生物学的および機能的な候補遺伝子/SNP のパネルで補完されました。

ジェノタイピングの結果は、WGS に由来するマーカーの良好なレベルの情報性を示し、高い割合のマーカーが中間頻度を示し (56% が 0.3 < MAF < 0.7 を示した)、低分離のために廃棄されたマーカーの割合はほとんどありません (10.6%)。これらの結果は、WGS データから実行される SNP 発見アプローチの有効性を裏付けています。逆に、文献目録および RNAseq データから得られたマーカーはあまり有益ではなく (33% が 0.3 < MAF < 0.7 を示した)、高い割合が単型またはほとんど有益ではない (50%) として検出されました。 WGS マーカーは関連研究に使用された同じイベリコ交雑種集団の動物からの配列データから得られたのに対し、候補遺伝子/SNP は異なるブタ品種および集団で以前に同定されていたため、この発見は予想外ではありませんでした。

WGS 由来のマーカー (102 対 903) と比較して、以前に知られていた候補 SNP の数が少なく、情報量が低いにもかかわらず、関連性または発現差の以前の証拠を持つ SNP について、はるかに多くの有意な関連性が観察されました。これらは、WGS データからの現在の研究で発見されたものよりも機能的な基準を満たしているため、これは驚くべきことではありません。実際、174 個の重要な関連性のうち、152 個は候補遺伝子の以前に同定された SNP に関係しており、WGS 由来のマーカーに関係しているのは 22 個だけでした。さらに、以前に特徴付けられたSNPに関連する152個の関連性のうち、62個は以前に関連性の証拠がある既知のマーカーに一致し、90個はイベリコ豚の遺伝子型間で差次的な発現を示す遺伝子内または潜在的な調節因子とのRNAseqデータから発見されたマーカーに一致した。役割。この発見は、表現型に興味深い影響を与える構造変異体を明らかにするためのディープシーケンシングのトランスクリプトームデータマイニングの有用性を裏付けています。

成長と発達の形質に関しては、ACACA と SCD という 2 つの主要な遺伝子が、それぞれ出生尺度と生涯のパフォーマンス形質に影響を与えることが判明しました。これらの遺伝子は脂肪の合成と代謝に関与しており、実際、FA 組成に対する両方の遺伝子の重大な影響も見つかりました。しかし、成長形質に対するそれらの潜在的な役割は今のところ明らかではありません。アセチル-CoA カルボキシラーゼ (ACACA) は、FA 合成の律速段階であるアセチル-CoA からマロニル-CoA へのカルボキシル化を触媒する多機能酵素系です。一致して、我々の結果は、背脂肪および筋肉中のパルミトレイン酸存在量に対するこの遺伝子の有意な影響と、多くの主要なFAおよびFA比に対する示唆的な影響を示しています。この主な役割に加えて、ACACA 遺伝子はミトコンドリア膜電位に影響を与えるため、細胞増殖と組織成長を調整する可能性があります 44。実際、ACACA 活性の欠損は、ヒトの初期成長と筋肉発達の欠陥に関連していました 45。ステアロイル CoA デサチュラーゼ (SCD) は、パルミトイル CoA とステアロイル CoA の Δ9-cis 不飽和化を触媒し、それぞれパルミトレオイル CoA とオレオイル CoA に変換されます。 MUFA は、組織脂質の成分であるだけでなく、シグナル伝達、細胞分化、食物摂取、アポトーシス、突然変異誘発のメディエーターとしても機能するため、哺乳動物における SCD 機能の変動は、さまざまな重要な機能に影響を与えると予想されます。成長を含む生理学的イベント46。したがって、BW の増加は SCD の活動の増加と関連している可能性があることが提案されています。実際、SCD 遺伝子の異なる多型が、痩せた豚の ADG および FCR と関連していました 47。それにもかかわらず、この研究で分析された 3 つの SNP (g.2108C>T; g.2228T>C; g.2281A>G) は、プロモーター領域に位置し、完全な連鎖不均衡にあり、これまではパフォーマンス特性と関連していませんでした 48,49 、50。この新たな発見は、特にBW形質に対する影響の大きさの点でも注目に値し、相加的影響は30kg近くと推定されている。ホモ接合遺伝子型の頻度は低く、効果が過大評価される可能性があるため、この発見は慎重に解釈される必要があります。興味深いことに、Estany ら 48 は、純血種デュロック系統において、CTA ハプロタイプが筋肉と皮下脂肪の両方で C18∶1/C18∶0 不飽和指数の増強と相加的に関連しているが、肝臓では関連していないことを発見しました。我々の結果によると、FA組成に関するSCD遺伝子の主な効果は肝臓の中性脂質画分で検出され、CTAハプロタイプは不飽和指数およびC16:1n-9含有量を増強します。また、この変異については IMF に対する潜在的な影響が報告されています 49 が、現在の研究では検証されていません。

成長形質の候補としてまだ確立されていない他の興味深い遺伝子および SNP も同定されました。それらの一部は RNAseq データから取得され、FREM2、EPHX1、SOWAHB、または TFRC5、13、14、15 など、イベリア豚とイベリア×デュロック交雑種またはデュロック純血種ブタの間で差次的に発現されました。したがって、これらの遺伝子は、イベリア遺伝子型とデュロック遺伝子型の間で観察される表現型の違いに関与するとみなされる機能的基準を満たしています。 FREM2 遺伝子 (FRAS1 関連細胞外マトリックスタンパク質 2) は、形態形成プロセスに役割を持つ細胞外マトリックスタンパク質をコードします 51。私たちの発見と一致して、FRAS/FREM複合体は開発の初期段階で特に重要です。ヒトでは、この遺伝子の変異は、皮膚上皮とその下にある間葉の構造的接着に関与する多系統疾患であるフレーザー症候群を引き起こし、ブタではこの遺伝子は、肉垂の存在などの発生形態学的形質の候補として提案されている53。 SOWAHB と EPHX1 は、いくぶん類似した関連結果を示し、どちらも生涯成長形質と一部の脂肪沈着形質に影響を与えました。遺伝子 SOWAHB (ソソンドワアンキリンリピートドメインファミリーメンバー B) はアンキリンファミリーの一部であり、膜骨格の組織化、イオン輸送、タンパク質認識、および細胞間接着の制御に関与しています 54。さらに、SOWAHB はインスリンに応答して活性化され、インスリン抵抗性に関与している可能性があります 55。成長形質に対するその効果に加えて、さまざまな年齢での肥満形質との関連性の検定（背脂肪の厚さ、コレステロール、補足表S6）で有意な名目p値が観察され、脂質代謝における役割を示唆しています。この遺伝子は、純血種イベリコ豚のロース肉で、交雑種の成長期豚よりも 4 倍高い発現を示しました 14。 EPHX1 (ミクロソームエポキシドヒドロラーゼ) 遺伝子は、成長および FA 組成形質、特に C20:1(n-9) および C20:3(n-6) 含量に関して有意な関連結果を示しました。この遺伝子によってコードされる酵素は、内因性多価不飽和FA由来のエポキシドを含む広範囲の基質を解毒または生物活性化することができ56、炎症や血管形成などのいくつかの生物学的プロセスを媒介し、細胞の恒常性、脂肪細胞の分化、およびインスリンの重要なシグナル伝達経路を調節します。応答57、58、59。この遺伝子は、デュロック遺伝子型と比較して、純粋なイベリア人の大腿二頭筋、背最長筋、背脂肪で上方制御されていることが示されました 13、14、15。純粋なイベリア人におけるさまざまな組織における SOWAHB および EPHX1 遺伝子の上方制御は、インスリン反応および肥満におけるそれらの潜在的な役割に関連している可能性があり、これは、それらの構造変異と肥満形質との示唆的な関連性および成長段階の体重との有意な関連性と一致します。イベリコ豚は、デュロック遺伝子型よりも代償的な成長とより激しい発育を発達させることが知られています。 TFRC 遺伝子 (トランスフェリン受容体) は、細胞の鉄取り込みに必要な細胞表面受容体をコードしており、疾患感受性の位置候補遺伝子として研究されています 60。さらに、JNK 経路の活性化を調節することでミトコンドリアの機能を調節する脂質センサーとしても機能します 61。興味深いことに、TFRCは出生時のデュロック種×イベリコ交雑種豚の腰筋で上方制御されていたが14、成長中のイベリコ純血種豚の大腿二頭筋と背脂肪では上方制御されていた5,15。私たちの結果によると、この遺伝子の SNP は BW と ADG に影響を与えます。したがって、出生期と成長期に観察される反対の規制は、イベリコとデュロックの遺伝子型で観察される筋肉発達プロセスの異なる不足に関連している可能性があり、後者はより高い出生前発育と出生時体重を示し、イベリコ純血種は成長期に筋肉発達の増加を示します。

新生児の発達形質は、WGS を通じて検出された、LRIG3、DENND1B、TMCC1 などの比較的未知の遺伝子に位置するいくつかの SNP によっても影響を受けました。これらの遺伝子はどれも生産的形質との関連について研究されていませんが、すべて初期の発生プロセスへの関与を示唆する役割または証拠を持っています。 LRIG3 (ロイシンリッチリピートおよび免疫グロブリン様ドメイン 3) は、頭蓋顔面の形態形成や神経堤の形成など、胚の発生に役割を果たしています 62。 KO マウスは野生型マウスよりも小さいため、LRIG3 と初期成長の直接的な関係も観察されました 63。一致して、我々の結果は、出生体重と胸郭周囲に重大な影響を及ぼし、頭部の発達に関連するものを含む新生児に記録される残りのすべての形質に示唆的な影響を与えることを示しています。 DENN Domain Containing 1B (DENND1B) 遺伝子は、サイトカイン産生と T 細胞受容体シグナル伝達の調節に重要な役割を果たしており、この因子の変異または欠損は新生児の免疫疾患に関連していました 64。興味深いことに、分娩時の静脈臍帯血中のDENND1B遺伝子のメチル化は、ヒトの出生体重と関連している65。我々の結果によると、この遺伝子は出生体重、腹囲、両頭頂径に影響を及ぼし、出生時と離乳時に記録された初期形質のほとんどに関連していることが示唆されています。 TMCC1 (膜貫通およびコイルドコイルドメインファミリー 1) は、中心脂肪蓄積およびウエスト周囲径と関連している内在性小胞体 (ER) 膜タンパク質であり 66、筋形成の調節に関与すると考えられています 67。我々の結果は、胸郭周囲だけでなく両頭頂部および胸郭の直径にも影響を及ぼし、また、他のいくつかの出生形質や背脂肪などのその後の肥満形質にも示唆的なレベルで影響を与えていることを示している。

肥育形質は主に、LEPR、LIPE、FTO などの既知の候補遺伝子にある SNP の影響を受け、LEPR で最も明確な結果が検出されました。レプチン受容体は、エネルギー恒常性に対する決定的な役割を持つことが知られており、脂肪沈着形質の機能的および位置的候補遺伝子として広く研究されています。イベリコ豚は、この遺伝子に固定変異を持っていることが知られており、これは、肥満と高い食欲、およびレプチン耐性表現型の特徴的な傾向に寄与する原因変異であると考えられています。この突然変異は、主に実験用ブタ集団だけでなく、他の品種や商業用ブタ集団など、さまざまな遺伝的背景で評価されています50、68、69、70。今回の結果によれば、このSNPの肥満に対する効果は、イベリコ×デュロックの商業集団においてさらに検証されており、さまざまな背脂肪の厚さの測定値と全体的なADGに大きな影響を及ぼし、BW、FA組成、コレステロール、HDL、LDLとの関連性が示唆されています。レベル。予想通り、イベリコ T 対立遺伝子は肥育を増加させ、背脂肪に対して 0.5 cm 近い相加効果がありました。これは、代替対立遺伝子のホモ接合動物間で推定約 1 cm の差があることを意味します (形質平均値の 20%)。

既知の候補遺伝子に加えて、NFE2L2 遺伝子 (核因子赤血球系 2 関連因子 2) も肥満と関連していました。この遺伝子は転写因子をコードしており、その主な役割は酸化損傷から保護する抗酸化タンパク質の発現の調節に関係します。さらに、この転写因子は恒常性の維持に決定因子であり、代謝、炎症、オートファジー、プロテオスタシス、ミトコンドリア生理学、免疫の調節に関与しています71。この遺伝子は、差次的に発現していないにもかかわらず、イベリコ豚の遺伝子型間で観察される差次的発現の予測制御因子としての潜在的な役割のため、SNP 発見のために RNAseq データから選択されました 14。我々の結果は、背脂肪の厚さに対する顕著な効果と、特にロースにおける筋肉内脂肪含量およびFA組成に対する示唆的な効果を示しており、これは筋肉のミトコンドリア生合成における既知の役割および筋肉代謝の調節因子としての予測される役割と一致している。

FA 組成に影響を与える主な遺伝子は KRT10 と NLE1 で、どちらも背脂肪とロース肉の SFA と MUFA の割合に影響を与えます。 KRT10 (ケラチン 10) は、交雑種豚と比較して純粋なイベリコ豚の筋肉で過剰発現されていたため研究されました 13 が、脂肪代謝に関連する役割はわかっていません。しかし、最近、KRT10 発現が血清脂質によって抑制されることが示され 72、興味深いことに、イベリコ豚とランドレース種の戻し交雑豚における 12 番染色体調査では、異なるケラチン遺伝子が FA 組成形質に関連していたことが示されました 73。実際、KRT10 遺伝子は SSC12 上にマッピングされており、FA 組成形質に関するいくつかの QTL が発見されています。さらに、KRT10は、ここでは分析されていない多型ALGA0066302Aに関するACACA遺伝子型に従って、最も有意に差次的に発現される遺伝子として見出され、これはFA組成と関連している。一方、NLE1 遺伝子 (Notchless ホモログ 1) は、脂肪細胞の分化と脂質代謝において中心的な役割を果たすことが知られている Wnt 経路の調節に役割を果たしています 74。実際、Wnt 経路は脂肪細胞において二重の機能を果たしており、これには、脂肪生成に対するよく知られた抑制効果と、成熟脂肪細胞におけるレプチン産生の刺激が含まれます 75。抗脂肪生成効果と一致して、純血種よりも痩せている交雑種イベリコ豚のロース肉でこの遺伝子の高い発現が観察されました 14。我々の関連性の結果は、この遺伝子を肥満形質と関連付けるものではないが、脂肪と筋肉という2つの主要な生産組織、特に最後の組織において、肉の品質に影響を与えるFA代謝とプロフィールに関連する関連形質にこの遺伝子が関与していることを明確に示している。オレイン酸含有量と一価不飽和脂肪から飽和脂肪の割合を増加させるこの対立遺伝子は、0.5～1SDの範囲の影響を及ぼし、成長後期のBWとADGに対して示唆的なプラスの効果を示しており、肉の品質を改善するための興味深いマーカーとなる可能性があることを意味している他の関連する形質にマイナスの多面発現効果を及ぼすことはありません。

他のいくつかの SNP は、少数の FA に限定して顕著な効果を示しました。 KCNH2 遺伝子 (カリウム電位依存性チャネルサブファミリー H メンバー 2) は背脂肪の主要な PUFA 含有量と関連していましたが、筋肉 FA プロファイルおよび一部の成長および枝肉形質に対して示唆的な効果が観察されました。この遺伝子は、交雑種やデュロック動物と比較して、純粋なイベリア人の腰と上腕二頭筋で過剰発現していることが以前に発見されました5,14。我々の結果は、KCNH2遺伝子の変異がインスリン恒常性およびグルコースおよび脂質代謝の変化に関連し76、そのメチル化が肥満に関連していることを示すヒトにおけるこれまでの知見と一致している。 AHNAK 遺伝子 (神経芽細胞分化関連タンパク質) は、脂肪生成と脂質代謝に関与する核タンパク質をコードします。 AHNAK KO マウスでは、脂肪の蓄積が減少し、血清トリグリセリドレベルが低下し、さらに脂肪分解と FA 酸化に関与する遺伝子の発現が増加しました 78。一致して、その発現は純粋なイベリコ動物でより高く14、交雑種よりも脂肪生成能が高いことを特徴としています。さらに、我々の結果は、背脂肪と肝臓のさまざまなFAに対するこのSNPの影響と、離乳時の背脂肪の厚さに対する示唆的な影響を示しています。

一方、遺伝子 FASN、LIPE、および TNF は、主に背脂肪と筋肉の C14:0 に対して同様の影響を及ぼしました。前述したように、LIPE 遺伝子（ホルモン感受性リパーゼ）は、屠殺時の背脂肪の深さに大きな影響を与える遺伝子の 1 つです。背脂肪にプラスの影響を与えるLIPE対立遺伝子は、背脂肪とロース肉のC14:0およびC16:0にもプラスの影響を与え、いくつかのMUFAパラメータおよびグルコース代謝指標に示唆的なマイナスの影響を及ぼしました。 LIPE は、アシルグリセロールからの FA の動員における重要な酵素であり 79、これは、種ごとに異なる脂肪沈着および FA プロファイル特性に関連しています。一方、FASN遺伝子（脂肪酸シンターゼ）はFAの合成を触媒する多機能酵素をコードしており、FA組成形質と繰り返し関連しており、我々の結果のように特に短鎖飽和FA21,80に影響を与えている。この遺伝子は、ACACA とともに、FA の新規合成を制御します。したがって、FASN、ACACA、およびLIPEは、脂肪生成と脂肪分解のバランスをとることによって脂肪沈着の調節に関与しており、これら3つすべてについて有意な関連結果が得られます。一方、腫瘍壊死因子（TNF）はインスリン抵抗性を促進するアディポカインであり、肥満誘発性の 2 型糖尿病と関連しています 81。一致して、我々の結果は背脂肪とロース肉のC14:0に対する有意な影響を示しているが、他の飽和FAや一価不飽和FA、枝肉や肉の品質形質（屠殺時の背脂肪の厚さ、ドリップロス）にも多くの示唆的な関連性が見出されている。そして成長に沿ったBWレコード。

代謝形質についてはほとんど有意な関連性が見出されなかった。 HDL 濃度のみが、2 つの異なる年齢における 2 つの異なる SNP と有意に関連しており、どちらも WGS データから得られました。生後 110 日で見つかった関連性は、膜輸送と分子輸送に役割を果たし、肥満では調節不全であることが知られている TBC1D16 遺伝子が関与しているため、興味深いものです 82。 HDL に対する有意な効果に加えて、この SNP は、腰部および肝臓の MUFA 含有量および MUFA/SFA 比に関連する多くの異なる FA 形質に対して示唆的な効果を示しました。これらは、体重や肥満指数などの肥満の物理的マーカーと相関していることが知られています 83。この遺伝子の役割と一致していますが、肥満には影響が観察されませんでした。

文献からの候補 SNP の回収と、WGS および RNAseq データからの SNP 発見を組み合わせたアプローチにより、これまで肥満、成長、FA 組成形質 (LEPR、SCD、ACACA、FASN など) に関連付けられていた候補遺伝子におけるさまざまな SNP の効果を検証することに成功しました。またはリペ）。さらに、私たちのアプローチにより、さまざまな年齢のBWとADG、および肥満とFA組成（NFE2L2、KRT10）に影響を与える、あまり知られていない遺伝子（LRIG3、DENND1B、NMNAT1、SOWAHB、EPHX1など）のSNPの興味深い新しい効果の発見が可能になりました。またはNLE1）。最も重要な関連性は、トランスクリプトームのディープシーケンシングデータの構造マイニングを通じて特定された SNP から得られました。 WGS 由来マーカーについて検出された有意な関連性は比較的少数であるにもかかわらず、ここで採用されたアプローチは、生産形質の未開発の遺伝的基盤を発見するための有用なゲノム規模の戦略として提案されています。私たちの結果は、成長と収量に関連する形質だけでなく、イベリコ豚の塩漬け製品に不可欠な肉の感覚的および技術的品質に関係する形質を含む、豚の関連形質の分子遺伝的基礎のより良い理解に貢献します。さらに、我々の結果は、異なるブタの品種および集団における関連研究および将来のマーカー支援選抜戦略の実施のための興味深いSNPおよび新しい候補遺伝子を提供する。 KRT10、NLE1、TFRC、NFE2L2、SCD、LEPR など、さまざまなマーカーが最も有望なものとして強調表示されます。これらのマーカーの中には、成長効率に悪影響を与えることなく肉の品質を向上させることができるものや、その逆の可能性があるものもあります。

著者らは、この研究の結果を裏付けるデータが論文およびその補足資料内で入手可能であることを確認しています。現在の研究中に生成された生の WGS データセットは、アクセッション番号 DRA014083 で DDBJ Sequence Read Archive (DRA) リポジトリに保管されています。

ロペス・ボーテ、CJ イベリコ豚品種の継続的な利用。肉科学。 49S1、S17–S27 (1998)。

記事 CAS Google Scholar

Rey, AI、Daza, A.、López Carrasco, C. & López-Bote, CJ 放し飼いまたは監禁状態でイベリコ豚にドングリと牧草を与えると、枝肉の特性、背最長筋と背脂肪の脂肪酸とトコフェロールの蓄積に影響を与える。肉科学。 73、66–74 (2006)。

記事 CAS Google Scholar

Juárez, M.、Clemente, I.、Polvillo, O.、Molina, A. 品種系統と交配による影響を受けたイベリコ豚のテンダーロインの肉質。肉科学。 81、573–579。 https://doi.org/10.1016/j.meatsci.2008.10.016 (2009)。

論文 PubMed Google Scholar

クレメンテ、I.ら。純イベリコ豚と交配イベリコ豚の豚肉品質の比較。第7回地中海の豚に関する国際シンポジウムにて。地中海オプション: シリーズ A 地中海セミナー、n. 101 (ペドロ、EJ およびカベサス、AB 編) 441–446 (CIHEAM、2012)。

Google スカラー

ベニテス、R. 他異なるエネルギー源を与えられたイベリコ豚の成長に伴う大腿二頭筋のトランスクリプトームの変化とデュロック種との比較分析。動物 11、3505 (2021)。

記事 Google Scholar

バスケス・ゴメス、M. 他子豚の出生体重と性別は、脂肪の多い豚の発育能力と脂肪酸組成に影響します。アニム。製品。科学。 60、573–583 (2020)。

記事 Google Scholar

オヴィロ、C.ら。ブタ第 6 染色体 QTL と LEPR の詳細なマッピングが、複数世代のイベリア人とランドレース種の交配における体組成に及ぼす影響。ジュネ研究所。 85、57–67 (2005)。

記事 CAS Google Scholar

オヴィロ、C. et al. ブタレプチン受容体 (LEPR)、神経ペプチド Y (NPY)、およびコカインおよびアンフェタミン調節転写 (CART) 遺伝子の視床下部発現は、LEPR 遺伝子型の影響を受けます。マム。ゲノム。 21、583–591 (2010)。

記事 Google Scholar

コロミナス、J. et al. ELOVL6 遺伝子の多型は、ブタの脂肪酸組成に対する主要な QTL 効果と関連しています。 PLoS One 8、e53687。 https://doi.org/10.1371/journal.pone.0053687 (2013)。

論文 ADS CAS PubMed PubMed Central Google Scholar

コロミナス、J. et al. 豚の肉品質形質の候補遺伝子としての豚ACSL4遺伝子の評価。アニム。ジュネット。 43、714–720。 https://doi.org/10.1111/j.1365-2052.2012.02335.x (2012)。

論文 CAS PubMed Google Scholar

メルカデ、A.ら。ブタのアシル-CoA シンテターゼ長鎖 4 遺伝子の特徴付けと、成長および肉の品質形質との関連。アニム。ジュネット。 37、219–224。 https://doi.org/10.1111/j.1365-2052.2006.01436.x (2006)。

論文 CAS PubMed Google Scholar

フェルナンデスバローゾ、M. Á. 他。野外放し飼いイベリコ豚の肉質特性に関する遺伝パラメーター推定と遺伝子関連分析。Ｊ．アニム．繁殖。ジュネット。 137、581–598。 https://doi.org/10.1111/jbg.12498 (2020)。

論文 CAS PubMed Google Scholar

Ayuso, M. et al. ブタの遺伝子型間の筋肉トランスクリプトームの比較分析により、成長、肥満、代謝に関連する遺伝子と制御機構が特定されます。 PLoS One 10、e0145162。 https://doi.org/10.1371/journal.pone.0145162 (2015)。

論文 CAS PubMed PubMed Central Google Scholar

Ayuso, M. et al. 発育段階、筋肉、および遺伝子型によりブタの筋肉トランスクリプトームが変化する: 成長と代謝に関与する遺伝子発現および調節因子への影響。 PLoS One 11、e0167858。 https://doi.org/10.1371/journal.pone.0167858 (2016)。

論文 CAS PubMed PubMed Central Google Scholar

ベニテス、R. 他異なるエネルギー源を与えられたイベリコ豚とデュロック豚の脂肪組織トランスクリプトームに対する品種、食事、および相互作用の影響。 Genes 10、589 (2019)。

記事 Google Scholar

ヴァン・ソン、M.ら。ゲノムワイド関連研究により、デュロック豚の SSC14 における背脂肪脂肪酸組成の主要な QTL が確認されました。 BMCジェノム。 18、369 (2017)。

記事 Google Scholar

Zhang、W.ら。 5 つの異なるブタ集団における脂肪酸代謝形質のゲノムワイド関連研究。科学。議員第 6 号、24718 (2016)。

記事 ADS CAS Google Scholar

ラマヨ・カルダス、Y.ら。イベリア×在来種交雑種における筋肉内脂肪酸組成のゲノムワイド関連研究。Ｊ．アニム．科学。 90、2883–2893 (2012)。

記事 CAS Google Scholar

ムニョス、M. 他高密度の遺伝子型解析および発現データを使用した、豚の背脂肪および筋肉内脂肪の脂肪酸組成のゲノムワイド分析。 BMCジェノム。 14、845 (2013)。

記事 Google Scholar

マルティネス・モンテス、Á. 金属。ゲノムワイド関連研究を使用して、イベリコ豚品種に基づく 3 つの異なる遺伝的背景における共通の QTL 領域を特定します。 PLoS One 13、e0190184。 https://doi.org/10.1371/journal.pone.0190184 (2018)。

論文 CAS PubMed PubMed Central Google Scholar

Crespo-Piazuelo、D. et al. イベリコ豚に基づく 3 つの交配における背脂肪および筋肉内脂肪酸組成の強力な候補遺伝子の同定。科学。議員 10、13962 (2020)。

記事 CAS Google Scholar

ペナ、RN et al. 5 つのゲノム領域がイベリコ豚の脂肪組成に大きな影響を与えます。科学。議員9、2031年。https://doi.org/10.1038/s41598-019-38622-7 (2019)。

論文 ADS CAS PubMed PubMed Central Google Scholar

ムニョス、M. 他高密度 SNP チップを使用して評価された、ヨーロッパの在来豚品種におけるゲノムの多様性、連鎖不均衡、および選択サイン。科学。議員9、13546。https://doi.org/10.1038/s41598-019-49830-62019 (2019)。

論文 ADS PubMed PubMed Central Google Scholar

Daetwyler、HD et al. 234 頭の雄牛の全ゲノム配列決定により、牛の単一遺伝子形質と複雑形質のマッピングが容易になります。ナット。ジュネット。 46、858–865 (2014)。

記事 CAS Google Scholar

de Blas, C.、Gasa, J. & Mateos, GG 豚の栄養要件。 FEDNA規格。スペインの動物栄養開発財団。 http://fundacionfedna.org/sites/default/files/Normas%20PORCINO_2013rev2_0.pdf (2013)。

Calvo, L.、Toldrá, F.、Aristoy, MC、Lopez-Bote, CJ & Rey, AI 豚肉の筋肉タンパク質分解活性と保水能力に対する食事性有機セレンの影響。肉科学。 121、1–11 (2016)。

記事 CAS Google Scholar

Segura, J. & Lopez-Bote, CJ 筋肉内脂肪分析のための実験室での効率的な方法。食品化学。 145、821–825 (2014)。

記事 CAS Google Scholar

Hulbert, AJ、Pamplona, R.、Buffenstein, R. & Buttemer, WA 生と死: 動物の代謝率、膜組成、寿命。生理。改訂 87、1175–1213 (2007)。

記事 CAS Google Scholar

ハルバー、MW et al. 骨格筋におけるステアロイル CoA デサチュラーゼ 1 発現の上昇は、肥満ヒトにおける異常な脂肪酸分配の一因となります。細胞メタブ。 2、251–261 (2005)。

記事 CAS Google Scholar

Bolger, AM、Lohse, M. & Usadel, B. Trimmomatic: Illumina シーケンスデータ用の柔軟なトリマー。バイオインフォマティクス 30、2114–2120 (2014)。

記事 CAS Google Scholar

ウォー、A.ら。ブタの遺伝学およびゲノム研究を可能にする改良されたブタ参照ゲノム配列。ギガサイエンス 9、giaa051 (2020)。

記事 Google Scholar

Li, H. & Durbin, R. Burrows-Wheeler 変換による高速かつ正確なショートリードアライメント。バイオインフォマティクス 25、1754–1760 (2009)。

記事 CAS Google Scholar

リー、Ｈら。配列アライメント/マップ形式と SAMtools。バイオインフォマティクス 25、2078–2079。 https://doi.org/10.1093/bioinformatics/btp352 (2009)。

論文 CAS PubMed PubMed Central Google Scholar

ポプリン、R. et al. 正確な遺伝的変異の発見を数万のサンプルに拡張します。バイオRxiv。 https://doi.org/10.1101/201178 (2018)。

記事 Google Scholar

Van der Auwera、GA & O'Connor、BD Genomics in the Cloud: using Docker、GATK、WDL in Terra 1st edn. （オライリーメディア、2020）。

Google スカラー

Danecek、P. et al. 1000ゲノムプロジェクト分析グループ。バリアント呼び出しフォーマットと VCFtools。バイオインフォマティクス 27、2156–2158。 https://doi.org/10.1093/bioinformatics/btr330 (2011)。

論文 CAS PubMed PubMed Central Google Scholar

マクラーレン、W. 他 Ensembl バリアント効果予測子。ゲノムバイオル。 17、122。https://doi.org/10.1186/s13059-016-0974-4 (2016)。

論文 CAS PubMed PubMed Central Google Scholar

Kumar、P.ら。 SIFT アルゴリズムを使用した、タンパク質機能に対する非同義バリアントのコーディングの影響の予測。ナット。プロトック。 4、1073–1081。 https://doi.org/10.1038/nprot.2009.86 (2009)。

論文 CAS PubMed Google Scholar

Chang, C. et al. 第 2 世代 PLINK: 大規模で豊富なデータセットの課題に対処します。ギガサイエンス 4, 559 (2015)。

記事 Google Scholar

Ayuso, M. et al. 純粋なイベリコ豚と交配したイベリコ豚の間のトランスクリプトームおよび表現型の違いに潜在的に関与する SNP の同定と優先順位付け。第 36 回国際動物遺伝学会 (ISAG) 会議、抄録集、p. 128 (2017)。

ヤンら。 GCTA: ゲノム全体にわたる複雑な形質解析のためのツール。午前。 J. ハム。ジュネット。 88、76–82 (2011)。

記事 CAS Google Scholar

Storyy, JD, Bass, AJ, Dabney, A. & Robinson, D. Qvalue: 誤検出率制御のための Q 値推定。 R パッケージバージョン 2.16.0。 http://github.com/jdstorey/qvalue (2019)。

RStudioチーム。 RStudio: R の統合開発。http://www.rstudio.com/ (RStudio、PBC、2020)。

Zhang、H.ら。 ACACAの下方制御は、ミトコンドリア電位の阻害を通じて前立腺がんの悪性進行を抑制します。 J. Cancer 12、232–243。 https://doi.org/10.7150/jca.49560 (2021)。

論文 CAS PubMed PubMed Central Google Scholar

Blom, W.、de MuinckKeizer、SMPF & Stolte、HR アセチル-CoA カルボキシラーゼ欠損症: 新しい脂肪酸合成の先天的エラー。Ｎ．Ｅｎｇｌ． J.Med. 305、465–466 (1981)。

記事 CAS Google Scholar

Ntambi, JM ステアロイル CoA デサチュラーゼ 1 は、エネルギー恒常性の生物学的調節因子です。脂質代謝におけるステアロイル CoA デサチュラーゼ遺伝子 (Ntambi, J. 編) (Springer、2013)。 https://doi.org/10.1007/978-1-4614-7969-7_3。

Google Scholar の章

Bartz, M.、Kociucka, B.、Salamon, S.、Jelen, H. & Switonski, M. ブタ SCD 遺伝子の転写レベルは筋肉および脂肪組織の脂肪酸組成に影響を及ぼさないが、推定上の脂肪酸組成には影響を与えない。 miRNA の標的部位は、毎日の増量および飼料変換率と関連しています。Ｊ．アニム．科学。 https://doi.org/10.2527/jas.2012-5380 (2012)。

論文 PubMed Google Scholar

Estany, J.、Ros-Freixedes, R.、Tor, M. & Pena, RN ステアロイル CoA デサチュラーゼ遺伝子プロモーターの機能的変異体は、豚肉の脂肪酸不飽和化を促進します。 PLoS One 9、e86177。 https://doi.org/10.1371/journal.pone.0086177 (2014)。

論文 ADS CAS PubMed PubMed Central Google Scholar

フェルナンデス、AI 他異なる遺伝的背景における脂肪酸不飽和化および脂肪沈着に対するブタ SCD ハプロタイプの影響を検証します。ライブベスト。科学。 205、98–105。 https://doi.org/10.1016/j.livsci.2017.09.021 (2017)。

記事 Google Scholar

Ros-Freixedes、R. et al. ゲノムワイド関連研究では、デュロック豚の筋肉内脂肪含量と脂肪酸組成に影響を与える 2 つの主要な遺伝子座として SCD と LEPR が挙げられています。 PLoS One 11、e0152496。 https://doi.org/10.1371/journal.pone.0152496 (2016)。

論文 CAS PubMed PubMed Central Google Scholar

ティマー、JR 他組織の形態形成と血管の安定性には、マウスの髄脳小胞遺伝子の産物である Frem2 タンパク質が必要です。手順国立アカド。科学。 USA 102、11746–11750 (2005)。

記事 ADS CAS Google Scholar

Petrou, PP、Makrygiannis, AK & Chalepakis, G. 細胞外マトリックスタンパク質の Fras1/Frem ファミリー: 構造、機能、およびフレーザー症候群およびマウスブレブ表現型との関連。接続する。組織検査 49、277–282 (2008)。

記事 CAS Google Scholar

Schiavo, G. et al. 土着のカセルタナ豚品種におけるいくつかの外部形質のゲノムワイド関連解析。ライブベスト。科学。 230、103842。https://doi.org/10.1016/j.livsci.2019.103842 (2019)。

記事 Google Scholar

Hryniewicz-Jankowska, A.、Czogalla, A.、Bok, E. & Sikorsk、AF アンキリン、多くの細胞経路に関与する多機能タンパク質。フォリア・ヒストケム。サイトビオール。 40、239–249 (2002)。

CAS PubMed Google Scholar

Haider, N.、Mann, M. & Kahn, CR 非糖尿病および糖尿病の個人におけるインスリン抵抗性および性別による変化に関連するシグナル伝達の欠陥。Ｊ．クリン．調査します。 131、e151818。 https://doi.org/10.1172/JCI151818 (2021)。

論文 CAS PubMed PubMed Central Google Scholar

Václavíková, R.、Hughes、DJ & Souček、P. ミクロソームエポキシド加水分解酵素 1 (EPHX1): ヒトの疾患における遺伝子、構造、機能、および役割。遺伝子 571、1 ～ 8。 https://doi.org/10.1016/j.gene.2015.07.071 (2015)。

論文 CAS PubMed PubMed Central Google Scholar

McReynolds, C.、Morisseau, C.、Wagner, K. & Hammock, B. エポキシ脂肪酸は、痛み、心血管疾患、およびミトコンドリア機能不全、小胞体ストレス、炎症を特徴とするその他の適応症の治療の有望な標的です。上級経験値医学。バイオル。 1274 年、71 ～ 99 年。 https://doi.org/10.1007/978-3-030-50621-6_5 (2020)。

論文 CAS PubMed PubMed Central Google Scholar

Gautheron、J. et al. EPHX1 変異は、エポキシドの加水分解障害と細胞老化の増加により、脂肪萎縮性糖尿病症候群を引き起こします。 Elife 10、e68445。 https://doi.org/10.7554/eLife.68445 (2021)。

論文 CAS PubMed PubMed Central Google Scholar

エディン、ML 他。エポキシド加水分解酵素 1 (EPHX1) はエポキシエイコサノイドを加水分解し、虚血後の心臓の回復を阻害します。Ｊ．Ｂｉｏｌ．化学。 293、3281–3292。 https://doi.org/10.1074/jbc.RA117.000298 (2018)。

論文 CAS PubMed PubMed Central Google Scholar

ヤコブセン、M.ら。ブタの ETEC F4ab/ac 候補領域内の 5 つの候補遺伝子の特性評価。 BMC 解像度注記 4、225 (2011)。

記事 CAS Google Scholar

Senyilmaz, D. et al. TFR1のステアロイル化によるミトコンドリアの形態と機能の調節。ネイチャー 525、124–128 (2015)。

記事 ADS CAS Google Scholar

アブレイラ、VE 他。 LRIG3 とネトリン 1 の間の交差抑制相互作用は、内耳の構造を形成します。開発 135、4091–4099。 https://doi.org/10.1242/dev.029330 (2008)。

論文 CAS PubMed Google Scholar

ヘルストロム、M. et al. Lrig3欠損マウスにおける心臓肥大と高密度リポタンパク質コレステロールの減少。午前。Ｊ．Ｐｈｙｓｉｏｌ．レギュラー。統合します。コンプ。生理。 310(11)、R1045–R1052 (2016)。

記事 Google Scholar

ヤン、CW et al. TH2 細胞およびアレルギー疾患における DENND1B による T 細胞受容体シグナル伝達の制御。セル 164、141 ～ 155。 https://doi.org/10.1016/j.cell.2015.11.052 (2016)。

論文 CAS PubMed Google Scholar

アガ、G.ら。在胎期間当たりの出生体重は、出生時および小児期の DNA メチル化と関連しています。クリン。エピジェネット。 8、118。 https://doi.org/10.1186/s13148-016-0285-3 (2016)。

記事 CAS Google Scholar

米山伸ほかヨーロッパ系の最大57,412人の中枢性肥満形質に関する遺伝子中心のメタ分析により、既知の遺伝子座が確認され、いくつかの新規関連性が明らかになった。ハム。モル。ジュネット。 23、2498–2510。 https://doi.org/10.1093/hmg/ddt626 (2014)。

論文 CAS PubMed Google Scholar

原田明ほか筋原性調節配列の空間的再構成は、遺伝子発現を時間的に制御します。核酸研究所 https://doi.org/10.1093/nar/gkv046 (2015)。

記事 PubMed PubMed Central Google Scholar

ペレス・モンタレロ、D. 他生産性と品質形質に対するブタレプチンとレプチン受容体多型の共同効果。アニム。ジュネット。 43、805–880 (2012)。

記事 Google Scholar

Galve, A. et al. デュロック×ランドレース/ラージホワイト交配における脂肪含量、脂肪分布、脂肪組成に対するレプチン受容体（LEPR）およびメラノコルチン-4受容体（MC4R）多型の影響。ライブベスト。科学。 145、145–152 (2012)。

記事 Google Scholar

植本裕也ほかデュロック純血種集団における背脂肪の厚さ、画像解析による脂肪面積比、および血清レプチン濃度に対するブタレプチン受容体遺伝子多型の影響。アニム。科学。 J. 83、375–385 (2012)。

記事 CAS Google Scholar

He、F.、Ru、X.、Wen、T. NRF2、ストレス応答以降の転写因子。内部。Ｊ．Ｍｏｌ．科学。 21、4777。https://doi.org/10.3390/ijms21134777 (2020)。

論文 CAS PubMed Central Google Scholar

Aldehlawi、H. et al. 血清脂質、レチノイン酸、およびフェノールレッドは、培養ケラチノサイトにおけるケラチン K1、K10、および K2 の発現を示差的に制御します。科学。議員 10、4829。https://doi.org/10.1038/s41598-020-61640-9 (2020)。

論文 ADS CAS PubMed PubMed Central Google Scholar

ムニョス、M. 他高密度SNPデータを用いた筋肉内脂肪の脂肪酸組成に影響を与えるSSC12領域の調査。フロント。ジュネット。 https://doi.org/10.3389/fgene.2011.00101 (2012)。

記事 PubMed PubMed Central Google Scholar

カプート、T.ら。白色脂肪組織における肥満関連炎症の開始時の抗脂肪生成シグナル。細胞。モル。生命科学。 78、227–247。 https://doi.org/10.1007/s00018-020-03485-z (2021)。

論文 CAS PubMed Google Scholar

チェン、ZLら。急性 Wnt 経路の活性化は、成熟脂肪細胞におけるレプチン遺伝子発現を正に制御します。セル信号。 27、587–597。 https://doi.org/10.1016/j.cellsig.2014.12.012 (2015)。

論文 CAS PubMed Google Scholar

Engelbrechtsen、L. et al. hERG (KCNH2) 電位依存性カリウムチャネルの一般的な変異は、空腹時およびグルコース刺激による血漿インクレチンおよびグルカゴン反応の変化に関連しています。 BMCジュネット。 19、15。https://doi.org/10.1186/s12863-018-0602-2 (2018)。

論文 CAS PubMed PubMed Central Google Scholar

ゾンネ、SB et al. 肥満は、脂肪細胞の DNA メチル化および遺伝子発現におけるデポー特異的な変化と関連しています。脂肪細胞 6、124–133。 https://doi.org/10.1080/21623945.2017.1320002 (2017)。

論文 CAS PubMed PubMed Central Google Scholar

キム、YN 他アーナク欠乏症は、FGF21 誘導を通じてマウスの高脂肪食誘発性脂肪肝を軽減します。経験値モル。医学。 53、468–482 (2021)。

記事 CAS Google Scholar

Lampidonis, AD、Rogdakis, E.、Voutsinas, GE & Stravopodis, DJ 哺乳類の脂肪分解におけるホルモン感受性リパーゼ (HSL) の復活。遺伝子 477、1–11。 https://doi.org/10.1016/j.gene.2011.01.007 (2011)。

論文 CAS PubMed Google Scholar

ザパテラ、M.ら。背脂肪脂肪酸組成と候補遺伝子のSNPとの関連研究により、大型白豚におけるFASN多型の影響が浮き彫りになった。肉科学。 156、75–84。 https://doi.org/10.1016/j.meatsci.2019.05.013 (2019)。

論文 CAS PubMed Google Scholar

Sethi、JK および Hotamisligil、GS 代謝メッセンジャー: 腫瘍壊死因子。ナット。メタブ。 3、1302–1312。 https://doi.org/10.1038/s42255-021-00470-z (2021)。

論文 CAS PubMed Google Scholar

ヤコブセン、MJ 他肥満ブタからの純粋な脂肪細胞画分のエピジェネティックおよびトランスクリプトーム特性評価により、代謝を制御する候補経路が特定されます。フロント。ジュネット。 17、1268 (2019)。

記事 Google Scholar

ジェヤクマール、SM 他脂肪酸不飽和指数は、Wistar NIN 肥満変異体ラット系統 WNIN/Ob および WNIN/GR-Ob における体重および肥満の肥満指数と相関します。ニュートル。メタブ。 (ロンドン) 11、27。https://doi.org/10.1186/1743-7075-6-27 (2009)。

記事 CAS Google Scholar

リファレンスをダウンロードする

この実験研究は、Agencia Estatal de Investigacion (プロジェクト AGL2013-48121-C3-R; AGL2016-79321-C2-1-R; PID2019-108695RB-C31/AEI/10.13039/501100011033) から資金提供を受けました。 MVG と CGC はスペイン政府の支援を受けました (MVG: FPU National Program Grant Number FPU014/01285. CGC: FPI National Program Grant Number BES-2014-070464)。動物素材の提供にご協力いただいたイベリコス・デ・アラウゾSLのスタッフに感謝いたします。

動物遺伝改善部門、INIA-CSIC、マドリード、スペイン

C. オビロ、Y. ヌニェス、F. ガルシア、JM ガルシア＝カスコ、M. ムニョス

家畜生物学研究所、FBN、ドゥマーシュトルフ、ドイツ

N. トラクールジュル、F. ハドリッヒ、E. ムラニ & K. ウィマーズ

CoPeD、獣医学部、アントワープ大学、ヴィルレイク、ベルギー

アユソさん

スペイン、グラナダのCSIC、栄養および持続可能な動物生産局

C.ガルシア＝コントレラス

ソルボンヌ大学、INSERM、NutiOmics、パリ、フランス

M・バスケス＝ゴメス

UCM獣医学部動物生産学科、マドリード、スペイン

AI レイ、C. ロペスボーテ、B. イザベル

カルデナル・エレーラ・CEU大学獣医学部、バレンシア、スペイン

A. ゴンザレス・ブルネス

PubMed Google Scholar でこの著者を検索することもできます

CO、AGB、BIR、AIR、CLB がこの研究を考案、設計し、資金を獲得しました。 YN、CGC、FG はウェットラボ作業を行いました。 CO、NT、FH、MM、MA は、バイオインフォマティクスと統計解析に貢献しました。 BIR、AIR、AGB、MVG、CGC がサンプルを提供しました。 CO、MM、MVG、JGC、CLB、AGB、EM、KW がデータ解釈に貢献しました。 CO が論文を書きました。すべての著者が提出されたバージョンを読んで承認しました。

C. オヴィロへの対応。

著者らは競合する利害関係を宣言していません。

シュプリンガーネイチャーは、発行された地図および所属機関における管轄権の主張に関して中立を保ちます。

オープンアクセスこの記事はクリエイティブコモンズ表示 4.0 国際ライセンスに基づいてライセンスされており、元の著者と情報源に適切なクレジットを表示する限り、あらゆる媒体または形式での使用、共有、翻案、配布、複製が許可されます。クリエイティブコモンズライセンスへのリンクを提供し、変更が加えられたかどうかを示します。この記事内の画像またはその他のサードパーティ素材は、素材のクレジットラインに別段の記載がない限り、記事のクリエイティブコモンズライセンスに含まれています。素材が記事のクリエイティブコモンズライセンスに含まれておらず、意図した使用が法的規制で許可されていない場合、または許可されている使用を超えている場合は、著作権所有者から直接許可を得る必要があります。このライセンスのコピーを表示するには、http://creativecommons.org/licenses/by/4.0/ にアクセスしてください。

転載と許可

オヴィロ、C.、トラクールジュル、N.、ヌニェス、Y. 他イベリコ交雑種豚の成長、肥満、肉質特性に関する SNP の発見と関連研究。 Sci Rep 12、16361 (2022)。 https://doi.org/10.1038/s41598-022-20817-0

引用をダウンロード

受信日: 2022 年 5 月 4 日

受理日: 2022 年 9 月 19 日

公開日: 2022 年 9 月 30 日

DOI: https://doi.org/10.1038/s41598-022-20817-0

次のリンクを共有すると、誰でもこのコンテンツを読むことができます。

申し訳ございませんが、現在この記事の共有リンクは利用できません。

Springer Nature SharedIt コンテンツ共有イニシアチブによって提供

コメントを送信すると、利用規約とコミュニティガイドラインに従うことに同意したことになります。虐待的なもの、または当社の規約やガイドラインに準拠していないものを見つけた場合は、不適切としてフラグを立ててください。