Skip to main content

Molecular and computational analysis of 45 samples with a serologic weak D phenotype detected among 132,479 blood donors in northeast China



RH1 is one of the most clinically important blood group antigens in the field of transfusion and in the prevention of fetal incompatibility. The molecular analysis and characterization of serologic weak D phenotypes is essential to ensuring transfusion safety.


Blood samples from a northeastern Chinese population were randomly screened for a serologic weak D phenotype. The nucleotide sequences of all 10 exons, adjacent flanking intronic regions, and partial 5′ and 3′ untranslated regions (UTRs) were detected for RHD genes. Predicted deleterious structural changes in missense mutations of serologicl weak D phenotypes were analyzed using SIFT, PROVEAN and PolyPhen2 software. The protein structure of serologic weak D phenotypes was predicted using Swiss-PdbViewer 4.0.1.


A serologic weak D phenotype was found in 45 individuals (0.03%) among 132,479 blood donors. Seventeen distinct RHD mutation alleles were detected, with 11 weak D, four partial D and two DEL alleles. Further analyses resulted in the identification of two novel alleles (RHD weak D 1102A and 399C). The prediction of a three-dimensional structure showed that the protein conformation was disrupted in 16 serologic weak D phenotypes.


Two novel and 15 rare RHD alleles were identified. Weak D type 15, DVI Type 3, and RHD1227A were the most prevalent D variant alleles in a northeastern Chinese population. Although the frequencies of the D variant alleles presented herein were low, their phenotypic and genotypic descriptions add to the repertoire of reported RHD alleles. Bioinformatics analysis on RhD protein can give us more interpretation of missense variants of RHD gene.


The number of blood group antigens currently recognized by the International Society of Blood Transfusion is 360, and 322 of them are clustered within 36 blood group systems [1]. The Rh blood group system is the most complex among all blood group systems [2]. The D (RH1) antigen is the most immunogenic and clinically significant antigen, which directly affects the hemolytic transfusion reaction and hemolytic diseases of fetuses and newborns [2, 3]. Besides D-positive and D-negative, RhD blood groups have multiple variants, including weak D, partial D and DEL phenotype [4]. A Working Group of the American Association of Blood Banks and College of American Pathologists published a proposal to use the term “serologic weak D phenotype” to distinguish the results of serological weak D testing using anti-human globulin with those of weak D genotyping based on molecular methods [5, 6].

The genetic alterations of RHD alleles differentially influence the RhD protein expression level and the number of RhD epitopes [7]. Weak D and most DEL have all D epitopes; partial D lack one or more D epitopes [3]. Anti-D production in RhD negative recipients transfused with blood from weak D, particle D and DEL donors has been reported [8,9,10,11]. Therefore, it is of great significance to accurately determine a serologic weak D phenotype. Nearly all serologic weak D phenotypes can be traced back to changes at the DNA level, including nonsense mutations, missense and synonymous mutations, frame shifts, unequal exchange, gene exchange and gene deletion, as well as others. This enables us to study the molecular mechanism of a serologic weak D phenotype at the gene level, and to determine the type of serologic weak D phenotype.

To date, more than 460 RHD alleles have been registered and nominated [12,13,14]. Serologic weak D phenotypes are often found in the blood samples of blood donors and patients, and molecular studies have been mainly conducted in Caucasian and African populations [15,16,17,18,19,20,21,22,23]. Corresponding research on the diversity of serologic weak D phenotypes has been reported in southern [24,25,26,27,28], but rarely in northeastern populations in China. Herein, we tested samples from a cohort of 132,479 blood donors in northeastern China for the serologic weak D phenotype. We subsequently sequenced the RHD gene of the 45 samples identified in this study as showing a serologic weak D phenotype. We built and optimized three-dimensional (3D) models of serologic weak D phenotypes identified in this study to explore their effect on RhD protein structure. Bioinformatics tools were employed to provide computational predictions on the RhD protein structure of serologic weak D phenotypes and enhance our understanding of how mutations affect phenotypes at the same time.

Materials and methods

Study participants

All 132,479 samples were collected from blood donors at the Blood Center of Liaoning Province, which is located in northeastern China, over a 5-year period (January 2012 to December 2016). Some donors may have donated repeatedly, which is common in similar large studies of the past and is known not to affect the statistics and conclusions. The study was approved by the Ethics Committee of the Liaoning Blood Center, Liaoning, China.

Serological studies

