Targeted imputation of sequence variants and gene expression profiling identifies twelve candidate genes associated with lactation volume,composition and calving interval in dairy cattle

Dairy cattle are an interesting model for gaining insights into the genes responsible for the large variation between and within mammalian species in the protein and fat content of their milk and their milk volume. Large numbers of phenotypes for these traits are available, as well as full genome sequence of key founders of modern dairy cattle populations. In twenty target QTL regions affecting milk production traits, we imputed full genome sequence variant genotypes into a population of 16,721 Holstein and Jersey cattle with excellent phenotypes. Association testing was used to identify variants within each target region, and gene expression data were used to identify possible gene candidates. There was statistical support for imputed sequence variants in or close to BTRC, MGST1, SLC37A1, STAT5A, STAT5B, PAEP, VDR, CSF2RB, MUC1, NCF4, and GHDC associated with milk production, and EPGN for calving interval. Of these candidates, analysis of RNA-Seq data demonstrated that PAEP, VDR, SLC37A1, GHDC, MUC1, CSF2RB, and STAT5A were highly differentially expressed in mammary gland compared to 15 other tissues. For nine of the other target regions, the most significant variants were in non-coding DNA. Genomic predictions in a third dairy breed (Australian Reds) using sequence variants in only these candidate genes were for some traits more accurate than genomic predictions from 632,003 common SNP on the Bovine HD array. The genes identified in this study are interesting candidates for improving milk production in cattle and could be investigated for novel biological mechanisms driving lactation traits in other mammals.     

A systematic review and integrative approach to decode the common molecular link between levodopa response and Parkinson’s disease

Background

PD is a progressive neurodegenerative disorder commonly treated by levodopa. The findings from genetic studies on adverse effects (ADRs) and levodopa efficacy are mostly inconclusive. Here, we aim to identify predictive genetic biomarkers for levodopa response (LR) and determine common molecular link with disease susceptibility. A systematic review for LR was conducted for ADR, and drug efficacy, independently. All included articles were assessed for methodological quality on 14 parameters. GWAS of PD were also reviewed. Protein-protein interaction (PPI) analysis using STRING and functional enrichment using WebGestalt was performed to explore the common link between LR and PD.

Results

From 37 candidate studies on levodopa toxicity, 18 genes were found associated, of which, CA_n STR 13, 14 (DRD2) was most significantly associated with dyskinesia, followed by rs1801133 (MTHFR) with hyper-homocysteinemia, and rs474559 (HOMER1) with hallucination. Similarly, 8 studies on efficacy resulted in 4 genes in which rs28363170, rs3836790 (SLC6A3) and rs4680 (COMT), were significant. To establish the molecular connection between LR with PD, we identified 35 genes significantly associated with PD. With 19 proteins associated with LR and 35 with PD, two independent PPI networks were constructed. Among the 67 nodes (263 edges) in LR, and 62 nodes (190 edges) in PD pathophysiology, UBC, SNCA, FYN, SRC, CAMK2A, and SLC6A3 were identified as common potential candidates.

Conclusion

Our study revealed the genetically significant polymorphism concerning the ADRs and levodopa efficacy. The six common genes may be used as predictive markers for therapy optimization and as putative drug target candidates.

     

Impact of variants on type-2 diabetes risk genes identified through genomewide association studies in polycystic ovary syndrome: a case–control study

Polycystic ovary syndrome (PCOS) is a common endocrine disorder in females, and is associated with altered metabolic processes in particular insulin resistance and diabetes mellitus. PCOS shares with type-2 diabetes (T2D) a number of features, including beta cell dysfunction, impaired glucose tolerance and dyslipidaemia. Recently, genomewide association studies (GWAS) have reported a number of genes with reproducible associations and susceptibilities to T2D. To address this, we examined the association between the T2D GWAS candidate genes (CDKAL1, CDKN2B, COL8A1, HHEX, IGF2BP2, KCNJ1, KCNQ1 and SLC30A8) and PCOS in Saudi women. A case–control study, includes 162 cases and 162 controls was enrolled. Genotyping was carried out by the allelic discrimination method. Our results showed that the variants including rs792837 of COL8A1, rs61873498 of KCNQ1 and rs13266634 of SLC30A8 genes to be significantly more frequent in PCOS patients than in controls. Our results suggest that COL8A1, KCNQ1 and SLC30A8, which are previously identified through GWAS as T2D-associated genes, are associated with PCOS.     

