Sistemik lupus eritematozus immun sisteminin insanın öz hüceyrə və toxumalarına qarşı reaksiya yaratdığı, genetik meyilliliklə ətraf mühit təsirlərinin birlikdə rol oynadığı mürəkkəb autoimmun xəstəlikdir. Xəstəliyin molekulyar xüsusiyyətlərindən biri tip I interferon siqnalının bir çox xəstədə həddindən artıq aktivləşməsidir. Lakin genom miqyaslı assosiasiya tədqiqatlarında lupusla əlaqəli tapılan genetik variantların əksəriyyəti zülal kodlayan bölgələrdə deyil. Buna görə bir variantın xəstəliklə statistik əlaqə göstərməsi onun hansı geni, hansı hüceyrədə və hansı molekulyar proses vasitəsilə təsir etdiyini təkbaşına izah etmir.

Bu tədqiqatda araşdırmaçılar Taiwan Precision Medicine Initiative çərçivəsində əldə edilən sistemik lupus GWAS xülasə statistikalarını müxtəlif funksional genomika mənbələri ilə birləşdirdilər. Analiz 5.184 lupus halı və 311.848 nəzarətdən əldə olunan GWAS nəticələrinə əsaslanırdı. Genom miqyasında əhəmiyyətli variantlardan 21 namizəd gen müəyyən edildi. Bu genlər içində xüsusilə IRF5, IRF8 və STAT4 interferon və sitokin siqnallaşması ilə əlaqəli olmaları və TNFSF10 geninin tənzimləyici bölgələrində özlərinə aid proqnozlaşdırılan bağlanma motivlərinin olması səbəbindən prioritetləşdirildi.

TNFSF10 immun tənzimlənməsi və hüceyrə siqnal prosesləri ilə əlaqələndirilən TRAIL zülalını kodlayır. Araşdırmanın mühüm nəticəsi TNFSF10 ətrafındakı genetik variantların çoxunun genin ümumi ifadə miqdarından daha çox, yaranan RNT-nin necə kəsilib-birləşdirilməsi ilə əlaqəli görünməsidir. Başqa sözlə, variantlar TNFSF10-un yalnız “nə qədər istehsal edildiyini” deyil, hansı RNT formalarının və ya izoformlarının istehsal edildiyini də təsirləndirə bilər.

GTEx məlumatlarında TNFSF10-un aşağı axınında və 3′ ucuna yaxın yerləşən çoxsaylı variantlar, xüsusilə EBV ilə transformasiya edilmiş limfositlərdə və tam qanda splicing QTL siqnalları göstərdi. EBV ilə transformasiya edilmiş limfositlərdə orta sQTL təsir ölçüsü tam qandakı orta təsirin təxminən 2,35 qatı idi. Variantların bir hissəsi terminal intronlara, ekzonlara və 3′UTR bölgəsinə yaxın olduğu halda, bir hissəsi genin bir neçə min baza cütü kənarında yerləşirdi.

Linkage disequilibrium analizləri bu variantların bir-birindən müstəqil hərəkət etmədiyini və bəzi populyasiyalarda güclü genetik bloklar yaratdığını göstərdi. Tayvan və Şərqi Asiya populyasiyalarında görünən əlaqə quruluşu Afrika, Amerika, Avropa və Cənubi Asiya populyasiyalarındakı quruluşlardan fərqli idi. Bu vəziyyət eyni bioloji bölgənin müxtəlif icmalarda fərqli marker variantlarla təmsil oluna biləcəyini göstərir.

ENCODE, ORegAnno, ReMap, H3K27ac və Micro-C məlumatları da TNFSF10-un aşağı axınındakı bəzi variantların tənzimləyici bölgələrlə və üçölçülü xromatin təmasları ilə üst-üstə düşdüyünü göstərdi. Xüsusilə rs62282578 TNFSF10-un 3′ ucundan təxminən 1.924 baza cütü aşağıda yerləşən və enhancer kimi işarələnən tənzimləyici bölgə ilə üst-üstə düşürdü. Bununla yanaşı, bu nəticələrin heç biri variantların TNFSF10-u birbaşa tənzimlədiyini və ya lupusa səbəb olduğunu sübut etmir. Tədqiqat eksperimental şəkildə doğrulanmalı namizəd variantları və mümkün mexanizmləri prioritetləşdirən hesablama yönümlü araşdırmadır.

Ətraflı İzah

Sistemik lupusda genetik mexanizmi anlamaq niyə çətindir?

Sistemik lupus eritematozus dəri, oynaq, böyrək, qan hüceyrələri, sinir sistemi və başqa orqanlara təsir edə bilən, klinik görünüşü insandan insana böyük dərəcədə dəyişən xroniki autoimmun xəstəlikdir. Xəstəliyin yaranmasını tək genin pozulması ilə izah etmək mümkün deyil. Çoxsaylı genetik variant, ətraf mühit tətikləyiciləri və immun sisteminin tənzimləyici mexanizmləri birlikdə rol oynayır.

Tədqiqatda ultrabənövşəyi şüalanma və virus infeksiyaları genetik meyilliliklə qarşılıqlı təsir göstərə bilən ətraf mühit amillərinə nümunə kimi verilir. Bu cür təsirlər anadangəlmə və qazanılmış immun hüceyrələri hərəkətə gətirən interferon sisteminə təsir edə bilər.

Tip I interferon siqnalının həddindən artıq işləməsi sistemik lupusun aydın molekulyar xüsusiyyətlərindən biridir. İnterferon siqnalları monositlərə, dendritik hüceyrələrə, makrofaqlara, T hüceyrələrinə və B hüceyrələrinə təsir edə bilər. Buna görə interferon istehsalını və ya interferona hüceyrə cavabını dəyişdirən irsi variantlar lupus meyilliliyinin anlaşılması üçün mühüm namizədlərdir.

Bununla yanaşı genom miqyaslı assosiasiya tədqiqatlarının tapdığı siqnalların əksəriyyəti genlərin zülal kodlayan hissələrində deyil. GWAS siqnalı çox vaxt aşağıdakı sualların heç birinə birbaşa cavab vermir:

• Statistik əlaqəni yaradan həqiqi funksional variant hansıdır?
• Variant hansı geni tənzimləyir?
• Genin ümumi ifadə səviyyəsini, yoxsa RNT işlənməsini dəyişdirir?
• Təsir hansı immun hüceyrəsində və ya toxumada ortaya çıxır?
• Eyni variant fərqli mənşədən gələn insanlarda eyni genetik quruluşun hissəsidirmi?
• Əlaqə birbaşa xarakter daşıyır, yoxsa yaxınlıqdakı başqa variantı təmsil edir?

