เผยแพร่เมื่อ | อัปเดตล่าสุดเมื่อ

คำตอบที่รวดเร็ว: การสุ่มแบบเมนเดล (Mendelian randomization หรือ MR) ใช้ตัวแปรทางพันธุกรรมที่ถ่ายทอดมาเป็นตัวแปรเครื่องมือเพื่อทดสอบว่าการได้รับปัจจัยเสี่ยงเป็นสาเหตุของผลลัพธ์หรือไม่ ซึ่งสามารถลดอคติจากตัวแปรแทรกซ้อนและสาเหตุย้อนกลับได้1,2 ผลลัพธ์จะน่าเชื่อถือได้ก็ต่อเมื่อเป็นไปตามข้อสมมติฐานสามประการ ได้แก่ ความเกี่ยวข้อง ความเป็นอิสระ และข้อจำกัดในการยกเว้น3 ในตัวอย่างที่เผยแพร่ภายใต้ลิขสิทธิ์ CC BY (ตัวแปรทางพันธุกรรม 185 ตัว, ไขมันในเลือด และโรคหลอดเลือดหัวใจ) การประมาณค่าถ่วงน้ำหนักแบบผกผันของความแปรปรวนทำให้ระดับคอเลสเตอรอล HDL ดูเหมือนมีผลป้องกันโรค (OR 0.78, 95% CI 0.68 ถึง 0.89) แต่ค่ามัธยฐานถ่วงน้ำหนัก (OR 0.93, 95% CI 0.81 ถึง 1.07) และ MR-Egger (OR 0.99, 95% CI 0.79 ถึง 1.24) กลับไม่เป็นเช่นนั้น4 เราคำนวณอัตราส่วนความเสี่ยงและช่วงความเชื่อมั่นใหม่ใน R จากค่าประมาณที่เผยแพร่แล้ว เราไม่ได้วิเคราะห์ข้อมูลระดับตัวแปรใหม่ ควรใช้ตัวประมาณค่าหลายตัว ไม่ใช่แค่ตัวเดียว

การสุ่มแบบเมนเดลคืออะไร?

การสุ่มแบบเมนเดล (Mendelian randomization) เป็นวิธีการใช้ความแปรผันทางพันธุกรรมที่ถ่ายทอดมาเพื่อถามว่าปัจจัยเสี่ยงที่สามารถปรับเปลี่ยนได้ เช่น คอเลสเตอรอล LDL ก่อให้เกิดผลลัพธ์ เช่น โรคหลอดเลือดหัวใจหรือไม่ วิธีนี้สร้างขึ้นจากกฎการถ่ายทอดทางพันธุกรรมของเมนเดลและการประมาณค่าตัวแปรเครื่องมือ ซึ่งทำให้สามารถประมาณผลกระทบเชิงสาเหตุได้แม้ว่าจะไม่เคยมีการวัดตัวแปรแทรกซ้อนบางอย่างก็ตาม5ข้อโต้แย้งดั้งเดิมคือ การจัดเรียงยีนแบบสุ่มจากพ่อแม่สู่ลูกในระหว่างการสร้างเซลล์สืบพันธุ์และการปฏิสนธิ หมายความว่าโพลีมอร์ฟิซึมที่เลียนแบบความเชื่อมโยงทางชีวภาพกับปัจจัยเสี่ยงนั้น โดยทั่วไปจะไม่ไวต่อสาเหตุย้อนกลับหรือตัวแปรแทรกซ้อนที่บิดเบือนการศึกษาเชิงสังเกตแบบดั้งเดิม1

กล่าวโดยง่ายคือ หากความแปรผันทางพันธุกรรมที่ส่งผลต่อระดับสารก่อภูมิแพ้ เช่น คอเลสเตอรอล LDL เกิดขึ้นโดยบังเอิญ การเปรียบเทียบผู้ที่มีและไม่มีความแปรผันทางพันธุกรรมนั้น ก็เหมือนกับการเปรียบเทียบสองกลุ่มที่ได้จากการโยนเหรียญ โดยมีเงื่อนไขว่าความแปรผันทางพันธุกรรมนั้นส่งผลต่อผลลัพธ์ผ่านทางสารก่อภูมิแพ้นั้นเท่านั้น

การวิเคราะห์ทุกครั้งสร้างขึ้นจากองค์ประกอบสามส่วนที่เชื่อมโยงกัน

  • เครื่องดนตรี การกลายพันธุ์ทางพันธุกรรมอย่างน้อยหนึ่งอย่างที่เกี่ยวข้องกับการสัมผัสสารดังกล่าว
  • การได้รับรังสี ปัจจัยเสี่ยงที่สามารถปรับเปลี่ยนได้ซึ่งอยู่ระหว่างการศึกษา เช่น ระดับคอเลสเตอรอลชนิด LDL
  • ผลลัพธ์. โรคหรือลักษณะที่สนใจ เช่น โรคหลอดเลือดหัวใจตีบ

วิธีการนี้มีข้อจำกัดที่การทบทวนเบื้องต้นได้ระบุไว้แล้ว ได้แก่ การรบกวนจากโพลีมอร์ฟิซึมที่อยู่ในภาวะสมดุลการเชื่อมโยงกับตัวแปร โพลีมอร์ฟิซึมที่มีผลต่อฟีโนไทป์หลายอย่าง การขาดโพลีมอร์ฟิซึมที่เหมาะสมสำหรับการสัมผัสหลายอย่าง และการปรับสมดุล ซึ่งเป็นการลดทอนผลกระทบทางพันธุกรรมในระหว่างการพัฒนา1

ข้อสมมติฐานหลักสามประการที่อยู่เบื้องหลังการศึกษาแบบสุ่มทางพันธุกรรมทุกครั้ง

ข้อ 5 ของ STROBE-MR ขอให้ผู้เขียนระบุสมมติฐานตัวแปรเครื่องมือหลักสามประการอย่างชัดเจนสำหรับการวิเคราะห์หลัก พร้อมทั้งสมมติฐานของการวิเคราะห์เพิ่มเติมหรือการวิเคราะห์ความไวใดๆ ข้อ 7 ขอให้พวกเขาอธิบายวิธีการหรือความรู้ก่อนหน้าที่ใช้ในการประเมินสมมติฐานเหล่านั้น3

เลื่อนไปด้านข้างเพื่อดูทุกคอลัมน์

ข้อสมมติพูดให้เข้าใจง่ายๆ ก็คือวิธีการตรวจสอบ
ติดต่อโฆษณาตัวแปรดังกล่าวมีความเกี่ยวข้องกับการสัมผัสสารนั้นรายงานค่าความแปรปรวนที่อธิบายได้ (R²) และโดยเฉพาะอย่างยิ่งในกรณีตัวอย่างเดียว ให้รายงานค่าสถิติ F เพื่อเป็นตัววัดความแข็งแกร่งของเครื่องมือ6
ความเป็นอิสระตัวแปรดังกล่าวไม่เกี่ยวข้องกับผลลัพธ์ผ่านกลไกที่ทำให้เกิดความสับสนความถูกต้องขึ้นอยู่กับความเป็นอิสระจากปัจจัยรบกวนในระดับประชากร การวิเคราะห์ภายในครอบครัวเป็นอีกทางเลือกหนึ่งหากข้อมูลและกำลังการวิเคราะห์เอื้ออำนวย6
การจำกัดการยกเว้นตัวแปรดังกล่าวไม่มีผลต่อผลลัพธ์ ยกเว้นอาจผ่านทางการสัมผัสการตรวจสอบภาวะพลีโอโทรปี: วิธี MR-Egger, ค่ามัธยฐานถ่วงน้ำหนัก, สถิติความแปรปรวน และวิธีการหาค่าผิดปกติ ซึ่งจะกล่าวถึงในหัวข้อถัดไป

