SMIM1 variants rs1175550 and rs143702418 independently modulate Vel blood group antigen expression

A Bioinformatically Initiated Approach to Evaluate GATA1 Regulatory Regions in Samples with Weak D, Del, or D- Phenotypes Despite Normal RHD Exons

Introduction

With over 360 blood group antigens in systems recognized, there are antigens, such as RhD, which demonstrate a quantitative reduction in antigen expression due to nucleotide variants in the non-coding region of the gene that result in aberrant splicing or a regulatory mechanism. This study aimed to evaluate bioinformatically predicted GATA1-binding regulatory motifs in the RHD gene for samples presenting with weak or apparently negative RhD antigen expression but showing normal RHD exons.

Methods

Publicly available open chromatin region data were overlayed with GATA1 motif candidates in RHD. Genomic DNA from weak D, Del or D- samples with normal RHD exons (n = 13) was used to confirm RHD zygosity by quantitative PCR. Then, RHD promoter, intron 1, and intron 2 regions were amplified for Sanger sequencing to detect potential disruptions in the GATA1 motif candidates. Electrophoretic mobility shift assay (EMSA) was performed to assess GATA1-binding. Luciferase assays were used to assess transcriptional activity.

Results

Bioinformatic analysis identified five of six GATA1 motif candidates in the promoter, intron 1 and intron 2 for investigation in the samples. Luciferase assays showed an enhancement in transcription for GATA1 motifs in intron 1 and for intron 2 only when the R 2 haplotype variant (rs675072G>A) was present. GATA1 motifs were intact in 12 of 13 samples. For one sample with a Del phenotype, a novel RHD c.1-110A>C variant disrupted the GATA1 motif in the promoter which was supported by a lack of a GATA1 supershift in the EMSA and 73% transcriptional activity in the luciferase assay. Two samples were D+/D- chimeras.

Conclusion

The bioinformatic predictions enabled the identification of a novel DEL allele, RHD c.1-110A>C, which disrupted the GATA1 motif in the proximal promoter. Although the majority of the samples investigated here remain unexplained, we provide GATA1 targets which may benefit future RHD regulatory investigations.

Epigenetic dissection of human blood group genes reveals regulatory elements and detailed characteristics of KEL and four other loci

BACKGROUND

Blood typing is essential for safe transfusions and is performed serologically or genetically. Genotyping predominantly focuses on coding regions, but non-coding variants may affect gene regulation, as demonstrated in the ABO, FY and XG systems. To uncover regulatory loci, we expanded a recently developed bioinformatics pipeline for discovery of non-coding variants by including additional epigenetic datasets.

METHODS

Multiple datasets including ChIP-seq with erythroid transcription factors (TFs), histone modifications (H3K27ac, H3K4me1), and chromatin accessibility (ATAC-seq) were analyzed. Candidate regulatory regions were investigated for activity (luciferase assays) and TF binding (electrophoretic mobility shift assay, EMSA, and mass spectrometry, MS).

RESULTS

In total, 814 potential regulatory sites in 47 blood-group-related genes were identified where one or more erythroid TFs bound. Enhancer candidates in CR1, EMP3, ABCB6, and ABCC4 indicated by ATAC-seq, histone markers, and co-occupancy of 4 TFs (GATA1/KLF1/RUNX1/NFE2) were investigated but only CR1 and ABCC4 showed increased transcription. Co-occupancy of GATA1 and KLF1 was observed in the KEL promoter, previously reported to contain GATA1 and Sp1 sites. TF binding energy scores decreased when three naturally occurring variants were introduced into GATA1 and KLF1 motifs. Two of three GATA1 sites and the KLF1 site were confirmed functionally. EMSA and MS demonstrated increased GATA1 and KLF1 binding to the wild-type compared to variant motifs.

DISCUSSION

