منصة موائع جزيئية أولتراهيغ الإنتاجية لتسلسل جينوم خلية واحدة
Summary
يكشف تسلسل خلية واحدة قد تغاير في النظم البيولوجية، ولكنها تفتقر إلى التكنولوجيات الحالية الإنتاجية اللازمة للتنميط العميقة وظيفة وتكوين المجتمع. هنا، يمكننا وصف سير عمل موائع جزيئية للتسلسل > خلية 50,000 جينوم خلية واحدة من تنوع السكان.
Abstract
تقنيات التسلسل شهدت تحولاً من الأكبر إلى خلية واحدة القرار استجابة لفهم دور التغايرية الخلوية متطورة في النظم البيولوجية. ومع ذلك، تسلسل خلية واحدة كبيرة من السكان قد أعاق القيود في تجهيز الجينوم للتسلسل. في هذه الورقة، ويصف لنا طريقة لتسلسل جينوم خلية واحدة (SiC-seq) الذي يستخدم ميكروفلويديكس الحبرية إلى عزل، وتضخيم، والرمز الشريطي الجينوم للخلايا المفردة. تغليف الخلية في ميكروجيلس يسمح بتنقية مجزأة وتاجمينتاتيون من الحمض النووي، بينما اندماج موائع جزيئية أزواج كفاءة كل الجينوم مع رمز شريطي فريد اليغنوكليوتيد خلية واحدة، مما يتيح > 50,000 الخلايا المفردة تكون متسلسلة كل تشغيل. هو ديمولتيبليكسيد البيانات التسلسل بالرمز الشريطي، توليد مجموعات من القراءات التي تنشأ من خلايا مفردة. كأسلوب الإنتاجية العالية والمنخفضة-التحيز في تسلسل خلية واحدة، سيمكن SiC seq طائفة أوسع من الدراسات الجينومية تستهدف السكان خلية المتنوعة.
Introduction
الجينوم بمثابة مخططا للهوية الخلوية والدالة، التي تحتوي على الكائن الحي بأكملها الترميز المحتملة. فهم البيولوجيا الخلوية على مستوى الجينوم يمكن أن تفسر تنوع المظهرية ملاحظتها داخل الخلية غير المتجانسة السكان. هذا التباين الواضح في الأنظمة البيولوجية، وله آثار واسعة النطاق على صحة الإنسان والمرض. على سبيل المثال، ترتبط الجينات نسخ عدد الاختلافات بين الخلايا السرطانية لتطور وانتشار السرطان1،2. في العدوى البكتيرية، الإمراضية جزر موجودة في جزء صغير من الجينوم يمكن نقلها أفقياً، ويؤدي إلى انتشار البكتيريا المقاومة للمضادات الحيوية3،4. هو تحدي الرئيسي في دراسة الجينوم على مستوى الخلايا المفردة كميات منخفضة من الحمض النووي المتاحة، فضلا عن الحاجة إلى تحليل آلاف خلايا العينة التنوع الكامل للأنماط الجينية. لهذه الأسباب، أعاقت القيود الإنتاجية التجريبية فعالية الدراسات خلية واحدة، يتحامل نتائج نحو الخلايا الأكثر وفرة. تقنيات عزل خلية واحدة مثل تدفق الفرز5،6، ملاقط بصرية7، امبيدمينت في معظم المواد الهلامية8، وميكروفلويديكس9 قادرون على تجهيز مئات خلايا لتسلسل؛ ومع ذلك، يمثل هذا جزء صغير فقط من معظم العينات. يسمح أسلوب لتسلسل جينوم خلية واحدة مع إنتاجية أعلى بكثير التنميط أكثر عمقاً وأكثر اكتمالا من السكان الخلية، وبالتالي توضيح دور أوضح التنوع داخل هذه المجتمعات.
ميكروفلويديكس الحبرية يتيح التلاعب الفائق بالخلايا والمواد الكاشفة البيولوجية داخل ملايين مفاعلات بيكوليتير-مقياس. وحتى الآن، ميكرودروبليت قد استخدمت تقنيات لدراسة أنماط التعبير التفاضلية بين خلايا من أنسجة غير متجانسة10،،من1112، عميق تسلسل جزيئات طويلة13،14 ،15، وسلوك الكروماتين إيمونوبريسيبيتيشن التسلسل (الرقائق-seq) تحليلات بشأن خلايا مفردة16. وفي الواقع، ميكرودروبليتس قادرون على عمليات الفائق، ومجزأة، وجعلها قابلة للتطبيقات في علم الجينوم خلية واحدة. تطوير هذه التكنولوجيا يعرض التحديات التكنولوجية الفريدة الخاصة به، ومع ذلك. يجب تفكيك الخلايا وتنقية وتضخيم مع التحيز الحد الأدنى، إلى شكل موحد عينة خلية السكان. بالإضافة إلى ذلك، خلافا للنصوص مرناً بوليادينيلاتيد في خلايا الثدييات، هناك لا فكرة الجزيئية قابلة للمقارنة في الجينوم لتسهيل القبض على الحمض النووي المستهدف. ولهذه الأسباب، قد يصعب تنفيذه في منصات ميكرودروبليت تسلسل جينوم خلية واحدة.
في هذا العمل، نحن نقدم بروتوكول مفصل لنهجنا ذكر سابقا موائع جزيئية خلية واحدة قادرة على تسلسل جينوم عشرات آلاف خلايا في تجربة واحدة17. مغلفة في النطاق ميكرون الهلاميات المائية مع هذه التكنولوجيا، ودعا SiC-seq، الخلايا البكتيرية وتفكيك فردياً، تاجمينتيد، واندمجت مع ميكرودروبليت الذي يحتوي على رمز شريطي اليغنوكليوتيد فريدة من نوعها، وهي تقسم على الحمض النووي للخلية عن طريق تداخل واحد ملحق تفاعل البوليميراز المتسلسل (PCR). الهلاميات المائية مثابة حاويات معزولة فيها عالية الوزن الجزيئي المجينية الحمض النووي هو ستيريكالي المغطى، السماح لجزيئات أصغر مثل المنظفات والأنزيمات الحال للوصول وتنقية الحمض النووي قبل المتوازية18. هذا البروتوكول العمليات > 50,000 الخلايا المفردة في غضون ساعات، أسفر عن المراقب مكتبة جاهزة للتسلسل. عقب التسلسل، ما يلي: ديمولتيبليكسيد وفقا لتسلسلها الباركود خلية واحدة، أدى إلى dataset تتألف من ملايين القراءات، كل مع فهرس خلوية.
Protocol
Representative Results
Discussion
سير العمل موائع جزيئية SiC seq تنتج بيانات تسلسل جينوم خلية واحدة من آلاف خلايا البكتيرية. الرموز الشريطية الرقمية تقسم على مورثات الخلايا مغلفة ميكروجيل تسمح deconvolution في السيليكون خ ع البيانات إلى مجموعات من يقرأ المراقب التي تنشأ من خلية واحدة. تجربة التحكم مع مجتمع الميكروبي لتكوين المعروفة أمر ضروري لتقييم نقاء المجموعات الباركود. جزء كبير من المجموعات المنخفضة-نقاء يشير إلى أن معدل تغليف خلية مرتفع للغاية، أو أن هناك التلوث عبر التجميعية الهامة التي تحدث أثناء خطوات المعالجة موائع جزيئية. وفقا لإحصاءات بواسون، ينبغي تغليف الرموز الشريطية والخلايا بنسبة مستهدفة من الجسيمات 1 لكل 10 قطرات الحد من معدل الأحداث تغليف متعددة إلى أقل من 5 في المائة من جميع قطرات غير فارغة. سعر تغليف أعلى من هذا يزيد معدلات doublets أضعافاً مضاعفة، حتى تحقق نسبة التغليف أثناء عملية دروبماكينج ذات أهمية حاسمة. ينبغي أن يكون المستخدمين الحذر خاصة لتغليف خلايا متعددة في ميكروجيل واحد لأنه لا يمكن أن يقرأ من خلايا مختلفة تقاسم نفس تسلسل الباركود بيوينفورماتيكالي مفصولة. في الحالة يتلقى تلك الخلية 1 2 الرموز الشريطية مختلفة، نقاء المجموعة رمز شريطي، لا تتأثر على الرغم من أن قياسات الوفرة هي منحرفة عند عد بتسلسل الباركود.
التلوث عبر الحبرية قد تنشأ أيضا بسبب ظروف الاندماج دون الحد الأمثل. أثناء عملية ناجحة، يمكن أن كنترولبلي زوج الجهاز الاندماج موائع جزيئية (الشكل 5) 1 الحبرية الباركود مع ميكروجيل 1 وحجم الكاشف PCR. سيؤدي إلى معدلات التدفق المثالي عدم معالجة تجميعية الأقران في نسب غير صحيحة: 1 الباركود يمكن أن يقترن ميكروجيلس 2، على سبيل المثال. المقصود أن تكون تقديرات معدلات تدفق جميع المدرجة في البروتوكول وقد تحتاج إلى تعديلها تبعاً لاختلافات طفيفة في أحجام التجميعية وهندسة الجهاز. المستخدمين الذين لديهم وصول للكاميرات مع قدرات تسجيل عالية السرعة (> إطارات 10,000/s) ينبغي التحقق من عملية الدمج الصحيح الحبرية في البداية وعلى مدى العملية موائع جزيئية. يمكن للمستخدمين دون الحصول على كاميرا عالية السرعة جمع كمية صغيرة من الإخراج المدمجة وقياس أحجام الحبرية تحت مجهر يدوياً. يجب أن يكون حجم الحبرية موحدة: فائض في الباركود غير مدمجة أو قطرات ميكروجيل يشير إلى أنه ينبغي خفض معدلات حقن تبعاً لذلك.
ينبغي اتخاذ عدة احتياطات عامة عند التعامل مع ميكروجيلس وميكرودروبليتس للحفاظ على سلامتهم. ميكروجيلس، على الرغم من أن قوي ميكانيكية، يجب أن يكون بما فيه الكفاية تبرد قبل كسر والغسيل خطوات لضمان جيلاتيون كاملة. ميكروجيلس عدم كروية هي إشارة إلى أن [اغروس] لا تعطي الوقت الكافي لترسيخ. عند غسل ميكروجيلس، تدور التعليق أسفل السرعات المطلوبة لتجنب خسارة المنتج. المائية [اغروس] إنكسار عن كثب مطابقة للمياه وقد يكون من الصعب على انظر في أنبوب22، حيث المستخدمين ينبغي أن تحدد بدقة الحدود هلام السائل قبل تطلع. قطرات الماء في النفط معرضة للتلاحم بحشد قوات ثابتة23 في مختبر القفازات وأنابيب. لهذا السبب، نوصي بتحميل الحقن حقن الحبرية بأيديهم العارية وعلاج كافة الأسطر مع بندقية مكافحة ساكنة قبل فتيلة مضخة حقن. يمكن إزالة قطرات كبيرة ملتئم بالتناوب المستحلبات في حقنه ببطء ويدويًا يسفط قطرات أكبر، التي تتراكم بالقرب من أعلى بسبب قوتها ازدهار أكبر.
Seq سيك هو أول تقنية لإظهار تسلسل جينوم خلية واحدة > 50,000 الخلايا البكتيرية. هذا البرنامج يوفر مزايا كبيرة في إنتاجية أكثر من النهج القائمة وتمكن أعمق لأخذ عينات من المجتمعات الميكروبية غير متجانسة. وحتى الآن استخدمت التكنولوجيات موائع جزيئية لتسلسل جينوم خلية واحدة ميكروتشامبيرس9 وميكروويلس24 لعزل الخلية والتضخيم، ولكن مع الوسيطة في نطاق فقط من عشرات إلى مئات خلايا. تدفق فرز الخلايا المفردة في ويلبلاتيس5،6 يتطلب لا الأجهزة المتخصصة موائع جزيئية ولكن يمتلك إنتاجية منخفضة على نحو مماثل. نظراً لأن عينات التربة والمياه من البيئة قد يشيع التنوعات ألفا > 1,000 في هذه الأنواع من المستوى25،26، SiC seq مفيد جداً بحكم قدرته على العينة عدد أكبر بكثير من الكائنات الحية. يكون سير العمل SiC seq قابلة للتكيف للمدخلات الخلية من مختبر الثقافة، والبيئة الطبيعية، أو مضيف معيشة. عينة خلية تحتاج فقط في وقف مائي وخالية من جزيئات كبيرة (> 10 ميكرون) تكون مناسبة لتغليف موائع جزيئية. على سبيل المثال، الأسلوب قد طبق سابقا على عينة من مياه البحر باستخدام سلسلة من المياه والصرف الصحي وتصفية خطوات لعملية ما قبل الخلايا قبل تغليف17.
البروتوكول SiC seq يولد كمية ضئيلة نسبيا من تسلسل البيانات من كل خلية مفردة، وقد لا تكون مناسبة لجميع التطبيقات. تتطلب بعض خوارزميات المعلوماتية الحيوية مثل الجمعية الجينوم حيثياته أو متغير واحد النوكليوتيدات (الهولندية) تطلق على أعماق تغطية أكبر للعمل بفعالية. بدلاً من ذلك، يمكن أن تكون مجموعات الرمز الشريطي متفاوت المسافات في السيليكون التصنيفية binning أساليب27 حيث أنه يمكن تطبيق الخوارزميات في مجموعات أكبر من ما يلي. يجوز كفاءة سير العمل SiC seq المتوازية الشاملة منخفضة نسبيا أيضا تحديات في الحالات التي يكون فيها توافر العينة الإدخال المنخفضة. Seq سيك ويعتمد على خطوة تغليف توزيع Poisson باركود، ولذلك تلقي رمز شريطي جزيئي حوالي 10% الخلايا وتتضخم خلال الخطوة إعداد مكتبة النهائي. بينما هذا يماثل الأخرى المستندة إلى ميكرودروبليت المتوازية مخططات10، المستخدمين الذين يعملون مع عينات الخلايا الثمينة قد يجدون صعوبة في تحقيق العائد مكتبة كافية للتسلسل وقد تحتاج إلى زيادة عدد دورات PCR في النهائي الخطوة التضخيم. حل محتمل آخر للمستخدمين ذوي الخبرة موائع جزيئية لفرز قطرات الباركود إيجابية بعد الخطوة PCR الرقمية، مما يجعل كفاءة المتوازية الشاملة > 85%28.
اتجاه مستقبل محتمل للتكنولوجيا SiC seq هو تكييف سير العمل للاستخدام مع خلايا الثدييات، مما يمهد الطريق للجديد من الدراسات السريرية في خلية واحدة. على سبيل مثال، خلايا تحليل التباين رقم نسخة بين السرطان واحد أيار زيادة فهمنا لدور عدم التجانس في أمراض السرطان2. وبدلاً من ذلك، سيمكن دمج SiC seq مع الأساليب القائمة التحقيق وإثراء تسلسل الحمض النووي لفائدة29 تسلسل وحيدة الخلية المستهدفة من الفئات السكانية الفرعية أو السلالات النادرة من الخلايا. مع العينات البيئية، يمكن استهداف الجينات من ضمن مسار التمثيل الغذائي معروفة وتحليلها سياقياً جنبا إلى جنب مع جينات للتعرف على رواية المجينية الجزر المجاورة. من داخل بيئة الإنسان مضيف، عينات عيار منخفض البكتيريا المسببة للأمراض يمكن أن تكون معزولة ومتسلسلة على مستوى خلية واحدة لدراسة المزيد عن كثب أصولهم أوضح من ضراوة.
Disclosures
The authors have nothing to disclose.
Acknowledgements
هذا العمل كان تدعمها “المؤسسة الوطنية للعلوم” من خلال جائزة مهنة (رقم المنحة DBI 1253293)؛ المعاهد الوطنية للصحة (NIH) (أرقام المنح HG007233-01، R01-EB019453-01، 1R21HG007233، DP2-AR068129-01، R01-HG008978)؛ و “الدفاع المتقدم البحث مشاريع الوكالة يعيشون المسابك البرنامج” (العقد رقم HR0011-12-ج-0065، N66001-12-ج-4211، HR0011-12-ج-0066).
Materials
| 3" silicon wafers, P type, virgin test grade | University Wafers | 447 | |
| SU-8 3025 photoresist | Microchem | 17030192 | |
| Spin coater | Specialty Coating Systems | G3P-8 | |
| Photomasks | CadArt Servcies | (custom) | See Supplemental Files for mask designs |
| PGMEA developer | Sigma-Aldrich | 484431 | |
| Isopropanol | Sigma-Aldrich | 109827 | |
| Sylgard 184 silicone elastomer kit | Krayden | 4019862 | |
| Degassing chamber | Bel-Art | 42025 | |
| 0.75 mm biopsy punch | World Precision Instruments | 504529 | |
| Glass microscope slides (75 mm x 50 mm) | Corning | 294775X50 | |
| Aquapel (hydrophobic glass treatment) | Pittsburgh Glass Works | 47100 | |
| PE-2 polyethylene tubing | Scientific Commodities | B31695-PE/2 | |
| 1 mL syringes | BD | 309628 | |
| 27 gauge needles | BD | 305109 | |
| Syringe pump | New Era Pump Systems | NE-501 | |
| Novec HFE-7500 fluorinated oil (HFE) | 3M | 98-0212-2928-5 | |
| FC-40 fluorinated oil | Sigma-Aldrich | F9755 | |
| PEG-PFPE surfactant | Ran Biotechnologies | 008-FluoroSurfactant | |
| Space heater | Lasko | CD09250 | |
| Agarose, low gelling temperature | Sigma-Aldrich | a9414 | |
| TE (10X) | Rockland | mb-007 | |
| PBS 1X, pH 7.4 | E&K Scientific Products | EK-65083 | |
| OptiPrep (density gradient medium) | Sigma-Aldrich | d1556 | |
| 1H,1H,2H-Perfluoro-1-Octanol (PFO) | Sigma-Aldrich | 370533 | |
| Span 80 (sorbitane monooleate) | Sigma-Aldrich | s6760 | |
| Hexane | Sigma-Aldrich | 139386 | |
| Tween 20 (polysorbate 20) | Sigma-Aldrich | p2287 | |
| Lysozyme Type IV | MP Biomedicals | 195303 | |
| Mutanolysin | Sigma-Aldrich | M9901 | |
| Zymolyase (yeast lytic enzyme) | Zymo Research | e1004 | |
| Lysostaphin | Sigma-Aldrich | L7386 | |
| Sodium chloride | Sigma-Aldrich | S9888 | |
| EDTA | Sigma-Aldrich | E6758 | |
| Tris-HCl, pH 7.5, 1M | Invitrogen | 15567-027 | |
| Dithiothreitol (DTT) | Teknova | d9750 | |
| Lithium dodecyl sulfate | Sigma-Aldrich | L9781 | |
| Proteinase K | New England Biosciences | P8107S | |
| Ethanol, 200 Proof (100%) | Koptec | V1001 | |
| SYBR Green I (nucleic acid stain) | Invitrogen | S7563 | |
| PEG 6k | Sigma-Aldrich | 81260 | |
| Triton X-100 (octylphenol ethoxylate) | Sigma-Aldrich | t8787 | |
| Nextera DNA Library Prep Kit | Illumina | FC-121-1030 | |
| Phusion Hot Start Flex Master Mix (High-Fidelity Hot Start Master Mix) | New England Biosciences | m05365 | |
| Platinum Multiplex PCR Master Mix (Taq Master Mix) | Applied Biosystems | 4464263 | |
| Warmstart 2.0 Bst Polymerase (isothermal polymerase) | New England Biosciences | m0538m | |
| NT buffer from Nextera XT kit (neutralization buffer) | Illumina | FC-131-1024 | |
| Cold cathode fluorescent inverter | (custom) | (custom) | |
| DC power supply | Mastech | HY1503D | |
| Zerostat 3 anti-static gun | Milty | 5036694022153 | |
| 3D-printed centrifuge syringe holder | (custom) | (custom) | See Supplemental Files for 3D print file |
| Zymo DNA Clean & Concentrator-5 | Zymo Research | D4003 |
References
- Navin, N., et al. Tumour evolution inferred by single-cell sequencing. Nature. 472 (7341), 90-94 (2011).
- Ni, X., et al. Reproducible copy number variation patterns among single circulating tumor cells of lung cancer patients. Proceedings of the National Academy of Sciences of the United States of America. 110 (52), 21083-21088 (2013).
- Schmidt, H., Hensel, M. Pathogenicity islands in bacterial pathogenesis. Clinical Microbiology Reviews. 17, 14-56 (2004).
- Martínez, J. L., Baquero, F. Interactions among strategies associated with bacterial infection: pathogenicity epidemicity, and antibiotic resistance. Clinical Microbiology Reviews. 15 (4), 647-679 (2002).
- Rinke, C., et al. Obtaining genomes from uncultivated environmental microorganisms using FACS-based single-cell genomics. Nature Protocols. 9 (5), 1038-1048 (2014).
- Rinke, C., et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 499 (7459), 431-437 (2013).
- Zhang, H., Liu, K. K. Optical tweezers for single cells. Journal of the Royal Society Interface. 5 (24), 671-690 (2008).
- Xu, L., Brito, I. L., Alm, E. J., Blainey, P. C. Virtual microfluidics for digital quantification and single-cell sequencing. Nature Methods. 13 (9), 759-762 (2016).
- Gawad, C., Koh, W., Quake, S. R. Dissecting the clonal origins of childhood acute lymphoblastic leukemia by single-cell genomics. Proceedings of the National Academy of Sciences of the United States of America. 111 (50), 17947-17952 (2014).
- Macosko, E. Z., et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 161 (5), 1202-1214 (2015).
- Klein, A. M., et al. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell. 161 (5), 1187-1201 (2015).
- Rotem, A., et al. High-throughput single-cell labeling (Hi-SCL) for RNA-Seq using drop-based microfluidics. PLoS One. 10 (5), 1-14 (2015).
- Amini, S., et al. Haplotype-resolved whole-genome sequencing by contiguity-preserving transposition and combinatorial indexing. Nature Reviews Genetics. 46 (12), 1343-1349 (2014).
- Zheng, G. X. Y., et al. Haplotyping germline and cancer genomes with high-throughput linked-read sequencing. Nature Biotechnology. 34 (3), 303-311 (2016).
- Lan, F., Haliburton, J. R., Yuan, A., Abate, A. R. Droplet barcoding for massively parallel single-molecule deep sequencing. Nature Communications. 7, 11784 (2016).
- Rotem, A., et al. Single-cell ChIP-seq reveals cell subpopulations defined by chromatin state. Nature Biotechnology. 33 (11), 1165-1172 (2015).
- Lan, F., Demaree, B., Ahmed, N., Abate, A. R. Single-cell genome sequencing at ultra-high-throughput with microfluidic droplet barcoding. Nature Biotechnology. 35 (7), 640-646 (2017).
- Novak, R., et al. Single-cell multiplex gene detection and sequencing with microfluidically generated agarose emulsions. Angewandte Chemie Internation Edition. 50 (2), 390-395 (2011).
- Gill, C., Van De Wijgert, J. H. H. M., Blow, F., Darby, A. C. Evaluation of lysis methods for the extraction of bacterial DNA for analysis of the vaginal microbiota. PLoS One. 11 (9), 1-16 (2016).
- Picelli, S., et al. Tn5 transposase and tagmentation procedures for massively scaled sequencing projects. Genome Research. 24 (12), 2033-2040 (2014).
- Abate, A. R., Chen, C. H., Agresti, J. J., Weitz, D. A. Beating Poisson encapsulation statistics using close-packed ordering. Lab on a Chip. 9 (18), 2628 (2009).
- Jain, A., Yang, A. H. J., Erickson, D. Gel-based optical waveguides with live cell encapsulation and integrated microfluidics. Optic Letters. 37 (9), 1472 (2012).
- Karbaschi, M., Shahi, P., Abate, A. R. Rapid chemical-free breaking of microfluidic emulsions with a hand-held antistatic gun. Biomicrofluidics. 11 (4), 1-6 (2017).
- Gole, J., et al. Massively parallel polymerase cloning and genome sequencing of single cells using nanoliter microwells. Nature Biotechnology. 31 (12), 1126-1132 (2013).
- Chao, Y., et al. Metagenomic analysis reveals significant changes of microbial compositions and protective functions during drinking water treatment. Scientific Reports. 3 (1), 3550 (2013).
- Fierer, N., et al. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proceedings of the National Academy of Sciences of the United States of America. 109 (52), 21390-21395 (2012).
- Mande, S. S., Mohammed, M. H., Ghosh, T. S. Classification of metagenomic sequences: methods and challenges. Briefings in Bioinformatics. 13 (6), 669-681 (2012).
- Eastburn, D. J., et al. Microfluidic droplet enrichment for targeted sequencing. Nucleic Acids Research. 43 (13), e86 (2015).
- Clark, I. C., Abate, A. R. Finding a helix in a haystack: nucleic acid cytometry with droplet microfluidics. Lab on a Chip. 17 (12), 2032-2045 (2017).
