Evaluation of targeted exome sequencing for 28 protein-based blood group systems, including the homologous gene systems, for blood group genotyping

Detecting serologically difficult ABO blood groups using single-molecule real-time sequencing technology

BACKGROUND AND OBJECTIVES

Recently, third-generation long-read sequencing technology has been increasingly applied to the detection of various blood group systems. Because of its long read length and use of single-molecule sequencing, it is capable of obtaining the sequences of blood group genes in their entirety as well as of distinguishing haplotypes. Therefore, here, we collected ABO blood group samples that were difficult to classify serologically and analysed the sequences of the coding regions of the ABO genes as well as the sequences upstream and downstream of the coding regions.

MATERIALS AND METHODS

Samples with ABO antigen typing and reverse serum typing discrepancies were screened in a total of 21 patients. All samples were subjected to serological testing and preliminary ABO genotyping (polymerase chain reaction with sequence-specific primers [PCR-SSP]), followed by single-molecule real-time (SMRT) sequencing to obtain complete ABO gene sequences. PCR sequence-based typing (PCR-SBT) was performed to validate the results.

RESULTS

Of the 21 samples, 15 had common ABO types, and 6 had rare ABO subtypes. One new allele, ABO*B.NEW (c.861C>T), and one allelic base recombination event was identified. Forty-two haplotype sequences were obtained via SMRT sequencing with intronic single-nucleotide variants (SNVs) specific to the ABO allele, and all of the exon region sequences were consistent with the PCR-SBT results.

CONCLUSION

SMRT sequencing is capable of accurately obtaining complete ABO gene sequences, distinguishing haplotypes and identifying allelic recombination.

A Bioinformatically Initiated Approach to Evaluate GATA1 Regulatory Regions in Samples with Weak D, Del, or D- Phenotypes Despite Normal RHD Exons

Introduction

With over 360 blood group antigens in systems recognized, there are antigens, such as RhD, which demonstrate a quantitative reduction in antigen expression due to nucleotide variants in the non-coding region of the gene that result in aberrant splicing or a regulatory mechanism. This study aimed to evaluate bioinformatically predicted GATA1-binding regulatory motifs in the RHD gene for samples presenting with weak or apparently negative RhD antigen expression but showing normal RHD exons.

Methods

Publicly available open chromatin region data were overlayed with GATA1 motif candidates in RHD. Genomic DNA from weak D, Del or D- samples with normal RHD exons (n = 13) was used to confirm RHD zygosity by quantitative PCR. Then, RHD promoter, intron 1, and intron 2 regions were amplified for Sanger sequencing to detect potential disruptions in the GATA1 motif candidates. Electrophoretic mobility shift assay (EMSA) was performed to assess GATA1-binding. Luciferase assays were used to assess transcriptional activity.

Results

Bioinformatic analysis identified five of six GATA1 motif candidates in the promoter, intron 1 and intron 2 for investigation in the samples. Luciferase assays showed an enhancement in transcription for GATA1 motifs in intron 1 and for intron 2 only when the R 2 haplotype variant (rs675072G>A) was present. GATA1 motifs were intact in 12 of 13 samples. For one sample with a Del phenotype, a novel RHD c.1-110A>C variant disrupted the GATA1 motif in the promoter which was supported by a lack of a GATA1 supershift in the EMSA and 73% transcriptional activity in the luciferase assay. Two samples were D+/D- chimeras.

Conclusion

The bioinformatic predictions enabled the identification of a novel DEL allele, RHD c.1-110A>C, which disrupted the GATA1 motif in the proximal promoter. Although the majority of the samples investigated here remain unexplained, we provide GATA1 targets which may benefit future RHD regulatory investigations.

Simultaneous high throughput genotyping of 36 blood group systems using NGS based on probe capture technology

Human blood group antigen has important biological functions, and transfusion of incompatible blood can cause alloimmunization and may lead to serious hemolytic reactions. Currently, serological methods are most commonly used in blood group typing. However, this technique has certain limitations and cannot fully meet the increasing demand for the identification of blood group antigens. This study describes a next-generation sequencing (NGS) technology platform based on exon and flanking region capture probes to detect full coding exon and flanking intron regions of the 36 blood group systems, providing a new high-throughput method for the identification of blood group antigens. The 871 capture probes were designed for the exon and flanking intron sequences of 36 blood group system genes, and synchronization analysis for 36 blood groups was developed. The library for NGS was tested using the MiSeq Sequencing Reagent Kit (v2, 300 cycles) by Illumina NovaSeq, and the data were analyzed by the CLC Genomics Workbench 21.0 software. A total of 199 blood specimens have been sequenced for the 41 genes from 36 blood groups. Among them, heterozygote genotypes were found in the ABO, Rh, MNS, Lewis, Duffy, Kidd, Diego, Gerbich, Dombrock, Globoside, JR, LAN, and Landsteiner-Wiene blood group systems. Only the homozygous genotype was found in the remaining 22 blood group systems. The obtained data in the NGS method shows a good correlation (99.98 %) with those of the polymerase chain reaction-sequence-based typing. An NGS technology platform for 36 blood group systems genotyping was successfully established, which has the characteristics of high accuracy, high throughput, and wide coverage.

A cross-sectional pilot study to estimate the frequency of minor blood group alleles and phenotypes in RhD-negative North Indian blood donors by DNA microarray analysis

INTRODUCTION

There are scarce data on Indian blood donors with respect to blood group phenotypes using molecular diagnostic modalities. Hence, we planned to estimate frequencies of blood group alleles/phenotypes using DNA microarray analysis in the north Indian RhD-negative blood donor population. With this initial pilot study, we plan to expand it to our entire donor population.

METHODOLOGY

The cross-sectional prospective study was conducted on 50 Indian blood donors (O RhD negative), to study the blood group genotype frequency. Genotyping for the most relevant red blood cell antigens (Rh, Kell, Duffy, Kidd, MNS, Lutheran, and Dombrock) was done using Bioarray Precise TypeHM Human Erythrocyte Antigen BeadChip kit containing probes directed to polymorphic sites.

RESULTS

In the Rh system, the most common alleles were RHCE*e/RHCE*e (98%) and RHCE*c/RHCE*c (80%). Phenotype K-k+ (genotype- KEL*02/KEL*02) was seen in 98% of samples, Js(a-b+) (KEL*02.07/KEL*02.07) was detected in 98% (49/50) of the samples tested. Jk(a + b+) (JK*01/JK*02) was the most common phenotype (48%) in the Kidd blood group system. In MNSs system, M+N+ (GYPA*01/GYPA*02) 44% and S+s+U+ (GYPB*03/GYPB*04) 34% were the most common phenotypes detected.

CONCLUSION

This pilot study shows the feasibility of genotyping a Northern Indian donor population. To the best of our knowledge, it is the first study on molecular blood grouping in Indian blood donors using the Bioarray platform.

Machine learning to optimize automated RH genotyping using whole-exome sequencing data

ABSTRACT

