Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer.

INTRODUCTION: Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex. METHODS: In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa). RESULTS: Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer. CONCLUSION: Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

A spatio-temporal model of multi-marker antimalarial resistance.

The emergence and spread of drug-resistant Plasmodium falciparum parasites have hindered efforts to eliminate malaria. Monitoring the spread of drug resistance is vital, as drug resistance can lead to widespread treatment failure. We develop a Bayesian model to produce spatio-temporal maps that depict the spread of drug resistance, and apply our methods for the antimalarial sulfadoxine-pyrimethamine. We infer from genetic count data the prevalences over space and time of various malaria parasite haplotypes associated with drug resistance. Previous work has focused on inferring the prevalence of individual molecular markers. In reality, combinations of mutations at multiple markers confer varying degrees of drug resistance to the parasite, indicating that multiple markers should be modelled together. However, the reporting of genetic count data is often inconsistent as some studies report haplotype counts, whereas some studies report mutation counts of individual markers separately. In response, we introduce a latent multinomial Gaussian process model to handle partially reported spatio-temporal count data. As drug-resistant mutations are often used as a proxy for treatment efficacy, point estimates from our spatio-temporal maps can help inform antimalarial drug policies, whereas the uncertainties from our maps can help with optimizing sampling strategies for future monitoring of drug resistance.

Analysis of 12 X-STRs loci in a population from southeastern Brazil.

Short tandem repeat (STR) markers on the X chromosome have a high potential for solving complex kinship analysis and individual identification cases. To achieve such purposes, allele and haplotype frequencies for the specific population are necessary. Nonetheless, such frequencies are not always available. Therefore, we obtained haplotypes from 520 unrelated males from four different geographic regions of Espírito Santo - Brazil, using the Investigator Argus X-12 kit. Forensic parameters for linked groups of four X-STR loci are reported. Genetic distance analyzes suggest that ES population is genetically closer to the Italian population and farther from the Mexican one, among the populations analyzed in this study.

Association between Human Leukocyte Antigen and End-Stage Renal Disease in Patients from Transylvania, Romania.

End-stage renal disease (ESRD) is the final stage of chronic kidney disease. This study explored the association between human leukocyte antigen (HLA) and ESRD. The interaction between genetic and environmental factors may also play a role in the development of ESRD. The study included 2392 ESRD patients who were awaiting renal transplantation. Blood samples were genotyped by SSOP and SSP-PCR methods. Multivariate logistic regression analysis showed that HLA-A*11 (p = 0.027), HLA-A*34 (p = 0.017), HLA-A*69 (p = 0.012), HLA-B*41 (p < 0.001), HLA-B*50 (p = 0.004), HLA-DRB1*10 (p = 0.027), and HLA-DRB1*14 (p = 0.004) were positively associated with ESRD (OR > 1); HLA-DRB1*07 (p < 0.001), HLA-DRB1*08 (p = 0.005), and HLA-DRB1*13 (p < 0.001) were protective against ESRD (OR < 1); and the three-locus haplotype HLA-A*02-B*41-DRB1*03, containing one susceptible allele, was strongly associated with ESRD (p < 0.001, OR = 3.15). In conclusion, this retrospective analysis of HLA typing in patients with ESRD of various etiologies suggests that molecular data on the HLA polymorphism should be collected in order to identify high-risk ESRD patients and to improve graft survival after kidney transplantation.

HLA allele and haplotype frequencies of registered stem cell donors in Chile.

Patients in need of hematopoietic stem cell transplantation often rely on unrelated stem cell donors matched in certain human leukocyte antigen (HLA) genes. Donor search is complicated by the extensive allelic variability of the HLA system. Therefore, large registries of potential donors are maintained in many countries worldwide. Population-specific HLA characteristics determine the registry benefits for patients and also the need for further regional donor recruitment. In this work, we analyzed HLA allele and haplotype frequencies of donors of DKMS Chile, the first Chilean donor registry, with self-assessed "non-Indigenous" (n=92,788) and "Mapuche" (n=1,993) ancestry. We identified HLA alleles that were distinctly more abundant in the Chilean subpopulations than in worldwide reference populations, four of them particularly characteristic for the Mapuche subpopulation, namely B*39:09g, B*35:09, DRB1*04:07g, and DRB1*16:02g. Both population subsamples carried haplotypes of both Native American and European origin at high frequencies, reflecting Chile's complex history of admixture and immigration. Matching probability analysis revealed limited benefits for Chilean patients (both non-Indigenous and Mapuche) from donor registries of non-Chilean donors, thus indicating a need for ongoing significant donor recruitment efforts in Chile.

