Introduction
集団内の個々の細胞は均一な生理的刺激1、2、3、4に大幅に異なる応答を示すことができます。集団中の細胞の遺伝的変異は、応答のこの品種のための1つのメカニズムであるが、それでも細胞のクローン集団で、応答の変動性を向上させることができ、いくつかの非遺伝的要因もあります。例えば、個々のタンパク質および他の重要なシグナル伝達分子のレベルは、下流の遺伝子発現プロファイルの変化を生じる、セル単位で変化させることができます。さらに、遺伝子活性化は、バースト7,8,9あたりの転写物の比較的少数に限定することができる転写5,6の短時間のバーストで発生することができます。そのような遺伝子活性化における偶然性が大幅に生物学的応答の変動に寄与することができ、生理的刺激に応答して、微生物10におよび哺乳動物細胞1に2を選択的利点を提供することができます。変化の遺伝的および非遺伝的供給源の両方に、刺激に応答する任意の細胞の遺伝子発現プロファイルは、バルク応答の測定から得られた平均遺伝子発現プロファイルとは大きく異なっていてもよいです。個々の細胞は、刺激に応答して可変性を示す程度を決定することは、個々の細胞の単離のための技術は、目的の転写物の発現レベルを測定し、得られた発現データのコンピュータ分析を必要とします。
コストの広い範囲をカバーし、単一細胞における遺伝子発現をアッセイするためのいくつかのアプローチがあり、転写産物の数はプローブ、及び定量化の精度。例えば、単一細胞RNA-のSeqは、転写物のカバレッジと、個々の細胞の中で最も高度に発現される遺伝子のための別個の転写産物の数千を定量する能力の広い深さを提供しています。コストが減少し続けるが、しかし、このようなシーケンシングの深さに関連するコストは、法外することができます。逆に、in situハイブリダイゼーション（FISH smRNA） における単一分子RNAの蛍光は、目的の遺伝子につき合理的なコストであっても低発現遺伝子のための転写物の正確な定量を提供しています。しかしながら、標的遺伝子の少数は、このアプローチによって、所定の細胞においてアッセイすることができます。このプロトコールに記載定量的PCRベースのアッセイは、これらの技術の間の妥協点を提供します。これらのアッセイは、最大96個の細胞で一度に関心の96転写物まで定量化するために、マイクロ流体リアルタイムPCR装置を採用しています。前述の方法のそれぞれは、必要なハードウェアのコストを有しているが、任意の個々のqPCRアッセイのコストは比較的です低いです。このプロトコルは、（プロトコルADP 41、フリュー）マイクロ流体リアルタイムPCR装置の製造元が推奨する1から適応されます。 PCRベースのアプローチで各転写産物の絶対数の推定を可能にするために、我々は、複数の実験で使用できる準備標的遺伝子アンプリコンの内部統制を利用するプロトコルを拡大してきました。
この技術の一例として、MCF-7ヒト乳癌細胞において腫瘍抑制因子p53によって調節される遺伝子の発現の定量化は、11に記載されています。細胞は、DNA二本鎖切断を誘発する化学物質でチャレンジされています。以前の研究は、DNA二重鎖切断に対するp53応答は、p53レベル12の観点と異なる標的遺伝子11の活性化の両方で、個々の細胞における異質性を大いに発揮することが示されています。さらに、p53は100以上の発現を調節します細胞周期停止、アポトーシス及び老化13、14を含む多くの下流経路に関与する標的遺伝子を十分に特徴付け。各細胞におけるp53媒介性応答は、複雑かつ可変でもあるので、アプローチからシステムの利点の分析とは、約100の標的遺伝子は、以下に記載されるような個々の細胞、同時にプローブすることができます。 （例えば、単一細胞の単離および溶解のための代替方法として）わずかな変更を加えて、プロトコルは、容易に、哺乳動物の細胞型、転写、および細胞応答の広い範囲を研究するために適合させることができます。
適切な事前準備で、細胞選別および遺伝子発現測定のラウンドは、3日間にわたって、このプロトコルに従って行うことができます。 、事前に関心の転写産物を選択し、それらのtranscrからcDNAを増幅するプライマー対を識別し、検証：以下のタイミングが示唆されていますIPTSは、および標準とそれらのプライマーを用いてプライマーミックスを準備します。細胞治療の後、1日目に、収穫と逆転写および特定の標的増幅を行い、細胞を選別し、取り込まれなかったプライマーを除去するためにエキソヌクレアーゼでサンプルを扱います。 2日目、定量PCRを使用して選別した細胞上の品質管理を行います。最後に、3日目に、マイクロ流体定量PCRを使用して選別した細胞における遺伝子発現を測定します。 図1は、必要な手順をまとめたものです。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1.事前準備
- その発現が測定される関心の96遺伝子まで選択。
注：これらの遺伝子の少なくとも1つは、ACTBまたはGAPDHなどの「ハウスキーピング遺伝子」である必要があり、それは、実験に用いた条件下で、相対的に高い及び一定のレベルで発現することが知られています。この遺伝子は、正にソートウェル（ステップ8.1）、増幅サンプル（ステップ10.1）を識別するために使用されます。
注：例えば、実験のために、公知の機能11、14およびハウスキーピングコントロールを選択したとしてGAPDHの様々なp53の、直接の標的を十分に特徴付けられました。本研究で用いた標的遺伝子およびプライマー配列の完全なリストについては、参考文献11を参照してください。 - 15（このようなプライマー-BLASTとして科学文献またはプライマーの設計ツールを使用して、例えば ）目的の遺伝子に特異的な潜在的なプライマー対を特定します。プライマーを持っています合成だとヌクレアーゼフリー水に、100μMのストック濃度でそれらを維持します。
注：これらのプライマーは15-25塩基の長さでなければなりません。アンプリコンを90〜130塩基対（bp）の長さを生成します。存在する場合、ゲノムDNAを増幅する可能性を低減するために、エクソンジャンクションにまたがります。 60±1℃の融解温度を有します。そして、最小限の二次構造16を有しています。これらのプライマーは、目的の転写物からcDNAを増幅し、染料ベースのqPCR結合DNAにおける遺伝子発現を測定し、プライム逆転写に使用されます。 - プライマーを検証します。
- まず、目的の細胞からcDNAを得ます。製造業者の取扱説明書または標準的な分子生物学的方法17に従ってRNA採取キットを用いて細胞から収穫RNA。製造元の指示または標準メトに従ってランダムヘキサマーを用いた逆転写キットを用いてRNAをcDNAに逆転写しますDS 17。
- 各プライマー対について、DNAが前方、染料のマスターミックスを結合し、逆方向プライマー、および前のステップからのcDNA、提案した反応条件のマスターミックスメーカーの推奨と次のように定量PCRを設定します。溶融曲線の取得を含んでマスターミックスのメーカーが推奨する熱サイクルプロトコルを用いたリアルタイムPCRマシン上でこれらの反応を実行します。
- 定量PCRの後、2％のw / vのアガロースゲル上で増幅されたDNAサンプルを実行し、標準的なプロトコル17以下のDNA挿入色素を介したDNAを可視化します。
- cDNAの連続希釈液から標準曲線を構築することによって、プライマー対の効率を決定する標準的なプロトコル16に従って（ステップ1.3.1で行われました）。
注：良いプライマー対は、ゲル16の期待位置に単一の、明確に定義されたピークと融解曲線と単一バンドを生じるであろう、
- 精製されたDNAアンプリコンから基準を準備します。
- 各検証プライマー対の場合：
- プライマーと鋳型DNA濃度、反応条件、及びPCRサイクル条件、または標準的なPCRプロトコールを使用するためのポリメラーゼを製造業者の推奨に従って、高忠実度DNAポリメラーゼを用いてcDNAから関心領域を増幅"外部参照"> 17。
- 2％のw / vのアガロースゲル上で増幅されたcDNAを実行し、DNA挿入色素17でバンドを可視化します。きれいなカミソリの刃を使用してアンプリコンを表すバンドを切除し、マイクロチューブにゲル片を配置します。
- 製造者の指示に従って、例えば、アガラーゼベース抽出17スピンカラムゲル抽出キットなどの標準的なゲル抽出法を用いてゲルからアンプリコンを精製します。
- 製造業者の指示に従って蛍光ベースのdsDNA定量キットを使用してアンプリコンの濃度を測定します。 1μl当たり単位複製配列分子の数を得るために、アンプリコンの配列に基づいて、単位複製配列の既知の分子量測定のdsDNAの濃度を分割します。
- 2.0ミリリットルの低結合マイクロチューブでは、DNAの懸濁緩衝液の体積に精製されたアンプリコンの2μlを添加するようにアンプの最終濃度liconは5×10 8分子/μlです。
- すべてのアンプリコンは精製されている場合には、2.0ミリリットルの低結合マイクロチューブに各精製されたアンプリコンの18μLを兼ね備えています。 1,800μLの総容量を作るためにDNAを懸濁緩衝液を加えます。このように、 各アンプリコンの濃度が5×10 6分子/μlです。
注：既知の等モル量で関心のあるすべての遺伝子からのcDNAアンプリコンが含まれている。この混合物は、単一セルのqPCRにおける標準となります。 - 渦この「5E6」標準徹底的に、2000×gで10秒間スピン、低結合PCRチューブ中の20μlのアリコートに分割します。 -80℃でアリコートを保管してください。
- 各検証プライマー対の場合：
- 特定の標的増幅（STA）プライマーミックス（10倍）を準備します。
- 15ミリリットルコニカルチューブでは、（192まで）、各100μMのストックプライマー25μlのを兼ね備えています。 5,000μLの総容量を作るためにDNAを懸濁緩衝液を加えます。
注記：この混合物中の各プライマーの濃度は500 nmです。選別された細胞の各プレートは、プライム逆転写および特定の標的増幅にこのプライマーミックスの105.6μLを使用します。ここで与えられたプライマーミックスのボリュームは、ソートに45版を作るのに十分です。一度だけ行われる必要があるようにプライマー混合物の十分に大きな容積を作成することが望ましいです。必要に応じてボリュームを拡張できます。 - 110μlのアリコートに分割し、ボルテックスで混合し、-20℃でアリコートを保存します。
- 15ミリリットルコニカルチューブでは、（192まで）、各100μMのストックプライマー25μlのを兼ね備えています。 5,000μLの総容量を作るためにDNAを懸濁緩衝液を加えます。
- マイクロ流体チップベースのqPCRでの使用のために、アッセイミックス、各プライマー対のための1を準備します。
- 96ウェルPCRプレートの各ウェルには、所定の遺伝子のために対応する100μMリバースプライマーの3.0μlを、DNAの懸濁緩衝液の24.0μlを、そして2×アッセイローディング試薬の30.0μLを100μMのフォワードプライマーの3.0μLを組み合わせます。
- 渦このプレート、500×gで10秒間スピンし、-20℃で保存。
- ;およびEXOI（エキソヌクレアーゼI処理; 表2）RTSTA（ 表1転写および特定の標的増幅を逆）：2熱サイクルプログラムを作成します。
2.治療
- それらは実験的研究の条件に応じて、処理の際に所望のコンフルエンスまで成長するように、組織培養皿上で哺乳動物細胞をプレート。この研究で提示され、例えば、10％ウシ胎児血清（FBS）を含むRPMI中で、プレート4×10 5個のMCF-7のp53ビーナス細胞18 6当たり10cmディッシュ、100 U / mlペニシリン、100ミリグラム/ mlのストレプトマイシン、および250 ngの/ 5％CO 2で37℃で2日間mlのアムホテリシンBをインキュベート細胞彼らは50％の密集度に達するまで。
- 実験的研究に基づいて、関心の条件を確立するために、細胞を処理します。本研究で提示例えば、400 ngの/ mlのneocarzinoでMCF-7のp53-金星細胞をインキュベート3時間、8.5時間、14時間、または24時間スタチン。
細胞選別のための3の溶解緩衝液の準備
注：溶解緩衝液の単板を作るには、約1時間かかります。細胞選別は非効率的であるとは検出されなかった細胞で多くの井戸を得ることができるように、複数のプレートを作成し、ソートすることをお勧めします。
- PCRの準備のためにDNAの除染液とベンチ、ピペット、チューブラック、コールドPCRプレートホルダー、手袋を清掃してください。
- ヌクレアーゼフリー水を1.5mlチューブに株式RNase阻害剤を添加することにより、2.64 U /μlのRNアーゼ阻害剤を希釈します。
注：これは、RNアーゼ阻害剤は、サブマイクロリットル容積をピペッティングすることなく、溶解緩衝液に添加されることを可能にします。 - 表3に列挙した成分を組み合わせることによって、細胞のための溶解緩衝液を作ります。
- 別のチューブへの細胞のための溶解バッファーの92.4μLを移します。これは、アンプリコン標準の溶解バッファーになります。
- LOを使用しましたワット結合ピペットチップやマイクロ遠心チューブ、DNA懸濁液緩衝液で、それぞれ、200μlの3.1 PG /μlの6.2 PG /μlに希釈した大腸菌 DNAの500μlのを準備します。
- 細胞のための溶解緩衝液に3.1 PG /μlの大腸菌 DNAの184.8μLを加えます。
注： 大腸菌キャリアDNAは、非特異的プライマー結合の可能性の範囲を広げることにより、プライマー二量体は、PCR産物につながる可能性を低減するのに役立ちます。 - 基準を準備します。
- 5E5、5E4、... 5e0として6低結合マイクロ遠心チューブにラベルを付けます。
- 「5E5」のチューブにDNAを懸濁緩衝液の90μLを加えます。
- 他の5つの管のそれぞれに6.2 PG /μlの大腸菌 DNAの90μLを加えます。氷の上ですべてのチューブを保管してください。
注：規格にキャリアDNAを組み込むには、1セルのウェル中に匹敵する標準ウェルにおけるDNAの総質量になります。これは、PRを低減します基準濃度はチューブ又はピペットチップへのDNA結合によって影響されることobability。 - ストレージから（ステップ1.4で調製した各アンプリコンの5×10 6分子/μL、）「5E6」標準のアリコートを取ります。簡単にボルテックスし、液体がチューブの底にあることを保証するために、簡単に（500×gで5秒）スピン。
- 低結合ピペットチップを使用して、「5E5」チューブに5E6標準の10μlを添加します。簡単にボルテックス、スピン（500×gで5秒）、そして氷上に戻ってチューブを置きます。
- 「5E4」チューブに「5E5」管からのピペットを10μl。簡単にボルテックス、スピン（500×gで5秒）、そして氷上に戻しました。
- すべてのチューブは、標準の希釈になるまで、各チューブが前のチューブから10μlのを受けて、これらの連続希釈液を繰り返します。
- キャップで12 PCRチューブの行を準備します。各チューブに、「細胞」の溶解緩衝液67μLを配布します。チューブに蓋をし、簡単にスピンダウン（5500×gで秒）必要であれば。
- 12チャンネルピペットを使用して、PCRプレートの最初の7列の各ウェルにPCRチューブの行から「細胞」の溶解緩衝液（ 図3）の9μLを分配します。
- プレートの最後の行（ 図3）の各ウェルに、「標準」の溶解緩衝液7μLを配布します。
- 低結合ピペットチップを使用して、プレートマップ（ 図3）に記載の最下行の各ウェルに標準の2μlを添加します。
注： など 10分子、5×10 0個/μlの標準量の2μlの - 10細胞および100細胞ウェルに3.1 PG /μlの大腸菌 DNAの2μLを加えます。無細胞ウェルに6.2 PG /μlの大腸菌 DNAの2μlを添加します。
注：各1セル、10セル、および100細胞はよく単一のヒト細胞では、ゲノムDNAの質量を近似、 大腸菌 DNAの6.2 pgのを持っています。各標準および無細胞ウェルが12.4 pgのO持ちますF 大腸菌 DNA、およそ2つのヒト細胞の価値があります。 - 細胞はソーティングのための準備が整うまで4℃で滅菌した熱シール、500×gで5秒間スピンし、ストアを使用してプレートをシール。
4.セルソーターのセットアップ
- プログラムプレートマップ（ 図3）によれば、96ウェルプレートに選別し、蛍光活性化セルソーター。何の細胞をウェルH1-H8またはH11-H12に分類されるべきではありません。
- セルソーターは、ウェルPCRプレートに並べ替えを目的としていることを確認します。
- 可能な限り最小にソートストリームの角度を設定します。この設定は、細胞が側面を打つのではなく、井戸の底に分類される可能性を増大させます。
- マシンを目指すためには、機械19の指示に従って、空のPCRプレートを密封し、空のプレートのシールにPBSの「並べ替え」の液滴にテストストリームを使用します。液滴が目に着陸する必要がありますウェルの中央より上の位置でのシールの電子表面。最良の結果を得るために、プレートの長さと幅を横切って目的を確認します。
- 液滴が正しい位置に上陸しない場合は、シールの表面を、それらを拭き取り、機械メーカーの指示に従って、セルソーターを再キャリブレーション、およびソート・ストリーム内の液滴を正確に堆積されるまで繰り返します。
5.細胞収穫と選別
- 細胞株のような適切な細胞をトリプシン処理します。例えば、MCF-7細胞について、0.05％トリプシン1mlで細胞をすすぎ、細胞から培地を吸引し、37℃、5％CO 2で5分間0.05％トリプシン2mLで細胞をインキュベート。
- トリプシン処理した細胞に細胞培養培地を8mlを加え、15mlチューブに細胞懸濁液の全10 mlで移します。 400×gで4分間遠心分離します。
- 上清を吸引除去し、PBS + 2％FBS中で細胞ペレットを再懸濁します。<BR />注：小さい再懸濁体積は細胞を集中し、細胞選別を容易にします。 500μlの再懸濁量は50％の集密度で細胞を6 10cmディッシュに適しています。
- 細胞の塊を除去するために、40ミクロンのフィルターを通して、チューブサイトメトリーピペッティングを流すmlの5に細胞懸濁液を転送します。
- 液体はウェルの底にあることを保証するために、細胞を選別する前に、4℃で30秒間、400×gでの溶解緩衝液プレートを遠心。
- セルソーターに細胞懸濁液をチューブを挿入します。溶解バッファーの板にシールを開き、慎重にプレートマップに従って、セルソーターの製造元の指示19を次のプレートに関心の細胞を選別。単一細胞が最大純度でソートされていることを確認するには、「単一セル」モードでマシンを実行し、前方および側方散乱光19に基づいて、標準的なフローサイトメトリーゲーティングを使用します。
注：中実験例、個々のセル11の最も明るい15％を単離するために、p53を金星の蛍光に基づいて、ゲートの細胞。 - 4℃で1分間、400×gで新しい、滅菌熱シールや遠心分離機でプレートをシール。サーマルサイクラーに板を置き、RTSTAプログラム（ステップ1.7を参照）を実行します。
注：これは分類された細胞からのmRNAの逆転写を行い、次いで、目的の転写産物のcDNAを増幅します。
6.エキソヌクレアーゼ処理
- 表4に列挙した成分を組み合わせることによって希釈エキソヌクレアーゼIを混ぜます。キャップで12 PCRチューブの行を準備します。各チューブに、希釈エキソヌクレアーゼIの30.5μLを配布します。チューブに蓋をし、液体がチューブの底にあることを保証するために、簡単に（500×gで5秒）スピン。
- RTSTAが終了すると、500×gで5秒間のサンプルプレートをスピン。 12チャンネルピペットを用い、希釈エキソヌクレアーゼの3.6μLを分配しますIプレートの各ウェルに。
- 滅菌熱シールでプレートをシールし、（500×gで5秒）を簡単にスピン。サーマルサイクラーに板を置き、EXOIプログラム（ステップ1.7を参照）を実行します。
注：彼らは将来のPCRに干渉しないように、このステップは、特定の標的増幅後に残っている任意のプライマーを劣化させます。
7.サンプル希釈
- エキソヌクレアーゼ処理が終了すると、サンプルプレートの各ウエルにTE緩衝液32.4μlを添加します。これは、任意の停止点であり、サンプルは数日間-20℃で保存することができます。
定量PCRを使用した8ソート品質管理
注：細胞選別が完全に効率的ではないので、このステップは、ソートプレートのウェルは、実際にセルを受信し識別することが必要です。これらのサンプルは、さらなる分析のために使用することができます。
- qPCRにより、各試料中のハウスキーピング遺伝子の発現を測定します。
- ステップ1.1で選択したハウスキーピング遺伝子のいずれかのプライマーを用いてポリメラーゼの製造元のプロトコルに従って96の反応のための染料のqPCRマスターミックスをDNA結合900μlのを準備します。この反応のため等の反応緩衝液、ポリメラーゼ、プライマーの正確な濃度は、使用される特定のポリメラーゼに依存するであろう。
- 96ウェルPCRプレートの各ウェルにマスターミックスの9μLを配布します。 PCRプレートの対応するウェルにサンプルプレートからサンプルを1μlを追加します。
- マスターミックス製造元が推奨する熱サイクルプロトコルを使用するか、または標準プロトコル20次のリアルタイムPCR装置にプレートを実行します。
- 無細胞ウェルからの測定値が低い（バックグラウンド）レベルを示していることを確認します遺伝子発現または全く発現。
- 無細胞ウェルに似た高い遺伝子発現および発現または全く発現のバックグラウンドレベルで陰性サンプルで陽性サンプル、：1 - 細胞ウェルからの測定値は、二つのグループに明確に分けることを確認します。
注：高い遺伝子発現を有するサンプルは、実際の選別された細胞を表します。バックグラウンド発現を有するものが、さらなる分析から除外すべきです。 - 10細胞および100-細胞ウェルからの測定値は、非効率的なソートがそれらのウェルが望まれるより少ない細胞を持たせることができる可能性を考慮して、正の1セルウェルよりも高い遺伝子発現を反映していることを確認します。
- 標準ウェルから測定値が均等に（C t値が均等に分布していること、すなわち ）対数空間に分布していることを確認します。
マイクロ流体定量PCRを使用した9遺伝子発現測定
注：このセクションのすべてのステップのために、試薬中の気泡の形成を最小限にするための唯一の最初のピットストップにピペット。
- 8 PCRチューブの行を開きます。サンプルの複数のプレートがある場合、8 PCRチューブの2行目を開き、チューブ1-8を標識します。この行には、複数のプレートから基準をプールするためになります。
- 1.5mlチューブにおいて、色素サンプルローディング試薬を結合するDNAの36μlのDNA結合色素定量PCRマスターミックスの360μLを混ぜます。渦簡単にと2000×gで10秒間スピン。各チューブに48μlずつ配置し、8 PCRチューブの標識されていない行にこの混合物（396μl）を分割します。
- 8チャンネルピペットを用いて、各ウェルにマスターミックス/サンプルローディング試薬混合物の3.3μLを配布新しい96ウェルPCRプレート。このプレート、ラベル「サンプルミックスを。」
- サンプルの各PCRプレート用：
- 渦500×gで10秒間10秒とスピンのためのプレート。カバーを取り外し、そのは、サンプルに対応するウェルにそれぞれ正1 - 細胞サンプルの2.7μLを移す1-細胞ウェルのためのプレートをミックスプレートマップに示されています。
- サンプルの一つだけのプレートがある場合は、そのは、サンプルに対応するウェルに移し、各増幅された標準の2.7μlをプレートマップに従ってプレートをミックスします。
- サンプルの複数のプレートがある場合、PCRチューブの標識された列内の対応する位置にそれぞれ増幅標準の2.7μLを移します。すでにそこに任意の規格にサンプルの整流板から基準を追加します。封止フィルムを用いて行ったときのサンプルプレートをシール。
- サンプルの複数のプレートがある場合：
- 基準を含むPCRチューブの列を渦2,000×gで10秒間10秒とスピン。よくサンプルに対応するに混合標準の各チューブから2.7μLを移し、プレートマップに従ってプレートをミックスします。
- プレートマップ上に表示ヌクレアーゼフリー水2.7μlのNTC井戸をロードします。必要になるまで、サンプルを4℃でプレートや店舗をミックスシール。
- ピールステッカー新しいマイクロ流体定量PCRチップの底オフ。まだ上にキャップと制御線流体の注射器を使用して、それを緩めるために、各アキュムレータスプリングを押し下げ、制御線液の150〜200μlの各アキュムレータを注入します。
- 搭載機にチップを挿入し、「プライム」スクリプトを実行します。これは〜20分かかります。
- プライミングは、渦を行って、アッセイ混合物のプレートをスピンダウンした場合、サンプルミックス（ステップ1.6で調製しました）。
- 8チャンネルピペットを用いて、図中によるマイクロ流体定量PCRチップの左側に各アッセイミックスの4.5μLを移します
- 図6によるマイクロ流体定量PCRチップの右側に、各サンプルミックスの4.5μLを移します。搭載機にチップを挿入し、「ロード・ミックス "スクリプトを実行します。これは〜90分かかります。
- 、マイクロ流体リアルタイムPCR装置の電源をオンにし、データ収集ソフトウェアを起動し、それをウォームアップするためにランプを点灯。これは〜20分かかります。
- ロードミックススクリプトが実行されると、チップ-この全体の何行がないことを確認すると、ロードの問題となります。慎重にスコッチテープやラボのテープを使用して、チップからほこり粒子を除去。マイクロ流体リアルタイムPCR機にチップを挿入し、データ収集スクリプトを実行します。結果を検査し、さらなる分析のためのデータをエクスポートする解析ソフトウェアを使用してください。
- 定量PCR増幅曲線は20（曲線は、標準的なS字型増幅曲線に似ていない、すなわち 、）低品質を示している測定値を削除します。
- 30以上の低品質またはゼロ測定またはのためのハウスキーピング遺伝子の発現を有するサンプルを削除する（ステップ1.1参照）他のサンプル（ 図7）の大部分よりもはるかに低いです。
- 各遺伝子について、各試料からの溶融ピーク位置の中央値として正確な融解ピーク位置を推定します。測定値が中央値から以上1°Cに位置する溶融ピークを持っている場合、考慮20からその測定値を削除します。
- 各遺伝子について、スプレッドシートやスクリプト言語20を使用して、各サンプル中のmRNAの数を推定するための標準曲線を作成します。
注：この実験での標準は10、100、千から成り、または10,000の分子dsDNAの関心の各転写物に対応します。逆転写は完全に効率的であった場合には、mRNAおよびcDNAは一本鎖であるため、基準は、20、200、2,000、およびmRNA 20,000の分子に対応します。実際には、これはそうではありませんので、我々は、E RTは逆転写の効率であるE RT、×分子単位の推定値として測定値を表現します。- 各標準ウェルについて、 メートル 、増幅前のssDNAの分子の数（dsDNAの分子の数の2倍）（ 例えば、20分子、200分子など ）を識別します。
- 各標準ウェルについて、定量PCR 20からのC t値を識別します。
- mおよびC tの値と、各遺伝子のパラメータaおよびbを有する標準曲線を生成するために、以下の式を用いて線形回帰を実行し、のようなR 2の値を見つけます適合度のn個の指標：
- 標準曲線から、各遺伝子のPCR効率Eを計算します ：
注：Eの値1.0よりもはるかに大きい（ 例えば、E> 1.1）は、物理的に現実的ではありません。 1.0よりもはるかに少ない線形回帰（ 例えば、R 2 <0.985）のためのR 2値は、基準が確実に動作していないことを意味します。これらの問題は、一般的に最も低いアンプリコンの標準は、信頼性の高い定量化の閾値を下回った、従って最低の標準は、線形回帰から除外されるべきであることを示しています。適切な線形回帰は、少なくとも3つのユニークな基準を使用して作製することができない場合には、遺伝子の発現の測定を捨てます。 - 線形回帰が適切である場合、得られた適合パラメータaおよびbを使用して<mRNAの数を推定するために、各遺伝子について/ em>の測定されたC tの値から、各単一細胞サンプル中の対応する遺伝子のために（E RT×分子の単位で） メートルの東側を分子：
- また、m個のEST = 0として任意の測定メートルの東側 <1を指定し、= 0メートルの東側として定量PCRで増幅し、したがってなしのC t値なしで任意の測定値を指定します。
注：estの最低未満または回帰で使用される最も高い標準よりも大きいmの測定は、標準曲線から外挿によって発見されました。これらの推定値は必ずしも無効ではありませんが、慎重にみなされるべきです。
- 単一細胞の遺伝子発現データを視覚化します。
- beeswarmプロット（ドーン、J.聖、T.、2012、http://www.mathworks.com/matlabcentral/fileeを作ります各点は、特定の細胞における遺伝子発現を表す、XCHANGE / 37105）;例えば、 図8に示されています。
- 「バイオリン」の幅がで特定のレベルのまわりの遺伝子発現測定の周波数を表しているバイオリンプロット（ドーン、J.、2012、http://www.mathworks.com/matlabcentral/fileexchange/23661）を作ります分析した細胞の集団。例えば、 図8に示されています。
注：より洗練された分析は、遺伝子21のペア間の発現の関係をプローブ22および遺伝子発現ネットワークを推測することができます。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
プロトコルの一般的な概要は、細胞治療、FACSによる単一細胞の単離、単一細胞溶解物からのcDNAライブラリーの生成および前増幅、単一細胞cDNAライブラリーを確認するステップを含む、 図1に示されていますソートされたウェル、および定量PCRにより遺伝子発現を測定します。
単一細胞の単離および遺伝子発現分析の準備のために、最初に、目的の各標的遺伝子のために有効なオリゴヌクレオチドプライマー対を特定する必要があります。 図2は、プライマーの品質管理方法の2つの例を示しています。溶融曲線は、プライマーペアがテストされていると定量PCRを以下に生成されます。有効なプライマーは、単一のピーク（ 図2A、緑の曲線）との融解曲線になります。複数のピーク（ 図2A、青の曲線）または顕著な肩の曲線（Figur場合E 2（a）、赤曲線）、これは、複数のPCR産物の存在を示し、存在し、プライマーを再設計する必要があります。第二の品質管理方法としては、アガロースゲル電気泳動は、視覚的に正しいサイズの単一バンドが（ 図2B）、試験した各プライマー対のために存在すること。確認するために使用することができます間違ったサイズの複数のバンドやバンドはプライマーが（ 図2B）が有効でないことを示しています。
プライマー検証し、正しいアンプリコンの定量PCR標準のライブラリーの生成に続いて、細胞選別を行うことができます。例えばソート方式は、図3に示されています。各ソートプレートは、各標的アンプリコンの10、100、1,000または10,000コピーして、アンプリコン標準の希釈系列を含むウェルを含める必要があります。また、何の細胞、10細胞、または100細胞は、単一細胞ウェルに加えて、ソートされていないそこに井戸を含めることをお勧めします。非効率性を選別するために、単一の細胞が正常に堆積されたプレートのウェルを同定することがしばしば必要です。この検証ステップのために、逆転写のための手順および特定の標的増幅後、定量PCRがハウスキーピング遺伝子（ 図4）の発現を測定するために行われるべきです。アンプリコン規格のウェル（ 図4、黒の曲線）は、等間隔の増幅曲線を示すべきです。 （;赤い曲線図4）、「10セル「ウェル（ 図4;青の曲線）、および「100細胞」ウェル（ 図4も「無細胞」ウェルに対応する曲線に明確な分離があるはずです;紫色の曲線）。明確な分離にも成功した細胞成長（ 図4に対応する「1セル「井戸で発生します; C t値と緑の増幅曲線少ない無細胞対照のためのものよりますがグラム無細胞対照用のものに近いC t値と緑の曲線）;失敗した細胞成長（ 図4対）10-細胞対照のためのものよりreater。単一細胞溶解物を有するウェルが同定されたら、（潜在的に複数のプレートから）をウェルからのcDNAは、ウェルあたり単一の細胞溶解物と、マイクロ流体のqPCRチップに転送することができます。例えば、 図5は、異なる3つのqPCR分析のための単一の新しいプレートに96ウェルプレートをソートから単一細胞からcDNAを組み合わせた仮想的な配列を示しています。
多くの欠落またはゼロ遺伝子発現測定（ 図7、矢印で示されている行）または異常に低いハウスキーピング遺伝子の発現によって示されるように、マイクロ流体定量PCRの結果を分析するために、可能性の高いセルを受信しなかったサンプルは、最初に、削除する必要があります。各アッセイのために、との任意の測定値は、パーソナルプラグインのものと一致しないピークを溶かします同じアッセイ中のRサンプルも削除する必要があります。不良サンプルとの測定値を除去した後、線形回帰が、各単一細胞サンプル中の目的の各転写物の絶対数を推定するために、アンプリコンの規格に対応した測定を行うことができます。 図8は beeswarmとバイオリンプロットとして、逆転写効率×分子の単位で測定し、これらの遺伝子発現の推定値の可視化を示しています。この例では、未処理の細胞は多くの細胞が全くCDKN1Aを発現していないとし、二桁に及ぶレベルで遺伝子を発現する他の人と、CDKN1A発現レベルの広い範囲を示しています。 3時間ネオカルチノスタチン（NCS）で処理した後、ほとんどの細胞がはるかに高いレベルでCDKN1Aを発現し 、ばらつきが大きさの約単一の注文に減少します。治療が増加する期間としては、CDKN1A発現の平均レベルが低下し、変動性の式が展開されます。
図1： 単一細胞の遺伝子発現プロファイリングのステップの概要。処理した細胞を回収し、直接、溶解緩衝液にFACSを介してソートされます。関心対象のmRNA種は、逆転写され、そして得られたcDNAをPCRにより増幅されます。各試料中のハウスキーピング遺伝子の発現は、サンプルが実際にセルを受信したかを決定するためにqPCRにより測定されます。この品質管理試験に合格した試料では、最大96の遺伝子の発現は、マイクロ流体のqPCRを使用して測定されます。この図は、以前の出版物11から変更されています。 この図の拡大版をご覧になるにはこちらをクリックしてください。
19 / 55219fig2.jpg "/>
図2： プライマー試験結果の一例。 （A）は、定量PCRから曲線を溶かします。青色の融解曲線（TNFRSF10D）は、複数のピークを有しており、赤色の融解曲線（TRIM22）は、最初のピークに近い第二のピークに達する」肩を、 "持っています。これらは、これらの反応に使用したプライマーは、複数の標的を増幅して再設計する必要があることを示しています。緑色の融解曲線（SCN3B）は、単一の標的を増幅するプライマーと一致し、単一のピークを有しています。 （B）増幅されたDNAをアガロースゲル上で実行されます。 SCN3B、TRIM22、およびTRPM2のためのレーンは複数のバンドが含まれています。これらのDNAアンプリコンを作製するために使用されるプライマーは、複数の標的を増幅して再設計する必要があります。 TNFRSF10Dのためのレーンは、単一の標的を増幅するプライマーと矛盾単一のバンドを持っています。それはTNFRSF10Dのものは、ゲル試験に合格するが、融解曲線を失敗しながらSCN3Bのためのプライマーは、融解曲線の試験に合格したがゲル試験を失敗していることは注目に値しますテスト;プライマーの検証のこれらの2つの方法は、相補的です。この例では、1つの4つすべてのプライマーセットを再設計する必要があると結論することができます。 この図の拡大版をご覧になるにはこちらをクリックしてください。
図3：細胞選別のためのプレートマップ。このプレートマップは、1セルのウェル（赤）、10-および100-細胞ウェル（濃い赤）、標準（緑）、および無細胞ウェル（白色）を含みます。各プレート上の各標準の2つの複製を含めると、規格の変動を補償するのに役立ちます。 10細胞および100-細胞ウェル含めると、遺伝子発現のための健全性検査として有用です。無細胞ウェルは陰性対照として役立ちます。 トンの拡大版をご覧になるにはこちらをクリックしてください。彼の姿。
図4：GAPDHの発現を測定するソート品質管理から定量PCR増幅曲線の例。無細胞サンプル（赤）の発現のバックグラウンドレベルを示しています。 1 - 細胞サンプル（緑）は、2つのグループ、高発現と無細胞サンプルと同様の低発現を有する1と1に分けます。高発現サンプルは、最も可能性の高いソート時に実際のセルを受け取りました。低発現サンプルは、ほとんどなかったです。 10細胞（青）と100細胞（紫）のサンプルは、1セルのサンプルよりも高い発現を示します。 10-細胞サンプルが原因ソート非効率性の最も高い発現1 - 細胞サンプルに近いです。期待通りに（黒）規格が均等に分散されています。 この図の拡大版をご覧になるにはこちらをクリックしてください。 </ P>
図5：サンプルミックスプレートマップの例。このプレートマップは、qPCRのためのすべての3つのプレート（H1-H8）、無細胞対照（H9-H10）、および無テンプレート対照から混合標準（H11-H12、三つの異なるプレート（行AG）から1セルのサンプルが含まれて）。 この図の拡大版をご覧になるにはこちらをクリックしてください。
図6： ロードの順序とマイクロ流体定量PCRチップの図（1、2、3 ...）。ノッチ（A1）は、プレートの左上隅です。 THIの拡大版をご覧になるにはこちらをクリックしてください。sのフィギュア。
図7： 定量PCR結果のヒートマップ。各行は、1-細胞試料、標準、または対照を表します。各列は、測定された転写物を表します。黒い四角またはxs（矢印）の異常な数を持つ行は、おそらく実際のセルを受信しなかったサンプルを示しています。 この図の拡大版をご覧になるにはこちらをクリックしてください。
図8：Beeswarmとバイオリンプロット。 beeswarmプロット（青い点）において、各点は、特定の細胞における遺伝子発現を表します。バイオリンプロット（灰色の形状）において、「バイオリン」の幅は、遺伝子発現の頻度を表すmeasu分析した細胞の集団中の特定のレベルのまわりrements。水色のバーは、最も低い標準を表します。青いバー内の測定値が外挿されている間、青色バーの上の測定は、規格間で補間されています。ここに示したデータは、DNA二重鎖切断を誘導するために示された時間のためのネオカルチノスタチン（NCS）で処理された金星タグ付きのp53を発現するMCF-7細胞でCDKN1A式の測定値です。この図は、以前の出版物11から変更されています。 この図の拡大版をご覧になるにはこちらをクリックしてください。
| 調子 | 温度 | 時間 |
| ホールド | 50°C | 15分 |
| ホールド | 95°C | 2分 |
| 20サイクル | 95°C | 15秒 |
| 60°C | 4分 | |
| ホールド | 4°C | ∞ |
表1：逆転写および特定の標的増幅（RTSTA）サーマルサイクリング・プログラム。
| 調子 | 温度 | 時間 |
| ホールド | 37°C | 30分 |
| ホールド | 80°C | 15分 |
| ホールド | 4°C | ∞ |
表2：エキソヌクレアーゼI処理（EXOI）サーマルサイクリング・プログラム。
| 成分 | ボリューム/ウェル（μL） | / 10％の超過ワット96容量（μL） |
| 2倍の反応ミックス | 5.00 | 528.00 |
| 逆転写酵素/ポリメラーゼミックス | 0.20 | 21.12 |
| 10倍STAプライマーミックス | 1.00 | 105.60 |
| 2.64 U /μlの（希）RNA分解酵素阻害剤 | 0.02 | 2.00 |
| ヌクレアーゼフリー水 | 0.78 | 82.48 |
| 合計 | 7.00 | 739.20 |
表3：細胞のための溶解緩衝液の成分。
| 成分 | ボリューム/ウェル（μL） | / 10％の超過ワット96容量（μL） |
| ヌクレアーゼフリー水 | 2.52 | 266.00 |
| エキソヌクレアーゼIの反応バッファー（10倍） | 0.36 | 38.00 |
| エキソヌクレアーゼI（20 U /μl）を | 0.72 | 76.00 |
| 合計 | 3.60 | 380.00 |
表4：エキソヌクレアーゼ混合物の成分が増幅されたcDNAサンプルを処理します。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
我々は、培養で増殖させた接着性細胞の集団から個々の哺乳動物細胞を単離するため、各セル内の約96の遺伝子の発現をアッセイするための方法を提示しています。この方法がうまく機能するために良好な事前準備が重要です。プライマーは、単セル測定の品質を決定するように、特に興味のある転写物（1.2〜1.3ステップ）に特異的な設計およびテストプライマーペアは、時間がかかるが、重要なステップです。信頼性のあるプライマー対が得られたら、それらは、目的の転写物からcDNAを増幅するために使用されます。アンプリコンは、その後の規格（ステップ1.4）を作るために等モル量で組み合わされています。標準は絶対転写産物数を推定するために必要とされるので、このステップは重要です。規格は、一度作っ等分して、細胞選別および遺伝子発現測定値の後続のすべてのラウンドのために使用されるべきです。同様に、細胞がソートされた日に、それは特にimportanです溶解バッファーのプレートを作る際に慎重に（3.7および3.11ステップ）低結合ピペットチップやマイクロ遠心チューブを使用して標準を希釈するトン。 （ステップ4.2）セルソーターを目指して、別の重要なステップです。これは、細胞が正常にPCRプレートに溶解緩衝液9μlに目標を達成するためには非常に慎重に行われなければなりません。
研究ニーズの多様に適合させることができるプロトコルにはいくつかの変更があります。ここでは、遺伝子発現を定量するための比較的手頃な方法を提供し、DNA結合染料ベースのqPCRアプローチに焦点を当てました。しかし、この方法は、任意の増幅されたDNAから、高いバックグラウンドシグナルを生成する可能性を有し、さらにそのオフターゲット配列は、蛍光シグナルが生じます。このような背景には、特定の標的を増幅するために使用されるプライマー対は、単一のピークを有する融解曲線と正確なサイズの単一のPCR産物（ 図2をもたらすことを確実にすることによって最小化することができます）。バックグラウンドは、このような品質管理方法を用いた後に依然として問題である場合には、プローブベースの定量PCRアプローチは、試薬の高いコストであるが、オフターゲット検出を最小化するために使用することができます。別のオプションは、プライマーのセットは、写しを逆転し、転写産物の大領域と、その大領域に含まれるより小さな領域を増幅するプライマーの第2のセットを、増幅するRTSTA工程で使用したネステッドPCRアプローチであろう、定量PCRにより遺伝子発現を測定するために使用されます。このアプローチは、染料ベースのqPCR 23 DNA結合の特異性を改善することが示されています。
いくつかの方法は、遺伝子発現分析のために個々の細胞を単離するために利用可能です。細胞分離方法の選択は、（例えば、蛍光活性化セルソーターなど）設備の利用可能性を含むいくつかの因子に依存する（例えば、組織からの初代細胞に対する樹立細胞株など）を評価するための細胞の供給源、または細胞を単離する必要があるとスピード。このプロトコルで示す例では、我々は明らかに検出可能な蛍光タグ化タンパク質を発現する細胞株のFACSを用いました。 FACSは、希少細胞の検出能力と比較的急速な細胞分離の利点があります。この方法の1つの課題は、96ウェルPCRプレート中で溶解緩衝液の必要な小体積に関心対象の細胞の堆積は非効率的であることができ、成功した分離を確実にするために、細胞選別のジオメトリの最適化を必要とし得ることです。低速および/またはより高いコストでの細胞分離でより高い精度をもたらす可能性が別のアプローチは、個々の細胞のマイクロピペット、腫瘍サンプルから特定の細胞のレーザーキャプチャーマイクロダイセクション、またはマイクロ流体システムの使用（ 例えば、フリューC1を含みますシステム）。
ここで説明されたプロトコルの1つの主要な利点は、それが内部統制、精製されたPCRアンプリコンの希釈系列、Eは各セル内の各転写物の絶対数の推定をNABLE。これらの基準なしに、遺伝子発現の唯一の相対的レベルは、C t値から導出さ計算に基づいて取得することができます。しかしながら、これらの規格で、絶対的な転写産物の数は、理想的には、正確に定量化されたアンプリコンの標準の範囲内挿（ 図8）を介して、推定することができます。遺伝子発現の定量化の任意のPCRベースの方法と同様に、正確な絶対的な定量化は逆転写を含む各プロセスの効率性についての知識を必要とします。
単一細胞内の転写レベルを定量化における現在の課題は、多くの方法の検出限界は、トランスクリプトームのかなりの割合の信頼性の定量化を防止し、セル24当たり10-100のmRNAであると推定されていることです。別のアプローチとして、1は、特に関心のある目標を定量化するためにsmRNA FISHを使用することができます非常に少数の転写産物を有するものを含むお尻= "外部参照"> 25、26、。 smRNA FISHの進歩は、かつて27で検出することができる別個の転写産物の数と適切な機器与え、一度28でプロファイリングすることができる細胞の数を拡大しています。 smRNA FISHは独自の制限がありますが（潜在的な偽陽性と偽陰性の転写産物の検出、細胞あたりわずか数標的遺伝子への制限を、機器や試薬のコストが）、それは結果を補完し、検証するための強力な方法を提供することができます関心対象の標的遺伝子のサブセットの定量PCRベースの分析。
Subscription Required. Please recommend JoVE to your librarian.
Acknowledgments
私たちは、このプロトコルの開発中に細胞選別を行うの彼女の援助のためにCCR ETIBフローサイトメトリーコアでV.カプールに感謝したいと思います。また、このプロトコルの開発中に定量PCRを行う際の援助のためにM. RaffeldとCCR LP分子診断ユニットとJ.朱とNHLBI DNA配列決定およびゲノミクスコアに感謝します。この研究は、NIHの学内プログラムによってサポートされていました。
Materials
| Name | Company | Catalog Number | Comments |
| RNeasy Plus Mini Kit | Qiagen | 74134 | |
| High Capacity cDNA Reverse Transcription Kit with RNase Inhibitor | ThermoFisher | 4374966 | |
| Phusion High-Fidelity DNA Polymerase | New England BioLabs | M0530S | |
| QIAquick Gel Extraction Kit | Qiagen | 28704 | |
| Quant-iT High-Sensitivity dsDNA Assay Kit | ThermoFisher | Q33120 | |
| 2.0 ml low adhesion microcentrifuge tubes | USA Scientific | 1420-2600 | |
| DNA Suspension Buffer | Teknova | T0221 | |
| Axygen 0.2 ml Maxymum Recovery Thin Wall PCR Tubes | Corning | PCR-02-L-C | |
| GE 96.96 Dynamic Array DNA Binding Dye Sample & Assay Loading Reagent Kit | Fluidigm | 100-3415 | |
| HyClone RPMI 1640 media | GE Healthcare Life Sciences | SH30027.01 | |
| Fetal Bovine Serum, Certified (US) | ThermoFisher | 16000-044 | |
| Antibiotic-Antimycotic Solution | Corning | 30-004-CI | |
| Neocarzinostatin | Sigma | N9162 | |
| ELIMINase | Decon Labs | 1101 | |
| SUPERase-In | ThermoFisher | AM2696 | |
| CellsDirect One-Step qRT-PCR Kit | ThermoFisher | 11753500 | |
| E. coli DNA | Affymetrix | 14380 10 MG | |
| ThermalSeal Sealing Film, Sterile | Excel Scientific | STR-THER-PLT | |
| BD FACSAria IIu | BD Biosciences | ||
| HyClone Trypsin 0.05% | GE Healthcare Life Sciences | SH30236.01 | |
| PBS, 1x | Corning | 21-040-CV | |
| Falcon 40 µm Cell Strainer | Corning | 352340 | |
| Exonuclease I | New England BioLabs | M0293S | |
| SsoFast EvaGreen Supermix with Low ROX | Bio-Rad | 172-5210 | |
| 96.96 Dynamic Array IFC for Gene Expression (microfluidic qPCR chip) | Fluidigm | BMK-M-96.96 | |
| IFC Controller HX (loading machine) | Fluidigm | ||
| BioMark or BioMark HD (microfluidic qPCR machine) | Fluidigm | ||
| Real-Time PCR Analysis software | Fluidigm | ||
| MATLAB software | MathWorks |
References
- Feinerman, O., Veiga, J., Dorfman, J. R., Germain, R. N., Altan-Bonnet, G. Variability and Robustness in T Cell Activation from Regulated Heterogeneity in Protein Levels. Science. 321 (5892), 1081-1084 (2008).
- Spencer, S. L., Gaudet, S., Albeck, J. G., Burke, J. M., Sorger, P. K. Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis. Nature. 459 (7245), 428-432 (2009).
- Geva-Zatorsky, N., Rosenfeld, N., et al. Oscillations and variability in the p53 system. Mol. Syst. Biol. 2 (1), (2006).
- Colman-Lerner, A., Gordon, A., et al. Regulated cell-to-cell variation in a cell-fate decision system. Nature. 437 (7059), 699-706 (2005).
- Chong, S., Chen, C., Ge, H., Xie, X. S. Mechanism of Transcriptional Bursting in Bacteria. Cell. 158 (2), 314-326 (2014).
- Raj, A., Peskin, C. S., Tranchina, D., Vargas, D. Y., Tyagi, S. Stochastic mRNA Synthesis in Mammalian Cells. PLoS Biol. 4 (10), e309+ (2006).
- Senecal, A., Munsky, B., et al. Transcription Factors Modulate c-Fos Transcriptional Bursts. Cell Rep. 8 (1), 75-83 (2014).
- Dey, S. S., Foley, J. E., Limsirichai, P., Schaffer, D. V., Arkin, A. P. Orthogonal control of expression mean and variance by epigenetic features at different genomic loci. Mol. Syst. Biol. 11 (5), 806+ (2015).
- Dar, R. D., Razooky, B. S., et al. Transcriptional burst frequency and burst size are equally modulated across the human genome. Proc. Natl. Acad. Sci. U.S.A. 109 (43), 17454-17459 (2012).
- Thattai, M., van Oudenaarden, A. Stochastic gene expression in fluctuating environments. Genetics. 167 (1), 523-530 (2004).
- Porter, J. R., Fisher, B. E., Batchelor, E. p53 Pulses Diversify Target Gene Expression Dynamics in an mRNA Half-Life-Dependent Manner and Delineate Co-regulated Target Gene Subnetworks. Cell Syst. 2 (4), 272-282 (2016).
- Lahav, G., Rosenfeld, N., et al. Dynamics of the p53-Mdm2 feedback loop in individual cells. Nat. Genet. 36 (2), 147-150 (2004).
- Levine, A. J., Oren, M. The first 30 years of p53: growing ever more complex. Nat. Rev. Cancer. 9 (10), 749-758 (2009).
- Riley, T., Sontag, E., Chen, P., Levine, A. Transcriptional control of human p53-regulated genes. Nat. Rev. Mol. Cell Biol. 9 (5), 402-412 (2008).
- Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., Madden, T. L. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinf. 13 (1), 134 (2012).
- PCR Technologies: A Technical Guide. , Sigma-Aldrich. St. Louis, MO. (2014).
- Sambrook, J., Russell, D. W. Molecular Cloning: A Laboratory Manual. , Cold Spring Harbor Laboratory Press. Cold Spring Harbor, NY. (2001).
- Batchelor, E., Mock, C. S., Bhan, I., Loewer, A., Lahav, G. Recurrent initiation: a mechanism for triggering p53 pulses in response to DNA damage. Mol. Cell. 30 (3), 277-289 (2008).
- Flow Cytometry: Principles and Applications. , Humana Press. Totowa, NJ. (2007).
- Real-time PCR. , Taylor & Francis. (2006).
- Song, L., Langfelder, P., Horvath, S. Comparison of co-expression measures: mutual information, correlation, and model based indices. BMC Bioinf. 13 (1), 328 (2012).
- Margolin, A. A., Nemenman, I., et al. ARACNE: An Algorithm for the Reconstruction of Gene Regulatory Networks in a Mammalian Cellular Context. BMC Bioinf. 7 (Suppl 1), (2006).
- Haff, L. A. Improved quantitative PCR using nested primers. PCR Methods Appl. 3 (6), 332-337 (1994).
- Hashimshony, T., Senderovich, N., et al. CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq. Genome Biol. 17, 77 (2016).
- Ronander, E., Bengtsson, D. C., Joergensen, L., Jensen, A. T. R., Arnot, D. E. Analysis of Single-cell Gene Transcription by RNA Fluorescent In Situ Hybridization (FISH). J. Vis. Exp. (68), e4073 (2012).
- Raj, A., van den Bogaard, P., Rifkin, S. A., van Oudenaarden, A., Tyagi, S. Imaging individual mRNA molecules using multiple singly labeled probes. Nat. Methods. 5 (10), 877-879 (2008).
- Lubeck, E., Cai, L. Single-cell systems biology by super-resolution imaging and combinatorial labeling. Nat. Methods. 9 (7), 743-748 (2012).
- Battich, N., Stoeger, T., Pelkmans, L. Image-based transcriptomics in thousands of single human cells at single-molecule resolution. Nat. Methods. 10 (11), 1127-1133 (2013).
