Using Whole Genome Sequencing to Characterize Clinically Significant Blood Groups Among Healthy Older Australians
=================================================================================================================

* Sudhir Jadhao
* Candice Davison
* Eileen V. Roulis
* Simon Lee
* Paul Lacaze
* Moeen Riaz
* John J McNeil
* David M Thomas
* Natalie M. Pecheniuk
* Catherine A. Hyland
* Robert L. Flower
* Shivashankar H. Nagaraj

## ABSTRACT

There have been no comprehensive studies of a full range of blood group polymorphisms within the Australian population. The problem is compounded by the absence of any databases carrying genomic information on chronically transfused patients and low frequency blood group antigens in Australia. Here, we use RBCeq, a web server-based blood group genotyping software, to identify unique blood group variants among Australians and compare the variation detected versus global data. Whole genome sequencing data was analysed from for 2796 healthy older Australians from the Medical Genome Reference Bank and compared with data from 1000G phase 3 (1KGP3) databases comprising 661 African, 347 American, 503 European, 504 East Asian, and 489 South Asian participants. There were 688 rare variants detected in this Australian sample population, including nine variants that had clinical associations. Notably, we identified 149 variants that were computationally predicted to be novel and deleterious. No clinically significant rare or novel variants were found associated with the genetically complex ABO blood group system. For the Rh blood group system, one novel and 16 rare variants were found. Our detailed blood group profiling results provide a starting point for the creation of an Australian blood group variant database.

**Key points**

*   We identified unique blood group variants among the healthy older Australian population compared with global data using RBCeq software.

*   Our detailed blood group profiling result may be a starting point for the creation of an Australian blood group variant database.

## INTRODUCTION

Modern transfusion medicine has highlighted the geographical and ethnic variability in blood group allele and genotype frequencies. The International Society of Blood Transfusion (ISBT) has recognised 43 red blood cell RBC membrane-associated blood group systems, involving over 360 antigen phenotypes defined by >1500 alleles across 45 genes and two transcription factors (TF).1-4 Blood group genes encode multiple structurally and functionally distinct molecules and exhibit varying degrees of polymorphic complexity (insertions/deletions [indels], single nucleotide variants [SNVs], copy number variations [CNVs], and structural variants [SVs]) within the population.5,6 Misidentification of any of these variants may contribute towards alloimmunizationof individuals requiring RBC transfusion support, which introduces a risk of adverse events such as Hemolytic Transfusion Reactions (HTRs), Hemolytic Disease of the Fetus and Newborn (HDFN), pregnancy complications, or more subtle allergic reactions of clinical significance 7.

Routinely in Australia, the Australian New Zealand Society of Blood Transfusion (ANZSBT) guidelines for Transfusion and Immunohematology laboratory practice advise to select ABO and RhD matched red cell products for transfusion to ensure transfusion safety.8,9,10 The RH (e.g., D, C/c, E/e) and MNS (e.g., M/N, S/s, U) blood group antigens are encoded by multiple complex alleles, genetic variations, and gene rearrangements between *RHD/RHCE* and *GYPA/GYPB/GYPE* with distinct population specific distributions. For example, the D negative blood type is prevalent in groups with Caucasian (15%)11 and African (8%) ancestries but rare amongst Asians (<0.1%)12. The variations in blood group antigen expression also may increase or decrease host susceptibility to infections in some diseases. The glycoprotein that carries the antigens of the Duffy blood group system (FY) is a typical chemokine receptor 1 (ACKR1), also previously known as Duffy Antigen Receptor for Chemokines (DARC). ACKR1 acts as a receptor for *Plasmodium vivax* and *Plasmodium knowlesi* malaria parasites. The absence of FY glycoprotein genotype is common among Africans and African-Americans, conferring these populations with resistance to malaria infection.13,14 Variations in *ACKR1* are also associated with a survival advantage in leukopenic HIV patients.15 The recessive African-specific *ACKR1* null allele increases the risk of HIV-1 infection.14 In the ABO blood group system, individuals with group O blood type are more prone to vaso-occlusive crises associated with sickle cell disease (SCD).16 However, the sickle hemoglobin (HbS) of SCD combined with O-type blood group confers some resistance to *Plasmodium falciparum* infection.17 O-blood type has also been found to confer resistance towards developing cardiovascular disease or type-2 diabetes.

Blood group typing by traditional serological, molecular or SNP microarray methods have limited functionality in characterizing blood group antigens which are rare, have weak expression or recombinant, partial, or novel.18 There is mounting evidence that clinically significant rare antigens and novel variants confound conventional serologic typing and SNP approaches. 19 To achieve extended RBC typing, a more comprehensive approach would be to apply Next-Generation Sequencing (NGS) to overcome the limitations of serological and molecular techniques.20-24 NGS approaches including targeted exome sequencing (TES), whole exome sequencing (WES), and whole genome sequencing (WGS) have the potential to provide a new basis of pre-transfusion testing by facilitating the accurate characterization of an individual’s complete blood group profile, supporting precision-based medicine.25,26 Curated and detailed DNA-phenotype-annotated databases storing blood group allele and antigen data has been developed by ISBT, providing tools for blood type calling from genetic data. As NGS is increasingly employed in immunohematology, Erythrogene, a custom-designed blood group allele database associated with the 1000 Genomes (1KGP3) project has been developed and applied.6,27