Rh phenotype matching reduces but does not eliminate alloimmunization in patients with sickle cell disease (SCD) due to RH genetic diversity that is not distinguishable by serological typing. RH genotype matching can potentially mitigate Rh alloimmunization but comprehensive and accessible genotyping methods are needed. We developed RHtyper as an automated algorithm to predict RH genotypes using whole-genome sequencing (WGS) data with high accuracy. Here, we adapted RHtyper for whole-exome sequencing (WES) data, which are more affordable but challenged by uneven sequencing coverage and exacerbated sequencing read misalignment, resulting in uncertain predictions for (1) RHD zygosity and hybrid alleles, (2) RHCE∗C vs. RHCE∗c alleles, (3) RHD c.1136C>T zygosity, and (4) RHCE c.48G>C zygosity. We optimized RHtyper to accurately predict RHD and RHCE genotypes using WES data by leveraging machine learning models and improved the concordance of WES with WGS predictions from 90.8% to 97.2% for RHD and 96.3% to 98.2% for RHCE among 396 patients in the Sickle Cell Clinical Research and Intervention Program. In a second validation cohort of 3030 cancer survivors (15.2% Black or African Americans) from the St. Jude Lifetime Cohort Study, the optimized RHtyper reached concordance rates between WES and WGS predications to 96.3% for RHD and 94.6% for RHCE. Machine learning improved the accuracy of RH predication using WES data. RHtyper has the potential, once implemented, to provide a precision medicine-based approach to facilitate RH genotype-matched transfusion and improve transfusion safety for patients with SCD. This study used data from clinical trials registered at ClinicalTrials.gov as #NCT02098863 and NCT00760656.

When and why is red blood cell genotyping applicable in transfusion medicine: a systematic review of the literature

This review aims to provide a better understanding of when and why red blood cell (RBC) genotyping is applicable in transfusion medicine. Articles published within the last 8 years in peer-reviewed journals were reviewed in a systematic manner. RBC genotyping has many applications in transfusion medicine including predicting a patient's antigen profile when serologic methods cannot be used, such as in a recently transfused patient, in the presence of autoantibody, or when serologic reagents are not available. RBC genotyping is used in prenatal care to determine zygosity and guide the administration of Rh immune globulin in pregnant women to prevent hemolytic disease of the fetus and newborn. In donor testing, RBC genotyping is used for resolving ABO/D discrepancies for better donor retention or for identifying donors negative for high-prevalence antigens to increase blood availability and compatibility for patients requiring rare blood. RBC genotyping is helpful to immunohematology reference laboratory staff performing complex antibody workups and is recommended for determining the antigen profiles of patients and prospective donors for accurate matching for C, E, and K in multiply transfused patients. Such testing is also used to determine patients or donors with variant alleles in the Rh blood group system. Information from this testing aides in complex antibody identification as well as sourcing rare allele-matched RBC units. While RBC genotyping is useful in transfusion medicine, there are limitations to its implementation in transfusion services, including test availability, turn-around time, and cost.

A cold case of hemolytic disease of the fetus and newborn resolved by genomic sequencing and population studies to define a new antigen in the Rh system

BACKGROUND

We report an obstetric case involving an RhD-positive woman who had developed a red blood cell (RBC) antibody that was not detected until after delivery of a newborn, who presented with a positive direct antiglobulin test result. Immunohematology studies suggested that the maternal antibody was directed against a low-prevalence antigen on the paternal and newborn RBCs.

RESULTS

Comprehensive blood group profiling by targeted exome sequencing revealed a novel nonsynonymous single nucleotide variant (SNV) RHCE c.486C>G (GenBank MZ326705) on the RHCE*Ce allele, for both the father and newborn. A subsequent genomic-based study to profile blood groups in an Indigenous Australian population revealed the same SNV in 2 of 247 individuals. Serology testing showed that the maternal antibody reacted specifically with RBCs from these two individuals.

DISCUSSION

The maternal antibody was directed against a novel antigen in the Rh blood group system arising from an RHCE c.486C>G variant on the RHCE*Ce allele linked to RHD*01. The variant predicts a p.Asn162Lys change on the RhCE protein and has been registered as the 56th antigen in the Rh system, ISBT RH 004063.

CONCLUSION

This antibody was of clinical significance, resulting in a mild to moderate hemolytic disease of the fetus and newborn (HDFN). In the past, the cause of such HDFN cases may have remained unresolved. Genomic sequencing combined with population studies now assists in resolving such cases. Further population studies have potential to inform the need to design population-specific red cell antibody typing panels for antibody screening in the Australian population.

Feasibility for non-invasive prenatal fetal blood group and platelet genotyping by massively parallel sequencing: A single test system for multiple atypical red cell, platelet and quality control markers

Non-invasive prenatal tests (NIPT) to predict fetal red cell or platelet antigen status for alloimmunised women are provided for select antigens. This study reports on massively parallel sequencing (MPS) using a red cell and platelet probe panel targeting multiple nucleotide variants, plus individual identification single nucleotide polymorphisms (IISNPs). Maternal blood samples were provided from 33 alloimmunised cases, including seven with two red cell antibodies. Cell-free and genomic DNA was sequenced using targeted MPS and bioinformatically analysed using low-frequency variant detection. The resulting maternal genomic DNA allele frequency was subtracted from the cell-free DNA counterpart. Outcomes were matched against validated phenotyping/genotyping methods, where available. A 2.5% subtractive allele frequency threshold was set after comparing MPS predictions for K, RhC/c, RhE/e and Fya /Fyb against expected outcomes. This threshold was used for subsequent predictions, including HPA-15a, Jka /Jkb , Kpa /Kpb and Lua . MPS outcomes were 97.2% concordant with validated methods; one RhC case was discordantly negative and lacked IISNPs. IISNPs were informative for 30/33 cases as controls. NIPT MPS is feasible for fetal blood group genotyping and covers multiple blood groups and control targets in a single test. Noting caution for the Rh system, this has the potential to provide a personalised service for alloimmunised women.

Prevalence of GP. Mur variant phenotype among Malaysian blood donors

BACKGROUND AND OBJECTIVE

A number of glycophorin variant phenotypes or hybrid glycophorin variants of the MNS blood group system bear multiple immunogenic antigens such as Mia, Mur, and MUT. In the East and Southeast Asian populations, glycoprotein (GP.) Mur is the most common glycophorin variant phenotype expressing those three immunogens. The aim of this study was to detect MNS system glycophorin variant phenotypes (GP. Mur, GP. Hop, GP. Bun, GP. HF, and GP. Hut) among Malaysian blood donors.

MATERIALS AND METHODS

In this cross-sectional study, 144 blood donors were selected under stratified random sampling. The deoxyribonucleic acid was extracted from whole blood samples, followed by a polymerase chain reaction assay. Sanger sequencing was used to identify the specific MNS variants and then validated by a serological crossmatch with known anti-Mur and anti-MUT.

RESULTS

GP. Mur was identified among Malaysian blood donors with a prevalence of 6.94%, and no other variants of the MNS system were found.

CONCLUSION

The present study substantiates that GP. Mur is the main variant of the MNS system glycophorin (B-A-B) hybrid in Malaysian blood donors. GP. Mur-negative red blood cells must therefore be considered in the current transfusion policy in order to prevent alloimmunization and immune-mediated transfusion reactions, particularly in transfusion-dependent patients.

Hemolytic disease of the fetus and newborn caused by anti-sD antibody in a GP.Mur/Mur Thai mother and review of the prevalence of sD in Thai blood donors

BACKGROUND

Low-prevalence antigen s^D (MNS23) is encoded by GYPB c.173C > G. Hemolytic disease of the fetus and newborn (HDFN) due to anti-s^D is rare. A mother delivered a newborn whose red blood cells (RBCs) were DAT-positive and was later diagnosed with HDFN. Serum from the mother was incompatible with the father's RBCs and was used to screen 184 Thai blood donors. This study aimed to investigate the cause of HDFN in a Thai family and determine the prevalence of s^D in Thai blood donors.