This combined bioinformatics and experimental approach revealed multiple candidate regulatory regions and predicted TF co-occupancy sites. The KEL promoter was characterized in detail, indicating that two adjacent GATA1 and KLF1 motifs are most crucial for transcription.

Elucidation of the low-expressing erythroid CR1 phenotype by bioinformatic mining of the GATA1-driven blood-group regulome

Genetic determinants underlying most human blood groups are now clarified but variation in expression levels remains largely unexplored. By developing a bioinformatics pipeline analyzing GATA1/Chromatin immunoprecipitation followed by sequencing (ChIP-seq) datasets, we identify 193 potential regulatory sites in 33 blood-group genes. As proof-of-concept, we aimed to delineate the low-expressing complement receptor 1 (CR1) Helgeson phenotype on erythrocytes, which is correlated with several diseases and protects against severe malaria. We demonstrate that two candidate CR1 enhancer motifs in intron 4 bind GATA1 and drive transcription. Both are functionally abolished by naturally-occurring SNVs. Erythrocyte CR1-mRNA and CR1 levels correlate dose-dependently with genotype of one SNV (rs11117991) in two healthy donor cohorts. Haplotype analysis of rs11117991 with previously proposed markers for Helgeson shows high linkage disequilibrium in Europeans but explains the poor prediction reported for Africans. These data resolve the longstanding debate on the genetic basis of inherited low CR1 and form a systematic starting point to investigate the blood group regulome.

Functional dissection of inherited non-coding variation influencing multiple myeloma risk

Thousands of non-coding variants have been associated with increased risk of human diseases, yet the causal variants and their mechanisms-of-action remain obscure. In an integrative study combining massively parallel reporter assays (MPRA), expression analyses (eQTL, meQTL, PCHiC) and chromatin accessibility analyses in primary cells (caQTL), we investigate 1,039 variants associated with multiple myeloma (MM). We demonstrate that MM susceptibility is mediated by gene-regulatory changes in plasma cells and B-cells, and identify putative causal variants at six risk loci (SMARCD3, WAC, ELL2, CDCA7L, CEP120, and PREX1). Notably, three of these variants co-localize with significant plasma cell caQTLs, signaling the presence of causal activity at these precise genomic positions in an endogenous chromosomal context in vivo. Our results provide a systematic functional dissection of risk loci for a hematologic malignancy.

Cataloguing experimentally confirmed 80.7 kb-long ACKR1 haplotypes from the 1000 Genomes Project database

Background

Clinically effective and safe genotyping relies on correct reference sequences, often represented by haplotypes. The 1000 Genomes Project recorded individual genotypes across 26 different populations and, using computerized genotype phasing, reported haplotype data. In contrast, we identified long reference sequences by analyzing the homozygous genomic regions in this online database, a concept that has rarely been reported since next generation sequencing data became available.

Study design and methods

Phased genotype data for a 80.6 kb region of chromosome 1 was downloaded for all 2,504 unrelated individuals of the 1000 Genome Project Phase 3 cohort. The data was centered on the ACKR1 gene and bordered by the CADM3 and FCER1A genes. Individuals with heterozygosity at a single site or with complete homozygosity allowed unambiguous assignment of an ACKR1 haplotype. A computer algorithm was developed for extracting these haplotypes from the 1000 Genome Project in an automated fashion. A manual analysis validated the data extracted by the algorithm.

Results

We confirmed 902 ACKR1 haplotypes of varying lengths, the longest at 80,584 nucleotides and shortest at 1,901 nucleotides. The combined length of haplotype sequences comprised 19,895,388 nucleotides with a median of 16,014 nucleotides. Based on our approach, all haplotypes can be considered experimentally confirmed and not affected by the known errors of computerized genotype phasing.

Conclusions

Tracts of homozygosity can provide definitive reference sequences for any gene. They are particularly useful when observed in unrelated individuals of large scale sequence databases. As a proof of principle, we explored the 1000 Genomes Project database for ACKR1 gene data and mined long haplotypes. These haplotypes are useful for high throughput analysis with next generation sequencing. Our approach is scalable, using automated bioinformatics tools, and can be applied to any gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04169-6.