Role of allelic genes of matrix metalloproteinases and their tissue inhibitors in the risk of peptic ulcer disease development

Peptic ulcer disease is a chronic disease of the gastrointestinal tract, mainly manifesting itself in the formation of the fairly persistent ulcer defect of the mucous membrane of the stomach and/or duodenum. Association analysis of common polymorphisms of matrix metalloproteinases genes MMP-1 (rs1799750, rs494379), MMP-2 (rs2285052), MMP-3 (rs3025058), MMP-9 (rs3918242, rs17576), and MMP-12 (rs2276109) and their tissue inhibitors TIMP-2 (rs8179090) and TIMP-3 (rs9619311) was carried out in 353 patients with a gastric ulcer or duodenal ulcer and in 325 unrelated healthy individuals from the Republic of Bashkortostan. Associations of polymorphic variants rs1799750 and rs494379 of gene MMP-1, rs3025058 of gene MMP-3, rs3918242 and rs17576 of gene MMP-9, and rs9619311 of gene TIMP-3 with the risk of peptic ulcer disease in Russians and Tatars were revealed.     

Association of several polymorphic loci of serotoninergic genes with unipolar depression

Serotoninergic system is one of the major brain neurotransmitter systems that is involved in the development of depressive spectrum disorders. Regulatory genes of this system are the principle candidate genes predisposing to unipolar depression. Using PCR-RFLP analysis, we have conducted a study of polymorphic loci of several genes of this system: C1019G of serotonin receptor 1A gene, (HTR1A); A-1438G of serotonin receptor 2A gene, (HTR2A); G861C of serotonin receptor 1B gene, (HTR1B); Stin2VNTR and 5-HTTLPR of serotonin transporter gene (SLC6A4) in patients with unipolar depression from Tatar and Russian population. The results of the study suggest that genotype 10/10 of the SLC6A4 gene as well as genotype G/G and allele G of the HTR2A gene can predispose to increased risk of unipolar depression development in ethnic Russians. In contrast, genotype 12/10 of the SLC6A4 gene is a marker of low risk of the disease in both groups.     

Association of genes of different functional classes with type 1 diabetes

The genetic structure of susceptibility to type 1 diabetes (T1D) in the population of Tomsk was studied. We had a group of T1D patients (N = 285) and a population sample (N = 300) and we studied 58 SNPs localized in the 47 genes which products are involved in various metabolic pathways and processes as fibrogenesis, endothelial dysfunction, and inflammation. Genotyping was performed by mass spectrometry using the Sequenom MassARRAY system (United States). We compared the group of T1D patients and the population sample and found an association with the predisposition to disease for seven markers: rs3765124 of the ADAMDEC1 gene, genotype AA (p = 0.004), allele A (p = 0.033); rs1007856 of the ITGB5 gene, genotype TT (p = 0.015), allele T (p = 0.036); rs20579 of the LIG1 gene, genotype CC (p = 0.004), allele C (p = 0.002); rs12980602 of the IFNL2 gene, allele C (p = 0.029); rs4986819 of the PARP4 gene, allele C (p = 0.044); rs1143674 of the ITGA4 gene genotype GG (p = 0.002); rs679620 of the MMP3 gene, genotype AA (p = 0.008). Thus, the products of genes associated with T1D belong to different molecular classes: metalloproteases (ADAMDEC1, MMP3), cytokines (IL28A), cell surface receptors (ITGA4), adhesion molecules (ITGB5), DNA ligases (LIG1), and ribosyltransferase enzymes (PARP4). The ADAMDEC1, ITGA4, and ITGB5 genes belong to two biological processes: cell communication and signal transduction. The LIG1 and PARP4 genes regulate the metabolism of nucleic acids, MMP3 is involved in the regulation of protein metabolism, and the IFNL2 is involved in the immune response.     

Genetic associations between ADHD and dopaminergic genes (<Emphasis Type="Italic">DAT1</Emphasis> and <Emphasis Type="Italic">DRD4</Emphasis>) VNTRs in Korean children

