A comprehensive characterization of blood group antigen variants in the Middle Eastern population genomes - Insights into genetic epidemiology

Population scale genomic datasets

Whole genome and whole exome sequencing datasets generated by population specific studies in the Greater Middle Eastern region were compiled and used for the analysis. This includes data from the Greater Middle East Variome (21), The Qatar genome and exomes (10) and the Iranome (9) genome projects. These accounted for a total of 2828 high quality exomes and 88 genomes (GME : 1111 exomes; Qatar : 88 genomes and 917 exomes ; Iranome : 800 exomes), encompassing 689299, 20159694 and 1575702 total genetic variants respectively. The population scale data is a collation of major ethnic groups of the middle eastern region including North Western Africa (NWA), North East Africa (NEA), Arabian Peninsula (AP), Turkish Peninsula (TP), Persia and Pakistan (PP), Syrian Desert (SD), Qatari subgroups and Iranian Arabs, Azeris, Baluchs, Kurds, Lurs, Turkmen and Persians. All the datasets correspond to human genome 19 (GRCh37/hg19) assembly. Genetic variations were retrieved from the variant call format (VCF) files and were used for further analyses.

Human blood group related alleles and variations

Clinically approved genomic coordinates of 50 genes associated with 43 human blood group systems and 2 erythroid specific transcription factors were fetched from Locus Genomic Reference (LRG) (22) and were used for analysis as described previously (23). There exists a handful of public resources that provide comprehensive collection of human blood alleles and antigens including the International Society of Blood Transfusion (ISBT) (24), The Blood Group Antigen Gene Mutation Database (dbRBC) (25), Erythrogene (26) RhesusBase (27), Blood Antigens (28,29) and BGvar (30). A comprehensive list of approved and predicted human blood group alleles were retrieved in a pre-formatted template (30) and was used as the reference data for further allele matching procedures.

Data processing and variant annotations

Primary analysis involved the filtering of all genetic variations spanning the LRG approved genomic coordinates of blood group genes and erythroid specific transcription factors. Subsequently, the filtered variants were systematically annotated for their functional consequences from a range of computational tools including SIFT (31), Polyphen (32), LRT, MutationTaster, Mutation Assessor (33), FATHMM (34), PROVEAN (35), CADD (36), GERP (37), PhyloP (38) and PhastCons using ANNOVAR (v. 2018-04-06) (39).

Known and novel blood group alleles/variants identification

Secondary stages of variant analyses involved the filtering and identification of well known and novel blood group alleles. Reference dataset was used to match all the variants which already had a known associated blood phenotype and reference genome nomenclature was conferred to those blood group systems with no blood group associated variants. In addition, potentially novel and rare blood group variants were also filtered and analyzed. A variant was deemed potentially novel if it lacked previous literature reports and SNP identification number (dbSNP ID). Variants with minor allele frequencies (MAF) < 5% were also filtered and checked for their association with blood group related phenotypes. A complete schematic representation of the workflow is shown in Figure 1.

• Download figure
• Open in new tab

Figure 1.

Schematic representation of the methodology followed for data analysis

Prediction of blood phenotypes in Qatari subpopulations

The publicly available genomic data of 1005 Qatari individuals comprising 5 major subpopulations were analyzed in the study to decipher the pattern of distribution of blood phenotypes among the Qatari subpopulations. The variant call file was systematically subsetted into each subpopulation using VCFTOOLS (40). Bespoke scripts were used to generate the number of samples with homozygous and heterozygous genotypes. Similarities and distinct differences in blood alleles and phenotype frequencies were studied. Minor allele frequencies were compared among the subpopulations and statistical significance was observed using Fisher’s exact test with a p-value < 0.05.

Estimation and systematic comparison of blood allele frequencies