Extending the discrete Laplace method: incorporating multi-copy loci, partial repeats and null alleles.

The discrete Laplace method can be used to estimate the frequency of a Y-chromosomal STR haplotype using a random sample from the population. Two limitations of the method are the assumptions that each profile has exactly one allele at every locus and that this allele has an integer repeat number. We relax these assumptions to allow for multi-copy loci, partial repeats and null alleles. We show how the parameters to the extension of the model can be estimated by numerical optimisation using an off-the-shelf solver. Concordance with the discrete Laplace method is obtained when the data satisfy the more stringent assumptions of the original method. We also investigate the performance of the (extended) discrete Laplace method when used to assign match probabilities for haplotypes. A simulation study shows that as more loci are used, match probabilities are underestimated more severely. This is consistent with the hypothesis that the discrete Laplace method cannot model the matches that arise by being identical by descent (IBD). As the number of loci increases the fraction of matches that are IBD increases. Simulation provides support that the discrete Laplace can model those matches that arise from identity by state (IBS) only.

Haplotype-dependent HLA-DRB1-DQB1 susceptibility to occult HBV infection in Xi'an Han population.

BACKGROUND: Occult hepatitis B virus (HBV) infection (OBI) is primarily characterized by the persistence of HBV-DNA in the liver tissues and/or in the serum without detectable HBsAg. Human leukocyte antigen (HLA) polymorphisms have been found to be strongly associated with HBV in different ethnic backgrounds. The association of HLA-DRB1-DQB1 haplotypes with OBI has not been previously reported in China. The aim of this study was to identify the potential association of HLA-DRB1-DQB1 haplotypes that may be involved in OBI genetic susceptibility. METHODS: A case-control study was conducted between 107 OBI subjects and 280 healthy controls from the blood donors in the Shaanxi Province Blood Center. The HLA-DRB1, DQB1 loci were genotyped using polymerase chain reaction-sequence based typing (PCR-SBT). Based on the genotype data of the two loci, haplotype estimation was performed. RESULTS: HLA-DRB1*07:01-DQB1*02:02 (pc = 0.344 × 10-3 , OR = 3.489, 95%CI = 2.000-6.088) and HLA-DRB1*09:01-DQB1*03:03 (pc = 0.02, OR = 2.370, 95%CI = 1.450-3.873) serve as the possible risk and susceptibility haplotypes for OBI in Xi'an Han after Bonferroni correction. CONCLUSIONS: This study demonstrated that HLA II haplotypes were significantly associated with OBI in the Xi'an Han population. To the best of our knowledge, this is the first study to associate HLA-DRB1-DQB1 haplotypes with OBI, which can provide valuable insights into the relationship between the various genetic factors and immune responses in the Xi'an population. The findings can also form the basis for future studies about the role of HLA in OBI.

Inferring Signatures of Positive Selection in Whole-Genome Sequencing Data: An Overview of Haplotype-Based Methods.

Signatures of positive selection in the genome are a characteristic mark of adaptation that can reveal an ongoing, recent, or ancient response to environmental change throughout the evolution of a population. New sources of food, climate conditions, and exposure to pathogens are only some of the possible sources of selective pressure, and the rise of advantageous genetic variants is a crucial determinant of survival and reproduction. In this context, the ability to detect these signatures of selection may pinpoint genetic variants that are responsible for a significant change in gene regulation, gene expression, or protein synthesis, structure, and function. This review focuses on statistical methods that take advantage of linkage disequilibrium and haplotype determination to reveal signatures of positive selection in whole-genome sequencing data, showing that they emerge from different descriptions of the same underlying event. Moreover, considerations are provided around the application of these statistics to different species, their suitability for ancient DNA, and the usefulness of discovering variants under selection for biomedicine and public health in an evolutionary medicine framework.

Spatial heterogeneity of knockdown resistance mutations in the dengue vector Aedes albopictus in Guangzhou, China.

