Molecular definition of an elusive third HLA-A9 molecule: HLA-A9.3

The HLA-A9 family has been characterized as possessing two well defined specificities; HLA-A23 and A24. Serological studies have suggested the presence of a third member of this family HLA-A9.3, however there is doubt surrounding the existence of this specificity. HLA-A23, A24, and the putative A9.3 proteins were analyzed biochemically by immunoprecipitation and isoelectric focusing. Both HLA-A24 and A9.3 have identical isoelectric points whereas A23 is different. We have sequenced cDNA encoding HLA-A23, A24, and A9.3. From the observed protein sequences, we found A9.3 to differ from A24 by two amino acid substitutions located in the ₂ helix of the class I molecule. These substitutions are expected to significantly change the shape of the peptide binding cleft.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M64740 (HLA-A ^*2402); M64741 (HLA-A ^*2403); M64742 (HLA-A ^*2301). Address correspondence and offprint requests to: P. Parham.     

Structure of the HLA-A*0204 antigen,found in South American Indians. Spatial clustering of HLA-A2 subtype polymorphism

The primary structure of the HLA-A2 subtype A^*0204 (isoelectric focusing variant A2.A) has been determined. cDNA encoding this subtype was amplified by the polymerase chain reaction. Four independent full-lenght cDNA clones encoding A^*0204 were analyzed to obtain a consensus sequence for this subtype. A^*0204 differs from A^*0201 by a single nucleotide change of G to T through the coding regions, resulting in an Arg to Met change at position 97. This substitution accounts for the isoelectric focusing pattern of the subtype. The same change occurs in other HLA-A specificities in association with other changes in its vicinity. The absence of additional substitutions in A^*0204 suggests that it could have arisen from A^*0201 by point mutation, and that recurrent mutations may take place during HLA diversification. The spatial location of this change implies that A^*0204 must be a functional variant. Comparison of its sequence with other HLA-A2 subtypes reveals that much of the HLA-A2 subtype polymorphism is generated by variations in four neighboring positions, including position 97, which are located in two adjacent -strands on the floor of the peptide binding site of the molecule.The nucleotide sequence data reported in this paper have been submitted to the EMBL nucleotide sequence database and have been assigned the accession number X57954. Address correspondence and offprint requests to: J. A. López de Castro.     

Two distinct HLA-A * 0101-specific submotifs illustrate alternative peptide binding modes

 Previous studies have defined two different peptide binding motifs specific for HLA-A ^* 0101. These motifs are characterized by the presence of tyrosine (Y) at the C-termini of 9-mer and 10-mer peptides, and either a small polar or hydrophobic (S, T, M) residue in position 2, or a negatively charged (D or E) residue in position 3. In this study, the structural requirements for peptide binding to A ^* 0101 have been further analyzed by examining the binding capacity of large sets of peptides corresponding to naturally occurring sequences which bore one or the other of these two A ^* 0101-specific motifs. By correlating the presence of specific residue types at each position along the peptide sequence with increased (or decreased) binding affinity, the prominent influence of secondary anchor residues was revealed. In most cases, the two anchors in positions 2 and 3 appear to act synergistically. With the exception of the DE₃ submotif in 9-mer peptides, a positive role for aromatic residues in position 1 and the center of the peptide (positions 4 or 5 of 9- or 10-mer peptides, respectively), and proline at C-3, were also consistently detected. However, secondary anchor residues also appear to differ significantly between the two different submotifs, demonstrating that A ^* 0101 can utilize alternative modes in binding its peptide ligands. According to these analyses, specific refined submotifs were also established, and their merit verified by independent sets of potential A ^* 0101 binding peptides. Besides providing useful insight into the nature of the interaction of the A ^* 0101 allele with its peptide ligands, such refined motifs should also facilitate accurate prediction of potential A ^* 0101-restricted peptide epitopes. Received: 16 July 1996 / Revised: 18 September 1996     

Role of strong anchor residues in the effective binding of 10-mer and 11-mer peptides to HLA-A*2402 molecules