Blood group associated variants filtered from all the datasets were systematically compiled in variant call or annotated formats along with their genotype information. Allele frequencies were estimated using PLINK {ref}. In addition, allele frequencies of the variants were fetched from major global population datasets including 1000 Genomes project {ref}, Exome Aggregation Consortium (ExAC v.0.3) {ref} and Genome Aggregation Database (gnomAD) {ref} and were used for comparison. Frequency of the filtered alleles among the subpopulations were also calculated.

Identification of clinically significant and actionable blood group alleles

Variants with known blood group associated phenotypes filtered from the datasets were further checked for their clinical relevance in transfusion procedures and in pregnancy settings. Literature databases and public resources like Pubmed and Google scholar were systematically queried to retrieve relevant clinically significant evidence of alleles.

Compendium of Greater Middle East blood group variants - Database architecture

The filtered variants were compiled into a preformatted mastersheet and annotations were duly mapped. The variant data and corresponding annotations were transformed to JavaScript Object Notation format and ported onto MongoDB 3.4.1. The user-friendly web interface for querying the database was coded in PHP 7.0, AngularJS, HTML, Bootstrap 4 and CSS. The web server was configured in Apache HTTP server.

Overview of blood group alleles in Greater Middle East

We retrieved a total 18717 unique variants mapping to 50 human blood group genes, out of which 1697, 15140 and 3496 variants were from the Greater Middle East Variome, The Qatar genome and the Iranome datasets respectively. Supplementary Table 1 provides a comprehensive summary of all the variant counts in blood group genes fetched from the above mentioned datasets. Of the variants, 2428 and 15 variants were found across exonic and splicing regions respectively. Of these, there were 1505 non synonymous SNVs, 766 synonymous SNVs, 50 stop gain and 3 stop loss variants. Variant summary and functional classifications across the datasets are schematically represented in Figure 2.

• Download figure
• Open in new tab

Figure 2.

Overview of the functional classification of blood group related variants in various datasets used in the study

Characterization of genetic variations in human blood group genes

Of the variants which mapped to the 50 blood group genes, a total of 147, 110 and 112 variants were previously shown to be associated with blood group related phenotypes in GME, Qatar and Iranome datasets respectively. Of these variants, a total of 77 unique variants belonging to 20 blood group systems were overlapping in all three datasets used in this study. The details of the variants along with their frequencies are tabulated in Table 1. Figure 3 illustrates the distribution of these variants and their corresponding allele frequencies in various datasets. Forty-six (46) of the total 77 variants were SNVs mapping back to 16 blood groups and the rest were combination mutations for which blood phenotypes could not be predicted. Phenotypes were reported the same as that of the reference genome (hg19) nomenclature for the rest blood group systems as previously described (41), (23). Table 2 summarizes the phenotypes of the variants and their corresponding zygosity information across the datasets. Phenotype frequencies of blood group systems with known blood group associated variants and those with reference genome nomenclature are detailed in Supplementary Table 2.

• Download figure
• Open in new tab

Figure 3.

Schematic illustration of the distribution of common variants with blood group associated phenotypes along with their corresponding allele frequencies. (A), Heatmap illustrating the similarities and differences in the allele frequencies among various datasets (B), Distribution of blood group related variants in various datasets used in the study (C) Distribution of common alleles among human blood group systems

Overview of potentially novel and rare blood group variants

Variants were systematically classified based on the dbSNP identifiers and were deemed as potentially novel if they lacked one. This analysis revealed a total of 1727 potentially novel unique variants across the datasets, of which 730 were found to span the exonic and splicing regions. Seventy (70) exonic variants which were not reported previously in global population frequencies including 1000 Genomes, gnomAD and EXAC were predicted to be deleterious by at least three or more computational tools. Distribution of these novel variants across different blood group genes along with their corresponding frequencies in the Middle Eastern datasets are summarized in Supplementary Table 3. A missense H blood group variant (p. Arg220Leu, chr19:49253880C>A) which was computationally predicted as deleterious, has a couple of alternate variants reported at the same position associated with weak H+ phenotype (p.Arg220Cys - FUT*01W.09 and p.Arg220His - FUT*01W.10). Yet another missense variant predicted from the same blood group (p. Try154Asp, chr19:49254079A>C) possess reported variants at the same amino acid position associated with weak and null H phenotypes (p.Tyr154Cys - FUT*01N.02, p.Tyr154His - FUT*01W.o5.01 and p.Tyr154Ter - FUT*01N.03).

Two missense variants belonging to JUNIOR blood group predicted as deleterious (p.Gly262Glu, chr4:89039317C>T and p.Arg246Leu, chr4:89039365C>A) were found to occur at the same amino acid sites which were reported to be associated with Jr(a-) phenotypes (p.Gly262X - ABCG2*01N.05 and p.Arg246X - ABCG2*01N.03 respectively).

Distribution of null and partial/weak antigens

Weak or partial expression of RBC antigens carries increased risk of hemolytic transfusion reactions as they are often mistyped as antigen-negative. Of the total variants analysed, our study reports 19 variants responsible for weak/partial or null blood group phenotypes. The RHD*weak D type 4.0 allele (rs1053355) was observed in 10 individuals in a homozygous state. This allele, responsible for the weak expression of D antigen, was found to be often mistyped as D- in serological tests (42). In addition, one of the rare low frequency antigens of the RH blood group system, C^w+ (Rh8 - rs138268848) was found in about 0.17% of Middle Eastern populations (0.13% in Iranome and 0.36% in GME). This antigen has been reported to have serious transfusion implications including hemolytic disease of the newborn when mismatched (43).

A weak JK*A allele (JK*01W.01 - rs2298720) of the KIDD blood group system, associated with weakened expression of Jk^a is found capable of causing alloimmunization when mistyped serologically as Jk_null (44), (45). Distribution of this variant is found to vary highly among various global populations with an allele frequency as high as 40.6% in East Asians, 26.8% in South Asians, 18-19% in Latino/Admixed Americans and African/African Americans, 12% in Europeans of Finnish origin, 10% in Ashkenazi Jews to as low as 7% in Middle Easterners, 6% in Europeans of non-Finnish origins and 3% in Amish population. Supplementary Table 4 enlists the complete details of weak/partial alleles identified in the study dataset.

Clinically significant blood group variants

The KELL blood group system is significant in the field of transfusion medicine owing to the severe transfusion incompatible reactions induced by its antigens (46). This system comprises over 20 different antigens, of which K1 (K) has been well studied as a strong immunogen capable of causing hemolytic disease of the newborn in sensitized mothers. Approximately 9% of the random RBC donor samples were found to react against K1 (46,47). Anti-K1 is ranked third most clinically significant antibody after ABO and RH(D). HDN induced by anti-K is characterized by immune destruction of K+ erythroid (47) progenitor cells by macrophages in the fetal liver (48). Frequency of this variant in Middle Easterners was found to be 3.5% which was comparable to Europeans of non-Finnish origin (4.1%) and differ from South Asians (<1%) (49), (50)

Anti-Cr^a is a rarely encountered antibody against the Cr(a) antigen of the CROMER blood group system. CROM:c.679G>C (rs60822372) variation in exon 6 of CROM gene is responsible for Cr(a-) phenotype which is characterized by the absence of Cr(a) antigen on the RBC membrane. Studies have shown negative crossmatching of individuals carrying the variant in homozygous state and positive while in heterozygous state (51). This variant was observed at a frequency of about 10% in Middle Eastern populations which was found to be comparable to African/African Americans and South Asians (52), (53),(54).

Blood group phenotyping of Qatari subpopulations

The genomic data of 1005 Qataris comprises 5 major subpopulations namely African (n=70), Arabian (n=193), Bedouin (n=490), Persian (n=170) and South Asian (n=76). Distribution of blood group related variants in different subpopulations was identified and corresponding phenotypes were predicted . In addition, systematic comparisons of allele and phenotype frequencies among the subpopulations were performed.

Comparison and statistical analysis of blood group variants/profiles among Middle Eastern population datasets

Distinct differences in the minor allele frequencies of blood group alleles between the Middle Easterners and various other global populations were statistically estimated using Fisher’s exact test. A total of 26 variants mapping back to 12 blood group systems were found statistically distinct in all the Middle Easterners (Qatar, GME and Iranome) in comparison to global datasets (1000 Genomes Project and gnomAD version?). Summary of these variants along with their p-value is provided in Supplementary Table 5a. In addition, differences in the pattern of distribution of these blood group variants were checked among the 5 major Qatari subpopulations. The Arabian subpopulation was found to carry the maximum number of variants whose frequencies were statistically distinct from the overall frequency in the Qatar dataset. Complete details of these variants are tabulated in Supplementary Table 5b.

User friendly data access

alnasab (Alleles and antigens in Arab and Persian populations associated with blood groups) is a user-friendly online search engine, composed of a comprehensive collection of blood group related variations mined from various population scale datasets of the Greater Middle Eastern region. The resource encompasses a total of 2149 blood group related variants fetched from about 47 blood group genes. The search interface enables the user to query the database based on gene name, variant name (chr-pos:ref>alt), nucleotide change and protein change. Variant information is provided in six different sections namely, (i) Blood group system, (ii) Variant details, (iii) Transcript details, (iv) Variant annotation, (v) Global allele frequencies and (vi) Allele frequencies - Greater Middle East. Figure 4a and 4b illustrate the search features and various sections of result display of the resource.

• Download figure
• Open in new tab

Figure 4.

Search query and sections of result display in alnasab. (A) Homepage of the web resource with query search examples. (B) Result display which provides an extensive summary of variant functional annotations in various sections

Clinical impact

An extended screening of blood group antigens between the donors and the recipients other than the major ABO, RH groups and minor DUFFY, KELL and KIDD groups is often triggered only after encountering clinical complications. Frequencies of a total of 26 variants belonging to 12 different blood group systems were found to differ significantly from other global populations. Fetal anaemia and Hemolytic disease of the fetus and newborn (HDFN) are the most serious consequences of transfusion in the majority of pregnancy settings (55). The clinically significant antibodies involved in such cases include Rh, C, E, c and e antigens with greater than 50% risk of mild to severe HDFN (56),(57). Our study reports a weak RHD allele (RHD*weak D type 4.0) with an overall frequency of 0.01% in the Middle eastern population. Often patients carrying this allele get serologically categorized to be transfused with D-RBCs. A recent report from a transfusion work group recommendation explained the necessity of RHD molecular genotyping in expectant mothers which would help categorize RHD alleles that can further allow safe D+ transfusion (42).

Literature comparison

There exists a handful of research studies and reports documenting the blood group profiles of donors across various parts of the Arabian peninsula. Phenotypes of ABO, RH and various other minor blood group systems predicted by standard serological methods were found to match with the genotyping profiling results of the study. (50), (58), (58,59), (60), (61). Js^a (KEL6) is a low-incidence antigen reported in less than 1% of the general global populations, whereas on other hand, could be prevalent upto 19.5% in African Americans (62). Js(a+b-) phenotype which was not previously reported in serological studies was found in one sample from North East Africa. In compliance with earlier evidence, frequencies of minor blood group phenotypes including Fy(a+b+), K-k+, Lu(a-b+), Jk(a+b+), Do(a+b+), Co(a+) were found comparable.

Limitations and future needs

Owing to the limitations in the available public datasets, our study is mainly limited by the fact that RBC antigenic expressions regulated by large deletions and insertions (especially in RH and MNS blood groups) have not been profiled. In addition, the level of concordance with serology based phenotype predictions to investigate the functional implications of novel and rare variants remains to be explored and validated in future studies.