A preliminary evaluation of next-generation sequencing as a screening tool for targeted genotyping of erythrocyte and platelet antigens in blood donors

Simultaneous genotyping for human platelet antigen systems and HLA-A and HLA-B loci by targeted next-generation sequencing

In order to treat the alloimmunization platelet transfusion refractoriness (PTR), human leukocyte antigen (HLA)-type and/or human platelet antigen (HPA)-type matched platelets between donors and patients are usually used. Therefore, genotyping of HLA-A and HLA-B loci, as well as HPA systems, for donors and patients, is of great significance. However, there is a rare report of genotyping for HLA-A and HLA-B loci as well as HPA systems at the same time. In this study, a high-throughput method for simultaneous genotyping of HLA-A and HLA-B loci, as well as HPA genotyping, was developed. A RNA capture probe panel was designed covering all exon sequences of the GP1BA, GP1BB, ITGA2, CD109, ITGB3, and ITGA2B genes and HLA-A and HLA-B loci. The HLA-A, HLA-B, and 34 HPA systems were genotyped using a targeted next-generation sequencing (NGS) method. The genotypes of the HLA-A and HLA-B loci, as well as the HPA, were assigned based on the nucleotides in the polymorphism sites. Using the NGS method, 204 unrelated blood specimens were successfully genotyped for all 34 HPA systems as well as HLA-A and HLA-B loci. The accuracy of the NGS method was 100%. Only HPA-2, HPA-3, HPA-5, HPA-6w, HPA-15, and HPA-21w showed polymorphism with frequencies of 0.9412, 0.6863, 0.9853, 0.9779, 0.4314, and 0.9951 for a allele, respectively. Thirty-two single nucleotide variants (SNVs) were detected. Of them, 12 SNVs can lead to amino acid change. HLA-A*11:01 and HLA-B*46:01 are the most common alleles for HLA-A and HLA-B loci. A targeted next-generation sequencing method for simultaneously genotyping HPA systems and HLA-A and HLA-B loci was first established, which could be used to create a database of HLA-typed and/or HPA-typed unrelated donors.

High-Throughput Next-Generation Sequencing of the Kidd Blood Group: Unexpected Antigen Expression Properties of Four Alleles and Detection of Novel Variants

Background

The blood supply for patients with foreign ethnic backgrounds can be challenging, as they often have blood group and HPA patterns that differ from the variants prevalent in the German population. In addition, hemoglobinopathies requiring regular blood transfusion may be more common in such populations. High-throughput genotyping tests can facilitate the identification of the most compatible blood products, thereby reducing the risk of transfusion reactions. The present study reports the results of a molecular study for the Kidd (JK) blood group. Allele frequencies and antigen prevalence data are presented for >8,000 individuals of various origins.

Material and Methods

More than 8,000 blood donors were genotyped for 22 blood group systems and 5 HPA genes using an amplicon-based next-generation sequencing (NGS) approach. As part of the test system, we focused on the JK system in more detail. Double-ARMS PCR analysis was performed for the haplotype phasing of the JK1/JK2 and two more common synonymous polymorphisms. We performed transcript analysis to detect potential alternative splice products. For a subset of samples, a comparison between serotype and red cell genotype was conducted. Allele frequencies were determined for geographically different panels of individuals.

Results

We successfully genotyped the JK blood group for 99.6% of the samples. Haplotype phasing revealed 96 different alleles. For several alleles that carry one of the synonymous SNVs c.588A>G and c.810G>A, we could not confirm the reported JK phenotypes. We found a higher frequency of JK:1 alleles for all populations except Iraqis. JK*01W.01 alleles were more common in the Asian groups and sub-Saharan Africans. A variant of the allele JK*02N.01 was present exclusively in Southeast Asians.

Conclusion

Genotyping for JK antigens with a targeted NGS assay can easily be performed in routine. The interpretation that c.588A>G leads to a weak phenotype and c.810G>A to a null phenotype is questionable. IDs as well as the descriptions of alleles carrying these SNVs should be revised in the ISBT JK table.

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.

Diversity of maize (Zea mays L.) rhizobacteria with potential to promote plant growth