BACKGROUND: The city of Guangzhou has been the epicenter of dengue fever in China since the 1990s, with Aedes albopictus being the primary vector. The main method used to control vectors and prevent dengue fever has been the application of chemical insecticides; however, this control strategy has resulted in the development of resistance to these insecticides in mosquitoes. Here we report our investigation of the patterns of knockdown resistance (kdr) mutations in 15 field populations of Ae. albopictus collected from 11 districts in Guangzhou. RESULTS: Four mutant alleles (V1016G, F1534S, F1534C, F1534L) were detected in domain II and III of the voltage-gated sodium channel (VGSC) gene. Various allele frequencies of kdr mutations were observed (3.1-25.9% for V1016G, 22.6-85.5% for F1534S, 0-29.0% for F1534L, 0.6-54.2% for F1534C). Seven kdr haplotypes (VF, VS, VL, VC, GF, GC, GS) were identified; the highest frequency of haplotypes was found for the single mutant haplotype VS (50.8%), followed by the wild-type VF haplotype (21.7%) and the single mutant haplotype VC (11.9%). Of the three double mutant haplotypes, GF was the most frequent (8.8%), followed by GC (1.2%) and GS (0.8%). Aedes albopictus showed spatial heterogeneity in deltamethrin resistance in populations collected in Guangzhou. We also observed significant differences in haplotype frequency. The frequency of the VC haplotype was significantly higher in high-risk dengue areas than in low-risk ones. CONCLUSIONS: The kdr allele V1016G was discovered for the first time in Guangzhou. Genetic isolation in mosquito populations and long-term insecticide selection seem to be responsible for the persistent, patchy distribution of kdr mutant alleles. The small-scale spatial heterogeneity in the distribution and frequency of kdr mutations may have important implications for vector control operations and insecticide resistance management strategies.

A multi-dimensional evaluation of the 'NIST 1032' sample set across four forensic Y-STR multiplexes.

This manuscript reports Y-chromosomal short tandem repeat (Y-STR) haplotypes for 1032 male U.S. population samples across 30 Y-STR loci characterized by three capillary electrophoresis (CE) length-based kits (PowerPlex Y23 System, Yfiler Plus PCR Amplification Kit, and Investigator Argus Y-28 QS Kit) and one sequence-based kit (ForenSeq DNA Signature Prep Kit): DYF387S1, DYS19, DYS385 a/b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS505, DYS518, DYS522, DYS533, DYS549, DYS570, DYS576, DYS612, DYS627, DYS635, DYS643, and Y-GATA-H4. The length-based Y-STR haplotypes include six loci that are not reported in the sequence-based kit (DYS393, DYS449, DYS456, DYS458, DYS518, and DYS627), whereas three loci included in the sequence-based kit are not present in length-based kits (DYS505, DYS522, and DYS612). For the latter, a custom multiplex was used to generate CE length-based data, allowing 1032 samples to be evaluated for concordance across the 30 Y-STR loci included in these four commercial Y-STR typing kits. Discordances between typing methods were analyzed further to assess underlying causes such as primer binding site mutations and flanking region insertions/deletions. Allele-level frequency and statistical information is provided for sequenced loci, excluding the multi-copy loci DYF387S1 and DYS385 a/b, for which locus-specific haplotype-level frequencies are provided instead. The resulting data reveals the degree of information gained through sequencing: 88% of sequenced Y-STR loci contain additional sequence-based alleles compared to length-based data, with the DYS389II locus containing the most additional alleles (51) observed by sequencing. Despite these allelic increases, only minimal improvement was observed in haplotype resolution by sequence, with all four commercial kits providing a similar ability to differentiate length-based haplotypes in this sample set. Finally, a subset of 369 male samples were compared to their corresponding additionally sequenced father samples, revealing the sequence basis for the 50 length-based changes observed, and no additional sequence-based mutations. GenBank accession numbers are reported for each unique sequence, and associated records are available in the STRSeq Y-Chromosomal STR Loci National Center for Biotechnology Information (NCBI) BioProject, accession PRJNA380347. Haplotype data is updated in the Y-STR Haplotype Reference Database (YHRD) for the 'NIST 1032' data set to now achieve the level of maximal haplotype of YHRD. All supplementary files including revisions to previously published Y-STR data are available in the NIST Public Data Repository: U.S. population data for human identification markers, DOI 10.18434/t4/1500024.

Human leukocyte antigen class I association with occult hepatitis B virus infection in the Shaanxi Han group: Analysis at the haplotype level.