The binding capacity of one-hundred-and-seventy-two 8-mer to 11-mer peptides carrying HLA-A24 anchor residues to HLA-A^*2402 molecules was analyzed by using a HLA class I stabilization assay. Most (76.2%) of these peptides bound to HLA-A^*2402 molecules. These results confirmed previous findings that Tyr and Phe at P2 as well as Phe, Trp, Ile, and Leu at the C-terminus were main anchor residues for HLA-A^*2402. Tyr at P2 was a stronger anchor residue than Phe, while bulky aromatic hydrophobic residues Phe and Trp at the C-terminus are stronger anchors than aliphatic hydrophobic residues Ile and Leu. These results were also supported by an analysis using a panel of mutated 9-mer peptides at P2 and P9. Taken together, these results suggest that HLA-A^*2402 molecules have deep B- and F-pockets because they favor peptides carrying bulky aromatic hydrophobic residues at P2 and the C-terminus. The affinity of 8-mer peptides was significantly lower than that of 9-mer to 11-mer peptides, while there was no difference in affinity between 9-mer, 10-mer, and 11-mer peptides. The affinity of peptides carrying bulky aromatic hydrophobic residues at the C-terminus was higher than that of peptides carrying aliphatic hydrophobic residues in each of the 8-mer to 11-mer peptides, though the greatest difference in affinity was observed in 11-mer peptides. The strong interaction of side chains of these anchor residues with the corresponding pockets may permit the effective binding of 10-mer and 11-mer peptides to HLA-A^*2402 molecules.     

Evaluation of Human Leukocyte Antigen-A (HLA-A), Other Non-HLA Markers on Chromosome 6p21 and Risk of Nasopharyngeal Carcinoma

Background

The association between human leukocyte antigen (HLA) genes (located in the Major Histocompatibility Complex [MHC] region of chromosome 6p21) and NPC has been known for some time. Recently, two genome-wide association studies (GWAS) conducted in Taiwan and China confirmed that the strongest evidence for NPC association was mapped to the MHC region. It is still unclear, however, whether these findings reflect direct associations with Human Leukocyte Antigen (HLA) genes and/or to other genes in this gene-rich region.

Methods

To better understand genetic associations for NPC within the MHC region of chromosome 6, we conducted an evaluation that pooled two previously conducted NPC case-control studies in Taiwan (N = 591 cases and N = 521 controls). PCR-based genotyping was performed for 12 significant SNPs identified within 6p21 in the Taiwan NPC GWAS and for the HLA-A gene (exons 2 and 3).

Findings

After confirming homogeneity between the two studies, pooled odds ratios (OR) and 95% confidence intervals (CI) were estimated by logistic regression. We found that HLA-A (p-trend = 0.0006) and rs29232 (within the GABBR1 gene; p-trend = 0.005) were independent risk factors for NPC after adjustment for age, gender, study and each other. NPC risk was highest among individuals who were homozygous for the HLA-A*0207 risk allele and carriers of the rs29232 risk allele (A).

Conclusion

Our study suggests that most of the SNPs significantly associated with NPC from GWAS reflect previously identified HLA-A associations. An independent effect of rs29232 (GABBR1), however, remained, suggesting that additional genes within this region might be associated with NPC risk.     

New Frame Shift Mutation in Puroindoline B in Indian Wheat Cultivars Hyb65 and NI5439

