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.

Disruption of the tumour-associated EMP3 enhances erythroid proliferation and causes the MAM-negative phenotype.

The clinically important MAM blood group antigen is present on haematopoietic cells of all humans except rare MAM-negative individuals. Its molecular basis is unknown. By whole-exome sequencing we identify EMP3, encoding epithelial membrane protein 3 (EMP3), as a candidate gene, then demonstrate inactivating mutations in ten known MAM-negative individuals. We show that EMP3, a purported tumour suppressor in various solid tumours, is expressed in erythroid cells. Disruption of EMP3 by CRISPR/Cas9 gene editing in an immortalised human erythroid cell line (BEL-A2) abolishes MAM expression. We find EMP3 to associate with, and stabilise, CD44 in the plasma membrane. Furthermore, cultured erythroid progenitor cells from MAM-negative individuals show markedly increased proliferation and higher reticulocyte yields, suggesting an important regulatory role for EMP3 in erythropoiesis and control of cell production. Our data establish MAM as a new blood group system and demonstrate an interaction of EMP3 with the cell surface signalling molecule CD44.

Missense mutations in the C-terminal portion of the B4GALNT2-encoded glycosyltransferase underlying the Sd(a-) phenotype.

Sda is a high-frequency carbohydrate histo-blood group antigen, GalNAcß1-4(NeuAcα2-3)Galß, implicated in pathogen invasion, cancer, xenotransplantation and transfusion medicine. Complete lack of this glycan epitope results in the Sd(a-) phenotype observed in 4% of individuals who may produce anti-Sda. A candidate gene (B4GALNT2), encoding a Sda-synthesizing ß-1,4-N-acetylgalactosaminyltransferase (ß4GalNAc-T2), was cloned in 2003 but the genetic basis of human Sda deficiency was never elucidated. Experimental and bioinformatic approaches were used to identify and characterize B4GALNT2 variants in nine Sd(a-) individuals. Homozygosity for rs7224888:T > C dominated the cohort (n = 6) and causes p.Cys466Arg, which targets a highly conserved residue located in the enzymatically active domain and is judged deleterious to ß4GalNAc-T2. Its allele frequency was 0.10-0.12 in different cohorts. A Sd(a-) compound heterozygote combined rs7224888:T > C with a splice-site mutation, rs72835417:G > A, predicted to alter splicing and occurred at a frequency of 0.11-0.12. Another compound heterozygote had two rare nonsynonymous variants, rs148441237:A > G (p.Gln436Arg) and rs61743617:C > T (p.Arg523Trp), in trans. One sample displayed no differences compared to Sd(a+). When investigating linkage disequilibrium between B4GALNT2 variants, we noted a 32-kb block spanning intron 9 to the intergenic region downstream of B4GALNT2. This block includes RP11-708H21.4, a long non-coding RNA recently reported to promote tumorigenesis and poor prognosis in colon cancer. The expression patterns of B4GALNT2 and RP11-708H21.4 correlated extremely well in >1000 cancer cell lines. In summary, we identified a connection between variants of the cancer-associated B4GALNT2 gene and Sda, thereby establishing a new blood group system and opening up for the possibility to predict Sd(a+) and Sd(a‒) phenotypes by genotyping.

GBGT1 is allelically diverse but dispensable in humans and naturally occurring anti-FORS1 shows an ABO-restricted pattern.

BACKGROUND: The FORS histo-blood group system was described in 2013 and much remains to be investigated regarding its genetic and immunohematologic characteristics, as well as its clinical importance. While presence of the c.887G>A-mutated GBGT1 gene, which results in FORS1 glycosphingolipid expression on human red blood cells (RBCs), is rare in the populations tested so far, naturally occurring anti-FORS1 in plasma appears common. STUDY DESIGN AND METHODS: The Erythrogene database was utilized to probe genetic variation in GBGT1 among 2504 individuals in the 1000 Genomes Project. We screened 1108 Swedish blood donors for three principally important single-nucleotide polymorphisms (c.363C>A, c.886C>T, and c.887G>A) and selected samples were analyzed further. Screening for naturally occurring anti-FORS1 in plasma from 100 donors was performed using antigen-positive RBCs. RESULTS: We identified 68 GBGT1 alleles, of which three were previously listed blood group alleles. Eight potential null alleles were observed, based on three different nonsense mutations. Four healthy donors were found homozygous for c.363C>A, which truncates the GBGT1-encoded Fs synthase prematurely. This is the first description of human knock-outs for GBGT1. The c.886C>T mutation that alters the same codon (p.Arg296Trp) changed by c.887G>A (p.Arg296Gln) was overexpressed to investigate if it induces the FORS1+ phenotype. However, c.886C>T did not result in synthesis of FORS1. We detected anti-FORS1 in 10% of all donors tested but none in the A1 or A1B groups. CONCLUSION: We have extended the knowledge of GBGT1 variants, allele frequencies, and the characteristics of naturally occurring antibodies in our newest carbohydrate blood group system, FORS. The finding of c.363C>A-homozygous donors indicates that GBGT1 is dispensable.

Disruption of a GATA1-binding motif upstream of XG/PBDX abolishes Xga expression and resolves the Xg blood group system.