BACKGROUND: Human leucocyte antigen (HLA) class I genes and haplotypes correlate with hepatitis B virus (HBV) infection. Occult HBV infection (OBI), a special type of chronic HBV infection, is defined as HBV surface antigen negative patients with or those without serologic markers by the means of HBV DNA detection in human plasma or in liver tissue by a diagnostic test. So far, the associations of HLA I haplotypes with OBI have not been reported previously in China. METHODS: A case-control study between 107 OBI subjects and 280 healthy controls from blood donors in the Blood Center of the Shaanxi Province was conducted in the present association analysis. The HLA-A, -B and -C loci of case-control subjects were detected and genotyped by polymerase chain reaction-sequence based typing. The HLA-A, -B and -C haplotypic frequencies were calculated by the maximum likelihood method. RESULTS: The HLA-A*33:03-C*07:01G (pc = 0.039, odds ratio [OR] = 8.996, 95% confidence interval [CI] = 1.825-44.338), B*44:03-C*07:01G (pc = 0.0069, OR = 12.000, 95% CI = 2.507-57.436) and A*33:03-B*44:03-C*07:01G (pc = 0.04, OR = 7.094, 95% CI = 1.387-36.288) haplotypes showed a a significant positive association with OBI. Independent effects demonstrated that HLA-B*44:03 and HLA-C*07:01G gave the main contribution to risk, whereas HLA-A*33:03 was associated only by linkage disequilibrium. CONCLUSIONS: This present study is the first to demonstrate that HLA I haplotypes are associated with OBI in the Shaanxi Han population. The present results suggest that HLA-B*44:03-C*07:01G might be a potential risk factor for OBI. Comparisons of the frequencies of HLA I haplotypes at high resolution were made between OBI from northern Chinese Han and controls. Frequencies of HLA-A-B, A-C, B-C, A-B-C between the OBI group (gray column) and the control group (blank column) were compared. Haplotypes with p < 0.05 in either group were shown. p values for multiple comparisons (pc ) were corrected by Bonferroni correction (*pc  < 0.05).

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps.

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Estimating HLA haplotype frequencies from homozygous individuals - A Technical Report.

We estimated HLA haplotype frequencies based on individuals homozygous for 4, 5 or 6 loci. Validation of our approach using a sample of over 3.4 million German individuals was successful. Compared to an expectation-maximization algorithm, the errors were larger. However, our approach allows the unequivocal detection of rare haplotypes.

Availability of HLA-allele-matched unrelated donors and registry size: Estimation from haplotype frequency in the Italian population.

In Italy, an HLA-matched unrelated donor is currently the primary donor when a HLA matched sibling is not found for allogeneic haematopoietic stem cell transplantation (HSCT). Better outcomes for transplantation require optimal matching between donor and recipient at least at the HLA-A, -B, -C, and -DRB1 loci; therefore, the availability of HLA-matched unrelated donors is important. The enormous HLA polymorphism has always necessitated registries with a large number of individuals in order to be able to provide well-matched donors to a substantial percentage of patients. In order to increase the efficiency of the Italian Bone Marrow Donor Registry (IBMDR) in providing Italian patients with a suitable donor, the probability of finding an HLA-A, -B, -C, and -DRB1 allele-matched (8/8) or a single mismatch unrelated donor (7/8) was estimated in this study according to IBMDR size. Using a biostatistical approach based on HLA haplotype frequencies of more than 100,000 Italian donors enrolled in the IBMDR and HLA-typed at high-resolution level, the probability of finding an 8/8 HLA-matched donor was 23.8%; 33.4%; and 41.4% in simulated registry sizes of 200,000; 500,000; and 1,000,000 donors; respectively. More than 2 million recruited donors are needed to increase the likelihood of identifying an HLA 8/8 matched donor for 50% of Italian patients. If one single mismatch at HLA I class loci was accepted, the probability of finding a 7/8 HLA-matched donor was 62.8%; 73.7%; and 80.3% in 200,000 donors; 500,000; and 1,000,000 donors; respectively. Using the regional haplotype frequencies of IBMDR donors, the probability of recruiting a donor with a new HLA phenotype, in the different Italian regions, was also calculated. Our findings are highly relevant in estimating the optimal size of the national registry, in planning a cost-effective strategy for donor recruitment in Italy, and determining the regional priority setting of recruitment activity in order to increase the phenotypic variability of IBMDR as well as its efficiency.