Plant growth-limiting factors, such as low nutrient availability and weak pathogen resistance, may hinder the production of several crops. Plant growth-promoting bacteria (PGPB) used in agriculture, which stimulate plant growth and development, can serve as a potential tool to mitigate or even circumvent these limitations. The present study evaluated the feasibility of using bacteria isolated from the maize rhizosphere as PGPB for the cultivation of this crop. A total of 282 isolates were collected and clustered into 57 groups based on their genetic similarity using BOX-PCR. A representative isolate from each group was selected and identified at the genus level with 16S rRNA sequencing. The identified genera included Bacillus (61.5% of the isolates), Lysinibacillus (30.52%), Pseudomonas (3.15%), Stenotrophomonas (2.91%), Paenibacillus (1.22%), Enterobacter (0.25%), Rhizobium (0.25%), and Atlantibacter (0.25%). Eleven isolates with the highest performance were selected for analyzing the possible pathways underlying plant growth promotion using biochemical and molecular techniques. Of the selected isolates, 90.9% were positive for indolepyruvate/phenylpyruvate decarboxylase, 54.4% for pyrroloquinoline quinine synthase, 36.4% for nitrogenase reductase, and 27.3% for nitrite reductase. Based on biochemical characterization, 9.1% isolates could fix nitrogen, 36.6% could solubilize phosphate, 54.5% could produce siderophores, and 90.9% could produce indole acetic acid. Enzymatic profiling revealed that the isolates could degrade starch (90.1%), cellulose (72.7%), pectin (81.8%), protein (90.9%), chitin (18.2%), urea (54.5%), and esters (45.4%). Based on the data obtained, we identified three Bacillus spp. (LGMB12, LGMB273, and LGMB426), one Stenotrophomonas sp. (LGMB417), and one Pseudomonas sp. (LGMB456) with the potential to serve as PGPB for maize. Further research is warranted to evaluate the biotechnological potential of these isolates as biofertilizers under field conditions.

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.

High-Throughput Screening of Blood Donors for Twelve Human Platelet Antigen Systems Using Next-Generation Sequencing Reveals Detection of Rare Polymorphisms and Two Novel Protein-Changing Variants

Background

Exposure to non-matching human platelet alloantigens (HPA) may result in alloimmunization. Antibodies to HPA can be responsible for post-transfusion purpura, refractoriness to donor platelets, and fetal and neonatal alloimmune thrombocytopenia. For the supply of compatible apheresis platelet concentrates, the HPA genotypes are determined in a routine manner.

Methods

Here, we describe a novel method for genotyping twelve different HPA systems simultaneously, including HPA-1 to HPA-5, HPA-9w, HPA-10w, HPA-16w, HPA-19w, HPA-27w, and the novel HPA-34w by means of amplicon-based next-generation sequencing (NGS). Blood donor samples of 757 individuals with a migration background and 547 of Western European ancestry were genotyped in a mass-screening setup. An in-house software was developed for fast and automatic analysis. TaqMan assay and Sanger sequencing results served for validation of the NGS workflow. Finally, blood donors were divided in several groups based on their country of origin and the allele frequencies were compared.

Results

For 1,299 of 1,304 samples (99.6%) NGS was successfully performed. The concordance with TaqMan assay and Sanger sequencing results was 99.8%. Allele-calling dropouts that were observed for two samples with the TaqMan assay caused by rare single nucleotide polymorphisms were resolved by NGS. Additionally, twenty rare and two novel variants in the coding regions of the genes ITGB3, GPB1A, ITGBA2, and CD109 were detected. The determined allele frequencies were similar to those published in the gnomAD database.

Conclusions

No significant differences were observed in the distribution of allele frequencies of HPA-1 through HPA-5 and HPA-15 throughout the analyzed groups except for a lower allele frequency for the HPA-1b allele in the group of donors with Southern Asian ancestry. In contrast, other nucleotide variants that have not yet been phenotypically characterized occurred three times more often in blood donors with a migration background. High-throughput amplicon-based NGS is a reliable method for screening HPA genotypes in a large sample cohort simultaneously. It is easily upgradeable for genotyping additional targets without changing the setup or the analysis pipeline. Mass-screening methods will help building up blood donor registries to provide matched blood products.

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.

Next-Generation Sequencing Technologies in Blood Group Typing

Sequencing of the human genome has led to the definition of the genes for most of the relevant blood group systems, and the polymorphisms responsible for most of the clinically relevant blood group antigens are characterized. Molecular blood group typing is used in situations where erythrocytes are not available or where serological testing was inconclusive or not possible due to the lack of antisera. Also, molecular testing may be more cost-effective in certain situations. Molecular typing approaches are mostly based on either PCR with specific primers, DNA hybridization, or DNA sequencing. Particularly the transition of sequencing techniques from Sanger-based sequencing to next-generation sequencing (NGS) technologies has led to exciting new possibilities in blood group genotyping. We describe briefly the currently available NGS platforms and their specifications, depict the genetic background of blood group polymorphisms, and discuss applications for NGS approaches in immunohematology. As an example, we delineate a protocol for large-scale donor blood group screening established and in use at our institution. Furthermore, we discuss technical challenges and limitations as well as the prospect for future developments, including long-read sequencing technologies.

