Identification of the hybrid gene LILRB5-3 by long-read sequencing and implication of its novel signaling function.

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.

Distinct Clinical Effects of Two RP1L1 Hotspots in East Asian Patients With Occult Macular Dystrophy (Miyake Disease): EAOMD Report 4.

Purpose: To characterize the clinical effects of two RP1L1 hotspots in patients with East Asian occult macular dystrophy (OMD). Methods: Fifty-one patients diagnosed with OMD harboring monoallelic pathogenic RP1L1 variants (Miyake disease) from Japan, South Korea, and China were enrolled. Patients were classified into two genotype groups: group A, p.R45W, and group B, missense variants located between amino acids (aa) 1196 and 1201. The clinical parameters of the two genotypes were compared, and deep learning based on spectral-domain optical coherence tomographic (SD-OCT) images was used to distinguish the morphologic differences. Results: Groups A and B included 29 and 22 patients, respectively. The median age of onset in groups A and B was 14.0 and 40.0 years, respectively. The median logMAR visual acuity of groups A and B was 0.70 and 0.51, respectively, and the survival curve analysis revealed a 15-year difference in vision loss (logMAR 0.22). A statistically significant difference was observed in the visual field classification, but no significant difference was found in the multifocal electroretinographic classification. High accuracy (75.4%) was achieved in classifying genotype groups based on SD-OCT images using machine learning. Conclusions: Distinct clinical severities and morphologic phenotypes supported by artificial intelligence-based classification were derived from the two investigated RP1L1 hotspots: a more severe phenotype (p.R45W) and a milder phenotype (1196-1201 aa). This newly identified genotype-phenotype association will be valuable for medical care and the design of therapeutic trials.

Chronomobility of international students under COVID-19 Australia.

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.

Positive association between serum resistin and smoking was strongest in homozygotes of the G-A haplotype at c.-420 C>G and c.-358 G>A in RETN promoter: the Toon Genome Study.

Resistin is mainly expressed in human monocytes/macrophages and is associated with insulin resistance, inflammation, and atherosclerosis. Serum resistin is strongly correlated with the G-A haplotype defined by single nucleotide polymorphisms (SNPs) c.-420 C>G (SNP-420) (rs1862513) and c.-358 G>A (SNP-358) (rs3219175) in the promoter region of the human resistin gene (RETN). Smoking is also associated with insulin resistance. We investigated the association between smoking and serum resistin and the effect of the G-A haplotype on this association. Participants were recruited under the Toon Genome Study (an observational epidemiology research in the Japanese population). Of these, 1975 subjects genotyped for both SNP-420 and SNP-358 were analyzed for serum resistin by grouping them based on smoking status and G-A haplotype status. RETN mRNA, isolated from whole blood cells, was evaluated in smokers (n = 7) and age-, sex-, and BMI-matched non-smokers (n = 7) with the G-A haplotype homozygotes. Serum resistin tended to be higher in current smokers who smoked more cigarettes per day (P for trend < 0.0001). The positive association between serum resistin and smoking was strongest in the G-A haplotype homozygotes, followed by heterozygotes and non-carriers (interaction P < 0.0001). This positive association was stronger in the G-A homozygotes than the C-G homozygotes (interaction P < 0.0001). RETN mRNA was 1.40-fold higher in smokers than non-smokers with the G-A homozygotes (P = 0.022). Therefore, the positive association between serum resistin and smoking was strongest in the G-A haplotype homozygotes defined by RETN SNP-420 and SNP-358.

Modern Japanese ancestry-derived variants reveal the formation process of the current Japanese regional gradations.

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.

Resistin G-A haplotype at SNP-420/-358 is associated with the latent sarcopenic obesity index in the toon genome study.

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-α.

Prediction Model with HLA-A*33:03 Reveals Number of Days to Develop Liver Cancer from Blood Test.

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.

Impaired protective role of HLA-B*57:01/58:01 in HIV-1 CRF01_AE infection: a cohort study in Vietnam.

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.

Super high-resolution single-molecule sequence-based typing of HLA class I alleles in HIV-1 infected individuals in Ghana.

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.

Genetic association of IL17 and the importance of ABO blood group antigens in saliva to COVID-19.

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.

Recurrent Hypoglycemia Due to a High Titer of Insulin Antibody in Response to Exogenous Insulin Administration in Two Cases of Type 1 Diabetes.