Buna görə GWAS-dan sonra tapılan statistik siqnalların funksional genomika məlumatları ilə izah edilməsi lazımdır. Bu proses ümumilikdə post-GWAS analizi adlanır.

Araşdırmanın əsas sualı

Tədqiqatın mərkəzində TNFSF10 geni dayanır. TNFSF10 şiş nekrozu faktoru ilə əlaqəli apoptoz induksiya edən ligand kimi tanınan TRAIL zülalını kodlayır. TRAIL yalnız hüceyrə ölümünü tənzimləyən molekul deyil; T hüceyrə aktivləşməsi, iltihab və immun tolerantlıqla əlaqəli kontekstlərdə də rol oynaya bilər.

Əvvəlki araşdırmalarda lupus xəstələrində TNFSF10 ifadəsinin və dövr edən TRAIL səviyyələrinin dəyişə bildiyi, bəzi nəticələrin xəstəlik aktivliyi və ya orqan tutulması ilə əlaqəli olduğu bildirilmişdir. İnterferon siqnalının TNFSF10 ifadəsini artıra bildiyinə dair tapıntılar da var.

Buna baxmayaraq interferonla əlaqəli lupus risk genlərinin TNFSF10-un tənzimlənməsi ilə hansı genetik arxitektura daxilində birləşdiyi aydın deyil. Buna görə tədqiqatçılar üçqat sual qoydular:

1. Lupus GWAS siqnalları içində hansı interferon genləri ön plana çıxır?
2. Bu genlər TNFSF10-un tənzimləyici bölgələri ilə potensial əlaqə daşıyırmı?
3. TNFSF10 ətrafındakı variantlar gen ifadəsindən daha çox RNT kəsilib-birləşdirilməsinə təsir göstərə bilərmi?

GWAS məlumat mənbəyi və tədqiqatın ölçüsü

GWAS xülasə statistikaları Tayvanda 16 tibb mərkəzindən iştirakçıların yer aldığı Taiwan Precision Medicine Initiative məlumat mənbəyindən əldə edildi. TPMI proqramında ümumilikdə 565.390 iştirakçının olduğu bildirilir.

Sistemik lupus xüsusiyyəti üçün istifadə olunan GWAS məlumat dəsti:

• 5.184 sistemik lupus halından,
• 311.848 nəzarətdən

ibarət idi. Hallar ICD-9 710.0 və ICD-10 M32 kodları vasitəsilə müəyyən edildi. Nəzarət qrupunda lupusla birlikdə seçilmiş autoimmun, iltihabi və dermatoloji diaqnoz kodlarına sahib şəxslər xaric edildi.

Genetik profillər Taiwan Biobank tərəfindən hazırlanmış TPMI-yə xüsusi SNP array-ləri ilə yaradıldı. Variant-xəstəlik əlaqələri böyük biobank tədqiqatlarında nümunə qohumluğunu və qeyri-bərabər hal-nəzarət nisbətlərini nəzərə almaq üçün istifadə olunan SAIGE üsulu ilə analiz edildi.

Genom miqyasında əhəmiyyətlilik nə deməkdir?

Tədqiqatçılar namizəd variantları seçərkən standart genom miqyasında əhəmiyyətlilik həddindən istifadə etdilər:

\[ P < 5 \times 10^{-8} \]

Bu hədd yüz minlərlə və ya milyonlarla genetik variantın eyni anda sınaqdan keçirildiyi GWAS tədqiqatlarında təsadüfi müsbət nəticələrin sayını azaltmaq üçün istifadə olunan çox sərt statistik həddir.

Variantın bu həddi keçməsi onun xəstəliyə birbaşa səbəb olduğunu göstərmir. Variant həqiqi funksional variantla güclü linkage disequilibrium içində olan genetik marker ola bilər. Post-GWAS analizinin məqsədi bu geniş siqnalı daha dar və bioloji baxımdan mənalı namizədlərə çevirməkdir.

Manhattan və Q-Q qrafikləri nə göstərir?

Genom miqyaslı assosiasiya nəticələrini göstərən Manhattan qrafikində xromosomlar boyunca müxtəlif yüksək siqnallar görünür. IRF5, STAT4, IRF8, HLA-DQA1, TNFAIP3, TNFSF4, ATG5, BLK, CDKN1B və başqa immun genlərinə yaxın bölgələr ön plana çıxır.

Qrafikdə hər nöqtə genetik variantı, şaquli ox isə P dəyərinin mənfi loqarifmini təmsil edir. Nöqtə nə qədər yuxarıdadırsa, xəstəliklə statistik əlaqə o qədər güclüdür. Lakin yüksək zirvə təkbaşına hansı variantın funksional olduğunu göstərmir; eyni əlaqə blokundakı çoxsaylı variant birlikdə yüksəlmiş ola bilər.

Q-Q qrafikində əhəmiyyətli variantların gözlənilən sıfır paylanmasından ayrıldığı, genomik inflyasiya faktorunun isə 1-ə yaxın olduğu bildirildi. Tədqiqatçılar bunu geniş və sistematik statistik şişmənin məhdud olduğu kimi şərh etdilər.

İyirmi bir namizəd genin ortaq xüsusiyyəti nə idi?

Genom miqyasında əhəmiyyətli variantlardan ümumilikdə 21 namizəd gen uyğunlaşdırıldı. Bu genlər DAVID istifadə olunaraq Gene Ontology bioloji prosesləri və Reactome yolları baxımından araşdırıldı.

Düzəldilmiş P dəyəri 0,05-dən aşağı olan zənginləşmiş başlıqlar arasında xüsusilə:

• Tip I interferon istehsalının tənzimlənməsi,
• İmmun sistemində sitokin siqnallaşması,
• İnterferonla əlaqəli anadangəlmə immun yolları,
• DDX58 və IFIH1 siqnalının mənfi tənzimləyiciləri

yer aldı.

Bu qruplara töhfə verən genlər arasında IRF5, IRF8, ISG15, STAT4, SH2B3, HLA-DQA1, TNFSF4, CDKN1B, ATG5 və TNFAIP3 var idi.

Enrichr-KG ilə yaradılan funksional şəbəkə də oxşar quruluş göstərdi. Şəbəkə düyünlərinin mühüm hissəsi interferon və sitokin tənzimlənməsi ətrafında birləşdi. Bu uyğunluq namizəd genlərin təsadüfi və bir-biri ilə əlaqəsiz siyahı yaratmadığını, lupus biologiyasının məlum interferon şəbəkəsində sıxlaşdığını düşündürdü.

Toxuma ifadə xəritəsi niyə vacib idi?

GTEx RNT ifadə profilləri araşdırıldıqda bir çox namizəd genin və TNFSF10-un immun baxımdan zəngin nümunələrdə nisbətən yüksək ifadə səviyyəsinə sahib olduğu görüldü. Xüsusilə:

• EBV ilə transformasiya edilmiş limfositlər,
• Tam qan,
• Dalaq və immunitetlə əlaqəli bəzi toxumalar

ön plana çıxdı.

Bu müşahidə variantların yalnız genomik mövqe baxımından deyil, lupus baxımından əlaqəli ola biləcək hüceyrə mühitlərində ifadə olunan genlərlə bağlı olduğunu dəstəkləyir. Lakin toxuma səviyyəsində ifadə xəritəsi konkret immun hüceyrə alt qrupunda tənzimləyici təsirin olduğunu sübut etmir.

IRF5, IRF8 və STAT4 niyə prioritetləşdirildi?

Funksional zənginləşdirmə nəticələrindən sonra tədqiqatçılar interferonla əlaqəli namizəd genlərin TNFSF10-un tənzimləyici arxitekturasına bağlanıb-bağlanmadığını araşdırdılar.

Transkripsiya faktoru motiv analizində:

• IRF5 üçün TNFSF10-un 5′ intronik bölgəsində,
• IRF8 üçün TNFSF10-un yuxarı axın tənzimləyici bölgəsində iki ayrı mövqedə,
• STAT4 üçün TNFSF10-un yuxarı axın bölgəsində

proqnozlaşdırılan motivlər tapıldı.

Mətndə verilən koordinatlar belədir:

• IRF5: chr3:172.520.205-172.520.249
• IRF8: chr3:172.523.606-172.523.619
• IRF8: chr3:172.524.863-172.524.876
• STAT4: chr3:172.527.465-172.527.474

Motiv şəklinin IRF5 alt panelində isə ilk aralığın son koordinatı 172.520.219 kimi yazılmışdır. Buna görə preprintin əsas mətni ilə şəkil etiketi arasında IRF5 motiv koordinatı baxımından fərq var. Bu vəziyyət motivin bioloji mövcudluğunu təkbaşına etibarsız etmir; lakin dəqiq koordinatın orijinal analiz faylları və ya sonrakı versiya vasitəsilə yoxlanması lazımdır.

Motivin tapılması müvafiq transkripsiya faktorunun həqiqətən həmin bölgəyə bağlandığını göstərmir. Motiv DNT ardıcıllığının IRF5, IRF8 və ya STAT4 tərəfindən tanına biləcək nümunəyə bənzədiyini göstərir.

Birbaşa bağlanmanı göstərmək üçün aşağıdakı tip təcrübələr lazımdır:

• ChIP-seq və ya ChIP-qPCR,
• DNT bağlanma təcrübələri,
• Reportyor gen təcrübələri,
• CRISPR əsaslı bölgə silmə və ya dəyişdirmə,
• IRF5, IRF8 və ya STAT4 basqılama və artırma təcrübələri.

Üç əsas GWAS variantı

IRF5, IRF8 və STAT4 ilə əlaqəli üç aparıcı variant ətraflı araşdırıldı.

Beta işarələrinin fərqli olması alternativ allellərin GWAS xüsusiyyəti ilə eyni istiqamətdə əlaqə göstərmədiyini ortaya qoyur. Müsbət beta analizdə alternativ allelin lupus xüsusiyyəti ilə artım istiqamətində, mənfi beta isə azalma istiqamətində əlaqəli olduğunu göstərir. Lakin bu işarələr birbaşa bioloji mexanizm, gen aktivləşməsi və ya klinik qoruyuculuq demək deyil.

Üç variant arasında yalnız rs3757387 üçün GTEx eQTL qeydləri tapıldı. Təsirlər günəşə məruz qalan dəri, EBV ilə transformasiya edilmiş limfositlər və tam qan kimi lupus baxımından əlaqəli nümunələrdə müşahidə edildi.

rs35703946 IRF8 ilə LINC01082 arasında yerləşən genlərarası variantdır. Tədqiqatçılar bu mövqenin IRF8-dən uzaq görünməsinə baxmayaraq uzun məsafəli xromatin qarşılıqlı təsirləri vasitəsilə tənzimləyici potensial daşıya biləcəyini düşünürlər. Bu ehtimal birbaşa təcrübə ilə göstərilməmişdir.

rs7568275-in STAT4-ün üçüncü intronunda yerləşməsi gen ifadəsi ilə yanaşı RNT işlənməsi ilə əlaqəli funksiya ehtimalını gündəmə gətirir. Lakin tədqiqat rs7568275-in STAT4 kəsilib-birləşdirilməsini birbaşa dəyişdirdiyini göstərməmişdir.

Allel tezlikləri populyasiyalar arasında niyə dəyişir?

Eyni variantın tezliyi müxtəlif coğrafi və genetik mənşəyə sahib populyasiyalarda dəyişə bilər. Bu vəziyyət variantın bir icmada vacib, başqa icmada isə əhəmiyyətsiz olması demək deyil. Xəstəliklə əlaqəli funksional bölgə eyni qala bilər, lakin həmin bölgəni təmsil edən yayılmış marker variant dəyişə bilər.

Tədqiqatda rs3757387-nin alternativ allel tezliyi:

• TPMI-də 0,151,
• Taiwan Biobank-da 0,143,
• 1000 Genomes Şərqi Asiya qrupunda 0,107,
• Avropa qrupunda 0,454,
• Amerika qrupunda 0,324

kimi bildirildi.

rs35703946-nın Tayvan tezliyi təxminən 0,095-0,096 idi. Tədqiqatçılar bunun Şərqi Asiya referensindən daha çox bəzi Amerika populyasiyası təxminlərinə yaxın göründüyünü qeyd etdilər.

rs7568275 isə Tayvanda təxminən 0,643-0,645 tezliyinə sahib idi; bəzi Afrika və Amerika mənbələrində daha yüksək tezliklər görünürdü.

Bu cür fərqlər dörd səbəbdən qaynaqlana bilər:

• Populyasiyaların tarixi ayrılması,
• Genetik drift,
• Founder effektləri və bottleneck-lər,
• Eyni bölgədə fərqli haplotiplərin yayılması.

Linkage disequilibrium nədir?

Yaxın genomik bölgələrdəki variantlar nəsillər boyunca birlikdə ötürülə bilər. İki variantın birlikdə irsən keçmə dərəcəsi linkage disequilibrium, yəni linkage disequilibrium ilə ifadə olunur.

Tədqiqatda istifadə olunan əsas ölçü r² dəyəridir:

\[ 0 \leq r^2 \leq 1 \]

• r² dəyəri 0-a yaxınlaşdıqda iki variant böyük ölçüdə müstəqildir.
• r² dəyəri 1-ə yaxınlaşdıqda variantlar populyasiyada güclü şəkildə birlikdə daşınır.
• Tədqiqatda r² > 0,8 dəyərləri yüksək linkage disequilibrium kimi şərh edilmişdir.