ข้อสมมติฐานเหล่านี้เหมือนกับที่เบอร์เจสและเพื่อนร่วมงานระบุไว้ในแนวทางปฏิบัติของพวกเขา คือ ตัวแปรนั้นเกี่ยวข้องกับการสัมผัส ไม่เกี่ยวข้องกับผลลัพธ์เนื่องจากกลไกการรบกวน และไม่มีผลต่อผลลัพธ์ ยกเว้นอาจจะผ่านทางการสัมผัส6

เครื่องมือที่อ่อนแอ

เครื่องมือที่อ่อนแอจะทำให้การประมาณค่าเอนเอียงไปในทิศทางของความสัมพันธ์เชิงสังเกตที่สับสน และขนาดของความเอนเอียงขึ้นอยู่กับค่าสถิติ F เบอร์เจสและทอมป์สันยังแสดงให้เห็นว่ากฎทั่วไปที่อ้างถึงกันว่าค่า F มากกว่า 10 จะหลีกเลี่ยงความเอนเอียงได้นั้นเป็นสิ่งที่ทำให้เข้าใจผิด7ในการจำลองด้วยตัวแปรที่ไม่สัมพันธ์กันและข้อมูลสรุป เบอร์เจสและเพื่อนร่วมงานพบว่าความเอนเอียงของเครื่องมือที่อ่อนแอมีขนาดเล็กเมื่อแต่ละตัวแปรมีค่า P น้อยกว่า 0.00001 สำหรับความสัมพันธ์กับปัจจัยเสี่ยง8บทความที่เกี่ยวข้องของเราเกี่ยวกับแนวทางตัวแปรเครื่องมือครอบคลุมตรรกะทั่วไปนอกเหนือจากพันธุศาสตร์

การออกแบบการทดลองแบบตัวอย่างเดียว สองตัวอย่าง และหลายตัวแปร รวมถึงตัวประมาณค่า

ในการวิเคราะห์ตัวอย่างเดียว ข้อมูลเกี่ยวกับปัจจัยเสี่ยงและผลลัพธ์มาจากแหล่งข้อมูลเดียว ในการวิเคราะห์สองตัวอย่าง ข้อมูลมาจากชุดข้อมูลที่ไม่ทับซ้อนกัน9ข้อมูลสรุปสาธารณะทำให้การวิเคราะห์สองตัวอย่างมีความโปร่งใส เนื่องจากบุคคลที่สามที่มีข้อมูลเดียวกันสามารถทำซ้ำได้ ในขณะที่ข้อมูลระดับบุคคลในการตั้งค่าตัวอย่างเดียวจำเป็นสำหรับการตรวจสอบผลกระทบที่ไม่เป็นเชิงเส้น6อคติจากเครื่องมือที่อ่อนแอแตกต่างกันเช่นกัน: ในการตั้งค่าตัวอย่างเดียว อคติจะมุ่งไปสู่ความสัมพันธ์ที่สับสน และในการตั้งค่าสองตัวอย่างที่ไม่ทับซ้อนกัน อคติจะมุ่งไปสู่ค่าว่าง6

โดยปกติแล้วจะมีการเรียกใช้โปรแกรมประมาณการหลายโปรแกรมพร้อมกัน

  • การถ่วงน้ำหนักแบบผกผันความแปรปรวน (IVW) เป็นวิธีที่ใช้กันทั่วไปมากที่สุด โดยมีหลายรูปแบบ เป็นวิธีที่มีประสิทธิภาพมากที่สุด แต่จะมีความคลาดเคลื่อนหากค่าเฉลี่ยของผลกระทบแบบหลายยีนแตกต่างจากศูนย์6 สามารถหาได้จากข้อมูลสรุปโดยใช้ค่าเฉลี่ยถ่วงน้ำหนักผกผันของความแปรปรวนของความสัมพันธ์เฉพาะตัวแปร8
  • การถดถอย MR-Egger วิธีนี้ใช้ค่าความชันเป็นค่าประมาณเชิงสาเหตุ และค่าจุดตัดแกนเป็นค่าเฉลี่ยของผลกระทบแบบหลายยีน วิธีนี้อนุญาตให้ทุกตัวแปรมีผลกระทบแบบหลายยีนได้ แต่ต้องอาศัยข้อสมมติฐานของ InSIDE และมักจะมีประสิทธิภาพน้อยกว่า10,6
  • ค่ามัธยฐานถ่วงน้ำหนัก มีความสอดคล้องแม้ว่าข้อมูลมากถึง 50% จะมาจากเครื่องมือที่ไม่ถูกต้องก็ตาม4
  • MR หลายตัวแปร วิธีนี้ใช้ตัวแปรที่เกี่ยวข้องกับการสัมผัสหลายอย่างที่สัมพันธ์กันเพื่อประเมินผลกระทบของการสัมผัสแต่ละครั้งต่อผลลัพธ์หนึ่งอย่าง และสามารถใช้เพื่อประเมินการเป็นตัวกลางภายในกรอบงาน MR ได้11

สำหรับการจำลองที่นำตัวประมาณค่าเหล่านี้ไปใช้กับภาวะพลีโอโทรปีประเภทต่างๆ โปรดดูบทความของ Editverse เกี่ยวกับสิ่งที่ตัวประมาณค่าทั้งสี่ตัวนี้ช่วยปกป้องคุณจากอะไรบ้าง

กรณีศึกษาจริง: คอเลสเตอรอล HDL และ LDL กับโรคหลอดเลือดหัวใจ

Bowden และเพื่อนร่วมงานได้ตีพิมพ์บทความแบบเปิดเผย (CC BY 4.0) ที่ทดสอบค่ามัธยฐานถ่วงน้ำหนักเทียบกับ IVW และ MR-Egger บนข้อมูลสรุปสำหรับตัวแปรทางพันธุกรรม 185 ตัวจาก Do และเพื่อนร่วมงาน โดยใช้ความสัมพันธ์ของไขมันจาก Global Lipids Genetics Consortium และความสัมพันธ์ของโรคหลอดเลือดหัวใจจาก CARDIoGRAM consortium 4ค่าประมาณคืออัตราส่วนลอการิทึมของอัตราต่อรองต่อไขมันที่สูงขึ้น 1 ส่วนเบี่ยงเบนมาตรฐาน เราใช้ค่าประมาณและข้อผิดพลาดมาตรฐานจากตารางที่ 4 ของพวกเขา (คอลัมน์สำหรับตัวแปรทั้งหมดที่เกี่ยวข้องกับการสัมผัสที่ระดับนัยสำคัญทั่วทั้งจีโนม) และคำนวณอัตราส่วนอัตราต่อรอง ช่วงความเชื่อมั่น 95% และค่า p ใหม่ใน R โดยใช้การประมาณค่าปกติ4

แผนภาพฟอเรสต์แสดงอัตราส่วนความเสี่ยงของโรคหลอดเลือดหัวใจต่อระดับไขมันที่สูงขึ้น 1 SD โดยใช้ตัวประมาณค่าแบบสุ่มทางพันธุกรรม (Mendelian randomization estimator) คอเลสเตอรอล LDL: IVW 1.62, ค่ามัธยฐานถ่วงน้ำหนัก 1.58, MR-Egger 1.85 คอเลสเตอรอล HDL: IVW 0.78, ค่ามัธยฐานถ่วงน้ำหนัก 0.93, MR-Egger 0.99 โดยค่าจุดตัดของ MR-Egger บ่งชี้ถึงภาวะพหุผล (pleiotropy)
อัตราส่วนความเสี่ยงของโรคหลอดเลือดหัวใจต่อระดับไขมันในเลือดที่สูงขึ้น 1 SD แหล่งที่มา: Bowden et al. 2016 (CC BY 4.0), ตารางที่ 4; อัตราส่วนความเสี่ยงและช่วงความเชื่อมั่นคำนวณใน R โดย Editverse

เลื่อนไปด้านข้างเพื่อดูทุกคอลัมน์