Puroindoline genes pinA and pinB are the main components of the 15 kD friabilin protein reported to be associated with kernel softness. However, grain hardness of Hyb65 and NI5439, the two Indian wheat varieties, could not be explained based on the earlier identified alleles in puroindolines in wheat. Hyb65 and NI5439 are hard but based on the earlier identified allelic forms of puroindolines both the varieties could have been soft. In this investigation, puroindolines (a and b) from Hyb65 and NI5439 were characterised to understand their role in determining grain hardness. The sequence of puroindoline genes from both the varieties indicated that there was no mutation in pinA. However, there was frame shift mutation in pinB generated by insertion of a guanine residue 126 bp downstream from the start codon in both the varieties. This created new hardness allele of pinB designated as pinb-D1h. Frame shift also culminated into stop codon (TGA) 231 bp downstream from the start codon terminating protein synthesis at 77^th amino acid position. Five more stop codons (4TGA & 1TAG) were also created to the downstream positions of the first stop codon because of frame shift. There was additional point mutation in NI5439 (transition from A to G) resulting into change of amino acid residue from thymine to arginine at 205^th nucleotide position. Thus single nucleotide change in pinB resulted into truncated pin B and consequently the harder texture.     

A tissue culture induced Adh1 null mutant of maize results from a single base change

Summary We have found a null Adh1 allele which arose as a somaclonal variant following tissue culture of maize embryos carrying Adh1-1S and Adh1-1F alleles. Cloning and sequencing shows that the mutant allele derives from Adh1-1S and that there has been a single base change in the coding region of the gene which converts and AAG lysine codon to a TAG stop codon. The rate of nucleotide substitution (two per 218 embryos cultured) is much greater than normal mutation rates.     

Molecular basis for subtypic differences of the “a” subunit of coagulation factor XIII with description of the genesis of the subtypes

The “a” subunit of human coagulation factor XIII (F13A) exhibits genetic polymorphism defined by four common alleles, F13A*1A, *1B, *2A, and *2B. We have previously suggested on the basis of the isoelectric focusing patterns of the four allele products that point mutations at two separate sites and one intragenic crossing over might be involved in the genes of F13A polymorphism. Here, we report nucleotide substitutions associated with F13A polymorphism. A C/T transition of the second nucelotide of codon 564 in exon 12 is responsible for the difference between F13A*1A and *1B and that between F13A*2A and *2B, and a set of two base changes in codons 650 and 651 in exon 14 leads to the differences between F13A*1A and *2A and those between F13A*1B and *2B. The four combinations of the point mutations at the two exons thus correspond to the four alleles, two of which were generated by the point mutations from ancestral monomorphic gene. The results suggest strongly that intragenic crossing over must be involved in the genesis of the fourth allele. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods discriminating these base changes in exons 12 and 14 are also presented.     

HLA-G gene polymorphism in a Japanese population

 Polymorphism of the HLA-G gene in a Japanese population was investigated employing polymerase chain reaction (PCR)-single-strand conformation polymorphism (SSCP) analysis, PCR sequence-specific oligonucleotide (SSO) analysis, and DNA direct sequencing. Nucleotide sequence variations in exons 2, 3, and 4 of the HLA-G gene in 54 healthy Japanese individuals were examined. In addition, seven Japanese samples carrying common HLA haplotypes were analyzed. In total, nine single-base substitutions compared with the sequence of G ^* 01011 were identified: one in intron 1 (nucleotide position 970), one in exon 2 (the third base of codon 57: G → A), three in intron 2 (1264, 1276, and 1292), three in exon 3 (the third base of codon 93: C → T, the third base of codon 107: A → T, and the first base of codon 110: C → A), and one in intron 3 (2334). The substitution at codon 110 was non-synonymous and led to an amino acid substitution from leucine to isoleucine. The other three nucleotide substitutions in exons were synonymous. Through analysis of combinations of the exon 2, 3, and 4 nucleotide sequences we identified four alleles, which we provisionally designated GJ1, GJ2, GJ3, and GJ4. The allele frequencies were estimated to be 0.33, 0.16, 0.45, and 0.06, respectively. Nucleotide sequences of GJ1, GJ2, and GJ4 were identical to G ^* 01011, the clone 7.0E, and G ^* 01013, respectively. GJ3 was a newly observed allele and was officially designated G ^* 0104 by the WHO Nomenclature Committee in January 1996. Strong positive associations were observed between HLA-G alleles and HLA-A, -B, or -DRB1 alleles. Received: 15 February 1996 / Revised: 26 March 1996     