Yüksək linkage disequilibrium əlaqə göstərən SNP-nin funksional variant olduğu demək deyil. Eyni blok daxilində çoxsaylı variant bir-birinə oxşar statistik siqnal göstərə bilər. Həqiqi funksional variantı müəyyən etmək üçün əlavə funksional sübut lazımdır.

IRF5, IRF8 və STAT4 bölgələrinin LocusZoom nəticələri

LocusZoom qrafiklərində aparıcı variantların ətrafındakı digər SNP-lər assosiasiya gücü və aparıcı variantla əlaqə dərəcəsinə görə göstərildi.

rs3757387 və rs7568275-in yaxınlıqdakı müxtəlif variantlarla birləşdirilmiş populyasiya analizində r² > 0,8 səviyyəsində əlaqə göstərdiyi bildirildi. Bununla birlikdə eyni bölgələrin populyasiyaya xas LD xəritələri fərqli nümunələr yaratdı.

Bu vəziyyət iki mühüm nəticə doğurur:

1. Bir populyasiyada aparıcı variant kimi görünən SNP başqa populyasiyada həqiqi funksional variantı daha zəif təmsil edə bilər.
2. Fərqli populyasiyaların birlikdə araşdırılması geniş əlaqə bloklarını parçalayaraq namizəd səbəbsəl variantların daraldılmasına kömək edə bilər.

Güc analizi nəyi göstərdi, nəyi göstərmədi?

GAS Power Calculator-un yüksək dərəcədə qeyri-bərabər hal-nəzarət dizaynlarındakı giriş məhdudiyyətləri səbəbindən tədqiqatçılar həqiqi TPMI nümunə ölçüsünü sözbəsöz yenidən hesablamadılar. Bunun əvəzinə 1.662 halı sabit saxlayaraq müxtəlif nəzarət sayları ilə həssaslıq analizi apardılar:

• 1.662 hal / 1.662 nəzarət,
• 1.662 hal / 3.324 nəzarət,
• 1.662 hal / 16.620 nəzarət,
• 1.662 hal / 100.000 nəzarət.

Xəstəlik prevalansı 0,0007, genom miqyasında əhəmiyyətlilik həddi 5 × 10^-8 və multiplikativ xəstəlik modeli istifadə edildi. Ən yüksək nəzarət sayına sahib ssenaridə üç variant üçün gücün təxminən yüzdə 80 səviyyəsinə çata bildiyi bildirildi.

Bu hesablama TPMI GWAS-ının həqiqi və dəqiq statistik gücü kimi oxunmamalıdır. Tədqiqatçılar da analizi nümunə ölçüsü replikasiyası deyil, nəzarət sayındakı artımın mümkün təsirini göstərən həssaslıq qiymətləndirməsi kimi təsvir edirlər.

eQTL ilə sQTL arasındakı fərq

Tədqiqatın ən mühüm nöqtələrindən biri genetik variantların gen ifadəsini və RNT işlənməsini fərqli formalarda təsirləndirə bilməsidir.

Tədqiqatın məntiqini izah etmək üçün əsas əlaqələr belə göstərilə bilər:

\[ \mathrm{eQTL}: \ Genetik\ varyant \rightarrow Toplam\ gen\ ifadesindeki\ değişim \]

\[ \mathrm{sQTL}: \ Genetik\ varyant \rightarrow RNA\ kesilip\ birleştirilmesindeki\ değişim \]

Bu göstərimlər məqalədə yayımlanmış riyazi tənliklər deyil; eQTL və sQTL anlayışlarının fərqini izah etmək məqsədilə verilmişdir.

eQTL müəyyən genetik variantın bir genin ümumi RNT miqdarı ilə əlaqəsini göstərir. Məsələn, bir alleli daşıyanlarda TNFSF10 RNT-si daha çox və ya daha az ola bilər.

sQTL isə ümumi RNT miqdarından daha çox RNT-nin hansı hissələrinin saxlanılacağı, hansı intronların çıxarılacağı və hansı ekzonların birləşdiriləcəyi ilə əlaqəlidir. Eyni gen müxtəlif kəsilib-birləşdirmə formaları ilə fərqli RNT izoformları yarada bilər.

Bu fərq fabrikin ümumi istehsalı ilə məhsul çeşidi arasındakı fərqə bənzədilə bilər. eQTL fabrikin neçə məhsul istehsal etdiyinə, sQTL isə eyni xammaldan hansı məhsul modelinin istehsal edildiyinə təsir göstərə bilər.

TNFSF10 bölgəsində niyə sQTL siqnalları önə çıxdı?

GTEx QTL məlumatları araşdırıldıqda TNFSF10 bölgəsindəki sQTL siqnallarının eQTL siqnallarından daha çox olduğu görüldü.

EBV ilə transformasiya edilmiş limfositlərdə yalnız bəzi variantlar TNFSF10-un ümumi ifadəsi üzərində mülayim mənfi eQTL təsirləri göstərdi. Bunlar:

• rs10936723,
• rs3136599,
• rs3181139,
• rs2241063,
• rs2270418,
• rs3136586

kimi siyahıya alındı. Təsir ölçüləri təxminən 0,2 səviyyəsində idi.

Bunun əvəzində daha çox variant TNFSF10 kəsilib-birləşdirilməsi üzərində müsbət sQTL təsirləri göstərdi. Siqnallar xüsusilə:

• EBV ilə transformasiya edilmiş limfositlərdə,
• Tam qanda,
• Ağciyərdə,
• Ürəkdə,
• Dəri toxumalarında

müşahidə edildi.

EBV ilə transformasiya edilmiş limfositlər və tam qan arasında ortaq sQTL variant bloku chr3:172.474.449 ilə chr3:172.522.815 arasında yerləşirdi.

Bu variantların mühüm hissəsi TNFSF10-un aşağı axınında, terminal intronlarında və ya 3′UTR ətrafında yerləşirdi. EBV ilə transformasiya edilmiş limfositlərdə orta sQTL təsiri tam qandakı orta təsirin təxminən 2,35 qatı idi.

Bu fərq eyni DNT variantının fərqli hüceyrə mühitlərində eyni böyüklükdə funksiya göstərməyə biləcəyini düşündürür. Limfositlərdə olan transkripsiya faktorları, RNT bağlayan zülallar və xromatin vəziyyəti TNFSF10 RNT işlənməsinə tam qandakı qarışıq hüceyrə populyasiyasından fərqli təsir edə bilər.

3′UTR və terminal intronlar niyə vacibdir?

Genlərin 3′ uc bölgəsi zülal kodlanması tamamlandıqdan sonra funksiyasız uzantı deyil. 3′UTR:

• RNT sabitliyinə,
• RNT-nin hüceyrə daxilində daşınmasına,
• MikroRNT bağlanmasına,
• Zülal istehsal sürətinə,
• Alternativ poliadenilasiyaya

təsir edə bilər.

Terminal intronlar və son ekzon ətrafı da RNT kəsilib-birləşdirilməsi baxımından vacibdir. Bu bölgədəki variant:

• Ekzonun atlanmasına,
• İntronun RNT daxilində saxlanılmasına,
• Alternativ kəsim nöqtəsinin istifadəsinə,
• Fərqli uzunluqda transkriptlərin yaranmasına

təsir göstərə bilər.

Araşdırmada TNFSF10 ətrafındakı sQTL variantları iki əsas qrupa bölündü. Bir qrup genin 3′ ucuna yaxın intronik, ekzonik və 3′UTR variantlarından, digər qrup isə TNFSF10-dan 4.000 baza cütündən daha uzaqda yerləşən aşağı axın variantlarından ibarət idi.

İki alt qrup arasında r² > 0,8 səviyyəsində güclü əlaqənin olması fərqli mövqelərdəki variantların eyni irsi haplotip blokunun hissələri ola biləcəyini göstərir. Buna görə sQTL siqnalını yaradan həqiqi funksional variant siyahıya alınmış SNP-lərdən biri və ya bu SNP-lərlə əlaqəli başqa variant ola bilər.

Əvvəlki TNFSF10 variantları ilə əlaqə

Tədqiqatçılar əvvəlki Misir və İran lupus tədqiqatlarında bildirilmiş rs1131568 və rs1131580 variantlarını da müzakirə etdilər. Bu variantlar TNFSF10-un 3′UTR bölgəsində yerləşir və lupus meyilliliyi və ya klinik xüsusiyyətlərlə əlaqələndirilmişdir.

GTEx analizində hər iki variant da EBV ilə transformasiya edilmiş limfositlərdə və tam qanda güclü TNFSF10 sQTL təsirləri göstərdi. Bu vəziyyət əvvəllər xəstəliklə əlaqələndirilən 3′UTR variantlarının ümumi gen ifadəsindən daha çox transkript işlənməsi vasitəsilə funksiya göstərə biləcəyi hipotezini dəstəkləyir.

Lakin mövcud tədqiqat bu variantların lupus riskini RNT kəsilib-birləşdirilməsi vasitəsilə artırdığını eksperimental şəkildə göstərməmişdir.

Populyasiyaya xas TNFSF10 əlaqə blokları

Taiwan Biobank məlumatları ilə yaradılan TNFSF10 LD xəritəsində variantların əksəriyyəti bir-biri ilə yüksək əlaqə göstərdi. Bəzi variantlar bu ümumi blokdan daha müstəqil görünürdü:

• rs2241063,
• rs3136594,
• rs3821821,
• rs3774313,
• rs3774315,
• rs1823227,
• rs2270417.

Şərqi Asiya populyasiyasında Tayvana bənzər əlaqə bloku görüldü. Afrika, Amerika, Avropa və Cənubi Asiya populyasiyalarında isə blok sərhədləri və variantlar arasındakı əlaqə dərəcələri dəyişirdi.

Tədqiqatçılar Şərqi Asiya və Cənubi Asiya qruplarında aşağı axın variantları ilə transkript bölgəsindəki variantlar arasında bəzi digər populyasiyalara nisbətən daha güclü əlaqə olduğunu qeyd etdilər.

Bu fərqlər variantın klinik test dəyərinin və ya xəstəliklə əlaqəsinin bütün populyasiyalara eyni şəkildə ümumiləşdirilə bilməyəcəyini göstərir. Bununla yanaşı, tədqiqat bu SNP-lərin heç birini klinik genetik biomarker kimi doğrulamamışdır.

Tənzimləyici genom anotasiya məlumatları nə göstərdi?

TNFSF10 ətrafındakı variantlar UCSC Genome Browser vasitəsilə müxtəlif tənzimləyici məlumat qatları ilə müqayisə edildi:

• ENCODE namizəd cis-tənzimləyici elementləri,
• ORegAnno tənzimləyici qeydləri,
• ReMap DNT bağlanma bölgələri,
• H3K27ac işarələri,
• Micro-C üçölçülü xromatin qarşılıqlı təsirləri.

H3K27ac aktiv enhancer və promotor bölgələri ilə tez-tez əlaqələndirilən histon işarəsidir. Variantın H3K27ac bölgəsi ilə üst-üstə düşməsi həmin mövqenin tənzimləyici ola biləcəyini düşündürür, lakin variantın funksiyasını təkbaşına sübut etmir.

ReMap çoxsaylı ChIP-seq təcrübələrindən əldə edilən transkripsiya faktoru və tənzimləyici zülal bağlanma bölgələrini birləşdirir.

ENCODE cCRE qeydləri xromatin əlçatanlığı, histon işarələri və başqa məlumatları birləşdirərək ehtimal olunan promotor və ya enhancer bölgələrini təsnif edir.

TNFSF10-un aşağı axınındakı çoxsaylı SNP-nin bu tənzimləyici qatların biri və ya bir neçəsi ilə üst-üstə düşdüyü görüldü. Üç əsas bölgədə H3K27ac və ReMap siqnallarının sıxlaşması potensial enhancer bölgələrinin mövcudluğunu düşündürdü.

rs62282578 niyə diqqət çəkdi?

rs62282578 TNFSF10-un 3′ ucundan təxminən 1.924 baza cütü aşağıda yerləşən variant kimi ön plana çıxarıldı. Variant ENCODE4 tərəfindən enhancer xüsusiyyəti daşıyan namizəd cis-tənzimləyici element kimi anotasiya edilən bölgə ilə üst-üstə düşürdü.

Mətndə koordinat ifadəsi “chr3:172.505.505 to 172.503.581” şəklində verilmişdir və başlanğıc-son koordinatları tərs ardıcıllıqda görünür. Buna görə bölgənin dəqiq koordinat düzülüşü preprintin sonrakı versiyasında və ya mənbə anotasiya faylında təsdiqlənməlidir. Araşdırmanın əsas iddiası rs62282578-in TNFSF10 3′ ucundan 1.924 baza cütü aşağıda bir enhancer anotasiya ilə üst-üstə düşməsidir.

Bu variantın seçilməsində üç sübut qatı birlikdə rol oynayır:

1. TNFSF10-un aşağı axınında yerləşməsi,
2. Tənzimləyici enhancer anotasiya ilə üst-üstə düşməsi,
3. Terminal intronik bölgələrlə üçölçülü xromatin təması siqnalı göstərməsi.