ประมาณการคอเลสเตอรอล LDL, OR (95% CI)คอเลสเตอรอล HDL, OR (95% CI)
ถ่วงน้ำหนักด้วยค่าผกผันของความแปรปรวน1.62 (1.44 ถึง 1.82)0.78 (0.68 ถึง 0.89)
ค่ามัธยฐานถ่วงน้ำหนัก1.58 (1.39 ถึง 1.80)0.93 (0.81 ถึง 1.07)
ความชัน MR-Egger1.85 (1.51 ถึง 2.27)0.99 (0.79 ถึง 1.24)
ค่าจุดตัด MR-Egger (SE) ที่คำนวณได้ p-0.009 (0.005), p = 0.072-0.014 (0.005), p = 0.005

ค่าคอเลสเตอรอล LDL ให้คำตอบที่สอดคล้องกันในตัวประมาณค่าทั้งสามตัว โดยไม่มีสัญญาณที่ชัดเจนของภาวะพลีโอโทรปีแบบมีทิศทาง (ค่าจุดตัด MR-Egger p = 0.072) สำหรับค่าคอเลสเตอรอล HDL ช่วง IVW ไม่รวม 1 ในขณะที่ช่วงค่ามัธยฐานถ่วงน้ำหนักและช่วง MR-Egger รวม 1 และค่าจุดตัด MR-Egger แตกต่างจากศูนย์ (p = 0.005) ผู้เขียนสรุปว่าผลกระทบเชิงสาเหตุของค่า IVW ของคอเลสเตอรอล HDL นั้นน่าสงสัย4เมื่อพวกเขาจำกัดการวิเคราะห์เฉพาะตัวแปรที่มีความสัมพันธ์ที่แข็งแกร่งที่สุดกับคอเลสเตอรอล HDL รูปแบบยังคงอยู่: IVW OR 0.87 (0.77 ถึง 0.99), ค่ามัธยฐานถ่วงน้ำหนัก OR 0.94 (0.82 ถึง 1.06), ค่าความชัน MR-Egger OR 1.10 (0.89 ถึง 1.35) (คำนวณโดย Editverse) 4

กรณีนี้แสดงให้เห็นว่า:รายงาน IVW เพียงอย่างเดียวจะแสดงผลว่าคอเลสเตอรอล HDL มีผลดีต่อสุขภาพ และตัวประมาณค่าเพิ่มเติมได้เปลี่ยนค่าที่ได้ แต่กรณีนี้ไม่ได้แสดงให้เห็นว่า:คอเลสเตอรอล HDL ไม่มีผลใดๆ เนื่องจากช่วงค่ามัธยฐานถ่วงน้ำหนักและช่วง MR-Egger นั้นกว้าง หรือว่าตัวประมาณค่าถูกต้องเสมอไป เราได้วาดกราฟใหม่จากตัวเลขที่ตีพิมพ์แล้ว และไม่ได้ทำการคำนวณข้อมูล 185 ตัวแปรซ้ำ เนื่องจากเราไม่ได้รับข้อมูลเหล่านั้น

พลีโอโทรปี ภัยคุกคามสำคัญต่อความถูกต้อง

ภาวะพลีโอโทรปีในแนวนอน ซึ่งตัวแปรหนึ่งส่งผลต่อผลลัพธ์โดยผ่านเส้นทางที่ข้ามการสัมผัส ทำให้ข้อจำกัดการยกเว้นถูกละเมิด Verbanck และคณะได้พัฒนา MR-PRESSO เพื่อระบุค่าผิดปกติของพลีโอโทรปีในแนวนอนใน MR ระดับสรุปแบบหลายเครื่องมือ และนำไปใช้กับลักษณะและโรคที่ซับซ้อน พวกเขาพบว่าสามารถตรวจพบพลีโอโทรปีในแนวนอนได้ในความสัมพันธ์เชิงสาเหตุที่มีนัยสำคัญมากกว่า 48% โดยมีการบิดเบือนของการประมาณค่าเชิงสาเหตุโดยเฉลี่ยตั้งแต่ -131% ถึง 201% และความสัมพันธ์เชิงสาเหตุที่เป็นผลบวกเท็จในความสัมพันธ์มากถึง 10% 12ตารางที่ 1 ของแนวทาง Burgess ระบุว่า MR-PRESSO มีอัตราผลบวกเท็จสูงมากเมื่อเครื่องมือหลายตัวไม่ถูกต้อง ดังนั้นจึงเป็นการตรวจสอบเพียงหนึ่งในหลายๆ วิธี6

การวิเคราะห์ที่น่าเชื่อถือจะต้องตรวจสอบหลายประเด็นร่วมกัน

  • การสกัดกั้น MR-Egger และความชันเป็นค่าประมาณความไว10
  • ค่ามัธยฐานถ่วงน้ำหนัก และวิธีการประมาณค่าที่มีประสิทธิภาพอื่นๆ4
  • ความแตกต่างหลากหลาย แนวทางดังกล่าวแนะนำให้รายงานค่า Cochran's Q หรือค่าสถิติ I² สำหรับการวิเคราะห์แบบพหุพันธุกรรม6
  • การกรองแบบ Leave-one-out และ Steiger STROBE-MR แสดงรายการกราฟแบบตัดออกหนึ่งรายการ (leave-one-out plots) ไว้ในกลุ่มกราฟเพิ่มเติม3 และการกรองแบบ Steiger จะกำจัดตัวแปรที่มีความสัมพันธ์กับผลลัพธ์มากกว่าความสัมพันธ์กับปัจจัยเสี่ยง6
  • ชุดตัวเลือกแบบเสรีนิยมและแบบอนุรักษ์นิยม โดยทั่วไปแล้ว ปรากฏการณ์พลีโอโทรปีมีแนวโน้มที่จะทำให้การประมาณค่าเบี่ยงเบนไปจากสมมติฐานว่าง ดังนั้นผลลัพธ์ที่เป็นศูนย์ในการวิเคราะห์แบบเสรีนิยมจึงเป็นหลักฐานที่น่าเชื่อถือมากกว่าว่าสมมติฐานว่างนั้นเป็นจริง6

การทับซ้อนของตัวอย่างและอคติจากเครื่องมือที่ไม่แข็งแรง

กลุ่มวิจัยขนาดใหญ่หลายกลุ่มมีผู้เข้าร่วมวิจัยร่วมกัน ดังนั้นชุดข้อมูลการสัมผัสและผลลัพธ์ในการวิเคราะห์แบบสองกลุ่มตัวอย่างจึงอาจทับซ้อนกันได้ เบอร์เจสและคณะได้จำลองสถานการณ์นี้สำหรับผลลัพธ์แบบต่อเนื่องและผลกระทบเชิงสาเหตุที่เป็นศูนย์: การวิเคราะห์แบบกลุ่มตัวอย่างเดียวที่มีอคติสัมพัทธ์ 10% (ซึ่งสอดคล้องกับพารามิเตอร์ F เท่ากับ 10) ลดลงเหลือ 5% เมื่อมีการทับซ้อน 50% และ 3% เมื่อมีการทับซ้อน 30% เนื่องจากอคติเป็นฟังก์ชันเชิงเส้นของสัดส่วนการทับซ้อน9ข้อ 10d ของ STROBE-MR ขอให้ผู้เขียนงานวิจัยแบบสองกลุ่มตัวอย่างรายงานจำนวนบุคคลที่ทับซ้อนกันระหว่างการศึกษาการสัมผัสและผลลัพธ์3

แผนภูมิแท่งแสดงความคลาดเคลื่อนสัมพัทธ์ของผลกระทบเชิงสาเหตุที่เป็นศูนย์ที่ระดับความแข็งแกร่งของเครื่องมือ F เท่ากับ 10: 3% เมื่อผู้เข้าร่วมซ้ำกัน 30%, 5% เมื่อซ้ำกัน 50% และ 10% เมื่อซ้ำกัน 100% (ตัวอย่างเดียว) จากการศึกษาจำลอง
ความเอนเอียงสัมพัทธ์ที่ F = 10 ในการจำลอง แหล่งที่มา: Burgess et al. 2016 (CC BY 4.0); ตรวจสอบความสัมพันธ์เชิงเส้นใน R โดย Editverse นี่เป็นผลลัพธ์จากการจำลอง ไม่ใช่ชุดข้อมูลจริง

