Ultra-high resolution HLA genotyping and allele discovery by highly multiplexed cDNA amplicon pyrosequencing

ABSTRACT: BACKGROUND: High-resolution HLA genotyping is a critical diagnostic and research assay. Current methods rarely achieve unambiguous high-resolution typing without making population-specific frequency inferences due to a lack of locus coverage and difficulty in exon-phase matching. Achieving high-resolution typing is also becoming more challenging with traditional methods as the database of known HLA alleles increases. RESULTS: We designed a cDNA amplicon-based pyrosequencing method to capture 94% of the HLA class I open-reading-frame with only two amplicons per sample, and an analogous method for class II HLA genes, with a primary focus on sequencing the DRB loci. We present a novel Galaxy server-based analysis workflow for determining genotype. During assay validation, we performed two GS Junior sequencing runs to determine the accuracy of the HLA class I amplicons and DRB amplicon at different levels of multiplexing. When 116 amplicons were multiplexed, we unambiguously resolved 99%of class I alleles to four- or six-digit resolution, as well as 100% unambiguous DRB calls. The second experiment, with 271 multiplexed amplicons, missed some alleles, but generated high-resolution, concordant typing for 93% of class I alleles, and 96% for DRB1 alleles. In a third, preliminary experiment we attempted to sequence novel amplicons for other class II loci with mixed success. CONCLUSIONS: The presented assay is higher-throughput and higher-resolution than existing HLA genotyping methods, and suitable for allele discovery or large cohort sampling. The validated class I and DRB primers successfully generated unambiguously high-resolution genotypes, while further work is needed to validate additional class II genotyping amplicons.     

Contemporary genetic structure and postglacial demographic history of the black scorpionfish,Scorpaena porcus,in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.     

Arachidonic acid metabolites and colchicine in Beh?et's disease (BD)

Vasculitis is accepted to be the basis of Beh?et's disease (BD) which is a multisystem disease, and the arachidonic acid(AA) metabolites acting as balancing mediators in the organism are accepted to be responsible for the vasculitis. In this study, we examined the prostaglandin E2 (PGE2) and leukotriene C4 (LTC4) levels of the patients with BD before and after colchicine therapy. We found a statistical decrease in the PGE2 and LTC4 levels after colchicine therapy compared to the previous levels, concluding that colchicine inhibits the inflammation and the polymorphonuclear leukocyte (PML) chemotaxis by inhibiting the cyclooxygenase and lipoxygenase pathways.     

Small and Large Ribosomal Subunit Deficiencies Lead to Distinct Gene Expression Signatures that Reflect Cellular Growth Rate

1. Download high-res image (149KB)
2. Download full-size image

     

Phylogenetic relationships among four species and a sub-species of Mullidae (Actinopterygii; Perciformes) based on mitochondrial cytochrome B, 12S rRNA and cytochrome oxidase II genes

DNA sequence comparisons of three mitochondrial DNA genes were used to reveal phylogenetic relationships among four species and a sub-species of Mullidae family. This is the first report using mitochondrial DNA sequence data to infer intraspecific relationship among different populations of Mullus barbatus and Mullus surmuletus; phylogenetic relationships between M. barbatus and its sub-species; M. barbatus ponticus. Cytochrome b, 12S ribosomal RNA, and cytochrome oxidase II regions of 242 individuals belonging to species M. barbatus, M. surmuletus, Upeneus moluccensis, Upeneus pori and sub-species M. barbatus ponticus were sequenced and phylogenetic trees were constructed using four different algorithms. The phylogenetic trees constructed support the existing taxonomical data of two mullid genera (Mullus, Upeneus). Molecular data shows no significant difference between same species of different geographical populations. The results suggest that the molecular difference is not large enough between M. barbatus and M. barbatus ponticus to consider them as sub-species.     

Analysis and network representation of hotspots in protein interfaces using minimum cut trees

