Skip to content

Advertisement

Reproductive Biology and Endocrinology

What do you think about BMC? Take part in

Open Access

THBD sequence variants potentially related to recurrent pregnancy loss

  • Paula Quintero-Ronderos1,
  • Eric Mercier2, 3,
  • Jean-Christophe Gris2, 3,
  • Clara Esteban-Perez4,
  • Harold Moreno-Ortiz4,
  • Dora Janeth Fonseca1,
  • Elkin Lucena4,
  • Daniel Vaiman5, 6 and
  • Paul Laissue1Email author
Reproductive Biology and Endocrinology201715:92

https://doi.org/10.1186/s12958-017-0311-0

Received: 20 September 2017

Accepted: 22 November 2017

Published: 1 December 2017

Abstract

Recurrent pregnancy loss (RPL) is a frequently occurring disease, which is classified as idiopathic in more than 50% of cases. THBD, the endothelial cell receptor for thrombin, has been associated with distinct biological processes and considered a coherent RPL-related candidate gene. In the present study, we have sequenced the complete coding region of THBD in 262 patients affected by RPL. Bioinformatics analysis and screening of controls strongly suggested that the THBD-p.Trp153Gly mutation might be related to RPL aetiology. It could be used, after its validation by functional assays, as a molecular marker for diagnostic/prognostic purposes.

Keywords

Recurrent pregnancy lossTHBDMolecular markerFemale infertility

Introduction

Recurrent pregnancy loss (RPL), which is clinically defined as at least two consecutive pregnancy losses prior to 20 weeks of gestation, is a frequently occurring disease. Although several RPL’s aetiologies have been described, the disease is still classified as idiopathic in more than 50% of cases, which underlines the potential involvement of genetic factors. However, only a few studies have revealed links between specific genetic variants and RPL’s origin.

THBD, the endothelial cell receptor for thrombin, has been associated with distinct biological processes and considered a coherent RPL-related candidate gene [1]. It has been largely studied regarding its involvement in the natural anticoagulation system. Its binding to thrombin leads to decreased fibrinogen-to-fibrin conversion and activation of the coagulation cascade (e.g. factors V, VIII, platelets). The THBD/thrombin complex catalyses protein C activation which, in turn, increases inhibition of the coagulation effect via cofactor Va and VIIIa inactivation. These features have been well documented in vivo as Thbd knock-out mice have displayed embryonic death secondary to coagulation dysfunction at the foetal-maternal interface [2]. A THBD truncating mutation in humans (p.Cys537X), which disturbs transmembrane domain conformation and has eliminated the cytoplasmic domain, has been functionally related to an autosomal dominant bleeding disorder [3]. Regarding human pregnancy dysfunction, it has been postulated that RPL might be associated with relative resistance to THBD as well as to its decrease in placental tissue. However, studying a potential association between some THBD SNP and RPL has led to contradictory results.

In the present study, we have sequenced the complete coding region of THBD in 262 patients affected by RPL. Bioinformatics analysis and screening of controls strongly suggested that the THBD-p.Trp153Gly mutation might be related to RPL aetiology. It could be used, after their validation by functional assays, as a molecular marker for diagnostic/prognostic purposes.

Methods

Patients and controls

Two hundred and sixty two women affected by RPL were included in the study. 233 and 29 were of French (Caucasian) and Colombian (Mestizo) origin, respectively. French women affected by RPL formed part of a previously-established group (the Nîmes Obstetricians and Haematologists patient cohort) of women suffering pregnancy loss [4]. Colombian women suffering RPL attended the Center For Research in Genetics and Genomics (Universidad del Rosario, Bogotá, Colombia). They presented normal karyotypes, and were affected by at least 2 pregnancy losses. They did not display antecedents of coagulation, autoimmune (e.g. antiphospholipid syndrome) and/or metabolic disorders. Autoimmune and coagulation disorders were excluded by biochemical tests. Uterine morphology was normal. All cases lack familial antecedents of consanguinity or RPL. No extensive genotyping initiatives at genome scale have been carried out in Colombian population. Thus, to screen interesting variants identified in RPL and Colombian patients, we have recruited 165 Colombian (Mestizo) women under 50 years old having at least one child and lacking antecedents of medical complications during pregnancy. We did not include Caucasian control women for THBD sequencing because public databases of SNPs are enriched by data of individuals from this ethnical origin. All the experimental steps of this study were approved by the Universidad del Rosario’s Ethics Committee, and was conducted in line with the Declaration of Helsinki (approval date: 1 February 2017. Institional Review Board reference number: CEI-DVN021–1-063).

THBD molecular and bioinformatics analysis

We extracted genomic DNA from whole blood samples by using the standard salting-out procedure. Then, DNA was used to amplify (PCR) the complete THBD open reading frame. Amplicons were purified by using shrimp alkaline phosphatase and exonuclease I. Sequencing involved using internal primers. Sequences were compared to the THBD wild type version (ENSG00000178726). Primer sequences, PCR and sequencing technical conditions are available upon request.

The frequency of THBD sequence variants in distinct healthy populations was screened using the Ensembl database. Allele frequencies of the c.1418 C > T (p.Ala473Val) variant in European controls were obtained from the Ensembl database and compared to that found in our French RPL women. Allele frequencies of variants identified in RPL Colombian patients were compared to those from a control group (n = 165) from the same ethnical origin. SIFT and PolyPhen bioinformatics tools were used for predicting potential deleterious effects of missense variants. Protein alignments between species were performed to assess the conservation during evolution of interchanged residues.

Statistical analysis

We performed chi-square tests for calculating significant differences of allele frequencies between cases and controls individuals.

Results

RPL patients were divided in different subgroups (G1 through G4) according to the number of miscarriages: G1 (2 RPLs), G2 (3 RPLs), G3 (4 RPLs) and G4 (> 4 RPLs). Most RPL patients (59.9%, n = 157) belonged to the G1. 26.3% (n = 69) and 6.9% (n = 18) of patients were classified into the G2 and G3 groups, respectively. G4 included 18 RPL patients (6.9%). Screening for mutations of the THBD gene displayed 2 non-synonymous sequence variants: c.457 T > G (p.Trp153Gly) and c.1418 C > T (p.Ala473Val) (Table 1). THBD-c.457 T > G and THBD-c.1418 C > T displayed, in the SNP Ensembl database, minor allele frequencies (MAF) of 0.00 and 0.16, respectively. For the c.457 T > G variant, which was found uniquely in Colombian patients, statistically significant differences in allele frequencies were identified between cases and controls from the same ethnical origin (p = 9 × 10−6). For the c.1418 C > T variant, a statistically significant difference (p = 0.02) was identified between French RPL patients and European controls recorded in the Ensembl database. Conversely, no differences in allelic frequencies for the c.1418 C > T variant were identified between Colombian RPL patients and control individuals from the same ethnical origin. The c.457 T > G and c.1418 C > T variants were most frequently found in the G1, in both French and Colombian RPL women (Table 1).
Table 1

THBD open reading frame sequencing in RPL patients and control individuals

French RPL patients

DNA

Protein

Zygocity (n)

Rs

RPL Group

Patients

European controls (Ensembl)

MAF

Allelic Frequency (cases)

Allelic Frequency (controls)

P value

SIFT

PolyPhen

Evolutive conservation

   c.1418 C > T

p.Ala473Val

Homozygote (3)

rs1042579

G1: 65,6%

32/233

173/503

0.16

C: 0,92 T: 0,08

C: 0,80 T: 0,20

0,02

Tolerated (0,42)

Bening (0)

Yes

Heterozygote (29)

G2: 25%

 

G3: 9,4%

Colombian RPL patients

DNA

Protein

Zygocity (n)

Rs

RPL Group

Patients

Colombian Controls

MAF

Allelic Frequency (cases)

Allelic Frequency (controls)

P value

SIFT

PolyPhen

Evolutive conservation

   c. 457 T > G

p.Trp153Gly

Homozygote (3)

 

G1: 50%

4/29

0/165

 

T: 0,88 G: 0,12

T: 1 G: 0

0,000009

Deleterious (0)

Probably damaging (0.976)

Yes

Heterozygote (1)

G3: 25%

 

G4: 25%

   c.1418 C > T

p.Ala473Val

Heterozygote (4)

rs1042579

G1: 50%

4/29

51/165

0.16

C: 0,93 T: 0,07

C: 0,95 T: 0,05

0,73

Tolerated (0,42)

Bening (0)

Yes

G2: 25%

G4: 25%

Discussion

In the present study, we focused our attention on the research of rare THBD non-synonymous mutations underlying theoretically putative moderate/strong functional effects contributing to RPL. We consider that p.Trp153Gly is particularly interesting as low MAF was reported in public databases of SNP and significant statistical differences were found in allele frequencies between cases and controls of Colombian origin. The p.Trp153Gly mutation involved the exchange of an aminoacid conserved during evolution (Trp153), and Trp to Gly exchange evoked a drastic modification of protein’s local physicochemical properties. In addition, SIFT and Polyphen scores were compatible with deleterious effects. These features might be related to an improper folding of the protein leading to functional disturbances.

