تصف هذه الدراسة سير العمل لإدارة بيانات مثيلة الحمض النووي التي تم الحصول عليها بواسطة تقنيات microarray . يوضح البروتوكول الخطوات من إعداد العينات إلى تحليل البيانات. يتم وصف جميع الإجراءات بالتفصيل ، ويعرض الفيديو الخطوات المهمة.
مقالة منهجية
* These authors contributed equally
تصف هذه الدراسة سير العمل لإدارة بيانات مثيلة الحمض النووي التي تم الحصول عليها بواسطة تقنيات microarray . يوضح البروتوكول الخطوات من إعداد العينات إلى تحليل البيانات. يتم وصف جميع الإجراءات بالتفصيل ، ويعرض الفيديو الخطوات المهمة.
ترتبط السمنة ارتباطا مباشرا بنمط الحياة وقد ارتبطت بتغيرات مثيلة الحمض النووي التي قد تسبب تغيرات في عمليات تكوين الدهون وتخزين الدهون التي تسهم في تطور المرض. نحن نظهر بروتوكولا كاملا من الاختيار إلى تحليل البيانات اللاجينية للمرضى الذين يعانون من السمنة أو بدونها. تم اختبار جميع الخطوات من البروتوكول والتحقق من صحتها في دراسة تجريبية. شاركت 32 امرأة في الدراسة ، حيث تم تصنيف 15 شخصا يعانون من السمنة وفقا لمؤشر كتلة الجسم (BMI) (45.1 ± 5.4 كجم / م 2) ؛ وتم تصنيف 17 فردا بدون سمنة وفقا لمؤشر كتلة الجسم (22.6 ± 1.8 كجم / م 2). في المجموعة المصابة بالسمنة ، تم تحديد 564 موقعا من مواقع CpG المتعلقة بكتلة الدهون من خلال تحليل الانحدار الخطي. كانت مواقع CpG في مناطق المروجين. وجد التحليل التفاضلي 470 CpGs hypomethylated و 94 موقعا مفرط الميثيل في الأفراد الذين يعانون من السمنة. كانت المسارات الأكثر تخصيبا بنقص ميثيل في RUNX ، وإشارات WNT ، والاستجابة لنقص الأكسجة. كانت مسارات فرط الميثيل مرتبطة بإفراز الأنسولين ، وإشارات الجلوكاجون ، و Ca2 +. نستنتج أن البروتوكول حدد بشكل فعال أنماط مثيلة الحمض النووي ومثيلة الحمض النووي المرتبطة بالسمات. يمكن أن ترتبط هذه الأنماط بالتعبير الجيني المتغير ، مما يؤثر على تكوين الدهون وتخزين الدهون. أكدت نتائجنا أن نمط الحياة البدينة يمكن أن يعزز التغيرات اللاجينية في الحمض النووي البشري.
تم استخدام تقنيات omics واسعة النطاق بشكل متزايد في دراسات الأمراض المزمنة. ميزة مثيرة للاهتمام لهذه الأساليب هي توافر كمية كبيرة من البيانات التي تم إنشاؤها للمجتمع العلمي. لذلك ، نشأ طلب على توحيد البروتوكولات للسماح بالمقارنة الفنية بين الدراسات. تقترح هذه الدراسة توحيد بروتوكول للحصول على بيانات مثيلة الحمض النووي وتحليلها ، باستخدام دراسة تجريبية كمثال مطبق.
يسود إنفاق الطاقة السلبي في أنماط الحياة البشرية الحديثة ، مما يؤدي إلى تراكم مفرط للأنسجة الدهنية ، وبالتالي تطور السمنة¹. أدت العديد من العوامل إلى زيادة معدلات السمنة ، مثل المستقرة ، والوجبات الغذائية عالية السعرات الحرارية ، والروتين المجهد. قدرت منظمة الصحة العالمية (WHO) أن 1.9 مليار بالغ يعانون من السمنة المفرطة في عام 2016 ، مما يعني أن أكثر من 20٪ من سكان العالم لديهم أكثر من 30 كجم / م 2 BMI2. كشف آخر تحديث لعام 2018 أن انتشار السمنة في الولايات المتحدة الأمريكية كان أعلى من 42٪ 3.
علم الوراثة اللاجينية هو التكيف الهيكلي للمناطق الكروموسومية لتسجيل حالات النشاط المتغيرة أو الإشارة إليها أو إدامتها4. مثيلة الحمض النووي هي تغيير كيميائي عكسي في مواقع ثنائي النوكليوتيدات السيتوزين - غوانوزين (مواقع CpG) ، وتشكيل 5-methylcytosine-pG (5mCpG). يمكنه تعديل التعبير الجيني عن طريق تنظيم وصول آلية النسخ إلى الحمض النووي5،6،7،8. في هذا السياق، من الضروري فهم مواقع CpG المرتبطة بالسمات المرتبطة بالسمنة9. يمكن للعديد من العوامل دعم أو منع مثيلة الحمض النووي الخاصة بالموقع. ويمكن للإنزيمات الضرورية لهذه العملية، مثل ميثيل ترانسفيراز الحمض النووي10 (DNTMs) وعشرة إلى أحد عشر عملية نقل (TETs)، أن تعزز مثيلة الحمض النووي أو إزالة الميثيل في ظل التعرض البيئي11.
بالنظر إلى الاهتمام المتزايد بدراسات مثيلة الحمض النووي على مدى السنوات الماضية ، كان اختيار استراتيجية التحليل الأكثر ملاءمة للإجابة بدقة على كل سؤال مصدر قلق أساسي للباحثين12،13،14. تعد مصفوفة مثيلة الحمض النووي 450K الطريقة الأكثر شيوعا ، وتستخدم في أكثر من 360 منشورا14 لتحديد ملف تعريف مثيلة الحمض النووي. يمكنه تحديد مثيلة ما يصل إلى 485000 CpGs الموجودة في 99٪ من الجينات المعروفة15. ومع ذلك ، تم إيقاف هذه الصفيف واستبدالها ب EPIC ، والتي تغطي مواقع 850,000 CpG. يمكن تطبيق البروتوكول الحالي على كل من 450K و EPIC16 و 17 و 18.
يتم عرض البروتوكول خطوة بخطوة في الشكل 1 ويشمل الخطوات التالية: اختيار السكان ، وأخذ العينات ، وإعداد التجارب ، وخط أنابيب مثيلة الحمض النووي ، وتحليل المعلوماتية الحيوية. يتم عرض دراسة تجريبية أجريت في مختبرنا هنا لتوضيح خطوات البروتوكول المقترح.

