RHD Genotypes in a Chinese Cohort of Pregnant Women

RHD variants in D¯ Chinese pregnant women arose difficulties in management during pregnancy. Therefore, this study aims to precisely manage D¯ pregnant women by evaluating the spectrum of RHD mutations in D¯ pregnant women and getting insight into the possible rare alleles of RHD. A total of 76 D¯ pregnant women were analyzed by performing polymerase chain reactions with sequence-specific primers (PCR-SSP), the 10 RHD exons Sanger sequencing, RHD zygosity detection, and mRNA sequencing (mRNA-seq). About 40% of alleles are variations of RHD, including RHD 1227A homozygous, RHD-CE(2-9)-D, et al. Therefore, we developed a molecular diagnostic strategy for Chinese women, and most D¯ pregnant women can be diagnosed with this simple decision tree. After RHD genotyping for D¯ pregnancy women, we eliminated at least 15% unnecessary ante- and postpartum injections of Rh immunoglobulin (RhIG). As the first pedigree study and the first functional analysis under physiological conditions, mRNA-seq revealed that c.336-1G>A mutation mainly led to the inclusion of the intron 2, which indirectly explained the D¯ phenotype in this family. We also developed a robust protocol for determining fetal RhD status from maternal plasma. All 31 fetuses were predicted as RhD positive and confirmed the RhD status after birth.

Introduction

The D antigen is one of the most immunogenic, diverse, and clinically crucial protein-based blood groups. Anti-D is still the leading cause of the hemolytic disease of the fetus and the newborn (HDN) (Flegel et al., 2009). Therefore, Rh immunoglobulin (RhIG) treatment was recommended for D¯ pregnant women for the prevention of HDN (Qureshi et al., 2014). However, RHD genotyping identified that approximately 40% of the D¯ pregnant women with weak, discrepant, or inconclusive D¯ phenotype are not candidates for RhIG treatment (Londero et al., 2020). Thus, RHD genotyping is critical for the precise management of D¯ pregnant women (Sandler et al., 2015). RHD genotyping is also crucial for avoiding alloimmunization in blood transfusion. Evidence has confirmed that some D¯ blood donors have variant RHD alleles, which might cause alloimmunization in D¯ recipients. Therefore, genotyping for RHD in D¯ donors was suggested as a routine procedure in blood centers (Krog et al., 2011; Perez-Alvarez et al., 2019).

RHD genotyping is more critical for D¯ Chinese people, for about 40% of them have RHD variants instead of RHD deletion (Colin et al., 1991; Zhang et al., 2019). However, few studies get insights into the mRNA product of RHD variations in China (Liu et al., 2010). Moreover, Most RHD variants studies rely on the bioinformatics tools in silico and plasmid construct (Fichou et al., 2015; El Wafi et al., 2017). Very little functional data under physiological conditions is available to characterize the effect of the mutation at the molecular level.

At the same time, due to the high frequency of RHD variants in D¯ Chinese pregnant women, the high-throughput method, automated real-time quantitative polymerase chain reaction (PCR), is not eligible for non-invasive prenatal testing (NIPT) for predicting fetal RhD status in D¯ Chinese people (Clausen et al., 2019). Therefore, a more flexible method needs to be developed for D¯ Chinese pregnant women.

Thus, to precisely manage D¯ pregnant women, this study elucidated the molecular basis of D¯ pregnant women in Shaoxing and established a simple molecular diagnostic strategy for Chinese D¯ pregnant women. In addition, we investigate the first family study and the first mRNA sequencing (mRNA-seq) for c.336-1G>A alleles under physiological conditions to analyze the mechanism of intronic mutation. We also developed a robust protocol for determining fetal RhD status from maternal plasma.

Materials and Methods

Test Subjects

This study was approved by the Shaoxing Maternal and Child Health Hospital (approval n. 277). A total of 76 RhD-negative pregnant women underwent polymerase chain reactions with sequence-specific primers (PCR-SSP), 10 RHD exons Sanger sequencing, RHD zygosity detection, and mRNA-seq between November 2018 and June 2020.

Test Method

DNA extraction

Nucleic acid extraction or purification was performed using commercial reagents (Tianjin Super Biotechnology Developing Co., Ltd., Tianjin, China; Lot No.: 201808001). Sample DNA was extracted according to the instructions of the DNA Extraction Kit, and the DNA concentration was measured and diluted to a final concentration of approximately 30 ng/μL (A260/280: 1.6–1.9). If the samples could not be tested immediately, they were stored at -20°C.

Polymerase Chain Reactions with Sequence-Specific Primers

The Human Erythrocyte RHD Genotyping Kit (PCR-SSP) (Tianjin Super Biotechnology Developing Co., Ltd., Tianjin, China; Lot No. 190830001) was used to detect eight common RHD genotypes: RHD-positive, RHD deletion, RHD-CE(2-9)-D, DVa (Hus), DVI III, weak D15, DEL RHD 1227A homozygous, and DEL RHD 1227A heterozygous. Reaction parameters were set according to the Human Erythrocyte RHD Genotyping Kit instructions, and amplification was performed on a Hema 9600 Gradient Thermal Cycler (Hangzhou Bioer Technology Co., Ltd., Zhejiang, China). The method is briefly described as follows. The PCR amplification system consisted of dNTP-Buffer working solution 80 μL, Taq enzyme 0.8 μL, DNA 10 μL, and 10 μL of the above-mixed solution was added into each of the eight wells coated with primers. PCR amplification procedure: pre-denaturation at 96°C for 2 min; 96°C for 20 s, 68°C for 1 min, 5 cycles; 96°C 20 s, 65°C 50 s, 72°C 45 s, 10 cycles; 96°C 20 s, 62°C 50 s, 72°C 45 s, 18 cycles; Finally, it was extended at 72°C for 5 min. PCR products were visualized on a gel imager after 2.5% agarose gel electrophoresis, and the results were interpreted according to the instructions provided with the kit.

Zygosity Detection and Sanger Sequencing

RHD zygosity detection and the 10 RHD exons Sanger sequencing were commissioned to Tianjin Super Biotechnology Developing Co. (China). RHD zygosity was assessed by the PCR method. Two pairs of primers were designed to amplify the hybrid Rhesus box of RHD- (2,700 bp) and internal control (1009 bp), respectively. If a hybrid Rhesus box is detected, there is a complete deletion of the RHD gene. The Sanger dideoxy method performed direct Sequencing of all the 10 RHD exons and flanking intron regions. Sequence analysis was performed by DNAMAN v9 (Lynnon Biosoft Co., United States) and ChromasPro v1.2 (Technelysium Pty. Ltd., Australia), and the reference allele was RHD*01 (NG_007494.1)

RHD mRNA-seq

RHD mRNA-seq was commissioned to Beijing Beikang Medical Laboratory Co. (China). First, the total RNA of whole blood was extracted, and the mRNA with polyA tail was enriched by Oligo (dT) magnetic column. Subsequently, the obtained mRNA was randomly interrupted with divalent cations in NEB Fragmentation Buffer, and the library was constructed by the standard NEB library construction method. NEBNext® UltraTM RNA Library Prep Kit for Illumina® was used for the library construction. The first strand of cDNA was synthesized in the M-MuLV reverse transcriptase system with the fragmented mRNA as templates and random oligonucleotides as primers, and then the RNA strands were degraded with RNaseH.

Moreover, dNTPs were used to synthesize the second strand of cDNA under the DNA polymerase Ⅰ system. After purifying the double-stranded cDNA, end-repair, A-tailing, and ligation of sequencing adapter were performed. Next, the 250–300 bp cDNA was screened with AMPure XP beads, PCR amplification was performed, and the PCR product was purified again with AMPure XP beads to obtain the library. Qubit® 2.0 Fluorometer (Thermo Fisher Scientific, Walsham, United States) was used for preliminary quantification after the library was constructed. Then the library was diluted to 1.5 ng/μL, and Agilent 2,100 bioanalyzer was used to detect the insert size of the library. Then qRT-PCR was used to measure the effective concentration of the library whose insert size was in line with the expectation. Accurate quantification by qRT-PCR (the effective concentration of the library was higher than 2 nM) was to ensure the quality of the library. Then the library was pooled and finally analyzed on the Novaseq 6,000 sequencer.

Non-Invasive Fetal RHD Genotyping

Cell-free fetal DNA (cffDNA) extraction from maternal plasma (200 μL) was performed using the BGISP-300 (BGI, Shenzhen, China) and the Nucleic Acid Extraction (BGI, Shenzhen, China) kits. First, the DNA amplification products were quantified on a Qubit® 2.0 Fluorometer (Thermo Fisher Scientific, Walsham, United States) using the QubitTM dsDNA HS Assay kit (Thermo Fisher Scientific, Walsham, United States). Then use the Human Erythrocyte RHD Genotyping Kit (PCR-SSP) (Tianjin Super Biotechnology Developing Co., Ltd., Tianjin, China; Lot No. 190830001) to detect exon 1, 5, 6, 7, 9. The PCR amplification system consisted of dNTP-Buffer working solution 80 μL, Taq enzyme 0.8 μL, DNA 10 μL, and 10 μL of the above-mixed solution was added into each of the eight wells coated with primers. PCR amplification procedure: pre-denaturation at 96°C for 2 min; 96°C for 20 s, 68°C for 1 min, 8 cycles; 96°C 20 s, 65°C 50 s, 72°C 45 s, 11 cycles; 96°C 20 s, 62°C 50 s, 72°C 45 s, 30 cycles; Finally, it was extended at 72°C for 5 min. PCR amplification products were imaged on a gel imager after 2.5% agarose gel electrophoresis, and the results were interpreted RhD+ when any of the exons (1, 5, 6, 7, 9) were positive for RHD deletion women, and any of the exons (5, 6, 7, 9) were positive for RHD-CE(2-9)-D women.

Results

RHD Genotypes in D¯ Pregnancy Women