A rapid method for quantitative prediction of high affinity CTL epitopes: QSAR studies on peptides having affinity with the class I MHC molecular HLA-A*0201

Summary It would be useful to develop a method to rapidly identify peptide epitopes for vaccine development. In this paper, empirical three-dimensional quantitative structure-affinity relationship (3D-QSAR) methods were used to study the relationship between the three dimensional structural parameters (the isotropic surface area, ISA, and the electronic charge index, ECI) of the HLA-A^*0201 binding peptide and the HLA-A^*0201/peptide binding affinities. A set of 102 peptides having affinity with the class I MHC HLA-A^*0201 molecule was used as training set. A test set of 40 peptides was used to determine the predictive value of the models. The 3D-QSAR models gave aq ²=0.5724 and highr _pred ² =0.6955. According to the standard regression coefficients, it is known that the hydrophobic interactions (in these studies, the ISA is highly correlative with the hydrophobic property) play a dominant role in peptide-MHC molecule binding, and also which amino acid residue with what property is needed at specific position of the peptide. The approach we have taken is highly complementary to the many excellent methods described in references and appears highly predictive. It is a rapid and convenient method for detecting high affinity peptide epitopes.     

Molecular characterisation of the inactive allele of the gene Glu-A1 and the development of a set of AS-PCR markers for HMW glutenins of wheat

The present work reports new PCR markers that amplify the complete coding sequence of the specific alleles of the high molecular weight (HMW) glutenin genes. A set of AS-PCR molecular markers was designed which use primers from nucleotide sequences of the Glu-A1 and Glu-D1 genes, making use of the minor diffeences between the sequences of the x1, x2* of Glu-A1, and the x5 and y10 of Glu-D1. These primers were able to distinguish between x2* and the x1 or xNull of Glu-A1. Also x5 was distinguishable from x2, and y10 from y12. The primers amplified the complete coding regions and corresponded to the upstream and downstream flanking positions of Glu-A1 and Glu-D1. Primers designed to amplify the Glu-A1 gene amplified a single product when used with genomic DNA of common wheats and the xNull allele of this gene. This work also describes the cloning and characterisation of the nucleotide sequence of this allele. It possesses the same general structure as x2* and x1 (previously determined) and differs from these alleles in the extension of the coding sequence for a presumptive mature protein with only 384 residues. This is due to the presence of a stop codon (TAA) 1215-bp downstream from the start codon. A further stop codon (TAG), 2280-bp downstream from the starting codon is also found. The open reading frame of xNull and x1 alleles has the same size in bp. Both are larger than x2* which shows two small deletions. The reduced size of the presumptive mature protein encoded by xNull could explain the negative effect of this allele on grain quality. Received: 16 May 1999 / Accepted: 16 September 1999     

Mitochondrial rps14 is a transcribed and edited pseudogene in Arabidopsis thaliana

We have isolated and analysed a 2 kb region of the mitochondrial genome of Arabidopsis thaliana (Columbia) showing a high level of nucleotide identity with the mitochondrial (mt) rps14 small-subunit ribosomal protein gene from Oenothera berteriana and Vicia faba, as well as with an open reading frame (ORF) located upstream of the nad3 locus in O. berteriana. The rps14 locus is present as a single copy in the A. thaliana mt genome and has a translational stop codon located near the initiation codon, as well as a deletion of one nucleotide that disturbs the coding sequence. The cloning and sequencing of nine amplified mt rps14 cDNAs clearly demonstrated that this gene is transcribed and that the mRNA precursors are edited at three positions, all involving C-to-U conversions. No editing events changing the stop codon and restoring the correct coding sequence were witnessed within the 9 individual cDNA clones. Therefore, we conclude that the single rps14 sequence of the mitochondrial genome from A. thealiana is in fact a pseudogene that is transcribed and edited but not translated.     

New insights into the incorporation of natural suppressor tRNAs at stop codons in Saccharomyces cerevisiae

