Search

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) on human chromosome 19q13.4 encode 11 immunoglobulin superfamily receptors, exhibiting genetic diversity within and between human populations. Among the LILR genes, the genomic region surrounding LILRB3 and LILRA6 has yet to be fully characterized due to their significant sequence homology, which makes it difficult to differentiate between them. To examine the LILRB3 and LILRA6 genomic region, a tool named JoGo-LILR CN Caller, which can call copy number from short-read whole genome sequencing (srWGS) data, was applied to an extensive international srWGS dataset comprising 2,504 samples. During this process, a previously unreported loss of both LILRB3 and LILRA6 was detected in three samples. Using long-read sequencing of these samples, we have discovered a novel large deletion (33,692 bp) in the LILRB3 and LILRA6 genomic regions in the Japanese population. This deletion spanned three genes, LILRB3, LILRA6, and LILRB5, resulting in LILRB3 exons 12-13 being located immediately downstream of LILRB5 exons 1-12 with the loss of LILRA6, suggesting the potential expression of a hybrid gene between LILRB5 and LILRB3 (LILRB5-3). Transcription and subsequent translation of the LILRB5-3 hybrid gene were also verified. The hybrid junction was located within the intracellular domain, resulting in an LILRB5 extracellular domain fused to a partial LILRB3 intracellular domain with three immunoreceptor tyrosine-based inhibitory motifs (ITIMs), suggesting that LILRB5-3 acquired a novel signaling function. Further application of the JoGo-LILR tool to srWGS samples suggested the presence of the LILRB5-3 hybrid gene in the CEU population. Our findings provide insight into the genetic and functional diversity of the LILR family.

ABSTRACT Aim/Introduction Resistin, which induces insulin resistance, is mainly expressed in monocytes/macrophages in humans. We reported previously that serum resistin was highest in the G–A haplotype defined by resistin single nucleotide polymorphisms (SNPs) at −420 (rs1862513) and − 358 (rs3219175). As sarcopenic obesity is associated with insulin resistance, we aimed to examine whether serum resistin and its haplotypes were associated with sarcopenic obesity at a latent stage. Materials and Methods We cross‐sectionally analyzed 567 community‐dwelling Japanese participants attending annual medical check‐ups in which the sarcopenic obesity index was evaluated. The age‐ and gender‐matched normal glucose tolerance subjects with G–A homozygotes and those with C–G homozygotes were examined via RNA‐sequencing and pathway analysis (each n = 3), and RT‐PCR (each n = 8). Results In multivariate logistic regression analyses, the fourth quartile (Q4) of serum resistin and G–A homozygotes were both associated with the latent sarcopenic obesity index defined by a visceral fat area of ≥ 100 cm2 and grip strength Q1 after adjustment for age and gender, with or without other confounding factors. RNA sequencing and pathway analysis showed that tumor necrosis factor (TNF) was involved in the top five pathways in the whole blood cells of G–A homozygotes compared with C–G homozygotes. RT‐PCR revealed that TNF mRNA was higher in G–A homozygotes than in C–G homozygotes. Conclusions The G–A haplotype was associated with the latent sarcopenic obesity index defined by grip strength in the Japanese cohort, could be mediated by TNF‐α.

This article investigates the chronomobility of international students in Australia going through COVID-19. Existing literature on international students approaches them largely in two manners: a market or victims. Using Shanti Robertson’s chronomobility, the study focuses on international students’ coping mechanisms and strategies for their next moves. Drawing from 15 in-depth interviews with international students formally enrolled in Australian institutions in Melbourne, the longest lockdown city during the pandemic, the authors find various ways of short-term coping mechanisms through meditation, physical exercises, virtual escapism and counselling. Furthermore, despite pandemic immobility, students presented a high level of resilience in making future decisions for post-pandemic mobilities. We conclude that family support and social networks are key to realise full potentials of international students as skilled migrants and valued members of society. Our manuscript contributes to the field of migration and mobility by enriching Robertson’s concept of chrono-mobility and adding the empirical case study from international students in Australia during the latest pandemic in 2020–2021.

