Mass Spectrometry Profiling of HLA-Associated Peptidomes in Mono-allelic Cells Enables More Accurate Epitope Prediction.

Identification of human leukocyte antigen (HLA)-bound peptides by liquid chromatography-tandem mass spectrometry (LC-MS/MS) is poised to provide a deep understanding of rules underlying antigen presentation. However, a key obstacle is the ambiguity that arises from the co-expression of multiple HLA alleles. Here, we have implemented a scalable mono-allelic strategy for profiling the HLA peptidome. By using cell lines expressing a single HLA allele, optimizing immunopurifications, and developing an application-specific spectral search algorithm, we identified thousands of peptides bound to 16 different HLA class I alleles. These data enabled the discovery of subdominant binding motifs and an integrative analysis quantifying the contribution of factors critical to epitope presentation, such as protein cleavage and gene expression. We trained neural-network prediction algorithms with our large dataset (>24,000 peptides) and outperformed algorithms trained on datasets of peptides with measured affinities. We thus demonstrate a strategy for systematically learning the rules of endogenous antigen presentation.

Centriolar kinesin Kif24 interacts with CP110 to remodel microtubules and regulate ciliogenesis.

We have identified a protein, Kif24, that shares homology with the kinesin-13 subfamily of motor proteins and specifically interacts with CP110 and Cep97, centrosomal proteins that play a role in regulating centriolar length and ciliogenesis. Kif24 preferentially localizes to mother centrioles. Loss of Kif24 from cycling cells resulted in aberrant cilia assembly but did not promote growth of abnormally long centrioles, unlike CP110 and Cep97 depletion. We found that loss of Kif24 leads to the disappearance of CP110 from mother centrioles, specifically in cycling cells able to form cilia. Kif24 is able to bind and depolymerize microtubules in vitro. Remarkably, ectopically expressed Kif24 specifically remodels centriolar microtubules without significantly altering cytoplasmic microtubules. Thus, our studies have identified a centriolar kinesin that specifically remodels a subset of microtubules, thereby regulating cilia assembly. These studies also suggest mechanistic differences between the regulation of microtubule elongation associated with centrioles and cilia.

Mapping the NPHP-JBTS-MKS protein network reveals ciliopathy disease genes and pathways.

Nephronophthisis (NPHP), Joubert (JBTS), and Meckel-Gruber (MKS) syndromes are autosomal-recessive ciliopathies presenting with cystic kidneys, retinal degeneration, and cerebellar/neural tube malformation. Whether defects in kidney, retinal, or neural disease primarily involve ciliary, Hedgehog, or cell polarity pathways remains unclear. Using high-confidence proteomics, we identified 850 interactors copurifying with nine NPHP/JBTS/MKS proteins and discovered three connected modules: "NPHP1-4-8" functioning at the apical surface, "NPHP5-6" at centrosomes, and "MKS" linked to Hedgehog signaling. Assays for ciliogenesis and epithelial morphogenesis in 3D renal cultures link renal cystic disease to apical organization defects, whereas ciliary and Hedgehog pathway defects lead to retinal or neural deficits. Using 38 interactors as candidates, linkage and sequencing analysis of 250 patients identified ATXN10 and TCTN2 as new NPHP-JBTS genes, and our Tctn2 mouse knockout shows neural tube and Hedgehog signaling defects. Our study further illustrates the power of linking proteomic networks and human genetics to uncover critical disease pathways.

Galectin-4 Antimicrobial Activity Primarily Occurs Through its C-Terminal Domain.