Stop codon readthrough may be promoted by the nucleotide environment or drugs. In such cases, ribosomes incorporate a natural suppressor tRNA at the stop codon, leading to the continuation of translation in the same reading frame until the next stop codon and resulting in the expression of a protein with a new potential function. However, the identity of the natural suppressor tRNAs involved in stop codon readthrough remains unclear, precluding identification of the amino acids incorporated at the stop position. We established an in vivo reporter system for identifying the amino acids incorporated at the stop codon, by mass spectrometry in the yeast Saccharomyces cerevisiae. We found that glutamine, tyrosine and lysine were inserted at UAA and UAG codons, whereas tryptophan, cysteine and arginine were inserted at UGA codon. The 5′ nucleotide context of the stop codon had no impact on the identity or proportion of amino acids incorporated by readthrough. We also found that two different glutamine tRNA^Gln were used to insert glutamine at UAA and UAG codons. This work constitutes the first systematic analysis of the amino acids incorporated at stop codons, providing important new insights into the decoding rules used by the ribosome to read the genetic code.     

Analysis of Low Frequency Protein Truncating Stop-Codon Variants and Fasting Concentration of Growth Hormone

Background

The genetic background of Growth Hormone (GH) secretion is not well understood. Mutations giving rise to a stop codon have a high likelihood of affecting protein function.

Objectives

To analyze likely functional stop codon mutations that are associated with fasting plasma concentration of Growth Hormone.

Methods

We analyzed stop codon mutations in 5451 individuals in the Malmö Diet and Cancer study by genotyping the Illumina Exome Chip. To enrich for stop codon mutations with likely functional effects on protein function, we focused on those disrupting >80% of the predicted amino acid sequence, which were carried by ≥10 individuals. Such mutations were related to GH concentration, measured with a high sensitivity assay (hs-GH) and, if nominally significant, to GH related phenotypes, using linear regression analysis.

Results

Two stop codon mutations were associated with the fasting concentration of hs-GH. rs121909305 (NP_005370.1:p.R93*) [Minor Allele Frequency (MAF) = 0.8%] in the Myosin 1A gene (MYO1A) was associated with a 0.36 (95%CI, 0.04 to 0.54; p=0.02) increment of the standardized value of the natural logarithm of hs-GH per 1 minor allele and rs35699176 (NP_067040.1:p.Q100*) in the Zink Finger protein 77 gene (ZNF77) (MAF = 4.8%) was associated with a 0.12 (95%CI, 0.02 to 0.22; p = 0.02) increase of hs-GH. The mutated high hs-GH associated allele of MYO1A was related to lower BMI (β-coefficient, -0.22; p = 0.05), waist (β-coefficient, -0.22; p = 0.04), body fat percentage (β-coefficient, -0.23; p = 0.03) and with higher HDL (β-coefficient, 0.23; p = 0.04). The ZNF77 stop codon was associated with height (β-coefficient, 0.11; p = 0.02) but not with cardiometabolic risk factors.

Conclusion

We here suggest that a stop codon of MYO1A, disrupting 91% of the predicted amino acid sequence, is associated with higher hs-GH and GH-related traits suggesting that MYO1A is involved in GH metabolism and possibly body fat distribution. However, our results are preliminary and need replication in independent populations.     

Different immunodominance of HIV-1-specific CTL epitopes among three subtypes of HLA-A*26 associated with slow progression to AIDS

