Genotype List String: a grammar for describing HLA and KIR genotyping results in a text string.

Knowledge of an individual's human leukocyte antigen (HLA) genotype is essential for modern medical genetics, and is crucial for hematopoietic stem cell and solid-organ transplantation. However, the high levels of polymorphism known for the HLA genes make it difficult to generate an HLA genotype that unambiguously identifies the alleles that are present at a given HLA locus in an individual. For the last 20 years, the histocompatibility and immunogenetics community has recorded this HLA genotyping ambiguity using allele codes developed by the National Marrow Donor Program (NMDP). While these allele codes may have been effective for recording an HLA genotyping result when initially developed, their use today results in increased ambiguity in an HLA genotype, and they are no longer suitable in the era of rapid allele discovery and ultra-high allele polymorphism. Here, we present a text string format capable of fully representing HLA genotyping results. This Genotype List (GL) String format is an extension of a proposed standard for reporting killer-cell immunoglobulin-like receptor (KIR) genotype data that can be applied to any genetic data that use a standard nomenclature for identifying variants. The GL String format uses a hierarchical set of operators to describe the relationships between alleles, lists of possible alleles, phased alleles, genotypes, lists of possible genotypes, and multilocus unphased genotypes, without losing typing information or increasing typing ambiguity. When used in concert with appropriate tools to create, exchange, and parse these strings, we anticipate that GL Strings will replace NMDP allele codes for reporting HLA genotypes.

Multi-locus HLA class I and II allele and haplotype associations with follicular lymphoma.

Follicular lymphoma (FL) is an indolent, sometimes, fatal disease characterized by recurrence at progressively shorter intervals and is frequently refractive to therapy. Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) in the human leukocyte antigen (HLA) region on chromosome 6p21.32-33 that are statistically significantly associated with FL risk. Low to medium resolution typing of single or multiple HLA genes has provided an incomplete picture of the total genetic risk imparted by this highly variable region. To gain further insight into the role of HLA alleles in lymphomagenesis and to investigate the independence of validated SNPs and HLA alleles with FL risk, high-resolution HLA typing was conducted using next-generation sequencing in 222 non-Hispanic White FL cases and 220 matched controls from a larger San Francisco Bay Area population-based case-control study of lymphoma. A novel protective association was found between the DPB1*03:01 allele and FL risk [odds ratio (OR) = 0.39, 95% confidence interval (CI) = 0.21-0.68]. Extended haplotypes DRB1*01:01-DQA1*01:01-DQB1*05:01 (OR = 2.01, 95% CI = 1.22-3.38) and DRB1*15-DQA1*01-DQB1*06 (OR = 0.55, 95% CI = 0.36-0.82) also influenced FL risk. Moreover, DRB1*15-DQA1*01-DQB1*06 was highly correlated with an established FL risk locus, rs2647012. These results provide further insight into the critical roles of HLA alleles and SNPs in FL pathogenesis that involve multi-locus effects across the HLA region.

A multi-site study using high-resolution HLA genotyping by next generation sequencing.

The high degree of polymorphism at human leukocyte antigen (HLA) class I and class II loci makes high-resolution HLA typing challenging. Current typing methods, including Sanger sequencing, yield ambiguous typing results because of incomplete genomic coverage and inability to set phase for HLA allele determination. The 454 Life Sciences Genome Sequencer (GS FLX) next generation sequencing system coupled with conexio atf software can provide very high-resolution HLA genotyping. High-throughput genotyping can be achieved by use of primers with multiplex identifier (MID) tags to allow pooling of the amplicons generated from different individuals prior to sequencing. We have conducted a double-blind study in which eight laboratory sites performed amplicon sequencing using GS FLX standard chemistry and genotyped the same 20 samples for HLA-A, -B, -C, DPB1, DQA1, DQB1, DRB1, DRB3, DRB4, and DRB5 (DRB3/4/5) in a single sequencing run. The average sequence read length was 250 base pairs and the average number of sequence reads per amplicon was 672, providing confidence in the allele assignments. Of the 1280 genotypes considered, assignment was possible in 95% of the cases. Failure to assign genotypes was the result of researcher procedural error or the presence of a novel allele rather than a failure of sequencing technology. Concordance with known genotypes, in cases where assignment was possible, ranged from 95.3% to 99.4% for the eight sites, with overall concordance of 97.2%. We conclude that clonal pyrosequencing using the GS FLX platform and CONEXIO ATF software allows reliable identification of HLA genotypes at high resolution.

High-resolution, high-throughput HLA genotyping by next-generation sequencing.

