Skip to main content

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



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

Study design and methods

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


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


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

Peer Review reports


Data generated by next generation sequencing (NGS) are often utilized in the emerging fields of precision and personalized medicine. This massively parallel processing chemistry can identify genetic factors that predict treatment and response to therapies. Reference nucleotide sequences are critical for analyzing NGS data, as exemplified by routine clinical diagnosis for HLA antigens [1].

Genotype phasing is the process to determine if genetic variants, often single nucleotide variations, called SNVs, belong to 2 separate chromosomes (in trans). If SNVs are located on the same chromosome (in cis), they constitute a haplotype or an allele. Genotype phasing has often been inferred using computational methods [2, 3], which are prone to certain types of error [4]. These errors are encountered in samples harboring novel variants, low frequency or rare variants, and structural variants [5]. Almost all of these errors can be precluded by laboratory based methods, such as sequencing the genomes of both parents and sibling offspring [6], physical separation of homologous chromosomes in diploid cells [7, 8], sequencing in sperm cells [9], allele specific PCR [10], single DNA molecule dilution [11] and single molecule sequencing chemistry [12, 13]. These laboratory based methods are, however, labor-intensive and time consuming, and thus infrequently applied in clinical diagnostics.

The human genome contains many regions that are known as long contiguous stretches of homozygosity (LCSH) [14,15]. Their presence in unrelated individuals across different populations is attributed to a lower average recombination rate in these regions of the human genome [14].