Accurate prediction of blood group phenotypes based on NGS data requires immuno-genetic knowledge since multiple genotypes can lead to the same phenotype (e.g. ABO, MNS and LE systems) and not all blood group antigens are direct products of primary genes like the ABO, LE and H systems.25 Even though various tools/algorithms have been developed applying statistical and machine learning approaches to process NGS data,28,29 there is no single tool which provides complete and comprehensive automation of blood group characterization in user-friendly manner. Only two tools are available to date, BOOGIE30 and BloodTyper,25 neither of which has the potential to identify novel blood group variants. To overcome this limitation, we developed a comprehensive and secure bioinformatics platform called RBCeq. RBCeq ([https://www.rbceq.org/](https://www.rbceq.org/)) is a web server-based blood group genotyping software. It is able to provide fast and accurate extended mass screening of blood groups, including from diverse populations with distinctive and complex blood group profiles.

Blood group antigen profiles are well characterized for European, North American, and some East Asian populations, but no extensive study has been carried out to date on the blood group antigen profile of the Australian population. Australia constitutes a highly heterogeneous multi-cultural and multi-ethnic population. 31,32 The patient population in Australia is more diverse than blood donors and there is a greater need in patients for blood types that are rare in a Caucasian population. Most donors are Caucasian, previous studies have indicated that immigrants and ethnic minorities are less involved in blood donation33. To date, there has not been any study on the blood group polymorphisms within the Caucasian Australian population as a whole. The problem is compounded by the absence of rare blood group antigen reports on a large number of Australian individuals. However, recently the Medical Genome Reference Bank (MGRB) was released comprising WGS data of 4000 healthy elderly Australians mostly of European descent.34,35 These individuals are participants of the ASPirin in Reducing Events in the Elderly (ASPREE) study, initiated to investigate whether the daily use of aspirin would prolong the healthy life span of older adults. The largest ongoing study of healthy ageing in the Southern Hemisphere36. In this study, we use genomic data from the MGRB ASPREE participants to identify unique blood group variants among a healthy older Australian general population and compared this with global data. To the best of our knowledge, this is the first study using WGS data for identifying the prevalence of blood group polymorphisms within a large Australian population dataset compared with global trends. Results of this study will help characterize key blood group systems and the frequency of rare phenotypes amongst the Australian general population, providing a tool to improve treatment accuracy and reduce the rate of hemolytic transfusion reactions. Also, providing an important proof-of-concept for using WGS for higher resolution RBC typing in the future.

## METHODS

### Study design

The design of this study included multiple interrelated components (Figure 1). We collated 2,796 WGS data from the ASPREE study, and 2504 from the 1KGP3 phase 3 (1KGP3) databases. Descriptions of the database ASPREE35 and 1KGP337 are included in Supplementary Table 1. The gVCF and BAM files were used to predict genotype and phenotype using RBCeq. We collated 96 709 high-quality coding variants from genomes of ASPREE participants. Allele frequencies and clinical relevance of the variant were obtained using ANNOVAR38 (gnomAD genome collection [v2.1.1]). Additionally, sequencing coverage and sequence identity of the blood group genes were extracted from the alignment files using BAMSTAT39 and IGV40, respectively.

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/T1)

Table 1: 
The Frequency of null phenotypes identified in ASPREE data compared with and 1KG dataset.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/22/2021.04.18.21255241/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/F1)

Figure 1: 
An overview of workflow to comprehensively characterise population specific blood group variants and phenotypes. The population stratification was performed using plink; blood group profiling using RBCeq and blood group gene coverage calculation using BAMstat and bedtools.

To evaluate the distinction between the ASPREE dataset and other global populations, we merged 789 linkage disequilibrium pruned coding variants of the ASPREE and 1KGP3 datasets, all of which had minor allele frequencies (MAF) ≥ 0.1. We then utilized PLINK v1.941 to generate a principal component analysis (PCA) plot from the resultant merged data and plotted using ggplot242 (Figure 2). Next, we created a Circos plot containing variant frequencies and the gene annotation using R libraries circlize43, tidyverse44, dplyr45, ComplexHeatmap46, and stringr44 (Figure 3). The prediction of RHD zygosity and C/c antigen was confirmed by RBCeq using the RH blood group gene coverage plot. Initially, all sample RH genes coverage was calculated using BedTools47 with bin size of 1. The coverage data was then smoothed and downscaled by averaging every 300 bases. The smoothed data were plotted using R libraries.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/22/2021.04.18.21255241/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/F2)

Figure 2. 
PCA plots showing 789 (MAF ≥ 0.1 and LD prune) markers in 5300 individuals from the ASPREE and 1KGP3 datasets. The X-axis denotes the value of PC1, while Y-axis denotes the value of PC2, with each dot in the figure representing one individual. The first two principal components shown here account for ∼80% of the observed variance in the combined dataset. (AFR: African, AMR: American, EAS: East Asian, EUR: European, SAS: South Asian, LD: linkage disequilibrium, MAF: minor allele frequency, PCA: principal component analysis.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/22/2021.04.18.21255241/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/F3)