Objectives: Human Leukocyte Antigen HLA-B*57:01 and B*58:01 are considered anti-HIV-1 protective alleles. HLA-B*57:01/58:01-restricted HIV-1 Gag TW10 (TSTLQEQIGW, Gag residues 240-249) epitope-specific CD8+ T cell responses that frequently select for a Gag escape mutation, T242N, with viral fitness cost are crucial for HIV-1 control. Although this finding has been observed in cohorts where HIV-1 subtype B or C predominates, the protective impact of HLA-B*57:01/58:01 has not been reported in Southeast Asian countries where HIV-1 CRF01_AE is the major circulating strain. Here, the effect of HLA-B*57:01/58:01 on CRF01_AE infection was investigated. Methods: The correlation of HLA-B*57:01/58:01 with viral load and CD4 counts were analyzed in the CRF01_AE-infected Vietnamese cohort (N = 280). The impact of the T242N mutation on CRF01_AE replication capacity was assessed. Results: HLA-B*57:01/58:01-positive individuals mostly had HIV-1 with T242N (62/63) but showed neither a significant reduction in viral load nor increased CD4 counts relative to B*57:01/58:01-negative participants. In vitro and in vivo analyses revealed a significant reduction in viral fitness of CRF01_AE with T242N. In silico analysis indicated reduced presentation of epitopes in the context of CRF01_AE compared to subtype B or C in 10/16 HLA-B*57:01/58:01-restricted HIV-1 epitopes. Conclusion: The protective impact of HLA-B*57:01/58:01 on CRF01_AE infection is impaired despite strong suppressive pressure by TW10-specific CD8+ T cells.

Summary: Modern Japanese people have two major ancestral populations: indigenous Jomon hunter-gatherers and continental East Asian farmers. To determine the formation process of the current Japanese population, we developed a detection method for variants derived from ancestral populations using a summary statistic, the ancestry marker index (AMI). We applied AMI to modern Japanese population samples and identified 208,648 single nucleotide polymorphisms (SNPs) that were likely derived from the Jomon people (Jomon-derived variants). Analysis of Jomon-derived variants in 10,842 modern Japanese individuals recruited from all over Japan revealed that the admixture proportions of the Jomon people varied between prefectures, probably owing to the prehistoric population size difference. The estimated allele frequencies of genome-wide SNPs in the ancestral populations of the modern Japanese suggested their adaptive phenotypic characteristics to their respective livelihoods. Based on our findings, we propose a formation model for the genotypic and phenotypic gradations of the current Japanese archipelago populations.

Abstract The outbreak of COVID-19 caused by infection with SARS-CoV-2 virus has become a worldwide pandemic, and the number of patients presenting with respiratory failure is rapidly increasing in Japan. An international meta-analysis has been conducted to identify genetic factors associated with the onset and severity of COVID-19, but these factors have yet to be fully clarified. Here, we carried out genomic analysis based on a genome-wide association study (GWAS) in Japanese COVID-19 patients to determine whether genetic factors reported to be associated with the onset or severity of COVID-19 in the international meta-GWAS are replicated in the Japanese population, and whether new genetic factors exist. Although no significant genome-wide association was detected in the Japanese GWAS, an integrated analysis with the international meta-GWAS identified for the first time the involvement of the IL17A/IL17F gene in the severity of COVID-19. Among nine genes reported in the international meta-GWAS as genes involved in the onset of COVID-19, the association of FOXP4-AS1, ABO, and IFNAR2 genes was replicated in the Japanese population. Moreover, combined analysis of ABO and FUT2 genotypes revealed that the presence of oral AB antigens was significantly associated with the onset of COVID-19. FOXP4-AS1 and IFNAR2 were also significantly associated in the integrated analysis of the Japanese GWAS and international meta-GWAS when compared with severe COVID-19 cases and the general population. This made it clear that these two genes were also involved in not only the onset but also the severity of COVID-19. In particular, FOXP4-AS1 was not found to be associated with the severity of COVID-19 in the international meta-GWAS, but an integrated analysis with the Japanese GWAS revealed an association with severity. Individuals with the SNP risk allele found between IL17A and IL17F had significantly lower mRNA expression levels of IL17F, suggesting that activation of the innate immune response by IL17F may play an important role in the severity of SARS-CoV-2 infection.