MATERIALS AND METHODS

Three family members and four blood donors were investigated in the study. Massively Parallel Sequencing (MPS) was used for genotyping. Standard hemagglutination techniques were used in titration studies, phenotyping, and enzyme/chemical studies. Anti-s, anti-Mi^a , anti-JENU, and anti-s^D reagents were used in serological investigations.

RESULTS

The mother was GYP*Mur/Mur. The father and the four donors were GYPB*s/s^D predicting S - s + s^D +. The baby was GYP*Mur/sD and his RBCs were Mi^a +, s + ^w with anti-s (P3BER) and JENU+^w . RBCs from two GYPB*s^D -positive blood donors reacted with anti-s^D (Dreyer). Proteolytic enzyme α-chymotrypsin-treated s^D + cells did not react with anti-s^D (Wat) produced by the GP.Mur/Mur mother but reacted with the original anti-s^D (Dreyer).

DISCUSSION

This is the first report of HDFN due to anti-s^D in the Asian population. The genotype frequency for GYPB*s^D in a selected Thai blood donor population is 2.2% (4/184). Anti-s^D should be considered in mothers with Southeast Asian or East Asian background when antibody identification is unresolved in pregnancies affected by HDFN.

Prediction of various blood group systems using Korean whole-genome sequencing data

Aims

This study established blood group analysis methods using whole-genome sequencing (WGS) data and conducted blood group analyses to determine the domestic allele frequency using public data from the Korean whole sequence analysis of the Korean Reference Genome Project conducted by the Korea Disease Control and Prevention Agency (KDCA).

Materials and methods

We analyzed the differences between the human reference sequences (hg19) and the conventional reference cDNA sequences of blood group genes using the Clustal Omega website, and established blood group analysis methods using WGS data for 41 genes, including 39 blood group genes involved in 36 blood group antigens, as well as the GATA1 and KLF1 genes, which are erythrocyte-specific transcription factor genes. Using CLC genomics Workbench 11.0 (Qiagen, Aarhus, Denmark), variant analysis was performed on these 41 genes in 250 Korean WGS data sets, and each blood group’s genotype was predicted. The frequencies for major alleles were also investigated and compared with data from the Korean rare blood program (KRBP) and the Erythrogene database (East Asian and all races).

Results

Among the 41 blood group-related genes, hg19 showed variants in the following genes compared to the conventional reference cDNA: GYPA, RHD, RHCE, FUT3, ACKR1, SLC14A1, ART4, CR1, and GCNT2. Among 250 WGS data sets from the Korean Reference Genome Project, 70.6 variants were analyzed in 205 samples; 45 data samples were excluded due to having no variants. In particular, the FUT3, GNCT2, B3GALNT1, CR1, and ACHE genes contained numerous variants, with averages of 21.1, 13.9, 13.4, 9.6, and 7.0, respectively. Except for some blood groups, such as ABO and Lewis, for which it was difficult to predict the alleles using only WGS data, most alleles were successfully predicted in most blood groups. A comparison of allele frequencies showed no significant differences compared to the KRBP data, but there were differences compared to the Erythrogene data for the Lutheran, Kell, Duffy, Yt, Scianna, Landsteiner-Wiener, and Cromer blood group systems. Numerous minor blood group systems that were not available in the KRBP data were also included in this study.

Conclusions

We successfully established and performed blood group analysis using Korean public WGS data. It is expected that blood group analysis using WGS data will be performed more frequently in the future and will contribute to domestic data on blood group allele frequency and eventually the supply of safe blood products.

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

There have been no comprehensive studies of a full range of blood group polymorphisms within the Australian population. This 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 vs global data. Whole-genome sequencing data were analyzed for 2796 healthy older Australians from the Medical Genome Reference Bank and compared with data from 1000 Genomes phase 3 (1KGP3) databases comprising 661 African, 347 American, 503 European, 504 East Asian, and 489 South Asian participants. There were 661 rare variants detected in this Australian sample population, including 9 variants that had clinical associations. Notably, we identified 80 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, 2 novel and 15 rare variants were found. Our detailed blood group profiling results provide a starting point for the creation of an Australian blood group variant database.

A new high-prevalence LW antigen detected by an antibody in an Indigenous Australian homozygous for LW*A c.309C>A variant

BACKGROUND AND OBJECTIVES

The LW gene encodes the LW glycoprotein that carries the antigens of the LW blood group system. LW antigens are distinct from D antigen, however, they are phenotypically related and anti-LW antibodies are often mistaken as anti-D. An antibody was detected in an Australian patient of Aboriginal descent who consistently typed as LW(a+b-). This study aimed to describe the antibody recognizing a high-prevalence antigen on the LW glycoprotein.

STUDY DESIGN AND METHODS

Samples from the patient and her four siblings were investigated. DNA was genotyped by single nucleotide polymorphism (SNP)-microarray and massively parallel sequencing (MPS) platforms. Red blood cells (RBCs) were phenotyped using standard haemagglutination techniques. Antibody investigations were performed using a panel of phenotyped RBCs from adults and cord blood cells.

RESULTS

SNP-microarray and MPS genotyped all family members as LW*A/A, (c.299A), predicting LW(a+b-). In addition, a novel LW*A c.309C>A single nucleotide variant was detected in all family members. The patient and one of her siblings (M4) were LW c.309C>A homozygous. Antibody from the patient reacted positive to all reagent panel RBCs and cord blood cells but negative with RBCs from LW(a-b-), Rh_null and sibling M4. Antibody failed to react with RBCs treated with dithiothreitol.

CONCLUSION

Antibody detected in the patient recognized a novel high-prevalence antigen, LWEM, in the LW blood group system. LWEM-negative patients who developed anti-LWEM can be safely transfused with D+ RBCs, however, D- is preferred. Accurate antibody identification can help better manage allocation of blood products especially when D- RBCs are in short supply.

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

Background

Methods and materials

Results

Conclusion

Application of Blood Group Genotyping by Next-Generation Sequencing in Various Immunohaematology Cases

BACKGROUND

Next-generation sequencing (NGS) technology has been recently introduced into blood group genotyping; however, there are few studies using NGS-based blood group genotyping in real-world clinical settings. In this study, we applied NGS-based blood group genotyping into various immunohaematology cases encountered in routine clinical practice.

METHODS

This study included 4 immunohaematology cases: ABO subgroup, ABO chimerism, antibody to a high-frequency antigen (HFA), and anti-CD47 interference. We designed a hybridization capture-based NGS panel targeting 39 blood group-related genes and applied it to the 4 cases.

RESULTS

NGS analysis revealed a novel intronic variant (NM_020469.3:c.29-10T>G) in a patient with an Ael phenotype and detected a small fraction of ABO*A1.02 (approximately 3-6%) coexisting with the major genotype ABO*B.01/O.01.02 in dizygotic twins. In addition, NGS analysis found a homozygous stop-gain variant (NM_004827.3:c.376C>T, p.Gln126*; ABCG2*01N.01) in a patient with an antibody to an HFA; consequently, this patient's phenotype was predicted as Jr(a-). Lastly, blood group phenotypes predicted by NGS were concordant with those determined by serology in 2 patients treated with anti-CD47 drugs.

CONCLUSION

