Search

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.

Emm is a high incidence red cell antigen with eight previously reported Emm− probands. Anti-Emm appears to be naturally occurring yet responsible for a clinically significant acute hemolytic transfusion reaction. Previous work suggests that Emm is located on a GPI-anchored protein, but the antigenic epitope and genetic basis have been elusive. We investigated samples from a South Asian Indian family with two Emm− brothers by whole genome sequencing (WGS). Additionally, samples from four unrelated Emm− individuals were investigated for variants in the candidate gene. Filtering for homozygous variants found in the Emm− brothers and by gnomAD frequency of PIGG [NM_001127178.3:c.2624_2625delTA, p.(Leu875*), rs771819481]. PIGG encodes for a transferase, GPI-ethanolaminephosphate transferase II, which adds ethanolamine phosphate (EtNP) to the second mannose in a GPI-anchor. The four additional unrelated Emm− individuals had various PIGG mutations; deletion of Exons 2–3, deletion of Exons 7–9, insertion/deletion (indel) in Exon 3, and new stop codon in Exon 5. The Emm− phenotype is associated with a rare deficiency of PIGG, potentially defining a new Emm blood group system composed of EtNP bound to mannose, part of the GPI-anchor. The results are consistent with the known PI-linked association of the Emm antigen, and may explain the production of the antibody in the absence of RBC transfusion. Any association with neurologic phenotypes requires further research.

BACKGROUND: Although P1 and Xg(a) are known to be associated with the A4GALT and XG genes respectively, the genetic basis of antigen expression has been elusive. Recent reports link both P1 and Xg(a) expression with nucleotide changes in the promotor regions and with antigen negative phenotypes due to disruption of transcription factor binding. STUDY DESIGN AND METHODS: Whole genome sequencing was performed on 113 individuals as part of the MedSeq Project with serologic RBC antigen typing for P1 (n=77) and Xg(a) (n=15). Genomic data were analyzed by two approaches, nucleotide frequency correlation or serologic correlation, to find A4GALT and XG changes associated with P1 and Xg(a) expression. RESULTS: For P1, the frequency approach identified 29 possible associated nucleotide changes, and the serologic approach revealed four among them correlating with the P1+/P1− phenotype: chr22:43,115,523_43,115,520AAAG/delAAAG (rs66781836); chr 22:43,114,551C/T (rs8138197); chr22:43,114,020T/G (rs2143918); and chr22:43,113,793G/T (rs5751348). For Xg(a), the frequency approach identified 82 possible associated nucleotide changes, and among these the serologic approach revealed one correlating with the Xg(a+)/Xg(a–) phenotype: chrX:2,666,384G/C (rs311103). CONCLUSION: A bioinformatics analysis pipeline was created to identify genetic changes responsible for RBC antigen expression. This study, in progress prior to the recently published reports, independently confirms the basis for P1 and Xg(a). Although, this enabled molecular typing of these antigens, the Y chromosome PAR1 region interfered with Xg(a) typing in males. This approach could be used to identify and confirm the genetic basis of antigens, potentially replacing the historical approach using family pedigrees as genomic sequencing becomes common place.

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 data-sets of peptides with measured affinities. We thus demonstrate a strategy for systematically learning the rules of endogenous antigen presentation.

BACKGROUND There are 346 serologically defined red blood cell (RBC) antigens and 33 serologically defined platelet (PLT) antigens, most of which have known genetic changes in 45 RBC or six PLT genes that correlate with antigen expression. Polymorphic sites associated with antigen expression in the primary literature and reference databases are annotated according to nucleotide positions in cDNA. This makes antigen prediction from next-generation sequencing data challenging, since it uses genomic coordinates. STUDY DESIGN AND METHODS The conventional cDNA reference sequences for all known RBC and PLT genes that correlate with antigen expression were aligned to the human reference genome. The alignments allowed conversion of conventional cDNA nucleotide positions to the corresponding genomic coordinates. RBC and PLT antigen prediction was then performed using the human reference genome and whole genome sequencing (WGS) data with serologic confirmation. RESULTS Some major differences and alignment issues were found when attempting to convert the conventional cDNA to human reference genome sequences for the following genes: ABO, A4GALT, RHD, RHCE, FUT3, ACKR1 (previously DARC), ACHE, FUT2, CR1, GCNT2, and RHAG. However, it was possible to create usable alignments, which facilitated the prediction of all RBC and PLT antigens with a known molecular basis from WGS data. Traditional serologic typing for 18 RBC antigens were in agreement with the WGS-based antigen predictions, providing proof of principle for this approach. CONCLUSION Detailed mapping of conventional cDNA annotated RBC and PLT alleles can enable accurate prediction of RBC and PLT antigens from whole genomic sequencing data.