The PCR-SSP results indicated 63.2% (48/76) RHD deletion homozygous, 15.8% (12/76) RHD 1227A homozygous, 14.5% (11/76) RHD-CE(2-9)-D, 2.6% (2/76) RHD VI III, and 3.9% (3/76) genotypes could not be identified by PCR-SSP (Table 1).

TABLE 1

TABLE 1. Genotyping results of RHD by PCR-SSP method in 76 pregnant women in Shaoxing.

Three Cases of RHD Point Mutation

The first case was RHD*10.08 partial D type (exon 3 c.340C>T) (Figure 1, Sample a). The second was RHD*01EL.02 type (exon 1 c.3G>A) (Figure 1, Sample b), and a mutation at the splicing site was found as RHD*01N.25 type (intron 2 c.336-1G>A) (Figure 1, Sample c) in the third case.

FIGURE 1

FIGURE 1. Three cases in which the genotype could not be confirmed by PCR-SSP method were found to have point mutations after RHD Sanger sequencing: Sample a. RHD*10.08 partial D type (exon 3 c.340C>T); Sample b. RHD*01EL.02 type (exon 1 c.3G>A); Sample c. RHD*01N.25 type (intron 2 c.336-1G>A).

Family Study of 336-1G>A Mutation

Two other RhD-negative family members (III-1, III-3) were identified in the family of the patient with RHD*01N.25 (intron 2 c.336-1G>A) (Figure 2). PCR-SSP and RHD zygosity analysis confirmed the proband's mother (II-4) and uncle-in-law (II-1) had the hybrid Rhesus box detected, with RHD+/RHD- genotype (Supplementary Figures S1–S2). In comparison, proband's father (II-3), aunt (II-2) (Supplementary Figure S3) have the same genotype RHD+/336-1G>A, which speculated the zygosity of the proband and her cousin (III-1) as RHD-/336-1G>A. RhD serological results and RHD genotype of this family are shown in Table 2.

FIGURE 2

FIGURE 2. Family pedigree of the patient with RHD*01N.25 type (intron 2 c.336-1G>A), with ABO blood group, RhD antigen phenotype, and genotype (some family members were tested), labeled sequentially below each case in the pedigree.

TABLE 2

TABLE 2. RhD serological results and RHD genotype of the family of the case with RHD*01N.25 (intron 2 c.336-1G>A).

mRNA Sequence Analysis

The sequencing results (Figure 3) verified 336-1G>A mutation site and revealed the multiple splicing products of different lengths with the intron 2 residue. In addition, the number of reads in intron 2 was significantly higher than other introns (Supplementary Figures S4–S6), suggesting that intron 2 was improperly spliced after c.336-1G>A mutation.

FIGURE 3

FIGURE 3. mRNA-seq results using the Integrative Genomics Viewer in a web browser (IGV-Web app version 1.6.3) (Robinson et al., 2011), showing multiple splicing products with different lengths (31 bp - 197 bp) of intron 2 sequences residue which are much longer than the 28 bp predicted by Fichou (Fichou et al., 2015), verifying that a cryptic site upstream of the constitutive acceptor site was activated in the presence of c.336-1G>A (green). To avoid the possibility of paralogs of transcript isotype, we BLAT the read's sequences in different regions and confirm that 197 bp intron 2 sequences residue near the mutation site is unique and has no isotypes.

Non-invasive Fetal RHD Genotyping

Fetal RhD status in 28 cases of maternal homozygous RHD deletion and 3 cases of maternal RHD-CE(2-9)-D were predicted (Table 3). All predicted results were consistent with the newborn RhD status. Two cases showed RhD negative in the first samples but confirmed positive in the second samples.

TABLE 3

TABLE 3. Fetus RhD status predicted results of 33 samples of NIPT plasma in 31 pregnant women.

Discussion

RHD Genotypes of D¯ Pregnant Women in Shaoxing and the Molecular Diagnostic Strategy for Chinese D¯ Pregnant Women

RhD-negative frequencies show vast racial differences (Xu et al., 2003; Flegel et al., 2009; Polin et al., 2009; Cruz et al., 2012; Crottet et al., 2014; Perez-Alvarez et al., 2019) (Table 4). Thus, for diverse ethnic populations, it is necessary to adopt different RHD genotyping strategies to the spectrum of prevalent alleles. In our study, the main genotypes of RHD in D¯ pregnant women in Shaoxing were RHD deletion, RHD 1227A homozygous type, and RHD-CE(2-9)-D type. The profile of our survey was in concordance with that of another study in China (Zhang et al., 2019), indicating minor differences in the RHD genotypes of D¯ pregnant women across different regions in China. Using the PCR-SSP method, designed for eight common Chinese RHD genotypes: RHD-positive, RHD deletion, RHD-CE(2-9)-D, DVa (Hus), DVI III, weak D15, DEL RHD 1227A homozygous, and DEL RHD 1227A heterozygous, we found that about 96.1% of D¯ women can identify RHD genotype, and the rate rises to above 99% when adding the 10 RHD exons sequencing. Therefore, we established a simple molecular diagnostic strategy for Chinese D¯ pregnant women (Figure 4).