NGS-based blood group genotyping can be used for identifying ABO subgroup alleles, low levels of blood group chimerism, and antibodies to HFAs. Furthermore, it can be applied to extended blood group antigen matching for patients treated with anti-CD47 drugs.

RBCeq: A robust and scalable algorithm for accurate genetic blood typing

Background

While blood transfusion is an essential cornerstone of hematological care, patients requiring repetitive transfusion remain at persistent risk of alloimmunization due to the diversity of human blood group polymorphisms. Despite the promise, user friendly methods to accurately identify blood types from next-generation sequencing data are currently lacking. To address this unmet need, we have developed RBCeq, a novel genetic blood typing algorithm to accurately identify 36 blood group systems.

Methods

RBCeq can predict complex blood groups such as RH, and ABO that require identification of small indels and copy number variants. RBCeq also reports clinically significant, rare, and novel variants with potential clinical relevance that may lead to the identification of novel blood group alleles.

Findings

The RBCeq algorithm demonstrated 99·07% concordance when validated on 402 samples which included 29 antigens with serology and 9 antigens with SNP-array validation in 14 blood group systems and 59 antigens validation on manual predicted phenotype from variant call files. We have also developed a user-friendly web server that generates detailed blood typing reports with advanced visualization (https://www.rbceq.org/).

Interpretation

RBCeq will assist blood banks and immunohematology laboratories by overcoming existing methodological limitations like scalability, reproducibility, and accuracy when genotyping and phenotyping in multi-ethnic populations. This Amazon Web Services (AWS) cloud based platform has the potential to reduce pre-transfusion testing time and to increase sample processing throughput, ultimately improving quality of patient care.

Funding

This work was supported in part by Advance Queensland Research Fellowship, MRFF Genomics Health Futures Mission (76,757), and the Australian Red Cross LifeBlood. The Australian governments fund the Australian Red Cross Lifeblood for the provision of blood, blood products and services to the Australian community.

Molecular genetic analysis of Mia -positive hybrid glycophorins revealed two novel alleles of GP.Vw and multiple variant transcripts of GYPB existing in both the homozygous GP.Mur and wild-type GPB individuals

BACKGROUND

The hybrid glycophorins of MNS blood group system express a series of low incidence antigens including Mi^a , which are commonly found in Southeast Asian populations. In this study, the molecular basis of Mi^a -positive hybrid glycophorins was firstly clarified in the Chinese Southern Han population. RNA transcripts of GYPB gene in the homozygous GP.Mur individuals were also analyzed.

STUDY DESIGN AND METHODS

DNAs were extracted from the whole blood samples of 111 Mi^a -positive donors. Then, high-resolution melting (HRM) analysis for GYP(B-A-B) was used to analyze the genotypes. Sequencing of GYPB pseudoexon 3 was conducted in the samples with variant melting curves. TA-cloning and subsequent sequencing of GYPA exons 2-4 were performed in the Mi^a -positive samples with normal GYPB/GYPB genotype by HRM. The transcript analysis of GYPB was conducted in homozygous GP.Mur and wild-type glycophorin B (GPB) individuals using RNA extracted from the cultured erythroblast.

RESULTS

The heterozygous GYP*Mur/GYPB (n = 101), homozygous GYP*Mur/GYP*Mur (n = 7) including one novel GYP*Mur allele with an extra GYPA/GYPE specific nucleotide substitution (c.229+110A>T), heterozygous GYP*Bun/GYPB (n = 1) and GYP*Vw/GYPA (n = 2) with two novel GYP*Vw alleles were identified. RNA transcript analysis revealed multiple transcripts of GYPB existing in both homozygous GP.Mur and normal GPB individuals.

CONCLUSION

The results showed the genetic diversity of hybrid glycophorins in the Chinese population. Besides, the successful analysis of GYPB transcripts indicates that the cultured erythroblast is a good source for RNA transcript analysis for the protein only expressed on the red blood cells.

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

Background

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.

Results

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.

Conclusions

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04169-6.

A novel algorithm comprehensively characterizes human RH genes using whole-genome sequencing data

RHD and RHCE genes encode Rh blood group antigens and exhibit extensive single-nucleotide polymorphisms and chromosome structural changes in patients with sickle cell disease (SCD). RH variation can drive loss of antigen epitopes or expression of new epitopes, predisposing patients with SCD to Rh alloimmunization. Serologic antigen typing is limited to common Rh antigens, necessitating a genetic approach to detect variant antigen expression. We developed a novel algorithm termed RHtyper for RH genotyping from existing whole-genome sequencing (WGS) data. RHtyper determined RH genotypes in an average of 3.4 and 3.3 minutes per sample for RHD and RHCE, respectively. In a validation cohort consisting of 57 patients with SCD, RHtyper achieved 100% accuracy for RHD and 98.2% accuracy for RHCE, when compared with genotypes obtained by RH BeadChip and targeted molecular assays and after verification by Sanger sequencing and independent next-generation sequencing assays. RHtyper was next applied to WGS data from an additional 827 patients with SCD. In the total cohort of 884 patients, RHtyper identified 38 RHD and 28 RHCE distinct alleles, including a novel RHD DAU allele, RHD* 602G, 733C, 744T 1136T. RHtyper provides comprehensive and high-throughput RH genotyping from WGS data, facilitating deconvolution of the extensive RH genetic variation among patients with SCD. We have implemented RHtyper as a cloud-based public access application in DNAnexus (https://platform.dnanexus.com/app/RHtyper), enabling clinicians and researchers to perform RH genotyping with next-generation sequencing data.

An open-source python library for detection of known and novel Kell, Duffy and Kidd variants from exome sequencing

BACKGROUND AND OBJECTIVES

Next generation sequencing (NGS) has promising applications in transfusion medicine. Exome sequencing (ES) is increasingly used in the clinical setting, and blood group interpretation is an additional value that could be extracted from existing data sets. We provide the first release of an open-source software tailored for this purpose and describe its validation with three blood group systems.

MATERIALS AND METHODS

The DTM-Tools algorithm was designed and used to analyse 1018 ES NGS files from the ClinSeq® cohort. Predictions were correlated with serology for 5 antigens in a subset of 108 blood samples. Discrepancies were investigated with alternative phenotyping and genotyping methods, including a long-read NGS platform.

RESULTS

Of 116 genomic variants queried, those corresponding to 18 known KEL, FY and JK alleles were identified in this cohort. 596 additional exonic variants were identified KEL, ACKR1 and SLC14A1, including 58 predicted frameshifts. Software predictions were validated by serology in 108 participants; one case in the FY blood group and three cases in the JK blood group were discrepant. Investigation revealed that these discrepancies resulted from (1) clerical error, (2) serologic failure to detect weak antigenic expression and (3) a frameshift variant absent in blood group databases.

CONCLUSION

DTM-Tools can be employed for rapid Kell, Duffy and Kidd blood group antigen prediction from existing ES data sets; for discrepancies detected in the validation data set, software predictions proved accurate. DTM-Tools is open-source and in continuous development.

Targeted exome sequencing designed for blood group, platelet, and neutrophil antigen investigations: Proof-of-principle study for a customized single-test system

BACKGROUND

Immunohematology reference laboratories provide red blood cell (RBC), platelet (PLT), and neutrophil typing to resolve complex cases, using serology and commercial DNA tests that define clinically important antigens. Broad-range exome sequencing panels that include blood group targets provide accurate blood group antigen predictions beyond those defined by serology and commercial typing systems and identify rare and novel variants. The aim of this study was to design and assess a panel for targeted exome sequencing of RBC, PLT, and neutrophil antigen-associated genes to provide a comprehensive profile in a single test, excluding unrelated gene targets.

STUDY DESIGN AND METHODS

An overlapping probe panel was designed for the coding regions of 64 genes and loci involved in gene expression. Sequencing was performed on 34 RBC and 17 PLT/neutrophil reference samples. Variant call outputs were analyzed using software to predict star allele diplotypes. Results were compared with serology and previous sequence genotyping data.

RESULTS

Average coverage exceeded 250×, with more than 94% of targets at Q30 quality or greater. Increased coverage revealed a variant in the Scianna system that was previously undetected. The software correctly predicted allele diplotypes for 99.5% of RBC blood groups tested and 100% of PLT and HNA antigens excepting HNA-2. Optimal throughput was 12 to 14 samples per run.

CONCLUSION

This single-test system demonstrates high coverage and quality, allowing for the detection of previously overlooked variants and increased sample throughput. This system has the potential to integrate genomic testing across laboratories within hematologic reference settings.

Overcoming the challenges of interpreting complex and uncommon RH alleles from whole genomes

BACKGROUND AND OBJECTIVES

Rh is one of the most diverse and complex blood group systems. Recently, next generation sequencing (NGS) has proven to be a viable option for RH genotyping. We have developed automated software (bloodTyper) for determining alleles encoding RBC antigens from NGS-based whole genome sequencing (WGS). The bloodTyper algorithm has not yet been optimized and evaluated for complex and uncommon RH alleles.

MATERIALS AND METHODS

Twenty-two samples with previous polymerase chain reaction (PCR) and Sanger sequencing-based RH genotyping underwent WGS. bloodTyper was used to detect RH alleles including those defined by structural variation (SV) using a combination of three independent strategies: sequence read depth of coverage, split reads and paired reads.

RESULTS

bloodTyper was programmed to identify D negative and positive phenotypes as well as the presence of alleles encoding weak D, partial D and variant RHCE. Sequence read depth of coverage calculation accurately determined RHD zygosity and detected the presence of RHD/RHCE hybrids. RHCE*C was determined by sequence read depth of coverage and by split read methods. RHD hybrid alleles and RHCE*C were confirmed by using a paired read approach. Small SVs present in RHCE*CeRN and RHCE*ceHAR were detected by a combined read depth of coverage and paired read approach.

CONCLUSIONS

The combination of several different interpretive approaches allowed for automated software based-RH genotyping of WGS data including RHD zygosity and complex compound RHD and RHCE heterozygotes. The scalable nature of this automated analysis will enable RH genotyping in large genomic sequencing projects.

Multiple GYPB gene deletions associated with the U- phenotype in those of African ancestry

BACKGROUND

The MNS blood group system is defined by three homologous genes: GYPA, GYPB, and GYPE. GYPB encodes for glycophorin B (GPB) carrying S/s and the "universal" antigen U. RBCs of approximately 1% of individuals of African ancestry are U- due to absence of GPB. The U- phenotype has long been attributed to a deletion encompassing GYPB exons 2 to 5 and GYPE exon 1 (GYPB*01N).

STUDY DESIGN AND METHODS

Samples from two U-individuals underwent Illumina short read whole genome sequencing (WGS) and Nanopore long read WGS. In addition, two existing WGS datasets, MedSeq (n = 110) and 1000 Genomes (1000G, n = 2535), were analyzed for GYPB deletions. Deletions were confirmed by Sanger sequencing. Twenty known U- donor samples were tested by a PCR assay to determine the specific deletion alleles present in African Americans.

RESULTS

Two large GYPB deletions in U- samples of African ancestry were identified: a 110 kb deletion extending left of GYPB (DEL_B_LEFT) and a 103 kb deletion extending right (DEL_B_RIGHT). DEL_B_LEFT and DEL_B_RIGHT were the most common GYPB deletions in the 1000 Genomes Project 669 African genomes (allele frequencies 0.04 and 0.02). Seven additional deletions involving GYPB were seen in African, Admixed American, and South Asian samples. No samples analyzed had GYPB*01N.

CONCLUSIONS

The U- phenotype in those of African ancestry is primarily associated with two different complete deletions of GYPB (with intact GYPE). Seven additional less common GYPB deletion backgrounds were found. GYPB*01N, long assumed to be the allele commonly encoding U- phenotypes, appears to be rare.

Defining Blood Group Gene Reference Alleles by Long-Read Sequencing: Proof of Concept in the ACKR1 Gene Encoding the Duffy Antigens

Background

In the novel era of blood group genomics, (re-)defining reference gene/allele sequences of blood group genes has become an important goal to achieve, both for diagnostic and research purposes. As novel potent sequencing technologies are available, we thought to investigate the variability encountered in the three most common alleles of ACKR1, the gene encoding the clinically relevant Duffy antigens, at the haplotype level by a long-read sequencing approach.

Materials and Methods

After long-range PCR amplification spanning the whole ACKR1 gene locus (∼2.5 kilobases), amplicons generated from 81 samples with known genotypes were sequenced in a single read by using the Pacific Biosciences (PacBio) single molecule, real-time (SMRT) sequencing technology.

Results

High-quality sequencing reads were obtained for the 162 alleles (accuracy >0.999). Twenty-two nucleotide variations reported in databases were identified, defining 19 haplotypes: four, eight, and seven haplotypes in 46 ACKR1*01, 63 ACKR1*02, and 53 ACKR1*02N.01 alleles, respectively.

Discussion

Overall, we have defined a subset of reference alleles by third-generation (long-read) sequencing. This technology, which provides a "longitudinal" overview of the loci of interest (several thousand base pairs) and is complementary to the second-generation (short-read) next-generation sequencing technology, is of critical interest for resolving novel, rare, and null alleles.

Frequency of Mia (MNS7) and Classification of Mia-Positive Hybrid Glycophorins in an Australian Blood Donor Population

Background

MNS blood group system genes GYPA and GYPB share a high degree of sequence homology and gene structure. Homologous exchanges between GYPA and GYPB form hybrid genes encoding hybrid glycophorins GP(A-B-A) and GP(B-A-B). Over 20 hybrid glycophorins have been characterised. Each has a distinct phenotype defined by the profile of antigens expressed including Mi^a. Seven hybrid glycophorins carry Mi^a and have been reported in Caucasian and Asian population groups. In Australia, the population is diverse; however, the prevalence of hybrid glycophorins in the population has never been determined. The aims of this study were to determine the frequency of Mi^a and to classify Mi^a-positive hybrid glycophorins in an Australian blood donor population.

Method

Blood samples from 5,098 Australian blood donors were randomly selected and screened for Mi^a using anti-Mi^a monoclonal antibody (CBC-172) by standard haemagglutination technique. Mi^a-positive red blood cells (RBCs) were further characterised using a panel of phenotyping reagents. Genotyping by high-resolution melting analysis and DNA sequencing were used to confirm serology.

Result

RBCs from 11/5,098 samples were Mi^a-positive, representing a frequency of 0.22%. Serological and molecular typing identified four types of Mi^a-positive hybrid glycophorins: GP.Hut (n = 2), GP.Vw (n = 3), GP.Mur (n = 5), and 1 GP.Bun (n = 1). GP.Mur was the most common.

Conclusion

This is the first comprehensive study on the frequency of Mi^a and types of hybrid glycophorins present in an Australian blood donor population. The demographics of Australia are diverse and ever-changing. Knowing the blood group profile in a population is essential to manage transfusion needs.

Banking with precision: transfusion medicine as a potential universal application in clinical genomics

PURPOSE OF REVIEW

To summarize the most recent scientific progress in transfusion medicine genomics and discuss its role within the broad genomic precision medicine model, with a focus on the unique computational and bioinformatic aspects of this emergent field.

RECENT FINDINGS

Recent publications continue to validate the feasibility of using next-generation sequencing (NGS) for blood group prediction with three distinct approaches: exome sequencing, whole genome sequencing, and PCR-based targeted NGS methods. The reported correlation of NGS with serologic and alternative genotyping methods ranges from 92 to 99%. NGS has demonstrated improved detection of weak antigens, structural changes, copy number variations, novel genomic variants, and microchimerism. Addition of a transfusion medicine interpretation to any clinically sequenced genome is proposed as a strategy to enhance the cost-effectiveness of precision genomic medicine. Interpretation of NGS in the blood group antigen context requires not only advanced immunohematology knowledge, but also specialized software and hardware resources, and a bioinformatics-trained workforce.

SUMMARY

Blood transfusions are a common inpatient procedure, making blood group genomics a promising facet of precision medicine research. Further efforts are needed to embrace transfusion bioinformatic challenges and evaluate its clinical utility.

The role of genomics in transfusion medicine

PURPOSE OF REVIEW

To summarize recent advances in red blood cell (RBC) blood group genotyping, with an emphasis on advances in the use of NGS next generation sequencing (NGS) to detect clinically relevant blood group gene variation.

RECENT FINDINGS

Genetic information is useful in predicting RBC blood group antigen expression in several clinical contexts, particularly, for patients at high-risk for allosensitization, such as multiple transfused patients. Blood group antigen expression is directed by DNA variants affecting multiply genes. With over 300 known antigens, NGS offers the attractive prospect of comprehensive blood group genotyping. Recent studies from several groups show that NGS reliably detects blood group gene single nucleotide variants (SNVs) with good correlation with other genetic methods and serology. Additionally, new custom NGS methods accurately detect complex DNA variants, including hybrid RH alleles. Thus, recent work shows that NGS detects known and novel blood group gene variants in patients, solves challenging clinical cases, and detects relevant blood group variation in donors.

SUMMARY

New work shows that NGS is particularly robust in identifying SNVs in blood group genes, whereas custom genomic tools can be used to identify known and novel complex structural variants, including in the RH system.

Automated typing of red blood cell and platelet antigens from whole exome sequences

BACKGROUND

Genotyping has expanded the number red blood cell (RBC) and platelet (PLT) antigens that can readily be typed, but often represents an additional testing cost. The analysis of existing genomic data offers a cost-effective approach. We recently developed automated software (bloodTyper) for determination of RBC and PLT antigens from whole genome sequencing. Here we extend the algorithm to whole exome sequencing (WES).

STUDY DESIGN AND METHODS

Whole exome sequencing was performed on samples from 75 individuals. WES-based bloodTyper RBC and PLT typing was compared to conventional polymerase chain reaction (PCR) RHD zygosity testing and serologic and single-nucleotide polymorphism (SNP) typing for 38 RBC antigens in 12 systems (17 serologic and 35 SNPs) and 22 PLT antigens (22 SNPs). Samples from the first 20 individuals were used to modify bloodTyper to interpret WES followed by blinded typing of 55 samples.

RESULTS

Over the first 20 samples, discordances were noted for C, M, and N antigens, which were due to WES-specific biases. After modification, bloodTyper was 100% accurate on blinded evaluation of the last 55 samples and outperformed both serologic (99.67% accurate) and SNP typing (99.97% accurate) reflected by two Fy^b and one N serologic typing errors and one undetected SNP encoding a Jk_null phenotype. RHD zygosity testing by bloodTyper was 100% concordant with a combination of hybrid Rhesus box PCR and PCR-restriction fragment length polymorphism for all samples.

CONCLUSION

The automated bloodTyper software was modified for WES biases to allow for accurate RBC and PLT antigen typing. Such analysis could become a routing part of future WES efforts.

Electron paramagnetic resonance oximetry as a novel approach to monitor the effectiveness and quality of red blood cell transfusions

BACKGROUND

The goal of red blood cell transfusion is to improve tissue oxygenation. Assessment of red blood cell quality and individualised therapeutic needs can be optimised using direct oxygen (O₂) measurements to guide treatment. Electron paramagnetic resonance oximetry is capable of accurate, repeatable and minimally invasive measurements of tissue pO₂. Here we present preclinical proof-of-concept of the utility of electron paramagnetic resonance oximetry in an experimental setting of acute blood loss, transfusion, and post-transfusion monitoring.

MATERIALS AND METHODS

Donor rat blood was collected, leucocyte-reduced, and stored at 4 °C in AS-3 for 1, 7 and 14 days. Red blood cell morphology, O₂ equilibrium, p50 and Hill numbers from O₂ binding and dissociation curves were evaluated in vitro. Recipient rats were bled and maintained at a mean arterial pressure of 30-40 mmHg and hind limb muscle (biceps femoris) pO₂ at 25-50% of baseline. Muscle pO₂ was monitored continuously over the course of experiments to assess the effectiveness of red blood cell preparations at different stages of blood loss and restoration.

RESULTS

Red blood cell morphology, O₂ equilibrium and p50 values of intra-erythrocyte haemoglobin were significantly altered by refrigerated storage for both 7 and 14 days. Transfusion of red blood cells stored for 7 or 14 days demonstrated an equivalently impaired ability to restore hind limb muscle pO₂, consistent with in vitro observations and transfusion with albumin. Red blood cells refrigerated for 1 day demonstrated normal morphology, in vitro oxygenation and in vivo restoration of tissue pO₂.

DISCUSSION

Electron paramagnetic resonance oximetry represents a useful approach to assessing the quality of red blood cells and subsequent transfusion effectiveness.

Blood group genotyping

Genomics is affecting all areas of medicine. In transfusion medicine, DNA-based genotyping is being used as an alternative to serological antibody-based methods to determine blood groups for matching donor to recipient. Most antigenic polymorphisms are due to single nucleotide polymorphism changes in the respective genes, and DNA arrays that target these changes have been validated by comparison with antibody-based typing. Importantly, the ability to test for antigens for which there are no serologic reagents is a major medical advance to identify antibodies and find compatible donor units, and can be life-saving. This review summarizes the evolving use and applications of genotyping for red cell and platelet blood group antigens affecting several areas of medicine. These include prenatal medicine for evaluating risk of fetal or neonatal disease and candidates for Rh-immune globulin; transplantation for bone marrow donor selection and transfusion support for highly alloimmunized patients and for confirmation of A2 status of kidney donors; hematology for comprehensive typing for patients with anemia requiring chronic transfusion; and oncology for patients receiving monoclonal antibody therapies that interfere with pretransfusion testing. A genomics approach allows, for the first time, the ability to routinely select donor units antigen matched to recipients for more than ABO/RhD to reduce complications. Of relevance, the growth of whole-genome sequencing in chronic disease and for general health will provide patients' comprehensive extended blood group profile as part of their medical record to be used to inform selection of the optimal transfusion therapy.

Molecular Detection of Glycophorins A and B Variant Phenotypes and their Clinical Relevance

Crossover or conversion between the homologous regions of glycophorin A (GYPA) and glycophorin B (GYPB) gives rise to several different hybrid glycophorin genes encoding a number of different glycophorin variant phenotypes which bear low prevalence antigens in the MNS blood group system. GP.Mur is the main glycophorin variant phenotype which causes hemolytic transfusion reaction (HTR) and hemolytic disease of the fetus and newborn (HDFN) in East and Southeast Asians. The detection of glycophorin variant phenotypes using serological methods is limited to phenotyping reagents that are not commercially available. Moreover, the red blood cells used for antibody identification are usually of the GP.Mur phenotype. The current Polymerase Chain Reaction (PCR)-based methods and loop-mediated isothermal amplification (LAMP) are available alternatives to phenotyping that allow for the specific detection of glycophorin variant phenotypes. This review highlights the molecular detection method for glycophorins A and B variant phenotypes and their clinical relevance.

Developments beyond blood group serology in the genomics era

Blood group serology and single nucleotide polymorphism-based genotyping platforms are accurate but do not provide a comprehensive cover for all 36 blood group systems and do not cover the antigen diversity observed among population groups. This review examines the extent to which genomics is shaping blood group serology. Resources for genomics include the Human Reference Genome Sequence assembly; curated blood group tables listing variants; public databases providing information on genetic variants from world-wide studies; and massively parallel sequencing technologies. Blood group genomic studies span the spectrum, from bioinformatic data mining of huge data sets containing whole genome and whole exome information to laboratory investigations utilising targeted sequencing approaches. Blood group predictions based on genome sequencing and genomic studies are proving accurate, and have shown utility in both research and reference settings. Overall, studies confirm the potential for blood group genomics to reshape donor and patient transfusion management strategies to provide more compatible blood transfusions.

Comprehensive blood group antigen profile predictions for Western Desert Indigenous Australians from whole exome sequence data

BACKGROUND

The distribution of RBC antigens, which define blood group types, differs among populations. In contrast to many world populations, blood group profiles for Indigenous Australians have not been well studied. As it is now possible to predict comprehensive blood group antigen profiles from genomic data sets, we aimed to apply this for Indigenous Australians and to provide a comparison to other major world populations.

STUDY DESIGN AND METHODS

Whole exome sequence data for 72 Western Desert Indigenous Australians was provided by the Telethon Kids Institute. Variants (against hg19) were annotated using computer software (ANNOVAR, Qiagen Bioinformatics) and filtered to include only variants in genes for 36 blood group systems, and the transcription factors KLF1 and GATA1. The RHCE*C allele and RHD zygosity were identified by copy number variant analysis of sequence alignments. The impact of missense variants was investigated in silico using a meta-predictor of disease-causing variants (Meta-SNP).

RESULTS

For 21 blood group systems the predicted blood group antigen frequencies were comparable to those for other major world populations. For 13 systems, interesting points of contrast were identified. Furthermore, we identified 12 novel variants, one novel D allele, and four rare variants with potential clinical significance.

CONCLUSION

This is the first systematic assessment of genomic data to elucidate blood group antigen profiles for Indigenous Australians who are linguistically and culturally diverse. Our study paves the way to understanding the geographic distribution of blood group variants in different Indigenous groups and the associated RBC phenotypes. This in turn is expected to guide transfusion practice for Indigenous individuals.

Genomic coordinates and continental distribution of 120 blood group variants reported by the 1000 Genomes Project

BACKGROUND

The 1000 Genomes Project provides a database of genomic variants from whole genome sequencing of 2504 individuals across five continental superpopulations. This database can enrich our background knowledge of worldwide blood group variant geographic distribution and identify novel variants of potential clinical significance.

STUDY DESIGN AND METHODS

The 1000 Genomes database was analyzed to 1) expand knowledge about continental distributions of known blood group variants, 2) identify novel variants with antigenic potential and their geographic association, and 3) establish a baseline scaffold of chromosomal coordinates to translate next-generation sequencing output files into a predicted red blood cell (RBC) phenotype.