Detection of somatic mutations in human leukocyte antigen (HLA) genes using whole-exome sequencing (WES) is hampered by the high polymorphism of the HLA loci, which prevents alignment of sequencing reads to the human reference genome. We describe a computational pipeline that enables accurate inference of germline alleles of class I HLA-A, B and C genes and subsequent detection of mutations in these genes using the inferred alleles as a reference. Analysis of WES data from 7,930 pairs of tumor and healthy tissue from the same patient revealed 298 nonsilent HLA mutations in tumors from 266 patients. These 298 mutations are enriched for likely functional mutations, including putative loss-of-function events. Recurrence of mutations suggested that these 'hotspot' sites were positively selected. Cancers with recurrent somatic HLA mutations were associated with upregulation of signatures of cytolytic activity characteristic of tumor infiltration by effector lymphocytes, supporting immune evasion by altered HLA function as a contributory mechanism in cancer.

Comprehensive multiple sequence alignments of the multisubunit DNA-dependent RNA polymerase (RNAP) large subunits, including the bacterial β and β′ subunits and their homologs from archaebacterial RNAPs, eukaryotic RNAPs I–III, nuclear–cytoplasmic large double-stranded DNA virus RNAPs, and plant plastid RNAPs, were created [Lane, W. J. and Darst, S. A. (2009). Molecular evolution of multisubunit RNA polymerases: sequence analysis. In press]. The alignments were used to delineate sequence regions shared among all classes of multisubunit RNAPs, defining common, fundamental RNAP features as well as identifying highly conserved positions. Here, we present a systematic, detailed structural analysis of these shared regions and highly conserved positions in terms of the RNAP structure, as well as the RNAP structure/function relationship, when known.

In the Rh blood group system, the RHD gene is bordered by two homologous DNA sequences called the upstream and downstream Rhesus boxes. The most common cause of the D− phenotype in people of European descent is a deletion of the RHD gene region, which results in a hybrid combination of the two Rhesus boxes. PCRbased testing can detect the presence or absence of the hybrid box to determine RHD zygosity. PCR hybrid box testing on fathers can stratify risk for haemolytic disease of the fetus and new born in mothers with anti-D antibodies. Red blood cells and genomic DNA were isolated from 37 individuals of European descent undergoing whole genome sequencing as part of the MedSeq Project. A whole genome sequence-based RHD sequence read depth analysis was used to determine RHD zygosity (homozygous, hemizygous, or null states) with 100% agreement (n=37) when compared to conventional RhD serology and PCR-based hybrid box assay.

Aim Accurate classification of serologically-distinct antigens (Ag) is critical in organ allocation, and nowadays it is possible to define anti-HLA antibody specificity to the allelic level. We examined class I antibodies (Ab) for allelic specificities that may be serologically distinct and evaluated its impact on crossmatch (XM) results and organ allocation. Methods 6726 sera were tested for anti-HLA Ab using LABScreen Single Ag assay (One Lambda), in which 17 class I HLA Ag are represented by more than one allele. Discordance in MFI values between the two allelic beads was examined as a categorical variable to mimic clinical practice, and as a continuous variable (MFI value) to better reflect biological plausibility. To evaluate the significance of an allele-specific Ab on XM results, we examined an unbiased cohort of deceased donor-candidate testing (125,828 CDC XMs between 2014 and 2017). In our lab, all candidates within a blood group have a CDC XM performed, regardless of whether they have donor-specific antibody (DSA). This uniquely allows us to examine whether a candidate with an allele-specific DSA, who has already been excluded from the match run, would indeed have a positive XM against a donor who has a different allele. Statistical analyses were performed using JMP12Pro (SAS Institute). Results 1. The incidence of discordance between the two allelic beads varied depending on whether MFI values were considered as continuous vs categorical variables. 2. The presence of allele-specific Ab for rare alleles was out of proportion to their frequencies in the population. 3. The presence of Ab against only the rare allele was a poor predictor of a positive CDC XM, with a positive predictive value of 0–7% compared with 52.5% if the candidate had DSA against an A or B locus Ag. Conclusions 9 out of 17 Class I HLA antigens had potential allele-specific reactivity, meaning that candidates with antibody selectivity against rare alleles may be unnecessarily excluded from many organ offers based on definition of unacceptable antigens by serological equivalence. Conversely, for centers relying on virtual-XM, prematurely defining alleles as new serological splits could unfortunately allow transplants for which the recipient has DSA. Further testing is needed to rigorously confirm these allele-specific reactivity patterns.