Figure 3: 
The distribution of gnomAD (genome and exome) datasets genetic variants and their frequency in RBC antigen encoding genes. The outer ring (red) represents the RBC antigen encoding genes; box length represents the number of variants observed; G denotes gnomAD genome frequency, and E denotes gnomAD exome frequency. The outer green (light/dark) circle indicates the distribution of variants frequency across different blood group genes from the gnomAD data. The red (light/dark) circle indicates the number of variants with ISBT associations relative to all gnomAD variants. The blue (light/dark) circle indicates the distribution of rare non-ISBT variants in all six populations. The dark grey circle indicates the number of non-ISBT variants annotated to the ClinVar database. The yellow circle shows the distribution of the number of novel variants. The data was visualised using Circos (reference)

### Genotyping using RBCeq

The RBCeq algorithm determined the genotype and predicted phenotype of the ASPREE and 1000 genome phase three cohort individuals([https://www.rbceq.org/](https://www.rbceq.org/)). RBCeq predicts blood group profiles based on both SNV and CNV data. The gVCF of each sample was used as input along with the blood group gene coverage calculated using BAMTrimmer ([https://github.com/MayurDivate/BamTrimmer](https://github.com/MayurDivate/BamTrimmer)). The RBCeq predicted known blood group profiles for 36 blood group systems and two transcription factors, and identified ClinVar, Rare (<=0.05 in gnomAD dataset) and deleterious novel variants.

## RESULTS

### The blood group gene variant landscape in healthy older Australians

The PCA plot revealed blood group genetics for the ASPREE and European participants from the 1KGP3 datasets exhibited a unique genomic identity distinct from other global population groups. The nearest correlation to ASPREE and European blood group genetics was found in the American population. Clear differentiation between the African, East Asian and South East Asian data was observed (Figure 2).

### Distribution of null phenotypes of high prevalence antigens

Caucasian people with rare null phenotypes are at risk of alloimmunisation if transfused. We found 0.21% ASPREE participants to express rare Colton Co(a-b+) phenotype. The rare Vel-negative blood type was found only among the ASPREE participants (0.04%).48 In the YT system, the Yt(a-b+) subtype was present in 0.14% of ASPREE, 0.20% of European, and 0.20% of South Asian participants. Transfusion of alloimmunised Co(a-) and Vel-individuals would require antigen negative blood to avoid HTR.

### Distribution of weak, partial and null antigens

Analysis of the distribution of the weak, partial and null antigens among the representative populations from the two databases showed a 93.34% distribution of the Duffy null phenotype Fy(a-b-) amongst African participants, which was higher than that reported in previous studies (68%-70%) (Table 2). Apart from African populations, Duffy null was only found in the American population (African-American). The Duffy null phenotype was not found amongst all other populations consistent with earlier reports.11,49 We report for the first time the distribution of a rare parabombay H+□ type (*FUT2*01W*.*02*.*01*), which was found only among the ASPREE (0.07%, n=2) and East Asian (19.05%, n=96) populations and was not represented in the European dataset. We also found a very low percentage of the extremely rare Kel 6,-7 antigen (0.15%) previously reported frequency in African populations is 19% and in Caucasians <1%50. The rare LU:-13 phenotype was found among the ASPREE participants (0.21%) previously reported only amongst European populations (0.40%).

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/T2)

Table 2: 
Comparing Weak, Partial and Null Blood group antigens distribution in the ASPREE and 1KG datasets.

Not a single individual with a very rare Jk(a–b–) phenotype, which is most frequent in Polynesian populations, was found among all six populations studied (Table 2). The Jk(a+wb-) phenotype which causes weak or partial expression of the Jka antigen, was represented at frequencies 0.54% (ASPREE), 4.54% (African), 4.61% (American), 1.19% (European), 15.87% (East Asians) and 8.38% (South Asian) within the six population datasets. We found the highest prevalence of the Jr(a+□) phenotype among East Asian populations (7.9%) while the African participants completely lacked the phenotype. Similar prevalence of Jr(a+) was observed among the other populations in our study (ASPREE [1.4%], European [1.2%], South Asian [0.6%]). 11 The A3 (0.15%) subtype was only observed in the African population, while the B3 phenotype was observed in both ASPREE (1.65%) and African (0.15%) populations.

### Prevalence of RHD/RHCE blood group phenotypes

We compared the frequency of Rh blood group phenotypes frequencies in ASPREE and 1KGP3 samples with previously reported data (Table 3). We analyzed the D gene percentage with respect to homozygous null expression within the six populations. Our results indicated that 99.8% of East Asian, 96.52% of African, 93.95% of American, and 94.48% of South Asian participants were homozygous for the RHD gene. Previous reports state approximately 85% of the Caucasian population is D-positive and RBCeq also predicted 83. 9% in European populations consistent with earlier reports (**Supplementary Figure 1**). In the 474 (16.92%) D-negative samples, 473 were due to deletion of the RHD gene and one sample was homozygous for the RHD*01N.20 allele. This is possibly the first to report for partial RHD antigen named DIIItype4 and DUC2 in Australian Caucasian cohorts. In Aspree two hemizygous calls for DIIItype4 and two homozygous calls for DUC2 were observed. Other than ASPREE, only in African populations DIIItype4 and only in American population DUC2 were observed. These two partial D types are rarely observed in Caucasians (0.1% European Americans), the same rarity we observed in other 1KG populations.

View this table:
[Table 3:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/T3)

Table 3: 
RHD/RHCE blood group phenotype frequency prevalence in ASPREE and 1KGP3 samples compared with previously reported data *RHD*DIII* ()0.04% and *RHD*DUC2*

The C antigen occurs when reads of exon 2 of RHCE are mapped to exon 2 of RHD. It results in an increase of aligned reads to RHD exon 2 with no alignment to RHCE exon 2 sequence51-53 (Supplementary Figure 2). Our data for C/c antigen distribution are similar to earlier reports except for the ASPREE, European and South Asian participants. Our predictions showed the prevalence of C/c antigen among the ASPREE and European participants was 40-50% lower than previous report. In South Asian, the C/c antigen prevalence was also 20-25% lower than previous reports. We found the prevalence of the Ece and CEce phenotypes was two-fold higher in the ASPREE and European populations when compared with previous data54. We also detected partial c and e phenotypes in the 1KGP3, only in African population. Two individuals from the European population exhibited the *RHCE*02*.*08*.*01* allele at a homozygous level which encodes low prevalence antigen Cw, a partial C and MAR null phenotype. *RHCE*02*.*08*.*01* is clinically significant and is associated with HDFN, has been reported more frequently in the Polish population. There is a demand for MAR-blood in Poland 55. The frequency of *RHCE*02*.*08*.*01* allele phenotype has not been previously reported in a large scale study.

### Comprehensive analysis of MNS phenotype

Table 4 shows the frequency of phenotype M+N-S+s- and M+N+S+s- in additional ethnic (EAS/SAS/EUR/AMR) data. The frequency of M+N-S+s+, M+N+S-s+, and M+N+S+s+ was previously reported for Caucasian and African population; we observed a two-fold lower prevalence of all three blood types in European participants compared with earlier reports (4.97% vs 14.0%, 7.75% vs 22.0%, and 12.92% vs 24.0% respectively). Conversely, the prevalence of M-N+S+s+ phenotype was found to be five-fold higher in the European participants (30.02%) in our study compared with previous data (6.0%). We found prevalence of all four MNS phenotypes in the ASPREE participants to be similar to those reported for Caucasian populations earlier.

View this table:
[Table 4:](http://medrxiv.org/content/early/2021/04/22/2021.04.18.21255241/T4)

Table 4: 
Comparison of MNS blood group phenotypes and frequencies observed in the ASPREE and 1KGP3 datasets.

### Investigation of novel and clinically significant blood group variants amongst healthy older Australians

We characterized variants within the ASPREE population dataset that have not previously been reported to be associated with blood group alleles but have the potential to affect the formation of antigenic structures (Figure 3). We detected nine variants that had clinical associations (Supplementary Table 2). There were 688 rare variants with frequencies of ≤ 0.05 (genomAD) amongst the ASPREE participants’ data (Supplementary Table 3). Notably, we identified 149 novel variants that were computationally predicted to be deleterious (Figure 3, Supplementary Table 4). Predictably, no clinically significant, rare and novel variants were found to be associated with the genetically complex and most studied ABO blood group system, and only one novel and 16 rare variants were associated with the RH blood group system, which is also amongst the most studied blood group systems.

## Discussion

In this study, we compared the blood group genotype profile of Australian participants from the ASPREE study with African, American, Asian, and European population data from the 1KGP3 and gnomAD databases. The population stratification analysis revealed the genetic makeup of the ASPREE and European participants distinct from the other populations and its likely because many ASPREE participants were of European descent (Figure 2). We also observed a clear differentiation between the African participants and other populations, consistent with previous reports on the evolution of human blood group genetics.

We analyzed the distribution of rare phenotypes amongst the six population datasets (Table 1). We found the Colton Co(a-b+) phenotype was high frequency amongst the ASPREE participants, consistent with European ancestry. The Co(a-b+) phenotype is extremely rare and we found this predicted phenotype dataset at a frequency of 0.21% (n=6) in the ASPREE dataset, also consistent with frequencies reported in Europe. Previously, seven separate studies on 13 460 donors from Northern Europe and North America showed 91.3% prevalence of Co(a+b-) and 8.5% prevalence of Co(a+b+) and only 0.2% prevalence of Co(a-b+).56 A study on 1706 African Americans showed a 100% prevalence of the Co(a+) phenotype.56.

In ASPREE, participants with the rare Vel-negative phenotype was also predicted. Antibodies to the Vel antigen are known to be associated with complement activation. Vel is a high frequency antigen, inherited as an autosomal recessive trait, that shows variable strength, ranging from strong to weak.21,57 Vel-negative individuals can develop anti-Vel antibodies after transfusion or pregnancy and may develop acute HTRs if transfused with Vel-positive blood.21,57 Therefore it is of value to screen the population to identify Vel-blood donors and add to the national frozen blood bank to fulfil future requests.

We also found the low prevalence ABO alleles A3 and B3 only among the African (A3:0.15%, B3:0.15%)and ASPREE (B3:1.65%) participants. The B3 prediction of ASPREE is with B genotype possibilities, and may require further serology work to confirm the genotype reflects phenotype. Earlier studies have indicated that distinct ABO variants were associated with different populations.58

It is important to type clinically significant RBC antigens of patients and blood donors in a pretransfusion setting to ensure blood components are matched appropriately for transfusion, to prevent transfusion related adverse events and reduce alloimmunisation.19 We typed clinically significant rare and weak blood group variants among the six populations (Table 2). The RBCeq prediction showed a significantly higher (93.34%) prevalence of the Duffy null Fy(a-b-) phenotype in the African dataset compared with earlier reports (68%-70%)11.

The observed Fy(a-b-) phenotype was predicted from genotype FY\*02N.01/FY\*02 by RBCeq. The phenotype FY(a-b-) sometimes may be genotypically FY\*01N.01/\*02, in such cases the patient is unlikely to form anti-Fy3 as the patient will express Fyb in tissues. We found para-Bombay phenotypes (H+□) amongst the ASPREE (0.07%) and East Asian (19.05%) participants. However, in the Chinese population, a greater frequency of para-Bombay than Bombay has been reported.11,59 Our results showed a very low percentage of the extremely rare KEL6,-7 antigen only amongst the African dataset (0.15%).

We also identified the extremely rare LU:-13 phenotype among the ASPREE participants (0.21%) which has only been observed previously among European populations. The Lutheran blood group system comprises 20 antigens defined by largely clinically benign antibodies. However, no specific data is available on the adverse events of transfusion of anti-LU13 due to the extreme rarity of the antibody. Therefore, there is a recommendation to transfuse Lu(a-b-) blood and do family studies to identify potential donors if a patient/donor is found to lack the LU:13 antigen.

We didn’t find any evidence of the Jk(a–b–) phenotype, which is found in Polynesians, amongst the six populations studied. Red Cross Lifeblood has appealed for blood donors of Polynesian heritage to donate in an attempt to identify more Jk(a-b-) donors as there is demand for blood of this phenotype. The weak or partial expression, Jk(a+w) antigen was observed in all six populations with the highest prevalence seen among the East Asian participants. The Jk(a+w) phenotype can be a basis of discordance between serological and molecular methods. A phenotype with weakened Jka expression may constitute a risk for HTRs if antigen-positive units are not identified during pre-transfusion/donation serology tests.

Anti-Jra has been reported to cause fatal cases of HDN and delayed HTRs60,61. We found the highest prevalence of the Jr(a+□) phenotype among the East Asian (7.9%) populations while the African population completely lacked the phenotype. Similar prevalence of prediction of the Jra antigen was observed among the other populations in our study. As with previous studies, we found the Kel6,-7phenotype only in the African dataset.11

Fetal anemia is one of the most serious consequences of paternal blood group incompatibility with the mother. The most clinically significant antibodies involved in HDFN are the RhD antigen, and other minor Rh antigens C, E, c, e antigens. Antigens of other blood group systems (Duffy, Kidd, Kell, M and S) rarely cause significant problems, however it is a requirement9 to test and monitor pregnant women for the formation of these antibodies, as some can cause HDN.62. For anti-D, anti-K, and anti-c, there is a >50% risk of mild to severe HDFN developing in case the fetus inherits the target paternal red cell antigen. We analyzed the RHD gene haplotypes percentage in all six populations (Table 3). Our results indicated the majority (99. 8%) of the East Asian participants, and most of the African (96.52%), American (93.95%) and South Asian (94.48%) participants were homozygous for the *RHD* gene. However, the European (83.9%) and ASPREE (83.09%) participants had lower homozygous *RHD* gene frequencies. Genetic variants that result in weakened or partial antigen expression are problematic in transfusion medicine and can result in mistyping of RBCs. We also identified rare partial RHD types *DUC2* (0.07%) and DIIItype4 (0.07%), observed in Caucasians (approx.0.1% European Americans). Partial D proteins lack certain D epitopes which need to be characterised in order to guide the selection of blood for transfusion and in antenatal care, as if these individuals are transfused with D+ blood (all D epitopes present) they can be alloimmunised to the epitope they lack. And individuals with these antigens should ideally receive D-negative donor RBCs to avoid alloimmunization. In some cases. individuals with weak D should receive D-negative RBCs to avoid alloimmunisation but are often mistyped as D positive with a risk of producing an alloanti-D antibody.25,53 Additionally, It is important to differentiate weak D from weak D/Partial and Partial D types in pregnant women to provide anti-D prophylaxis during pregnancy and in women of childbearing potential to prevent D alloimmunisation..25

RBCeq predictions for the C-antigen prevalence among the six population datasets followed similar patterns to earlier studies for the predominantly African and East Asian (Table 3). However, C-antigen prevalence in our predictions for the ASPREE, European and South Asian datasets was significantly lower than previously reported values63.

The points of contrast presented for the blood group profiles of the ASPREE cohort are of potential clinical significance. In our study, the ASPREE and European participants had two-fold higher prevalence of C-E+c+e+ and C+E+c+e+ compared with previous reports for Caucasian populations. We could also found two instances of Cw, a partial C and MAR null within the European population never reported before using genomics data.

To the best of our knowledge, this is the first study describing the frequency of M+N-S+s- and M+N+S+s- in different ethnic backgrounds using genomic data (Table 4). We found the prevalence of all MNS phenotypes among the ASPREE participants to be similar to those reported earlier for Caucasian populations.

The prevalence of blood group antigen variants based on exon and splice-site variants were compared to global population groups (Figure 3). The results reflected a high level of conservation in the majority of blood group systems. We detected nine variants that had pathogenic clinical associations (Supplementary table2). We also detected 688 rare variants (238 with greater than two-fold difference from the gnomAD exome) with minimum allele frequencies of ≤ 0.05 in any of the five populations (Supplementary Table 3). Most importantly, we identified 149 variants that were computationally predicted to be novel and deleterious (Supplementary Table 4). Predictably, no clinically significant, rare, and novel variants were found to be associated with the most studied genetically complex ABO blood group system, and only one novel and 16 rare variants was associated with the RH blood group system. This observation is expected from these two most commonly studied systems.

While this is possibly the first comprehensive multi-ethnic study where genetic data from six different populations was compared, certain limitations need to be acknowledged. The ASPREE dataset, though heterogeneous, predominantly included participants of European descent with Caucasian population values. The dataset also did not include investigation of the hybrid RH and MNS changes. Serological testing and exon level detections were not performed to back up the RBCeq phenotype predictions. However, despite the limitations, our results are providing new insights into blood group genetics and antigen frequency.

The frequency of inherited disorders like SCD and thalassaemia are increasing in Australia due to increased migration between the populations.64 RH alloimmunisation is not completely prevented in patients with SCD although matched for extended blood groups due to antigen polymorphism in the RH system. These are examples of the significant differences between Australia’s blood donors, who are mostly Caucasian and African SCD transfusion recipients.64 Alpha thalassaemia has been identified in Australian patients in Aboriginal and Torres Strait Islander communities in the Northern Territory and northern Western Australia.65.65 Rare blood group identification will improve rare blood unit management. Since language and culture are very diverse in Australian communities, additional research is needed to guide transfusion practice and improve transfusion outcomes in Australia.

## Supporting information

Supplementary data [[supplements/255241_file02.xlsx]](pending:yes)

## Data Availability

The Medical Genome Reference Bank

*   Received April 18, 2021.
*   Revision received April 18, 2021.
*   Accepted April 22, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Schoeman EM, Roulis EV, Liew YW, et al. Targeted exome sequencing defines novel and rare variants in complex blood group serology cases for a red blood cell reference laboratory setting. Transfusion. 2018;58(2):284–293.
    
    

2.  2.Lane WJ, Westhoff CM, Uy JM, et al. Comprehensive red blood cell and platelet antigen prediction from whole genome sequencing: proof of principle. Transfusion. 2016;56(3):743–754.
    
    

3.  3.Storry JR, Clausen FB, Castilho L, et al. International Society of Blood Transfusion Working Party on Red Cell Immunogenetics and Blood Group Terminology: Report of the Dubai, Copenhagen and Toronto meetings. Vox Sanguinis. 2019;114(1):95–102.
    
    

4.  4.Transfusion ISoB. Red Cell Immunogenetics and Blood Group Terminology 2021.
    
    

5.  5.Guo Y, Busch MP, Seielstad M, et al. Development and evaluation of a transfusion medicine genome wide genotyping array. Transfusion. 2019;59(1):101–111.
    
    

6.  6.Möller M, Jöud M, Storry JR, Olsson ML. Erythrogene: a database for in-depth analysis of the extensive variation in 36 blood group systems in the 1000 Genomes Project. Blood Adv. 2016;1(3):240–249.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmxvb2RvYSI7czo1OiJyZXNpZCI7czo3OiIxLzMvMjQwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMjIvMjAyMS4wNC4xOC4yMTI1NTI0MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

7.  7.Fung MK, Eder A, Spitalnik SL, Westhoff CM, American Association of Blood B. Technical manual, 2017.
    
    

8.  8.Gleadall NS, Veldhuisen B, Gollub J, et al. Development and validation of a universal blood donor genotyping platform: a multinational prospective study. Blood Adv. 2020;4(15):3495–3506.
    
    

9.  9. Ltd AaNZSoBT. FINAL-Guideline\_-for\_Transfusion\_and\_Immunohaematology\_Laboratory\_Practice. 2020.
    
    

10. 10.Flegel WA. ABO genotyping: the quest for clinical applications. Blood Transfus. 2013;11(1):6–9.
    
    

11. 11.Reid ME L-FC, Olsson ML. The Blood Group Antigen FactsBook. Vol 3. Elsevier. 2012;Vol 3.
    
    

12. 12.Chou ST, Westhoff CM. The Rh and RhAG blood group systems. Immunohematology. 2010;26(4):178–186.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22356455&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

13. 13.de Mattos LC. Molecular polymorphisms of human blood groups: a universe to unravel. Rev Bras Hematol Hemoter. 2011;33(1):6–7.
    
    

14. 14.da Silva-Malta MCF, Sales CC, Guimarães JC, et al. The Duffy null genotype is associated with a lower level of CCL2, leukocytes and neutrophil count but not with the clinical outcome of HTLV-1 infection. J Med Microbiol. 2017;66(8):1207–1216.
    
    

15. 15.Kulkarni H, Marconi VC, He W, et al. The Duffy-null state is associated with a survival advantage in leukopenic HIV-infected persons of African ancestry. Blood. 2009;114(13):2783–2792.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTE0LzEzLzI3ODMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNC8yMi8yMDIxLjA0LjE4LjIxMjU1MjQxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

16. 16.Al Huneini M, Alkindi S, Panjwani V, et al. Increased Vasoocclusive Crises in “O” Blood Group Sickle Cell Disease Patients: Association with Underlying Thrombospondin Levels. Mediterr J Hematol Infect Dis. 2017;9(1):e2017028.
    
    

17. 17.Amodu OK, Olaniyan SA, Adeyemo AA, Troye-Blomberg M, Olumese PE, Omotade OO. Association of the sickle cell trait and the ABO blood group with clinical severity of malaria in southwest Nigeria. Acta Trop. 2012;123(2):72–77.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22503377&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

18. 18.McBean RS, Hyland CA, Flower RL. Approaches to determination of a full profile of blood group genotypes: single nucleotide variant mapping and massively parallel sequencing. Comput Struct Biotechnol J. 2014;11(19):147–151.
    
    

19. 19.Kulkarni S, Maru H. Extended phenotyping of blood group antigens: Towards improved transfusion practices. Global Journal of Transfusion Medicine. 2020;5(2):120–125.
    
    

20. 20.McBean R, Liew YW, Wilson B, et al. Genotyping confirms inheritance of the rare At(a-) type in a case of haemolytic disease of the newborn. J Pathol Clin Res. 2016;2(1):53–55.
    
    

21. 21.Storry JR, Joud M, Christophersen MK, et al. Homozygosity for a null allele of SMIM1 defines the Vel-negative blood group phenotype. Nat Genet. 2013;45(5):537–541.
    
    

22. 22.Cvejic A, Haer-Wigman L, Stephens JC, et al. SMIM1 underlies the Vel blood group and influences red blood cell traits. Nat Genet. 2013;45(5):542–545.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2603&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23563608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

23. 23.Stenfelt L, Hellberg A, Moller M, Thornton N, Larson G, Olsson ML. Missense mutations in the C-terminal portion of the B4GALNT2-encoded glycosyltransferase underlying the Sd(a-) phenotype. Biochem Biophys Rep. 2019;19:100659.
    
    

24. 24.Omae Y, Ito S, Takeuchi M, et al. Integrative genome analysis identified the KANNO blood group antigen as prion protein. Transfusion. 2019;59(7):2429–2435.
    
    

25. 25.Lane WJ, Westhoff CM, Gleadall NS, et al. Automated typing of red blood cell and platelet antigens: a whole-genome sequencing study. Lancet Haematol. 2018;5(6):e241–e251.
    
    

26. 26.Roulis E, Schoeman E, Hobbs M, et al. Targeted exome sequencing designed for blood group, platelet, and neutrophil antigen investigations: Proof-of-principle study for a customized single-test system. Transfusion. 2020.
    
    

27. 27.Abecasis GR, Auton A, Brooks LD, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11632&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23128226&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310434500030&link_type=ISI) 

28. 28.Mills RE, Walter K, Stewart C, et al. Mapping copy number variation by population-scale genome sequencing. Nature. 2011;470(7332):59–65.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09708&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21293372&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286886400033&link_type=ISI) 