The human leukocyte antigen (HLA) class I and class II loci are the most polymorphic genes in the human genome. Hematopoietic stem cell transplantation requires allele-level HLA typing at multiple loci to select the best matched unrelated donors for recipient patients. In current methods for HLA typing, both alleles of a heterozygote are amplified and typed or sequenced simultaneously, often making it difficult to unambiguously determine the sequence of the two alleles. Next-generation sequencing methods clonally propagate in parallel millions of single DNA molecules, which are then also sequenced in parallel. Recently, the read lengths obtainable by one such next-generation sequencing method (454 Life Sciences, Inc.) have increased to >250 nucleotides. These clonal read lengths make possible setting the phase of the linked polymorphisms within an exon and thus the unambiguous determination of the sequence of each HLA allele. Here we demonstrate this capacity as well as show that the throughput of the system is sufficiently high to enable a complete, 7-locus HLA class I and II typing for 24 or 48 individual DNAs in a single GS FLX sequencing run. Highly multiplexed amplicon sequencing is facilitated by the use of sample-specific internal sequence tags (multiplex identification tags or MIDs) in the primers that allow pooling of samples yet maintain the ability to assign sequences to specific individuals. We have incorporated an HLA typing software application developed by Conexio Genomics (Freemantle, Australia) that assigns HLA genotypes for these 7 loci (HLA-A, -B, -C, DRB1, DQA1, DQB1, DPB1), as well as for DRB3, DRB4, and DRB5 from 454 sequence data. The potential of this HLA sequencing system to analyze chimeric mixtures is demonstrated here by the detection of a rare HLA-B allele in a mixture of two homozygous cell lines (1/100), as well as by the detection of the rare nontransmitted maternal allele present in the blood of a severe combined immunodeficiency disease syndrome (SCIDS) patient.

Non-conservative substitutions distinguish previously uncharacterized HLA-A molecules.

The extent of class I HLA polymorphism is not yet realized, and to provide a glimpse of the HLA-A polymorphism which remains undetected, we have analyzed approximately 3,700 National Marrow Donor Program (NMDP) Donor/Recipient Pair Retrospective Study Samples with HLA-A DNA sequence-based typing (SBT). Seventeen new HLA-A alleles were detected, with a total of 19 nucleotide substitutions distinguishing these new alleles from their closest HLA-A relatives. Nearly all of the new alleles differ by single nucleotide substitutions; a majority of these substitutions can be explained by gene conversion events but 6 alleles likely originated by point mutation. Fifteen of the 19 nucleotide substitutions translate into amino acid differences in the molecule. Structurally, the inferred amino acid alterations were non-conservative in terms of chemical property, and most substitutions were positioned in 1 or more of the specificity pockets which determine peptide binding. Although these new alleles were identified in a primarily Caucasian sample population, 9 of the 17 new HLA-A alleles were found in samples of non-Caucasoid origin. A new allele detection rate of 1 in approximately 200 individuals in our data set would, therefore, be higher in a non-Caucasoid sample population. In summary, the single nucleotide substitutions that distinguish undetected HLA-A alleles translate into functionally distinct HLA-A molecules. Further studies of the role of HLA-A in transplantation, in disease association, and in evolution must therefore accommodate the discovery of new alleles differing by single nucleotides.

Evolution of Pacific/Asian populations inferred from HLA class II allele frequency distributions.

The allele frequency distributions for the HLA class II loci, DRB1, DQB1 and DPB1, in eight Pacific/Asian populations: Hawaiian, Samoan, Malay, Papua New Guinea (PNG) Highlands, and two Indonesian and PNG Lowland groups, were determined using high-resolution polymerase chain reaction/sequence-specific oligonucleotide probe (PCR/SSOP) typing methods. The allele frequency distributions for the HLA-DRB1 locus were determined for a third Indonesian population as well as for an additional Filipino population. DRB1 alleles in the DR2 serogroup (or allelic lineage) are very common in this region; in some populations, more than 50% of the alleles belong to this serogroup. The DRB1*1502 allele is frequent in nine of the ten populations studied, reaching a frequency of 0.48 in one Indonesian population and among Filipinos. Extensive DR-DQ haplotype diversity was detected in these populations. Seven different DR2-DQB1 haplotypes were observed in the Indonesian and PNG Lowland populations, eight in the PNG Highlands and ten in Malays and Filipinos. The DRB1*0410 allele, commonly observed in Australia, is observed in the PNG Highlands at a low frequency (f=0.03) and is absent in the other populations. Two additional DRB1 alleles commonly observed in Australia, DRB1*0405 and *1407, are also observed in the PNG Highlands at high frequencies (f=0.132 and 0.126), while they are rare in the PNG Lowlands (f=0.039 and 0.013). These alleles are generally rare or absent in the other populations. The DPB1*0501 allele, common in Chinese and Japanese populations, is most frequent in the Samoan, Hawaiian, Indonesian, and Malay populations, and the *0401 allele is the most frequent DPB1 allele in the PNG Lowlands. Both of these alleles have the same very high frequency (f=0.34) in the PNG Highlands. Analyses of homozygosity (the Ewens-Watterson F statistic) in these and other populations indicate that, while most allele frequency distributions are consistent with balancing selection, values of F for the Indonesian and Javan populations may reflect positive directional selection. Phylogenetic trees constructed using the allele frequencies at the DRB1 locus of the populations reported here, as well as those for additional Pacific, Asian, and Australian populations, indicate that the PNG Highland population is more closely related to Australian populations than to PNG Lowland populations, while the PNG Lowlands are more closely related to other Melanesian populations.