The D antigen was serologically determined using a monoclonal anti-D reagent (IgM, Clone BS226, Bio-Rad Medical Diagnostics GmbH Industriestrabe, Germany) using a microplate test protocol and a fully automated blood grouping instrument (Hemo-Type automatic blood group analyzer; GSG Robotix, Milan, Italy). For the microplate test protocol and testing in a fully automatic blood grouping instrument, 6 μL of whole blood from a sample tube was added to 333 μL of 0.9% saline in a test tube and mixed. An erythrocyte suspension (35 μL) was absorbed to a micropore and mixed with 25 μL of anti-D reagent in accordance with the manufacturer’s instructions.

We retested all samples showing negative or equivocal agglutination by tube method with three commercially available anti-D reagents (1-IgM/IgG blend, Clone D175-2 and D415 1E4, Dominion Biologicals Ltd., Dartmouth, Nova Scotia, Canada; 2-IgM/IgG blend, Clone P3X61, P3X21223B10, P3X290 and P3X35, Diagast Ltd, Loos, France; 3-IgM/IgG blend, Clone TH-28 and MS-26, Merck Millipore Ltd, Livingston, UK). An indirect antiglobulin test (IAT) was performed in the case of negative reactions. RhD-negative phenotyping was performed for all samples by hemagglutination using two techniques (microtiter plate and tube); samples were not identified by serologic adsorption–elution techniques. We performed an antibody screening test for blood group alloantibodies on all D-negative and serologic weak D phenotype samples using a gel method (DiaMed-ID, microtyping system, DiaMed China Limited, Hong Kong, China). Routine Rh typing for C, c, E, and e antigens was performed by the tube method with commercial monoclonal immunoglobulin (IgM) reagents (anti-C, Clone MS-24; anti-c, Clone MS-33; anti-E, MS-80 + 258; and anti-e, Clone MS-16 + 21 + 63, Shanghai Hemo-Pharmaceutical & Biological, Inc., Shanghai, China).

Molecular analysis of genomic DNA