Although immune tolerance evolved to reduce reactivity with self, it creates a gap in the adaptive immune response against microbes that decorate themselves in self-like antigens. This is particularly apparent with carbohydrate-based blood group antigens, wherein microbes can envelope themselves in blood group structures similar to human cells. In this study, we demonstrate that the innate immune lectin, galectin-4 (Gal-4), exhibits strain-specific binding and killing behavior towards microbes that display blood group-like antigens. Examination of binding preferences using a combination of microarrays populated with ABO(H) glycans and a variety of microbial strains, including those that express blood group-like antigens, demonstrated that Gal-4 binds mammalian and microbial antigens that have features of blood group and mammalian-like structures. Although Gal-4 was thought to exist as a monomer that achieves functional bivalency through its two linked carbohydrate recognition domains, our data demonstrate that Gal-4 forms dimers and that differences in the intrinsic ability of each domain to dimerize likely influences binding affinity. While each Gal-4 domain exhibited blood group-binding activity, the C-terminal domain (Gal-4C) exhibited dimeric properties, while the N-terminal domain (Gal-4N) failed to similarly display dimeric activity. Gal-4C not only exhibited the ability to dimerize but also possessed higher affinity toward ABO(H) blood group antigens and microbes expressing glycans with blood group-like features. Furthermore, when compared to Gal-4N, Gal-4C exhibited more potent antimicrobial activity. Even in the context of the full-length protein, where Gal-4N is functionally bivalent by virtue of Gal-4C dimerization, Gal-4C continued to display higher antimicrobial activity. These results demonstrate that Gal-4 exists as a dimer and exhibits its antimicrobial activity primarily through its C-terminal domain. In doing so, these data provide important insight into key features of Gal-4 responsible for its innate immune activity against molecular mimicry.

Neoantigen vaccine generates intratumoral T cell responses in phase Ib glioblastoma trial.

Neoantigens, which are derived from tumour-specific protein-coding mutations, are exempt from central tolerance, can generate robust immune responses1,2 and can function as bona fide antigens that facilitate tumour rejection3. Here we demonstrate that a strategy that uses multi-epitope, personalized neoantigen vaccination, which has previously been tested in patients with high-risk melanoma4-6, is feasible for tumours such as glioblastoma, which typically have a relatively low mutation load1,7 and an immunologically 'cold' tumour microenvironment8. We used personalized neoantigen-targeting vaccines to immunize patients newly diagnosed with glioblastoma following surgical resection and conventional radiotherapy in a phase I/Ib study. Patients who did not receive dexamethasone-a highly potent corticosteroid that is frequently prescribed to treat cerebral oedema in patients with glioblastoma-generated circulating polyfunctional neoantigen-specific CD4+ and CD8+ T cell responses that were enriched in a memory phenotype and showed an increase in the number of tumour-infiltrating T cells. Using single-cell T cell receptor analysis, we provide evidence that neoantigen-specific T cells from the peripheral blood can migrate into an intracranial glioblastoma tumour. Neoantigen-targeting vaccines thus have the potential to favourably alter the immune milieu of glioblastoma.

Corrigendum: An immunogenic personal neoantigen vaccine for patients with melanoma.

This corrects the article DOI: 10.1038/nature22991.

Impact of allele-specific anti-human leukocyte antigen class I antibodies on organ allocation.

Technological advances in the field of histocompatibility have allowed us to define anti-human leukocyte antigen (HLA) antibody specificity at the allelic level. However, how allele-specific antibodies affect organ allocation is poorly studied. We examined allelic specificities of class I HLA antibodies in 6726 consecutive serum samples from 2953 transplant candidates and evaluated their impact on the corresponding crossmatch and organ allocation. Out of 17 class I HLA antigens represented by >1 allele in the LABScreen single antigen bead assay, 12 had potential allele-specific reactivity. Taking advantage of our unbiased cohort of deceased donor-candidate testing (123,135 complement-dependent cytotoxicity crossmatches between 2014 and 2017), we estimated that the presence of allele-specific antibody detected using a single antigen bead assay (median fluorescence intensity, >3000) against only the rare allele was a poor predictor of a positive complement-dependent cytotoxicity crossmatch, with a positive predictive value of 0% to 7%, compared with 52.5% in allele-concordant class I HLA antibodies against A or B locus antigens. Further, we confirmed allele-specific reactivity using flow crossmatch in 3 scenarios: A11:01/A11:02, A68:01/A68:02, and B44:02/B44:03. Our results suggest that allele-specific antibodies may unnecessarily exclude transplant candidates (up to 10%) from organ offers by overcalling unacceptable antigens; incorporation of selective reactivity pattern in allocation may promote precision matching and more equitable allocation.

ABO Genotyping finds more A2 to B kidney transplant opportunities than lectin-based subtyping.