ตัวอย่างที่ตีพิมพ์เพิ่มเติม

คอเลสเตอรอล LDL และโรคหลอดเลือดหัวใจ

เฟเรนซ์และคณะได้รวมการศึกษาแบบสุ่มทางพันธุกรรมของโพลีมอร์ฟิซึม 9 ตำแหน่งใน 6 ยีน โดยไม่พบหลักฐานของความแตกต่างกัน ในข้อมูลที่ไม่ซ้ำซ้อนจากผู้เข้าร่วม 312,321 คน การจัดสรรแบบสุ่มตามธรรมชาติให้มีระดับคอเลสเตอรอล LDL ต่ำในระยะยาวมีความสัมพันธ์กับการลดความเสี่ยงของโรคหลอดเลือดหัวใจลง 54.5% (ช่วงความเชื่อมั่น 95% 48.8% ถึง 59.5%) ต่อคอเลสเตอรอล LDL ที่ลดลง 1 มิลลิโมล/ลิตร ซึ่งประมาณ 3 เท่าของการลดลงที่พบระหว่างการรักษาด้วยยาสแตตินที่เริ่มในภายหลัง13เบอร์เจสและคณะใช้เพียง 5 ตัวแปรและข้อมูลสรุป ประมาณการว่าการลดคอเลสเตอรอล LDL ลง 30% จะช่วยลดความเสี่ยงของโรคหลอดเลือดหัวใจลง 67% (ช่วงความเชื่อมั่น 95% 54% ถึง 76%) 8

การสุ่มแบบเมนเดลของเป้าหมายยา

Gagnon และเพื่อนร่วมงานได้ทำการวิเคราะห์จีโนไทป์และลำดับอาร์เอ็นเอของตัวอย่างตับมนุษย์ 504 ตัวอย่าง และระบุตัวแปร (rs472495) ที่อธิบายการแสดงออกของPCSK9 ในตับได้ 5.6% เพื่อเลียนแบบการรบกวนอาร์เอ็นเอ และใช้PCSK9 R46L เพื่อจำลองการยับยั้งโดยใช้แอนติบอดี สำหรับการลดลงหนึ่งค่าเบี่ยงเบนมาตรฐานของอะโพลิโปโปรตีน B อัตราส่วนความเสี่ยงของโรคหลอดเลือดหัวใจจะอยู่ที่ 0.40 (95% CI 0.31 ถึง 0.51) และ 0.48 (0.43 ถึง 0.55) ตามลำดับ14

วิตามิน D

ฟางและคณะได้รวบรวมงานวิจัยแบบสุ่มทางพันธุกรรม (Mendelian randomization) จำนวน 133 ฉบับ ครอบคลุมผลลัพธ์ 275 รายการ จนถึงเดือนกันยายน 2022 พบว่าระดับ 25-ไฮดรอกซีวิตามินดีที่สูงในเชิงพันธุกรรมมีความสัมพันธ์กับความเสี่ยงที่ลดลงของโรคปลอกประสาทเสื่อมแข็ง (multiple sclerosis) และผลลัพธ์อื่นๆ รวมถึงความเสี่ยงที่สูงขึ้นของบางโรค เช่น โรคนิ่วในไต ความสัมพันธ์กับโรคปลอกประสาทเสื่อมแข็งมีความสอดคล้องกันในชุดเครื่องมือต่างๆ ในขณะที่ความสัมพันธ์อื่นๆ ส่วนใหญ่ปรากฏเฉพาะกับตัวแปรที่ไม่จำกัดเฉพาะวิถีทางของวิตามินดี และคุณภาพของการศึกษามีความหลากหลาย15จ้าวและคณะใช้ผู้เข้าร่วมจาก UK Biobank จำนวน 332,984 คน โดย 23,089 คนมีโรคภูมิต้านตนเองอย่างน้อยหนึ่งโรค และรายงานอัตราส่วนความเสี่ยง (odds ratio) ที่ 0.95 (95% CI 0.91 ถึง 0.99) ต่อ 25-ไฮดรอกซีโคลแคลซิเฟอรอลที่สูงขึ้น 10 ng/ml สำหรับกลุ่มที่มีการอักเสบจากภูมิคุ้มกันตนเอง แต่ไม่ใช่กลุ่มที่มีภูมิคุ้มกันตนเอง16

แอลกอฮอล์และโรคหลอดเลือดในประเทศจีน

มิลล์วูดและเพื่อนร่วมงานได้ทำการวิเคราะห์จีโนไทป์ของผู้เข้าร่วม 161,498 คนจากธนาคารชีวภาพจีนคาดูรี สำหรับ ยีน ALDH2 -rs671 และADH1B -rs1229984 การบริโภคที่รายงานด้วยตนเองมีความสัมพันธ์แบบรูปตัวยูกับโรคหลอดเลือดสมองและกล้ามเนื้อหัวใจตายในผู้ชาย แต่การบริโภคที่ทำนายจากจีโนไทป์ไม่มีความสัมพันธ์ดังกล่าว ความเสี่ยงสัมพัทธ์ต่อการบริโภค 280 กรัมต่อสัปดาห์อยู่ที่ 1.58 (1.36 ถึง 1.84) สำหรับภาวะเลือดออกในสมอง 1.27 (1.13 ถึง 1.43) สำหรับโรคหลอดเลือดสมองตีบ และ 0.96 (0.78 ถึง 1.18) สำหรับกล้ามเนื้อหัวใจตาย ผู้เขียนสรุปว่าผลที่เห็นได้ชัดว่าเป็นการป้องกันโรคหลอดเลือดสมองจากการดื่มแอลกอฮอล์ในปริมาณปานกลางนั้นส่วนใหญ่ไม่ใช่สาเหตุโดยตรง17

ควรใช้การสุ่มแบบเมนเดลเมื่อใด?

  • การทดลองแบบสุ่มไม่สามารถทำได้จริง วิธีการนี้มีประโยชน์มากที่สุดเมื่อไม่สามารถทำการทดลองแบบสุ่มได้ และการศึกษาเชิงสังเกตมีอคติเนื่องจากตัวแปรแทรกซ้อนหรือสาเหตุย้อนกลับ2
  • คุณจำเป็นต้องตรวจสอบเป้าหมายของยา มีการนำ MR มาใช้ในการทำนายประสิทธิภาพและความปลอดภัยของยาที่มุ่งเป้าไปที่ปัจจัยเสี่ยงของโรคหัวใจและหลอดเลือด รวมถึงการสำรวจการนำยาที่มีอยู่แล้วมาใช้ใหม่18
  • มีข้อมูลสรุปที่มีขนาดเหมาะสมอยู่แล้ว ข้อมูลสรุป GWAS ที่เผยแพร่สู่สาธารณะช่วยให้สามารถใช้การออกแบบสองกลุ่มตัวอย่างได้โดยไม่ต้องระบุจีโนไทป์ของแต่ละบุคคล6

วิธีการนี้ไม่เหมาะสมเมื่อไม่สามารถหาเหตุผลสนับสนุนตัวแปรทางพันธุกรรมใดๆ มาใช้เป็นเครื่องมือในการประเมินความเสี่ยงได้ เนื่องจากสมมติฐานเรื่องความเกี่ยวข้องจะล้มเหลว หรือเมื่อข้อมูลที่มีอยู่มีน้อยเกินไป เพราะการศึกษาแบบ MR มักมีกำลังการทดสอบต่ำ6หากคุณกำลังเลือกระหว่างการออกแบบทางพันธุกรรมกับการวิเคราะห์แบบสำรวจหรือทะเบียนข้อมูล ให้ตรวจสอบว่าควรใช้ฐานข้อมูลสุขภาพใดในการศึกษาของคุณ