Bununla yanaşı, tədqiqat rs62282578-i səbəbsəl variant kimi müəyyən etmir. Variant yalnız eksperimental araşdırma üçün prioritet namizədlərdən biridir.

Micro-C xromatin qarşılıqlı təsiri nə deməkdir?

DNT hüceyrə nüvəsində düz xətt kimi dayanmır. Xromosomlar üçölçülü şəkildə qatlanır. Xətti DNT ardıcıllığında bir-birindən minlərlə baza cütü uzaqda olan iki bölgə nüvədə fiziki olaraq bir-birinə yaxınlaşa bilər.

Micro-C üsulu xromatin daxilindəki bu təmasları yüksək rezolyusiyada ölçməyə çalışır. Araşdırmada TNFSF10-un 3′UTR bölgəsinə yaxın aşağı axın SNP-ləri ilə terminal intronlar arasında təmas siqnalları görüldü. Bəzi uzaq aşağı axın bölgələrinin qarşılıqlı təsir skorları 50-dən yuxarı idi.

Bu müşahidə aşağıdakı modeli dəstəkləyə bilər:

\[ Uzak\ düzenleyici\ bölge \rightarrow Üç\ boyutlu\ kromatin\ teması \rightarrow TNFSF10\ RNA\ işlenmesinde\ değişim \]

Bu əlaqə tədqiqatın mexanizmini izah etmək üçün verilmiş konseptual göstərimdir; araşdırmada birbaşa sübut edilmiş səbəbsəl zəncir deyil.

Mühüm məhdudiyyət istifadə olunan Micro-C məlumatının H1 insan embrionik kök hüceyrələrindən gəlməsidir. Bu hüceyrələr lupus xəstələrinin limfositləri, monositləri və ya başqa immun hüceyrələri deyil. Buna görə xromatin təması xəritəsi bölgənin struktur potensialını göstərə bilər, lakin eyni təmasların lupus immun hüceyrələrində baş verdiyi hələ göstərilməmişdir.

Şəkillər birlikdə necə oxunmalıdır?

Tədqiqatın vizual axını tək sübut əvəzinə bir-birini tamamlayan qatlar təqdim edir:

• Manhattan qrafiki: Lupusla genom miqyasında əlaqəli bölgələri göstərir.
• Zənginləşdirmə şəbəkəsi: Namizəd genlərin interferon və sitokin yollarında qruplaşdığını göstərir.
• Motiv xəritəsi: IRF5, IRF8 və STAT4 üçün TNFSF10 ətrafında mümkün DNT bağlanma ardıcıllıqlarını göstərir.
• LocusZoom qrafikləri: Aparıcı variantlarla yaxınlıqdakı SNP-lərin assosiasiya və LD quruluşunu göstərir.
• QTL və LD xəritəsi: TNFSF10 ətrafındakı siqnalların xüsusilə kəsilib-birləşdirmə və əlaqə blokları ilə bağlı olduğunu göstərir.
• Populyasiya LD qrafikləri: Eyni genomik bölgənin müxtəlif mənşə qruplarında fərqli haplotip quruluşları yaratdığını göstərir.
• UCSC və Micro-C xəritəsi: Variantların tənzimləyici anotasiya və üçölçülü xromatin təmasları ilə üst-üstə düşməsini göstərir.
• Son model sxemi: Yuxarı axın interferon genləri, TNFSF10 tənzimlənməsi, aşağı axın sQTL variantları və populyasiyaya xas təsirləri vahid konseptual modeldə birləşdirir.

Bu şəkillərin heç biri təkbaşına birbaşa mexanizm sübutu deyil. Tədqiqatın gücü müxtəlif məlumat qatlarında eyni bölgənin təkrar-təkrar prioritet qazanmasıdır.

Tədqiqatın təklif etdiyi inteqrativ model

Araşdırmanın nəticə modeli üç əsas səviyyədən ibarətdir.

Birinci səviyyə: İnterferon risk genləri

IRF5, IRF8 və STAT4 ətrafındakı lupusla əlaqəli variantlar bu genlərin ifadəsini və ya RNT işlənməsini təsirləndirə bilər. Bu ehtimal rs3757387 üçün eQTL məlumatı ilə qismən dəstəkləndiyi halda, digər variantlar üçün daha çox genomik mövqe və motiv çıxarışına əsaslanır.

İkinci səviyyə: TNFSF10-un transkripsional tənzimlənməsi

IRF5, IRF8 və STAT4 zülallarının tanıya biləcəyi proqnozlaşdırılan DNT motivləri TNFSF10-un yuxarı axın və ya 5′ intronik bölgələrində yerləşir. Bu, interferon siqnalı ilə TNFSF10 arasında mümkün tənzimləmə xətti təqdim edir.

Üçüncü səviyyə: TNFSF10 RNT izoformları

TNFSF10-un terminal intronları, 3′UTR bölgəsi və aşağı axınındakı SNP-lər xüsusilə immunitetlə əlaqəli toxumalarda sQTL siqnalları göstərir. Tənzimləyici anotasiya və xromatin təmasları bəzi uzaq variantların RNT işlənməsinə təsir göstərə biləcək bölgələrlə fiziki və ya funksional əlaqə qura biləcəyini düşündürür.

Bu model belə xülasə edilə bilər:

\[ Lupus\ ile\ ilişkili\ interferon\ varyantları \rightarrow IRF5/IRF8/STAT4\ düzenlenmesi \rightarrow TNFSF10\ düzenleyici\ mimarisi \]

\[ TNFSF10\ aşağı\ akım\ varyantları \rightarrow Alternatif\ RNA\ kesilip\ birleştirilmesi \rightarrow Farklı\ TNFSF10\ transkriptleri \]

Bu iki göstərim tədqiqatın inteqrativ hipotezini izah edir. Tədqiqatda zəncirin bütün addımları eksperimental şəkildə doğrulanmamışdır.

Tədqiqatın güclü tərəfləri

• GWAS siqnalları tək anotasiya mənbəyi ilə şərh edilməmiş, çoxsaylı müstəqil məlumat qatı birləşdirilmişdir.
• 5.184 hal və 311.848 nəzarəti əhatə edən böyük Tayvan GWAS məlumat mənbəyindən istifadə edilmişdir.
• İnterferonla əlaqəli genlər yalnız bioloji ilkin fərziyyə ilə deyil, funksional zənginləşdirmə analizləri ilə prioritetləşdirilmişdir.
• Gen ifadəsi ilə RNT kəsilib-birləşdirilməsi bir-birindən ayrılmışdır.
• GTEx məlumatları çoxsaylı toxuma üzrə müqayisə edilmişdir.
• Tayvan və Şərqi Asiya əlaqə quruluşu digər populyasiyalarla müqayisə edilmişdir.
• ENCODE, ORegAnno, ReMap, H3K27ac və Micro-C kimi müxtəlif tənzimləyici məlumat mənbələri birlikdə qiymətləndirilmişdir.
• Müəlliflər motiv, QTL və xromatin nəticələrini səbəbsəl sübut kimi təqdim etməmiş, eksperimental doğrulama ehtiyacını açıq şəkildə qeyd etmişdir.