RESULTS

Forty-two genes were investigated. A total of 604 known variants were mapped to the GRCh37 assembly; 120 of these were reported by 1000 Genomes in at least one superpopulation. All queried variants, including the ACKR1 promoter silencing mutation, are located within exon pull-down boundaries. The analysis yielded 41 novel population distributions for 34 known variants, as well as 12 novel blood group variants that warrant further validation and study. Four prediction algorithms collectively flagged 79 of 109 (72%) known antigenic or enzymatically detrimental blood group variants, while 4 of 12 variants that do not result in an altered RBC phenotype were flagged as deleterious.

CONCLUSION

Next-generation sequencing has known potential for high-throughput and extended RBC phenotype prediction; a database of GRCh37 and GRCh38 chromosomal coordinates for 120 worldwide blood group variants is provided as a basis for this clinical application.

Investigation of the variable In(Lu) phenotype caused by KLF1 variants

INTRODUCTION

KLF1 is an essential transcriptional activator that drives erythropoiesis. KLF1 variants can result in the Inhibitor of Lutheran, or In(Lu), phenotype where red blood cells (RBCs) have reduced BCAM (LU) and CD44 (IN). Other RBC surface molecules also have changed expression; however, there is controversy in the literature regarding which are truly impacted. We aimed to investigate KLF1 variants in the Australian population.