Interestingly, this mutation was located at the end of the C-type lectin-like (CLL) module, which has been associated with protective properties against inflammatory damage [5]. Due to its location at the N-terminal extracellular region of the protein, it would seem proper for the CLL domain to interact with other molecules and cells. Transgenic mice lacking the Thbd-CLL region have displayed enhanced synthesis of cytokines and a significant inflammation reaction when treated with lipopolysaccharide [5]. It has also been proposed that the THBD-CLL domain participates in cell-cell adhesion by maintaining the integrity of endothelial junctions and contributing to maintaining blood vessels’ quiescent state [6, 7]. CLL’s anti-inflammatory properties have been related to HMGB1 (high mobility group box 1) functions, such molecule originating in necrotic and inflammatory cells and displaying pro-inflammatory mediator properties [5]. The CLL domain of THBD binds directly to HMGB1, interfering with its binding to specific receptors and decreasing inflammation. We thus hypothesise that the THBD-p.Trp153Gly mutation might perturb the interaction between THBD and HMGB1, thereby leading to local enhanced inflammation at the foetal-placental interface and contributing towards miscarriage aetiology. We cannot discard a decrease of mutant THBD bioavailability at the cell surface secondary to abnormal protein folding and cytoplasmic trafficking. We consider therefore that the p.Trp153Gly mutation might be a relevant molecular marker for RPL in Colombian population. Concerning the c.1418 C > T (p.Ala473Val) variant (located in the epidermal growth-factor–like repeats region, EGF-LR), it is difficult to associate it with a potential deleterious effect since, although statistical differences were identified in allelic frequencies between cases and controls from French origin, no differences were observed in individuals from Colombian origin. A lack of association of this variant with the RPL phenotype was reported by Kaare et al., 2007 and Guerra-Shinohara et al., 2012 [8, 9]. Conversely, Cao et al., 2013 reported that this variant was significantly associated with the increased risk of RPL in Chinese population [10]. These differences might be due to a small population size in some studies and/or secondary to genetic ethnically specific differences. Thus, we cannot discard that the c.1418 C > T variant confers susceptibility to RPL. To note, this variant has been associated with thrombosis [1115].

The mutations found in the present study might alter different mechanisms related to complement activation, coagulation, fibrinolysis, inflammation and cell/cell adhesion. Dysfunction of the CLL domain might be related to enhanced susceptibility to inflammation and alterations in the trophoblast/endometrium adhesion. Mutations in the EGF-LR domain might lead to expression disturbances of key molecules such as the activated protein C and TAFI, promoting an hypercoagulation state.

In our study the c.457 T > G and c.1418 C > T variants were most frequently identified in women affected by 2 miscarriages. However, due to the complexity of the RPL phenotype, in which numerous genes (and mutations) are probably involved, it was difficult to establish an accurate genotype-phenotype correlation. For clinicians, these variants might be used in the near future to establish more accurate treatment for women affected by RPL and genetic counselling.

Conclusions

In synthesis our findings argues in favour of an aetiological link between THBD sequence variants and RPL aetiology. We consider that the THBD-p.Trp153Gly mutation might be related to RPL aetiology in Colombian RPL women. It could be used, after its validation by functional assays, as a molecular marker for diagnostic/prognostic purposes.

Abbreviations

CLL: 

C-type lectin-like

EGF-LR: 

epidermal growth-factor–like repeats region.

HMGB1: 

high mobility group box 1

MAF: 

minor allelic frequency

RPL: 

recurrent pregnancy loss

SNP: 

single nucleotide polymorphism

THBD: 

thrombomodulin

Declarations

Acknowledgments

Not applicable

Funding

This study was supported by the Universidad del Rosario (grant: CS/ABN062/GENIUROS 017).

Availability of data and materials

The data analyzed during the current study are available from the corresponding author on reasonable request.

Authors’ contributions

JCG and EM performed the clinical activities. PQ, EM, DV, JCG, DF, CE, HM, EL, PL analysed the data and revised the paper. PL designed and directed the study, and wrote the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

All of this study’s clinical and experimental steps were approved by the Institutional Ethics Committee of both participating institutions (number: DVN021–1-063). The clinical investigation was performed according to the Helsinki Declaration of 1975, as revised in 1996. All of the women had given their informed consent to participate.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Center For Research in Genetics and Genomics-CIGGUR, GENIUROS Research Group, School of Medicine and Health Sciences, Universidad del Rosario
(2)
Department of Haematology, University Hospital
(3)
Faculty of Pharmacy and Biological Sciences and Research Team EA 2992, University of Montpellier
(4)
Department of Reproductive Genetics, Fertility and Sterility Colombian Center
(5)
Institut Cochin, Université Paris Descartes, CNRS (UMR 8104)
(6)
Inserm, U1016