29. 29.Lofling JC, Hauzenberger E, Holgersson J. Absorption of anti-blood group A antibodies on P-selectin glycoprotein ligand-1/immunoglobulin chimeras carrying blood group A determinants: core saccharide chain specificity of the Se and H gene encoded alpha1,2 fucosyltransferases in different host cells. Glycobiology. 2002;12(3):173–182.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/glycob/12.3.173&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11971861&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175501800005&link_type=ISI) 

30. 30.Giollo M, Minervini G, Scalzotto M, Leonardi E, Ferrari C, Tosatto SC. BOOGIE: Predicting Blood Groups from High Throughput Sequencing Data. PLoS One. 2015;10(4):e0124579.
    
    

31. 31.Delaney M, Harris S, Haile A, Johnsen J, Teramura G, Nelson K. Red blood cell antigen genotype analysis for 9087 Asian, Asian American, and Native American blood donors. Transfusion. 2015;55(10):2369–2375.
    
    

32. 32.Volken T, Crawford RJ, Amar S, Mosimann E, Tschaggelar A, Taleghani BM. Blood Group Distribution in Switzerland - a Historical Comparison. Transfus Med Hemother. 2017;44(4):210–216.
    
    

33. 33.Tanya Davison JD, Robert Flower Australia’s ethnic face is changing, and so are our blood types. 2019.
    
    

34. 34.Pinese M, Lacaze P, Rath EM, et al. The Medical Genome Reference Bank contains whole genome and phenotype data of 2570 healthy elderly. Nat Commun. 2020;11(1):435.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-14079-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

35. 35.Lacaze P, Pinese M, Kaplan W, et al. The Medical Genome Reference Bank: a whole-genome data resource of 4000 healthy elderly individuals. Rationale and cohort design. Eur J Hum Genet. 2019;27(2):308–316.
    
    

36. 36.Banks E, Redman S, Jorm L, et al. Cohort profile: the 45 and up study. Int J Epidemiol. 2008;37(5):941–947.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dym184&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17881411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259771500008&link_type=ISI) 

37. 37.Genomes Project C, Auton A, Brooks LD, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26432245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

38. 38.Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq603&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20601685&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

39. 39.Gan HM, Falk S, Morales HE, Austin CM, Sunnucks P, Pavlova A. Genomic evidence of neo-sex chromosomes in the eastern yellow robin. Gigascience. 2019;8(12).
    
    

40. 40.Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–192.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bib/bbs017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22517427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

41. 41.Purcell S. PLINK v1.9 2019.
    
    

42. 42.H W. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. 2016, ISBN 978-3-319-24277-4.
    
    

43. 43.Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize Implements and enhances circular visualization in R. Bioinformatics. 2014;30(19):2811–2812.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24930139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000343082900017&link_type=ISI) 

44. 44.ickham H AM, Bryan J, Chang W, McGowan LD, François R, Grolemund G, Hayes A, Henry L, Hester J, Kuhn M, Pedersen TL, Miller E, Bache SM, Müller K, Ooms J, Robinson D, Seidel DP, Spinu V, Takahashi K, Vaughan D, Wilke C, Woo K, Yutani H. Welcome to the tidyverse. Journal of Open Source Software. 2019;4(43), 1686. doi: 10.21105/joss.01686.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21105/joss.01686&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15461798&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

45. 45.Hadley Wickham RF, Lionel Henry and Kirill Müller. dplyr: A Grammar of Data Manipulation. 2021.
    
    

46. 46.Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–2849.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw313&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27207943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

47. 47.Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–842.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20110278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275243500019&link_type=ISI) 

48. 48.<Human blood groups. 3rd ed.pdf>.
    
    

49. 49.Hyland CA, Roulis EV, Schoeman EM. Developments beyond blood group serology in the genomics era. British Journal of Haematology. 2019;184(6):897–911.
    
    

50. 50. Kristin Stendahl Cat,  Ian L. Baine. Chapter2: Methods of RBC Alloimmunization toABO and Non-ABO Antigens, and TestMethodologies, Immunologic Concepts in Transfusion Medicine. 2020.
    
    