STUDY DESIGN AND METHODS

In(Lu) samples were sourced through screening and through the RBC reference laboratory. Blood donor samples (8036) were screened to identify weakened/absent Lu^b antigen. Samples were genotyped by massively parallel sequencing, while surface carbohydrates and blood group molecules were assessed by flow cytometry. Hemoglobin (Hb) types were analyzed by high-performance liquid chromatography.

RESULTS

Four of 8036 donors were identified to be In(Lu), and two previously identified In(Lu) samples were provided from the RBC reference laboratory. Five different KLF1 variants were identified; two were novel: c.954G>C/p.Trp318Cys and c.421C>T/p.Arg141*. BCAM and CD44 were reduced in all samples, consistent with previous reports. As a group, In(Lu) RBCs had reduced CD35 (KN), ICAM4 (LW), and CD147 (OK), and demonstrated increased binding of lectins ECA and SNAI. One In(Lu) sample had elevated HbF and another elevated HbA2.

CONCLUSION

Different KLF1 variants may potentially produce variable phenotypes. A framework for investigating KLF1 variants and their phenotypic impact has been provided. In the future, given available international databases, further testing algorithms (as advocated here) will allow for correlation of phenotype with genotype and therefore accurately document this variability between KLF1 variants.

Severe hemolytic disease of the fetus and newborn due to allo-anti-D in a patient with a partial DEL phenotype arising from the variant allele described as RHD*148+1T (RHD*01EL.31)