HLA class II haplotype associations with inflammatory bowel disease in Jewish (Ashkenazi) and non-Jewish caucasian populations.

Ulcerative colitis (UC) and Crohn's disease (CD) are the clinical entities comprising idiopathic inflammatory bowel disease (IBD). Previous studies on the association of IBD and human leukocyte antigen (HLA) class II genes suggested a role for HLA in this disease. Here we present HLA class II (DRB1, DQB1, DQA1, DPB1) allele and haplotype distributions determined using the polymerase chain reaction and sequence-specific oligonucleotide probe methods. A total of 578 UC and CD Caucasian patients and controls from Jewish (Ashkenazi) and non-Jewish populations was examined. Our previously reported association of DR1-DQ5 with CD was attributable to DRB1*0103. A dramatic association with IBD and the highly unusual DRB1*0103-DQA1*0501-DQB1*0301 haplotype (OR = 6.6, p = 0.036) was found. The more common DR1 haplotype, DRB1*0103-DQA1*0101-DQB1*0501, was also associated with IBD (OR = 3.1, p = 0.014), a result suggesting that interaction between DR and DQ may determine the extent of disease risk. Our previously reported association of DR2 with UC was attributable to DRB1*1502 (OR = 2.6, p = 0.006). At the DPB1 locus, a significant association of DPB1*0401 with CD was observed for the combined populations (OR = 1.85, p = 0.007). These observations indicate that some class II alleles and haplotypes confer susceptibility to both UC and CD, implying common immunogenetic mechanisms of pathogenesis, while others confer risk to only one of these diseases, and illustrate the value of DNA HLA typing in disease susceptibility analyses.

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations.

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.

Hardy-Weinberg testing for HLA class II (DRB1, DQA1, DQB1, and DPB1) loci in 26 human ethnic groups.

Testing the fit of population data to Hardy-Weinberg proportions is crucial in the validation of many current approaches in population genetic studies. In this paper, we tested fit to Hardy-Weinberg proportions using exact approaches for both the overall and individual heterozygote genotype data of four HLA Class II loci: DRB1, DQA1, DQB1, and DPB1, from 26 human populations. Eighty of 99 overall tests fit the Hardy-Weinberg expectation (73% for DRB1, 89% for DQA1, 81% for DQB1 and 81% for DPB1). Deviations from Hardy-Weinberg proportions were both locus and group specific. Although we could not rule out other mechanisms at work, the individual test results indicated that the departure was possibly partly due to recent admixture. Evidence for selection and other sources of deviation are also discussed.

Results of Expedicion Humana. I. Analysis of HLA class II (DRB1-DQA1-DPB1) alleles and DR-DQ haplotypes in nine Amerindian populations from Colombia.