The Xga blood group is differentially expressed on erythrocytes from men and women. The underlying gene, PBDX, was identified in 1994, but the molecular background for Xga expression remains undefined. This gene, now designated XG, partly resides in pseudoautosomal region 1 and encodes a protein of unknown function from the X chromosome. By comparing calculated Xga allele frequencies in different populations with 2612 genetic variants in the XG region, rs311103 showed the strongest correlation to the expected distribution. The same single-nucleotide polymorphism (SNP) had the most significant impact on XG transcript levels in whole blood (P = 2.0 × 10-22). The minor allele, rs311103C, disrupts a GATA-binding motif 3.7 kb upstream of the transcription start point. This silences erythroid XG messenger RNA expression and causes the Xg(a-) phenotype, a finding corroborated by SNP genotyping in 158 blood donors. Binding of GATA1 to biotinylated oligonucleotide probes with rs311103G but not rs311103C was observed by electrophoretic mobility shift assay and proven by mass spectrometry. Finally, a luciferase reporter assay indicated this GATA motif to be active for rs311103G but not rs311103C in HEL cells. By using an integrated bioinformatic and molecular biological approach, we elucidated the underlying genetic basis for the last unresolved blood group system and made Xga genotyping possible.

Identification of human glycosyltransferase genes expressed in erythroid cells predicts potential carbohydrate blood group loci.

Glycans are biologically important structures synthesised by glycosyltransferase (GT) enzymes. Disruptive genetic null variants in GT genes can lead to serious illness but benign phenotypes are also seen, including antigenic differences on the red blood cell (RBC) surface, giving rise to blood groups. To characterise known and potential carbohydrate blood group antigens without a known underlying gene, we searched public databases for human GT loci and investigated their variation in the 1000 Genomes Project (1000 G). We found 244 GT genes, distributed over 44 families. All but four GT genes had missense variants or other variants predicted to alter the amino acid sequence, and 149 GT genes (61%) had variants expected to cause null alleles, often associated with antigen-negative blood group phenotypes. In RNA-Seq data generated from erythroid cells, 155 GT genes were expressed at a transcript level comparable to, or higher than, known carbohydrate blood group loci. Filtering for GT genes predicted to cause a benign phenotype, a set of 30 genes remained, 16 of which had variants in 1000 G expected to result in null alleles. Our results identify potential blood group loci and could serve as a basis for characterisation of the genetic background underlying carbohydrate RBC antigens.

Allele-selective RUNX1 binding regulates P1 blood group status by transcriptional control of A4GALT.

P1 and Pk are glycosphingolipid antigens synthesized by the A4GALT-encoded α1,4-galactosyltransferase, using paragloboside and lactosylceramide as acceptor substrates, respectively. In addition to the compatibility aspects of these histo-blood group molecules, both constitute receptors for multiple microbes and toxins. Presence or absence of P1 antigen on erythrocytes determines the common P1 (P1+Pk+) and P2 (P1-Pk+weak) phenotypes. A4GALT transcript levels are higher in P1 individuals and single-nucleotide polymorphisms (SNPs) in noncoding regions of A4GALT, particularly rs5751348, correlate with P1/P2 status. Despite these recent findings, the molecular mechanism underlying these phenotypes remains elusive. The In(Lu) phenotype is caused by Krüppel-like factor 1 (KLF1) haploinsufficiency and shows decreased P1 levels on erythrocytes. We therefore hypothesized KLF1 regulates A4GALT expression. Intriguingly, P1 -specific sequences including rs5751348 revealed potential binding sites for several hematopoietic transcription factors, including KLF1. However, KLF1 binding did not explain P1 -specific shifts in electrophoretic mobility-shift assays and small interfering RNA silencing of KLF1 did not affect A4GALT transcript levels. Instead, protein pull-down experiments using P1 but not P2 oligonucleotide probes identified runt-related transcription factor 1 (RUNX1) by mass spectrometry. Furthermore, RUNX1 binds P1 alleles selectively, and knockdown of RUNX1 significantly decreased A4GALT transcription. These data indicate that RUNX1 regulates A4GALT and thereby the expression of clinically important glycosphingolipids implicated in blood group incompatibility and host-pathogen interactions.

Erythrogene: a database for in-depth analysis of the extensive variation in 36 blood group systems in the 1000 Genomes Project.

Blood group genotyping has recently developed into a clinical tool to improve compatibility of blood transfusions and management of pregnancies. Next-generation sequencing (NGS) is rapidly moving toward routine practice for patient and donor typing and has the potential to remedy some of the limitations of currently used platforms. However, a large-scale investigation into the blood group genotypes obtained by NGS in a multiethnic cohort is lacking. The 1000 Genomes Project provides information on genome variation among 2504 individuals representing 26 populations worldwide. We extracted their NGS data for all 36 blood group systems to a custom-designed database. In total, 210 412 alleles from 43 blood group-related genes were imported and curated. Matching algorithms were developed to compare them to blood group variants identified to date. Of the 1241 non-synonymous variants identified in the coding regions, 241 are known blood group polymorphisms. Interestingly, 357 of the remaining 1000 variants are predicted to occur on extracellular portions of 31 different blood group-carrying proteins and some may represent undiscovered antigens. Of the alleles analyzed, 1504 were not previously described. The ABO/GBGT1/FUT2/FUT3 and GYPB/GYPC genes showed the highest degree of variation per kilobase coding sequence, and ACKR1 variants had the most skewed distribution across 5 continental superpopulations in the dataset. Results were exported to an online search engine, www.erythrogene.com, which presents data according to the allele nomenclature developed for clinical reporting by the International Society of Blood Transfusion. The established database deepens our knowledge on blood group polymorphism globally and provides a long-sought platform for future research.