It is well known that dopaminergic genes affect the development of attention deficit hyperactivity disorder (ADHD) in various populations. Many studies have shown that variable number tandem repeats (VNTRs) located within the 3′-untranslated region of DAT1 and in exon 3 of DRD4 are associated with ADHD development; however, these results were inconsistent. Therefore, we investigated the genetic association between two VNTRs and ADHD in Korean children. We determined the VNTRs using PCR. We examined genotype and allele frequency differences between the experimental and control groups, along with the odds ratios, using Chi square and exact tests. We observed a significant association between the children with ADHD and the control group in the 10R/10R genotype of DAT1 VNTRs (p?=?0.025). In addition, the 11R allele of DAT1 VNTRs showed a higher frequency in the control group than in the ADHD group (p?=?0.023). Also, the short repeat (without 11R) and long repeat alleles (including 11R) were associated with ADHD (p?<?0.05). The analysis of DRD4 VNTRs revealed that the 2R allele is associated with ADHD (p?=?0.025). A significant result was also observed in long and short repeats (p?<?0.05). Additionally, ADHD subtypes showed that the DRD4 VNTRs are associated with combined and hyperactive-impulsive subtype groups (p?<?0.05). Therefore, our results suggest that DAT1 VNTRs and DRD4 VNTRs play a role in the genetic etiology of ADHD in Korean children.     

Association of six CpG-SNPs in the inflammation-related genes with coronary heart disease

Background

Chronic inflammation has been widely considered to be the major risk factor of coronary heart disease (CHD). The goal of our study was to explore the possible association with CHD for inflammation-related single nucleotide polymorphisms (SNPs) involved in cytosine-phosphate-guanine (CpG) dinucleotides. A total of 784 CHD patients and 739 non-CHD controls were recruited from Zhejiang Province, China. Using the Sequenom MassARRAY platform, we measured the genotypes of six inflammation-related CpG-SNPs, including IL1B rs16944, IL1R2 rs2071008, PLA2G7 rs9395208, FAM5C rs12732361, CD40 rs1800686, and CD36 rs2065666). Allele and genotype frequencies were compared between CHD and non-CHD individuals using the CLUMP22 software with 10,000 Monte Carlo simulations.

Results

Allelic tests showed that PLA2G7 rs9395208 and CD40 rs1800686 were significantly associated with CHD. Moreover, IL1B rs16944, PLA2G7 rs9395208, and CD40 rs1800686 were shown to be associated with CHD under the dominant model. Further gender-based subgroup tests showed that one SNP (CD40 rs1800686) and two SNPs (FAM5C rs12732361 and CD36 rs2065666) were associated with CHD in females and males, respectively. And the age-based subgroup tests indicated that PLA2G7 rs9395208, IL1B rs16944, and CD40 rs1800686 were associated with CHD among individuals younger than 55, younger than 65, and over 65, respectively.

Conclusions

In conclusion, all the six inflammation-related CpG-SNPs (rs16944, rs2071008, rs12732361, rs2065666, rs9395208, and rs1800686) were associated with CHD in the combined or subgroup tests, suggesting an important role of inflammation in the risk of CHD.

     

Role of Gene-Gene Interactions in the Chromosomal Instability in Workers at Coal Thermal Power Plants

The 23 polymorphic variants in genes encoding the enzymes of xenobiotics biotransformation (CYP1A1 (rs4646903), CYP1A2 (rs762551), GSTP1 (rs1138272, rs1695), GSTM1 (del), and GSTT1 (del)), DNA repair (XRCC1 (rs25489, rs25487), APEX1 (rs1130409), hOGG1 (rs1052133), ADPRT (rs1136410), XPD (rs13181), XPG (rs17655), XPC (rs2228001), ATM (rs1801516), NBS (rs1805794), XRCC2 (rs3218536), and XRCC3 (rs861539)), antioxidant system (MnSOD (rs4880) and GPx1 (rs1050450)), cell cycle control and apoptosis (TP53 (rs1042522)), DNA methylation (MTHFR (rs1801133) and MTR (rs1805087)), and chromosomal aberrations in lymphocytes in the workers at thermal power plants were analyzed. We found that allelic variants in the CYP1A1 (rs4646903), hOGG1 (rs1052133), XRCC1 (rs25487), and APEX1 (rs1130409) genes were associated with an increased level of chromosomal aberrations in workers. Informative models of gene-gene interactions including CYP1A1 (rs4646903, T>C), CYP1A2 (rs762551, C>A), GSTT1 (del); XRCC1 (rs25487, G>A), MTHFR (rs1801133, C>T), GSTT1 (del); XRCC1 (rs25487, G>A), APEX1 (rs1130409, T>G), TP53 (rs1042522, G>C) determining the formation of the increased frequency of chromosomal aberrations in the workers at coal thermal power plants were discovered.     

Association of polymorphic markers of chemokine genes,their receptors,and <Emphasis Type="Italic">CD14</Emphasis> gene with coronary atherosclerosis

Atherosclerosis represents an inflammatory response to the disturbance of the endothelial layer in the arterial bloodstream. In the present study, an analysis of associations of polymorphic markers for the genes controlling synthesis of proteins involved in atherosclerosis pathogenesis in coronary atherosclerosis (CA) patients (217 subjects) and in a control group (250 subjects) was conducted. The following genes were examined: rs991804 (CCL2 gene), rs1126579 (CXCR2 gene), rs4074 (CXCL1 gene), rs4073 (CXCL8 gene), rs333 (CCR5 gene), rs2471859 (CXCR4 gene), rs1801157 (CXCL12 gene), and rs2569190 (CD14 gene). Using the Monte Carlo and Markov chain (APSampler) method, allele/genotype combinations associated with both low and high CA risk were revealed. The most important findings included the following: CXCR4*T/T + CCL2*C + CCR5*I/I (P_perm = 1 × 10^–6, OR = 0.44, 95% CI 0.3–0.63), CXCR2*C + CD14*C + CXCL12*G + CCL2*C + CCR5*D (P_perm = 4 × 10^–6, OR = 5.78, 95% CI 2.34–14.28), CD14*C + CCL2*C/C + CCR5*D (P_perm = 6.3 × 10^–6, OR = 5.81, 95% CI 2.17–15.56), CXCL8*A + CXCR2*C + CD14*T + CXCR4*C (P_perm = 0.01, OR = 3.21, 95% CI 1.63–6.31).     

Combined effect between two functional polymorphisms of SLC6A12 gene is associated with temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is the most common epilepsy subtype with complex genetic structure. A recent study in four populations (Ireland, UK, Australia and Finland) reported an allelic association between betaine/GABA transporter-1 (BGT-1 or SLC6A12) and mesial temporal lobe epilepsy with hippocampal sclerosis. To demonstrate the association between SLC6A12 gene polymorphisms and TLE, TaqMan method was used to genotype five single-nucleotide polymorphisms of SLC6A12 gene in 358 TLE patients and 596 nonepileptic control subjects of Chinese Han origin. Real-time PCR was used to detect the effects of variations on gene expression associated with TLE. Though, the single-marker analysis did not demonstrate allelic association with TLE, rs542736–rs557881 interaction showed significant association. The SLC6A12 expression levels in peripheral blood mononuclear cells were significantly higher in TLE patients than in control subjects and were correlated to rs542736 G–rs557881 A haplotypes. Our preliminary results suggested combined effect of two common polymorphisms on SLC6A12 gene may be associated with TLE, but the precise mechanism needs further investigation.     

Polymorphic Imprinting of <Emphasis Type="Italic">SLC38A4</Emphasis> Gene in Bovine Placenta

Imprinted genes are characterized by monoallelic expression that is dependent on parental origin. Comparative analysis of imprinted genes between species is a powerful tool for understanding the biological significance of genomic imprinting. The slc38a4 gene encodes a neutral amino acid transporter and is identified as imprinted in mice. In this study, the imprinting status of SLC38A4 was assessed in bovine adult tissues and placenta using a polymorphism-based approach. Results indicate that SLC38A4 is not imprinted in eight adult bovine tissues including heart, liver, spleen, lung, kidney, muscle, fat, and brain. It was interesting to note that SLC38A4 showed polymorphic status in five heterogeneous placentas, with three exhibiting paternal monoallelic expression and two exhibiting biallelic expression. Monoallelic expression of imprinted genes is generally associated with allele-specific differentially methylation regions (DMRs) of CpG islands (CGIs)-encompassed promoter; therefore, the DNA methylation statuses of three CGIs in the SLC38A4 promoter and exon 1 region were tested in three placentas (two exhibiting paternal monoallelic and one showing biallelic expression of SLC38A4) and their corresponding paternal sperms. Unexpectedly, extreme hypomethylation (<?3%) of the DNA was observed in all the three detected placentas and their corresponding paternal sperms. The absence of DMR in bovine SLC38A4 promoter region implied that DNA methylation of these three CGIs does not directly or indirectly affect the polymorphic imprinting of SLC38A4 in bovine placenta. This suggested other epigenetic features other than DNA methylation are needed in regulating the imprinting of bovine SLC38A4, which is different from that of mouse with respect to a DMR existence at the mouse’s slc38a4 promoter region. Although further work is needed, this first characterization of polymorphic imprinting status of SLC38A4 in cattle placenta provides valuable information on investigating the genomic imprinting phenomenon itself.     