BACKGROUND

RhD DEL variants may show complete or partial expression of RhD epitopes. There have been only rare reports of anti-D causing hemolytic disease of the fetus and newborn (HDFN) in this context. We report a case of severe HDFN associated with a recently described DEL variant.

CASE REPORT

A multiparous woman presented with an allo-anti-D and showed incongruent phenotyping and genotyping results on initial study. Further investigations identified the RHD mutation, defined as RHD*148+1T and named RHD*01EL.31, which had been previously associated with a DEL phenotype. Extended RhD phenotyping by adsorption-elution showed that there was reactivity with four of nine monoclonal anti-D antibodies, suggesting a partial DEL phenotype. The first child showed no clinical evidence of HDFN, although the cord direct antiglobulin test was positive. The second child developed fetal anemia treated with intrauterine transfusion, and neonatal hyperbilirubinemia requiring exchange transfusion.

CONCLUSION

The RHD allele, RHD*148+1T, results in a partial Del phenotype, and the anti-D formed in pregnant women with this phenotype is capable of causing severe HDFN.

Massively parallel and multiplex blood group genotyping using next-generation-sequencing

OBJECTIVES

Thirty-six blood group systems are listed by the International Society of Blood Transfusion, containing almost 350 antigens. Most of these result from a single nucleotide polymorphism (SNP). Serology is the standard method for blood group typing. However, this technique has some limitations and cannot respond to the growing demand of blood product typing for a large number of antigens. Here we describe a blood group genotyping assay directly from whole blood samples using Next-Generation Sequencing (NGS), allowing the simultaneous identification of 15 SNPs associated with the blood group systems of 95 patients in a single run.