TABLE 4

TABLE 4. Frequency of RHD alleles in serologic RhD-negative blood donors.

FIGURE 4

FIGURE 4. A simple decision tree for the molecular diagnosis of D^– pregnancy women in China.

The Precise Guide of RhIG use for D¯ Pregnant Women

The occurrence of alloimmunity in D¯ pregnant women is related to several factors, such as maternal RHD genotype and RhD antigen epitope (Flegel, 2011); quantity of D-positive fetal red blood cells entering the body (Woodrow, 1971); and other factors such as ABO blood type, maternal HLA class, or fetal sex (Hadley and Soothill, 2002). However, the main factor is maternal RHD genotype and RhD antigen epitope. Therefore, our study guides RhIG use according to the different RHD genotypes in D¯ pregnant women (Table 5).

TABLE 5

TABLE 5. RHD alleles and guidance for managing transfusion or RhIG administration.

The RHD (c.1227G>A) allele, named the DEL type or Asian DEL type, is the most prevalent DEL allele in Asians. Asian DEL is always mistyped as D-negative by routine serological assays (Kim et al., 2020). Therefore, RHD genotyping is the gold standard for detecting DEL (Nuchnoi et al., 2014). The RHD (c.1227G>A) should be managed as D-positive for RhIG administration or selection of blood components for transfusion (Chun et al., 2020). In 2015, a study revealed that erythrocytes carrying the 1227A mutation might express very low “normal” D antigen levels (Fichou et al., 2015). This theory explains why “Asian DEL types” do not cause hemolytic disease in newborns and fetuses. Thus, at least 15% of D¯ pregnant women in China, typed as the RHD (c.1227G>A) allele, will benefit from RHD genotyping and be free from unnecessary RhIG use. At the same time, to avoid transfusion of DEL RBC units to D¯ recipients, RBC units carrying c.1227G>A mutation should be transferred into the D+ pool (Flegel et al., 2009).

RHD-CE(2-9)-D is a hybrid allele that replaces exons 2 to 9 by RHCE. Lacking D antigen expression, RBCs of RHD-CE(2-9)-D will not risk D¯ recipients (Flegel et al., 2009). At the same time, women with RHD-CE(2-9)-D should be managed as D¯ and be given RhIG for D immunoprophylaxis in pregnancy (Table 5).

DVI is the most common partial D that produces anti-D, and newborns born to DVI mothers with anti-D may develop hemolytic disease. DVI III is named a D-Ce(3-6)-D hybrid and DVI III erythrocytes carry relatively high RhD antigen densities. Therefore, DVI III recipients should be transfused with RhD negative blood, while DVI III donor should be managed as D⁺ (Wagner et al., 1998).

Our study characterized two point mutations in exons: c.340C>T and c.3G>A. According to the guideline (Flegel et al., 2020), both should be managed as D¯ and given RhIG for D immunoprophylaxis. More interesting, a recent study showed that weak-D and Asia-type DEL alleles’ coexistence would completely express the D-antigen (Chun et al., 2020), which means free from RhIG administration.

For our study’s 336-1G>A mutation, the mRNA products of the proband showed long intronic segment retention, which will induce DEL or D-negative phenotypes (Fichou et al., 2015). The serological study for family members reveals D¯ phenotypes in RHD-/336-1G>A heterozygote and D positive in RHD+/336-1G>A heterozygote (Table 2). However, lacking the detailed evidence for no D-antigen in red blood cells, we still suggested these individuals with RHD-/336-1G>A heterozygote should be managed as RhD positive donors in transfusion in an abundance of caution.

Thus, in our study, RHD genotyping for D¯ pregnant women eliminates at least 15% unnecessary ante- and postpartum injections of RhIG and gave a detailed guide of RhIG use for these D¯ pregnant women.

c.336-1G>A Mutation Analysis and Family Study

More than 300 variant alleles of RHD have been reported (Kawano et al., 1998; Ye et al., 2007; Fichou et al., 2015; Chen et al., 2016; Ogasawara et al., 2016; El Wafi et al., 2017; Chun et al., 2018; Raud et al., 2019), including single-nucleotide polymorphisms (SNPs), small or large fragment deletions, gene rearrangements, and complete RHD deletions. The intron mutation disrupts a constitutive splice site, resulting in improper retention of an intron or activation of a cryptic splice site in the vicinity of the mutant. Thus, this mutation will induce DEL or D-negative phenotypes by minute to no expression of the D antigen at the surface of RBCs. As many intronic mutated alleles of RHD were identified in Sequencing, several bioinformatics tools were used to predict the defect resulting from a genetic variation in silico. However, very few functional data are available to confirm these predictions. Furthermore, lacking fresh blood samples for RNA extraction, many functional studies mainly relied on recombinant plasmids. Hence, clinical transcriptome data and protein analysis are urgently needed to investigate genotype-phenotype mechanisms (Liu et al., 2010; Fichou et al., 2015; El Wafi et al., 2017).