ABO compatibility is important for kidney transplantation, with longer waitlist times for blood group B kidney transplant candidates. However, kidneys from non-A1 (eg, A2) subtype donors, which express less A antigen, can be safely transplanted into group B recipients. ABO subtyping is routinely performed using anti-A1 lectin, but DNA-based genotyping is also possible. Here, we compare lectin and genotyping testing. Lectin and genotype subtyping was performed on 554 group A deceased donor samples at 2 transplant laboratories. The findings were supported by 2 additional data sets of 210 group A living kidney donors and 124 samples with unclear lectin testing sent to a reference laboratory. In deceased donors, genotyping found 65% more A2 donors than lectin testing, most with weak lectin reactivity, a finding supported in living donors and samples sent for reference testing. DNA sequencing and flow cytometry showed that the discordances were because of several factors, including transfusion, small variability in A antigen levels, and rare ABO∗A2.06 and ABO∗A2.16 sequences. Although lectin testing is the current standard for transplantation subtyping, genotyping is accurate and could increase A2 kidney transplant opportunities for group B candidates, a difference that should reduce group B wait times and improve transplant equity.

Genome screening, reporting, and genetic counseling for healthy populations.

Rapid advancements of genome sequencing (GS) technologies have enhanced our understanding of the relationship between genes and human disease. To incorporate genomic information into the practice of medicine, new processes for the analysis, reporting, and communication of GS data are needed. Blood samples were collected from adults with a PCR-confirmed SARS-CoV-2 (COVID-19) diagnosis (target N = 1500). GS was performed. Data were filtered and analyzed using custom pipelines and gene panels. We developed unique patient-facing materials, including an online intake survey, group counseling presentation, and consultation letters in addition to a comprehensive GS report. The final report includes results generated from GS data: (1) monogenic disease risks; (2) carrier status; (3) pharmacogenomic variants; (4) polygenic risk scores for common conditions; (5) HLA genotype; (6) genetic ancestry; (7) blood group; and, (8) COVID-19 viral lineage. Participants complete pre-test genetic counseling and confirm preferences for secondary findings before receiving results. Counseling and referrals are initiated for clinically significant findings. We developed a genetic counseling, reporting, and return of results framework that integrates GS information across multiple areas of human health, presenting possibilities for the clinical application of comprehensive GS data in healthy individuals.

Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations.

Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.

An immunogenic personal neoantigen vaccine for patients with melanoma.

Effective anti-tumour immunity in humans has been associated with the presence of T cells directed at cancer neoantigens, a class of HLA-bound peptides that arise from tumour-specific mutations. They are highly immunogenic because they are not present in normal tissues and hence bypass central thymic tolerance. Although neoantigens were long-envisioned as optimal targets for an anti-tumour immune response, their systematic discovery and evaluation only became feasible with the recent availability of massively parallel sequencing for detection of all coding mutations within tumours, and of machine learning approaches to reliably predict those mutated peptides with high-affinity binding of autologous human leukocyte antigen (HLA) molecules. We hypothesized that vaccination with neoantigens can both expand pre-existing neoantigen-specific T-cell populations and induce a broader repertoire of new T-cell specificities in cancer patients, tipping the intra-tumoural balance in favour of enhanced tumour control. Here we demonstrate the feasibility, safety, and immunogenicity of a vaccine that targets up to 20 predicted personal tumour neoantigens. Vaccine-induced polyfunctional CD4+ and CD8+ T cells targeted 58 (60%) and 15 (16%) of the 97 unique neoantigens used across patients, respectively. These T cells discriminated mutated from wild-type antigens, and in some cases directly recognized autologous tumour. Of six vaccinated patients, four had no recurrence at 25 months after vaccination, while two with recurrent disease were subsequently treated with anti-PD-1 (anti-programmed cell death-1) therapy and experienced complete tumour regression, with expansion of the repertoire of neoantigen-specific T cells. These data provide a strong rationale for further development of this approach, alone and in combination with checkpoint blockade or other immunotherapies.

USP13 interacts with cohesin and regulates its ubiquitination in human cells.