DESIGN AND METHOD

After an automated DNA extraction, targets are amplified by multiplex polymerase chain reaction (PCRm). Two panels addressing 9 groups have been developed (MNS, Lutheran, Kell, Duffy, Kidd, Diego, Yt, Dombrock, and Colton), one for 8 SNPs, the other for 7 SNPs. For each sample, both panels corresponding to 14 amplicons (1 amplicon containing 2 SNPs) are pooled. Then a dual-indexed library is generated from each pool by linking Illumina adaptors directly onto amplicons, followed by sequencing using the MiSeq platform (Illumina).

RESULTS

In a single experiment, 95 blood donor samples have been sequenced for the genes of interest. Among the 1425 targeted single nucleotide polymorphisms, 1420 were identified by sequencing, reflecting a coverage of 99.65%. The obtained data shows a good correlation (99% for all SNPs) with other blood group typing methods. Depending on the allele pairs analyzed, correlations vary between 97.12 and 100%.

CONCLUSION

Next-Generation sequencing would supplement serological and molecular techniques and, in the near future, could replace it with complete and fast results acquisition for pre-screening and identification of rare blood bags.

A DEL phenotype attributed to RHD Exon 9 sequence deletion: slipped-strand mispairing and blood group polymorphisms

BACKGROUND

The RhD blood group antigen is extremely polymorphic and the DEL phenotype represents one such class of polymorphisms. The DEL phenotype prevalent in East Asian populations arises from a synonymous substitution defined as RHD*1227A. However, initially, based on genomic and cDNA studies, the genetic basis for a DEL phenotype in Taiwan was attributed to a deletion of RHD Exon 9 that was never verified at the genomic level by any other independent group. Here we investigate the genetic basis for a Caucasian donor with a DEL partial D phenotype and compare the genomic findings to those initial molecular studies.

STUDY DESIGN AND METHODS

The 3'-region of the RHD gene was amplified by long-range polymerase chain reaction (PCR) for massively parallel sequencing. Primers were designed to encompass a deletion, flanking Exon 9, by standard PCR for Sanger sequencing. Targeted sequencing of exons and flanking introns was also performed.

RESULTS

Genomic DNA exhibited a 1012-bp deletion spanning from Intron 8, across Exon 9 into Intron 9. The deletion breakpoints occurred between two 25-bp repeat motifs flanking Exon 9 such that one repeat sequence remained.

CONCLUSION

Deletion mutations bordered by repeat sequences are a hallmark of slipped-strand mispairing (SSM) event. We propose this genetic mechanism generated the germline deletion in the Caucasian donor. Extensive studies show that the RHD*1227A is the most prevalent DEL allele in East Asian populations and may have confounded the initial molecular studies. Review of the literature revealed that the SSM model explains some of the extreme polymorphisms observed in the clinically significant RhD blood group antigen.

Non-invasive fetal RHD genotyping for RhD negative women stratified into RHD gene deletion or variant groups: comparative accuracy using two blood collection tube types

Non-invasive fetal RHD genotyping in Australia to reduce anti-D usage will need to accommodate both prolonged sample transport times and a diverse population demographic harbouring a range of RHD blood group gene variants. We compared RHD genotyping accuracy using two blood sample collection tube types for RhD negative women stratified into deleted RHD gene haplotype and RHD gene variant cohorts. Maternal blood samples were collected into EDTA and cell-free (cf)DNA stabilising (BCT) tubes from two sites, one interstate. Automated DNA extraction and polymerase chain reaction (PCR) were used to amplify RHD exons 5 and 10 and CCR5. Automated analysis flagged maternal RHD variants, which were classified by genotyping. Time between sample collection and processing ranged from 2.9 to 187.5 hours. cfDNA levels increased with time for EDTA (range 0.03-138 ng/μL) but not BCT samples (0.01-3.24 ng/μL). For the 'deleted' cohort (n=647) all fetal RHD genotyping outcomes were concordant, excepting for one unexplained false negative EDTA sample. Matched against cord RhD serology, negative predictive values using BCT and EDTA tubes were 100% and 99.6%, respectively. Positive predictive values were 99.7% for both types. Overall 37.2% of subjects carried an RhD negative baby. The 'variant' cohort (n=15) included one novel RHD and eight hybrid or African pseudogene variants. Review for fetal RHD specific signals, based on one exon, showed three EDTA samples discordant to BCT, attributed to high maternal cfDNA levels arising from prolonged transport times. For the deleted haplotype cohort, fetal RHD genotyping accuracy was comparable for samples collected in EDTA and BCT tubes despite higher cfDNA levels in the EDTA tubes. Capacity to predict fetal RHD genotype for maternal carriers of hybrid or pseudogene RHD variants requires stringent control of cfDNA levels. We conclude that fetal RHD genotyping is feasible in the Australian environment to avoid unnecessary anti-D immunoglobulin prophylaxis.

Targeted exome sequencing defines novel and rare variants in complex blood group serology cases for a red blood cell reference laboratory setting

BACKGROUND

We previously demonstrated that targeted exome sequencing accurately defined blood group genotypes for reference panel samples characterized by serology and single-nucleotide polymorphism (SNP) genotyping. Here we investigate the application of this approach to resolve problematic serology and SNP-typing cases.

STUDY DESIGN AND METHODS

The TruSight One sequencing panel and MiSeq platform was used for sequencing. CLC Genomics Workbench software was used for data analysis of the blood group genes implicated in the serology and SNP-typing problem. Sequence variants were compared to public databases listing blood group alleles. The effect of predicted amino acid changes on protein function for novel alleles was assessed using SIFT and PolyPhen-2.

RESULTS

Among 29 unresolved samples, sequencing defined SNPs in blood group genes consistent with serologic observation: 22 samples exhibited SNPs associated with varied but known blood group alleles and one sample exhibited a chimeric RH genotype. Three samples showed novel variants in the CROM, LAN, and RH systems, respectively, predicting respective amino acid changes with possible deleterious impact. Two samples harbored rare variants in the RH and FY systems, respectively, not previously associated with a blood group allele or phenotype. A final sample comprised a rare variant within the KLF1 transcription factor gene that may modulate DNA-binding activity.

CONCLUSION

Targeted exome sequencing resolved complex serology problems and defined both novel blood group alleles (CD55:c.203G>A, ABCB6:c.1118_1124delCGGATCG, ABCB6:c.1656-1G>A, and RHD:c.452G>A) and rare variants on blood group alleles associated with altered phenotypes. This study illustrates the utility of exome sequencing, in conjunction with serology, as an alternative approach to resolve complex cases.