In the first case, a 60-year-old man who was using continuous subcutaneous insulin infusion (CSII), developed recurrent hypoglycemia due to insulin antibodies. This is the first report of such a case using CSII. In the second case, a 70-year-old man was follow-up case who developed hypoglycemia while using human insulin. In both cases, the hypoglycemia subsided after switching to multiple daily insulin injection and/or insulin preparation. The results of Scatchard analyses of the two cases were similar to those of cases of insulin autoimmune syndrome (IAS) that improved after recovery from hypoglycemia.The clinical characteristics and Scatchard analysis data were essentially the same as those for IAS, except for the presence of insulin administration.

Predictors associated with a better response to the Japanese aluminum-free hepatitis A vaccine, Aimmugen® , for people living with HIV.

AIM: After the hepatitis A virus (HAV) outbreak among men who have sex with men (MSM) around 2018, the importance of HAV vaccination was emphasized, especially for MSM-living with human immunodeficiency virus (MSM-LWHIV). Aimmugen® is licensed and distributed exclusively in Japan. While administration of three doses is recommended, 85% of recipients in the general population were reported to acquire seroprotection after the second dose. In this study, we evaluated the efficacy of two or three vaccine doses along with predictors associated with the response to Aimmugen® in MSM-LWHIV. METHODS: We retrospectively examined anti-HA-IgG titers of MSM-LWHIV vaccinated with Aimmugen® in our hospital. Patients' data were collected from medical records. RESULTS: Between January 2018 and October 2019, 141 subjects whose median age was 46 years old, were examined. All the subjects were on antiretroviral therapy (ART) and the median CD4 count was 615/µL. The acquisition rate of protectable anti-HA-IgG titers after the second and third dose was 71.1% and 98.6%, respectively. In 114 subjects whose anti-HA-IgG titers were tested after the second-dose, factors significantly associated with better response were prolonged ART duration and higher CD4 count. The titers of anti-HA-IgG after the third dose were higher in those who became seropositive after the second-dose than those who did not. CONCLUSIONS: Three-dose of Aimmugen® for MSM-LWHIV was effective while two-dose was less effective compared to non-HIV-infected people. People-LWHIV with shorter duration of ART and lesser CD4 cell count achieved lower titers of anti-HA-IgG and might require an additional vaccination.

Cytochrome b-245 Alpha Chain Gene Variants and Arterial Function in Indonesian Short Stature Children.

BACKGROUND: The association between short stature, undernutrition and the risk to cardiovascular disease has been clinically established. Genetic factor, particularly the variants in cytochrome b-245 alpha chain (CYBA) gene, which alter the formation of nicotinamide adenine dinucleotide phosphate (NADPH) oxidase level, might affect arterial function. This study aimed to observe the association between single nucleotide variants (SNVs) of the CYBA gene and the arterial function of short stature children to understand the reason why some people with short stature develop cardiovascular disease. METHODS: A total of 142 genomic deoxyribonucleic acid (DNA) samples have been collected from short stature children in Brebes, Central Java, Indonesia. Four common single-nucleotide polymorphisms (SNPs): C242T (rs4673), A640G (rs1049255), -930A>G (rs9932581) and *49A>G (rs7195830) in the CYBA gene were examined using TaqMan allelic discrimination assay. The arterial function was measured using transthoracic echocardiography and described as aortic stiffness and distensibility index. Statistical analysis was done to find a significant difference in arterial function between genotypes of each SNV. RESULTS: A P-value of < 0.05 was considered significant. In rs9932581 (-930A>G) of CYBA gene, the subjects with GG genotype were found to have significantly lower arterial stiffness and higher distensibility compared to AA and AG genotypes. No significant difference was found in the other SNVs. CONCLUSION: The GG genotype in rs9932581 of the CYBA gene might have a protective effect on cardiovascular disease in short stature children.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians.

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.

Geographic variation in the polygenic score of height in Japan.

A geographical gradient of height has existed in Japan for approximately 100 years. People in northern Japan tend to be taller than those in southern Japan. The differences in annual temperature and day length between the northern and southern prefectures of Japan have been suggested as possible causes of the height gradient. Although height is well known to be a polygenic trait with high heritability, the genetic contributions to the gradient have not yet been explored. Polygenic score (PS) is calculated by aggregating the effects of genetic variants identified by genome-wide association studies (GWASs) to predict the traits of individual subjects. Here, we calculated the PS of height for 10,840 Japanese individuals from all 47 prefectures in Japan. The median height PS for each prefecture was significantly correlated with the mean height of females and males obtained from another independent Japanese nation-wide height dataset, suggesting genetic contribution to the observed height gradient. We also found that individuals and prefectures genetically closer to continental East Asian ancestry tended to have a higher PS; modern Japanese people are considered to have originated as result of admixture between indigenous Jomon people and immigrants from continental East Asia. Another PS analysis based on the GWAS using only the mainland Japanese was conducted to evaluate the effect of population stratification on PS. The result also supported genetic contribution to height, and indicated that the PS might be affected by a bias due to population stratification even in a relatively homogenous population like Japanese.