The human atypical chemokine receptor 1 gene (ACKR1, MIM #613,665) [16] encodes a multi-pass trans-membrane glycoprotein. It is a receptor for pro-inflammatory cytokines [17] and malaria Plasmodium parasites (P. vivax and P. knowlesi) [18]. The ACKR1 glycoprotein carries the five antigens of the Duffy blood group system (Fy) [19, 20]. Recent sequencing studies in the ACKR1 gene have identified approximately 30 haplotypes, albeit at limited lengths of 2.1 kb [21], 2.5 kb [22], 5.2 kb [23], and 5.6 kb [24], respectively. We previously applied these ACKR1 haplotypes to predict the Duffy phenotype in Neanderthal samples [21]. Later, high-coverage genome sequences of Neanderthals were established [25,26,27], which confirmed our prediction [21]. A recent similar comparative study, involving long genomic segments, identified a 50 kb segment in humans, which was inherited from Neanderthals and represented a genetic risk factor in SARS-CoV-2 infection [28].

The 1000 Genomes Project (1000GP) provides a comprehensive database of genotypes and haplotypes in 2,504 unrelated individuals across 26 populations worldwide [29, 30]. As a proof of principle using data from the 1000GP for the ACKR1 gene, we establish a list of 902 haplotypes, some more than 80 kb long. Our scalable approach can be applied to any gene in any population.

Materials and methods

Algorithm workflow

A Python algorithm was developed (Supplementary Information, File S1) to download and analyze genotype data for 80.6 kb region of chromosome 1 (between positions NC_000001.11: 159,203,314–159,283,887) flanked between 2 genes, CADM3 and FCER1A, and encompassing the ACKR1 gene (Fig. 1) for all 2,504 unrelated individuals of the final release 1000GP panel (Phase 3; GRCh38) using Bcftools [31]. The SNV data was downloaded from the dbSNP database [32]. Individual sequences with heterozygosity at a single site or with complete homozygosity were automatically extracted as an unambiguous ACKR1 haplotype that can be considered experimentally confirmed, which applied a time-proven concept [4]. The algorithm outputs three files: a sequence file containing the distinct haplotypes, a meta-data file containing information about the population in which the haplotypes are found, and a folder containing graphical representations of the population distribution of the distinct haplotypes.

Fig. 1
figure 1

Schematic representation of chromosome 1 region analyzed. The ACKR1 gene is bordered by the 2 genes CADM3 in centromeric and FCER1A in telomeric direction at chromosomal position 1q23.2 (a). The structure of the ACKR1 gene (b) comprises 2 exons (closed boxes) and include the coding sequence (CDS,  black) and the 5’- and 3’-untranslated region (UTR, grey). The intron 1 joins the 2 exons (black line). The number of SNVs observed in the for the dbSNP (b) and 1000GP databases (c) are shown for the 5’-UTR, CDS, intron, CDS and 3’-UTR


Phased haplotype data for 80.6 kb region of chromosome 1 (between positions NC_000001.11: 159,203,314–159,283,887) was manually downloaded for all 2535 individuals of the 1000GP panel (Phase 3; GRCh37) from the 1000 Genomes browser. After removing 31 related individuals, haplotype data from 2504 unrelated individuals was imported into Microsoft Excel. Individuals with heterozygosity at a single site or with complete homozygosity in the 1,626 nucleotide-long ACKR1 gene (NG_011626.3; NC_000001.11:159,204,875–159,206,500) allowed unambiguous assignment of an ACKR1 haplotype. These unambiguous ACKR1 haplotypes were further analyzed individually using Excel spreadsheets, and their sequences were extended in both 5'- and 3'-directions until a heterozygous SNV was encountered. The region between 2 SNVs was catalogued as a haplotype and compared with the previous automated results. The manual analysis was performed and thus a validation dataset generated before the Python algorithm was developed.

Neanderthal genome

The published DNA sequence of the Neanderthal genome (Chagyrskaya, Altai, and Vindija 33.19, [2527] was analyzed (Integrative genomics viewer version 2.3.20) [33] and aligned to the human genome (NCBI Build GRCh38/hg38). We searched for the longest match, if any, with the haplotypes in the 1000GP.


Using the 1000GP database and a Python algorithm, we extracted and catalogued long haplotypes that encompassed the ACKR1 gene and were flanked between 2 SNVs (Fig. 1). Among 2,504 individuals included in the 1000GP database, 1,520 individuals were homozygous for the 1,626 nucleotide-long ACKR1 gene or heterozygous with only 1 SNV. The ACKR1 sequences for these individuals were further analyzed both upstream and downstream of ACKR1 gene until SNVs were encountered. The extension in both directions allowed us to identify long ACKR1 haplotypes that can be considered experimentally verified. The results obtained with our computational approach were validated by a manual method, performed in a blinded fashion.

ACKR1 and SNVs

For the ACKR1 gene (Fig. 1), the dbSNP database [32] lists 549 SNVs spread over 1,626 nucleotides (Fig. 1b). We encountered, however, only 43 SNVs of the ACKR1 gene in the 1000GP database (Fig. 1c) out of the 549 known SNVs.

ACKR1 haplotypes

We identified 31 distinct haplotypes with ≥ 10 observations (Table 1). They ranged in length from 2,383 nucleotides to 17,739 nucleotides. A total of 902 haplotypes were observed, ranging in length from 1,901 nucleotides to 80,584 nucleotides, some extending into the adjacent CADM3 and FCER1A genes (Fig. 2). The combined length of haplotype sequences comprised 19,895,388 nucleotides with a median of 16,014 nucleotides (Quartile 1 – Quartile 3: 7,588 – 30,729 nucleotides; Interquartile Range: 23,141 nucleotides). The length of the haplotypes was inversely proportional to the number of observations (Fig. 3). Most of the common haplotypes (70.13%) were small (< 10 kb; Table 2) and ranged in length between 1,901 to 9,927 nucleotides. The most common ACKR1 allele observed was the Duffy-null allele (FY*02 N.01) followed by FY*A (FY*01) and FY*B (FY*02), respectively (Table 3). For each of these 3 common ACKR1 alleles, we were able to identify reference sequences longer than 80 kb (Table 3).

Table 1 Experimentally confirmed ACKR1 haplotypes with ≥ 10 observations in the 1000GP database*
Fig. 2
figure 2

ACKR1 haplotypes observed in the 1000GP. A total of 902 unique haplotypes were observed and sorted according to their length (bars). All haplotypes comprise the ACKR1 gene (shaded column), their positions in the ACKR1 gene locus (top, see Fig. 1) is indicated. The cumulative number is listed (right). Haplotypes of similar lengths are grouped together (for exact lengths see Supplementary Information, Excel files S1 and S2)

Fig. 3
figure 3

Correlation between length and observations of ACKR1 haplotypes. The length of the ACKR1 haplotypes (x-axis) observed in the 1000GP was inversely proportional to the number of observations (y-axis)

Table 2 ACKR1 haplotypes and length distribution in the 1000GP database among 1520 individuals
Table 3 Length distribution of the 3 common ACKR1 alleles observed in the 1000GP

ACKR1 alleles in the Neanderthal samples

The 3 Neanderthal samples were GATA box negative (-67 T) and represented the ancestral FY*B allele (Table 4). None of the 3 Neanderthal ACKR1 sequences (Chagyrskaya, Altai, and Vindija 33.19) fully matched any of the 902 haplotypes. The 2 haplotypes closest to the Neanderthal sequences had 1 mismatch in the GATA box (Table 4).

Table 4 ACKR1 alleles in the 1000GP and 3 Neanderthal samples


In the current study, we identified 902 experimentally confirmed reference haplotypes for the ACKR1 gene, using only publicly available data from the large scale 1000GP study database. Our approach is easily scalable. It can be applied to similar databases, including the UK10K Consortium [34], the African Genome Variation Project [35] and the upcoming All of Us Research Program [36]. For proof of principle, we demonstrated the application using a Python algorithm for one gene. The approach can, however, define reference sequences for any segment of the genome, with genes or without.

We showed that reference sequences can be obtained from databases and verified without ambiguity at lengths exceeding 80 kb. Such reference sequences can be catalogued inexpensively for use in clinical diagnostics. The catalogue comprised the set of the longest unique haplotypes that can be distinguished by the gene’s nucleotide sequence. In clinical diagnostics with molecular-based assays, common and well documented (CWD) [37] reference haplotypes are routinely applied, for example in HLA typing [1]. Exact matching at the haplotype level improves survival following bone marrow transplantation [38] and reduces alloimmunization in chronically transfused patients [39,40,41]. A limited number of common haplotypes represented the majority in the population [42], and identifying haplotypes from databases is an economical way to obtain such reference sequences.

Apart from clinical diagnostics, long-range haplotypes are also useful to understand the influence of environment on positive selection of genes in human populations [43], for association mapping of genes that contribute to disease and other phenotypes [44], for correlating the geographical distribution of haplotypes with endemicity of disease [45], for identifying evolutionarily conserved elements and regulatory elements [46], and for improving the reliability of genotype imputation [47]. Long haplotypes identified by using SNV data from high-density oligonucleotide arrays and the International HapMap Project [48] have been shown to be population dependent and can provide important insights into human evolutionary history [49]. These studies may also identify regions of positive selection with important roles in human health and disease [50].

Next generation sequencing is increasingly used for blood group genes [51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78]. In contrast to HLA [79], most blood group genes lack well documented long reference sequences associated with them [80]. Hence, a comprehensive reference database for blood group genes will facilitate blood group genotyping by NGS. The Erythrogene database [59] contains the complete coding region sequence of many different blood group alleles obtained from the 1000GP. However, it lacks information for sequence variants in the non-coding regions, such as promoter, splice sites and long intronic regions, which can also affect the expression of antigens and helps to ascertain the allele and its coding sequence [81,82,83,84].

A large number of haplotypes were more than 50 kb long with some extending at least to 80.5 kb in length (Fig. 2). Our observations are consistent with previous reports suggesting that most of the human genome is contained in blocks of a few kb to more than 100 kb [85, 86]. However, most of the ACKR1 haplotypes in the 1000GP were small and concentrated closely around the ACKR1 gene. The number of haplotypes decreased as their length increased and extended into the intergenic regions (Fig. 3). This is explained because most of the variants in the dbSNP database resides in the intergenic regions [87].

Our 2 haplotypes HAP897 and HAP899 (Additional file 4:Table S3), observed once each in African populations, were closest to the 3 Neanderthal samples. Both haplotypes carried the GATA box mutation (c.-67C), which all Neanderthal samples lacked (c.-67T). Individuals homozygous for the GATA box mutation (c.-67C) do not express the Duffy glycoprotein on the red cell surface [81] making them resistant to invasion by the malarial parasite P. vivax [88,89,90]. This similarity in alleles, discrepant at nucleotide position c.-67 only, was consistent with the fact that the GATA box mutation (c.-67C) started to spread in Africa only around 30,000 years ago [91], while the 3 Neanderthals Vindija, Altai and Chagyrskaya are 50,000, 120,000 and 50,000 years old, respectively [25,26,27].

In clinical diagnostics for patients, long-range haplotypes harboring novel or rare SNVs can only be detected when the haplotype is sequenced at full-length [92]. Using Sanger sequencing, we have previously characterized the ERMAP [93], ICAM4 [94], and ACKR1 [23] blood group genes at the haplotype level and identified prevalent long-range reference alleles, a time consuming and low throughput approach. We showed in this study how long contiguous stretches of homozygosity (LCSH) can serve to generate a database of long haplotypes, as defined by full length nucleotide sequences rather than the concatenation of known SNVs. Relying on SNV data would miss patients carrying novel or rare alleles with possible clinical relevance, which are not identical to the reference sequences. Features of the 1000GP allowed us to catalogue these extended nucleotide sequences with population specific frequencies. Our approach will enable the positive identification of patients carrying these reference sequences.

We plan to extend this approach to all blood group systems recognized by the International Society of Blood Transfusion (ISBT) [95]. A tool under development will allow researchers the customized online extraction of long haplotypes from databases and genes or genomic regions of their choice. Eventually, our approach can be applied to any region of a chromosome. For now, the 902 ACKR1 alleles identified through our novel approach will be useful as templates for analyzing data from NGS, thus enhancing the reliability of clinical diagnostics.

Web Resources

1000 Genomes browser ( accessed on Aug 05, 2019. ISBT ( Max Planck Institute for Evolutionary Anthropology (

Availability of data and materials

The datasets analyzed and generated during the current study are available as supplementary tables and at 1000 Genomes browser (


  1. Robinson J, et al. The IPD and IMGT/HLA database: allele variant databases. Nucleic Acids Res. 2015;43:D423-431.

    Article  CAS  PubMed  Google Scholar 

  2. Halldórsson, B. V. et al. A survey of computational methods for determining haplotypes. In: Istrail S., Waterman M., Clark A. (eds) Computational methods for SNPs and haplotype inference. RSNPsH 2002. Lecture Notes in Computer Science. Springer, Berlin, Heidelberg. 2983, 26–47, (2004).

  3. Al Bkhetan Z, Zobel J, Kowalczyk A, Verspoor K, Goudey B. Exploring effective approaches for haplotype block phasing. BMC Bioinform. 2019;20:540.

    Article  Google Scholar 

  4. Clark AG. Inference of haplotypes from PCR-amplified samples of diploid populations. Mol Biol Evol. 1990;7:111–22.

    Article  CAS  PubMed  Google Scholar 

  5. Glusman G, Cox HC, Roach JC. Whole-genome haplotyping approaches and genomic medicine. Genome Med. 2014;6:73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Roach JC, et al. Analysis of genetic inheritance in a family quartet by whole-genome sequencing. Science. 2010;328:636–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ma L, et al. Direct determination of molecular haplotypes by chromosome microdissection. Nat Methods. 2010;7:299–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Yang H, Chen X, Wong WH. Completely phased genome sequencing through chromosome sorting. Proc Natl Acad Sci U S A. 2011;108:12–7.

    Article  PubMed  Google Scholar 

  9. Kirkness EF, et al. Sequencing of isolated sperm cells for direct haplotyping of a human genome. Genome Res. 2013;23:826–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Arbeithuber, B., Heissl, A. & Tiemann-Boege, I. in Haplotyping: Methods and Protocols (eds Irene Tiemann-Boege & Andrea Betancourt) 3–22 (Springer New York, 2017).

  11. Zheng GX, et al. Haplotyping germline and cancer genomes with high-throughput linked-read sequencing. Nat Biotechnol. 2016;34:303–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Rhoads A, Au KF. PacBio sequencing and its applications. Genom Proteom Bioinform. 2015;13:278–89.

    Article  Google Scholar 

  13. Jain M, Olsen HE, Paten B, Akeson M. The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community. Genome Biol. 2016;17:239.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Li LH, et al. Long contiguous stretches of homozygosity in the human genome. Hum Mutat. 2006;27:1115–21.

    Article  CAS  PubMed  Google Scholar 

  15. Gibson J, Morton NE, Collins A. Extended tracts of homozygosity in outbred human populations. Hum Mol Genet. 2006;15:789–95.

    Article  CAS  PubMed  Google Scholar 

  16. Nibbs RJB, Graham GJ. Immune regulation by atypical chemokine receptors. Nat Rev Immunol. 2013;13:815–29.

    Article  CAS  PubMed  Google Scholar 

  17. Horuk, R. The Duffy antigen receptor for chemokines DARC/ACKR1. Front Immunol 6, doi: (2015).

  18. Miller LH, Mason SJ, Dvorak JA, Mcginniss MH, Rothman IK. Erythrocyte receptors for (Plasmodium-Knowlesi) malaria - duffy blood-group determinants. Science. 1975;189:561–3.

    Article  CAS  PubMed  Google Scholar 

  19. Meny GM. The Duffy blood group system: a review. Immunohematology. 2010;26:51–6.

    Article  CAS  PubMed  Google Scholar 

  20. Meny GM. An update on the Duffy blood group system. Immunohematology. 2019;35:11–2.

    Article  PubMed  Google Scholar 

  21. Schmid P, Ravenell KR, Sheldon SL, Flegel WA. DARC alleles and Duffy phenotypes in African Americans. Transfusion. 2012;52:1260–7.

    Article  CAS  PubMed  Google Scholar 

  22. Fichou Y, et al. Defining blood group gene reference alleles by long-read sequencing: proof of concept in the ACKR1 gene encoding the duffy antigens. Transfusion Med Hemotherapy. 2020;47:23–32.

    Article  Google Scholar 

  23. Yin Q, Srivastava K, Gebremedhin A, Makuria AT, Flegel WA. Long-range haplotype analysis of the malaria parasite receptor gene ACKR1 in an East-African population. Hum Genome Var. 2018;5:26.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Srivastava K, et al. ACKR1 alleles at 5.6 kb in a well-characterized renewable US Food and Drug Administration (FDA) reference panel for standardization of blood group genotyping. J Mol Diagn. 2020;22:1272-1279. doi:

  25. Prüfer K, et al. The complete genome sequence of a Neanderthal from the Altai mountains. Nature. 2014;505:43–9.

    Article  CAS  PubMed  Google Scholar 

  26. Prüfer K, et al. A high-coverage Neandertal genome from Vindija Cave in Croatia. Science. 2017;358:655–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mafessoni F, et al. A high-coverage Neandertal genome from Chagyrskaya Cave. Proc Natl Acad Sci U S A. 2020;117:15132–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zeberg H, Pääbo S. The major genetic risk factor for severe COVID-19 is inherited from Neanderthals. Nature. 2020.

    Article  PubMed  Google Scholar 

  29. Genomes Project, C. et al. A global reference for human genetic variation. Nature. 2015;526:68–74.

  30. Sudmant PH, et al. An integrated map of structural variation in 2,504 human genomes. Nature. 2015;526:75–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics (Oxford, England). 2011;27:2987–93.

    Article  CAS  Google Scholar 

  32. Sherry ST, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Robinson JT, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Walter K, et al. The UK10K project identifies rare variants in health and disease. Nature. 2015;526:82–90.

    Article  CAS  PubMed  Google Scholar 

  35. Gurdasani D, et al. The African genome variation project shapes medical genetics in Africa. Nature. 2015;517:327–32.

    Article  CAS  PubMed  Google Scholar 

  36. Denny JC, et al. The “All of Us” research program. N Engl J Med. 2019;381:668–76.

    Article  PubMed  Google Scholar 

  37. Mack SJ, et al. Common and well-documented HLA alleles: 2012 update to the CWD catalogue. Tissue Antigens. 2013;81:194–203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Tay GK, et al. Matching for MHC haplotypes results in improved survival following unrelated bone marrow transplantation. Bone Marrow Transpl. 1995;15:381–5.

    CAS  Google Scholar 

  39. Chou ST, Liem RI, Thompson AA. Challenges of alloimmunization in patients with haemoglobinopathies. Br J Haematol. 2012;159:394–404.

    Article  CAS  PubMed  Google Scholar 

  40. Tournamille C, et al. Partial C antigen in sickle cell disease patients: clinical relevance and prevention of alloimmunization. Transfusion. 2010;50:13–9.

    Article  CAS  PubMed  Google Scholar 

  41. Allen ES, et al. Immunohaematological complications in patients with sickle cell disease after haemopoietic progenitor cell transplantation: a prospective, single-centre, observational study. Lancet Haematol. 2017;4:e553–61.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Slater N, et al. Power laws for heavy-tailed distributions: modeling allele and haplotype diversity for the national marrow donor program. PLoS Comput Biol. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Vallender EJ, Lahn BT. Positive selection on the human genome. Hum Mol Genet. 2004.

    Article  PubMed  Google Scholar 

  44. Gibson G, Muse SV. A primer of genome science. Sunderland, MA: Sinauer Associates; 2009.

    Google Scholar 

  45. Filosa S, et al. G6PD haplotypes spanning Xq28 from F8C to red/green color vision. Genomics. 1993;17:6–14.

    Article  CAS  PubMed  Google Scholar 

  46. Li MJ, Yan B, Sham PC, Wang J. Exploring the function of genetic variants in the non-coding genomic regions: approaches for identifying human regulatory variants affecting gene expression. Brief Bioinform. 2015;16:393–412.

    Article  CAS  PubMed  Google Scholar 

  47. Gudbjartsson DF, et al. Large-scale whole-genome sequencing of the Icelandic population. Nat Genet. 2015;47:435–44.

    Article  CAS  PubMed  Google Scholar 

  48. The International HapMap Project. Nature. 2003;426:789–96.

    Article  CAS  Google Scholar 

  49. Gusev A, et al. The architecture of long-range haplotypes shared within and across populations. Mol Biol Evol. 2012;29:473–86.

    Article  CAS  PubMed  Google Scholar 

  50. Zhang C, et al. A whole genome long-range haplotype (WGLRH) test for detecting imprints of positive selection in human populations. Bioinformatics (Oxford, England). 2006;22:2122–8.

    Article  CAS  Google Scholar 

  51. Stabentheiner S, et al. Overcoming methodical limits of standard RHD genotyping by next-generation sequencing. Vox Sang. 2011;100:381–8.

    Article  CAS  PubMed  Google Scholar 

  52. Rieneck K, et al. Next-generation sequencing: proof of concept for antenatal prediction of the fetal Kell blood group phenotype from cell-free fetal DNA in maternal plasma. Transfusion. 2013;53:2892–8.

    Article  CAS  PubMed  Google Scholar 

  53. Fichou Y, Audrézet MP, Guéguen P, Le Maréchal C, Férec C. Next-generation sequencing is a credible strategy for blood group genotyping. Br J Haematol. 2014;167:554–62.

    Article  CAS  PubMed  Google Scholar 

  54. Wieckhusen C, Bugert P. 454-sequencing for the KEL, JR, and LAN blood groups. Methods Mol Biol. 2015;1310:123–133. doi:

  55. Giollo M, et al. BOOGIE: predicting blood groups from high throughput sequencing data. PLoS ONE. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Lane WJ, et al. Comprehensive red blood cell and platelet antigen prediction from whole genome sequencing: proof of principle. Transfusion. 2016;56:743–54.

    Article  CAS  PubMed  Google Scholar 

  57. Lang K, et al. ABO allele-level frequency estimation based on population-scale genotyping by next generation sequencing. BMC Genomics. 2016;17:374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Fichou Y, Mariez M, Le Maréchal C, Férec C. The experience of extended blood group genotyping by next-generation sequencing (NGS): investigation of patients with sickle-cell disease. Vox Sang. 2016;111:418–24.

    Article  CAS  PubMed  Google Scholar 

  59. 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:240–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Baronas J, Westhoff C, Vege S, Mah H, Aguad M. RHD zygosity determination from whole genome sequencing data. J Blood Disord Transfus. 2016;7:1–5.

    Article  Google Scholar 

  61. Schoeman EM, 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:1078–88.

    Article  CAS  PubMed  Google Scholar 

  62. Dezan MR, et al. RHD and RHCE genotyping by next-generation sequencing is an effective strategy to identify molecular variants within sickle cell disease patients. Blood Cells Mol Dis. 2017;65:8–15.

    Article  CAS  PubMed  Google Scholar 

  63. Chou ST, et al. Whole-exome sequencing for RH genotyping and alloimmunization risk in children with sickle cell anemia. Blood Adv. 2017;1:1414–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Jakobsen MA, Dellgren C, Sheppard C, Yazer M, Sprogøe U. The use of next-generation sequencing for the determination of rare blood group genotypes. Transfus Med. 2019;29:162–8.

    Article  CAS  PubMed  Google Scholar 

  65. Schoeman EM, 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:284–93.

    Article  CAS  PubMed  Google Scholar 

  66. Orzińska A, et al. A preliminary evaluation of next-generation sequencing as a screening tool for targeted genotyping of erythrocyte and platelet antigens in blood donors. Blood Transf. 2018;16:285–292.

  67. Lane WJ, et al. Automated typing of red blood cell and platelet antigens: a whole-genome sequencing study. Lancet Haematol. 2018;5:e241–51.

  68. Wheeler MM, et al. Genomic characterization of the RH locus detects complex and novel structural variation in multi-ethnic cohorts. Genet Med. 2019;21:477–86.

    Article  CAS  PubMed  Google Scholar 

  69. Wu PC, et al. ABO genotyping with next-generation sequencing to resolve heterogeneity in donors with serology discrepancies. Transfusion. 2018;58:2232–42.

    Article  CAS  PubMed  Google Scholar 

  70. Montemayor-Garcia C, et al. Genomic coordinates and continental distribution of 120 blood group variants reported by the 1000 Genomes Project. Transfusion. 2018;58:2693–704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Tounsi WA, Madgett TE, Avent ND. Complete RHD next-generation sequencing: establishment of reference RHD alleles. Blood Adv. 2018;2:2713–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Schoeman EM, Roulis EV, Perry MA, Flower RL, Hyland CA. Comprehensive blood group antigen profile predictions for Western Desert Indigenous Australians from whole exome sequence data. Transfusion. 2019;59:768–78.

    Article  CAS  PubMed  Google Scholar 

  73. Orzińska A, et al. Prediction of fetal blood group and platelet antigens from maternal plasma using next-generation sequencing. Transfusion. 2019;59:1102–7.

    Article  CAS  PubMed  Google Scholar 

  74. Lane WJ, et al. Automated typing of red blood cell and platelet antigens from whole exome sequences. Transfusion. 2019;59:3253–63.

    Article  CAS  PubMed  Google Scholar 

  75. Halls JBL, et al. Overcoming the challenges of interpreting complex and uncommon RH alleles from whole genomes. Vox Sang. 2020.

    Article  PubMed  Google Scholar 

  76. Fürst D, et al. Next-generation sequencing technologies in blood group typing. Transf Med Hemother. 2020;47:4–13.

  77. Wu PC, Pai S-C, Chen P-L. Blood group genotyping goes next generation: featuring ABO, RH and MNS. ISBT Sci Ser. 2018;13:290–7.

    Article  Google Scholar 

  78. Orzinska A, Guz K, Brojer E. Potential of next-generation sequencing to match blood group antigens for transfusion. Int J Clin Transfus Med. 2019;7:11–22.

    Article  Google Scholar 

  79. Barone JC, et al. HLA-genotyping of clinical specimens using Ion Torrent-based NGS. Hum Immunol. 2015;76:903–9.

    Article  CAS  PubMed  Google Scholar 

  80. Reid ME. Transfusion in the age of molecular diagnostics. Hematol Am Soc Hematol Educ Program. 2009;2009:171–7.

  81. Tournamille C, Colin Y, Cartron JP, Le Van Kim C. Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy-negative individuals. Nat Genet. 1995;10;224–8.

  82. Lucien N, et al. Characterization of the gene encoding the human Kidd blood group/urea transporter protein. Evidence for splice site mutations in Jknull individuals. J Biol Chem. 1998;273:12973–80.

  83. Lomas-Francis C, Reid ME. The Dombrock blood group system: a review. Immunohematology. 2010;26:71–8.

    Article  CAS  PubMed  Google Scholar 

  84. Christophersen MK, et al. SMIM1 variants rs1175550 and rs143702418 independently modulate Vel blood group antigen expression. Sci Rep. 2017;7:40451.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Gabriel SB, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–9.

    Article  CAS  PubMed  Google Scholar 

  86. Wall JD, Pritchard JK. Haplotype blocks and linkage disequilibrium in the human genome. Nat Rev Genet. 2003;4:587–97.

    Article  CAS  PubMed  Google Scholar 

  87. Jin Y, Wang J, Bachtiar M, Chong SS, Lee CGL. Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes. Hum Genomics. 2018;12:43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Miller LH, Mason SJ, Clyde DF, McGinniss MH. The resistance factor to Plasmodium vivax in blacks. The Duffy-blood-group genotype, FyFy. N Engl J Med. 1976;295:302–304.

  89. Chaudhuri A, et al. Purification and characterization of an erythrocyte membrane protein complex carrying Duffy blood group antigenicity. Possible receptor for Plasmodium vivax and Plasmodium knowlesi malaria parasite. J Biol Chem. 1989;264:13770–13774.

  90. Hadley TJ, Peiper SC. From malaria to chemokine receptor: the emerging physiologic role of the Duffy blood group antigen. Blood. 1997;89:3077–91.

    Article  CAS  PubMed  Google Scholar 

  91. Hamblin MT, Di Rienzo A. Detection of the signature of natural selection in humans: evidence from the Duffy blood group locus. Am J Hum Genet. 2000;66:1669–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Suk EK, et al. A comprehensively molecular haplotype-resolved genome of a European individual. Genome Res. 2011;21:1672–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Srivastava K, Lee E, Owens E, Rujirojindakul P, Flegel WA. Full-length nucleotide sequence of ERMAP alleles encoding Scianna (SC) antigens. Transfusion. 2016;56:3047–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Yin Q, et al. Molecular analysis of the ICAM4 gene in an autochthonous East African population. Transfusion. 2019;59:1880–1.

    Article  PubMed  PubMed Central  Google Scholar 

  95. (2020).

Download references


Bo Lan participated in the study during his Summer Internship Program at NIH in 2019.


Open Access funding provided by the National Institutes of Health (NIH). This work was supported in part by the Intramural Research Program (projects ZIC CL002128 and  RASCL#727301) of the NIH Clinical Center at the National Institutes of Health. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



WAF and KS conceived the study; KS designed the analysis and downloaded the data; ASF programmed the algorithm; WAF, KS, ASF and BL analyzed the data; WAF, ASF and KS wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Willy Albert Flegel.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors consent to the publication of this manuscript.

Competing interests


Statement of disclaimer

The views expressed do not necessarily represent the view of the National Institutes of Health, the Department of Health and Human Services, or the U.S. Federal Government.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

File S1. Python algorithm.

Additional file 2

Table S1. Populations in the 1000GP database.

Additional file 3

Table S2. Sequence data for the 902 long range ACKR1 haplotypes in the 1000GP.

Additional file 4

Table S3. Metadata file for the ACKR1 long range haplotypes in the 1000GP.

Additional file 5

Table S4. Exonic SNV distribution in the 902 experimentally confirmed ACKR1 haplotypes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Srivastava, K., Fratzscher, AS., Lan, B. et al. Cataloguing experimentally confirmed 80.7 kb-long ACKR1 haplotypes from the 1000 Genomes Project database. BMC Bioinformatics 22, 273 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: