Search

Polymorphisms in human leukocyte antigen (HLA) class I loci are known to have a great impact on disease progression in HIV-1 infection. Prevailing HIV-1 subtypes and HLA genotype distribution are different all over the world, and the HIV-1 and host HLA interaction could be specific to individual areas. Data on the HIV-1 and HLA interaction have been accumulated in HIV-1 subtype B- and C-predominant populations but not fully obtained in West Africa where HIV-1 subtype CRF02_AG is predominant. In the present study, to obtain accurate HLA typing data for analysis of HLA association with disease progression in HIV-1 infection in West African populations, HLA class I (HLA-A, -B, and -C) four-digit allele typing was performed in treatment-naïve HIV-1 infected individuals in Ghana (n = 324) by a super high-resolution single-molecule sequence-based typing (SS-SBT) using next-generation sequencing. Comparison of the SS-SBT-based data with those obtained by a conventional sequencing-based typing (SBT) revealed incorrect assignment of several alleles by SBT. Indeed, HLA-A*23:17, HLA-B*07:06, HLA-C*07:18, and HLA-C*18:02 whose allele frequencies were 2.5%, 0.9%, 4.3%, and 3.7%, respectively, were not determined by SBT. Several HLA alleles were associated with clinical markers, viral load and CD4+ T-cell count. Of note, the impact of HLA-B*57:03 and HLA-B*58:01, known as protective alleles against HIV-1 subtype B and C infection, on clinical markers was not observed in our cohort. This study for the first time presents SS-SBT-based four-digit typing data on HLA-A, -B, and -C alleles in Ghana, describing impact of HLA on viral load and CD4 count in HIV-1 infection. Accumulation of these data would facilitate high-resolution HLA genotyping, contributing to our understanding of the HIV-1 and host HLA interaction in Ghana, West Africa.

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Sequences homologous to human herpesvirus 6 (HHV-6) are integrated within the nuclear genome of about 1% of humans, but it is not clear how this came about. It is also uncertain whether integrated HHV-6 can reactivate into an infectious virus. HHV-6 integrates into telomeres, and this has recently been associated with polymorphisms affecting MOV10L1. MOV10L1 is located on the subtelomere of chromosome 22q (chr22q) and is required to make PIWI-interacting RNAs (piRNAs). As piRNAs block germline integration of transposons, piRNA-mediated repression of HHV-6 integration has been proposed to explain this association. In vitro, recombination of the HHV-6 genome along its terminal direct repeats (DRs) leads to excision from the telomere and viral reactivation, but the expected "solo-DR scar" has not been described in vivo. Here we screened for integrated HHV-6 in 7,485 Japanese subjects using whole-genome sequencing (WGS). Integrated HHV-6 was associated with polymorphisms on chr22q. However, in contrast to prior work, we find that the reported MOV10L1 polymorphism is physically linked to an ancient endogenous HHV-6A variant integrated into the telomere of chr22q in East Asians. Unexpectedly, an HHV-6B variant has also endogenized in chr22q; two endogenous HHV-6 variants at this locus thus account for 72% of all integrated HHV-6 in Japan. We also report human genomes carrying only one portion of the HHV-6B genome, a solo-DR, supporting in vivo excision and possible viral reactivation. Together these results explain the recently-reported association between integrated HHV-6 and MOV10L1/piRNAs, suggest potential exaptation of HHV-6 in its coevolution with human chr22q, and clarify the evolution and risk of reactivation of the only intact (non-retro)viral genome known to be present in human germlines.

Parasite virulence, an important factor contributing to the severity of Plasmodium falciparum infection, varies among P. falciparum strains. Relatively little is known regarding markers of virulence capable of identifying strains responsible for severe malaria. We investigated the effects of genetic variations in the P.f. merozoite surface protein 2 gene (msp2) on virulence, as it was previously postulated as a factor. We analyzed 300 msp2 sequences of single P. falciparum clone infection from patients with uncomplicated disease as well as those admitted for severe malaria with and without cerebral disease. The association of msp2 variations with disease severity was examined. We found that the N allele at codon 8 of Block 2 in the FC27-like msp2 gene was significantly associated with severe disease without cerebral complications (odds ratio = 2.73, P = 0.039), while the K allele at codon 17 of Block 4 in the 3D7-like msp2 gene was associated with cerebral malaria (odds ratio = 3.52, P = 0.024). The data suggests possible roles for the associated alleles on parasite invasion processes and immune-mediated pathogenicity. Multiplicity of infection was found to associate with severe disease without cerebral complications, but not cerebral malaria. Variations in the msp2-FC27-block 2-8N and 3D7-block 4-17K allele appear to be parasite virulence markers, and may be useful in determining the likelihood for severe and cerebral malaria. Their interactions with potential host factors for severe diseases should also be explored.