Aim ABO typing of recipients and donors is important for optimal transplant outcomes. ABO typing by serology requires a fresh blood sample, whereas living donor evaluations as well as collection of samples for the donor registry often involve a buccal swab for DNA isolation. ASHI and CAP have recently released guidelines that allow DNA-based genetic ABO determination as a screening tool for initial donor evaluation or for listing donors on the registry. We sought to develop a targeted ABO next generation sequencing (NGS) assay, along with companion interpretive software, for automated data processing and ABO prediction. Methods To develop a target approach, a series of PCR primer combinations spanning the ABO gene were evaluated on blood and buccal swab isolated DNA. The DNA products were run on agarose gel and Agilent bioanalyzer to determine product size and quality. PCR products underwent Illumina TruSight NGS library preparation and sequencing using Illumina MiSeq. A curated ABO allele table was created, and customize software developed to detect and phase ABO alleles from NGS data. Results PCR primer evaluation showed reproducible PCR products spanning exons 6–7 of the ABO gene, which yielded good quality NGS data. These exons encode the major ABO specificities and variations. To verify the ABO NGS assay, a combination of 50 blood and buccal swab samples were tested and found to be 100% concordant with serologic ABO typing which included 19 Group A, 18 Group O, 12 Group B, and 1 Group AB. Conclusions A targeted ABO NGS based assay was developed and optimized for blood and buccal swab isolated DNA. Companion interpretive software, which automates NGS ABO analysis, was developed and correctly determined the ABO type when compared with serology. This assay and software interpretation provides a screening tool for prediction of ABO blood type when no blood sample is available.

Personalized neoantigen therapies for cancer require accurate epitope selection. However, most Class I prediction algorithms in common use today are based on biochemical binding assays that are difficult to scale and do not address processing steps upstream of peptide-MHC binding. Here we present an alternative approach based on the LC-MS/MS identification of MHC Class I-bound peptides. While MS-based profiling is not new, we optimized the system for rule learning by focusing on cell lines expressing only a single HLA-A or HLA-B allele and by collecting parallel transcriptomic and proteomic measurements. Identifying over 24,000 peptides across 16 individual alleles, we were able to discover novel binding motifs, which were validated biochemically, and develop novel neural network algorithms. Furthermore, we systematically interrogated processing rules - discovering a novel motif conserved across multiple cell types - and developed a principled framework for integrating epitope cleavability, expression, and MHC binding potential into an overall ranking. Validating on external datasets, we saw a doubling in positive predictive value with respect to standard approaches. We thus demonstrate a scalable strategy for systematically learning the rules of endogenous antigen presentation that can be deployed for the optimal selection of patient-specific cancer neoantigens. Citation Format: Michael S. Rooney, Jenn Abelin, Siranush Sarkizova, Derin Keskin, Christine Hartigan, Wandi Zhang, John Sidney, William Lane, Jonathan Stevens, Guang L. Zhang, Karl Clauser, Nir Hacohen, Steve Carr, Cathy Wu. Next-generation epitope prediction using mass spectrometry and integrative genomics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-179. doi:10.1158/1538-7445.AM2017-LB-179