We found a family with 336-1G>A intron mutation in our study. The mutation site of 336-1G>A is an intron 2 acceptor. The 336-1G>A mutation was first reported in Korea in 2005 (Kim et al., 2005) and first identified in Chinese blood donors in 2009 (Ye et al., 2009), named as the RHD*01N.25 type (IVS2-1G>A) by ISBT. In 2015, The first functional analysis by plasmid recombination experiments revealed that the cryptic splicing site of c.336-1G>A was activated, and mRNA product was either increased by 28 base pairs (bp) (from intron 2) or decreased by 21 bp (exon 3) (Fichou et al., 2015).

Our serological D-negative proband was an RHD-deletion/336-1G>A heterozygote in our study. SpliceAI showed a 94% decrease in the probability of acting as a splice acceptor in the 336-1G>A site. Our mRNA products of the 336-1G>A mutation under physiological conditions confirmed the prediction. A cryptic site upstream of the constitutive acceptor site was activated in the presence of c.336-1G>A, which resulted in insertion of an intronic sequence within the mature transcript. The retentional intronic segments (31 bp–197 bp) are much longer in our study than the 28 bp predicted by Fichou (Fichou et al., 2015), which suggests that the mRNA product under physiological conditions is different from the plasmid recombination experiments. This mature RHD transcript with intronic segments might be translated into a protein with an additional 10–66 amino-acid sequence, which may alter the structure of the second extracellular loop of RhD. We didn't identify the active cryptic splice site in intron 2. However, we supposed the intronic retention inducing in-frame shift because the intron 3 seems correctly splicing in our study. All reported RHD-/336-1G>A cases are completely serological D-negative, which means no antigens can be detected (Kim et al., 2005; Ye et al., 2009; Fichou et al., 2015). Our results first revealed that the 336-1G>A intron mutation induces the retention of intron 2, which will generate improper peptides.

Fetal RHD Genotyping in Maternal Plasma

Prenatal detection of the fetal RhD status is helpful to assess the risk of hemolytic disease in fetus. Non-invasive prenatal testing for predicting fetal RhD status is available in many countries (Finning et al., 2008; Akolekar et al., 2011; Wikman et al., 2012; Chitty et al., 2014; Clausen et al., 2014). The standard quantitative PCR method is based on the assumption of homozygous RHD deletion in D¯ pregnant women. RHD variants in D¯ women will yield high false-positive results (Yang et al., 2019). Therefore, Non-invasive fetal RHD genotyping is only feasible for D¯ women with total or partial deletion of the RHD (Clausen et al., 2019; Yang et al., 2019).

Unlike the high-throughput quantitative PCR method widely used in European countries, we use the PCR byproduct from routine NIPT for aneuploidy as DNA template, then amplified by PCR-SSP method. Because most Chinese pregnant women choose NIPT for aneuploidy after 12 gestational weeks, our NIPT + PCR-SSP protocol for D¯ pregnant women is cost-effective and suitable for the small number of D¯ samples in China. In addition, the results of routine NIPT for aneuploidy can provide good quality control data, such as fetal fraction concentration.

Many research chooses the exons 5, 7, 10 as the targeted exons to predict fetal RhD status (Finning et al., 2008; Akolekar et al., 2011; Chitty et al., 2014; Clausen et al., 2014), while Wikman (Wikman et al., 2012) choose exon 4 as target exon. However, our study found that any one of exon 1, 5, 6, 7, and 9 positive can predict fetal D+ status. Furthermore, all our predicted results were consistent with the newborn RhD status.

Our study also emphasized the importance of resampling for the first negative results. Two cases showed negative results in the first sample but were later confirmed as positive in the second sample, which was likely due to the lower fetal fraction content or PCR bias (Sabina and Leamon, 2015). Thus, we argued the second sampling after a few weeks to confirm the first negative results (Table 3).

No fetuses were predicted as RhD negative in our study. Asian D-negative pregnant women have a 96% possibility of having a D-positive fetus (Chun et al., 2020). Predicting fetal RhD status is not cost-effective in China, whatever methods are used.

Limitation of Our Study

To our knowledge, our study used the mRNA-seq to analyze the splicing products of the 336-1G>A under physiological conditions at the first time. However, there has a significant limitation by using an inappropriate sample type, whole peripheral blood, instead of enriched reticulocytes. Hence the reads numbers are relatively small in mRNA-seq. The active cryptic site upstream of the constitutive acceptor site is not identified in our study.

Conclusion

In conclusion, we established a pathway to precisely manage Chinese D¯ pregnant women, including using the PCR-SSP method and Sanger sequencing to provide precise RhIG administration and promote transfusion safety. The robust protocol for determining fetal RhD status from maternal plasma is suitable for Chinese D¯ women. In addition, the detailed study on rare alleles will help us investigate the molecular mechanisms between genotype and phenotype of RHD variations.

Data Availability Statement