Cohesin is a multiprotein ring complex that regulates 3D genome organization, sister chromatid cohesion, gene expression, and DNA repair. Cohesin is known to be ubiquitinated, although the mechanism, regulation, and effects of cohesin ubiquitination remain poorly defined. We previously used gene editing to introduce a dual epitope tag into the endogenous allele of each of 11 known components of cohesin in human HCT116 cells. Here we report that mass spectrometry analysis of dual-affinity purifications identified the USP13 deubiquitinase as a novel cohesin-interacting protein. Subsequent immunoprecipitation/Western blots confirmed the endogenous interaction in HCT116, 293T, HeLa, and RPE-hTERT cells; demonstrated that the interaction occurs specifically in the soluble nuclear fraction (not in the chromatin); requires the ubiquitin-binding domains (UBA1/2) of USP13; and occurs preferentially during DNA replication. Reciprocal dual-affinity purification of endogenous USP13 followed by mass spectrometry demonstrated that cohesin is its primary interactor in the nucleus. Ectopic expression and CRISPR knockout of USP13 showed that USP13 is paradoxically required for both deubiquitination and ubiquitination of cohesin subunits in human cells. USP13 was dispensable for sister chromatid cohesion in HCT116 and HeLa cells, whereas it was required for the dissociation of cohesin from chromatin as cells transit through mitosis. Together these results identify USP13 as a new cohesin-interacting protein that regulates the ubiquitination of cohesin and its cell cycle regulated interaction with chromatin.

Recommendation for validation and quality assurance of non-invasive prenatal testing for foetal blood groups and implications for IVD risk classification according to EU regulations.

BACKGROUND AND OBJECTIVES: Non-invasive assays for predicting foetal blood group status in pregnancy serve as valuable clinical tools in the management of pregnancies at risk of detrimental consequences due to blood group antigen incompatibility. To secure clinical applicability, assays for non-invasive prenatal testing of foetal blood groups need to follow strict rules for validation and quality assurance. Here, we present a multi-national position paper with specific recommendations for validation and quality assurance for such assays and discuss their risk classification according to EU regulations. MATERIALS AND METHODS: We reviewed the literature covering validation for in-vitro diagnostic (IVD) assays in general and for non-invasive foetal RHD genotyping in particular. Recommendations were based on the result of discussions between co-authors. RESULTS: In relation to Annex VIII of the In-Vitro-Diagnostic Medical Device Regulation 2017/746 of the European Parliament and the Council, assays for non-invasive prenatal testing of foetal blood groups are risk class D devices. In our opinion, screening for targeted anti-D prophylaxis for non-immunized RhD negative women should be placed under risk class C. To ensure high quality of non-invasive foetal blood group assays within and beyond the European Union, we present specific recommendations for validation and quality assurance in terms of analytical detection limit, range and linearity, precision, robustness, pre-analytics and use of controls in routine testing. With respect to immunized women, different requirements for validation and IVD risk classification are discussed. CONCLUSION: These recommendations should be followed to ensure appropriate assay performance and applicability for clinical use of both commercial and in-house assays.

International Society of Blood Transfusion Working Party on Red Cell Immunogenetics and Blood Group Terminology Report of Basel and three virtual business meetings: Update on blood group systems.