High-resolution HLA allele and haplotype frequencies in several unrelated populations determined by next generation sequencing: 17th International HLA and Immunogenetics Workshop joint report.

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

Diversities of HLA-A, -B, -C, -DRB1 and -DQB1 loci in Chinese Kazak population and its genetic relatedness dissection with multiple populations: a comparative study.

Studying the allele and haplotype distributions of human leukocyte antigen (HLA) loci at 2nd-field level in different populations was important. Allele and haplotype frequencies of HLA-A, -B, -C, -DRB1 and -DQB1 loci in 110 unrelated healthy Kazak individuals living in Xinjiang (China) were analyzed using polymerase chain reaction sequence based typing. Thirty HLA-A, 48 HLA-B, 24 HLA-C, 34 HLA-DRB1 and 18 HLA-DQB1 alleles were detected at the 2nd-field level in the Kazak population. Frequencies of HLA alleles, genotypes, and haplotypes were calculated, and some exhibited significantly different distributions among different populations. A neighbor-joining (NJ) tree, heatmap, multidimensional scaling (MDS) and principal component analysis (PCA) were used to explore the genetic relationships between the Kazak population and 32 reference populations distributed in Asia, Africa, America and Europe using frequency data of HLA-A, -B, -C and -DRB1 loci. The NJ tree, heatmap, and MDS of the 33 populations were constructed based on pairwise DA values of populations obtained by the HLA-A, -B, -C and -DRB1 allele frequencies. Different PCA plots were constructed based on the allele frequencies of HLA-A, -B, -C and -DRB1 or estimated haplotypic frequencies of HLA-A, -B, -C loci. The data obtained in the present research can be used for research on HLA-related diseases or paternity relationships, and aid to finding the best matched donors in stem cell transplantation for Kazak individuals.

High-resolution HLA allele and haplotype frequencies of the Saudi Arabian population based on 45,457 individuals and corresponding stem cell donor matching probabilities.

We estimated HLA allele and haplotype frequencies of the Saudi Arabian population from a sample of 45,457 registered stem cell donors. The most frequent HLA alleles were A*02:01g (18.5%), C*06:02g (16.1%), B*51:01g (14.1%), DRB1*07:01g (16.2%), DQB1*02:01g (30.5%), and DPB1*04:01g (33.6%). The most frequent 5-locus haplotypes were A*02:05g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.73%), A*02:01g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.66%), and A*26:01g~C*07:02g~B*08:01g~DRB1*03:01g~DQB1*02:01g (1.38%). Furthermore, we used the calculated haplotype frequencies to estimate stem cell donor matching probabilities for Saudi Arabian donor and patient populations under various matching requirements. These results are relevant for strategic donor registry planning in the Kingdom of Saudi Arabia.

Study on the polymorphisms of HLA-ABCDQB1DRB1 alleles and haplotypes in Hubei Han population of China.

The present study aimed to analyse the frequencies of human leukocyte antigen HLA-ABCDQB1 and HLA-DRB1 alleles and haplotypes in a subset of 3,732 Han population from Hubei of China. All samples were typed in the HLA-ABCDQB1 and HLA-DRB1 loci using the sequence-based typing method; subsequently, the HLA polymorphisms were analysed. A total of 47 HLA-A, 89 HLA-B, 43 HLA-C, 49 HLA-DRB1 and 24 HLA-DQB1 alleles were identified in the Hubei Han population. The top three most frequent alleles in the HLA-ABCDQB1 and HLA-DRB1 were A*11:01 (0.2617), A*24:02 (0.1590), A*02:07 (0.1281); B*46:01 (0.1502), B*40:01 (0.1409) and B*58:01 (0.0616); C*01:02 (0.2023), C*07:02 (0.1691) and C*03:04 (0.1175); and DQB1*03:01 (0.2000), DQB1*03:03 (0.1900), DQB1*06:01 (0.1187); DRB1*09:01 (0.1790), DRB1*15:01 (0.1062) and DRB1*12:02 (0.0841), respectively. Meanwhile, the three most frequent two-loci haplotypes were A*02:07-C*01:02 (0.0929), B*46:01-C*01:02 (0.1366) and DQB1*03:03-DRB1*09:01 (0.1766). The three most frequent three-loci haplotypes were A*02:07-B*46:01-C*01:02 (0.0883), B*46:01-DQB1*03:03-DRB1*09:01 (0.0808) and C*01:02-DQB1*03:03-DRB1*09:01 (0.0837). The three most frequent four-loci haplotypes were A*02:07-B*46:01-C*01:02-DQB1*03:03 (0.0494), B*46:01-DRB1*09:01-C*01:02-DQB1*03:03 (0.0729) and A*02:07-B*46:01-DQB1*03:03-DRB1*09:01 (0.0501). The most frequent five-loci haplotype was A*02:07-B*46:01-C*01:02-DQB1*03:03-DRB1*09:01 (0.0487). Heat maps and multiple correspondence analysis based on the frequencies of HLA specificity indicated that the Hubei Han population might be described into Southern Chinese populations. Our results lay a certain foundation for future population studies, disease association studies and donor recruitment strategies.