1) The data presented in the study are deposited in the China National GeneBank DataBase (CNGBdb) repository, accession number CNP0002279, could be founded with https://db.cngb.org/cnsa/project/CNP0002279/reviewlink/. 2) The genetic data could not be publicly available because of (Regulation of the People’s Republic of China on the Administration of Human Genetic Resources), any inquries could contact the corresponding author for details.

Ethics Statement

The studies involving human participants were reviewed and approved by the Shaoxing Maternal and Child Health Hospital (approval n. 277). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

This work was supported by the Shaoxing Public Welfare Technology Application Research Program Project (2018C30039).

Conflict of Interest

Authors XP, ZZ, LA, and YZ were employed by the company Tianjin Super Biotechnology Developing Co., Ltd. Authors JP and MW were employed by the company Zhejiang Biosan Biotechnology Co., Ltd. Author JZ was employed by the company BGI-Shenzhen.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank Professor Choy Kwongwai of the Chinese University of Hong Kong and Zhang Yanyan of Shenzhen BGI Co., Ltd for their guidance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.752485/full#supplementary-material

References

Akolekar, R., Finning, K., Kuppusamy, R., Daniels, G., and Nicolaides, K. H. (2011). Fetal RHD Genotyping in Maternal Plasma at 11-13 Weeks of Gestation. Fetal Diagn. Ther. 29 (4), 301–306. doi:10.1159/000322959

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D.-P., Sun, C.-F., Ning, H.-C., Wang, W.-T., and Tseng, C.-P. (2016). Comprehensive Analysis ofRHDsplicing Transcripts Reveals the Molecular Basis for the Weak Anti-D Reactivity of Del-Red Blood Cells. Transfus. Med 26 (2), 123–129. doi:10.1111/tme.12270

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitty, L. S., Finning, K., Wade, A., Soothill, P., Martin, B., Oxenford, K., et al. (2014). Diagnostic Accuracy of Routine Antenatal Determination of Fetal RHD Status across Gestation: Population Based Cohort Study. BMJ 349, g5243. doi:10.1136/bmj.g5243

PubMed Abstract | CrossRef Full Text | Google Scholar

Chun, S., Yun, J. W., Park, G., and Cho, D. (2018). The Synonymous Nucleotide Substitution RHD 1056C>G Alters mRNA Splicing Associated with Serologically Weak D Phenotype. J. Clin. Lab. Anal. 32 (4), e22330. doi:10.1002/jcla.22330

PubMed Abstract | CrossRef Full Text | Google Scholar

Chun, S., Kim, H., Yun, J. W., Yu, H., Seo, J. Y., and Cho, D. (2020). RHD Genotyping Is Recommended for All Patients with Serological Weak-D Phenotypes in Asian Populations - Cases with Coexistence of Weak-D and Asia Type DEL Alleles Results in Complete Expression of D-Antigen. Transfus. Apher. Sci. 59 (4), 102807. doi:10.1016/j.transci.2020.102807

PubMed Abstract | CrossRef Full Text | Google Scholar

Banch Clausen, F., Steffensen, R., Christiansen, M., Rudby, M., Jakobsen, M. A., Jakobsen, T. R., et al. (2014). Routine Noninvasive Prenatal Screening for fetalRHDin Plasma of RhD-Negative Pregnant Women-2 Years of Screening Experience from Denmark. Prenat Diagn. 34 (10), 1000–1005. doi:10.1002/pd.4419

PubMed Abstract | CrossRef Full Text | Google Scholar

Clausen, F. B., Rieneck, K., Krog, G. R., Bundgaard, B. S., and Dziegiel, M. H. (2019). “Noninvasive Antenatal Screening for Fetal RHD in RhD Negative Women to Guide Targeted Anti-D Prophylaxis,” in Prenatal Diagnosis. Editor B. Levy. (New York, NY: Springer), 347–359. doi:10.1007/978-1-4939-8889-1_23

PubMed Abstract | CrossRef Full Text | Google Scholar

Colin, Y., Chérif-Zahar, B., Le Van Kim, C., Raynal, V., Van Huffel, V., and Cartron, J. (1991). Genetic Basis of the RhD-Positive and RhD-Negative Blood Group Polymorphism as Determined by Southern Analysis. Blood 78 (10), 2747–2752. doi:10.1182/blood.v78.10.2747.bloodjournal78102747

PubMed Abstract | CrossRef Full Text | Google Scholar

Crottet, S. L., Henny, C., Meyer, S., Still, F., Stolz, M., Gottschalk, J., et al. (2014). Implementation of a Mandatory Donor RHD Screening in Switzerland. Transfus. Apher. Sci. 50 (2), 169–174. doi:10.1016/j.transci.2014.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz, B. R., Chiba, A. K., Moritz, E., and Bordin, J. O. (2012). RHD Alleles in Brazilian Blood Donors with Weak D or D-Negative Phenotypes. Transfus. Med. 22 (2), 84–89. doi:10.1111/j.1365-3148.2011.01129.x

PubMed Abstract | CrossRef Full Text | Google Scholar