References

  1. Vaiman D. Genetic regulation of recurrent spontaneous abortion in humans. Biom J. 2015;38:11.Google Scholar
  2. Isermann B, Sood R, Pawlinski R, Zogg M, Kalloway S, Degen JL, et al. The thrombomodulin–protein C system is essential for the maintenance of pregnancy. Nat Med. 2003;9:331–7.View ArticlePubMedGoogle Scholar
  3. Maclachlan A, Dolan G, Grimley C, Watson SP, Morgan NV, on behalf of the UK GAPP Study Group. Whole exome sequencing identifies a mutation in thrombomodulin as the genetic cause of a suspected platelet disorder in a family with normal platelet function. Platelets. 2017;6:611-3.Google Scholar
  4. Lissalde-Lavigne G, Fabbro-Peray P, Cochery-Nouvellon E, Mercier E, Ripart-Neveu S, Balducchi J-P, et al. Factor V Leiden and prothrombin G20210A polymorphisms as risk factors for miscarriage during a first intended pregnancy: the matched case-control “NOHA first” study. J Thromb Haemost. 2005;3:2178–84.View ArticlePubMedGoogle Scholar
  5. Li Y-H, Kuo C-H, Shi G-Y, Wu H-L. The role of thrombomodulin lectin-like domain in inflammation. J Biomed Sci. 2012;19:34.View ArticlePubMedPubMed CentralGoogle Scholar
  6. Huang H-C, Shi G-Y, Jiang S-J, Shi C-S, C-M W, Yang H-Y, et al. Thrombomodulin-mediated cell adhesion: involvement of its lectin-like domain. J Biol Chem. 2003;278:46750–9.View ArticlePubMedGoogle Scholar
  7. Li Y-H, Shi G-Y, Wu H-L. The role of thrombomodulin in atherosclerosis: from bench to bedside. Cardiovasc Hematol Agents Med Chem. 2006;4:183–7.View ArticlePubMedGoogle Scholar
  8. Kaare M. Ulander V-MM, painter JN, Ahvenainen T, Kaaja R, Aittomäki K. Variations in the thrombomodulin and endothelial protein C receptor genes in couples with recurrent miscarriage. Hum Reprod. 2007;22:864–8.View ArticlePubMedGoogle Scholar
  9. Guerra-Shinohara EM, Bertinato JF, Tosin Bueno C, Cordeiro da Silva K, Burlacchini de Carvalho MH. Pulcineli Vieira Francisco R, et al. polymorphisms in antithrombin and in tissue factor pathway inhibitor genes are associated with recurrent pregnancy loss. Thromb Haemost. 2012;108:693–700.View ArticlePubMedGoogle Scholar
  10. Cao Y, Zhang Z, Xu J, Yuan W, Wang J, Huang X, et al. The association of idiopathic recurrent pregnancy loss with polymorphisms in hemostasis-related genes. Gene. 2013;530:248–52.View ArticlePubMedGoogle Scholar
  11. Navarro S, Medina P, Bonet E, Corral J, Martinez-Sales V, Martos L, et al. Association of the Thrombomodulin Gene c.1418C>T Polymorphism With Thrombomodulin Levels and With Venous Thrombosis Risk. Arterioscler Thromb Vasc Biol. 2013;33:1435–40.View ArticlePubMedGoogle Scholar
  12. Aleksic N, Folsom AR, Cushman M, Heckbert SR, Tsai MY, Prospective WKK. Study of the A455V polymorphism in the thrombomodulin gene, plasma thrombomodulin, and incidence of venous thromboembolism: the LITE study. J Thromb Haemost. 2003;1:88–94.View ArticlePubMedGoogle Scholar
  13. Heit JA, Petterson TM, Owen WG, Burke JP, DE Andrade M, Melton LJ. Thrombomodulin gene polymorphisms or haplotypes as potential risk factors for venous thromboembolism: a population-based case-control study. J Thromb Haemost. 2005;3:710–7.View ArticlePubMedGoogle Scholar
  14. Sugiyama S, Hirota H, Kimura R, Kokubo Y, Kawasaki T, Suehisa E, et al. Haplotype of thrombomodulin gene associated with plasma thrombomodulin level and deep vein thrombosis in the Japanese population. Thromb Res. 2007;119:35–43.View ArticlePubMedGoogle Scholar
  15. Tang L, Wang H-F, Lu X, Jian X-R, Jin B, Zheng H, et al. Common genetic risk factors for venous thrombosis in the Chinese population. Am J Hum Genet. 2013;92:177–87.View ArticlePubMedPubMed CentralGoogle Scholar

Copyright

© The Author(s). 2017

Advertisement