SMIM1, carrier of the Vel blood group, is a tail-anchored transmembrane protein and readily forms homodimers in a cell-free system

Antibodies to the Vel blood group antigen can cause adverse hemolytic reactions unless Vel-negative blood units are transfused. Since the genetic background of Vel-negativity was discovered in 2013, DNA-based typing of the 17-bp deletion causing the phenotype has facilitated identification of Vel-negative blood donors. SMIM1, the gene underlying Vel, encodes a 78-amino acid erythroid transmembrane protein of unknown function. The transmembrane orientation of SMIM1 has been debated since experimental data supported both the N- and C-termini being extracellular. Likewise, computational predictions of its orientation were divided and potential alternatives such as monotopic or dual-topology have been discussed but not investigated.

We used a cell-free system to explore the topology of SMIM1 when synthesized in the endoplasmic reticulum (ER). SMIM1 was tagged with an opsin-derived N-glycosylation reporter at either the N- or C-terminus and synthesized in vitro using rabbit reticulocyte lysate supplemented with canine pancreatic microsomes as a source of ER membrane. SMIM1 topology was then determined by assessing the N-glycosylation of its N- or C-terminal tags. Complementary experiments were carried out by expressing the same SMIM1 variants in HEK293T/17 cells and establishing their membrane orientation by immunoblotting and flow cytometry.

Our data consistently indicate that SMIM1 has its short C-terminus located extracellularly and that it most likely belongs to the tail-anchored class of membrane proteins with the bulk of the polypeptide located in the cytoplasm. Having established its membrane orientation in an independent model system, future work can now focus on functional aspects of SMIM1 as a potential regulator of erythropoiesis.

SMIM1 missense mutations exert their effect on wild type Vel expression late in erythroid differentiation

BACKGROUND

Vel expression on erythrocytes is variable due to polymorphisms, complicating Vel typing. Weak Vel expression can be caused by mutations within SMIM1 in a heterozygous setting, suggesting a dominant negative effect of SMIM1 mutants on wild type (wt)SMIM1 expression. Here we report how SMIM1 expression is regulated during erythropoiesis, to understand its variable expression on erythrocytes.

STUDY DESIGN AND METHODS

Peripheral blood reticulocytes at different stages, cultured erythroid precursors and HEK293T cells were used to investigate expression and putative competition between wtSMIM1 and mutated SMIM1 VEL*01W.01, (c.152T>A (p.Met51Lys)), VEL*01W.02 (c.152T>G (p.Met51Arg)), and VEL*01W.03 (c.161T>C (p.Leu54Pro)).

RESULTS

Depending on the mutations in SMIM1 an effect on total and membrane expression of SMIM1 was observed in transfected HEK293T cells, but co-expression of wtSMIM1 and mutatedSMIM1 did not have an effect on wtSMIM1 membrane expression. During differentiation of donors expressing VEL*01W.01, VEL*01W.03, Vel-positive, Vel-negative (homozygote SMIM1*64_80del), and Vel-heterozygote SMIM1*64_80del primary human erythroblasts no overt defect was found in Vel expression dynamics or total SMIM1 expression levels when compared with wtSMIM1 erythroblasts. However, during enucleation, total Vel expression was significantly lower on reticulocytes of Vel-weak donors expressing heterozygote mutated SMIM1 compared to Vel-positive or Vel-heterozygote SMIM1*64_80del donors, while Vel expression on extruded nuclei was maintained. In addition, reticulocyte maturation in vivo showed further loss of Vel expression in these individuals and nearly absent on erythrocytes.

CONCLUSION

These results suggest that SMIM1 mutations exert a dominant negative effect on wtSMIM1 probably by affecting SMIM1 multimerization and thereby Vel epitope presentation at the latest stages of erythroid differentiation.

