Comparative validation of computer programs for haplotype frequency estimation from donor registry data.

Estimation of human leukocyte antigen (HLA) haplotype frequencies from unrelated stem cell donor registries presents a challenge because of large sample sizes and heterogeneity of HLA typing data. For the 14th International HLA and Immunogenetics Workshop, five bioinformatics groups initiated the 'Registry Diversity Component' aiming to cross-validate and improve current haplotype estimation tools. Five datasets were derived from different donor registries and then used as input for five different computer programs for haplotype frequency estimation. Because of issues related to heterogeneity and complexity of HLA typing data identified in the initial phase, the same five implementations, and two new ones, were used on simulated datasets in a controlled experiment where the correct results were known a priori. These datasets contained various fractions of missing HLA-DR modeled after European haplotype frequencies. We measured the contribution of sampling fluctuation and estimation error to the deviation of the frequencies from their true values, finding equivalent contributions of each for the chosen samples. Because of patient-directed activities, selective prospective typing strategies and the variety and evolution of typing technology, some donors have more complete and better HLA data. In this setting, we show that restricting estimation to fully typed individuals introduces biases that could be overcome by including all donors in frequency estimation. Our study underlines the importance of critical review and validation of tools in registry-related activity and provides a sustainable framework for validating the computational tools used. Accurate frequencies are essential for match prediction to improve registry operations and to help more patients identify suitably matched donors.

16(th) IHIW: global analysis of registry HLA haplotypes from 20 million individuals: report from the IHIW Registry Diversity Group.

This project has the goal to validate bioinformatics methods and tools for HLA haplotype frequency analysis specifically addressing unique issues of haematopoietic stem cell registry data sets. In addition to generating new methods and tools for the analysis of registry data sets, the intent is to produce a comprehensive analysis of HLA data from 20 million donors from the Bone Marrow Donors Worldwide (BMDW) database. This report summarizes the activity on this project as of the 16IHIW meeting in Liverpool.

The impact of frequent HLA haplotypes in high linkage disequilibrium on donor search and clinical outcome after unrelated haematopoietic SCT.

The MHC region on chromosome 6 contains a large number of non-HLA genes next to the HLA genes. Matching for HLA in unrelated hematopoietic SCT (HSCT) does not necessarily mean that these non-HLA genes are also matched. We selected 348 Northwest European patients transplanted with an HLA-A-, -B-, -C-, -DRB1-, -DQB1-matched unrelated donor (MUD) between 1987 and 2008. Patients' haplotypes were identified via descend. We were unable to determine the haplotypes of the donor; therefore we used frequent haplotypes (FH) in high linkage disequilibrium (LD) as a proxy for haplotype matching. Presence of a FH in a patient positively affected the probability and speed of identifying a matched unrelated donor. Competing risk survival analysis showed that patients with one or two FH have a statistically significantly decreased probability of developing ≥ grade II acute GVDH (aGVHD) without increased risk of relapse compared to patients without FH (HR (95% CI): 0.53 (0.31-0.91)). This association was strongest for those FH with the highest LD between both HLA-A and -C or -B, and HLA-C or -B and -DRB1 (HR (95% CI): 0.49 (0.26-0.92)). These results extend evidence that non-HLA allele coding regions have a significant impact on development of ≥ grade II aGVHD. We conclude that there is more to successful HSCT than matching for HLA genes.

Analysis of HLA class I-II haplotype frequency and segregation in a cohort of patients with advanced stage ovarian cancer.

In solid tumors, human leucocyte antigen (HLA)-A2 has been suggested to be a risk factor and a negative prognostic factor. The HLA-A2 allele in Scandinavia has a high prevalence; it decreases with latitude and also with ovarian cancer mortality in Europe. Furthermore, an association of the HLA-A2 allele with severe prognosis in serous adenocarcinoma of the ovary in stages III-IV was found. Thirty-two unrelated Swedish women with relapsing or progressive ovarian cancer were analysed for the genotypes at the HLA-A, HLA-B, HLA-Cw, and HLA-DRB1 loci by the polymerase chain reaction/sequence-specific primer method. The frequencies of HLA alleles of healthy Swedish bone marrow donors provided by the coordinating centre of the Bone Marrow Donors Worldwide Registries, Leiden, the Netherlands were used as controls. When this cohort of epithelial ovarian cancer patients was compared with healthy Swedish donors, the frequency of HLA-A1 and HLA-A2 gene/phenotype appears, although not statistically significant, to be increased in patients with ovarian carcinoma, while HLA-A3 was decreased. HLA-A2 homozygotes were twofold higher in patients. The A2-B8 haplotype was significantly increased (corrected P value). A2-B5, A2-B15, A2-DRB1*03, A2-DRB1*04, A2-B15-Cw3, and A2-B8-DRB1*03 had odds ratio as well as the level of the lower confidence interval above 1 and significant P value only when considered as single, non-corrected analysis. HLA-B15 and HLA-Cw3 were only present in HLA-A2-positive patients showing that the HLA-A2-HLA-Cw3 and HLA-B15 haplotypes were segregated. In this selected cohort with advanced disease, there are indications of an unusual overrepresentation of HLA class I and II genes/haplotypes as well as segregation for the HLA-A2-HLA-Cw3 and HLA-B15 haplotypes. These findings are presented as a descriptive analysis and need further investigations on a larger series of ovarian cancer patients to establish prognostic associations.

Bone Marrow Donors Worldwide and cord blood stem cell transplantation.

Stem cell transplantation using cord blood donors mismatched for one or two HLA antigens can be successful. A simulation study was performed to assess the percentage of patients for whom a matched and/or a one or two antigen mismatched donor is available. Using either bone marrow donors worldwide (BMDW) (n = 925109) or the cord blood donor file in BMDW (n = 11049). As a donor source, for 97% of the Dutch patients a 0 or 1 antigen mismatched cord blood donor was available, for patients from the minority groups this percentage was 55%. These percentages compare favourably with the fraction of 0 mismatched donors found in BMDW.

Validation of large data sets, an essential prerequisite for data analysis: an analytical survey of the Bone Marrow Donors Worldwide.

Large data sets like the Bone Marrow Donors Worldwide (BMDW) data set can be used for population genetic analyses. The qualities of such data sets are unique. To be able to use the BMDW data for analyses, several problems, like limited size and selective DR typing, of the data have to be solved and the quality of the registry data subsets has to be examined. We describe these problems and methods to overcome them. Also, we give an overview of the qualities of the different registry subsets. Sixteen of the twenty-nine examined subsets contain data that can be used for population genetic analysis. We will deal with these analyses in the future. Additionally, we present a method to calculate the minimum number of individuals required for reliable haplotype frequency estimation.