Abstract Background Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors). Methods In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors. Results The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors. Conclusions Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Abstract Hepatitis B (HB) vaccines (Heptavax-II and Bimmugen) designed based on HBV genotypes A and C are mainly used for vaccination against HB in Japan. To determine whether there are differences in the genetic background associated with vaccine responsiveness, genome-wide association studies were performed on 555 Heptavax-II and 1193 Bimmugen recipients. Further HLA imputation and detailed analysis of the association with HLA genes showed that two haplotypes, DRB1*13:02-DQB1*06:04 and DRB1*04:05-DQB1*04:01, were significantly associated in comparison with high-responders (HBsAb > 100 mIU/mL) for the two HB vaccines. In particular, HLA-DRB1*13:02-DQB1*06:04 haplotype is of great interest in the sense that it could only be detected by direct analysis of the high-responders in vaccination with Heptavax-II or Bimmugen. Compared with healthy controls, DRB1*13:02-DQB1*06:04 was significantly less frequent in high-responders when vaccinated with Heptavax-II, indicating that high antibody titers were less likely to be obtained with Heptavax-II. As Bimmugen and Heptavax-II tended to have high and low vaccine responses to DRB1*13:02, 15 residues were found in the Heptavax-II-derived antigenic peptide predicted to have the most unstable HLA-peptide binding. Further functional analysis of selected hepatitis B patients with HLA haplotypes identified in this study is expected to lead to an understanding of the mechanisms underlying liver disease.

The development of liver cancer in patients with hepatitis B is a major problem, and several models have been reported to predict the development of liver cancer. However, no predictive model involving human genetic factors has been reported to date. For the items incorporated in the prediction model reported so far, we selected items that were significant in predicting liver carcinogenesis in Japanese patients with hepatitis B and constructed a prediction model of liver carcinogenesis by the Cox proportional hazard model with the addition of Human Leukocyte Antigen (HLA) genotypes. The model, which included four items—sex, age at the time of examination, alpha-fetoprotein level (log10AFP) and presence or absence of HLA-A*33:03—revealed an area under the receiver operating characteristic curve (AUROC) of 0.862 for HCC prediction within 1 year and an AUROC of 0.863 within 3 years. A 1000 repeated validation test resulted in a C-index of 0.75 or higher, or sensitivity of 0.70 or higher, indicating that this predictive model can distinguish those at high risk of developing liver cancer within a few years with high accuracy. The prediction model constructed in this study, which can distinguish between chronic hepatitis B patients who develop hepatocellular carcinoma (HCC) early and those who develop HCC late or not, is clinically meaningful.

Abstract Background Dengue patients develop different disease severity ranging from mild (dengue fever [DF]) to severe forms (dengue hemorrhagic fever [DHF] and the fatal dengue shock syndrome [DSS]). Host genetics are considered to be one factor responsible for the severity of dengue outcomes. To identify genes associated with dengue severity that have not been studied yet, we performed genetic association analyses of interferon lambda 3 (IFNL3), CD27, and human leukocyte antigen-DPB1 (HLA-DPB1) genes in Thai dengue patients. Methods A case–control association study was performed in 877 children (age ≤ 15 years) with dengue infection (DF, n = 386; DHF, n = 416; DSS, n = 75). A candidate single nucleotide polymorphism of each of IFNL3, CD27, and HLA-DPB1 was selected to be analyzed. Genotyping was performed by TaqMan real-time PCR assay, and the association with dengue severity was examined. Results The rs9277534 variant of HLA-DPB1 was weakly associated with DHF. The genotype GG and G allele conferred protection against DHF (p = 0.04, odds ratio 0.74 for GG genotype, p = 0.03, odds ratio 0.79 for G allele). The association became borderline significant after adjusting for confounders (p = 0.05, odds ratio 0.82). No association was detected for IFNL3 or CD27. Conclusions The present study demonstrated the weak association of the rs9277534 variant of HLA-DPB1 with protection against DHF. This variant is in the 3′ untranslated region and affects HLA-DPB1 surface protein expression. Our finding suggests that HLA-DPB1 may be involved in DHF pathogenesis.