Targeted exome sequencing designed for blood group, platelet, and neutrophil antigen investigations: Proof-of-principle study for a customized single-test system

BACKGROUND

Immunohematology reference laboratories provide red blood cell (RBC), platelet (PLT), and neutrophil typing to resolve complex cases, using serology and commercial DNA tests that define clinically important antigens. Broad-range exome sequencing panels that include blood group targets provide accurate blood group antigen predictions beyond those defined by serology and commercial typing systems and identify rare and novel variants. The aim of this study was to design and assess a panel for targeted exome sequencing of RBC, PLT, and neutrophil antigen-associated genes to provide a comprehensive profile in a single test, excluding unrelated gene targets.

STUDY DESIGN AND METHODS

An overlapping probe panel was designed for the coding regions of 64 genes and loci involved in gene expression. Sequencing was performed on 34 RBC and 17 PLT/neutrophil reference samples. Variant call outputs were analyzed using software to predict star allele diplotypes. Results were compared with serology and previous sequence genotyping data.

RESULTS

Average coverage exceeded 250×, with more than 94% of targets at Q30 quality or greater. Increased coverage revealed a variant in the Scianna system that was previously undetected. The software correctly predicted allele diplotypes for 99.5% of RBC blood groups tested and 100% of PLT and HNA antigens excepting HNA-2. Optimal throughput was 12 to 14 samples per run.

CONCLUSION

This single-test system demonstrates high coverage and quality, allowing for the detection of previously overlooked variants and increased sample throughput. This system has the potential to integrate genomic testing across laboratories within hematologic reference settings.

Molecular Blood Group Screening in Donors from Arabian Countries and Iran Using High-Throughput MALDI-TOF Mass Spectrometry and PCR-SSP

Background and Aims

Only little is known about blood groups other than ABO blood groups and Rhesus factors in Arabian countries and Iran. During the last years, increased migration to Central Europe has put a focus on the question how to guarantee blood supply for patients from these countries, particularly because hemoglobinopathies with the need of regular blood support are more frequent in patients from that region. Therefore, blood group allele frequencies should be determined in individuals from Arabian countries and Iran by molecular typing and compared to a German rare donor panel.

Methods

1,111 samples including 800 individuals from Syria, 147 from Iran, 123 from the Arabian Peninsula, and 41 from Northern African countries were included in a MALDI-TOF MS assay to detect polymorphisms coding for Kk, Fy(a/b), Fy_null, C_w, Jk(a/b), Jo(a+/a-), Lu(a/b), Lu(8/14), Ss, Do(a/b), Co(a/b), In(a/b), Js(a/b), Kp(a/b), and variant alleles RHCE*c.697C>G and RHCE *c.733C>G. Yt(a/b), S-s-U-, Vel_null, Co_null, and RHCE *c.667G>T were tested by PCR-SSP.

Results

Of the Arabian donors, 2% were homozygous for the FY *02.01N allele (Fy_null), and 15.7% carried the heterozygous mutation. However, 0.8% of the German donors also carried 1 copy of the allele. 3.6% of all and 29.3% of Northern African donors were heterozygous for the RHCE *c.733C>G substitution, 0.4% of the Syrian probands were heterozygous for DO *01/DO *01.-05, a genotype that was lacking in German donors. Whereas the KEL *02.06 allele coding for the Js(a) phenotype was missing in Germans; 0.8% of the Syrian donors carried 1 copy of this allele. 1.8% of the Syrian but only 0.3% of the German donors were negative for YT *01. One donor from Northern Africa homo-zygously carried the GYPB *270+5g>t mutation, inducing the S-s-U+^w phenotype, and in 2 German donors a GYPB *c.161G>A exchange, which induces the Mit+ phenotype, caused a GYPB *03 allele dropout in the MALDI assay. The overall failure rate of the Arabian panel was 0.4%.

Conclusions