[This corrects the article DOI: 10.1371/journal.pone.0190418.].

The people of the Solomon Islands represent an Austronesian (AN)-speaking population's adaptation to a humid tropical environment and subsistence of tuberous crops. Genome-wide association studies (GWASs) of other populations (e.g. the Human Genome Diversity Project [HGDP]) have suggested the existence of genotypes adaptive to ecoregion, diet, and subsistence, and that those genotypes are also associated with metabolic and cardiovascular diseases. Recently, the incidence of non-communicable diseases has been increasing in the Solomon Islands. In the present study, we explored the association of genotypes adaptive to a tropical environment and tuberous crop diet with metabolic and cardiovascular conditions in rural and urban AN-speaking Melanesian and Micronesian populations of the Solomon Islands. A total of 561 participants were genotyped for single nucleotide polymorphisms (SNPs) potentially associated with a tropical environment (rs174570 and rs2237892) and a tuberous crop diet (rs162036, rs185819, and rs2722425). The results showed that the allele frequencies of the Solomon Islands populations adopted patterns similar to those in populations from other hot, tropical areas with a tuberous crop diet in previous studies. Furthermore, rs162036, rs185819, rs2237892, and rs2722425 were all strongly associated with one or more metabolic and cardiovascular conditions. The derived allele of rs2722425 (i.e. rs2722425-G) was significantly associated with an elevated LDL level (P = 0.000264) even after the significance level was adjusted for multiple testing (i.e., α = 0.0005). Our results suggest that the inhabitants of the Solomon Islands exhibit the effects of the tropical environment and tuberous crop diet on their allele frequencies, and that their susceptibility to metabolic and cardiovascular diseases is therefore considered to be associated with their environment and diet.

Previous studies have revealed the association between SNPs located on human leukocyte antigen (HLA) class II genes, including HLA-DP and HLA-DQ, and chronic hepatitis B virus (HBV) infection, mainly in Asian populations. HLA-DP alleles or haplotypes associated with chronic HBV infection or disease progression have not been fully identified in Asian populations. We performed trans-ethnic association analyses of HLA-DPA1, HLA-DPB1 alleles and haplotypes with hepatitis B virus infection and disease progression among Asian populations comprising Japanese, Korean, Hong Kong, and Thai subjects. To assess the association between HLA-DP and chronic HBV infection and disease progression, we conducted high-resolution (4-digit) HLA-DPA1 and HLA-DPB1 genotyping in a total of 3,167 samples, including HBV patients, HBV-resolved individuals and healthy controls. Trans-ethnic association analyses among Asian populations identified a new risk allele HLA-DPB1*09 ∶ 01 (P = 1.36 × 10(-6); OR= 1.97; 95% CI, 1.50-2.59) and a new protective allele DPB1*02 ∶ 01 (P = 5.22 × 10(-6); OR = 0.68; 95% CI, 0.58-0.81) to chronic HBV infection, in addition to the previously reported alleles. Moreover, DPB1*02 ∶ 01 was also associated with a decreased risk of disease progression in chronic HBV patients among Asian populations (P = 1.55 × 10(-7); OR = 0.50; 95% CI, 0.39-0.65). Trans-ethnic association analyses identified Asian-specific associations of HLA-DP alleles and haplotypes with HBV infection or disease progression. The present findings will serve as a base for future functional studies of HLA-DP molecules in order to understand the pathogenesis of HBV infection and the development of hepatocellular carcinoma.