BACKGROUND AND OBJECTIVES: Under the ISBT, the Working Party (WP) for Red Cell Immunogenetics and Blood Group Terminology is charged with ratifying blood group systems, antigens and alleles. This report presents the outcomes from four WP business meetings, one located in Basel in 2019 and three held as virtual meetings during the COVID-19 pandemic in 2020 and 2021. MATERIALS AND METHODS: As in previous meetings, matters pertaining to blood group antigen nomenclature were discussed. New blood group systems and antigens were approved and named according to the serologic, genetic, biochemical and cell biological evidence presented. RESULTS: Seven new blood group systems, KANNO (defined numerically as ISBT 037), SID (038), CTL2 (039), PEL (040), MAM (041), EMM (042) and ABCC1 (043) were ratified. Two (039 and 043) were de novo discoveries, and the remainder comprised reported antigens where the causal genes were previously unknown. A further 15 blood group antigens were added to the existing blood group systems: MNS (002), RH (004), LU (005), DI (010), SC (013), GE (020), KN (022), JMH (026) and RHAG (030). CONCLUSION: The ISBT now recognizes 378 antigens, of which 345 are clustered within 43 blood group systems while 33 still have an unknown genetic basis. The ongoing discovery of new blood group systems and antigens underscores the diverse and complex biology of the red cell membrane. The WP continues to update the blood group antigen tables and the allele nomenclature tables. These can be found on the ISBT website (http://www.isbtweb.org/working-parties/red-cell-immunogenetics-and-blood-group-terminology/).

Utility of Routine Computed Tomography Angiogram Screening for Blunt Carotid Injury in Isolated Mandible Fractures.

PURPOSE: Determine the incidence of blunt carotid injury (BCAI) in the blunt trauma induced isolated mandible fracture population and determine the utility of computed tomographic angiography (CTA) screening for BCAI in this population. METHODS: A retrospective data review was performed on patients presenting to Boston Medical Center from January 2008 to January 2019 with a diagnosis of mandible fracture. Population selected utilizing ICD-9 and ICD-10 diagnosis codes for mandible fracture and BCAI. Excluded populations were pediatric (less than 18 years) and penetrating carotid injuries. RESULTS: A total of 1,508 mandible fractures were identified, with 73% (n = 1,103) being isolated injuries. Five BCAIs were identified, and of these, 1 was associated with an isolated mandible fracture (incidence <0.01%). One of 110 isolated mandible fractures screened for BCAI with the use of CTA was positive (incidence 0.9%). CONCLUSIONS: BCAI is a rare finding in isolated mandible fractures. The inclusion of this population as an independent risk factor for BCAI should be questioned. Routine screening with CTA is not warranted.

Systematic proteomics of endogenous human cohesin reveals an interaction with diverse splicing factors and RNA-binding proteins required for mitotic progression.

The cohesin complex regulates sister chromatid cohesion, chromosome organization, gene expression, and DNA repair. Cohesin is a ring complex composed of four core subunits and seven regulatory subunits. In an effort to comprehensively identify additional cohesin-interacting proteins, we used gene editing to introduce a dual epitope tag into the endogenous allele of each of 11 known components of cohesin in cultured human cells, and we performed MS analyses on dual-affinity purifications. In addition to reciprocally identifying all known components of cohesin, we found that cohesin interacts with a panoply of splicing factors and RNA-binding proteins (RBPs). These included diverse components of the U4/U6.U5 tri-small nuclear ribonucleoprotein complex and several splicing factors that are commonly mutated in cancer. The interaction between cohesin and splicing factors/RBPs was RNA- and DNA-independent, occurred in chromatin, was enhanced during mitosis, and required RAD21. Furthermore, cohesin-interacting splicing factors and RBPs followed the cohesin cycle and prophase pathway of cell cycle-regulated interactions with chromatin. Depletion of cohesin-interacting splicing factors and RBPs resulted in aberrant mitotic progression. These results provide a comprehensive view of the endogenous human cohesin interactome and identify splicing factors and RBPs as functionally significant cohesin-interacting proteins.

A cloning and expression system to probe T-cell receptor specificity and assess functional avidity to neoantigens.

Recent studies have highlighted the promise of targeting tumor neoantigens to generate potent antitumor immune responses and provide strong motivation for improving our understanding of antigen-T-cell receptor (TCR) interactions. Advances in single-cell sequencing technologies have opened the door for detailed investigation of the TCR repertoire, providing paired information from TCRα and TCRß, which together determine specificity. However, a need remains for efficient methods to assess the specificity of discovered TCRs. We developed a streamlined approach for matching TCR sequences with cognate antigen through on-demand cloning and expression of TCRs and screening against candidate antigens. Here, we first demonstrate the system's capacity to identify viral-antigen-specific TCRs and compare the functional avidity of TCRs specific for a given antigen target. We then apply this system to identify neoantigen-specific TCR sequences from patients with melanoma treated with personalized neoantigen vaccines and characterize functional avidity of neoantigen-specific TCRs. Furthermore, we use a neoantigen-prediction pipeline to show that an insertion-deletion mutation in a putative chronic lymphocytic leukemia (CLL) driver gives rise to an immunogenic neoantigen mut-MGA, and use this approach to identify the mut-MGA-specific TCR sequence. This approach provides a means to identify and express TCRs, and then rapidly assess antigen specificity and functional avidity of a reconstructed TCR, which can be applied for monitoring antigen-specific T-cell responses, and potentially for guiding the design of effective T-cell-based immunotherapies.

Multiple GYPB gene deletions associated with the U- phenotype in those of African ancestry.

BACKGROUND: The MNS blood group system is defined by three homologous genes: GYPA, GYPB, and GYPE. GYPB encodes for glycophorin B (GPB) carrying S/s and the "universal" antigen U. RBCs of approximately 1% of individuals of African ancestry are U- due to absence of GPB. The U- phenotype has long been attributed to a deletion encompassing GYPB exons 2 to 5 and GYPE exon 1 (GYPB*01N). STUDY DESIGN AND METHODS: Samples from two U-individuals underwent Illumina short read whole genome sequencing (WGS) and Nanopore long read WGS. In addition, two existing WGS datasets, MedSeq (n = 110) and 1000 Genomes (1000G, n = 2535), were analyzed for GYPB deletions. Deletions were confirmed by Sanger sequencing. Twenty known U- donor samples were tested by a PCR assay to determine the specific deletion alleles present in African Americans. RESULTS: Two large GYPB deletions in U- samples of African ancestry were identified: a 110 kb deletion extending left of GYPB (DEL_B_LEFT) and a 103 kb deletion extending right (DEL_B_RIGHT). DEL_B_LEFT and DEL_B_RIGHT were the most common GYPB deletions in the 1000 Genomes Project 669 African genomes (allele frequencies 0.04 and 0.02). Seven additional deletions involving GYPB were seen in African, Admixed American, and South Asian samples. No samples analyzed had GYPB*01N. CONCLUSIONS: The U- phenotype in those of African ancestry is primarily associated with two different complete deletions of GYPB (with intact GYPE). Seven additional less common GYPB deletion backgrounds were found. GYPB*01N, long assumed to be the allele commonly encoding U- phenotypes, appears to be rare.

Overcoming the challenges of interpreting complex and uncommon RH alleles from whole genomes.

BACKGROUND AND OBJECTIVES: Rh is one of the most diverse and complex blood group systems. Recently, next generation sequencing (NGS) has proven to be a viable option for RH genotyping. We have developed automated software (bloodTyper) for determining alleles encoding RBC antigens from NGS-based whole genome sequencing (WGS). The bloodTyper algorithm has not yet been optimized and evaluated for complex and uncommon RH alleles. MATERIALS AND METHODS: Twenty-two samples with previous polymerase chain reaction (PCR) and Sanger sequencing-based RH genotyping underwent WGS. bloodTyper was used to detect RH alleles including those defined by structural variation (SV) using a combination of three independent strategies: sequence read depth of coverage, split reads and paired reads. RESULTS: bloodTyper was programmed to identify D negative and positive phenotypes as well as the presence of alleles encoding weak D, partial D and variant RHCE. Sequence read depth of coverage calculation accurately determined RHD zygosity and detected the presence of RHD/RHCE hybrids. RHCE*C was determined by sequence read depth of coverage and by split read methods. RHD hybrid alleles and RHCE*C were confirmed by using a paired read approach. Small SVs present in RHCE*CeRN and RHCE*ceHAR were detected by a combined read depth of coverage and paired read approach. CONCLUSIONS: The combination of several different interpretive approaches allowed for automated software based-RH genotyping of WGS data including RHD zygosity and complex compound RHD and RHCE heterozygotes. The scalable nature of this automated analysis will enable RH genotyping in large genomic sequencing projects.

Automated typing of red blood cell and platelet antigens from whole exome sequences.

BACKGROUND: Genotyping has expanded the number red blood cell (RBC) and platelet (PLT) antigens that can readily be typed, but often represents an additional testing cost. The analysis of existing genomic data offers a cost-effective approach. We recently developed automated software (bloodTyper) for determination of RBC and PLT antigens from whole genome sequencing. Here we extend the algorithm to whole exome sequencing (WES). STUDY DESIGN AND METHODS: Whole exome sequencing was performed on samples from 75 individuals. WES-based bloodTyper RBC and PLT typing was compared to conventional polymerase chain reaction (PCR) RHD zygosity testing and serologic and single-nucleotide polymorphism (SNP) typing for 38 RBC antigens in 12 systems (17 serologic and 35 SNPs) and 22 PLT antigens (22 SNPs). Samples from the first 20 individuals were used to modify bloodTyper to interpret WES followed by blinded typing of 55 samples. RESULTS: Over the first 20 samples, discordances were noted for C, M, and N antigens, which were due to WES-specific biases. After modification, bloodTyper was 100% accurate on blinded evaluation of the last 55 samples and outperformed both serologic (99.67% accurate) and SNP typing (99.97% accurate) reflected by two Fyb and one N serologic typing errors and one undetected SNP encoding a Jknull phenotype. RHD zygosity testing by bloodTyper was 100% concordant with a combination of hybrid Rhesus box PCR and PCR-restriction fragment length polymorphism for all samples. CONCLUSION: The automated bloodTyper software was modified for WES biases to allow for accurate RBC and PLT antigen typing. Such analysis could become a routing part of future WES efforts.