Tədqiqatın məhdudiyyətləri

GWAS xülasə məlumatları: Analiz əsasən xülasə statistikalarına əsaslanır. Eyni şəxslərdə genotip, RNT izoformu, zülal səviyyəsi və klinik aktivlik birlikdə ölçülməmişdir.

Səbəbsəl variant müəyyən edilmədi: Güclü LD bloklarında çoxsaylı SNP oxşar siqnal göstərə bilər. Tədqiqat hansı variantın həqiqi funksional dəyişikliyi yaratdığını dəqiqləşdirə bilmir.

Motivlər yalnız proqnozdur: IRF5, IRF8 və STAT4 motivləri birbaşa zülal-DNT bağlanması sübutu deyil.

GTEx əsasən sağlam toxumalara əsaslanır: GTEx QTL məlumatları lupus xəstəliyi zamanı aktivləşmiş immun hüceyrələri tam təmsil etməyə bilər.

Hüceyrə tipi rezolyusiyası məhduddur: Tam qan çoxsaylı hüceyrə növünün qarışığıdır. Müşahidə olunan sQTL siqnalının B hüceyrəsindən, T hüceyrəsindən, monositdən və ya başqa alt qrupdan gəldiyi müəyyən edilə bilmir.

EBV ilə transformasiya edilmiş hüceyrələr təbii limfosit deyil: EBV transformasiyası hüceyrə gen ifadəsini və RNT işləmə proqramlarını dəyişə bilər.

Micro-C hüceyrə konteksti fərqlidir: Xromatin təması məlumatı H1 embrionik kök hüceyrələrindən alınmışdır; lupus immun hüceyrələrindən deyil.

Müstəqil Tayvan replikasiyası yoxdur: Populyasiyaya xas variant və LD nəticələri müstəqil Tayvan lupus kohortunda yenidən yoxlanılmamışdır.

İzoformlar birbaşa ölçülməyib: Tədqiqat sQTL siqnallarını bildirir, lakin hansı TNFSF10 izoformunun artdığını və ya azaldığını funksional səviyyədə göstərmir.

Zülal və hüceyrə nəticəsi ölçülməyib: TRAIL zülal səviyyəsi, reseptor siqnalı, T hüceyrə aktivləşməsi və ya lupusla əlaqəli hüceyrə fenotipləri variantlara görə eksperimental şəkildə müqayisə edilməmişdir.

Klinik əlaqə qurulmayıb: Variantlar lupus aktivliyi, böyrək tutulması, müalicə cavabı və ya proqnozla birbaşa sınaqdan keçirilməmişdir.

Tədqiqat nə deyir, nə demir?

Keçmiş, bu gün və gələcək baxımından əhəmiyyəti

Keçmiş baxımından: Ənənəvi GWAS şərhləri çox vaxt ən yaxın geni xəstəlik geni kimi qəbul etməyə meylli idi. Bu tədqiqat variantların fiziki olaraq yaxın olduqları gendən başqa uzaq tənzimləyici bölgələr, RNT işlənməsi və üçölçülü xromatin təmasları vasitəsilə təsir göstərə biləcəyini vurğulayır.

Bu gün baxımından: Lupus genetikasında yalnız gen ifadəsinə fokuslanmağın yetərli olmaya biləcəyini göstərir. TNFSF10 bölgəsində dominant sQTL siqnalları alternativ RNT izoformlarının lupus araşdırmalarında daha çox öyrənilməli olduğunu düşündürür.

Gələcək baxımından: Nəticələr aşağıdakı təcrübələrin dizaynına kömək edə bilər:

• TNFSF10 izoformlarının uzun oxunuşlu RNT sekvensləşdirilməsi ilə müəyyən edilməsi,
• Lupus xəstələrinin təmizlənmiş B hüceyrəsi, T hüceyrəsi və monositlərində sQTL analizi,
• rs62282578 və digər namizəd variantların CRISPR ilə dəyişdirilməsi,
• IRF5, IRF8 və STAT4 bağlanmasının ChIP üsulları ilə sınaqdan keçirilməsi,
• Allel-spesifik enhancer reportyor təcrübələri,
• Lupus immun hüceyrələrində xromatin təması xəritələnməsi,
• Fərqli TNFSF10 izoformlarının TRAIL istehsalı və immun hüceyrə funksiyasına təsirinin ölçülməsi,
• Nəticələrin müstəqil Tayvan və digər Şərqi Asiya lupus kohortlarında doğrulanması.

Gündəlik həyat və klinik tətbiq baxımından mənası

Bu araşdırma birbaşa xəstələrə tətbiq edilə biləcək test və ya müalicə təqdim etmir. Sistemik lupus diaqnozu mövcud klinik qiymətləndirmələr, autoanticisimlər, orqan tutulmaları və mütəxəssis həkim müayinəsi əsasında qoyulur. Tədqiqatda önə çıxarılan SNP-lər hazırda diaqnoz qoymaq, xəstəlik şiddətini proqnozlaşdırmaq və ya dərman seçmək üçün təsdiqlənmiş deyil.

Tədqiqatın gündəlik həyat baxımından dolayı əhəmiyyəti eyni xəstəlik diaqnozuna sahib insanların genetik və molekulyar mexanizmlərinin fərqli ola biləcəyini göstərməsidir. Gələcəkdə interferon aktivliyi, TNFSF10 izoformları və populyasiyaya xas variantlar birlikdə doğrulanarsa, lupus xəstələrinin molekulyar alt qruplara ayrılmasına töhfə verə bilər. Lakin mövcud tədqiqat bu fərdiləşdirilmiş tibb tətbiqinin klinik uğur sübutunu təqdim etmir.

Tədqiqatın Metodu və Nəticələri

Texniki tədqiqat dizaynı

Analiz iş axını