Genomic DNA was extracted from a 0.2 mL blood sample using a DNA whole blood isolation kit (Tiangen Biotech, Beijing, China) in accordance with the manufacturer’s instructions. The RHD gene was sequenced in all serologic weak D phenotypes and 117 D antigen negative samples by IAT as previously described [24]. The nucleotide sequences of all 10 exons as well as adjacent flanking intronic regions, including partial 5′ and 3′ untranslated regions (UTRs), were determined (Table 1). Genomic DNA (50 to 100 ng) was used in a 25 μL reaction mix containing 200 mM dNTPs, 0.1 mM of each specific primer, 1.5 mM MgCl2, 1× PCR buffer, and 1 unit of GoTaq polymerase (Promega, Madison, WI, USA), supplemented with ddH2O. The following PCR program was used: 5 min of denaturation at 95 °C, 35 cycles of 30 s at 94 °C, 30 s at 62 °C (exons 1, 3, 4, 6–10), 30 s at 58 °C (exons 2, 5), and 1 min at 72 °C, followed by a final 10-min extension at 72 °C. The PCR procedure was carried out in a PE-9700 thermal cycler (Applied Biosystems, Foster City, CA, USA). Sequencing data were analyzed with FinchTV software (Geospiza Inc., Seattle, WA, USA) and all results compared to a NCBI Reference Sequence (RefSeq) database number NG_007494.1. The amino acid alignment of RhD was analyzed by CLUSTAL X (version 2.1) and the amino acid sequences used in the analysis were obtained from a protein database ( RHD zygosity was determined on all sequencing samples by the presence or absence of a hybrid Rhesus box as described [29].

Table 1 Primers for RHD gene amplification and sequencing

Statistical analysis

Allele frequencies were calculated from corresponding genotype counts. According to Hardy–Weinberg equilibrium, genotype frequency of D negative homozygote is equal to the square of the D negative allele frequency, and genotype frequency of the heterozygote (D variant and D negative) is equal to twice the product of the two allele frequencies. Allele frequencies for each molecular background of serologic weak D and D negative phenotypes were calculated.

Computational modeling of RhD protein and amino acid substitutions

The 3D structure of the RhD protein was visualized using Swiss-Pdb Viewer 4.1.0 (, which was used to generate models of the selected protein structure for the corresponding amino acid substitutions [30, 31]. Sorting Intolerant From Tolerant (SIFT) [32], Polymorphism Phenotyping algorithmV2 (PolyPhen-2) [33] and Protein Variation Effect Analyzer (PROVEAN) [34] software were used to predict the impact of amino acid substitutions on RhD protein structure.


Serological studies

Using routine methods, we screened 132,479 blood donors for D antigen, 131,939 of whom were found to be D + (99.592%), 495 were D − (0.374%), and 45 (0.034%) had a serologic weak D phenotype [5]. They were sorted based on the anti-D agglutination strength using the two routine techniques, and were also tested with three monoclonal anti-D reagents. No blood group alloantibody was detected in the plasmas of 45 serologic weak D phenotype samples.

Molecular characterization of D variants

We determined the RHD sequence of the coding region, adjacent flanking intronic regions, and the partial 5′ and 3′ UTRs in all 45 serologic weak D phenotype samples and detected 17 distinct alleles (Table 2). Twenty-nine individuals were weak D types representing 11 distinct alleles. Twelve individuals were partial D phenotypes representing four distinct alleles. Four individuals were DEL phenotypes representing two distinct alleles. Two novel single-nucleotide missense mutations (c.399G > C and c.1102G > A) were identified in three individuals. The DNA sequences of two novel alleles have been submitted to the GenBank Database (accession numbers KM282169 and KM282170, respectively). Six individuals carried complex alleles: one individual carried weak D 101G (c.101A > G) and weak D type15 (c.845 G > A) alleles, while another five individuals had a weak D (c.399G > C, c.845 G > A, c.763G > C, c.779A > G, c.1212C > A, respectively) and RHD1227A allele (c.1227G > A). One each of weak D type 72 (c.1212C > A), weak D type 18 (c.19C > T), weak D type 25 (c.341G > A), weak D type 31 (c.17C > T) and weak D type 54 (c.365C > T) were identified.

Table 2 RHD alleles found among donors with a serologic weak D phenotype

Calculation of allele frequencies

In our study, phenotype frequencies of serologic weak D and D negative phenotypes were 0.00034 and 0.00374, respectively. According to calculation, the D negative allele frequency was 0.06113, and frequency of a single D variant allele observation was 0.00005. In 117 D negative samples, 66.67% of the genetically tested RhD negatives seemed to be homozygous for the RHD deletion, and 33.33% had a deletion of RHD plus a variant RHD or a RHD-CE-D hybrid gene. Cumulative RHD 1227A allele frequency among weak D and D negative phenotypes was 0.00853. Expected homozygous deletion of RHD had a genotype frequency of 0.00244. Allele frequencies of serologic weak D phenotypes and D negative alleles were calculated as shown in Tables 3 and 4. Raw data for calculation of serologic weak D and RhD negative allele frequencies are shown in Additional file 1: Table S1.

Table 3 Allele frequencies of serologic weak D phenotype alleles
Table 4 Allele frequencies of RhD negative alleles

Predicted effect of nonsynonymous substitutions

SIFT, PolyPhen-2, and PROVEAN bioinformatic software programs were used to predict deleterious structural changes induced by nonsynonymous p.P6L, p.R7W, p.R10W, p.Y34C, p.R114Q, p.S122L, p.K133N, p.G255R, H260R, p.G282D, p.G368R, and p.D404E substitutions (Additional file 2: Table S2). The 12 nonsynonymous substitutions were distributed along the whole length of the RHD coding sequence without any apparent clustering (Fig. 1 ) [4, 45].

Fig. 1

Positions of single amino acid substitution of the RhD protein (adapted from Flegel [4] and Srivastava [45]). There are 417 amino acids in the RhD protein, shown here as circles. The mature protein in the cell membrane lacks the first amino acid. The nine exon boundaries in the RHD cDNA as reflected in the amino acid sequence are labeled as black bars. All detected amino acid substitutions encoding D variant alleles are labeled as colored circles. A synonymous single nucleotide polymorphism (SNP) caused no amino acid change (gray). The other SNPs are nonsynonymous and cause amino acid changes that are predicted to affect the RhD protein structure (red) or to be neutral (blue)

Bioinformatics analysis of RhD protein structure model

The template of RhD protein homology model was in accordance with the model based on computational hydropathy map [45, 46]. The model comprised 409 amino acids from Ser3 to Pro411 and lacked nine residues (two in the N terminus and six in the C terminus, Additional file 3: Table S3). The RhD 3D protein structures of the wild-type and 16 serologic weak D phenotypes highlighted the change in structure with altered amino acids. A 3D structure analysis of 16 serologic weak D phenotypes predicted amino acid position shifts in intracellular and exofacial loops, and the transmembraneous domain. The 3D structure model also demonstrated that p.P6L, p.R7W, p.R10W, p.Y34C, p.R114Q, p.S122L, p.K133 N, p.G280D, p.G368R and p.D404E mutations led to the disappearance of beta sheets, and position changes due to p.G255R and p.H260R mutations in beta sheets. The 3D structure model of four partial D types displayed the disappearance of beta sheets in DVI type 3 and DV type 2, an increase in beta sheets at amino acids 38–40 and 42–44 in DVI type 4, and no change in beta sheets in DFR type 2 only, with an amino acid position shift in intracellular and exofacial loops, and the transmembraneous domain. The 3D structure model of four partial D types displayed the disappearance of beta sheets in DVI type 3 and DV type 2, an increase in beta sheets at amino acids 38–40 and 42–44 in DVI type 4, and no change in beta sheets in DFR type 2 only, with an amino acid position shift in intracellular and exofacial loops, and the transmembraneous domain (Fig. 2).

Fig. 2

Comparison of the tertiary structure of modeled wild-type and mutant RhD proteins. RHD and D variant tertiary structures were modeled using Deep View-Swiss-Pdb Viewer 4.0.1. Helices are shown as deep pink and rose pink ribbons, sheets as yellow and orange, and coils as green. a1, b1, c1, and d1 are wild-type RhD protein. The weak D proteins of weak D type 31 (a2), weak D type 18 (a3), weak D type 61 (a4), weak D 101G (b2), weak D type 25 (b3), weak D 399C (b4), weak D type 54 (c2), weak D type 763C (c3) , weak D type15 (c4), weak D779G (d2), weak D 1102A (d3), and weak D type 72 (d4) compared with the left wild-type RhD proteins, respectively. The partial D proteins of DVI type 3 (e1), DVI type 4 (e2), D V type 2 (e3) and DFR type2 (e4) compared with the wild-type RhD protein (a1), respectively. The mutant weak D proteins (p.P6L, p.R7W, p.R10W, p.Y34C, p. R114Q, p.K133N, p.S122L, p.G255R, p.G280D, p.H260R, p.G368R and p.D404E) are highlighted


In the present study, we investigated the molecular characteristics of serologic weak D phenotypes in a northeastern population in China. The bioinformatics of 17 variant RHD alleles for serologic weak D phenotypes, including two novel alleles, were analyzed. The RHD allele distribution varies widely between Asian [24,25,26,27,28, 47,48,49] and European centers [15,16,17,18,19,20,21]. Differences in populations and routine serologic screening procedures employed, as well as in the molecular examinations used, may account for such differences to date, highlighting the need for standardization. In this study, Weak D type 15, DVI Type 3 and DEL (RHD1227A) were the most prevalent D variant alleles measured in the northeastern Chinese population, which were consistent with those reported in southern Chinese population [24, 25, 27]; however, they were rare in other populations. Therefore, the frequency of distribution of serologic weak D phenotypes varies among populations and ethnic groups.

Our tests detected two mutation types for DEL variants. One type was RHD1227A (c.1227G > A), and the other was weak D type 61 (c.28C > T). Weak D type 61 was determined on the basis of weak agglutination in the IAT procedure; it was first reported in the Chinese population [40]. The DEL (RHD1227A) variant with very low levels of D antigen detectable only by the adsorption-elution method accounts for 10% to 33% of apparent D negative phenotypes in eastern Asia [26, 40, 50]. As estimated, the maximum antigen site density per red cell was 36 and often no more than 22 [51]. In this study, RHD1227A was detected eight in 45 serologic weak D phenotypes and 28 in 117 RhD negative individuals by sequencing. Primary and second immunization of RhD negatives by RHD1227A blood have been shown to occur [10, 11, 52]. First, as the measure for improvement of transfusion safety in China, RhD negative individuals should be RHD genotyped, in order to reduce the number of immunizations of RhD negatives with RHD1227A positive blood, not identified by standard serological techniques. Second, testing of “serological weak D phenotype” donors would be of interest among donors, of higher importance in recipients. For blood samples of patients or donors with serological weak D phenotype, some are hard to determine serologically. Patients with several serologic weak D phenotypes (DFR, DV, DVI, and Weak D type 15) reported in this study have been found to develop alloanti-D [3]. Therefore, it is necessary to identify different RHD alleles and their frequencies in different populations. More practical investigations of Rh-related transfusion and obstetrics in China and other Asian populations are encouraged.

The 12 nonsynonymous variant mutations were dispersed throughout RhD protein, with no clustering at specific sites. They occurred in the intracellular, exofacial, and transmembraneous red blood cell membrane (Fig. 1). While weak D phenotypes derived mainly from amino acid substitutions in intracellular or transmembrane segments of RHD, partial D is located in extracellular portions of the RHD polypeptide [11, 35]. In this study, 12 weak D mutations were found in the intracellular and transmembrane region, except weak D101G (c.101 A > G) [25]. The possible reason is that the precise locations of the amino acid residues of RhD protein in the membrane is not yet clear; different models may predict the different locations of some amino acids [3].The substitutions may also affect the tertiary interactions and stabilization of the RhD protein. The prediction of 3D structures showed that the space conformation of the protein was disrupted in 16 serologic weak D phenotypes. These all affect the normal assembly of the tertiary structure, resulting in an activity change of the D antigen. These results indicate that bioinformatics analysis on RhD protein can give us more interpretation of missense variants of RHD gene.

The RHD gene coding region, splicing sites, partial introns, and 5′ and 3′ UTRs were detected in 45 samples with serologic weak D phenotypes in this study. The mutation sites of the 45 samples were all located in the coding region. At present, most studies on the serologic weak D phenotype are at the DNA level, and relevant available RNA information is not comprehensive. Therefore, the molecular mechanism(s) underlying serologic weak D phenotypes need to be further investigated. In addition, due to the relative scarcity of RhD negative samples in the Chinese population, especially that of serologic weak D phenotype samples, data about the overall characteristics of various ethnic groups in China are still relatively lacking at the present time. Therefore, increased specimen collection is an urgent problem that remains unresolved.


This study describes two novel and 15 rare RHD alleles by variant screening of large doubtful D phenotypes and provides a brief overview of serologic weak D phenotypes with respect to their underlying mutational mechanisms. We also applied bioinformatics analysis to predicted deleterious structural changes of serologic weak D phenotypes. These data extend our knowledge of serologic weak D phenotypes in blood donors and clinical transfusion recipients, which underlies the safety of blood transfusion and which may reduce the risk of anti-D immunization.

Availability of data and materials

Not applicable.



an indirect antiglobulin test


untranslated regions


single nucleotide polymorphism




Sorting Intolerant From Tolerant


Polymorphism Phenotypingv2


protein Variation Effect Analyzer


  1. 1.

    Storry JR, Clausen FB, Castilho L, Chen Q, Daniels G, Denomme G, Flegel WA, Gassner C, de Haas M, Hyland C, et al. International society of blood transfusion working party on red cell immunogenetics and blood group terminology: report of the Dubai. Copenhagen and Toronto meetings. Vox Sang. 2019;114:95–102.

    Article  Google Scholar 

  2. 2.

    Fichou Y, Parchure D, Gogri H, Gopalkrishnan V, Le Maréchal C, Chen J-M, Férec C, Madkaikar M, Ghosh K, Kulkarni S. Molecular basis of weak D expression in the Indian population and report of a novel, predominant variant RHD allele. Transfusion. 2018;58:1540–9.

    CAS  Article  Google Scholar 

  3. 3.

    Daniels G. Variants of RhD–current testing and clinical consequences. Br J Haematol. 2013;161:461–70.

    Article  Google Scholar 

  4. 4.

    Flegel WA. Molecular genetics and clinical applications for RH. Transfus Apher Sci. 2011;44:81–91.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Sandler SG, Flegel WA, Westhoff CM, Denomme GA, Delaney M, Keller MA, Johnson ST, Katz L, Queenan JT, Vassallo RR, Simon CD. It’s time to phase in RHD genotyping for patients with a serologic weak D phenotype. College of American Pathologists Transfusion Medicine Resource Committee Work Group. Transfusion. 2015;55:680–9.

    CAS  Article  Google Scholar 

  6. 6.

    Sandler SG, Chen LN, Flegel WA. Serological weak D phenotypes: a review and guidance for interpreting the RhD blood type using the RHD genotype. Br J Haematol. 2017;179:10–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Stegmann TC, Veldhuisen B, Bijman R, Thurik FF, Bossers B, Cheroutre G, Jonkers R, Ligthart P, de Haas M, Haer-Wigman L, van der Schoot CE. Frequency and characterization of known and novel RHD variant alleles in 37 782 Dutch D-negative pregnant women. Br J Haematol. 2016;173:469–79.

    CAS  Article  Google Scholar 

  8. 8.

    Flegel WA, Khull SR, Wagner FF. Primary anti-D immunization by weak D type 2 RBCs. Transfusion. 2000;40:428–34.

    CAS  Article  Google Scholar 

  9. 9.

    Yazer MH, Triulzi DJ. Detection of anti-D in D-recipients transfused with D + red blood cells. Transfusion. 2007;47:2197–201.

    Article  Google Scholar 

  10. 10.

    Kim KH, Kim KE, Woo KS, Han JY, Kim JM, Park KU. Primary anti-D immunization by DEL red blood cells. Korean J Lab Med. 2009;29:361–5.

    Article  Google Scholar 

  11. 11.

    Sandler SG, Flegel WA. Does transfusion of Asian-type DEL red blood cells to D- recipients cause D alloimmunization? Transfusion. 2019;59:2455–8.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Wagner FF, Flegel WA. The Rhesus Site. Transfus Med Hemother. 2014;41:357–63.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Patnaik SK, Helmberg W, Blumenfeld OO. BGMUT: NCBI dbRBC database of allelic variations of genes encoding antigens of blood group systems. Nucleic Acids Res. 2012;40:D1023–9.

    CAS  Article  Google Scholar 

  14. 14.

    Storry JR, Castilho L, Daniels G, Flegel WA, Garratty G, de Haas M, Hyland C, Lomas-Francis C, Moulds JM, Nogues N, et al. International Society of Blood Transfusion Working Party on red cell immunogenetics and blood group terminology: cancun report (2012). Vox Sang. 2014;107:90–6.

    CAS  Article  Google Scholar 

  15. 15.

    Fichou Y, Le Marechal C, Bryckaert L, Guerry C, Benech C, Dupont I, Jamet D, Ferec C, Chen JM. Variant screening of the RHD gene in a large cohort of subjects with D phenotype ambiguity: report of 17 novel rare alleles. Transfusion. 2012;52:759–64.

    CAS  Article  Google Scholar 

  16. 16.

    Van Sandt VS, Gassner C, Emonds MP, Legler TJ, Mahieu S, Kormoczi GF. RHD variants in Flanders, Belgium. Transfusion. 2015;55:1411–7.

    Article  CAS  Google Scholar 

  17. 17.

    Fichou Y, Le Marechal C, Jamet D, Bryckaert L, Ka C, Audrezet MP, Le Gac G, Dupont I, Chen JM, Ferec C. Establishment of a medium-throughput approach for the genotyping of RHD variants and report of nine novel rare alleles. Transfusion. 2013;53:1821–8.

    CAS  Article  Google Scholar 

  18. 18.

    McGowan EC, Lopez GH, Knauth CM, Liew YW, Condon JA, Ramadi L, Parsons K, Turner EM, Flower RL, Hyland CA. Diverse and novel RHD variants in Australian blood donors with a weak D phenotype: implication for transfusion management. Vox Sang. 2017;112:279–87.

    CAS  Article  Google Scholar 

  19. 19.

    Silvy M, Simon S, Gouvitsos J, Di Cristofaro J, Ferrera V, Chiaroni J, Bailly P. Weak D and DEL alleles detected by routine SNaPshot genotyping: identification of four novel RHD alleles. Transfusion. 2011;51:401–11.

    CAS  Article  Google Scholar 

  20. 20.

    Denomme GA, Wagner FF, Fernandes BJ, Li W, Flegel WA. Partial D, weak D types, and novel RHD alleles among 33,864 multiethnic patients: implications for anti-D alloimmunization and prevention. Transfusion. 2005;45:1554–60.

    CAS  Article  Google Scholar 

  21. 21.

    Chen Q, Flegel WA. Random survey for RHD alleles among D + European persons. Transfusion. 2005;45:1183–91.

    CAS  Article  Google Scholar 

  22. 22.

    Ouchari M, Jemni-Yaacoub S, Chakroun T, Abdelkefi S, Houissa B, Hmida S. RHD alleles in the Tunisian population. Asian J Transfus Sci. 2013;7:119–24.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Ba A, Beley S, Chiaroni J, Bailly P, Silvy M. RH diversity in Mali: characterization of a new haplotype RHD*DIVa/RHCE*ceTI(D2). Transfusion. 2015;55:1423–31.

    CAS  Article  Google Scholar 

  24. 24.

    Yan L, Wu J, Zhu F, Hong X, Xu X. Molecular basis of D variants in Chinese persons. Transfusion. 2007;47:471–7.

    CAS  Article  Google Scholar 

  25. 25.

    He J, Ying Y, Hong X, Xu X, Zhu F, Lv H. Molecular basis and zygosity determination of D variants including identification of four novel alleles in Chinese individuals. Transfusion. 2015;55:137–43.

    CAS  Article  Google Scholar 

  26. 26.

    Shao CP, Maas JH, Su YQ, Köhler M, Legler TJ. Molecular background of Rh D-positive, D-negative, Del and weak D phenotypes in Chinese. Vox Sang. 2002;83:156–61.

    CAS  Article  Google Scholar 

  27. 27.

    Ye L, Wang P, Gao H, Zhang J, Wang C, Li Q, Han S, Guo Z, Yang Y, Zhu Z. Partial D phenotypes and genotypes in the Chinese population. Transfusion. 2012;52:241–6.

    CAS  Article  Google Scholar 

  28. 28.

    Ji YL, Luo H, Wen JZ, Haer-Wigman L, Veldhuisen B, Wei L, Wang Z, Ligthart P, Loden-van Straaten M, Fu YS, et al. RHD genotype and zygosity analysis in the Chinese Southern Han D+, D− and D variant donors using the multiplex ligation-dependent probe amplification assay. Vox Sang. 2017;112:660–70.

    CAS  Article  Google Scholar 

  29. 29.

    Perco P, Shao CP, Mayr WR, Panzer S, Legler TJ. Testing for the D zygosity with three different methods revealed altered Rhesus boxes and a new weak D type. Transfusion. 2003;43:335–9.

    CAS  Article  Google Scholar 

  30. 30.

    Kalia N, Sharma A, Kaur M, Kamboj SS, Singh J. A comprehensive in silico analysis of non-synonymous and regulatory SNPs of human MBL2 gene. Springerplus. 2016;5:811.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    Johansson MU, Zoete V, Michielin O, Guex N. Defining and searching for structural motifs using DeepView/Swiss-PdbViewer. BMC Bioinform. 2012;13:173.

    Article  Google Scholar 

  32. 32.

    Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11:863–74.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS ONE. 2012;7:e46688.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Wagner FF, Gassner C, Müller TH, Schönitzer D, Schunter F, Flegel WA. Molecular basis of weak D phenotypes. Blood. 1999;93:385–93.

    CAS  Article  Google Scholar 

  36. 36.

    Wagner FF, Frohmajer A, Flegel WA. RHD positive haplotypes in D negative Europeans. BMC Genet. 2001;2:10.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Arnoni CP, Latini FR, Muniz JG, Gazito D, Person RD, Vendrame TA, Barreto JA, Castilho L. How do we identify RHD variants using a practical molecular approach? Transfusion. 2014;54(4):962–9.

    CAS  Article  Google Scholar 

  38. 38.

    Schmid P, von Zabern I, Scharberg EA, Wagner FF, Flegel WA. Specific amino acid substitutions cause distinct expression of JAL (RH48) and JAHK (RH53) antigens in RhCE and not in RhD. Transfusion. 2010;50:267–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Kormoczi GF, Forstemann E, Gabriel C, Mayr WR, Schonitzer D, Gassner C. Novel weak D types 31 and 32: adsorption-elution-supported D antigen analysis and comparison to prevalent weak D types. Transfusion. 2005;45:1574–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Li Q, Hou L, Guo ZH, Ye LY, Yue DQ, Zhu ZY. Molecular basis of the RHD gene in blood donors with DEL phenotypes in Shanghai. Vox Sang. 2009;97:139–46.

    CAS  Article  Google Scholar 

  41. 41.

    Fichou Y, Parchure D, Gogri H, Gopalkrishnan V, Le Maréchal C, Chen J-M, Férec C, Madkaikar M, Ghosh K, Kulkarni S. Molecular basis of weak D expression in the Indian population and report of a novel, predominant variant RHD allele. Transfusion. 2018;58:1540–9.

    CAS  Article  Google Scholar 

  42. 42.

    Rouillac C, Colin Y, Hughes-Jones NC, Beolet M, D’Ambrosio AM, Cartron JP, Le Van Kim C. Transcript analysis of D category phenotypes predicts hybrid Rh D-CE-D proteins associated with alteration of D epitopes. Blood. 1995;85:2937–44.

    CAS  Article  Google Scholar 

  43. 43.

    Wagner FF, Gassner C, Muller TH, Schonitzer D, Schunter F, Flegel WA. Three molecular structures cause rhesus D category VI phenotypes with distinct immunohematologic features. Blood. 1998;91:2157–68.

    CAS  Article  Google Scholar 

  44. 44.

    Esteban R, Montero R, Flegel WA, Wagner FF, Subirana L, Parra R, Ribera A, Nogués N. The D category VI type 4 allele is prevalent in the Spanish population. Transfusion. 2006;46:616–23.

    CAS  Article  Google Scholar 

  45. 45.

    Srivastava K, Polin H, Sheldon SL, Wagner FF, Grabmer C, Gabriel C, Denomme GA, Flegel WA. The DAU cluster: a comparative analysis of 18 RHD alleles, some forming partial D antigens. Transfusion. 2016;56:2520–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Chérif-Zahar B, Bloy C, Le Van Kim C, Blanchard D, Bailly P, Hermand P, Salmon C, Cartron JP, Colin Y. Molecular cloning and protein structure of a human blood group Rh polypeptide. Proc Natl Acad Sci USA. 1990;87:6243–7.

    Article  Google Scholar 

  47. 47.

    Luettringhaus TA, Cho D, Ryang DW, Flegel WA. An easy RHD genotyping strategy for D- East Asian persons applied to Korean blood donors. Transfusion. 2006;46:2128–37.

    CAS  Article  Google Scholar 

  48. 48.

    Ogasawara K, Sasaki K, Isa K, Tsuneyama H, Uchikawa M, Satake M, Tadokoro K. Weak D alleles in Japanese: a c960G > A silent mutation in exon 7 of the RHD gene that affects D expression. Vox Sang. 2016;110:179–84.

    CAS  Article  Google Scholar 

  49. 49.

    Liu HC, Eng HL, Yang YF, Wang YH, Lin KT, Wu HL, Lin TM. Aberrant RNA splicing in RHD 7-9 exons of DEL individuals in Taiwan: a mechanism study. Biochim Biophys Acta. 2010;1800:565–73.

    CAS  Article  Google Scholar 

  50. 50.

    Okuda H, Kawano M, Iwamoto S, Tanaka M, Seno T, Okubo Y, Kajii E. The RHD gene is highly detectable in RhD-negative Japanese donors. J Clin Invest. 1997;100:373–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Kormoczi GF, Gassner C, Shao CP, Uchikawa M, Legler TJ. A comprehensive analysis of DEL types: partial DEL individuals are prone to anti-D alloimmunization. Transfusion. 2005;45:1561–7.

    Article  Google Scholar 

  52. 52.

    Yasuda H, Ohto H, Sakuma S, Ishikawa Y. Secondary anti-D immunization by Del red blood cells. Transfusion. 2005;45:1581–4.

    Article  Google Scholar 

Download references


We acknowledge the contributions of Prof. Yujie Zhao of the China Medical University for his technical assistance with bioinformatics and blood donors for providing samples and histories.


This work was supported by Shenyang Science and Technology Foundation (F10-206-1-00; 18-014-4-50).

Author information




XZ, HF, and JL designed the experiments. XZ, GL, ZZ, XH, and LL carried out the experiments. XZ wrote the manuscript. HF and JL supervised the project. XZ, GL, CS, and YL analyzed the data, ZZ, XH, LL, and XL collected the samples. All authors provided critical comments and editorial modifications. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Hua Fan or Jianping Li.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the Liaoning Blood Center, Liaoning, China.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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: Table S1.

The calculation of allele frequencies for serologic weak D phenotype and RhD negative alleles.

Additional file 2: Table S2.

Predicted effects on RhD protein function based on SNP missense mutations.

Additional file 3: Table S3.

Predicted 3D structure of 16 serologic weak D phenotypes with the amino acid position shifts in the intracellular loop, transmembraneous domain and exofacial loop.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, X., Li, G., Zhou, Z. et al. Molecular and computational analysis of 45 samples with a serologic weak D phenotype detected among 132,479 blood donors in northeast China. J Transl Med 17, 393 (2019).

Download citation


  • RHD variant
  • Serological weak D phenotype
  • Molecular and computational analysis
  • Weak D
  • Partial D
  • DEL