الشكل 1: مخطط البروتوكول المعروض. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
وافقت لجنة الأخلاقيات في مستشفى جامعة ريبيراو بريتو كلية الطب بجامعة ساو باولو (HCRP-USP) على الدراسة (CAAE:14275319.7.0000.5440). ووقع المشاركون على استمارة الموافقة، وأجريت جميع الإجراءات عقب إعلان هلسنكي.
1. السكان وأخذ العينات
2. الأنثروبومتري وتكوين الجسم
الجدول 1: الخصائص السكانية. كانت جميع المتغيرات بارامترية (p > 0.05 ، Shapiro Wilk) ، وتم تقييم الاختلافات باستخدام اختبار t مستقل ، حيث اعتبر p < 0.05 معنويا. * ص < 0.0001. يرجى النقر هنا لتنزيل هذا الجدول.
3. جمع المواد البيولوجية
4. استخراج الحمض النووي
5. التحضير قبل تحليل المثيلة
6. خط أنابيب مثيلة الحمض النووي
ملاحظة: تنقسم تجربة مثيلة الحمض النووي إلى 4 أيام (الشكل 1). اتبع توصيات الشركة المصنعة ومواصفاتها للحصول على نتائج دقيقة.
7. تحليل المعلوماتية الحيوية
ملاحظة: من الضروري تغيير بعض سمات نوع المصفوفة الدقيقة (على سبيل المثال، يجب استبدال نوع الصفيف = "450k" بنوع الصفيف = "EPIC"). يصف مؤلف خط أنابيب ChAMP بالتفصيل جميع الحقول التي يجب تعديلها لتحليل المصفوفات على صفحة الموصلات الحيوية للحزمة30.
بعد استخدام خط أنابيب ChAMP ، تم النظر في 409,887 مسبارا في التحليل بعد جميع المرشحات (CpGs غير المؤهلة ، تحقيقات غير CG ، CpG بالقرب من SNPs ، مواقع متعددة الضربات ، و CpGs المتعلقة ب XY) ؛ ويمثل الشكل 2 مخططا.