Some blood group alleles that are largely lacking in Europeans but had been described in African individuals are present in Arabian populations at a somewhat lower frequency. In single cases, it could be challenging to provide immunized Arabian patients with compatible blood.

Dimerization of small integral membrane protein 1 promotes cell surface presentation of the Vel blood group epitope

The Vel blood group antigen is carried on the short extracellular segment of the 78-amino-acid-long, type II transmembrane protein SMIM1 of unknown function. Here, using biochemical analysis and flow cytometry of cells expressing wild-type and mutant alleles of SMIM1, we demonstrate that dimerization of SMIM1 promotes cell surface display of the Vel epitope. We show that SMIM1 dimerization is mediated both by an extracellular Cys77-dependent, homomeric disulfide linkage and via a GxxxG helix-helix interaction motif in the transmembrane domain. These results provide important context for the observed variability in reactivity patterns of clinically important anti-Vel identified in patient sera.

SMIM1 at a glance; discovery, genetic basis, recent progress and perspectives

Recent elucidation of the genetic basis of the Vel blood group system has offered the field of blood transfusion medicine an additional consideration in determining the causes of hemolytic reactions after a patient is transfused. The identification of the SMIM1 gene to be responsible for the Vel blood group allows molecular based tools to be developed to further dissect the function of this antigen. Genetic signatures such as the homozygous 17 bp deletion and the heterozygous 17 bp deletion in combination with other single nucleotide polymorphisms (SNPs) and insertion sequences regulate the expression level of the gene. With this knowledge, it is now possible to study this antigen in-depth.

SMIM1 intron 2 gene variations leading to variability in Vel antigen expression among Brazilian blood donors

BACKGROUND

There is a significant inter-individual heterogeneity of Vel antigen expression which can lead to inaccuracies on Vel phenotyping of blood donors and, potentially, to hemolytic post-transfusion reactions. Our aim was to evaluate the impact of genetic variants in the SMIM1 intron 2 on the expression of Vel antigen among Brazilian blood donors harboring the c.64_80del17 deletion in heterozygosity.

METHODS

Donors presenting the SMIM1 c.64_80del17 in heterozygosity were included in the study and subjected to SMIM1 intron 2 direct sequencing aiming to genotype the following polymorphisms: rs143702418, rs1181893, rs191041962, rs6673829, rs1175550 and rs9424296.

RESULTS

SMIM1 intron 2 sequencing was performed on two hundred donors presenting one c.64_80del17 allele. The rs1175550 polymorphism significantly impacted on Vel antigen expression. Variations in the strength of agglutination on Vel phenotyping were also observed according to the rs6673829 genotype, but this difference did not persist with statistical relevance after multivariate analysis.

CONCLUSION

The presence of the rs1175550A allele of SMIM1 is significantly and independently associated with a decrease in Vel antigen expression. Even though the population in Brazil is intensely mixed, the allele frequencies obtained in the current study were very similar to that reported for Europeans.

Development of a recombinant anti-Vel immunoglobulin M to identify Vel-negative donors

BACKGROUND

Alloimmunization against the high-frequency Vel blood group antigen may result in transfusion reactions or hemolytic disease of fetus and newborn. Patients with anti-Vel alloantibodies require Vel-negative blood but Vel-negative individuals are rare (1:4000). Identification of Vel-negative donors ensures availability of Vel-negative blood; however, accurate Vel blood group typing is difficult due to variable Vel antigen expression and limited availability of anti-Vel typing sera. We report the production of a recombinant anti-Vel that also identifies weak Vel expression.

STUDY DESIGN AND METHODS

A recombinant anti-Vel monoclonal antibody was produced by cloning the variable regions from an anti-Vel-specific B cell isolated from an alloimmunized patient into a vector harboring the constant regions of immunoglobulin (Ig)G1-kappa or IgM-kappa. Antibody Vel specificity was tested by reactivity to SMIM1-transfected HEK293T cells and by testing various red blood cells (RBCs) of donors with normal, weak, or no Vel expression. High-throughput donor screening applicability was tested using an automated blood group analyzer.