Objectives: Patients with dengue exhibit a range of symptoms from an acute febrile illness (dengue fever, DF), to dengue hemorrhagic fever (DHF), and to the most severe outcome, dengue shock syndrome (DSS). This study was performed to determine the host genetic factors responsible for dengue severity. Two single nucleotide polymorphisms (SNPs) of the interferon lambda 1 (IFNL1) gene (rs30461 and rs7247086) were analyzed for their association with dengue severity in a Thai population. Methods: This was a case–control association study involving 877 patients under the age of 15 years (DF, n = 386; DHF, n = 416; DSS, n = 75). Genotyping was performed by TaqMan real-time PCR assay. Results: It was found that the rs7247086 variant of IFNL1 was associated with DHF, but not DSS. Genotypes CT and TT and the T allele were protective against DHF (p = 0.03, odds ratio 0.62 for CT, odds ratio 0.13 for TT; and p = 0.01, odds ratio 0.54 for the T allele). The other SNP tested was not associated with DHF or DSS. Conclusions: The rs7247086 variant of IFNL1 (the T allele) was found to be protective against DHF, suggesting that IFNL1 may play a role in the pathogenesis of DHF. Keywords: Dengue, IFNL1, Genetic association, Disease severity

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.

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.

Abstract Aim/Introduction Insulin administration was found to trigger type 1 diabetes in six Japanese type 2 diabetes patients with type 1 diabetes high‐risk human leukocyte antigen class II and the class I allele of the insulin gene variable number tandem repeat genotype. The objective of the present study was to assess the contribution of non‐human leukocyte antigen single‐nucleotide polymorphisms (SNPs) to the risk of developing insulin‐triggered type 1 diabetes. Materials and Methods We genotyped 13 type 1 diabetes susceptible SNPs in six patients and compared them with those in Japanese controls (Hap Map3‐JPT). The SNPs that showed statistically significant results were further analyzed using non‐diabetic control participants and participants with type 2 diabetes at the Ehime University Hospital. Results The risk allele frequency of BACH2 rs3757247 in the six patients was significantly more frequent than that in 86 Japanese controls (P = 0.038). No significant difference in the allele frequency was observed in the other SNPs. This result was confirmed by the findings that the risk allele frequency of BACH2 in the six patients was significantly higher than that in the non‐diabetic control participants (n = 179) and type 2 diabetes with or without insulin treatment (n = 154 or n = 152; P = 0.035, 0.034 or 0.037, respectively). Despite being statistically not significant, the six patients were all homozygous for the CLEC16A rs12708716 risk allele and five were homozygous for the CLEC16A rs2903692 risk allele. Conclusions In addition to type 1 diabetes high‐risk human leukocyte antigen class II and the class I allele of the insulin gene variable number tandem repeat genotype, the possibility that the risk variants of BACH2 and CLEC16A could contribute to the development of insulin‐triggered type 1 diabetes cannot be excluded.