51. 51.Wagner FF, Flegel WA. Review: the molecular basis of the Rh blood group phenotypes. Immunohematology. 2004;20(1):23–36.
    
    

52. 52.Chou ST, Flanagan JM, Vege S, et al. Whole-exome sequencing for RH genotyping and alloimmunization risk in children with sickle cell anemia. Blood Adv. 2017;1(18):1414–1422.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmxvb2RvYSI7czo1OiJyZXNpZCI7czo5OiIxLzE4LzE0MTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNC8yMi8yMDIxLjA0LjE4LjIxMjU1MjQxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. 53.Schoeman EM, Lopez GH, McGowan EC, et al. Evaluation of targeted exome sequencing for 28 protein-based blood group systems, including the homologous gene systems, for blood group genotyping. Transfusion. 2017;57(4):1078–1088.
    
    

54. 54.Mbalibulha Y, Muwanguzi E, Mugyenyi G. Rhesus blood group haplotype frequencies among blood donors in southwestern Uganda. J Blood Med. 2018;9:91–94.
    
    

55. 55.G. D. Human blood groups. 3rd ed. Oxford: Wiley-Blackwell. 2013;3.
    
    

56. 56.Colton Blood Group System. Human Blood Groups, 2013:384–390.
    
    

57. 57.Cvejic A, Haer-Wigman L, Stephens JC, et al. SMIM1 underlies the Vel blood group and influences red blood cell traits. Nat Genet. 2013;45(5):542–545.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2603&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23563608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