الشكل 2: خط أنابيب تحليل المعلوماتية الحيوية. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
وعلاوة على ذلك، كشفت خطة الكثافة أن جميع العينات لها كثافات مماثلة لتوزيع بيتا تتعلق بخطوات مراقبة الجودة. يقيم هذا التحليل توزيع قيم بيتا ويشير إلى ما إذا كانت هناك أي عينات تحتاج إلى استبعادها. في هذه الدفعة، لم يتم استبعاد أي عينات بناء على هذا التحليل.

الشكل 3: مخطط كثافة القيم التجريبية الذي تم الحصول عليه باستخدام حزمة ChAMP. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
تم استخدام تحليل تحلل القيمة المفردة (SVD) للتحقق من المكونات الرئيسية التي أثرت بشكل كبير على تباين بيانات مثيلة الحمض النووي. وكشفت البيانات أن مؤشر كتلة الجسم وWC (p < 1e10-5) وFM (p < 0.05) كان لها تأثير كبير على تباين البيانات (الشكل 4).

الشكل 4: تحليل تحلل القيمة المفردة. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
كشف تقدير نوع الخلية أن كلا من الخلايا القاتلة الطبيعية (NK) والخلايا البائية كانت أعلى في النساء البدينات (الشكل 5).

الشكل 5: كسور الخلايا المقدرة بطريقة هاوسمان. غران: الخلايا المحببة. CD4T: مساعد الخلايا التائية ، الخلايا الليمفاوية. CD8T: الخلايا التائية السامة للخلايا ، الخلايا الليمفاوية. أحادي: مونوسيت. الخلية البائية: الخلايا الليمفاوية البائية. NK: الخلايا القاتلة الطبيعية ، الخلايا الليمفاوية. * ص < 0.05. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
اختلفت مستويات مثيلة الحمض النووي بين النساء البدينات وغير البدينات قبل وبعد تصحيح نوع الخلية. قبل تصحيح بيانات مثيلة الحمض النووي لأنواع الخلايا ، لوحظ وجود 43,463 موقعا ميثيليا تفاضليا (DMPs) ، وظل 3,329 CpGs كبيرا بعد ذلك. وكان 445 موقعا من مواقع CpG في المناطق بين الجينات (IGR)، و 2,884 موقعا في المناطق الوراثية، ومعظمها في منطقة المروج (n = 1,438). كان التوزيع على طول جميع المناطق TSS1500 (n = 612) ، TSS200 (n = 826) ، 5'UTR (n = 390) ، الإكسون الأول (n = 273) ، الجسم (n = 724) ، و 3'UTR (n = 59). وبالنظر إلى قيم Δβ <-0.05 و >0.05 ، كان 162 CpGs هيبوميثيل و 576 مفرط الميثيل في البدناء مقارنة بالأفراد غير البدناء (الشكل 6). والبيانات متاحة في قاعدة بيانات توقعات البيئة العالمية تحت رمز السجل GSE166611.

الشكل 6: المواضع الميثيلية التفاضلية. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
من الممكن تقييم المثيلة المختلفة بين المجموعات والعثور على مواقع CpG المرتبطة بسمات محددة في دراسات المثيلة. بالنسبة لدراسات كتلة الدهون ، تم العثور على مواقع 13,222 CpG. وارتبط 6,159 CpGs في المنطقة المروجة بكتلة الدهون، مع 470 فرط ميثيل و 94 هيبوميثيل (الشكل 7)، وتم إثراء الجينات المعنية (الجدول 2).