Importance of HBsAg recognition by HLA molecules as revealed by responsiveness to different hepatitis B vaccines.

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.

Association of HLA-DRB1*09:01 with severe COVID-19.

HLA-A, -C, -B, and -DRB1 genotypes were analyzed in 178 Japanese COVID-19 patients to investigate the association of HLA with severe COVID-19. Analysis of 32 common HLA alleles at four loci revealed a significant association between HLA-DRB1*09:01 and severe COVID-19 (odds ratio [OR], 3.62; 95% CI, 1.57-8.35; p = 0.00251 [permutation p value = 0.0418]) when age, sex, and other common HLA alleles at the DRB1 locus were adjusted. The DRB1*09:01 allele was more significantly associated with risk for severe COVID-19 compared to preexisting medical conditions such as hypertension, diabetes, and cardiovascular diseases. These results indicate a potential role for HLA in predisposition to severe COVID-19.

Characterization of LILRB3 and LILRA6 allelic variants in the Japanese population.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) are encoded by members of a human multigene family, comprising 11 protein-coding genes and two pseudogenes. Among the LILRs, LILRB3 and LILRA6 show the highest homology with each other, along with high allelic and copy number variations. Therefore, it has been difficult to discriminate between them, both genetically and functionally, precluding disease association studies of LILRB3 and LILRA6. In this study, we carefully performed variant screening of LILRB3 and LILRA6 by cDNA cloning from Japanese individuals and identified four allelic lineages showing significantly high non-synonymous-to-synonymous ratios in pairwise comparisons. Furthermore, the extracellular domains of the LILRB3*JP6 and LILRA6*JP1 alleles were identical at the DNA level, suggesting that gene conversion-like events diversified LILRB3 and LILRA6. To determine the four allelic lineages from genomic DNA, we established a lineage typing method that accurately estimated the four allelic lineages in addition to specific common alleles from genomic DNA. Analysis of LILRA6 copy number variation revealed one, two, and three copies of LILRA6 in the Japanese-in-Tokyo (JPT) population. Flow cytometric analysis showed that an anti-LILRB3 antibody did not recognize the second most common lineage in the Japanese population, indicating significant amino acid differences across the allelic lineages. Taken together, our findings indicate that our lineage typing is useful for classifying the lineage-specific functions of LILRB3 and LILRA6, serving as the basis for disease association studies.

The number of SNPs required for distinguishing Japanese from other East Asians.

In some cases, it is necessary to estimate the national origin of an unknown subject in forensic medicine. The use of single nucleotide polymorphism (SNP) markers appears to be very effective for this purpose, since genome-wide SNP genotype data of many human populations are publicly available. In this study, we examined the number of SNPs that could objectively and accurately distinguish Japanese subjects (1KG-JPT) from the other East Asians (1KG-CDX, -CHB, -CHS, and -KHV) using the combination of principal component analysis and hierarchical cluster analysis. A computer simulation showed that approximately 3000 randomly selected SNPs were enough for the discrimination. Our results suggest that at least a 0.024% coverage is needed in the next generation sequencing experiment to objectively determine whether an unknown person is Japanese or not if the amount of DNA sample from him/her is insufficient or the quality is low.

Prefecture-level population structure of the Japanese based on SNP genotypes of 11,069 individuals.

We analyzed genome-wide single-nucleotide polymorphism data of 11,069 Japanese individuals recruited from all 47 prefectures of Japan to clarify their genetic structure. The principal component analysis at the prefectural level enabled us to study the relationship between geographical location and genetic differentiation. The results revealed that the mainland Japanese were not genetically homogeneous, and the genetic structure could be explained mainly by the degree of Jomon ancestry and the geographical location. One of the interesting findings was that individuals in the Shikoku region (i.e., Tokushima Prefecture, Kagawa Prefecture, Ehime Prefecture, and Kochi Prefecture) were genetically close to Han Chinese. Therefore, the genetic components of immigrants from continental East Asia in the Yayoi period may have been well maintained in Shikoku. The present results will be useful for understanding the peopling of Japan, and also provide suggestions for recruiting subjects in genetic association studies.