RESULTS

A Vel-specific IgM class antibody was produced. The antibody was able to distinguish between Vel-negative and very weak Vel antigen-expressing RBCs by direct agglutination and in high-throughput settings using a fully automated blood group analyzer and performed better than currently used human anti-Vel sera. High-throughput screening of 13,288 blood donations identified three new Vel-negative donors.

CONCLUSION

We generated a directly agglutinating recombinant anti-Vel IgM, M3F5S-IgM, functional in manual, automated agglutination assays and flow cytometry settings. This IgM anti-Vel will improve diagnostics by facilitating the identification of Vel-negative blood donors.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry analysis of 36 blood group alleles among 396 Thai samples reveals region-specific variants

BACKGROUND

Blood group phenotype variation has been attributed to potential resistance to pathogen invasion. Variation was mapped in blood donors from Lampang (northern region) and Saraburi (central region), Thailand, where malaria is endemic. The previously unknown blood group allele profiles were characterized and the data were correlated with phenotypes. The high incidence of the Vel-negative phenotype previously reported in Thais was investigated.

STUDY DESIGN AND METHODS

DNA from 396 blood donors was analyzed by matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. Outliers were investigated by serology and DNA sequencing. Allele discrimination assays for SMIM1 rs1175550A/G and ACKR1 rs118062001C/T were performed and correlated with antigen expression.

RESULTS

All samples were phenotyped for Rh, MNS, and K. Genotyping/phenotyping for RhD, K, and S/s showed 100% concordance. Investigation of three RHCE outliers revealed an e-variant antigen encoded by RHCE*02.22. Screening for rs147357308 (RHCE c.667T) revealed a frequency of 3.3%. MN typing discrepancies in 41 samples revealed glycophorin variants, of which 40 of 41 were due to Mi^a . Nine samples (2.3%) were heterozygous for FY*01W.01 (c.265C > T), and six samples (1.5%) were heterozygous for JK*02N.01. All samples were wildtype SMIM1 homozygotes with 97% homozygosity for rs1175550A.

CONCLUSIONS

Matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry is an efficient method for rapid routine genotyping and investigation of outliers identified novel variation among our samples. The expected high prevalence of the Mi(a+) phenotype was observed from both regions. Of potential clinical relevance in a region where transfusion-dependent thalassemia is common, we identified two RHCE*02 alleles known to encode an e-variant antigen.

The use of next-generation sequencing for the determination of rare blood group genotypes

OBJECTIVES

Next-generation sequencing (NGS) for the determination of rare blood group genotypes was tested in 72 individuals from different ethnicities.

BACKGROUND

Traditional serological-based antigen detection methods, as well as genotyping based on specific single nucleotide polymorphisms (SNPs) or single nucleotide variants (SNVs), are limited to detecting only a limited number of known antigens or alleles. NGS methods do not have this limitation.

METHODS

NGS using Ion torrent Personal Genome Machine (PGM) was performed with a customised Ampliseq panel targeting 15 different blood group systems on 72 blood donors of various ethnicities (Caucasian, Hispanic, Asian, Middle Eastern and Black).

RESULTS

Blood group genotypes for 70 of 72 samples could be obtained for 15 blood group systems in one step using the NGS assay and, for common SNPs, are consistent with previously determined genotypes using commercial SNP assays. However, particularly for the Kidd, Duffy and Lutheran blood group systems, several SNVs were detected by the NGS assay that revealed additional coding information compared to other methods. Furthermore, the NGS assay allowed for the detection of genotypes related to VEL, Knops, Gerbich, Globoside, P1PK and Landsteiner-Wiener blood group systems.

CONCLUSIONS

The NGS assay enables a comprehensive genotype analysis of many blood group systems and is capable of detecting common and rare alleles, including alleles not currently detected by commercial assays.