Polymorphic variants of glutamate receptor (<Emphasis Type="Italic">GRIK5</Emphasis>, <Emphasis Type="Italic">GRIN2B</Emphasis>) and serotonin receptor (<Emphasis Type="Italic">HTR2A</Emphasis>) genes are associated with chronic

Chronic obstructive pulmonary disease (COPD) is a complex chronic inflammatory disease of the respiratory system that affects primarily distal respiratory pathways and lung parenchyma. Smoking tobacco is a major risk factor for COPD. The relationship of HTR4 (rs3995090), HTR2A (rs6313), GRIK5 (rs8099939), GRIN2B (rs2268132), and CHRNB4 (rs1948) gene polymorphisms and COPD, as well as the contribution of these polymorphisms to the variations in quantitative characteristics that describe respiratory function, smoking behavior, and nicotine dependence was assessed in an ethnically homogeneous Tatar population. The polymorphisms of HTR2A (rs6313) (P = 0.026, OR = 1.42 for the CC genotype) and GRIN2B (rs2268132) (P = 0.0001, OR = 2.39 for the TT genotype) were significantly associated with increased risk of COPD. The AA genotype of GRIK5 (rs8099939) had a protective effect (P = 0.02, OR = 0.61). Importantly, the HTR2A (rs6313), GRIN2B (rs2268132), and GRIK5 (rs8099939) polymorphisms were only associated with COPD in smokers. Smoking index (pack-years) was significantly higher in carriers of the GRIK5 genotype AC (rs8099939) (P = 0.0027). The TT genotype of GRIN2B (rs2268132) was associated with COPD in subjects with high nicotine dependence according to the Fagerström test (P = 0.002, OR = 2.98). The TT genotype of HTR2A (rs6313) was associated with a reduced risk of the disease in the group with moderate nicotine dependence (P = 0.02, OR = 0.22). The CC genotype of HTR2A (rs6313) and the TT genotype of GRIN2B (rs2268132) were associated with higher levels of nicotine dependence according to the Fagerström test (P = 0.0011 and P = 0.037). Our results may provide insight into potential molecular mechanisms that involve the glutamate (GRIK5, GRIN2B) and serotonin (HTR2A) receptor genes in the pathogenesis of COPD.     

No genetic association between attention-deficit/hyperactivity disorder (ADHD) and Parkinson’s disease in nine ADHD candidate SNPs

Attention-deficit/hyperactivity disorder (ADHD) and Parkinson’s disease (PD) involve pathological changes in brain structures such as the basal ganglia, which are essential for the control of motor and cognitive behavior and impulsivity. The cause of ADHD and PD remains unknown, but there is increasing evidence that both seem to result from a complicated interplay of genetic and environmental factors affecting numerous cellular processes and brain regions. To explore the possibility of common genetic pathways within the respective pathophysiologies, nine ADHD candidate single nucleotide polymorphisms (SNPs) in seven genes were tested for association with PD in 5333 cases and 12,019 healthy controls: one variant, respectively, in the genes coding for synaptosomal-associated protein 25 k (SNAP25), the dopamine (DA) transporter (SLC6A3; DAT1), DA receptor D4 (DRD4), serotonin receptor 1B (HTR1B), tryptophan hydroxylase 2 (TPH2), the norepinephrine transporter SLC6A2 and three SNPs in cadherin 13 (CDH13). Information was extracted from a recent meta-analysis of five genome-wide association studies, in which 7,689,524 SNPs in European samples were successfully imputed. No significant association was observed after correction for multiple testing. Therefore, it is reasonable to conclude that candidate variants implicated in the pathogenesis of ADHD do not play a substantial role in PD.     

Dysfunctional protection against advanced glycation due to thiamine metabolism abnormalities in gestational diabetes