1. TPMI lupus GWAS xülasə statistikalarından genom miqyasında əhəmiyyətli SNP-lər seçildi.
2. Variantların yaxınlığındakı və ya uyğunlaşdırılan namizəd genlər müəyyən edildi.
3. DAVID ilə Gene Ontology və Reactome yol zənginləşdirməsi aparıldı.
4. Oxşar funksional terminlər Kappa statistikasına əsaslanan təsnifatla qruplaşdırıldı.
5. Enrichr-KG ilə gen-yol funksional şəbəkəsi yaradıldı.
6. TNFSF10-un yuxarı axın və 5′ intronik bölgələrində transkripsiya faktoru motivləri skan edildi.
7. IRF5, IRF8 və STAT4 aparıcı variantlarının allel tezlikləri populyasiyalar arasında müqayisə edildi.
8. LocusZoom və PLINK istifadə edilərək linkage disequilibrium quruluşları qiymətləndirildi.
9. GTEx vasitəsilə TNFSF10 ətrafındakı eQTL və sQTL siqnalları araşdırıldı.
10. UCSC Genome Browser vasitəsilə ENCODE cCRE, ORegAnno, ReMap və H3K27ac anotasiya qatları əlavə edildi.
11. Micro-C məlumatları ilə uzun məsafəli xromatin təmasları qiymətləndirildi.
12. Bütün sübut qatları TNFSF10 üçün namizəd tənzimləyici və splicing variantlarını prioritetləşdirmək məqsədilə birləşdirildi.

İstifadə olunan məlumat mənbələri və alətlər

Əsas funksional zənginləşdirmə nəticələri

IRF5, IRF8 və STAT4 nəticələrinin texniki xülasəsi

TNFSF10 QTL tapıntılarının texniki xülasəsi

Tənzimləyici və xromatin tapıntılarının xülasəsi

Araşdırmanın yekun texniki şərhi

Tədqiqatda bir-birini tamamlayan üç əsas sübut qatı yaradıldı:

1. GWAS və yol qatı: Tayvan lupus GWAS namizədləri interferon və sitokin siqnallaşmasında qruplaşdı.
2. Transkripsional tənzimləmə qatı: IRF5, IRF8 və STAT4 üçün TNFSF10 ətrafında proqnozlaşdırılan motivlər tapıldı.
3. RNT işlənməsi və xromatin qatı: TNFSF10-un aşağı axınındakı variantlar sQTL siqnalları, tənzimləyici anotasiya, populyasiyaya xas LD və mümkün uzun məsafəli xromatin təmasları ilə üst-üstə düşdü.

Bu qatların birləşməsi TNFSF10 bölgəsindəki bəzi variantların lupus biologiyası baxımından daha ətraflı təcrübələrə götürülməsini dəstəkləyir. Bununla belə, tədqiqat heç bir variantı səbəbsəl, diaqnostik, proqnostik və ya müalicə müəyyən edən kimi doğrulamamışdır.

Mənbə və Metod Qeydi

Bu məzmun Kuo-Sheng Hung, Shan-Wen Lui, Jia-Huei Lyu, Hsiang-Cheng Chen, Ya-Ting Yang, Yu-Tien Chang, Feng-Cheng Liu, Tzu-Sheng Hsu və I-Tsu Chyuan tərəfindən hazırlanmış “Integrative post-GWAS analysis prioritizes interferon-linked regulatory and splicing-associated variants at the TNFSF10 locus in systemic lupus erythematosus” başlıqlı tədqiqata əsaslanır.

Tədqiqat SSRN-də yerləşən preprint tədqiqat məqaləsidir. Mətndə açıq şəkildə “This preprint research paper has not been peer reviewed” ifadəsi var. Buna görə nəticələr rəyçi qiymətləndirməsindən keçməmişdir. Mətndə jurnal qəbulu və ya məqaləyə aid DOI bildirilməmişdir.

Araşdırma klinik tədqiqat, dərman sınağı, hüceyrə kulturası təcrübəsi və ya heyvan təcrübəsi deyil. Taiwan Precision Medicine Initiative lupus GWAS xülasə statistikalarının ictimai və institusional funksional genomika mənbələri ilə birləşdirildiyi hesablama yönümlü post-GWAS analizidir.

Tədqiqat yeni xəstə nümunələrində TNFSF10 RNT izoformu, TRAIL zülalı, transkripsiya faktoru bağlanması və ya immun hüceyrə funksiyasını ölçməmişdir. IRF5, IRF8 və STAT4 motivləri hesablama proqnozudur. GTEx eQTL və sQTL nəticələri toxuma səviyyəsindəki əlaqələrdir. Micro-C məlumatı lupus xəstələrinin immun hüceyrələrindən deyil, H1 insan embrionik kök hüceyrələrindən əldə edilmişdir.

Nəticələr variantların funksional və səbəbsəl olduğunu sübut etmir. Tədqiqat yalnız gələcək ChIP, reportyor gen, CRISPR, allel-spesifik ifadə, xromatin profilləmə və transkript izoform təcrübələri üçün namizəd bölgələri prioritetləşdirir.

Tədqiqatda önə çıxarılan rs3757387, rs35703946, rs7568275, rs62282578 və ya digər TNFSF10 variantları hazırda klinik lupus diaqnozu, xəstəlik aktivliyi, proqnoz və ya müalicə seçimi üçün doğrulanmış biomarkerlər deyil. Araşdırma yeni dərman, müalicə zəmanəti və ya klinik qərar vasitəsi təqdim etmir.

TPMI tədqiqatının müvafiq etik şuralar tərəfindən təsdiqləndiyi və iştirakçıların genomik analizlə elektron sağlamlıq qeydlərinin əlaqələndirilməsinə razılıq verdiyi bildirilmişdir. Mövcud analiz GWAS xülasə statistikalarından istifadə etmişdir. Fərdi TPMI və Taiwan Biobank məlumatları institusional məlumat idarəçiliyi və məxfilik məhdudiyyətləri səbəbindən ictimai deyil.

Müəlliflər maraqlar toqquşması bildirməmişdir. Tədqiqat Tayvandakı Ministry of Science and Technology və National Science and Technology Council tərəfindən göstərilən layihələrlə dəstəklənmişdir; maliyyələşdirən qurumların dizayn, analiz, şərh və ya nəşr qərarında rol oynamadığı bəyan edilmişdir.

Müəlliflərin bəyanına görə ChatGPT yalnız dil redaktəsi və məqalə formatlaşdırmasına yardım məqsədilə istifadə edilmiş; məzmun sonradan müəlliflər tərəfindən yoxlanılmış və redaktə edilmişdir. Müəlliflər tədqiqatın bütün elmi məzmununun məsuliyyətini daşıdıqlarını bildirmişlər.

Bu məqalə yalnız yüklənmiş tədqiqatda təqdim edilən mətn, metod, cədvəl, şəkil və müəllif şərhlərinə əsaslanaraq hazırlanmışdır. PDF-də olmayan səbəbsəl, klinik və ya müalicəyə yönəlik iddialar əlavə edilməmişdir. eQTL, sQTL və xromatin tənzimlənməsini izah etmək üçün istifadə olunan ox göstərimləri tədqiqatın konseptual məntiqini sadələşdirmək məqsədilə verilmiş, məqalədə yayımlanmış riyazi tənliklər olmadığı açıq şəkildə qeyd edilmişdir.