วิธีการดำเนินการศึกษาแบบสุ่มทางพันธุกรรม (Mendelian randomization)

  1. 1ระบุสมมติฐานเชิงสาเหตุ ข้อ 2 และ 3 ของแบบสอบถาม STROBE-MR ถามถึงเหตุผลว่าทำไม MR จึงเป็นวิธีการที่มีประโยชน์ และวัตถุประสงค์เฉพาะ รวมถึงสมมติฐานเชิงสาเหตุที่กำหนดไว้ล่วงหน้า3
  2. 2เลือกแบบและข้อมูล ตัดสินใจว่าจะใช้ข้อมูลตัวอย่างเดียวหรือสองตัวอย่าง และข้อมูลระดับบุคคลหรือข้อมูลสรุป6
  3. 3เลือกตัวเลือกต่างๆ แนวทางทั่วไปคือการรวมตัวแปรที่เกี่ยวข้องกับการสัมผัสที่ระดับนัยสำคัญทั่วทั้งจีโนม โดยทั่วไปค่า P ต่ำกว่า 5 10 ×-8การเลือกใช้ชุดข้อมูลเดียวกันกับการวิเคราะห์อาจก่อให้เกิดความเสี่ยงต่อปรากฏการณ์ "คำสาปของผู้ชนะ" และทำให้ความลำเอียงจากเครื่องมือที่อ่อนแอแย่ลงไปอีก6
  4. 4ปรับสมดุลอัลลีล จัดเรียงอัลลีลที่เกี่ยวข้องกับการสัมผัสและผลลัพธ์ สำหรับตัวแปรแบบพาลินโดรม (A/T หรือ C/G) ที่มีความถี่ของอัลลีลใกล้เคียง 50% นักวิเคราะห์อาจตัดตัวแปรนั้นทิ้งหากไม่สามารถตรวจสอบทิศทางได้6
  5. 5ดำเนินการวิเคราะห์ข้อมูลเบื้องต้นและการวิเคราะห์ความไว ใช้ IVW เป็นตัวประมาณค่าหลักสำหรับตัวแปรหลายตัว จากนั้นจึงใช้ตัวประมาณค่าที่แข็งแกร่ง สถิติความแปรปรวน และแผนภาพ6
  6. 6รายงานถึง STROBE-MR. รายการตรวจสอบประกอบด้วยรายการหลัก 20 รายการ และรายการย่อย 30 รายการ รวมถึงเวอร์ชันซอฟต์แวร์และแพ็กเกจ (รายการ 9a) สถานะการลงทะเบียนล่วงหน้า (รายการ 9b) และการแบ่งปันข้อมูลและรหัส (รายการ 19)19,3

คำแถลง STROBE-MR ขยายความตรรกะทั่วไปของแนวทางการรายงาน STROBE

แหล่งข้อมูลและซอฟต์แวร์สำหรับการวิเคราะห์ MR

คำอธิบายด้านล่างนี้คัดมาจากหน้าเว็บของแต่ละแหล่งข้อมูล โดยอ่านเมื่อวันที่ 4 พฤศจิกายน 2026

เลื่อนไปด้านข้างเพื่อดูทุกคอลัมน์

ทรัพยากรมันคืออะไรการใช้งานทั่วไป
OpenGWAS (MRC IEU)20ฐานข้อมูลที่มีชุดข้อมูลสรุป GWAS มากกว่า 50,000 ชุด สามารถสืบค้นหรือดาวน์โหลดได้สถิติสรุปการได้รับและผลลัพธ์สำหรับการออกแบบการทดลองแบบสองกลุ่มตัวอย่าง
ฟินน์เจน21โครงการที่รวบรวมและวิเคราะห์ข้อมูลทางพันธุกรรมและข้อมูลสุขภาพจากผู้บริจาคไบโอแบงก์ชาวฟินแลนด์จำนวน 500,000 คนความสัมพันธ์ทางพันธุกรรมระหว่างผลลัพธ์ของโรค
แคตตาล็อก NHGRI-EBI GWAS22แคตตาล็อกการศึกษาความสัมพันธ์ทั่วทั้งจีโนมมนุษย์ของ NHGRI-EBI ส่วนสถิติสรุปมีเอกสารและสามารถเข้าถึงสถิติสรุปฉบับเต็มสำหรับการศึกษาต่างๆ ได้ (หากมี)การค้นหาและตรวจสอบความสัมพันธ์ระหว่างลักษณะทางพันธุกรรมกับลักษณะเด่น
TwoSampleMR (แพ็กเกจ R)23แพ็กเกจ R สำหรับการวิเคราะห์ MR ด้วยข้อมูลสรุป GWAS ที่สามารถดึงข้อมูลจาก OpenGWAS ได้การประสานงาน, IVW, MR-Egger, ค่ามัธยฐานถ่วงน้ำหนัก
การสุ่มแบบเมนเดล (แพ็กเกจ R)24แพ็กเกจ CRAN ที่รวบรวมวิธีการ MR หลายวิธีสำหรับข้อมูลสรุปการนำตัวประมาณค่าหลักไปใช้ในรูปแบบอื่น
mrpresso (แพ็กเกจ R)6ระบุไว้ในแนวทางของเบอร์เจสว่าเป็นซอฟต์แวร์สำหรับ MR-PRESSOการตรวจจับค่าผิดปกติ โปรดดูข้อควรระวังด้านบน

ข้อผิดพลาดทั่วไปในการสุ่มแบบเมนเดล

  • ไม่ได้รายงานความแข็งแรงของเครื่องมือ รายงานค่า R² และค่าสถิติ F และโปรดจำไว้ว่ากฎทั่วไปที่ว่าค่า F มากกว่า 10 นั้นอาจทำให้เข้าใจผิดได้6,7
  • รายงานเฉพาะค่าประมาณ IVW เท่านั้น กรณีคอเลสเตอรอล HDL ข้างต้นแสดงให้เห็นว่าเครื่องมือประเมินเพียงตัวเดียวอาจทำให้เกิดความเข้าใจผิดได้4
  • ไม่สนใจการทับซ้อนของตัวอย่าง รายงานจำนวนบุคคลที่มีข้อมูลซ้ำซ้อน และคาดว่าอคติจะแตกต่างกันไปตามระดับความซ้ำซ้อน9,3
  • การผสมผสานประชากรโดยไม่ระมัดระวัง หากตัวอย่างทั้งสองมาจากกลุ่มชาติพันธุ์ที่แตกต่างกัน รูปแบบของภาวะไม่สมดุลของการเชื่อมโยงทางพันธุกรรมก็อาจแตกต่างกันไปด้วย ดังนั้นตัวแปรทางพันธุกรรมอาจมีความสัมพันธ์กับปัจจัยเสี่ยงในชุดข้อมูลผลลัพธ์น้อยลง6
  • การอ่านช่วงค่าว่างที่กว้างหมายถึงไม่มีผลกระทบใดๆ ค่าศูนย์ในกรณีของ HDL มาพร้อมกับช่วงตั้งแต่ประมาณ 0.8 ถึง 1.2 ซึ่งสอดคล้องกับผลกระทบเล็กน้อยในทั้งสองทิศทาง4
  • ใช้ภาษาเชิงสาเหตุอย่างไม่เคร่งครัดนัก ข้อที่ 16 ของ STROBE-MR ขอให้ผู้เขียนใช้ภาษาที่แสดงถึงความสัมพันธ์เชิงสาเหตุอย่างระมัดระวัง และชี้แจงให้ชัดเจนว่าการประมาณค่าอาจแสดงผลกระทบเชิงสาเหตุได้ภายใต้สมมติฐานบางประการเท่านั้น3