El Wafi, M., El Housse, H., Zaid, N., Zouine, S., Nourichafi, N., Bouisk, K., et al. (2017). Novel Intronic RHD Variants Identified in Serologically D-Negative Blood Donors. Vox Sang 112 (8), 796–802. doi:10.1111/vox.12570

PubMed Abstract | CrossRef Full Text | Google Scholar

Fichou, Y., Gehannin, P., Corre, M., Le Guern, A., Le Maréchal, C., Le Gac, G., et al. (2015). Extensive Functional Analyses of RHD Splice Site Variants: Insights into the Potential Role of Splicing in the Physiology of Rh. Transfusion 55 (6 Pt 2), 1432–1443. doi:10.1111/trf.13083

PubMed Abstract | CrossRef Full Text | Google Scholar

Finning, K., Martin, P., Summers, J., Massey, E., Poole, G., and Daniels, G. (2008). Effect of High Throughput RHD Typing of Fetal DNA in Maternal Plasma on Use of Anti-RhD Immunoglobulin in RhD Negative Pregnant Women: Prospective Feasibility Study. BMJ 336 (7648), 816–818. doi:10.1136/bmj.39518.463206.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Flegel, W. A., von Zabern, I., and Wagner, F. F. (2009). Six Years' Experience performingRHDgenotyping to Confirm D− Red Blood Cell Units in Germany for Preventing Anti-D Immunizations. Transfusion 49 (3), 465–471. doi:10.1111/j.1537-2995.2008.01975.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Flegel, W. A., Denomme, G. A., Queenan, J. T., Johnson, S. T., Keller, M. A., Westhoff, C. M., et al. (2020). It's Time to Phase Out "serologic Weak D Phenotype" and Resolve D Types with RHD Genotyping Including Weak D Type 4. Transfusion 60 (4), 855–859. doi:10.1111/trf.15741

PubMed Abstract | CrossRef Full Text | Google Scholar

Flegel, W. A. (2011). Molecular Genetics and Clinical Applications for RH. Transfus. Apher. Sci. 44 (1), 81–91. doi:10.1016/j.transci.2010.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hadley, A., and Soothill, P. (2002). Alloimmune Disorders of Pregnancy. Anaemia, Thrombocytopenia and Neutropenia in the Fetus and Newborn. Transfus. Med. 12 (5), 331–332. doi:10.1046/j.1365-3148.2002.00397.x

CrossRef Full Text | Google Scholar

Kawano, M., Iwamoto, S., Okuda, H., Fukuda, S., Hasegawa, N., and Kajii, E. (1998). A Splicing Mutation of the RHAG Gene Associated with the Rhnull Phenotype. Ann. Hum. Genet. 62 (Pt 2), 107–113. doi:10.1046/j.1469-1809.1998.6220107.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. Y., Kim, S. Y., Kim, C.-A., Yon, G. S., and Park, S. S. (2005). Molecular Characterization of D- Korean Persons: Development of a Diagnostic Strategy. Transfusion 45 (3), 345–352. doi:10.1111/j.1537-2995.2005.04311.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, T. Y., Hong, Y. J., Kim, M. J., Kim, H., Kim, T. S., Park, J. S., et al. (2020). Recommendations Regarding Practical DEL Typing Strategies for Serologically D-Negative Asian Donors. Transfus. Med. Hemother 47 (1), 88–93. doi:10.1159/000500098

PubMed Abstract | CrossRef Full Text | Google Scholar

Krog, G. R., Clausen, F. B., Berkowicz, A., Jørgensen, L., Rieneck, K., Nielsen, L. K., et al. (2011). Is Current Serologic RhD Typing of Blood Donors Sufficient for Avoiding Immunization of Recipients? (CME). Transfusion 51 (11), 2278–2285. doi:10.1111/j.1537-2995.2011.03156.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H.-C., Eng, H.-L., Yang, Y.-F., Wang, Y.-H., Lin, K.-T., Wu, H.-L., et al. (2010). Aberrant RNA Splicing in RHD 7-9 Exons of DEL Individuals in Taiwan: a Mechanism Study. Biochim. Biophys. Acta (Bba) - Gen. Subjects 1800 (6), 565–573. doi:10.1016/j.bbagen.2010.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Londero, D., Monge, J., and Hellberg, A. (2020). A Multi‐Centre Study on the Performance of the Molecular Genotyping Platform ID RHD XT for Resolving Serological Weak RhD Phenotype in Routine Clinical Practice. Vox Sang 115 (3), 241–248. doi:10.1111/vox.12886

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuchnoi, P., Thongbus, J., Srisarin, A., Kerdpin, U., and Prachayasittikul, V. (2014). Clinical and Laboratory Update on the DEL Variant. Lab. Med. 45 (4), 285–290. doi:10.1309/lmtuz00o7vftgceb

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogasawara, K., Sasaki, K., Isa, K., Tsuneyama, H., Uchikawa, M., Satake, M., et al. (2016). Weak D Alleles in Japanese: a c.960G>A Silent Mutation in Exon 7 of theRHDgene that Affects D Expression. Vox Sang 110 (2), 179–184. doi:10.1111/vox.12322

