クローニングのための遺伝学
(前編)
Akifumi Shimizu1. 遺伝のしくみ
○ 遺伝子(gene、DNA 配列)は、染色体(chromosome)に乗って遺伝する ○ 遺伝子が染色体に座乗する位置は決まっている(遺伝子座 locus) ○ 同一遺伝子座に座乗できる遺伝子(対立遺伝子allele)は、複数種が存在しうる ○ 染色体は、減数分裂時に乗り換える(crossing over)ことがある ○ 染色体が乗り換えることで、染色体上の対立遺伝子の組み合わせは変化する 一対の染色体からなる架空の2倍体生物の場合の遺伝の様子をみてみましょう(図 1)。こ の生物の体細胞は2n=2 本の染色体を持ちます。子孫へと遺伝情報を伝えるのは、卵と精子 という特別な細胞(配偶子)のみです。配偶子は減数分裂をして n=1 本の染色体を持っていま す(半数体、一倍体と呼ぶ)。図 1 ではアルファベットで遺伝子座の種類を、下付き数字で対 立遺伝子の種類を区別しています。 減数分裂時には相同染色体の姉妹染色体分体間で乗り換え(Crossing over)が、少なくとも 一回起こります。乗り換えにより、遺伝子座の位置は変わりませんが、同一染色体上の対 立遺伝子の組合せが変わります。生じる組合せの頻度は遺伝子座間の位置関係(遺伝距離) と関係します。図1 の例では、4つの配偶子のうち 2 つで対立遺伝子の組合せが変わって います(新たに A1B1C1D1E1F1G2H2とA2B2C2D2E2F2G1H1が生じた)。 実際には配偶子は一度に何種類もつくられます。沢山の子孫の対立遺伝子セットを調べ ることで、遺伝子が一緒に伝わる(連鎖する)か、伝わらない(独立する)かの頻度を調べその 比率をできるようになります。この比率は、遺伝子の染色体上の座乗位置と関係します。 イネなど二倍体の自殖性植物の場合、両親を交配して得た交雑当代(F1)を自殖した F2集図1 減数分裂のモデル図(2n=2)
○はセントロメア, セントロメアを中心として短い方を短腕(A~C)、長い方(D~H)を長腕という a)普通の細胞 b)染色体の倍加 c)乗り換え d)4個の配偶子へ 2n = 2 (減数分裂中) (減数分裂中) n = 1 P1 P2 P1 P1 P2 P2 P1 P1 P2 P2 A1 A2 A1 A1 A2 A2 B1 B2 B1 B1 B2 B2 C1 C2 C1 C1 C2 C2 ○ ○→
○ ○→
○ ○→
○ ○ ○ ○ D1 D2 D1 D1 D2 D2 E1 E2 E1 E1 E2 E2 F1 F2 F1 F1 F2 F2 G1 G2 G1 G2 G1 G2 H1 H2 H1 H2 H1 H22. DNA マーカー
染色体の乗換えにより遺伝情報がどのように継承されるかを追跡するための便利なツー ルに、DNA マーカーがあります。手軽に使える DNA マーカーの1つに SSR マーカーがあ ります。SSR マーカーは、”GAGAGAGA”のように単純な塩基単位の繰り返し(Simple Sequence Repeats, SSR)によって構成される DNA 配列を元にしたマーカーで、SSR は真 核生物のゲノム中に普遍的に散在していることが知られています。SSR の反復数は変異し 易いことが知られており、反復数の変異を断片長多型として検出することで複数の対立遺 伝子を識別することができます。SSR マーカーは、電気泳動法によって断片長多型を識別 するのが一般的です(図2)。 図2 電気泳動法による SSR 断片長多型の検出例 或るSSR マーカーについて断片長が異なる 2 倍体の自殖品種 A と B を考えます。品種 A のSSR 断片長を A 型とし、品種 B のものを B 型としたとき、品種 A と品種 B を交配した F A 型 B 型両方の断片長を持つヘテロ遺伝子型になるはずです。図 3 は、品種 A の (GA)6 6反復 (GA)10 10反復 GAGAGAGAGAGA GAGAGAGAGAGAGAGAGAGA CTCTCTCTCTCT CTCTCTCTCTCTCTCTCTCT PCRで増幅 PCRで増幅 電気泳動による断片長多型の検出
図3 SSR マーカーの断片長多型と遺伝子型の実例
他のDNA マーカーとしては、SNP(Single Nucleotides Polymorphisms、一塩基多型)マ ーカーがあります。SNP マーカーは最小単位の DNA 多型であり、もっとも高密度な多型 情報を得られます。
3. マーカーの連鎖
例としてイネの染色体4 に座乗する SSR マーカー(RM252 と RM241)の連鎖を見てみま しょう。品種G×品種 K の自殖 F2 86 系統について遺伝子型を調べたところ、表 1 のよう になりました。各マーカーの遺伝子型は、品種G ホモ型(GG), 品種 K ホモ型(KK)、ヘテロ 型(GK)の 3 種類で、2 種の組み合わせは 3 種類×3 種類の 9 種類になります(表 1)。表1 この 2 マーカー間の組換え価を計算してみます。 マーカー間の組換え価がrであると仮定すると、RM252 と RM241 について、品種 G と 品種K を交配した雑種当代 F1個体の遺伝構成は次のようになるはずです。 この F1を自殖して F2分離集団をつくります。このとき、F1からできる配偶子は RM252 とRM241 について 4 種類に分類でき、 の頻度で非組換え型と組換え型の配偶子が生じます。卵配偶子(4 種類)と花粉配偶子(4 種類) RM252 RM241 86 GG GG 24 GG GK 4 GG KK 0 GK GG 4 GK GK 31 GK KK 3 KK GG 0 KK GK 7 KK KK 13 F1 組換え価
r
RM252 G K | |r
RM241 G K 非組換え型配偶子 組換え型配偶子頻度 (1-
r
)/2
r
/2
G
K
G
K
|
|
|
|
G
K
K
G
①
②
③
④
それぞれ、配偶子の頻度に応じて9 種類の遺伝子型クラスの期待頻度が記されています(自 分で紙に書き出して組み合わせ数がこの通りになることを確認しよう!!)。 RM252 と RM241 で実際に観測できたデータ数と上図モデルの結果をまとめてみると次 のようになります。
表2 F2 1 GG-GG (1-
r
)2/4 2 GG-GK 2(1-r
)r
/4 3 GG-KKr
2/4 RM252 G G RM252 G G G G RM252 G G | | | | | | | | RM241 G G RM241 G K K G RM241 K K ① ① ① ③ ③ ① ③ ③ 4 GK-GG 2(1-r
)r
/4 RM252 G K K G | | | | RM241 G G G G ① ④ ④ ① 5 GK-GK 2(1-r
)2/4+2r
2/4 RM252 G K K G G K K G | | | | | | | | RM241 G K K G K G G K ① ② ② ① ③ ④ ④ ③ 6 GK-KK 2(1-r
)r
/4 RM252 G K K G | | | | RM241 K K K K ③ ② ② ③ 7 KK-GGr
2/4 8 KK-GK 2(1-r
)r
/4 9 KK-KK (1-r
)2/4 RM252 K K RM252 K K K K RM252 K K | | | | | | | | RM241 G G RM241 G K K G RM241 K K ④ ④ ④ ② ② ④ ② ②組換え価 r の2 マーカーで、実際に観測される遺伝子型の個数は、確率関数としてあらわ すことができます。9 面のサイコロを振って(各面が表になる確率が fa~fi で与えられてい る)出る目の数を計測し(a~i)、観測値の出る確率を考えるのです。その確率関数は多項分布 の形であらわされ(N=a+b+c+d+e+f+g+h+i)、
( ) ( ) ( ) ( ) ( )
( ) ( )
( ) ( )
i i h h g g f f e e d d c c b b a af
f
f
f
f
f
f
f
f
i
h
g
f
e
d
c
b
a
N
r
P
!
!
!
!
!
!
!
!
!
!
)
(
=
のようになります。!記号は階乗を表し、N!=N×(N-1)×‥×2×1 です。 期待頻度の重複を考慮すると上式はもう少し簡単になって、( ) ( )
( ) ( )
e e g c c h f d b b i a af
f
f
f
i
h
g
f
e
d
c
b
a
N
r
P
=
+ + + + +!
!
!
!
!
!
!
!
!
!
)
(
式
1
のようになります。 多項分布は、高校の数学で習う二項分布の拡張版です。二項分布はコインの表裏など 2 つ の背反な事象を表す確率関数で、N 回コインを投げて表が出る回数が k 回(裏の出る数は N-k回)であり表になる確率がrのとき、( ) (
k)
N kr
r
k
N
k
N
r
P
−
−−
=
1
)!
(
!
!
)
(
のように表せます。 もし多項分布の意味がすんなりと理解できなくても、考え込まないようにしましょう。式 1 を見直してください。a~i、Nはただの観測値の数です。fa,fb,fc,fe も変数はたった一つr の関数です。rは0.0~0.5 までの値をとります。Excel のシート上でa~iを記入しておき (RM252-RM241 の場合の値は表 2 を参照)、あるrに対するP(r)を Excel の式で計算させ、 rの値を0.001~0.500 まで 0.001 刻みで計算させると、次のようになります。 図 組換え価の推定 RM252-RM241の組換え価 0 0.0000001 0.0000002 0.0000003 0.0000004 0.0000005 0.0000006 P(r)を最尤法といいます。マーカー間の真の組換え価は神様しか知らないので、人間が組換え 価の情報を得るには、実際にマーカー遺伝子型の分離を観測し、観測データから組換え価 を推定するしかないことに注意してください。 最尤法の使用方法をより単純な確率モデル(二項分布)で説明しましょう。コインは表裏が 均等にできていれば表がでる確率は0.5 であると考えます。今、コインの表面の重さを裏面 より10%重くしてみました。表の出る確率はいくらになるでしょうか?、、、、この場合、実 際にそのコインを投げて表裏を計測し、観察値から予測をします。50 回コインを投げて表 は30 回出たとします。表の出る確率をrとして計測値のようになる確率を考えてみると、