الشكل 7: المناطق المروجة للجينات الهيبوميثيل وفرط ميثيل ، المرتبطة بشكل مستقل بكتلة الدهون.
الجدول 2: الإثراء الوظيفي للجينات من مواقع CpG الميثيلية بشكل مختلف. يرجى النقر هنا لتنزيل هذا الجدول.
المواد التكميلية 1: يرجى النقر هنا لتنزيل هذا الملف.
صفائف مثيلة الحمض النووي هي الطرق الأكثر استخداما للوصول إلى مثيلة الحمض النووي بسبب نسبة التكلفة إلى العائد14. وصفت هذه الدراسة بروتوكولا مفصلا باستخدام منصة microarray متاحة تجاريا لتقييم مثيلة الحمض النووي في دراسة تجريبية أجريت في مجموعة برازيلية. وأكدت النتائج التي تم الحصول عليها من الدراسة التجريبية فعالية البروتوكول. ويبين الشكل 3 قابلية العينة للمقارنة والتحويل الكامل للبيسلفيت32.
كخطوة لمراقبة الجودة ، أوصت خوارزمية ChAMP باستبعاد مواقع CpGs أثناء عملية التصفية. الهدف من استبعاد التحقيقات هو تحسين تحليل البيانات والقضاء على التحيز. تمت إزالة CpGs منخفضة الجودة (قيم p أقل من 0.05) للقضاء على الضوضاء التجريبية في مجموعة البيانات. ظلت الأهداف مرتازة في تحليل مخطط الكثافة. وصف Zhou33 أهمية تصفية CpGs بالقرب من SNPs لتجنب عدم التطابق ، وسوء تفسير مثيلة السيتوسيبين متعدد الأشكال ، والتسبب في تبديل لون تصميم المسبار من النوع الأول34. أيضا ، نظرا لأن كروموسومات XY تتأثر بشكل مختلف بالبصمة ، فقد عزز Heiss و Just35 أهمية تصفية هذه المجسات لأنه ، في الإناث ، قد تكون مشاكل التهجين عوامل مربكة 35.
يعتبر تاريخ انتهاء صلاحية DMAPs ، وتاريخ افتتاح الفورماميد ، والجودة التحليلية للإيثانول المطلق ، وإجمالي عدد الكريات البيض خطوات حاسمة في البروتوكول.
علاوة على ذلك ، وفقا لملاحظاتنا ، فإن تقدير نوع الخلية ضروري في إجراء تحليل المعلوماتية الحيوية. تقوم طريقة Houseman بإجراء تقدير نوع الخلية كما هو موضح في دراسة تيان 30. تعتمد هذه الطريقة على 473 موقعا محددا ل CpG يمكنها التنبؤ بالنسب المئوية لأهم أنواع الخلايا ، مثل الخلايا المحببة والخلايا الوحيدة والخلايا البائية والخلايا التائية36. استخدمنا الوظيفة الموصى بها "myRefbase" من حزمة ChAMP. بعد التقدير ، تقوم خوارزمية ChAMP بضبط قيم بيتا وإزالة هذا التحيز من مجموعة البيانات. هذه الخطوة حاسمة في الدراسات التي تركز على السمنة لأن هذه المجموعة لديها اختلاف كبير في خلايا الدم البيضاء بسبب حالتها الالتهابية المزمنة.
لقد غيرنا فقط خريطة الغطاء الأصلية لختم PCR الشائع فيما يتعلق بتعديل الطريقة واستكشاف الأخطاء وإصلاحها. بعد كل عملية طرد مركزي ، تم تغيير الختم لآخر جديد. لم نتمكن من استخدام الختم الحراري القياسي وتكييفه باستخدام رقائق الألومنيوم حول اللوحة.
على الرغم من أن الفحوصات التجارية قد اعتبرت معيارا ذهبيا للدراسات اللاجينية ، إلا أن أحد قيود البروتوكول يمكن أن يكون خصوصية الكواشف والمعدات من علامة تجارية فريدة من نوعها 37،38،39،40. وثمة قيد آخر يتمثل في عدم وجود مؤشرات تسمح بتحديد التقدم الصحيح للتجربة(41).
يمثل توحيد هذا البروتوكول دليلا رائعا للبحوث اللاجينية ، مما يقلل من الأخطاء البشرية أثناء العملية ويسمح بتحليل البيانات بنجاح وقابلية المقارنة بين الدراسات المختلفة.
وفقا لنتائجنا، فإن تجارب مثيلة الحمض النووي مناسبة للدراسات التي تقارن الأفراد الذين يعانون من السمنة وبدونها43. كما أن تحليل المعلوماتية الحيوية المقترح يوفر بيانات عالية الجودة ويمكن النظر فيه في الدراسات واسعة النطاق.
باستخدام تحليل SVD ، حددنا أن السمات المرتبطة بالسمنة (مؤشر كتلة الجسم ، WC ، و FM) أثرت على التباين في بيانات مثيلة الحمض النووي. ونتيجة لذلك، يشير تقدير نوع الخلية إلى أن كلا من الخلايا القاتلة الطبيعية (NK) والخلايا البائية كانت أعلى لدى النساء المصابات بالسمنة منها في النساء غير المصابات بالسمنة (الشكل 5). ويمكن تفسير ارتفاع عدد هذه الخلايا بالحالة الالتهابية المنخفضة الدرجة لهؤلاء الأفراد44. لاحظنا أن المرضى الذين يعانون من السمنة لديهم CpGs نقص وفرط ميثيل في المناطق المروجة للجينات المرتبطة بكتلة الدهون. كانت معظم المواقع هيبوميثيل ، والتي يمكن أن تكون مرتبطة بالزيادة الطبيعية في مستويات أنواع الأكسجين التفاعلية (ROS) في هؤلاء الأفراد. قد تعزز حالة الإجهاد التأكسدي هذه الغوانين perturbance في موقع ثنائي النوكليوتيد ، وتشكل 8-hydroxy-2'-deoxyguanosine (8-OHdG) ، مما يؤدي إلى موقع ثنائي النوكليوتيد 5mCp-8-OHdG ، ويسبب تجنيد إنزيمات. كل هذه الأحداث يمكن أن تكون مسؤولة عن تعزيز نقص ميثيل الحمض النووي وفرط المثيلة من خلال آليات عمل مختلفة45.
بالإضافة إلى ذلك ، يبدو أن معدل تكوين الدهون يزداد لدى الأفراد الذين يعانون من السمنة ، مع ما يقرب من 10 ٪ من الخلايا الجديدة إلى الخلايا القديمة46،47. يمكن للمساهمات اللاجينية، التي تؤكد على البيئة المسببة للسمنة، أن تغير معدلات انتشار الخلايا وتمايزها، مما يفضل تطور كتلة الدهون48. يمكن أن تؤثر التغيرات اللاجينية أيضا على البرامج الشحمية، مما يسهل أو يقيد تطورها. عوامل النسخ الأولية (PPARγ أو C / EBPα) أو تجميع مجمعات البروتينات المتعددة، المتمركزة في مناطق المروج النهائية التي تعمل عن طريق تضمين أو استبعاد الإنزيمات المعدلة فوق الجينية، تنظم التعبير الجيني من خلال فرط أو نقص الميثيل45. تم وصف مسار PPARγ سابقا لتغيير مسار WNT ، الذي كان يحتوي على جينات غنية في هذه الدراسة. على الرغم من أنه لا يزال من غير المعروف كيف تحدث إشارات WNT أثناء تكوين الدهون ، فقد ذكرت الدراسات الحديثة أنه قد يكون لها أدوار أساسية في استقلاب الخلايا الشحمية ، خاصة في ظل الظروف المسببة للسمنة49.
ليس لدى مؤلفي هذه المقالة أي تضارب في المصالح وليس لديهم إفصاح مالي للإعلان عنه.
نود أن نشكر يوان تيان ، دكتوراه (tian.yuan@ucl.ac.uk) لكونه متاحا للإجابة على جميع الشكوك حول حزمة ChAMP. ونشكر أيضا غيليرمي تيليس، السيدة على مساهمته في المسائل التقنية والعلمية الواردة في هذه الورقة؛ قدم اعتبارات مهمة فيما يتعلق بعلم الوراثة اللاجينية وتقنيات التقاط الفيديو وتنسيقه (guilherme.telles@usp.br). تمويل المواد الاستهلاكية: مؤسسة ساو باولو للأبحاث (FAPESP) (#2018/24069-3) والمجلس الوطني للتنمية العلمية والتكنولوجية (CNPq: #408292/2018-0). التمويل الشخصي: (FAPESP: #2014/16740-6) وبرنامج التميز الأكاديمي من التنسيق لتطوير موظفي التعليم العالي (CAPES:88882.180020/2018-01). سيتم إتاحة البيانات للجمهور ومجانا دون قيود. مراسلات العنوان إلى NYN (البريد الإلكتروني: nataliayumi@usp.br) أو CBN (البريد الإلكتروني: carla@fmrp.usp.br).
| الاسم | الشركة | رقم فهرسي | التعليقات |
|---|---|---|---|
| الإيثانول المطلق | JT Baker | B5924-03 | |
| Agarose gel | Kasvi | K9-9100 | |
| المعاوقة الحيوية الكهربائية | Quantum BIA 450 Q - نظام | ||
| حمض الإيثيلين ديامين تترا أسيتيك (EDTA) Corning | 46-000-CI | ||
| EZ DNA مجموعة مثيلة الحمض | النووي الذهبية ZymoResearch ، إيرفين ، كاليفورنيا ، الولايات المتحدة الأمريكية | D5001 | |
| فورماميد | سيجما | F9037 | |
| FMS & [مدش] ؛ حل التجزئة | Illumina | 11203428 | الكواشف الموردة |
| HumanMethylation450 BeadChip | Illumina | ||
| Maxwell Instrument | Promega ، البرازيل | AS4500 | |
| MA1 & [مدش] ؛ تضخيم متعدد العينات 1 مزيج | Illumina | 11202880 | الكواشف الموردة |
| MicroAmp فيلم لاصق بصري | Thermo Fisher Scientific | 201703982 | |
| MSM & [مدش] ؛ متعدد العينات تضخيم رئيسي ميكس | Illumina | 11203410 | الكواشف الموردة |
| NaOH | F. MAIA | 114700 | |
| PB1 &[مدش] ؛ الكاشف المستخدم لإعداد رقائق الخرز للتهجين | Illumina | 11291245 | الكواشف الموردة |
| PB2 و [مدش] ؛ عازلة الترطيب المستخدمة أثناء التهجين | Illumina | 11191130 | الكواشف الموردة |
| 2-بروبانول | Emsure | 10،96،34،01،000 | |
| RA1 &[مدش] ؛ إعادة التعليق والتهجين والغسيل الحل | Illumina | 11292441 | الكواشف الموردة |
| RPM & [مدش] ؛ مزيج التمهيدي العشوائي | Illumina | 15010230 | الكواشف الموردة |
| STM & [مدش] ؛ متفوقة ثنائية اللون الرئيسية مزيج | Illumina | 11288046 | الكواشف الموردة |
| TEM & [مدش] ؛ لوحان تمديد رئيسي ميكس | Illumina | 11208309 | الكواشف الموردة |
| فائقة النقاء EDTA | Invitrogen | 155576-028 | |
| 96 بئر لوحة تفاعل مع الباركود (0.1 مل) | ByoSystems | 4346906 | |
| لوحة تفاعل 96 بئر مع الباركود (0.8 مل) Thermo | Fisher Scientific | AB-0859 | |
| XC1 & [مدش] ؛ حل XStain BeadChip 1 | الكواشف | الموردة Illumina | 11208288 |
| XC2 و [مدش] ؛ حل XStain BeadChip 2 | كواشف Illumina | 11208296 | الموردة |
| XC3 و [مدش] ؛ حل XStain BeadChip 3 | كواشف Illumina | 11208392 | الموردة |
| XC4 و [مدش] ؛ حل XStain BeadChip 4 | كواشف | Illumina | 11208430 |
طلب إذن لإعادة استخدام النص أو الأشكال في مقالة JoVE هذه
طلب إذن