Abstract Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical Mycobacterium tuberculosis (Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring katG-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, PPE18/19, gid, emrB, Rv1588c, and pncA, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to PPE46/47-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and 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.

Background: A missense variant (rs373863828:G > A; p.Arg457Gln) of the CREBRF gene is strongly associated with a higher body mass index (BMI; kg/m2) in Polynesian populations. This variant has also been reported to be associated with lower total cholesterol in Samoans. Aim: The aim of this study is to examine the association of rs373863828:G > A with levels of serum lipids in four Pacific populations. Methods: A total of 613 adult subjects were recruited from Tonga (Polynesians) and the Solomon Islands (Melanesians and Micronesians). Multiple regression analyses adjusted for age and sex were performed to examine the association of rs373863828 with levels of serum lipids in each population. Results: A significant association of rs373863828:G > A with lower level of HDL-cholesterol was detected in the Tonga population (β = –3.32 and p-value = 0.030). The expected change in HDL-cholesterol with respect to a single copy of the rs373863828-A allele was 3.32 mg/dL. However, the association between rs373863828-A and lower levels of HDL-cholesterol was not significant after further adjustment for BMI in the Tonga population (β = –2.32 and p-value = 0.13). Conclusions: The rs373863828-A allele may not directly affect the level of serum HDL-cholesterol independent of BMI. To confirm the present findings, association studies with large sample sizes and functional analyses are required.

Because ≈90% of malaria cases occur in Africa, emergence of artemisinin-resistant Plasmodium falciparum in Africa poses a serious public health threat. To assess emergence of artemisinin-resistant parasites in Uganda during 2014–2016, we used the recently developed ex vivo ring-stage survival assay, which estimates ring-stage–specific P. falciparum susceptibility to artemisinin. We conducted 4 cross-sectional surveys to assess artemisinin sensitivity in Gulu, Uganda. Among 194 isolates, survival rates (ratio of viable drug-exposed parasites to drug-nonexposed controls) were high (>10%) for 4 isolates. Similar rates have been closely associated with delayed parasite clearance after drug treatment and are considered to be a proxy for the artemisinin-resistant phenotype. Of these, the PfKelch13 mutation was observed in only 1 isolate, A675V. Population genetics analysis suggested that these possibly artemisinin-resistant isolates originated in Africa. Large-scale surveillance of possibly artemisinin-resistant parasites in Africa would provide useful information about treatment outcomes and help regional malaria control.

Objective: To know whether the effect of interferon-induced protein with tetratricopeptide repeats (IFIT) 1 polymorphism influences the susceptibility of cerebral malaria outcome. Methods: Case-control association study was performed among 314 Thai patients (110 with cerebral malaria and 204 with uncomplicated malaria) infected with Plasmodium falciparum. Genotyping for five tag-single nucleotide polymorphisms of IFIT1 was performed by endpoint genotyping. Results: Genotype frequencies of all tag-SNPs (single nucleotide polymorphisms) showed no association with malaria outcome. However, C allele of rs11203109 was associated with the protection from cerebral malaria (OR=0.62, 95% CI=0.38-0.99, P=0.048). Two single nucleotide polymorphisms (rs5786868 and rs57941432) were in linkage disequilibrium with rs11203109. Conclusions: This suggests that our associated single nucleotide polymorphism (rs11203109) might be a genetic marker of cerebral malaria progression in the Thai population.

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

This paper, drawing on Pratt’s notion of ‘transculturation’ and Bhabha’s ‘third space’, presents an example of language learning tasks that empower learners’ agency and promote their cross-cultural awareness and sensitivities to a different set of cultural expectations, using a naturally occurred Japanese thanking episodes. The paper discusses the merits of Natural Conversation Reconstruction Tasks (NCRTs) as a practical method for helping L2 learners develop this ‘intercultural competence’. It is based on a qualitative study of the results of one NCRT created for use in the context of teaching Japanese as a L2 in a multicultural society. It suggests the NCRT encourages the learners to explore the intersection where language use, speaker intention and L1 and L2 cultural norms meet. Such a process helps the learners become aware of socially expected patterns of communication in L1 and L2 in terms of the choices of speech act, formulaic expressions, sequential organization and politeness orientation. The learners’ comments suggest that the NCRT helps learners transcend their cultural boundaries by overcoming their narrow understanding of ‘thanking’ as ‘expressions of gratitude and appreciation’ and by cross-culturally widening their views of what counts as thanking. The NCRT with rich contextual information promotes the learners’ intercultural awareness, sensitivity to context and intercultural exploration in the space between L1 and L2, where they have authority and freedom of making sense of conversations, and pragmatics is fully integrated into language pedagogy.

This paper, drawing on Pratt’s notion of ‘transculturation’ and Bhabha’s ‘third space’, presents an example of language learning tasks that empower learners’ agency and promote their cross-cultural awareness and sensitivities to a different set of cultural expectations, using a naturally occurred Japanese thanking episodes. The paper discusses the merits of Natural Conversation Reconstruction Tasks (NCRTs) as a practical method for helping L2 learners develop this ‘intercultural competence’. It is based on a qualitative study of the results of one NCRT created for use in the context of teaching Japanese as a L2 in a multicultural society. It suggests the NCRT encourages the learners to explore the intersection where language use, speaker intention and L1 and L2 cultural norms meet. Such a process helps the learners become aware of socially expected patterns of communication in L1 and L2 in terms of the choices of speech act, formulaic expressions, sequential organization and politeness orientation. The learners’ comments suggest that the NCRT helps learners transcend their cultural boundaries by overcoming their narrow understanding of ‘thanking’ as ‘expressions of gratitude and appreciation’ and by cross-culturally widening their views of what counts as thanking. The NCRT with rich contextual information promotes the learners’ intercultural awareness, sensitivity to context and intercultural exploration in the space between L1 and L2, where they have authority and freedom of making sense of conversations, and pragmatics is fully integrated into language pedagogy.

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.