Background The MedSeq Project is a randomized clinical trial developing approaches to assess the impact of integrating genome sequencing into clinical medicine. To facilitate the return of results of potential medical relevance to physicians and patients participating in the MedSeq Project, we sought to develop a reporting approach for the effective communication of such findings. Methods Genome sequencing was performed on the Illumina HiSeq platform. Variants were filtered, interpreted, and validated according to methods developed by the Laboratory for Molecular Medicine and consistent with current professional guidelines. The GeneInsight software suite, which is integrated with the Partners HealthCare electronic health record, was used for variant curation, report drafting, and delivery. Results We developed a concise 5–6 page Genome Report (GR) featuring a single-page summary of results of potential medical relevance with additional pages containing structured variant, gene, and disease information along with supporting evidence for reported variants and brief descriptions of associated diseases and clinical implications. The GR is formatted to provide a succinct summary of genomic findings, enabling physicians to take appropriate steps for disease diagnosis, prevention, and management in their patients. Conclusions Our experience highlights important considerations for the reporting of results of potential medical relevance and provides a framework for interpretation and reporting practices in clinical genome sequencing. Electronic supplementary material The online version of this article (doi:10.1186/s12881-014-0134-1) contains supplementary material, which is available to authorized users.

Objective The mechanism for anagrelide's potent platelet lowering activity in human subjects is not well defined. Studies related to anagrelide function have been hampered by its lack of activity in nonhuman primates and water insolubility. In an effort to define the mechanism whereby anagrelide exerts its therapeutic effect, we identified a water-soluble metabolite (anagrelide.met). The availability of anagrelide.met allowed, for the first time, parallel in vitro and in vivo animal studies centered on the mechanisms by which anagrelide lowers platelet levels. Materials and Methods The effects of anagrelide.met on proliferation and maturation of mega-karyocytes (MKs) as well as platelet production were studied both in vitro and in vivo. Results Anagrelide.met is capable of blocking in vitro MK migration by 20% to 40%. At 100 ng/mL, anagrelide.met selectively blocked in vitro MK maturation, resulting in a 50% decrease in the total number of CD41a + MKs, corresponding with a 30% decrease in MK ploidy by day 10 and a 60% decrease by day 20. Daily intraperitoneal injections of anagrelide.met 100 μg into BALB/c mice was sufficient to significantly decrease platelet counts within 24 to 48 hours, stabilizing to 40 to 50% of normal levels by day 5. This was associated with a 45% decrease in the number of developing MKs and an increase in thrombopoietin levels. Anagrelide.met did not alter WBC counts, hematocrit, or bleeding time, or lead to any apparent signs of toxicity. Furthermore, unlike the parent anagrelide compound, anagrelide.met did not inhibit ADP-induced platelet aggregation even at high concentrations (10 μg/mL). Conclusion We describe a cross-species reactive anagrelide metabolite that selectively inhibits MK maturation and migration, lowering platelet levels without influencing platelet aggre-gation.

Cancer mutations yield neo-antigens, which are instrumental to immune-mediated recognition and control of cancer. Vaccine-based therapies targeting neo-antigens will require accurate prediction of which mutations yield peptides presented on polymorphic HLA class I. While in vitro methods have produced increasingly accurate predictors of peptide:MHC binding, there remains a need to define rules for endogenous antigen presentation. Here, we use rapid, high-resolution liquid chromatography mass spectrometry (LC-MS/MS) to identify >24,000 peptides associated with 16 HLA alleles in B cell lines that each express a single HLA allele. The elution of peptides from single HLA alleles allowed us to develop improved rules for endogenous peptide presentation based on the physicochemical properties of binding peptides, patterns of peptide cleavage and abundance of cognate transcripts. Finally, we trained models that integrated MS-derived peptide data and gene expression and demonstrate improved prediction of endogenous peptide presentation in independent datasets. Our strategy thus improves the performance of current predictive algorithms and provides a rapid and scalable method to generate rules for the massive and diverse set of human HLA alleles. Citation Format: Michael S. Rooney, Jennifer G. Abelin, Derin B. Keskin, Siranush Sarkizova, Christina Hartigan, Wandi Zhang, John Sidney, Jonathan Stevens, William J. Lane, Guang L. Zhang, Karl R. Clauser, Nir Hacohen, Steven A. Carr, Catherine J. Wu. High-throughput profiling of HLA allele-specific peptides by MS for improved epitope prediction [abstract]. In: Proceedings of the Second CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; 2016 Sept 25-28; New York, NY. Philadelphia (PA): AACR; Cancer Immunol Res 2016;4(11 Suppl):Abstract nr B089.