We propose a novel approach to analyze and visualize residue contact networks of protein interfaces by graph‐based algorithms using a minimum cut tree (mincut tree). Edges in the network are weighted according to an energy function derived from knowledge‐based potentials. The mincut tree, which is constructed from the weighted residue network, simplifies and summarizes the complex structure of the contact network by an efficient and informative representation. This representation offers a comprehensible view of critical residues and facilitates the inspection of their organization. We observed, on a nonredundant data set of 38 protein complexes with experimental hotspots that the highest degree node in the mincut tree usually corresponds to an experimental hotspot. Further, hotspots are found in a few paths in the mincut tree. In addition, we examine the organization of hotspots (hot regions) using an iterative clustering algorithm on two different case studies. We find that distinct hot regions are located on specific sites of the mincut tree and some critical residues hold these clusters together. Clustering of the interface residues provides information about the relation of hot regions with each other. Our new approach is useful at the molecular level for both identification of critical paths in the protein interfaces and extraction of hot regions by clustering of the interface residues. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

A Conserved E7-derived Cytotoxic T Lymphocyte Epitope Expressed on Human Papillomavirus 16-transformed HLA-A2+ Epithelial Cancers

Human Papillomavirus 16 (HPV-16) has been identified as the causative agent of 50% of cervical cancers and many other HPV-associated tumors. The transforming potential/tumor maintenance capacity of this high risk HPV is mediated by two viral oncoproteins, E6 and E7, making them attractive targets for therapeutic vaccines. Of 21 E6 and E7 peptides computed to bind HLA-A*0201, 10 were confirmed through TAP-deficient T2 cell HLA stabilization assay. Those scoring positive were investigated to ascertain which were naturally processed and presented by surface HLA molecules for CTL recognition. Because IFNγ ELISpot frequencies from healthy HPV-exposed blood donors against HLA-A*0201-binding peptides were unable to identify specificities for tumor targeting, their physical presence among peptides eluted from HPV-16-transformed epithelial tumor HLA-A*0201 immunoprecipitates was analyzed by MS³ Poisson detection mass spectrometry. Only one epitope (E7_11–19) highly conserved among HPV-16 strains was detected. This 9-mer serves to direct cytolysis by T cell lines, whereas a related 10-mer (E7_11–20), previously used as a vaccine candidate, was neither detected by MS³ on HPV-transformed tumor cells nor effectively recognized by 9-mer specific CTL. These data underscore the importance of precisely defining CTL epitopes on tumor cells and offer a paradigm for T cell-based vaccine design.     

Elicitation from virus-naive individuals of cytotoxic T lymphocytes directed against conserved HIV-1 epitopes

Cytotoxic T lymphocytes (CTL) protect against viruses including HIV-1. To avoid viral escape mutants that thwart immunity, we chose 25 CTL epitopes defined in the context of natural infection with functional and/or structural constraints that maintain sequence conservation. By combining HLA binding predictions with knowledge concerning HLA allele frequencies, a metric estimating population protection coverage (PPC) was computed and epitope pools assembled. Strikingly, only a minority of immunocompetent HIV-1 infected individuals responds to pools with PPC >95%. In contrast, virus-naive individuals uniformly expand IFNγ producing cells and mount anti-HIV-1 cytolytic activity. This disparity suggests a vaccine design paradigm shift from infected to normal subjects.     

Protein-protein interactions: hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: implications for docking

Energetic hot spots account for a significant portion of the total binding free energy and correlate with structurally conserved interface residues. Here, we map experimentally determined hot spots and structurally conserved residues to investigate their geometrical organization. Unfilled pockets are pockets that remain unfilled after protein-protein complexation, while complemented pockets are pockets that disappear upon binding, representing tightly fit regions. We find that structurally conserved residues and energetic hot spots are strongly favored to be located in complemented pockets, and are disfavored in unfilled pockets. For the three available protein-protein complexes with complemented pockets where both members of the complex were alanine-scanned, 62% of all hot spots (DeltaDeltaG>2kcal/mol) are within these pockets, and 60% of the residues in the complemented pockets are hot spots. 93% of all red-hot residues (DeltaDeltaG>/=4kcal/mol) either protrude into or are located in complemented pockets. The occurrence of hot spots and conserved residues in complemented pockets highlights the role of local tight packing in protein associations, and rationalizes their energetic contribution and conservation. Complemented pockets and their corresponding protruding residues emerge among the most important geometric features in protein-protein interactions. By screening the solvent, this organization shields backbone hydrogen bonds and charge-charge interactions. Complemented pockets often pre-exist binding. For 18 protein-protein complexes with complemented pockets whose unbound structures are available, in 16 the pockets are identified to pre-exist in the unbound structures. The root-mean-squared deviations of the atoms lining the pockets between the bound and unbound states is as small as 0.9A, suggesting that such pockets constitute features of the populated native state that may be used in docking.