Sporozoites of Plasmodium falciparum are transmitted to human hosts by Anopheles mosquitoes. Thrombospondin-related adhesive protein (TRAP) is expressed in sporozoites and plays a crucial role in sporozoite gliding and invasion of human hepatocytes. A previous study showed that the TRAP gene has been subjected to balancing selection in the Gambian P. falciparum population. To further study the molecular evolution of the TRAP gene in Plasmodium falciparum, we investigated TRAP polymorphisms in P. falciparum isolates from Suan Phueng District in Ratchaburi Province, Thailand. The analysis of the entire TRAP coding sequences in 32 isolates identified a total of 39 single nucleotide polymorphisms (SNPs), which comprised 37 nonsynonymous and two synonymous SNPs. McDonald-Kreitman test showed that the ratio of the number of nonsynonymous to synonymous polymorphic sites within P. falciparum was significantly higher than that of the number of nonsynonymous to synonymous fixed sites between P. falciparum and P. reichenowi. Furthermore, the rate of nonsynonymous substitution was significantly higher than that of synonymous substitution within Thai P. falciparum. These results indicate that the TRAP gene has been subject to diversifying selection in the Thai P. falciparum population as well as the Gambian P. falciparum population. Comparison of our P. falciparum isolates with those from another region of Thailand (Tak province, Thailand) revealed that TRAP was highly differentiated between geographically close regions. This rapid diversification seems to reflect strong recent positive selection on TRAP. Our results suggest that the TRAP molecule is a major target of the human immune response to pre-erythrocytic stages of P. falciparum.

The human leukocyte antigen (HLA) genes exhibit the highest degree of polymorphism in the human genome. This high degree of variation at classical HLA class I and class II loci has been maintained by balancing selection for a long evolutionary time. However, little is known about recent positive selection acting on specific HLA alleles in a local population. To detect the signature of recent positive selection, we genotyped six HLA loci, HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQB1, and HLA-DPB1 in 418 Japanese subjects, and then assessed the haplotype homozygosity (HH) of each HLA allele. There were 120 HLA alleles across the six loci. Among the 80 HLA alleles with frequencies of more than 1%, DPB1*04∶01, which had a frequency of 6.1%, showed exceptionally high HH (0.53). This finding raises the possibility that recent positive selection has acted on DPB1*04∶01. The DPB1*04∶01 allele, which was present in the most common 6-locus HLA haplotype (4.4%), A*33∶03-C*14∶03-B*44∶03-DRB1*13∶02-DQB1*06∶04-DPB1*04∶01, seems to have flowed from the Korean peninsula to the Japanese archipelago in the Yayoi period. A stochastic simulation approach indicated that the strong linkage disequilibrium between DQB1*06∶04 and DPB1*04∶01 observed in Japanese cannot be explained without positive selection favoring DPB1*04∶01. The selection coefficient of DPB1*04∶01 was estimated as 0.041 (95% credible interval 0.021-0.077). Our results suggest that DPB1*04∶01 has recently undergone strong positive selection in Japanese population.

Cerebral malaria is a major, life-threatening complication of Plasmodium falciparum malaria, and has very high mortality rate. In murine malaria models, natural killer (NK) cell responses have been shown to play a crucial role in the pathogenesis of cerebral malaria. To investigate the role of NK cells in the developmental process of human cerebral malaria, we conducted a case-control study examining genotypes for killer immunoglobulin-like receptors (KIR) and their human leukocyte antigen (HLA) class I ligands in 477 malaria patients. We found that the combination of KIR2DL3 and its cognate HLA-C1 ligand was significantly associated with the development of cerebral malaria when compared with non-cerebral malaria (odds ratio 3.14, 95% confidence interval 1.52-6.48, P = 0.00079, corrected P = 0.02). In contrast, no other KIR-HLA pairs showed a significant association with cerebral malaria, suggesting that the NK cell repertoire shaped by the KIR2DL3-HLA-C1 interaction shows certain functional responses that facilitate development of cerebral malaria. Furthermore, the frequency of the KIR2DL3-HLA-C1 combination was found to be significantly lower in malaria high-endemic populations. These results suggest that natural selection has reduced the frequency of the KIR2DL3-HLA-C1 combination in malaria high-endemic populations because of the propensity of interaction between KIR2DL3 and C1 to favor development of cerebral malaria. Our findings provide one possible explanation for KIR-HLA co-evolution driven by a microbial pathogen, and its effect on the global distribution of malaria, KIR and HLA.