HLA class II variation was analyzed in nine Native American populations of Colombia using PCR/SSOP typing methods. Under the auspices of the Expedition Humana, approximately 30 unrelated native Colombia Indian samples each from the Tule (NW Pacific Coast), Kogui (Sierra Nevada). Ijka (Sierra Nevada), Ingano (Amazonas), Coreguaje (Amazonas), Nukak (Amazonas), Waunana (Pacific), Embera (Pacific) and Sikuani (Northeastern Plains) were collected and analyzed at the DRBI, DQA1, DQB1 and DPB1 loci. The number of different DRB1, DQA1, DQB1 and DPB1 alleles in the Colombian Indians is markedly reduced in comparison with neighboring African Colombian populations, which exhibit a very high degree of class II variability, as discussed in an accompanying paper. In the Colombian Amerindian groups, DR2 (DRB1*1602), DR4 (DRB1*0407, *0404, *0403 AND *0411), DR6 (DRB1*1402) and DR8 (DRB1*0802) comprise > 95% of all DRB1 alleles. We also found an absence of DR3 in all populations, and DR1, DR7 and DR9 allelic groups were either very rare or absent. Each Colombian Amerindian population has a predominant DRB1 allele (f = approximately 0.22-0.65) and DRB1-DQA1-DQB1 haplotype. Several novel DR-DQ haplotypes were also found. At the DPB1 locus, DPB1*0402 (f = 0.28-0.82), *1401 (f = 0.03-0.45), and *3501 (f = 0.03-0.27), were the three most prevalent alleles, each population maintaining one of these three alleles as the predominant (f > 0.26) DPB1 allele. The reduction of diversity for the HLA class II alleles in the Colombian Indians is suggestive of a population bottleneck during the colonization of the Americans, with little to no subsequent admixture with neighboring African Colombian populations in the last approximately 300 years.

Results of Expedicion Humana. II. Analysis of HLA class II alleles in three African American populations from Colombia using the PCR/SSOP: identification of a novel DQB1*02 (*0203) allele.

PCR/SSOP typing methods were used to analyze the HLA Class II DRB1, DQA1, DQB1 and DPB1 loci of samples from three African American populations of Colombia. Forty samples from the Cauca (Pacific), and twenty samples each from the Choco (North Pacific Coast) and the Providencia (Caribbean island) populations, were collected and the Class II loci analyzed under the auspices of the Expedicion Humana. Despite the limited number of samples analyzed, the African Colombian populations exhibit a very high degree of class II polymorphism. A great diversity of DRB1 alleles was found, with representatives from all serological classes, including 19 DRB1 alleles in the Providencia, 16 in the Cauca and 14 in the Choco groups. In addition, a novel DQB1*02 allele (*0203) was found in two individuals from the Cauca population of the Pacific Coast. The sequence of the DQB1*0203 allele, associated with DR3, differs from DQB1*0201 by only one nucleotide substitution (C-->A) in the second position of codon 57, resulting in an Ala to Asp change. The addition of DQB1*0203 brings the total number of DQB1 alleles identified to date to 26. HLA class II diversity is much greater in these African Colombian populations than that seen in nearby Amerindian populations. Analysis of regional Colombian African American HLA population genetics is discussed with respect to the Colombian Amerindian HLA genetics described in an accompanying paper.

DNA-based HLA typing for cord blood stem cell transplantation.

Molecular analysis of the HLA loci has revealed a pattern of extensive sequence polymorphism. For the class II loci, the polymorphism is localized to the second exon, whereas for the class I loci, both the second and third exons are polymorphic. These polymorphic regions encode the peptide binding groove and appear to be functionally significant in terms of disease susceptibility and transplantation. However, much of this polymorphism cannot be detected by serologic HLA typing methods. DNA typing methods based on PCR amplification and hybridization with sequence-specific oligonucleotide (SSO) probes can distinguish the many allelic sequence variants identified at these loci. The use of arrays of immobilized SSO probes allows genetic typing at many polymorphic sequence motifs in a single PCR and single hybridization reaction, making possible the development of simple, robust, and automated tests. PCR-SSO probe typing of the HLA loci requires very little sample material, is capable of either general or fine discrimination of alleles, and can be used to detect maternal contamination of cord blood. The application of this approach to typing HLA class I and II loci is discussed with regard to hematopoietic transplantation therapy.

HLA class II linkage disequilibrium and haplotype evolution in the Cayapa Indians of Ecuador.

DNA-based typing of the HLA class II loci in a sample of the Cayapa Indians of Ecuador reveals several lines of evidence that selection has operated to maintain and to diversify the existing level of polymorphism in the class II region. As has been noticed for other Native American groups, the overall level of polymorphism at the DRB1, DQA1, DQB1, and DPB1 loci is reduced relative to that found in other human populations. Nonetheless, the relative evenness in the distribution of allele frequencies at each of the four loci points to the role of balancing selection in the maintenance of the polymorphism. The DQA1 and DQB1 loci, in particular, have near-maximum departures from the neutrality model, which suggests that balancing selection has been especially strong in these cases. Several novel DQA1-DQB1 haplotypes and the discovery of a new DRB1 allele demonstrate an evolutionary tendency favoring the diversification of class II alleles and haplotypes. The recombination interval between the centromeric DPB1 locus and the other class II loci will, in the absence of other forces such as selection, reduce disequilibrium across this region. However, nearly all common alleles were found to be part of DR-DP haplotypes in strong disequilibrium, consistent with the recent action of selection acting on these haplotypes in the Cayapa.