While the pathogenic role of dicarbonyl stress and accelerated formation of advanced glycation end products (AGEs) to glucose intolerance and to the development of diabetic complications is well established, little is known about these processes in gestational diabetes mellitus (GDM), a condition pathogenically quite similar to type 2 diabetes. The aims of the present study were (i) to determine plasma thiamine and erythrocyte thiamine diphosphate (TDP) and transketolase (TKT) activity in pregnant women with and without GDM, (ii) to assess relationships between thiamine metabolism parameters and selected clinical, biochemical and anthropometric characteristics and, finally, (iii) to analyse relationship between variability in the genes involved in the regulation of transmembrane thiamine transport (i.e. SLC19A2 and SLC19A3) and relevant parameters of thiamine metabolism. We found significantly lower plasma BMI adjusted thiamine in women with GDM (P = 0.002, Mann-Whitney) while levels of erythrocyte TDP (an active TKT cofactor) in mid-trimester were significantly higher in GDM compared to controls (P = 0.04, Mann-Whitney). However, mid-gestational TKT activity - reflecting pentose phosphate pathway activity - did not differ between the two groups (P > 0.05, Mann-Whitney). Furthermore, we ascertained significant associations of postpartum TKT activity with SNPs SLC19A2 rs6656822 and SLC19A3 rs7567984 (P = 0.03 and P = 0.007, resp., Kruskal-Wallis). Our findings of increased thiamine delivery to the cells without concomitant increase of TKT activity in women with GDM therefore indicate possible pathogenic role of thiamine mishandling in GDM. Further studies are needed to determine its contribution to maternal and/or neonatal morbidity.     

MicroRNA Biogenesis Pathway Gene Polymorphisms Are Associated with Breast Cancer Risk

MicroRNAs (miRNAs) play an important role as epigenetic regulators in cancer initiation and progression. One of the mechanisms of miRNA dysregulation is altered functioning of proteins involved in miRNA processing machinery. It has been suggested that single nucleotide polymorphisms (SNPs) within miRNA gene regions, miRNA target genes, and miRNA machinery genes may affect the miRNAs regulation. We selected 25 SNPs in the key genes of miRNA biosynthesis, including DROSHA/RNASEN, DGCR8, DICER1, XPO5, RAN, PIWIL1/HIWI, AGO1/EIF2C1, AGO2, GEMIN4, GEMIN3/DDX20, and DDX5, and investigated the association between these SNPs and the risk of breast cancer. The total number of breast cancer cases and cancer-free controls enrolled in the investigation were 778 (417 breast cancer patients and 361 healthy women). We found that rs11060845 and rs10773771 in the PIWIL1 gene, rs3809142/RAN, rs10719/DROSHA, rs1640299/DGCR8, rs563002/DDX20, rs595055/AGO1, and rs2740348/GEMIN4 were associated with breast cancer risk in Russians.     

Search for osteoarthritis genetic markers in women with undifferentiated connective tissue dysplasia

We conducted an association study of ten polymorphisms in six candidate genes of OA (rs1799750 (MMP1), rs35068180 (MMP3), rs2252070 (MMP13), rs63118460 and rs2276455 (COL2A1), rs143383 (GDF5), rs1544410, rs7975232, rs731236, and rs2228570 (VDR)) with the development of osteoarthritis (OA) in 333 women taking into account the localization of the pathological process, the age of disease manifestation, and ethnicity and investigated the presence of signs of undifferentiated connective tissue dysplasia (UCTD). On the basis of clinical and genetic data, we revealed statistically significant models to predict the development of osteoarthritis of various localizations (knee, hip, and generalized forms).     

Association of gene polymorphisms of matrix metalloproteinases with reproductive losses in the first trimester of pregnancy

In the present study, the frequencies of genotypes and alleles of candidate genes with respect to polymorphisms associated with increased pregnancy loss in the first trimester of pregnancy, including MMP1–1607insG, MMP9 A–8202G, and TIMP1 С536T, were reported. The frequency of homozygotes for allele MMP9 A–8202 was increased by a factor of two among women with miscarriage in the first trimester compared to the control. Significant models of interaction of genes MMPs and TIMP1 were revealed. The genotypes of genes MMP1 (rs1799750), MMP9 (rs11697325), and TIMP1 (rs11551797) increasing the risk of pregnancy loss in the first trimester were determined.