Complete RHD next-generation sequencing: establishment of reference RHD alleles

The Rh blood group system (ISBT004) is the second most important blood group after ABO and is the most polymorphic one, with 55 antigens encoded by 2 genes, RHD and RHCE This research uses next-generation sequencing (NGS) to sequence the complete RHD gene by amplifying the whole gene using overlapping long-range polymerase chain reaction (LR-PCR) amplicons. The aim was to study different RHD alleles present in the population to establish reference RHD allele sequences by using the analysis of intronic single-nucleotide polymorphisms (SNPs) and their correlation to a specific Rh haplotype. Genomic DNA samples (n = 69) from blood donors of different serologically predicted genotypes including R₁R₁ (DCe/DCe), R₂R₂ (DcE/DcE), R₁R₂ (DCe/DcE), R₂R_Z (DcE/DCE), R₁r (DCe/dce), R₂r (DcE/dce), and R₀r (Dce/dce) were sequenced and data were then mapped to the human genome reference sequence hg38. We focused on the analysis of hemizygous samples, as these by definition will only have a single copy of RHD For the 69 samples sequenced, different exonic SNPs were detected that correlate with known variants. Multiple intronic SNPs were found in all samples: 21 intronic SNPs were present in all samples indicating their specificity to the RHD*DAU0 (RHD*10.00) haplotype which the hg38 reference sequence encodes. Twenty-three intronic SNPs were found to be R₂ haplotype specific, and 15 were linked to R₁, R₀, and R_Z haplotypes. In conclusion, intronic SNPs may represent a novel diagnostic approach to investigate known and novel variants of the RHD and RHCE genes, while being a useful approach to establish reference RHD allele sequences.

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.

Automated typing of red blood cell and platelet antigens: a whole-genome sequencing study

BACKGROUND

There are more than 300 known red blood cell (RBC) antigens and 33 platelet antigens that differ between individuals. Sensitisation to antigens is a serious complication that can occur in prenatal medicine and after blood transfusion, particularly for patients who require multiple transfusions. Although pre-transfusion compatibility testing largely relies on serological methods, reagents are not available for many antigens. Methods based on single-nucleotide polymorphism (SNP) arrays have been used, but typing for ABO and Rh-the most important blood groups-cannot be done with SNP typing alone. We aimed to develop a novel method based on whole-genome sequencing to identify RBC and platelet antigens.

METHODS

This whole-genome sequencing study is a subanalysis of data from patients in the whole-genome sequencing arm of the MedSeq Project randomised controlled trial (NCT01736566) with no measured patient outcomes. We created a database of molecular changes in RBC and platelet antigens and developed an automated antigen-typing algorithm based on whole-genome sequencing (bloodTyper). This algorithm was iteratively improved to address cis-trans haplotype ambiguities and homologous gene alignments. Whole-genome sequencing data from 110 MedSeq participants (30 × depth) were used to initially validate bloodTyper through comparison with conventional serology and SNP methods for typing of 38 RBC antigens in 12 blood-group systems and 22 human platelet antigens. bloodTyper was further validated with whole-genome sequencing data from 200 INTERVAL trial participants (15 × depth) with serological comparisons.

FINDINGS

We iteratively improved bloodTyper by comparing its typing results with conventional serological and SNP typing in three rounds of testing. The initial whole-genome sequencing typing algorithm was 99·5% concordant across the first 20 MedSeq genomes. Addressing discordances led to development of an improved algorithm that was 99·8% concordant for the remaining 90 MedSeq genomes. Additional modifications led to the final algorithm, which was 99·2% concordant across 200 INTERVAL genomes (or 99·9% after adjustment for the lower depth of coverage).

INTERPRETATION

By enabling more precise antigen-matching of patients with blood donors, antigen typing based on whole-genome sequencing provides a novel approach to improve transfusion outcomes with the potential to transform the practice of transfusion medicine.

FUNDING

National Human Genome Research Institute, Doris Duke Charitable Foundation, National Health Service Blood and Transplant, National Institute for Health Research, and Wellcome Trust.