การสุ่มแบบเมนเดลเทียบกับวิธีการอนุมานเชิงสาเหตุอื่นๆ

ตารางนี้แสดงคุณลักษณะเด่นของแต่ละการออกแบบสำหรับปัญหาทั้งสองที่ MR มุ่งหมายจะแก้ไข

เลื่อนไปด้านข้างเพื่อดูทุกคอลัมน์

วิธีสับสนการก่อให้เกิดผลย้อนกลับ
สุ่มทดลองควบคุมการสุ่มจัดสรรการแทรกแซงตามหลักการแล้ว การแทรกแซงจะถูกกำหนดไว้ก่อนที่จะมีการวัดผลลัพธ์
การสุ่มแบบเมนเดเลียนจะลดลงหากตัวแปรนั้นเป็นอิสระจากปัจจัยรบกวน6จีโนไทป์ถูกกำหนดตั้งแต่การปฏิสนธิ ดังนั้นโดยทั่วไปจึงไม่เสี่ยงต่อการติดเชื้อ1
การจับคู่คะแนนความโน้มเอียงสมดุลจะวัดเฉพาะตัวแปรควบคุมเท่านั้นปัญหาดังกล่าวไม่ได้รับการแก้ไขโดยการจับคู่ตามตัวแปรพื้นฐาน
การถดถอยหลายตัวแปรปรับค่าเฉพาะตัวแปรที่วัดได้เท่านั้นไม่ได้พิจารณาโดยการปรับตัวแปรควบคุม

เนื่องจากการออกแบบแต่ละแบบมีแหล่งที่มาของอคติที่แตกต่างกัน การใช้สามเหลี่ยมในการตรวจสอบจะช่วยเสริมความมั่นใจเมื่อวิธีการเหล่านั้นสอดคล้องกัน และ MR สามารถใช้เป็นส่วนหนึ่งของการตรวจสอบได้25 , 5สำหรับข้อมูลพื้นฐานที่เกี่ยวข้อง โปรดดูคู่มือของเราเกี่ยวกับการจับคู่คะแนนความโน้มเอียงและการอนุมานเชิงสาเหตุในระบาดวิทยา

ต้องการดำเนินการวิเคราะห์หรือตรวจสอบผลหรือไม่? โปรดติดต่อ Editverse

ผู้เชี่ยวชาญเฉพาะด้านของ Editverse ทำงานโดยยึดหลักแนวทางของ COPE, ICMJE และ GPP 2022 พร้อมทั้งปฏิบัติตามหลักจริยธรรมเฉพาะสถาบันและเฉพาะประเทศด้วย

โลโก้ Editverse

แก้ไข

การศึกษาแบบสุ่มเมนเดล

ต้องการงานวิจัย MR แบบสองกลุ่มตัวอย่างฉบับสมบูรณ์ที่สร้างขึ้นตามมาตรฐานวารสารวิชาการหรือไม่?

คุณส่งข้อมูลต่อไปนี้มาให้เรา:ข้อมูลเกี่ยวกับปัจจัยเสี่ยง ผลลัพธ์ และคำถามเกี่ยวกับสาเหตุ สถิติสรุปเป็นข้อมูลสาธารณะ เราเป็นผู้จัดหาให้

คุณจะได้รับ:การศึกษา MR สองกลุ่มตัวอย่างพร้อมการวัดความแรงของเครื่องมือ การกรองแบบ Steiger ตัวประมาณค่าที่ทนทานต่อภาวะพลีโอโทรปีสี่ตัว การตรวจจับค่าผิดปกติ MR-PRESSO การวิเคราะห์ PheWAS ของเครื่องมือ และการทำซ้ำอย่างอิสระ รวมถึงต้นฉบับที่สมบูรณ์ แผนภาพกระจายและแผนภาพกรวย โค้ด TwoSampleMR ที่สามารถใช้งานได้ และรายการตรวจสอบ STROBE-MR

กำหนดเวลา: เมื่อคุณต้องการให้ดำเนินการศึกษาทั้งหมดและจัดทำรายงาน 15 วันทำการ

การศึกษาแบบสุ่มทางพันธุกรรม (Mendelian Randomization Study) · 2,900 ดอลลาร์สหรัฐ

วางแผนการศึกษา MR ของฉัน →

โลโก้ Editverse

แก้ไข

การตรวจสอบความเป็นไปได้

ไม่แน่ใจว่าการออกแบบ MR จะสามารถตอบคำถามของคุณได้หรือไม่?

คุณส่ง:คำถามวิจัยของคุณ

คุณจะได้รับ:คำแนะนำเป็นลายลักษณ์อักษรว่าคำถามนั้นได้รับการตอบแล้วหรือไม่ ข้อมูลสนับสนุนได้หรือไม่ และคุณมีกำลังการวิเคราะห์เพียงพอหรือไม่ รวมถึงความแข็งแกร่งของเครื่องมือสำหรับการออกแบบการสุ่มแบบเมนเดล (Mendelian randomization design) โดยจะได้รับเครดิตเต็มจำนวนสำหรับการศึกษาใดๆ ที่คุณว่าจ้างภายในหกเดือน

ช่วงเวลา:ก่อนที่คุณจะตัดสินใจเข้าร่วมการศึกษา 5 วันทำการ

การประเมินความเป็นไปได้และขั้นตอนการดำเนินงาน · 450 ดอลลาร์สหรัฐ

ตรวจสอบความเป็นไปได้ก่อน →

โลโก้ Editverse

แก้ไข

การวิเคราะห์ทางสถิติ

มีข้อมูลอยู่แล้วและต้องการให้ทำการวิเคราะห์ใน R ใช่หรือไม่?

คุณส่งข้อมูลของคุณ คำถามวิจัย และซอฟต์แวร์ที่คุณต้องการใช้มาได้เลย

คุณจะได้รับ:การวิเคราะห์ทางสถิติโดยใช้ซอฟต์แวร์ที่คุณต้องการ (SPSS, R, STATA หรือ SAS) พร้อมโค้ดที่มีคำอธิบายประกอบ ไฟล์ผลลัพธ์ที่สมบูรณ์ และคู่มือการตีความ

เมื่อไหร่:คุณต้องการผลลัพธ์ที่พร้อมสำหรับการตีพิมพ์ภายใน 4-6 วันทำการ

การวิเคราะห์ทางสถิติด้วย SPSS/R/STATA · 400 ดอลลาร์สหรัฐ

เริ่มการวิเคราะห์ของฉัน →

เครื่องมือวิจัยฟรีจาก Editverse

คำถามที่พบบ่อยเกี่ยวกับการสุ่มแบบเมนเดล

การสุ่มแบบเมนเดล (Mendelian randomization) เหมือนกับการทดลองแบบสุ่มที่มีกลุ่มควบคุม (randomized controlled trial) หรือไม่?

ไม่ การสุ่มแบบเมนเดล (Mendelian randomization) ใช้ได้กับข้อมูลจากการสังเกต และถือว่าตัวแปรทางพันธุกรรมที่ถ่ายทอดทางกรรมพันธุ์เป็นเครื่องมือ อัลลีลถูกส่งต่อจากพ่อแม่ไปยังลูกหลานโดยการแยกตัวแบบสุ่ม ซึ่งคล้ายกับการจัดสรรการรักษาแบบสุ่มในการทดลอง และควรลดปัจจัยรบกวนลง การประมาณค่านี้ยังขึ้นอยู่กับสมมติฐานที่การทดลองไม่จำเป็นต้องใช้ ดังนั้นจึงเป็นการเสริมหลักฐานจากการทดลอง ไม่ใช่การแทนที่หลักฐานจากการทดลอง

ข้อสมมติฐานหลักสามประการของการสุ่มแบบเมนเดลคืออะไร?

