The genomic landscape of blood groups in Indigenous Australians in remote communities.

BACKGROUND: Red blood cell (RBC) membrane-associated blood group systems are clinically significant. Alloimmunisation is a persistent risk associated with blood transfusion owing to the antigen polymorphisms among these RBC-associated blood groups. Next-generation sequencing (NGS) offers an opportunity to characterize the blood group variant profile of a given individual. Australia comprises a large multiethnic population where most blood donors are Caucasian and blood group variants remain poorly studied among Indigenous Australians. In this study, we focused on the Tiwi Islanders, who have lived in relative isolation for thousands of years. METHODS AND MATERIALS: We predicted the blood group phenotype profiles in the Tiwi (457) and 1000 Genomes Phase 3 (1KGP3-2504) cohort individuals using RBCeq (https://www.rbceq.org/). The predicted phenotype prevalence was compared with the previous literature report.
RESULTS: We report, for the first time, comprehensive blood group profiles corresponding to the 35 known blood group systems among the Indigenous Tiwi islander population and identify possible novel antigen variants therein. Our results demonstrate that the genetic makeup of the Tiwi participants is distinct from that of other populations, with a low prevalence of LU (Au[a-b+]) and ABO (A2) and D+C+c+E+e- phenotype, an absence of Diego blood group variants, and a unique RHD (DIII type4) variant.
CONCLUSION: Our results may contribute to the development of a database of predicted phenotype donors among the Tiwi population and aid in improving transfusion safety for the ~2800 Tiwi people and the ~800,000 other Indigenous Australians throughout the nation.

chronic kidney disease

copy number variations

haemolytic transfusion reactions

Insertions/deletions

Next‐generation sequencing

Red blood cell

1000 genomes phase 3

Single nucleotide polymorphisms

INTRODUCTION

Blood transfusion has historically been an important facet of hematological care for which there is no viable alternative. There is immense diversity among human blood group systems, with blood group genes displaying a range of antigenic polymorphisms, including single nucleotide polymorphisms (SNP), insertions/deletions (indels), and copy number variations (CNVs). Misidentification or inappropriate matching of these variants may pose a risk of alloimmunisation that can result in serious hemolytic transfusion reactions (HTRs). , Next‐generation sequencing (NGS) is the most comprehensive approach to achieving extended red blood cell (RBC) typing, enabling the accurate characterization of the complete blood group variant profiles for a given individual. NGS approaches include whole‐exome sequencing (WES), whole‐genome sequencing (WGS), and targeted exome sequencing (TES). While these technologies are powerful, there is currently a lack of comprehensive software packages with the storage and processing capability necessary to handle large volumes of NGS data and to accurately predict blood group phenotypes. Recently, our group developed RBCeq, a web server‐based blood group genotyping software designed to improve blood type compatibility testing, and to facilitate the accurate mass screening of blood groups in diverse populations with distinctive and complex blood group profiles.

Australia is a multiethnic nation, with roughly one‐third of the population having been born overseas. However, most Australian blood donors are Caucasian. It is well established that the prevalence of antigens can differ among populations, potentially hindering efforts to maintain a safe and robust blood supply in multiethnic countries. In an earlier study, we carried out extended blood group molecular typing to detect rare antigens in a large Australian Caucasian cohort. Blood group antigen profiling for Western Desert Indigenous Australians is available. However, the Tiwi Islanders are culturally distinct among these Australian Indigenous communities. They lived in relative isolation for 7000–15,000 years and consider themselves to be distinct from mainland Aboriginal communities, with their own origin stories, language, and customs. Although isolated from mainland Indigenous groups for thousands of years, they had intermittent contact with visitors from Southeast Asia and Europe beginning in the seventeenth century. The Tiwi population, along with Aboriginal people from Arnhem Land, has a higher mortality rate than the Australian mainstream. Tiwi people are at high risk of kidney, cardiovascular, and metabolic diseases that significantly reduce their life expectancy relative to the non‐Indigenous population. , Alpha thalassaemia has also been identified in some Aboriginal and Torres Strait Islander communities in the Northern Territory and Northwestern Australia, requiring regular transfusion for affected individuals. Studies of blood group variants among the Tiwi people have never been performed. Developing an approach to careful blood group matching is thus vital for Indigenous individuals who required blood transfusion. Here, we characterize blood samples from Tiwi people with the aim of developing a database of blood type polymorphisms among Indigenous Australians in order to facilitate accurate donor matching and thereby reduce the incidence of transfusion reactions.

METHODS

Study population and datasets

Blood samples were collected between 2013 and 2014 as detailed in a previous study investigating the prevalence of chronic kidney disease (CKD) in the Tiwi population. Whole‐genome sequencing of 189 samples was performed in three batches by three different research groups. Additionally, 268 Tiwi samples were genotyped using Infinium Global Screening Array v3.0. Publicly available variant dataset from the 1000 Genomes Phase 3 (1KGP3) project was downloaded and used for comparative analyses. Allele frequencies for gnomAD populations were obtained from ANNOVAR (gnomAD genome collection [v2.1.1]).

Variant calling and annotation

Variant calling was performed according to GATK best practices v4.0.2.0 and Sentieon DNASeq v.201808.08 using the GRCh37 human reference genome. The resulting variant calls were filtered based on variant quality score recalibration (VQSR). Autosomal variants for 39 samples aligned to hg38 in the form of a processed variant call format (VCF) files were obtained from the National Centre for Indigenous Genomics. To achieve comprehensive blood group allele annotation, all variant calls were converted to the GRCh37 assembly for autosomal variants.

Prediction of blood group alleles and their phenotype using RBCeq

We predicted the blood group antigen profiles for the Tiwi (457) and 1KGP3 (2504) cohort individuals using RBCeq (https://www.rbceq.org/). The genomic variant call format (gvcf) of each sample was given as input with the blood group gene coverage calculated using the BAMTrimmer tool within RBCeq. RBCeq considers SNV, indels from 45 genes and calculates CNV to profile 36 blood groups and two transcription factors (Table S1). RBCeq first, scans each blood group for reference and alternative allele in the input vcf file. The zygosity of detected reference‐alternate allele guides the sub allele selection. The selected alleles are paired against each other and are assigned with a score. The score helps to select allele pair which considers most of the known present variants from the input vcf file. The detailed method is described in Jadhao et al., 2022.

The RBCeq output included the prediction of International Society of Blood Transfusion (ISBT) known blood group antigen alleles and phenotype. RBCeq also classified the filtered high‐quality deleterious variants which are not mapped to known ISBT alleles. First, it classified the variant as clinically significant by querying against ClinVar, then rare variants which have gnomAD MAF≤0.05 and potential novel variants which are without any public database annotation (Figure 1). The Rh copy number variation (D zygosity and C/c) prediction accuracy by RBCeq was confirmed using the RH blood group gene coverage plot. All samples RH genes coverage were considered using BedTools (bin size: 1), which were further normalized for every 300 bases pairs and plotted using R libraries.

FIGURE 1

An overview of the workflow to comprehensively characterize population‐specific blood group variants and phenotypes [Color figure can be viewed at wileyonlinelibrary.com]

To do the population stratification between the Tiwi dataset and other global populations, we carried out principal component analysis (PCA) of 695 variants (LD pruned and MAF≥0.1)) from 42 autosomal blood groups encoding gene sequences of TIWI and 1KGP3 dataset using PLINK v1.9 (Figure 1). The PCA results were plotted using ggplot2. The blood group encoding gene variants annotation/categorization was done using ANNOVAR. The Circos plot was created to show variants frequency and the gene annotation using R libraries like circlize, tidyverse, dplyr, ComplexHeatmap, and stringr.

RESULTS

The blood group gene variant landscape in the Tiwi population

Results from a PCA plot analysis revealed that the Tiwi participants exhibited a unique blood group genomic identity that was not closely correlated with that of the European population, who account for most blood donors in Australia. While clear differentiation was observed among the Tiwi, African, and East Asian participants, all three groups also distinctly differed from the rest of the global populations (Figure 2).

FIGURE 2

PCA plots showing 695 (MAF ≥0.1 and LD prune) markers in 454 individuals from the Tiwi 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 showed here account for ~80% of the observed variance in the combined dataset. AMR, American; AFR, African; EAS, east Asian; EUR, European; LD, linkage disequilibrium; MAF, minor allele frequency; PCA, principal component analysis; SAS, south Asian [Color figure can be viewed at wileyonlinelibrary.com]

Distribution of weak and partial antigens

We observed that the distribution of ABO antigens in Tiwi was distinct from that in other populations. Four sub‐phenotypes were predicted in ABO blood groups among this population: O (81.18%), A1 (18.38%), A2 (0.22%), and B (0.22%). The distribution of the O phenotype was very high compared with the other populations (Caucasian: 44%, African: 49%, and Asian: 43%). A1 and A2 phenotypes account for 99% of all groups A individuals. The A2 subgroup is considered a weak A subgroup that is most identified via the presence of an unexpected anti‐A or anti‐A1. This is the first study to report the prevalence of A2 in an Indigenous population, revealing it to be very rare (0.22%) among this Tiwi population (Table 1). The rare Jk(a−b−) was missing in all six populations, but Jk(a+wb−) predicted from JK*01 W.01 (c.130G > A) had a much higher frequency (16.19%) in Tiwi as compared to others 4.39% (African), 4.74% (American), 1.19% (European), and 8.38% (South Asian) except in East Asians (16.47%) (Table 1). Jk(a+wb−) antigen causes weak or partial expression of Jk(a+b−) and can be mistyped as Jk(a−).

TABLE 1

Comparison of blood group profiles with previously reported prevalence in other populations

Blood group	Predicted phenotype	RBCeq prediction (%)						Previously published (%) 27
		Tiwi	AFR	AMR	EUR	EAS	SAS	Caucasians	Africa	Asian
		n = 457	n = 661	n = 347	n = 503	n = 504	n = 489
ABO	A2	0.22 (n = 1)	11.50 (n = 76)	5.64 (n = 19)	13.92 (n = 70)	1.19 (n = 6)	5.32 (n = 26)	10	8	rare
ABO	B	0.22 (n = 1)	23.90 (n = 158)	5.93 (n = 20)	42.86 (n = 216)	10.54 (n = 53)	47.44 (n = 233)	9	20	27
RHD	DIII type4	1.75 (n = 8)	2.87 (Four homozygous and 15 heterozygous)	‐	‐	‐	‐	‐	‐	‐
LU	LU:18‐,19 or Au(a−b+)	0.22 (n = 1)	19.21 (n = 127)	7.12 (n = 24)	10.34 (n = 52)	0.99 (n = 5)	4.91 (n = 24)	‐	‐	‐
FY	Fy(a−b+)	0.88 (n = 3)	‐	24.04 (n = 81)	0.4 (n = 2)	33 (n = 166)	12.68 (n = 62)	34	22	0.3–3
JK	Jk(a−b+)	5.69 (n = 26)	5.14 (n = 34)	30.27 (n = 102)	25.79 (n = 130)	25.45 (n = 128)	14.31 (n = 70)	23	8	27
JK	Jk(a+w b−)	16.19 (n = 76)	4.39 (n = 29)	4.74 (n = 16)	1.19 (n = 6)	16.47 (n = 83)	8.38 (41)	‐	‐	‐
DO	Do(a+b−)	50.32 (n = 230)	6.51 (n = 43)	14.54 (n = 49)	1.19 (n = 6)	13.92 (n = 70)	12.88 (n = 63)	18	11	0.5–2
DO	Do(a−b+)	10.5 (n = 48)	53.86 (n = 356)	40.35 (n = 136)	81.55 (n = 411)	39.36 (n = 198)	35.99 (n = 176)	33	45	77–87

Abbreviations: AMR, American; AFR, African; EAS, East Asian; EUR, European; SAS, South Asian.

Genetic variants that lead to weakened or partial antigen expression can be difficult to type. D category IV is one of the most clinically relevant and phylogenetically diverse partial D. We detected a DIII type 4 frequency of 1.75% (n = 8) at a homozygous level among the Tiwi population (Table 1). Phenotyping a patient with partial D as D– is safe for their management and enables the prevention of anti‐D alloimmunisation. However, these patients are often misidentified as being D+ and they are at risk for alloanti‐D development. The mother with anti‐D and fetus with partial D are at risk for hemolytic disease of the fetus and newborn. Antibodies to Lutheran antigens are rarely clinically significant; most of its antigens have high prevalence, but Au(a−b+) phenotype is polymorphic and rare. A multiethnic population‐level frequency for the Au(a−b+) phenotype has not been reported previously. The Tiwi distribution of Au(a−b+) phenotype was 0.22% (n = 1), which was very low relative to the frequency in other populations (African, American, East Asian, European, South Asian) (Table 1).

Predicted blood group profiles contrasting with previous reports in other population

The protein sequence of the A and B glycosyltransferases (GTA and GTB) are encoded by the ABO gene, which differs by four amino acid residues, and a mutation in these sequences resulting in protein truncation is what leads to the O phenotype. B antigen is common in most of the population and we also observed the same through the 1KG population dataset analysis. However, the prevalence of B antigen was very low among the Tiwi population such that only one individual harboring the B antigen was detected (Table 1).

Low immunogenic antigens from the Dombrock blood group systems can cause mild‐to‐moderate and acute‐to‐delayed HTRs. Dombrock antigens are difficult to detect, and typing is not routinely performed in the context of routine pretransfusion practice. Additionally, their antigenic frequency represents a significant ethnic marker. Both Do(a+) and Do(b+) antigens are common among most global populations. However, we detected at least two times higher Do(a+b−) prevalence (50.32%) among Tiwi as compared with the 1KGP3 population and previous reports. Conversely, the Do(a−b+) predicted phenotype prevalence among the Tiwi population was less than half of that reported for other global populations and detected in the 1KGP3 population in the current study (Table 1).

It is well established that the Duffy (Fy) antigens can cause maternal immunization and subsequent rare HDFN. The Fy(a−b−) null phenotype is more common among African populations. An anti‐Fya African sickle cell anemia patient may require multiple blood transfusions. The Fy(a−b+) phenotype is observed in about 34% of Caucasians, 22% of African Americans, and 3% of the Asian population. Here, we detected a low prevalence of Fy(a−b+) among the Tiwi population (0.88%) which was comparable to that in the East Asian population (0.4%) (Table 1).

The clinically significant anti‐Jka and anti‐Jkb are weak antigens that may trigger HTRs and delayed transfusion reactions. Anti‐Jka and anti‐Jkb are difficult to detect in routine serology tests and notorious for their evanescence. We predicted a three times lower prevalence (5.69%) of the Jk(a−b+) phenotype among the Tiwi population compared with that predicted in the 1KGP3 population, and a 1.3 times lower prevalence than reported in previous studies (Table 1).

Frequency of RHD/RHCE blood group phenotypes

We also equated the RHD/RHCE blood group predicted phenotype frequencies of Tiwi with those from 1KGP3 dataset and previous reports. In the comparison of the prevalence of the D with respect to homozygous and null expression, our results revealed that 99.80% of East Asian, 96.52% of African, 93.95% of American, and 94.48% of South Asian participants were homozygous for the RHD‐gene (Table 2). However, we found a very high (100%) prevalence of homozygous RHD‐gene (Table 2; Figure S1) in 151 WGS Tiwi samples. The previous indigenous study also reported only two D‐negative individuals (Table 2).

TABLE 2

RHD blood group phenotype frequency prevalence in Tiwi and 1KGP3 dataset

Blood group Gene	Genotype	Predicted phenotype	RBCeq prediction (%)						Previously published (%) 27
			Tiwi n = 151	AFR n = 661	AMR n = 347	EUR n = 503	EAS n = 504	SAS n = 489	Western Desert Indigenous Australians	Caucasians	Africans	Asian
RHD	DD	D+	100 (n = 151)	96.52 (n = 628)	93.95 (n = 326)	83.90 (n = 422)	99.80 (n = 503)	94.48 (n = 462)	97.22 (n = 72)	85.00	92.00	99.00
RHD	D‐	D‐	0 (n = 0)	3.48 (n = 23)	6.05 (n = 21)	16.10 (n = 81)	0.20 (n = 1)	5.52 (n = 27)	2.77 (n = 2)	15.00	8.00	1.00

Abbreviations: AMR, American; AFR, African; EAS, East Asian; EUR, European; SAS, South Asian.

The highly polymorphic C antigen expression is driven by the number of exon 2‐ RHCE reads which are mapped to exon 2 of RHD , , (Figure S2). The prevalence of the D+C−c+E+e+ phenotype was very low in Tiwi and East Asian populations as compared to other populations (Table 3). The prevalence of the D+C+c+E−e+ phenotype was two times higher (10.28) as compared to other populations other than American and African populations (Table 3). The rare D+C+c+E+e− phenotype was predicted in two individuals of the Tiwi population however it was absent in the 1KGP3 dataset.

TABLE 3

RHCE blood group phenotype frequency prevalence in Tiwi and 1KGP3 dataset

Phenotype	RBCeq prediction (%)						Previously published (%) 27 , 42
	TIWI (457)	AFR (n = 661)	AMR (n = 347)	EUR (n = 503)	SAS (n = 489)	EAS (n = 504)	Caucasians
D+C−c+E+e+	0.44 (n = 2)	8.02 (n = 53)	5.48 (n = 19)	12.33 (n = 62)	4.7 (n = 23)	0.6 (n = 3)	10
D+C+c+ E−e	10.28 (n = 47)	13.77 (n = 91)	8.65 (n = 30)	4.57 (n = 23)	5.32 (n = 26)	5.16 (n = 26)	34.9
D+C+c+E+e−	0.44 (n = 2)	‐	‐	‐	‐	‐	<0.1

Abbreviations: AMR, American; AFR, African; EAS, East Asian; EUR, European; SAS, South Asian.

The MNS phenotype of the Tiwi people

We observed at least a 1.5 times higher prevalence of the M−N+S−s+ among the Tiwi participants as compared with the 1KGP3 population (Table S2). Conversely, the prevalence of M+N−S−s+ phenotype was predicted to be two times lower (2.41%) in the Tiwi participants as compared with other populations except in European (Table S2).

Investigation of novel blood group variants with potential clinical relevance

In a secondary analysis, we characterized variants with the potential to affect the antigen specificity but with no previous blood group allele associations within the Tiwi population dataset (Figure 3). We detected 13 variants that had clinical associations (Table S3), 5 of which had two times higher frequency in Tiwi. And 22 variants had Human Gene Mutation Database (HGMD) annotation (Table S4), five were associated with blood group antigen structure rs3088190 (Dombrock), rs5036 (Diego), rs8138197 (P1PK), rs5751348 (P1PK), CD44*c.255C > G (Indian). The rs201915844 (FUT2) was also detected which impairs the H antigen activity. Other 17 were associated with enhancer activity, porphyria, non‐secretory phenotype, and sedimentation (Table S4). There were 28 rare variants with frequencies of ≤0.05 among the gnomAD dataset (Table S5). Most importantly, 13 potential novel variants were identified that were predicted to be deleterious by any of the SIFT, Polyphen2, MutationTaster2, LRT, FATHMM, or PROVEAN in silico methods (Figure 3, Table S6).

FIGURE 3

Tiwi blood group gene encoding variants distribution as known (ISBT), GnomAD mapped, rare (<=0.05 in gnomAD), ClinVar annotated and potential novel variants. The RBC antigen encoding genes are shown by the outer ring (red); the number of variations found is represented by the box length; G denotes gnomAD genome frequency, and E denotes gnomAD exome frequency. The outside green (light/dark) circle represents the Tiwi variants mapped to gnomAD data's across distinct blood group genes. The red (light/dark) circle indicates variants with ISBT annotation and their frequency in gnomAD dataset. The blue (light/dark) circle indicates the distribution of Tiwi variants with no ISBT annotation and is rare in gnomAD dataset. The dark gray circle indicates the number of non‐ISBT variants annotated to the ClinVar database with gnomAD frequency. The yellow circle shows the distribution of the number of potential novel variants [Color figure can be viewed at wileyonlinelibrary.com]

Comparison of Tiwi blood group antigens with other indigenous populations

We compared the blood group antigen profiles of the Tiwi participants with those of other Indigenous Australian populations reported previously, , , , , , including Western Desert Indigenous Australians for whom a comprehensive blood group antigen profile prediction was performed using WES (Table 4). Overall, 10% of Tiwi samples in this dataset lacked the high‐frequency Yk(a+) antigen belonging to the KN blood group system as compared to 38% of the Western Desert population. Yk(a+) is common in most populations other than the African population. The Fy(a−b+) negative phenotype common among Caucasians, was predicted to be rare in the Tiwi population as well as in other Indigenous Australian populations. No Diego blood group antigen was detected among the Tiwi dataset except reference antigen; the same finding was observed in the western desert report. Interestingly, the rare KYO+ phenotype predicted in the Western Desert population, was absent among the Tiwi people. KYO+ was previously only been reported in two Japanese individuals. K+k− was completely missing in Tiwi but present in Caucasians (0.2%) (Table 4). The Mc+ antigen from the MNS blood group system was present in Tiwi but at the heterozygous level which was also observed in the western desert.

TABLE 4

Blood group antigen prevalence comparisons for the Tiwi and other previously reported populations

Blood group	Tiwi (%)	Western Desert (%) 9	Indigenous literature (%) 50 , 51 , 52 , 53 , 54 , 55	Caucasian (%) 27
Yk(a−)	10.07%	38%	Missing	8%
Fy(a−)	0.88	1.39	0%	34%
Di(a+)	0%	0%	0%	<0.1%
KYO+	Missing	<0.1	Missing	0%
K+k−	0%	0%	0%	0.2

DISCUSSION

Blood group antigens are poorly characterized for Indigenous Australians. the most recent serological blood typing for Indigenous populations was performed in the 1960s. Here, we compared the blood group genetic profile of the Indigenous Tiwi participants with those of African, American, Asian, and European populations from the 1KGP3 dataset and other previously published reports. Our population stratification analysis revealed the genetic makeup of the Tiwi participants to be distinct from that of other populations, as was that of the African participants, consistent with prior reports. We analyzed the distribution of the low‐prevalence phenotype among these six populations. The ABO phenotype A2 has been present primarily among European and African populations. These ABO subgroups exhibit distinct frequencies among different populations, thus requiring data from and experience with individual populations to ensure safe blood transfusion practices.

We identified the polymorphic and rare Au(a−b+) phenotype, belonging to the Lutheran blood group system, in all six populations with different frequencies among these groups. Although the Lutheran blood group system comprises 20 clinically well‐defined antigens, no specific data is available regarding the clinical significance of anti‐Au(a−b+). Our results predicted the highest prevalence of the weak Jk(a+wb−) phenotype, among the East Asian and Tiwi participants.

The most clinically significant blood group variants involved in HDFN are the RHD antigens, and other RH antigens including the c, E, e and C antigens. The anti‐D, and anti‐c have >50% risk of developing mild to severe HDFN if the fetus inherits the target paternal red cell antigen. Our copy number variation analysis for the homozygosity and null percentages of the RHD gene in all six populations indicated that the majority of the Tiwi (100%), East Asian (90.28%), African (96.52%), American (99.35%) and South Asian (99.48%) participants were homozygous for RHD gene (D+ phenotype). Only one Tiwi individual was hemizygous for the RHD gene. We also detected DIII type4 in eight Tiwi individuals at the homozygous level. DIII type4 is often mistyped as D+ and may lead to the risk of alloanti‐D production. , Additionally, maternal anti‐D with a partial D‐type leads to the risk of HDFN in a D‐positive fetus. We also report the prevalence of the rare D+C+c+E+e− phenotype in the Tiwi population, which was completely absent in 1KGP3 populations.

In the secondary analysis, blood group gene variants that are not previously assigned to known blood group antigens and are predicted to be deleterious using in silico methods were scanned for clinical significance, rarity, and potential novelty with global population frequency. We uncovered 13 variants that had pathogenic clinical association mostly involved in leukocyte adhesion deficiency, hemolytic anemia, and secretor/nonsecretor polymorphism. We also detected 29 rare variants with minimum allele frequencies of ≤0.05 in any of the five populations (African, American, East Asian, South Asian, and European) and 11 of them had two times higher frequency in Tiwi. Most interestingly, we identified 13 potential novel variants using in silico prediction techniques, of which three were common, and eight were rare in the Tiwi dataset.

Although one of the first comprehensive blood group studies of a Tiwi‐Indigenous Australian population, the current study has certain limitations. The RH copy number variation analysis was restricted to 151 samples due to the availability of only WGS data, and it did not include an investigation of hybrid RH and MNS systems. The predicted clinically significant and novel variant was done solely by applying the in‐silico tools method from RBCeq, further investigation is needed for their antigenic properties. Blood group calls could not be made for the XG, GATA1 and XK systems because the variants were called for only autosomal chromosomes. Additionally, serology testing and exon‐level detection were not performed to validate these RBCeq output data. Despite the limitations, our results have the potential to help identify donors from Indigenous Australian populations and to aid in blood group matching for patients with rare phenotypes. Our study also suggests that, in comparison to other ethnicities, East Asian populations have the next closest blood group profile to Tiwi. In addition to the associated anthropological interest, the clarification of well‐known blood group polymorphisms in diverse populations represents a key step towards ensuring the safety of blood transfusion. This work will provide critical data necessary for healthcare facilities to implement regional and national blood transfusion programs aimed at securing transfusion safety for the Tiwi population and other Indigenous Australian communities.

CONFLICT OF INTEREST

The authors declare that there is no conflict of interest.

ETHICS STATEMENT

The study entitled “To confirm the genetic regions found to cause renal disease in the Tiwi Island Aboriginal Population” received the full support of the Tiwi Island Land Council and was approved by the human research ethics committees of The Northern Territory Department of Health (2012–1767), The Australian National University (2014–663), The University of Queensland (2012001146), The University of Tasmania (H0012832), and Queensland University of Technology (1800000895). The blood group study was approved by Tiwi Island Land Council on 22nd June 2017 through personal communication to Dr Nagaraj.

Data S1 Supporting Information.

Click here for additional data file.