D.YU. KLYUCHNIKOV ET AL., 2022
analysis of associations of genetic markers with the development of congenital scoliosis
D.Yu. Klyuchnikov1, E.Yu. Filatov2, I.V. Tyumin3, O.V. Tyumina1'3
1Samara Regional Medical Centre "Dynasty", Samara, Russia, 2National Ilizarov Medical Research Centre for Traumatology and Orthopaedics, Kurgan, Russia 3Samara State Medical University, Samara, Russia
Objective. To study the associations of single-nucleotide polymorphisms: rs6570507 in GPR126 gene, rs1800795 in IL-6 gene, rs1800469 in TGFB1 gene, rs731236 in VDR gene, rs625039 and rs11598564 polymorphisms in LBX1 gene, and rs12946942 in SOX9 gene with congenital scoliosis.
Material and Methods. The study included 90 patients with verified congenital anomalies of the spine (single and multiple malformations of the spine, ICD-10 Code: Q76.3) and 157 clinically healthy volunteers without diagnosed spinal deformity and without family history of spinal malformations or osteoarticular system diseases. Molecular genetic testing was performed by PCR with real-time registration of a signal from the developed oligonucleotides used to determine rs6570507, rs1800795, rs1800469, rs625039, rs11598564, rs12946942, and rs731236 polymorphisms. Reference sequences were selected from the dbSNP database, and sequence design was performed on the BLAST platform. Data analysis was performed using the R free software computing environment. Data were compared using Pearson's X2 test, and 95 % confidence interval limits were calculated to assess the significance of OR.
Results. Statistically significant association of the G allele and GG genotype of the rs1800795 polymorphism in the interleukin-6 gene with congenital scoliosis was found in group of Russian patients (p < 0.001). No significant association of alleles and genotypes of polymorphic variants of rs6570507, rs1800469, rs625039, rs11598564, rs12946942, and rs731236 with congenital scoliosis was found. Conclusion. The rs1800795 polymorphism can be considered as a promising marker for molecular genetic diagnostics of congenital scoliosis. Key Words: congenital scoliosis, single-nucleotide polymorphisms, SNPs, association of genetic markers.
Please cite this paper as: Klyuchnikov DYu, Filatov EYu, Tyumin IV, Tyumina OV. Analysis of associations of genetic markers with the development of congenital scoliosis. Hir. Pozvonoc. 2022;19(2):33—39. In Russian. DOI: http://dx.doi.org/10.14531/ss2022.2.33-39.
Congenital scoliosis is frontal or sagittal spine curvature, caused by abnormal vertebral segmentation and development [1].
Congenital scoliosis accounts for about 10% of spinal deformities [2]. The incidence of congenital spine anomalies ranges from 0.5 to 1.0 per 1,000 of liveborns [3]. The cause of congenital vertebral malformations may be disorders that occur at the stage of somitogenesis. [4].
Congenital spinal deformities are progressive in 50 % of cases, and patients with this disease need constant monitoring.
Currently, congenital anomalies of spine are detected with ultrasound screening. However, the possibilities of ultrasound in prenatal diagnostics of such anomalies are limited. Spine curvature may not be obvious at birth and diagnostics can fail. Therefore, developmental anomalies are often detected in the postnatal period much later, when
an extended primary scoliotic curve, quite often with structural compensatory curves, have already been formed [5].
In addition to ultrasound diagnostics, congenital scoliosis can be detected using X-ray methods, but they create undesirable radiation exposure on patients.
In this regard, studying the genetic nature of the onset of congenital scolio-sis is of great interest in order to develop early molecular genetic screening tests identifying the risks of the disease onset and progression. It is known that abnormal vertebral segmentation can be associated with more than 150 genetic disorders [6]. Data given by different authors, however, are contradictory and have not been reproduced in other studies.
The objective is to study the association of polymorphisms array with the risk of the development of congenital scoliosis in order to assess the possibilities of early molecular genetic testing of this pathology.
33
Despite the fact that congenital and idiopathic scolioses are clinically different, their pathogenetic mechanisms are similar [7]. Therefore, SNPs with the described association with the development of congenital and idiopathic scoliosis in different groups of patients were selected for analysis: polymorphism rs6570507 in GPR126 gene [8], rs1800795 in IL-6 gene [9], rs1800469 in TGFB1 gene [10], rs731236 in VDR gene [11], polymorphisms rs625039 and rs11598564 in LBX1 gene [12], as well as rs12946942 in SOX9 gene [13].
Materials and Methods
The study included 90 patients aged from 10 months up to 16 years, who underwent treatment at the National Ilizarov Medical Research Centre for Traumatology and Orthopaedics in 2019-2020, with verified congenital spine anomalies (single and multiple malformations of the spine, ICD-10 code: Q76.3 congenital
scoliosis caused by bone malformation). Also, the study included 157 clinically healthy volunteers aged 19-49 years without a diagnosed spinal deformity who did not have family history of spinal developmental pathology and diseases of the osteoarticular system. All patients or their legal representatives as well as volunteers from the comparison group, gave written informed consent to collecting biomaterial and conducting research. This study was performed on the basis of the immunological typing laboratory of the Dynasty Medical Center, the study protocol was approved by the Ethics Committee of the Centre as of January 11, 2019. The research was carried out in accordance with the WMA Declaration of Helsinki - Ethical Principles for Medical Research Involving Human Subjects with amendments as of 2000.
Venous blood samples and buccal epithelium scrapings were collected as a DNA source. Real-time PCR was performed on CFX96, BioRad and iQ5, Biorad (USA) amplifiers using the designed oligonucleotide sequences to determine rs6570507, rs1800795, rs1800469, rs625039, rs11598564, rs12946942, and rs731236. Reference sequences were selected from dbSNP database [14], sequence design was performed on Basic Local Alignment Search Tool (BLAST) platform [15].
The data analysis was performed with R free software computing environment (v.3.5.1) [16], the nominal data were compared using Pearson's 2 test, and 95 % confidence interval limits were calculated to assess the significance of OR (95 % CI).
Results
Allele and genotype frequency distributions in the study (n = 90) and the control (n = 157) groups by the analyzed polymorphisms rs11598564, rs12946942, rs1800469, rs625039, rs1800795, rs6570507, and rs731236 are given in Tables 1, 2.
rs11598564. Comparison of the allele carriership of the rs1 1598564 polymorphism in the study and control
groups showed that the odds ratio (OR) for carriership of the G allele against the A allele is equal to 1.118 (95 CI 0.7691.624). An analysis of a recessive model GG against AA + AG gives OR for the GG genotype as a risk factor at 1.235 (95 CI 0.688-2.215).
rs12946942. When comparing the study and the control groups, OR of the T allele as a risk factor is equal to 1.086 (95 CI 0.615- 1.917), and the TT genotype against GG + GT - 1.144 (95 CI 0.188-6.979).
rs1800469. OR of the G allele against an alternative A allele in the sample analysis is equal to 1.152 (95 CI 0.7701.724), and the homozygote genotype GG against AA + AG - 1.260 (95 CI 0.749-2.120).
rs1800795. Analyzing this polymorphism, it is found out that OR of the G risk allele against C risk allele is equal to 2.152 (95 CI 1.456-3.182). At the same time, 2 = 15.011, higher than the critical value, the significance of 2 < 0.001. OR of the genotype GG against CC + CG is equal to 2.764 (95 CI 1.602-4.770), at the same time the confidence interval limits, 2 = 13.686 and the significance level of 2 < 0.001 indicate the confidence of the GG genotype association.
rs625039. OR of the G allele as a risk factor against the A allele is equal to 1.368 (95 CI 0.76-2.464), and the GG genotype against AA + AG is equal to 1.344 (95 CI 0.704-2.564).
rs6570507. An analysis of this polymorphism shows that the odds of the G allele against an alternative A allele is 1.324 (95 CI 0.860-2.038), the GG genotype against AG + AA - 1.247 (95 CI 0.729-2.132).
rs731236. An analysis demonstrates that the odds of the A allele variant against G is 1.247 (95 CI 0.835-1.861), the AA homozygote genotype against AG + GG - 1.201 (95 CI 0.713-2.024). The results of OR of alleles and genotypes are given in Tables 3, 4.
The determined alleles and genotypes in the analyzed polymorphic loci were more frequent in the group of patient with congenital spine pathology: the G allele (OR 1.324) and the GG genotype (OR 1.247) for rs6570507 in GPR126
gene; the G allele (OR 1.152) and the GG genotype (OR1.26) for rs1800469 in TGFB1 gene; the A allele (OR 1.247) and the AA genotype (OR 1.201) for rs731236 in VDR gene; the G allele (OR 1.368) and the GG genotype (OR1.344) for rs625039 in LBX1 gene; the G allele (OR 1.118) and the GG genotype (OR 1.235) for rs11598564 of LBX1 gene; the T allele (OR 1.086) and the TT genotype (OR 1.144) for rs12946942 in SOX9 gene. However, confident association was found only for the G allele and the GG genotype for rs1800795 in IL-6 gene. Thus, it is shown that carriership of the GG genotype can increase the risk of congenital scoliosis by 2.7 times.
Discussion
Among the studied polymorphisms, a confident association was found for the GG genotype of the rs1800795 polymorphism. This polymorphism was identified on the short arm of the 7th chromosome (chr7:22727026 (GRCh38.p13) in the promoter region of IL-6 gene. The presence of the G allele in the genotype is associated with the enhanced expression of IL-6 gene [17, 18]. It is known that vertebral bodies are formed from somites as a result of budding from presomitic mesoderm, mediated by expression of the FGF, Wnt genes and genes of the Notch signaling pathway [19], while downregulation of Notch signaling pathway is observed at an elevated level of IL-6 [20]. The evidence suggests that a signal pathway disorder, which decreases Notch pathway signaling and regulatory failure of this process, are quite important in the pathogenesis of congenital spine deformities [21].
In view of the association with sco-liosis, rs1800795 has been analyzed in several studies with contradictory results. It is detected that carriers of the GG genotype have a higher risk of developing adolescent idiopathic scoliosis among the representatives of the European race [9] which was not detected in the Asian population [22]. The literature also points out that such polymorphism is reliably associated with the severity of the disease [23]. The data obtained as a result of our
34
Table 1
Allele frequencies
Polymorphism Allele Study group Control group
(n = 90) (n = 157)
frequency number frequency number
rs11598564 A 0.46 79 0.49 150
G 0.54 93 0.51 158
rs12946942 G 0.88 158 0.89 273
T 0.12 22 0.11 35
rs1800469 A 0.28 51 0.31 97
G 0.72 129 0.69 213
rs1800795 C 0.31 55 0.49 143
G 0.69 125 0.51 151
rs625039 A 0.10 18 0.13 41
G 0.90 158 0.87 263
rs6570507 A 0.23 42 0.29 81
G 0.77 138 0.71 201
rs731236 A 0.72 129 0.67 209
G 0.28 51 0.33 103
Table 2
Genotype frequencies
Polymorphism Genotype Study group Control group
(n = 90) (n = 157)
frequency number frequency number
rs11598564 AA 0.221 19 0.234 36
AG 0.477 41 0.506 78
GG 0.302 26 0.260 40
rs12946942 GG 0.778 70 0.792 122
GT 0.200 18 0.188 29
TT 0.022 2 0.020 3
rs1800469 AA 0.089 8 0.090 14
AG 0.389 35 0.445 69
GG 0.522 47 0.465 72
rs1800795 CC 0.144 13 0.265 39
CG 0.322 29 0.442 65
GG 0.533 48 0.292 43
rs625039 AA 0.011 1 0.026 4
AG 0.182 16 0.217 33
GG 0.807 71 0.757 115
rs6570507 AA 0.067 6 0.121 17
AG 0.333 30 0.333 47
GG 0.600 54 0.546 77
rs731236 AA 0.489 44 0.449 70
AG 0.456 41 0.422 69
GG 0.056 5 0.109 17
research confirm the described studies in a group of the Russian patients and allow us to attribute the GG genotype of the
rs1800795 polymorphism to the risk factors for developing congenital scoliosis.
35
The results on the frequencies of alleles and genotypes of rs1800469 differ from the results obtained on a sample of the patients from Moscow region, where the OR for the G allele is 1.73 and the GG genotype is 4.82 in comparison with the results obtained in the experiment -1.152 and 1.260, respectively [10]. The results on the frequencies of the GG genotype of the rs731236 polymorphism in VDR gene differ from those described in the other group of Russian patients [11].
Some differences were also noted in the frequencies of the rs1800795 and rs731236 alleles in the control group and the frequencies given in the databases of the ALFA, gnomAD-Genomes, 1000 Genomes projects and described for Cau-casoids. Allele frequencies in the control group are presented in Table 5.
On the one hand, a statistically unconfirmed association may be due to the small sample of the study and relative heterogeneity of the patients. On the other hand, the data of some studies are not reproduced in other populations due to the difference of ethnic origin that is confirmed by experimental data on the difference in the allele frequencies of the rs1800795 and rs731236 polymorphisms between the control group (inhabitants of the Middle Volga region) and the data in the European populations. The difference in the frequencies of alleles and genotypes in populations speaks in favor of forming relevant comparison groups based on population features.
Congenital scoliosis can be a polygen-ic disease, and polymorphisms can have different severity or affect only in combination. Introduction of risks calculation for development of congenital or idiopathic scoliosis using genetic markers requires additional research, particularly, of intergenic interactions and non-genetic environmental factors. The risk of congenital spinal disorders is also associated with pregnancy course, diabetes and maternal smoking [27, 28].
Conclusion
A frequency association of the G allele and the GG genotype of the rs1800795 polymorphism with congenital scoliosis was revealed on the example of a group
Table 4 Genotype odds ratio
Polymorphism Model OR (95 CI) x2 P
rs11598564 GG AA + AG 1.235 (0.688-2.215) 0.502 0.479
rs12946942 TT GT + GG 1.144 (0.188-6.979) 0.021 0.885
rs1800469 GG AG + AA 1.260 (0.749-2.120) 0.759 0.384
rs1800795 GG CG + CC 2.764 (1.602-4.770) 13.686 0.0002
rs625039 GG AG + AA 1.344 (0.704-2.564) 0.807 0.370
rs6570507 GG AG + AA 1.247 (0.729-2.132) 0.650 0.421
rs731236 AA AG + GG 1.201 (0.713-2.024) 0.475 0.491
Table 3 Allele odds ratio
Polymorphism Allele OR (95 CI) x2 P
rs11598564 A 0.895 (0.616-1.301) 0.340 0.560
G 1.118 (0.769-1.624)
rs12946942 G 0.921 (0.522-1.625) 0.081 0.776
T 1.086 (0.615-1.917)
rs1800469 A 0.868 (0.580-1.300) 0.472 0.492
G 1.152 (0.770-1.724)
rs1800795 C 0.465 (0.314-0.687) 15.011 0.0001
G 2.152 (1.456-3.182)
rs625039 A 0.731 (0.406-1.316) 1.099 0.295
G 1.368 (0.760-2.464)
rs6570507 A 0.755 (0.491-1.162) 1.634 0.202
G 1.324 (0.860-2.038)
rs731236 A 1.247 (0.835-1.861) 1.162 0.281
G 0.802 (0.537-1.198)
of the Russian patients (n = 90). The obtained results show that rs1800795 can be considered as one of the promising candidates for diagnostic markers of congenital scoliosis development, and require validation.
Molecular genetic screening tests in combination with other methods for early diagnostics of congenital scoliosis development will allow early risks diagnosis and the patients stratification. It will give an opportunity to assign timely and effective treatment and will enhance the patients' quality of life.
Phenotypic heterogeneity within the groups with congenital spinal deformities and a limited number of studies on small samples remain the main difficulties in determining the genetic basis for the development of congenital and idio-pathic scoliosis. Therefore, it is important and expedient to develop not only a system for registering patients with congenital spinal deformities, but also biobanking the genetic material of patients for testing and validating research results.
The authors are grateful to Natalya Sergeevna Osinovskaya, PhD in Biology, Senior Researcher at the Genomics Laboratory with the Bioresource Collection Group of the Scientific and Research Institute of Obstetrics, Gynecology and Reproductology named after D. O. Ott, for the support in preparing the manuscript for publication.
The study was carried out with the financial support of the Fund for the Promotion of Small Forms of Enterprises in the Scientific and Technical Field within the framework of the 3741TC1/63150 project.
The authors declare that they have no conflict of interest.
36
Table 5 Occurrence frequency of polymorphism alleles rs6570507, rs1800795, rs1800469, rs625039, rs11598564, rs12946942, rs731236 [24-26]
Polymorphism ALFA: European gnomAD-Genomes, European 1000 Genome, European Experimental data (comparison group)
rs6570507 G = 0.71 G = 0.71 G = 0.68 G = 0.71
A = 0.29 A = 0.29 A = 0.32 A = 0.29
rs1800795 C = 0.44 C = 0.46 C = 0.42 C = 0.49
G = 0.56 G = 0.54 G = 0.58 G = 0.51
rs1800469 A = 0.32 A = 0.30 A = 0.31 A = 0.31
G = 0.68 G = 0.70 G = 0.69 G = 0.69
rs625039 G = 0.86 G = 0.88 G = 0.86 G = 0.87
A = 0.14 A = 0.12 A = 0.14 A = 0.13
rs11598564 G = 0.52 G = 0.53 G = 0.54 G = 0.51
A = 0.48 A = 0.47 A = 0.46 A = 0.49
rs12946942 G = 0.93 G = 0.92 G = 0.91 G = 0.89
T = 0.07 T = 0.08 T = 0.09 T = 0.11
rs731236 A = 0.60 A = 0.63 A = 0.60 A = 0.67
G = 0.40 G = 0.37 G = 0.40 G = 0.33
37
spine deformities
References
1. Erol B, Kusumi K, Lou J, Dormans JP. Etiology of congenital scoliosis. Univ Penn Orthop J. 2002;15:37-42.
2. Birnbaum K, Weber M, Lorani A, Leiser-Neef U, Niethard FU. Prognostic significance of the Nasca classification for the long-term course of congenital scoliosis. Arch Orthop Trauma Surg. 2002;122:383-389. DOI: 10.1007/s00402-002-0401-z.
3. Wynne-Davies R. Congenital vertebral anomalies: aetiology and relationship to spina bifida cystica. J Med Genet. 1975;12:280-288. DOI: 10.1136/jmg.12.3.280.
4. Kusumi K, Turnpenny PD. Formation errors of the vertebral column. J Bone Joint Surg. 2007; 89(Suppl 1):64-71. DOI: 10.2106/JBJS.F.00486.
5. Ryabykh SO, Savin DM, Filatov Elu. A clinical case of surgical treatment of severe congenital kyphoscoliosis in an 11-year old child. Genij Ortopedii. 2017;23(2):216-219. DOI: 10.18019/1028-4427-2017-23-2-216-219.
6. Shifley ET, Cole SE. The vertebrate segmentation clock and its role in skeletal birth defects. Birth Defects Res. C Embryo Today. 2007;81:121-133. DOI: 10.1002/ bdrc.20090.
7. Giampietro PF. Genetic aspects of congenital and idiopathic scoliosis. Scientifica (Cairo). 2012;2012:152365. DOI: 10.6064/2012/152365.
8. Kou I, Watanabe K, Takahashi Y, Momozawa Y, Khanshour A, Grauers A, Zhou H, Liu G, Fan YH, Takeda K, Ogura Y, Zhou T, Iwasaki Y, Kubo M, Wu Z, Matsumoto M, Einarsdottir E, Kere J, Huang D, Qiu G, Qiu Y, Wise CA, Song YQ, Wu N, Su P, Gerdhem P, Ikegawa S. A multi-ethnic meta-analysis confirms the association of rs6570507 with adolescent idiopathic scoliosis. Sci Rep. 2018;8:11575. DOI: 10.1038/s41598-018-29011-7.
9. Sobhan MR, Mahdinezhad-Yazdi M, Dastgheib SA, Ahrar H, Aghili K, Neamat-zadeh H. Association of the IL-6-174G > C (rs1800795) polymorphism with adolescent idiopathic scoliosis: evidence from a case-control study and meta-analysis. Rev Bras Ortop. 2020;55:17-26. DOI: 10.1055/s-0039-1700813.
10. Ryzhkov II, Borzilov EE, Churnosov MI, Ataman AV, Dedkov AA, Poloni-kov AV. Transforming growth factor beta 1 is a novel susceptibility gene for adolescent idiopathic scoliosis. Spine. 2013;12:E699-704. DOI: 10.1097/BRS.0b013e31828de9e1.
11. Vissarionov SV, Larionova VI, Kazarian IV, Filippova AN, Kostik MM, Voitov-ich AN, Rotchev EV. The gene polymorphisms of COL1A1 and VDR in children with scoliosis. Pediatric Traumatology, Orthopaedics and Reconstructive Surgery. 2017;5(1):5-12. DOI: 10.17816/PTORS515-12.
12. Jiang H Yang Q, Liu Y, Guan Y, Zhan X, Xiao Z, Wei Q. Association between ladybird homeobox 1 gene polymorphisms and adolescent idiopathic scoliosis: A MOOSE-compliant meta-analysis. Medicine (Baltimore). 2019; 98:e16314. DOI: 10.1097/ MD.0000000000016314.
13. Takeda K, Kou I, Otomo N, Grauers A, Fan YH, Ogura Y, Takahashi Y, Momozawa Y, Einarsdottir E, Kere J, Matsumoto M, Qiu Y, Song Y-Q, Gerd-hem P, Watanabe K, Ikegawa S. A multiethnic meta-analysis defined the association of rs12946942 with severe adolescent idiopathic scoliosis. J Hum Genet. 2019;64:493-498. DOI: 10.1038/s10038-019-0575-7.
14. Sherry ST, Ward M, Sirotkin K. DbSNP - database for single nucleotide polymorphisms and other classes of minor genetic variation. Genome Res. 1999;9:677-679.
15. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-410. DOI: 10.1016/S0022-2836(05)80360-2.
16. Team R Core. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, 2019. Available at https://www.R-project. org/.
17. Guan Y, Wang S, Wang J, Meng D, Wu H, Wei Q, Jiang H. Gene polymorphisms and expression levels of interleukin-6 and interleukin-10 in lumbar disc disease: a meta-
analysis and immunohistochemical study. J Orthop Surg Res. 2020;15:54. DOI: 10.1186/ s13018-020-01588-8.
18. Gonzalez-Castro TB, Hernandez-Diaz Y, Perez-Hernandez N, Tovilla-Zarate CA, Juarez-Rojop IE, Lopez-Narvaez ML, Blachman-Braun R, Posadas-Sanchez R, Vargas-Alarcon G, Garcia-Flores E, Cazarin-Santos BG, Borgonio-Cuadra VM, Reyes-Lopez PA, Rodriguez-Perez JM. Interleukin 6 (rs1800795) gene polymorphism is associated with cardiovascular diseases: a metaanalysis of 74 studies with 86,229 subjects. EXCLI J. 2019;18:331-355. DOI: 10.17179/ excli2019-1248.
19. Pourquie O. Vertebrate segmentation: from cyclic gene networks to scoliosis. Cell. 2011;145:650-663. DOI: 10.1016/j.cell.2011.05.011.
20. Milinkovic I, Djinic Krasavcevic A, Nikolic N, Aleksic Z, Carkic J, Jezdic M, Jankovic S, Milasin J. Notch down-regulation and inflammatory cytokines and RANKL overexpression involvement in peri-implant mucositis and peri-implantitis: A cross-sectional study. Clin Oral Implants Res. 2021;32:1496-1505. DOI: 10.1111/ clr.13850.
21. Barhoumi T Nashabat M, Alghanem B, Alhallaj A, Boudjelal M, Umair M, Alarifi S, Alfares A, Mohrij SAA, Alfadhel M. Delta Like-1 gene mutation: a novel cause of congenital vertebral malformation. Front Genet. 2019;10:534. DOI: 10.3389/ fgene.2019.00534.
22. Liu Z, Tang NL, Cao XB, Liu WJ, Qiu XS, Cheng JC, Qiu Y. Lack of association between the promoter polymorphisms of MMP-3 and IL-6 genes and adolescent idio-pathic scoliosis: a case-control study in a Chinese Han population. Spine. 2010;35: 1701-1705. DOI: 10.1097/BRS.0b013e3181c6ba13.
23. Nikolova ST, Yablanski VT, Vlaev EN, Stokov LD, Savov AS, Kremensky IM, Loukanov AR. Association between IL-6 and MMP3 common genetic polymorphisms and idiopathic scoliosis in Bulgarian patients: a case-control study. Spine. 2016;41: 785-791. DOI: 10.1097/BRS.0000000000001360.
24. Phan L, Jin Y, Zhang H, Qiang W, Shekhtman E, Shao D, Revoe D, Villamarin R, Ivanchenko E, Kimura M, Wang ZY, Hao L, Sharopova N, Bihan M, Sturcke A, Lee M, Popova N, Wu W, Bastiani C, Ward M, Holmes JB, Lyo-shin V, Kaur K, Moyer E, Feolo M, Kattman BL. ALFA: Allele Frequency Aggregator. National Center for Biotechnology Information, U.S. National Library of Medicine. 2020. Available at: www.ncbi.nlm.nih.gov/snp/docs/gsr/alfa/.
25. Fairley S, Lowy-Gallego E, Perry E, Flicek P. The International Genome Sample Resource (IGSR) collection of open human genomic variation resources. Nucleic Acids Res. 2020;48:D941-D947. DOI: 10.1093/nar/gkz836.
26. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alfoldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP, Gauthier LD, Brand H, Solomonson M, Watts NA, Rhodes D, Singer-Berk M, England EM, Seaby EG, Kosmicki JA, Walters RK, Tashman K, Farjoun Y, Banks E, Poter-ba T, Wang A, Seed C, Whiffin N, Chong JX, Samocha KE, Pierce-Hoffman E, Zappala Z, O'Donnell-Luria AH, Minikel EV, Weisburd B, Lek M, Ware JS, Vittal C, Armean IM, Bergelson L, Cibulskis K, Connolly KM, Covarrubias M, Donnelly S, Ferriera S, Gabriel S, Gentry J, Gupta N, Jeandet T, Kaplan D, Llanwarne C, Munshi R, Novod S, Petrillo N, Roazen D, Ruano-Rubio V, Saltzman A, Schleicher M, Soto J, Tibbetts K, Tolonen C, Wade G, Talkows-ki ME, Neale BM, Daly MJ, MacArthur DG. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434-443. DOI: 10.1038/s41586-020-2308-7.
27. Tinker SC, Gilboa SM, Moore CA, Waller DK, Simeone RM, Kim SY, Jamie-son DJ, Botto LD, Reefhuis J. Specific birth defects in pregnancies of women with
38
28. Hackshaw A, Rodeck C, Boniface S. Maternal smoking in pregnancy and birth defects, a systematic review based on 173 687 malformed cases and 11.7 million controls. Hum Reprod Update. 2011;17:589-604. DOI: 10.1093/humupd/dmr022.
diabetes: National Birth Defects Prevention Study, 1997-2011. Am J Obstet Gynecol. 2020;222:176.e1-176.e11. D01:10.1016/j.ajog.2019.08.028.
Address correspondence to:
Klyuchnikov Dmitry Yurievich Medical Centre «Dynasty», Tissue Typing Lab, 159 Tashkentskaya str., Samara, 443095, Russia, [email protected]
Received 19.04.2022 Review completed 06.06.2022 Passed for printing 08.06.2022
Dmitry Yurievich Klyuchnikov, biologist in Tissue Typing Lab, Medical Centre "Dynasty", 159 Tashkentskaya str., Samara, 443095, Russia, ORCID: 0000-0003-4934-5619, dmklyu@gmail. com;
Egor Yurievich Filatov, MD, PhD, orthopedic surgeon, National Ilizarov Medical Research Centre for Traumatology and (Orthopaedics, 6 Marii Ulyanovoj str., Kurgan, 640014, Russia, ORCID: 0000-0002-3390-807X, [email protected];
Ivan Valerievich Tyumin, clinical resident, Department of Oncology, Samara State Medical University, 89 Chapayevskaya str., Samara, 443099, Russia, ORCID: 0000-0002-7053-9285, [email protected];
Olga Vladimirovna Tyumina, DMSc, director, Medical Centre "Dynasty", 159 Tashkentskaya str., Samara, 443095, Russia; professor, Samara State Medical University, 89 Chapayevskaya str., Samara, 443099, Russia, ORCID: 0000-0002-5608-1925, [email protected].
39