PubMed Abstract | CrossRef Full Text | Google Scholar

Perez-Alvarez, I., Hayes, C., Hailemariam, T., Shin, E., Hutchinson, T., and Klapper, E. (2019). RHD Genotyping of Serologic RhD-Negative Blood Donors in a Hospital-Based Blood Donor center. Transfusion 59 (7), 2422–2428. doi:10.1111/trf.15325

PubMed Abstract | CrossRef Full Text | Google Scholar

Polin, H., Danzer, M., Gaszner, W., Broda, D., St-Louis, M., Pröll, J., et al. (2009). Identification ofRHDalleles with the Potential of Anti-D Immunization Among Seemingly D− Blood Donors in Upper Austria. Transfusion 49 (4), 676–681. doi:10.1111/j.1537-2995.2008.02046.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Qureshi, H., Massey, E., Kirwan, D., Davies, T., Robson, S., White, J., et al. (2014). BCSH Guideline for the Use of Anti-D Immunoglobulin for the Prevention of Haemolytic Disease of the Fetus and Newborn. Transfus. Med 24 (1), 8–20. doi:10.1111/tme.12091

PubMed Abstract | CrossRef Full Text | Google Scholar

Raud, L., Ka, C., Gourlaouen, I., Callebaut, I., Férec, C., Le Gac, G., et al. (2019). Functional Analysis of NovelRHDvariants: Splicing Disruption Is Likely to Be a Common Mechanism of Variant D Phenotype. Transfusion 59 (4), 1367–1375. doi:10.1111/trf.15210

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, J. T., Thorvaldsdóttir, H., Winckler, W., Guttman, M., Lander, E. S., Getz, G., et al. (2011). Integrative Genomics Viewer. Nat. Biotechnol. 29 (1), 24–26. doi:10.1038/nbt.1754

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabina, J., and Leamon, J. H. (2015). Bias in Whole Genome Amplification: Causes and Considerations. Methods Mol. Biol. 1347, 15–41. doi:10.1007/978-1-4939-2990-0_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandler, S. G., Flegel, W. A., Westhoff, C. M., Denomme, G. A., Delaney, M., Keller, M. A., et al. (2015). It's Time to Phase inRHDgenotyping for Patients with a Serologic Weak D Phenotype. Transfusion 55 (3), 680–689. doi:10.1111/trf.12941

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, F. F., Gassner, C., Müller, T. H., Schönitzer, D., Schunter, F., and Flegel, W. A. (1998). Three Molecular Structures Cause Rhesus D Category VI Phenotypes with Distinct Immunohematologic Features. Blood 91 (6), 2157–2168. doi:10.1182/blood.v91.6.2157

PubMed Abstract | CrossRef Full Text | Google Scholar

Wikman, A. T., Tiblad, E., Karlsson, A., Olsson, M. L., Westgren, M., and Reilly, M. (2012). Noninvasive Single-Exon Fetal RHD Determination in a Routine Screening Program in Early Pregnancy. Obstet. Gynecol. 120 (2 Pt 1), 227–234. doi:10.1097/aog.0b013e31825d33d9

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodrow, J. G. (1971). Rh-Immunisation and its Prevention. Nord Med. 85 (22), 704–705.

PubMed Abstract | Google Scholar

Xu, Q., Zhang, J., Wang, Q., Zhang, S., and Si, G. (2003). RHD Gene Polymorphism Among RhD-Negative Han Chinese. Chin. Med. J. (Engl) 116 (10), 1539–1543. doi:10.3760/cma.j.issn.0366-6999.2003.10.123

CrossRef Full Text | Google Scholar

Yang, H., Llewellyn, A., Walker, R., Harden, M., Saramago, P., Griffin, S., et al. (2019). High-throughput, Non-invasive Prenatal Testing for Fetal Rhesus D Status in RhD-Negative Women: a Systematic Review and Meta-Analysis. BMC Med. 17 (1)–37. doi:10.1186/s12916-019-1254-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, L.-Y., Guo, Z.-H., Li, Q., and Zhu, Z.-Y. (2007). Molecular and Family Analyses Revealed Two Novel RHD Alleles in a Survey of a Chinese RhD-Negative Population. Vox Sang 92 (3), 242–246. doi:10.1111/j.1423-0410.2006.00886.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, L., Yue, D., Wo, D., Ding, X., Guo, S., Li, Q., et al. (2009). Molecular Bases of unexpressedRHDalleles in Chinese D− Persons. Transfusion 49 (8), 1655–1660. doi:10.1111/j.1537-2995.2009.02181.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Li, G., Zhou, Z., Shao, C., Huang, X., Li, L., et al. (2019). Molecular and Computational Analysis of 45 Samples with a Serologic Weak D Phenotype Detected Among 132,479 Blood Donors in Northeast China. J. translational Med. 17 (1), 393. doi:10.1186/s12967-019-02134-9

CrossRef Full Text | Google Scholar