High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

HLA-E typing of more than 2.5 million potential hematopoietic stem cell donors: Methods and population-specific allele frequencies

We developed a cost-efficient workflow for genotyping HLA-E by NGS and applied it for genotyping more than 2.5 million potential stem cell donors. The data obtained were used to determine HLA-E allele frequency distributions for 104 populations. Our results confirm the known dominance of the alleles E*01:01 and E*01:03, which have a combined frequency of more than 0.99 in 97 of the 104 populations. E*01:01 is more frequent in Africa and the western part of South America, E*01:03 in Southeast and East Asia. E*01:03 shows a pronounced regional substructure at the high-resolution level with E*01:03:01G being particularly common in a large connected region extending from Turkey to China, E*01:03:02G in Northwestern Europe and E*01:03:03 in Central and Eastern Europe as well as Central Asia. The presented results are relevant both as a basis for further population genetics studies and for optimizing stem cell donor searches.

Haplotype Motif-Based Models for KIR-Genotype Informed Selection of Hematopoietic Cell Donors Fail to Predict Outcome of Patients With Myelodysplastic Syndromes or Secondary Acute Myeloid Leukemia

Results from registry studies suggest that harnessing Natural Killer (NK) cell reactivity mediated through Killer cell Immunoglobulin-like Receptors (KIR) could reduce the risk of relapse after allogeneic Hematopoietic Cell Transplantation (HCT). Several competing models have been developed to classify donors as KIR-advantageous or disadvantageous. Basically, these models differ by grouping donors based on distinct KIR–KIR–ligand combinations or by haplotype motif assignment. This study aimed to validate different models for unrelated donor selection for patients with Myelodysplatic Syndromes (MDS) or secondary Acute Myeloid Leukemia (sAML). In a joint retrospective study of the European Society for Blood and Marrow Transplantation (EBMT) and the Center for International Blood and Marrow Transplant Research (CIBMTR) registry data from 1704 patients with secondary AML or MDS were analysed. The cohort consisted mainly of older patients (median age 61 years) with high risk disease who had received chemotherapy-based reduced intensity conditioning and anti-thymocyte globulin prior to allogeneic HCT from well-matched unrelated stem cell donors. The impact of the predictors on Overall Survival (OS) and relapse incidence was tested in Cox regression models adjusted for patient age, a modified disease risk index, performance status, donor age, HLA-match, sex-match, CMV-match, conditioning intensity, type of T-cell depletion and graft type. KIR genes were typed using high-resolution amplicon-based next generation sequencing. In univariable and multivariable analyses none of the models predicted OS and the risk of relapse consistently. Our results do not support the hypothesis that optimizing NK-mediated alloreactivity is possible by KIR-genotype informed selection of HLA-matched unrelated donors. However, in the context of allogeneic transplantation, NK-cell biology is complex and only partly understood. KIR-genes are highly diverse and current assignment of haplotype motifs based on the presence or absence of selected KIR genes is over-simplistic. As a consequence, further research is highly warranted and should integrate cutting edge knowledge on KIR genetics, and NK-cell biology into future studies focused on homogeneous groups of patients and treatment modalities.

High-resolution HLA allele and haplotype frequencies of the Saudi Arabian population based on 45,457 individuals and corresponding stem cell donor matching probabilities

We estimated HLA allele and haplotype frequencies of the Saudi Arabian population from a sample of 45,457 registered stem cell donors. The most frequent HLA alleles were A*02:01g (18.5%), C*06:02g (16.1%), B*51:01g (14.1%), DRB1*07:01g (16.2%), DQB1*02:01g (30.5%), and DPB1*04:01g (33.6%). The most frequent 5-locus haplotypes were A*02:05g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.73%), A*02:01g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.66%), and A*26:01g~C*07:02g~B*08:01g~DRB1*03:01g~DQB1*02:01g (1.38%). Furthermore, we used the calculated haplotype frequencies to estimate stem cell donor matching probabilities for Saudi Arabian donor and patient populations under various matching requirements. These results are relevant for strategic donor registry planning in the Kingdom of Saudi Arabia.

Estimation of German KIR Allele Group Haplotype Frequencies

The impact of the highly polymorphic Killer-cell immunoglobulin-like receptor (KIR) gene cluster on the outcome of hematopoietic stem cell transplantation (HCST) is subject of current research. To further understand the involvement of this gene family into Natural Killer (NK) cell-mediated graft-versus-leukemia reactions, knowledge of haplotype structures, and allelic linkage is of importance. In this analysis, we estimate population-specific KIR haplotype frequencies at allele group resolution in a cohort of n = 458 German families. We addressed the polymorphism of the KIR gene complex and phasing ambiguities by a combined approach. Haplotype inference within first-degree family relations allowed us to limit the number of possible diplotypes. Structural restriction to a pattern set of 92 previously described KIR copy number haplotypes further reduced ambiguities. KIR haplotype frequency estimation was finally accomplished by means of an expectation-maximization algorithm. Applying a resolution threshold of ½ n, we were able to identify a set of 551 KIR allele group haplotypes, representing 21 KIR copy number haplotypes. The haplotype frequencies allow studying linkage disequilibrium in two-locus as well as in multi-locus analyses. Our study reveals associations between KIR haplotype structures and allele group frequencies, thereby broadening our understanding of the KIR gene complex.

High-Throughput MICA/B Genotyping of Over Two Million Samples: Workflow and Allele Frequencies