ความเกี่ยวข้อง ความเป็นอิสระ และข้อจำกัดในการยกเว้น กล่าวโดยง่ายคือ ตัวแปรนั้นเกี่ยวข้องกับการสัมผัส ไม่ได้เกี่ยวข้องกับผลลัพธ์ผ่านกลไกการรบกวน และไม่ส่งผลต่อผลลัพธ์ ยกเว้นอาจจะผ่านการสัมผัส ข้อ 5 ของ STROBE-MR ขอให้ผู้เขียนระบุสิ่งเหล่านี้อย่างชัดเจน

ฉันจำเป็นต้องมีข้อมูลจีโนไทป์ระดับบุคคลเพื่อใช้ในการวิเคราะห์แบบเมนเดลแรนติเซชันหรือไม่?

ไม่เหมาะสำหรับการวิเคราะห์แบบสองกลุ่มตัวอย่าง ซึ่งสามารถใช้ข้อมูลสรุปที่เผยแพร่แล้วเกี่ยวกับความสัมพันธ์ของแต่ละตัวแปรกับปัจจัยเสี่ยงและผลลัพธ์ได้ จำเป็นต้องใช้ข้อมูลระดับบุคคลในการวิเคราะห์แบบกลุ่มตัวอย่างเดียวเพื่อตรวจสอบผลกระทบที่ไม่เป็นเชิงเส้น ข้อดีในทางปฏิบัติของข้อมูลสรุปที่เผยแพร่ต่อสาธารณะคือบุคคลที่สามสามารถทำการวิเคราะห์ซ้ำได้

การสุ่มแบบเมนเดลโดยใช้กลุ่มตัวอย่างเดียวกับการสุ่มแบบเมนเดลโดยใช้กลุ่มตัวอย่างสองกลุ่มแตกต่างกันอย่างไร?

ในการวิเคราะห์แบบกลุ่มตัวอย่างเดียว ข้อมูลเกี่ยวกับปัจจัยเสี่ยงและผลลัพธ์มาจากแหล่งข้อมูลเดียว ในการวิเคราะห์แบบสองกลุ่มตัวอย่าง ข้อมูลมาจากชุดข้อมูลที่ไม่ทับซ้อนกัน ทิศทางของอคติจากเครื่องมือที่ไม่แข็งแรงจึงแตกต่างกัน โดยจะเอนเอียงไปทางความสัมพันธ์ที่ซับซ้อนในกรณีการวิเคราะห์แบบกลุ่มตัวอย่างเดียว และเอนเอียงไปทางค่าว่างในกรณีการวิเคราะห์แบบสองกลุ่มตัวอย่างโดยไม่มีการทับซ้อนกันของกลุ่มตัวอย่าง

เหตุใดจึงต้องรายงานค่าประมาณที่มากกว่าค่าประมาณถ่วงน้ำหนักตามส่วนกลับของความแปรปรวน?

เนื่องจากวิธีการถ่วงน้ำหนักแบบผกผันความแปรปรวนจะให้ค่าประมาณที่สอดคล้องกันก็ต่อเมื่อตัวแปรทั้งหมดถูกต้องหรือภาวะพลีโอโทรปีสมดุลกัน ในตัวอย่างไขมันที่ตีพิมพ์ข้างต้น อัตราส่วนความน่าจะเป็นแบบ IVW สำหรับคอเลสเตอรอล HDL คือ 0.78 ในขณะที่ค่ามัธยฐานแบบถ่วงน้ำหนักให้ค่า 0.93 และ MR-Egger ให้ค่า 0.99 โดยมีช่วงความเชื่อมั่นที่รวม 1 การรายงานค่าประมาณหลายค่าช่วยให้ผู้อ่านเห็นว่าผลลัพธ์ขึ้นอยู่กับวิธีการใดวิธีการหนึ่ง

ค่า MR-Egger intercept บอกอะไรฉันบ้าง?

ในการวิเคราะห์การถดถอยแบบ MR-Egger ค่าจุดตัดจะประมาณค่าเฉลี่ยของผลกระทบแบบหลายยีน (pleiotropic effect) ในทุกตัวแปร และค่าความชันจะให้ค่าประมาณเชิงสาเหตุ ค่าจุดตัดที่แตกต่างจากศูนย์บ่งชี้ถึงผลกระทบแบบหลายยีนที่มีทิศทาง วิธีนี้จำเป็นต้องใช้สมมติฐาน InSIDE ที่ว่าผลกระทบแบบหลายยีนเป็นอิสระจากความสัมพันธ์ระหว่างตัวแปรและการสัมผัส และมักจะมีประสิทธิภาพน้อยกว่าวิธีการประมาณค่าอื่นๆ

STROBE-MR ต้องการอะไรบ้าง?

STROBE-MR เป็นแนวทางการรายงานที่มีรายการหลัก 20 รายการและรายการย่อย 30 รายการ ตัวอย่างเช่น การระบุข้อสมมติฐานตัวแปรเครื่องมือสามประการ (รายการที่ 5) การอธิบายวิธีการประเมินข้อสมมติฐาน (รายการที่ 7) การรายงานการวิเคราะห์ความไว (รายการที่ 8 และ 13) และสำหรับแบบแผนสองกลุ่มตัวอย่าง จำนวนบุคคลที่ซ้ำกันระหว่างการศึกษาการสัมผัสและการศึกษาผลลัพธ์ (รายการที่ 10d)

การสุ่มแบบเมนเดลสามารถพิสูจน์ความสัมพันธ์เชิงสาเหตุได้หรือไม่?

การทดสอบแบบสุ่มเมนเดล (Mendelian randomization) เพียงอย่างเดียวไม่เพียงพอ การทดสอบนี้ใช้สมมติฐานเชิงสาเหตุภายใต้ข้อสมมติเฉพาะบางประการ และข้อสมมติเหล่านั้นบางส่วนไม่สามารถตรวจสอบได้อย่างครบถ้วน จึงควรใช้ภาษาที่แสดงถึงสาเหตุอย่างระมัดระวัง รายงานการวิเคราะห์ความไว และถือว่าผลลัพธ์เป็นเพียงหลักฐานส่วนหนึ่งในการตรวจสอบความถูกต้องร่วมกับการออกแบบการวิจัยอื่นๆ