HLA-B alleles of the Cayapa of Ecuador: new B39 and B15 alleles.

Recent data suggest that HLA-B locus alleles can evolve quickly in native South American populations. To investigate further this phenomenon of new HLA-B variants among Amerindians, we studied samples from another South American tribe, the Cayapa from Ecuador. We selected individuals for HLA-B molecular typing based upon their HLA class II typing results. Three new variants of HLA-B39 and one new variant of HLA-B15 were found in the Cayapa: HLA-B*3905, HLA-B*3906, HLA-B*3907, and HLA-B*1522. A total of thirteen new HLA-B alleles have now been found in the four South American tribes studied. Each of these four tribes studied, including the Cayapa, had novel alleles that were not found in any of the other tribes, suggesting that many of these new HLA-B alleles may have evolved since the Paleo-Indians originally populated South America. Each of these 13 new alleles contained predicted amino acid replacements that were located in the peptide binding site. These amino acid replacements may affect the sequence motif of the bound peptides, suggesting that these new alleles have been maintained by selection. New allelic variants have been found for all common HLA-B locus antigenic groups present in South American tribes with the exception of B48. In spite of its high frequency in South American tribes, no evidence for variants of B48 has been found in all the Amerindians studied, suggesting that B48 may have unique characteristics among the B locus alleles.

Analysis of HLA class II haplotypes in the Cayapa Indians of Ecuador: a novel DRB1 allele reveals evidence for convergent evolution and balancing selection at position 86.

PCR amplification, oligonucleotide probe typing, and sequencing were used to analyze the HLA class II loci (DRB1, DQA1, DQB1, and DPB1) of an isolated South Amerindian tribe. Here we report HLA class II variation, including the identification of a new DRB1 allele, several novel DR/DQ haplotypes, and an unusual distribution of DPB1 alleles, among the Cayapa Indians (N = 100) of Ecuador. A general reduction of HLA class II allelic variation in the Cayapa is consistent with a population bottle-neck during the colonization of the Americas. The new Cayapa DRB1 allele, DRB1*08042, which arose by a G-->T point mutation in the parental DRB1*0802, contains a novel Val codon (GTT) at position 86. The generation of DRB1*08042 (Val-86) from DRB1*0802 (Gly-86) in the Cayapa, by a different mechanism than the (GT-->TG) change in the creation of DRB1*08041 (Val-86) from DRB1*0802 in Africa, implicates selection in the convergent evolution of position 86 DR beta variants. The DRB1*08042 allele has not been found in > 1,800 Amerindian haplotypes and thus presumably arose after the Cayapa separated from other South American Amerindians. Selection pressure for increased haplotype diversity can be inferred in the generation and maintenance of three new DRB1*08042 haplotypes and several novel DR/DQ haplotypes in this population. The DPB1 allelic distribution in the Cayapa is also extraordinary, with two alleles, DPB1*1401, a very rare allele in North American Amerindian populations, and DPB1*0402, the most common Amerindian DPB1 allele, constituting 89% of the Cayapa DPB1.(ABSTRACT TRUNCATED AT 250 WORDS)

Blackout detection as a multiobjective optimization problem.

We study new fast computational procedures for a pilot blackout (total loss of vision) detection in real time. Their validity is demonstrated by data acquired during experiments with volunteer pilots on a human centrifuge. A new systematic class of very fast suboptimal group filters is employed. The utilization of various inherent group invariancies of signals involved allows us to solve the detection problem via estimation with respect to many performance criteria. The complexity of the procedures in terms of the number of computer operations required for their implementation is investigated. Various classes of such prediction procedures are investigated, analyzed and trade offs are established. Also we investigated the validity of suboptimal filtering using different group filters for different performance criteria, namely: the number of false detections, the number of missed detections, the accuracy of detection and the closeness of all procedures to a certain bench mark technique in terms of dispersion squared (mean square error). The results are compared to recent studies of detection of evoked potentials using estimation. The group filters compare favorably with conventional techniques in many cases with respect to the above mentioned criteria. Their main advantage is the fast computational processing.

Visual evoked potential processing under acceleration stress.

A study of the use of visual evoked potentials to detect acceleration induced blackout is presented. Decisions concerning the presence or absence of the visual evoked potential within the measured EEG were made using detection and estimation-based approaches. The performance of each technique is evaluated in two ways: (a) via a performance index (PI) which is introduced in this work, and (b) by the accuracy of pure detection. Variations of the main processing techniques such as output averaging and decision thresholding were also investigated.