INTRODUCTION: Interstitial Lung Disease (ILD) is frequently associated with Rheumatoid Arthritis (RA) as one of extra-articular manifestations. Many studies for Human Leukocyte Antigen (HLA) allelic association with RA have been reported, but few have been validated in an RA subpopulation with ILD. In this study, we investigated the association of HLA class II alleles with ILD in RA. METHODS: An association study was conducted on HLA-DRB1, DQB1, and DPB1 in 450 Japanese RA patients that were or were not diagnosed with ILD, based on the findings of computed tomography images of the chest. RESULTS: Unexpectedly, HLA-DRB1*04 (corrected P [Pc] = 0.0054, odds ratio [OR] 0.57), shared epitope (SE) (P = 0.0055, OR 0.66) and DQB1*04 (Pc = 0.0036, OR 0.57) were associated with significantly decreased risk of ILD. In contrast, DRB1*16 (Pc = 0.0372, OR 15.21), DR2 serological group (DRB1*15 and *16 alleles) (P = 0.0020, OR 1.75) and DQB1*06 (Pc = 0.0333, OR 1.57, respectively) were significantly associated with risk of ILD. CONCLUSION: HLA-DRB1 SE was associated with reduced, while DR2 serological group (DRB1*15 and *16) with increased, risk for ILD in Japanese patients with RA.

The human dopamine receptor D4 (DRD4) gene contains a 48-bp variable number of tandem repeat (VNTR) in exon 3, encoding the third intracellular loop of this dopamine receptor. The DRD4 7R allele, which seems to have a single origin, is commonly observed in various human populations and the nucleotide diversity of the DRD4 7R haplotype at the DRD4 locus is reduced compared to the most common DRD4 4R haplotype. Based on these observations, previous studies have hypothesized that positive selection has acted on the DRD4 7R allele. However, the degrees of linkage disequilibrium (LD) of the DRD4 7R allele with single nucleotide polymorphisms (SNPs) outside the DRD4 locus have not been evaluated. In this study, to re-examine the possibility of recent positive selection favoring the DRD4 7R allele, we genotyped HapMap subjects for DRD4 VNTR, and conducted several neutrality tests including long range haplotype test and iHS test based on the extended haplotype homozygosity. Our results indicated that LD of the DRD4 7R allele was not extended compared to SNP alleles with the similar frequency. Thus, we conclude that the DRD4 7R allele has not been subjected to strong recent positive selection.

Insulin resistance is a feature of type 2 diabetes. Resistin, secreted from adipocytes, causes insulin resistance in mice. We previously reported that the G/G genotype of single nucleotide polymorphism (SNP) at -420 (rs1862513) in the human resistin gene (RETN) increased susceptibility to type 2 diabetes by enhancing its promoter activity. Plasma resistin was highest in Japanese subjects with G/G genotype, followed by C/G, and C/C. In this study, we cross-sectionally analyzed plasma resistin and SNPs in the RETN region in 2,019 community-dwelling Japanese subjects. Plasma resistin was associated with SNP-638 (rs34861192), SNP-537 (rs34124816), SNP-420, SNP-358 (rs3219175), SNP+299 (rs3745367), and SNP+1263 (rs3745369) (P

Phenotypic divergences between modern human populations have developed as a result of genetic adaptation to local environments over the past 100,000 years. To identify genes involved in population-specific phenotypes, it is necessary to detect signatures of recent positive selection in the human genome. Although detection of elongated linkage disequilibrium (LD) has been a powerful tool in the field of evolutionary genetics, current LD-based approaches are not applicable to already fixed loci. Here, we report a method of scanning for population-specific strong selective sweeps that have reached fixation. In this method, genome-wide SNP data is used to analyze differences in the haplotype frequency, nucleotide diversity, and LD between populations, using the ratio of haplotype homozygosity between populations. To estimate the detection power of the statistics used in this study, we performed computer simulations and found that these tests are relatively robust against the density of typed SNPs and demographic parameters if the advantageous allele has reached fixation. Therefore, we could determine the threshold for maintaining high detection power, regardless of SNP density and demographic history. When this method was applied to the HapMap data, it was able to identify the candidates of population-specific strong selective sweeps more efficiently than the outlier approach that depends on the empirical distribution. This study, confirming strong positive selection on genes previously reported to be associated with specific phenotypes, also identifies other candidates that are likely to contribute to phenotypic differences between human populations.