อ้างอิง

  1. Davey Smith G, Ebrahim S. 'การสุ่มแบบเมนเดล': ระบาดวิทยาทางพันธุกรรมสามารถช่วยให้เข้าใจปัจจัยด้านสิ่งแวดล้อมที่ก่อให้เกิดโรคได้หรือไม่? วารสารระบาดวิทยานานาชาติ 2003;32(1):1-22. PMID 12689998. doi.org/10.1093/ije/dyg070
  2. Sekula P, Del Greco MF, Pattaro C, Kottgen A. การสุ่มแบบเมนเดลเป็นแนวทางในการประเมินความสัมพันธ์เชิงสาเหตุโดยใช้ข้อมูลจากการสังเกต วารสารสมาคมโรคไตแห่งอเมริกา 2016;27(11):3253-3265. PMID 27486138. doi.org/10.1681/ASN.2016010098
  3. Skrivankova VW, Richmond RC, Woolf BAR และคณะ การเสริมสร้างความเข้มแข็งในการรายงานการศึกษาเชิงสังเกตในระบาดวิทยาโดยใช้การสุ่มแบบเมนเดล (STROBE-MR): คำอธิบายและการขยายความ BMJ 2021;375:n2233. การเข้าถึงแบบเปิด, CC BY 4.0. doi.org/10.1136/bmj.n2233
  4. Bowden J, Davey Smith G, Haycock PC, Burgess S. การประมาณค่าที่สอดคล้องกันในการสุ่มแบบเมนเดลด้วยเครื่องมือที่ไม่ถูกต้องบางส่วนโดยใช้ตัวประมาณค่ามัธยฐานแบบถ่วงน้ำหนัก Genet Epidemiol. 2016;40(4):304-314. การเข้าถึงแบบเปิด, CC BY 4.0; อัตราส่วนความน่าจะเป็น ช่วงความเชื่อมั่น และค่า p คำนวณใหม่โดย Editverse ใน R จากตารางที่ 4 doi.org/10.1002/gepi.21965
  5. Sanderson E, Glymour MM, Holmes MV และคณะ การสุ่มแบบเมนเดล คู่มือเบื้องต้นเกี่ยวกับวิธีการ Nat Rev 2022;2:6. PMID 37325194. doi.org/10.1038/s43586-021-00092-5
  6. Burgess S, Davey Smith G, Davies NM และคณะ แนวทางการดำเนินการวิจัยแบบสุ่มทางพันธุกรรม: ฉบับปรับปรุงสำหรับฤดูร้อนปี 2023 Wellcome เปิด Res. 2019;4:186 (เวอร์ชัน 3). การเข้าถึงแบบเปิด, CC BY 4.0. doi.org/10.12688/wellcomeopenres.15555.3
  7. Burgess S, Thompson SG. การหลีกเลี่ยงอคติจากเครื่องมือที่ไม่แข็งแรงในการศึกษาแบบสุ่มทางพันธุกรรม (Mendelian randomization) วารสารระบาดวิทยานานาชาติ 2011;40(3):755-764. PMID 21414999. doi.org/10.1093/ije/dyr036
  8. Burgess S, Butterworth A, Thompson SG. การวิเคราะห์แบบสุ่มทางพันธุกรรมด้วยตัวแปรทางพันธุกรรมหลายตัวโดยใช้ข้อมูลสรุป Genet Epidemiol. 2013;37(7):658-665. PMID 24114802. doi.org/10.1002/gepi.21758
  9. Burgess S, Davies NM, Thompson SG. อคติเนื่องจากการทับซ้อนของผู้เข้าร่วมในการสุ่มแบบเมนเดลสองกลุ่มตัวอย่าง Genet Epidemiol. 2016;40(7):597-608. การเข้าถึงแบบเปิด, CC BY 4.0; ตรวจสอบความสัมพันธ์เชิงเส้นโดย Editverse ใน R. doi.org/10.1002/gepi.21998
  10. Bowden J, Davey Smith G, Burgess S. การสุ่มแบบเมนเดลด้วยเครื่องมือที่ไม่ถูกต้อง: การประมาณผลกระทบและการตรวจจับอคติผ่านการถดถอยของ Egger วารสารระบาดวิทยานานาชาติ 2015;44(2):512-525. PMID 26050253. doi.org/10.1093/ije/dyv080
  11. Sanderson E. การสุ่มแบบเมนเดลหลายตัวแปรและการไกล่เกลี่ย Cold Spring Harb Perspect Med 2021;11(2):a038984. PMID 32341063. doi.org/10.1101/cshperspect.a038984
  12. Verbanck M, Chen CY, Neale B, Do R. การตรวจจับพลีโอโทรปีแนวนอนที่แพร่หลายในความสัมพันธ์เชิงสาเหตุที่อนุมานได้จากการสุ่มแบบเมนเดลระหว่างลักษณะที่ซับซ้อนและโรคต่างๆ Nat Genet 2018;50(5):693-698. PMID 29686387. doi.org/10.1038/s41588-018-0099-7
  13. Ference BA, Yoo W, Alesh I และคณะ ผลกระทบของการได้รับคอเลสเตอรอลชนิดไลโปโปรตีนความหนาแน่นต่ำในระดับต่ำเป็นระยะเวลานานตั้งแต่ช่วงต้นของชีวิตต่อความเสี่ยงของโรคหลอดเลือดหัวใจ: การวิเคราะห์แบบสุ่มทางพันธุกรรม (Mendelian randomization analysis) J Am Coll Cardiol 2012;60(25):2631-2639. PMID 23083789. doi.org/10.1016/j.jacc.2012.09.017
  14. Gagnon E, Gill D, Bourgault J และคณะ การแทรกแซง RNA เทียบกับการยับยั้ง PCSK9 โดยใช้แอนติบอดีเพื่อป้องกันโรคหัวใจและหลอดเลือด: การศึกษาแบบสุ่มทางพันธุกรรมเป้าหมายยา Cardiovasc Res 2025;121(7):1066-1075. PMID 40347490. doi.org/10.1093/cvr/cvaf078
  15. Fang A, Zhao Y, Yang P, Zhang X, Giovannucci EL. วิตามินดีและสุขภาพของมนุษย์: หลักฐานจากการศึกษาแบบสุ่มทางพันธุกรรม (Mendelian randomization studies) เออร์ เจ เอพิเดมิออล. 2024;39(5):467-490. PMID 38214845. doi.org/10.1007/s10654-023-01075-4
  16. Zhao SS, Mason A, Gjekmarkaj E, Yanaoka H, ​​Burgess S. ความสัมพันธ์ระหว่างวิตามินดีและโรคภูมิต้านตนเอง: การวิเคราะห์แบบสุ่มทางพันธุกรรม (Mendelian randomization analysis) โรคข้ออักเสบเซมิน 2023;62:152238. PMID 37437450. doi.org/10.1016/j.semarthrit.2023.152238
  17. Millwood IY, Walters RG, Mei XW และคณะ หลักฐานแบบดั้งเดิมและทางพันธุกรรมเกี่ยวกับสาเหตุของโรคหลอดเลือดที่เกิดจากแอลกอฮอล์: การศึกษาเชิงพยากรณ์ในผู้ชายและผู้หญิง 500,000 คนในประเทศจีน มีดหมอ 2019;393(10183):1831-1842. การเข้าถึงแบบเปิด, CC BY 4.0 doi.org/10.1016/S0140-6736(18)31772-0
  18. Larsson SC, Butterworth AS, Burgess S. การสุ่มแบบเมนเดลสำหรับโรคหัวใจและหลอดเลือด: หลักการและการประยุกต์ใช้ Eur Heart J. 2023;44(47):4913-4924. PMID 37935836. doi.org/10.1093/eurheartj/ehad736
  19. Skrivankova VW, Richmond RC, Woolf BAR และคณะ การเสริมสร้างความเข้มแข็งในการรายงานการศึกษาเชิงสังเกตในระบาดวิทยาโดยใช้การสุ่มแบบเมนเดล: คำแถลง STROBE-MR JAMA 2021;326(16):1614-1621. PMID 34698778. doi.org/10.1001/jama.2021.18236
  20. หน่วยระบาดวิทยาเชิงบูรณาการ MRC. OpenGWAS (หน้าหลัก อ่านเมื่อวันที่ 4 พฤศจิกายน 2026) gwas.mrcieu.ac.uk/
  21. FinnGen. หน้าหลักโครงการ (อ่านเมื่อ 4 พฤศจิกายน 2026) www.finngen.fi/en
  22. แคตตาล็อก NHGRI-EBI GWAS หน้าหลัก (อ่านเมื่อ 4 พฤศจิกายน 2026) www.ebi.ac.uk/gwas/
  23. หน่วยระบาดวิทยาเชิงบูรณาการ MRC เอกสารประกอบการใช้งาน TwoSampleMR (เวอร์ชัน 0.7.13 ณ วันที่ 4 พฤศจิกายน 2026) mrcieu.github.io/TwoSampleMR/
  24. CRAN. หน้าแพ็กเกจ MendelianRandomization (เวอร์ชัน 0.10.0 เผยแพร่เมื่อ 12 เมษายน 2024 อ่านเมื่อ 4 พฤศจิกายน 2026) cran.r-project.org/package=MendelianRandomization
  25. Lawlor DA, Tilling K, Davey Smith G. การใช้สามเหลี่ยมในการวิเคราะห์ระบาดวิทยาเชิงสาเหตุ วารสารระบาดวิทยานานาชาติ 2016;45(6):1866-1886. PMID 28108528. doi.org/10.1093/ije/dyw314

การอ่านที่เกี่ยวข้อง

ทีมวิจัย Editverse | อัปเดต พฤศจิกายน 2026