A combination of structural approaches yields a complete atomic model of the highly biochemically characterized Escherichia coli RNA polymerase, enabling fuller exploitation of E. coli as a model for understanding transcription., The Escherichia coli transcription system is the best characterized from a biochemical and genetic point of view and has served as a model system. Nevertheless, a molecular understanding of the details of E. coli transcription and its regulation, and therefore its full exploitation as a model system, has been hampered by the absence of high-resolution structural information on E. coli RNA polymerase (RNAP). We use a combination of approaches, including high-resolution X-ray crystallography, ab initio structural prediction, homology modeling, and single-particle cryo-electron microscopy, to generate complete atomic models of E. coli core RNAP and an E. coli RNAP ternary elongation complex. The detailed and comprehensive structural descriptions can be used to help interpret previous biochemical and genetic data in a new light and provide a structural framework for designing experiments to understand the function of the E. coli lineage-specific insertions and their role in the E. coli transcription program., Author Summary Transcription, or the synthesis of RNA from DNA, is one of the most important processes in the cell. The central enzyme of transcription is the DNA-dependent RNA polymerase (RNAP), a large, macromolecular assembly consisting of at least five subunits. Historically, much of our fundamental information on the process of transcription has come from genetic and biochemical studies of RNAP from the model bacterium Escherichia coli. More recently, major breakthroughs in our understanding of the mechanism of action of RNAP have come from high resolution crystal structures of various bacterial, archaebacterial, and eukaryotic enzymes. However, all of our high-resolution bacterial RNAP structures are of enzymes from the thermophiles Thermus aquaticus or T. thermophilus, organisms with poorly characterized transcription systems. It has thus far proven impossible to obtain a high-resolution structure of E. coli RNAP, which has made it difficult to relate the large collection of genetic and biochemical data on RNAP function directly to the available structural information. Here, we used a combination of approaches—high-resolution X-ray crystallography of E. coli RNAP fragments, ab initio structure prediction, homology modeling, and single-particle cryo-electron microscopy—to generate complete atomic models of E. coli RNAP. Our detailed and comprehensive structural models provide the heretofore missing structural framework for understanding the function of the highly characterized E. coli RNAP.

Transcription in all cellular organisms is performed by multi-subunit, DNA-dependent RNA polymerases that synthesize RNA from DNA templates. Previous sequence and structural studies have elucidated the importance of shared regions common to all multi-subunit RNA polymerases. In addition RNA polymerases contain multiple lineage-specific domain insertions involved in protein-protein and protein-nucleic acid interactions. We have created comprehensive multiple sequence alignments using all available sequence data for the multi-subunit RNA polymerase large subunits, including the bacterial β and β′ subunits and their homologues from archaebacterial RNA polymerases, the eukaryotic RNA polymerases I, II, and III, the nuclear-cytoplasmic large double-stranded DNA Virus RNA polymerases, and plant plastid RNA polymerases. In order to overcome technical difficulties inherent to the large subunit sequences, including large sequence length, small and large lineage-specific insertions, split subunits, and fused proteins, we created an automated and customizable sequence retrieval and processing system. In addition, we used our alignments to create a more expansive set of shared sequence regions and bacterial lineage-specific domain insertions. We also analyzed the intergenic gap between the bacterial β and β′ genes.

The control of bacterial transcription initiation depends on a primary sigma factor for housekeeping functions, as well as alternative sigma factors that control regulons in response to environmental stresses. The largest and most diverse subgroup of alternative sigma factors, the group IV extracytoplasmic function sigma factors, directs the transcription of genes that regulate a wide variety of responses, including envelope stress and pathogenesis. We determined the 2.3-A resolution crystal structure of the -35 element recognition domain of a group IV sigma factor, Escherichia coli sigma(E)4, bound to its consensus -35 element, GGAACTT. Despite similar function and secondary structure, the primary and group IV sigma factors recognize their -35 elements using distinct mechanisms. Conserved sequence elements of the sigma(E) -35 element induce a DNA geometry characteristic of AA/TT-tract DNA, including a rigid, straight double-helical axis and a narrow minor groove. For this reason, the highly conserved AA in the middle of the GGAACTT motif is essential for -35 element recognition by sigma(E)4, despite the absence of direct protein-DNA interactions with these DNA bases. These principles of sigma(E)4/-35 element recognition can be applied to a wide range of other group IV sigma factors.