MICA and MICB are ligands of the NKG2D receptor and thereby influence NK and T cell activity. MICA/B gene polymorphisms, expression levels and the amount of soluble MICA/B in the serum have been linked to autoimmune diseases, infections, and cancer. In hematopoietic stem cell transplantation, MICA matching between donor and patient has been correlated with reduced acute and chronic graft-vs.-host disease and improved survival. Hence, we developed an extremely cost-efficient high-throughput workflow for genotyping MICA/B for newly registered potential stem cell donors. Since mid-2017, we have genotyped over two million samples using NGS amplicon sequencing for MICA/B exons 2–5. In donors of German origin, MICA*008 is the most common MICA allele with a frequency of 42.3%. It is followed by MICA*002 (11.7%) and MICA*009 (8.8%). The three most common MICB alleles are MICB*005 (43.9%), MICB*004 (21.7%), and MICB*002 (18.9%). In general, MICB is less diverse than MICA and only 6 alleles, instead of 15, account for a cumulative allele frequency of 99.5%. In 0.5% of the samples we observed at least one allele of MICA or MICB which has so far not been reported to the IPD/IMGT-HLA database. By providing MICA/B typed voluntary donors, clinicians become empowered to include MICA/B into their donor selection process to further improve unrelated hematopoietic stem cell transplantation.

The biochemical properties of the Arabidopsis ecto-nucleoside triphosphate diphosphohydrolase AtAPY1 contradict a direct role in purinergic signaling

The Arabidopsis E-NTPDase (ecto-nucleoside triphosphate diphosphohydrolase) AtAPY1 was previously shown to be involved in growth and development, pollen germination and stress responses. It was proposed to perform these functions through regulation of extracellular ATP signals. However, a GFP-tagged version was localized exclusively in the Golgi and did not hydrolyze ATP. In this study, AtAPY1 without the bulky GFP-tag was biochemically characterized with regard to its suggested role in purinergic signaling. Both the full-length protein and a soluble form without the transmembrane domain near the N-terminus were produced in HEK293 cells. Of the twelve nucleotide substrates tested, only three--GDP, IDP and UDP--were hydrolyzed, confirming that ATP was not a substrate of AtAPY1. In addition, the effects of pH, divalent metal ions, known E-NTPDase inhibitors and calmodulin on AtAPY1 activity were analyzed. AtAPY1-GFP extracted from transgenic Arabidopsis seedlings was included in the analyses. All three AtAPY1 versions exhibited very similar biochemical properties. Activity was detectable in a broad pH range, and Ca(2+), Mg(2+) and Mn(2+) were the three most efficient cofactors. Of the inhibitors tested, vanadate was the most potent one. Surprisingly, sulfonamide-based inhibitors shown to inhibit other E-NTPDases and presumed to inhibit AtAPY1 as well were not effective. Calmodulin stimulated the activity of the GFP-tagless membranous and soluble AtAPY1 forms about five-fold, but did not alter their substrate specificities. The apparent Km values obtained with AtAPY1-GFP indicate that AtAPY1 is primarily a GDPase. A putative three-dimensional structural model of the ecto-domain is presented, explaining the potent inhibitory potential of vanadate and predicting the binding mode of GDP. The found substrate specificity classifies AtAPY1 as a nucleoside diphosphatase typical of N-terminally anchored Golgi E-NTPDases and negates a direct function in purinergic signaling.

The Arabidopsis apyrase AtAPY1 is localized in the Golgi instead of the extracellular space

BACKGROUND

The two highly similar Arabidopsis apyrases AtAPY1 and AtAPY2 were previously shown to be involved in plant growth and development, evidently by regulating extracellular ATP signals. The subcellular localization of AtAPY1 was investigated to corroborate an extracellular function.

RESULTS

Transgenic Arabidopsis lines expressing AtAPY1 fused to the SNAP-(O(6)-alkylguanine-DNA alkyltransferase)-tag were used for indirect immunofluorescence and AtAPY1 was detected in punctate structures within the cell. The same signal pattern was found in seedlings stably overexpressing AtAPY1-GFP by indirect immunofluorescence and live imaging. In order to identify the nature of the AtAPY1-positive structures, AtAPY1-GFP expressing seedlings were treated with the endocytic marker stain FM4-64 (N-(3-triethylammoniumpropyl)-4-(p-diethylaminophenyl-hexatrienyl)-pyridinium dibromide) and crossed with a transgenic line expressing the trans-Golgi marker Rab E1d. Neither FM4-64 nor Rab E1d co-localized with AtAPY1. However, live imaging of transgenic Arabidopsis lines expressing AtAPY1-GFP and either the fluorescent protein-tagged Golgi marker Membrin 12, Syntaxin of plants 32 or Golgi transport 1 protein homolog showed co-localization. The Golgi localization was confirmed by immunogold labeling of AtAPY1-GFP. There was no indication of extracellular AtAPY1 by indirect immunofluorescence using antibodies against SNAP and GFP, live imaging of AtAPY1-GFP and immunogold labeling of AtAPY1-GFP. Activity assays with AtAPY1-GFP revealed GDP, UDP and IDP as substrates, but neither ATP nor ADP. To determine if AtAPY1 is a soluble or membrane protein, microsomal membranes were isolated and treated with various solubilizing agents. Only SDS and urea (not alkaline or high salt conditions) were able to release the AtAPY1 protein from microsomal membranes.

CONCLUSIONS

AtAPY1 is an integral Golgi protein with the substrate specificity typical for Golgi apyrases. It is therefore not likely to regulate extracellular nucleotide signals as previously thought. We propose instead that AtAPY1 exerts its growth and developmental effects by possibly regulating glycosylation reactions in the Golgi.