Identification of 8-Digit HLA-A, -B, -C, and -DRB1 Allele and Haplotype Frequencies in Koreans Using the One Lambda AllType Next-Generation Sequencing Kit.

BACKGROUND: Recent studies have successfully implemented next-generation sequencing (NGS) in HLA typing. We performed HLA NGS in a Korean population to estimate HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies up to an 8-digit resolution, which might be useful for an extended application of HLA results. METHODS: A total of 128 samples collected from healthy unrelated Korean adults, previously subjected to Sanger sequencing for 6-digit HLA analysis, were used. NGS was performed for HLA-A, -B, -C, and -DRB1 using the AllType NGS kit (One Lambda, West Hills, CA, USA), Ion Torrent S5 platform (Thermo Fisher Scientific, Waltham, MA, USA), and Type Steam Visual NGS analysis software (One Lambda). RESULTS: Eight HLA alleles showed frequencies of ≥10% in the Korean population, namely, A*24:02:01:01 (19.5%), A*33:03:01 (15.6%), A*02:01:01:01 (14.5%), A*11:01:01:01 (13.3%), B*15:01:01:01 (10.2%), C*01:02:01 (19.9%), C*03:04:01:02 (11.3%), and DRB1*09:01:02 (10.2%). Nine previous 6-digit HLA alleles were further identified as two or more 8-digit HLA alleles. Of these, eight alleles (A*24:02:01, B*35:01:01, B*40:01:02, B*55:02:01, B*58:01:01, C*03:02:02, C*07:02:01, and DRB1*07:01:01) were identified as two 8-digit HLA alleles, and one allele (B*51:01:01) was identified as three 8-digit HLA alleles. The most frequent four-loci haplotype was HLA-A*33:03:01-B*44:03:01:01-C*14:03-DRB1*13:02:01. CONCLUSIONS: We identified 8-digit HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a healthy Korean population using NGS. These new data can be used as a representative Korean data for further disease-related HLA type analysis.

Distribution and haplotype polymorphism of HLA-A、-B、-C、-DRB1、-DQB1 alleles in Anhui Han population / 中国输血杂志

【Objective】 To study the distribution and haplotype polymorphism of HLA-A, -B, -C, -DRB1, -DQB1 alleles in Anhui Han population. 【Methods】 The HLA-A, -B, -C, -DRB1 and -DQB1 genotyping of 3 169 random unrelated stem cell donors was performed by PCR-SBT. The allele frequency, haplotype frequency and linkage imbalance parameters were calculated by counting method, maximum expectation algorithm and PyPop software. 【Results】 A total of 411 HLA alleles were detected in the population, of which 67, 143, 65, 75 and 64 alleles were detected for HLA-A, -B, -C, -DRB1 and -DQB1, respectively. The alleles with frequency >0.1 were HLA-A*11∶01, A*11∶01, A*24∶02, A*02∶01, C*01∶02, C*07∶02, C*06∶02, DRB1*09∶01, DRB1*15∶01, DRB1*07∶01, DQB1* 03∶01, DQB1* 03∶03, and DQB1*02∶01. 1426 HLA-A~HLA-B, 1 772 HLA-B~HLA-DRB1, 798 HLA-B~HLA-C, and 446 HLA-DRB1~HLA-DQB1 haplotypes were detected. The haplotypes showed linkage imbalance, and 19 of them showed strong linkage imbalance (RLD>0.80). 【Conclusion】 The frequency and haplotype distribution of HLA-A, -B, -C, -DRB1 and -DQB1 alleles in Anhui Han population were obtained. The distribution of those alleles and haplotypes have their own characteristics.