It is speculated that HLA-A^∗26-restricted HIV-1-specific CTLs can control HIV-1, since HLA-A^∗26 is associated with a slow progression to AIDS. In three major HLA-A^∗26 subtypes, HLA-A^∗2601-restricted, and HLA-A^∗2603-restricted HIV-1 epitopes have been identified, but HLA-A^∗2602-restricted ones have not. We here identified HLA-A^∗2602-restricted HIV-1 epitopes by using reverse immunogenetics and compared the immunodominance of the epitopes among the three subtypes. Out of 110 HIV-1 peptides carrying HLA-A^∗26 anchor residues, only the Gag169-177 peptide, which had been previously identified as an HLA-A^∗2601- and HLA-A^∗2603-restricted immunodominant epitope, induced Gag169-177-specific CD8⁺ T cells from only two of six HLA-A^∗2602⁺ HIV-1-infected individuals. No difference in affinity of this epitope peptide was found among these three HLA-A^∗26 subtypes, indicating that Gag169-177 was effectively presented by HLA-A^∗2602 but recognized as a subdominant epitope in HIV-1-infected HLA-A^∗2602⁺ individuals. These findings indicate different immunodominance of Gag169-177 epitope among 3 HLA-A^∗26 subtypes.     

Genetic Integrity and Microgeographic Population Structure of Westslope Cutthroat Trout, Oncorhynchus clarki lewisi, in the Pend Oreille Basin in Washington

Five microsatellite DNA loci (Ots-101 ^*,Ots-107 ^*,Oki-10 ^*, Ogo-3 ^*, and FGT-3 ^*) were screened to evaluate the genetic characteristics and population structure for cutthroat trout from eight tributaries of the Pend Oreille River in northeastern Washington and to compare these collections with two hatchery stocks of westslope cutthroat trout, Oncorhynchus clarki lewisi, Yellowstone cutthroat trout, Oncorhynchus clarki bouvieri and a hatchery rainbow trout, Oncorhynchus mykiss, strain that have been stocked in northeastern Washington. Relatively high levels of variation (numbers of alleles and heterozygosity) were observed in all collections and allele frequencies were quite variable among collections. Evidence of limited introgression by rainbow and/or Yellowstone cutthroat was found at several locations. Both F_ST values and tests of genetic differentiation indicated the existence of numerous, reproductively isolated populations. The population in Slate Creek was very similar to the Kings Lake Hatchery strain, and we conclude that this similarity is the result of historical introductions of this hatchery strain into what was presumably a stream without a native cutthroat population. In one stream, differences in introgression and allele frequencies were found above and below a barrier falls. Because of the substantial level of population differentiation observed among the various collections, we recommend that management and conservation actions be focused at the level of individual streams in order to maintain the productivity and genetic character of the existing populations of cutthroat trout.     

Extensive polymorphism of the FUT2 gene in an African (Xhosa) population of South Africa

The human secretor type α(1,2)fucosyltrans-ferase gene (FUT2) polymorphism was investigated in Xhosa and Caucasian populations of South Africa by polymerase chain reaction–restriction fragment length polymorphism and DNA sequencing. Six new base substitutions were found in the coding region of FUT2. A single base (C) deletion at nucleotide 778, which led to a frame shift and produced a stop codon at codon 275, was responsible for the enzyme inactivation. Three nonsynonymous base substitutions, A40G (Ile¹⁴Val), C379T (Arg¹²⁷Cys), and G481A (Asp¹⁶¹Asn), and two synonymous base substitutions, A375G (Glu¹²⁵) and C480T (His¹⁶⁰), were also identified in functional alleles. As a result, seven new alleles, Se ⁴⁰ , Se ⁴⁸¹ , Se ^40,481 , Se ^357,480 , Se ^357,379,480 , Se ³⁷⁵ , and se ^357,480,778 were identified. Population studies revealed that an allele containing a nonsense mutation G428A (Trp¹⁴³stop) (se ⁴²⁸ ) was the common null allele in both Xhosa and Caucasian populations, whereas an allele containing a missense A385T (Ile¹²⁹Phe) mutation (se ^357,385 ), which is the common null allele in Orientals, was found to be absent from both populations. The heterozygosity rates of FUT2 genotypes were as high as 0.75 in the Xhosa population and 0.65 in the Caucasian population. Therefore, the extensive polymorphism and race specificity of the FUT2 gene make it suitable for application as a new tool in genetic studies of modern human evolutionary history. Received: 23 March 1998 / Accepted: 9 May 1998     

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.