58. 58.Ying YL, Hong XZ, Xu XG, et al. Molecular Basis of ABO Variants Including Identification of 16 Novel ABO Subgroup Alleles in Chinese Han Population. Transfus Med Hemother. 2020;47(2):160–166.
    
    

59. 59.Patel JN, Donta AB, Patel AC, Pandya AN, Kulkarni SS. Para-Bombay phenotype: A case report from a tertiary care hospital from South Gujarat. Asian J Transfus Sci. 2018;12(2):180–182.
    
    

60. 60.Kim H, Park MJ, Sung TJ, et al. Hemolytic disease of the newborn associated with anti-Jra alloimmunization in a twin pregnancy: the first case report in Korea. Korean J Lab Med. 2010;30(5):511–515.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20890084&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F22%2F2021.04.18.21255241.atom) 

61. 61.Peyrard T, Pham BN, Arnaud L, et al. Fatal hemolytic disease of the fetus and newborn associated with anti-Jr. Transfusion. 2008;48(9):1906–1911.
    
    

62. 62.Johnsen JM. Using red blood cell genomics in transfusion medicine. Hematology Am Soc Hematol Educ Program. 2015;2015:168–176.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiYmxvb2Rib29rIjtzOjU6InJlc2lkIjtzOjEwOiIyMDE1LzEvMTY4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMjIvMjAyMS4wNC4xOC4yMTI1NTI0MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. 63.L. D. Blood Groups and Red Cell Antigens. Bethesda (MD): National Center for Biotechnology Information (US). 2005.
    
    

64. 64.Crighton G, Wood E, Scarborough R, Ho PJ, Bowden D. Haemoglobin disorders in Australia: where are we now and where will we be in the future? Intern Med J. 2016;46(7):770–779.
    
    

65. 65.Gaff C NJ, Saleh M <genetics\_in\_family_medicine.pdf>. National Health and Medical Research Council (NHMRC). 2007.