A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis.

BACKGROUND AIMS: Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-non-specific) and non-shared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC. APPROACH RESULTS: Protein tyrosine phosphatase non-receptor type 2 (PTPN2) was identified as a novel PBC susceptibility gene locus through a GWAS and subsequent genome-wide meta-analysis involving 2,181 cases and 2,699 controls from the Japanese population (GWAS-lead variant: rs8098858, p=2.6×10-8). In-silico and in-vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells (Tfh) and plasmacytoid dendritic cells (pDCs). Infiltration of PTPN2-positive T-cells and pDCs were confirmed in the portal area of the PBC-liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in-vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758. CONCLUSIONS: PTPN2, a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC via an insufficient negative feedback loop caused by the PTPN2 risk allele of rs2292758 in IFN signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases.

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Dried blood spot-based host genome analysis technique targeting pathological associations with hepatitis B: Development and clinical application in the Cambodian population.

AIM: Reports of patients with hepatitis B have highlighted associations between polymorphisms in the human leukocyte antigen (HLA)-DPB1, CXCL13, and CXCR5 genes and disease pathology. Owing to its potential to contribute to the development of new diagnostic and therapeutic methods, we aimed to establish a reliable host genome analysis technique that can be used in countries with inadequate infrastructure. METHOD: We compared multiple commercially available kits for dried blood spot (DBS)-based sample collection to develop a basic DBS-based host genome analysis technique. We then collected blood samples from Cambodian patients with hepatitis  B and performed single-nucleotide polymorphism genotyping and HLA allele typing by the DBS system. RESULT: We were able to perform single-nucleotide polymorphism genotyping and HLA allele typing with host DNA samples obtained using a combination of a HemaSpot™ filter paper-based device and a SMITEST® EX-R&D DNA extraction kit. The accuracy of genotyping using samples obtained by this method was not inferior to one using samples obtained by venipuncture. In the Cambodian population, significant associations of HLA-DPB1*04:01 with protection against chronic hepatitis B virus (HBV) infection, and HLA-DPB1*05:01 and HLA-DPB1*13:01 with susceptibility to chronic HBV infection were identified. CONCLUSION: Based on the DBS system, we clarified the associations of HLA-DPB1 alleles with chronic HBV infection in the Cambodian population for the first time. Because the DBS is a low-cost, durable, transportable, and easy-to-handle modality, genetic analysis based on the DBS system is a feasible strategy for obtaining a deeper understanding of HBV epidemiology, especially in middle- or low-income countries.

Novel HLA allele associations with susceptibility, staging, symptomatic state, autoimmune hepatitis and hepatocellular carcinoma events for primary biliary cholangitis in the Japanese population.

Primary biliary cholangitis (PBC) is a rare autoimmune disease with a clear predisposition for human leukocyte antigen (HLA)-DR/DQ-associated loss of immune tolerance for the E2 component of the pyruvate dehydrogenase complex. Three-field-resolution HLA imputation of 1,670 Japanese PBC patients and 2,328 healthy controls was conducted using Japanese population-specific HLA reference panels. Eighteen previously reported Japanese PBC-associated HLA alleles were confirmed and extended to 3-field-resolution, including HLA-DRB1*08:03 to HLA-DRB1*08:03:02, HLA-DQB1*03:01 to HLA-DQB1*03:01:01, HLA-DQB1*04:01 to HLA-DQB1*04:01:01 and HLA-DQB1*06:04 to HLA-DQB1*06:04:01. In addition, additional significant novel HLA alleles were identified, including 3 novel susceptible HLA-DQA1 alleles: HLA-DQA1*03:03:01, HLA-DQA1*04:01:01, HLA-DQA1*01:04:01 and 1 novel protective HLA-DQA1 allele, HLA-DQA1*05:05:01. In addition, PBC patients carrying HLA-DRB1*15:01:01 and HLA-DQA1*03:03:01 would have a higher predisposition toward developing concomitant autoimmune hepatitis (AIH). Further, late-stage and symptomatic PBC shared the same susceptible HLA alleles of HLA-A*26:01:01, HLA-DRB1*09:01:02 and HLA-DQB1*03:03:02. Lastly, HLA-DPB1*05:01:01 was identified as a potential risk HLA allele for development of hepatocellular carcinoma (HCC) in PBC patients. In conclusion, we have extended the current knowledge of HLA allele associations to 3-field resolution and identified novel HLA allele associations with predisposition risk, staging, symptomatic state, and AIH and HCC events for Japanese PBC patients.

Development of a novel microarray data analysis tool without normalization for genotyping degraded forensic DNA.

Since the arrest of the Golden State Killer in the US in April 2018, forensic geneticists have been increasingly interested in the investigative genetic genealogy (IGG) method. While this method has already been in practical use as a powerful tool for criminal investigation, we have yet to know well the limitations and potential risks. In this current study, we performed an evaluation study focusing on degraded DNA using the Affymetrix Genome-Wide Human SNP Array 6.0 platform (Thermo Fisher Scientific). We revealed one of the potential problems that occur during SNP genotype determination using a microarray-based platform. Our analysis results indicated that the SNP profiles derived from degraded DNA contained many false heterozygous SNPs. In addition, it was confirmed that the total amount of probe signal intensity on microarray chips derived from degraded DNA decreased significantly. Because the conventional analysis algorithm performs normalization during genotype determination, we concluded that noise signals could be genotype-called. To address this issue, we proposed a novel microarray data analysis method without normalization (nMAP). Although the nMAP algorithm resulted in a low call rate, it substantially improved genotyping accuracy. Finally, we confirmed the usefulness of the nMAP algorithm for kinship inferences. These findings and the nMAP algorithm will make a contribution to the advance of the IGG method.

Regulation of MEK inhibitor selumetinib sensitivity by AKT phosphorylation in the novel BRAF L525R mutant.

BACKGROUND: Oncogenic mutations in BRAF genes are found in approximately 5-10% of colorectal cancers. The majority of BRAF mutations are located within exons 11-15 of the catalytic kinase domains, with BRAF V600E accounting for more than 80% of the observed BRAF mutations. Sensitivity to BRAF- and mitogen-activated protein kinase (MEK) inhibitors varies depending on BRAF mutations and tumor cell types. Previously, we newly identified, BRAF L525R-mutation, in the activation segment of the kinase in colorectal cancer patient. Here, we characterized the function of the BRAF L525R mutation. METHODS: HEK293 cells harboring a BRAF mutation (V600E or L525R) were first characterized and then treated with cetuximab, dabrafenib, and selumetinib. Cell viability was measured using WST-1 assay and the expression of proteins involved in the extracellular signal-regulated kinase (ERK) and protein kinase B (AKT) signaling pathways was evaluated using western blot analysis. RESULTS: The MEK inhibitor selumetinib effectively inhibited cell proliferation and ERK phosphorylation in BRAF L525R cells but not in BRAF V600E cells. Further studies revealed that AKT phosphorylation was reduced by selumetinib in BRAF L525R cells but not in BRAF V600E cells or selumetinib-resistant BRAF L525R cells. Moreover, the AKT inhibitor overcame the selumetinib resistance. CONCLUSIONS: We established a model system harboring BRAF L525R using HEK293 cells. BRAF L525R constitutively activated ERK. AKT phosphorylation caused sensitivity and resistance to selumetinib. Our results suggest that a comprehensive network analysis may provide insights to identify effective therapies.

A common variant rs2054564 in ADAMST17 is associated with susceptibility to lumbar spondylosis.

The molecular pathophysiology underlying lumbar spondylosis development remains unclear. To identify genetic factors associated with lumbar spondylosis, we conducted a genome-wide association study using 83 severe lumbar spondylosis cases and 182 healthy controls and identified 65 candidate disease-associated single nucleotide polymorphisms (SNPs). Replication analysis in 510 case and 911 control subjects from five independent Japanese cohorts identified rs2054564, located in intron 7 of ADAMTS17, as a disease-associated SNP with a genome-wide significance threshold (P = 1.17 × 10-11, odds ratio = 1.92). This association was significant even after adjustment of age, sex, and body mass index (P = 7.52 × 10-11). A replication study in a Korean cohort, including 123 case and 319 control subjects, also verified the significant association of this SNP with severe lumbar spondylosis. Immunohistochemistry revealed that fibrillin-1 (FBN1) and ADAMTS17 were co-expressed in the annulus fibrosus of intervertebral discs (IVDs). ADAMTS17 overexpression in MG63 cells promoted extracellular microfibrils biogenesis, suggesting the potential role of ADAMTS17 in IVD function through interaction with fibrillin fibers. Finally, we provided evidence of FBN1 involvement in IVD function by showing that lumbar IVDs in patients with Marfan syndrome, caused by heterozygous FBN1 gene mutation, were significantly more degenerated. We identified a common SNP variant, located in ADAMTS17, associated with susceptibility to lumbar spondylosis and demonstrated the potential role of the ADAMTS17-fibrillin network in IVDs in lumbar spondylosis development.

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.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms.

BACKGROUND: Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. RESULTS: Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms. CONCLUSION: The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Practical forensic use of kinship determination using high-density SNP profiling based on a microarray platform, focusing on low-quantity DNA.

Instead of traditional short tandem repeat (STR) profiling, the genetic genealogy method, which uses hundreds of thousands of single nucleotide polymorphisms (SNPs) spread across genome-wide, has emerged as a powerful kinship determination tool and recently attracted great attention in forensic genetics. In this study, we explored the tolerance and viability of kinship discrimination based on a high-density SNP profile for forensic DNA, especially focusing on low-quantity DNA. Using the Affymetrix Genome-Wide Human SNP Array 6.0 platform (Thermo Fisher Scientific), the influence of low-quantity DNA on SNP genotype determination was evaluated. The low-quantity DNA samples failed once every few samples, the generated SNP profile had low data quality. Our investigation revealed that the SNP profile with low data quality contained many genotyping errors in which the SNP genotype changed from homozygote to heterozygote. The kinship discrimination analysis using KING software was directly influenced by these genotyping errors, which was confirmed that some unrelated pairs were mis-specified as 4th-degree relatives. We confirmed that the false heterozygous SNPs resulted in an inflation of kinship coefficient and a decrease of non-shared allele between a tested pair. To eliminate the influence of these genotyping errors and acquire an accurate kinship discrimination result, we developed a novel method to select only the robust SNPs, which stably give the genotype determination with high accuracy even in SNP profiles with low data quality. The application of our novel method led to the improved results of kinship discrimination up to the same level as in the SNP profile with high data quality. In addition, this study demonstrated the advantage of kinship analysis using a high-density SNP profile in the forensic field. It is well known that likelihood ratio calculation based on autosomal STR profile, which is the most commonly applied approach, has difficulty in gaining true kinship analysis results, especially when the relationship between the tested two individuals is more biologically distant. We showed the kinship discrimination analysis with a high-density SNP profile is more suitable for the case without close relatives, using the real case data. Although further study with larger samples will be necessary, this study indicated that practical forensic use of kinship determination with a high-density SNP profile would bring benefits to the forensic field.

Genome-Wide Association Study for Chronic Hepatitis B Infection in the Thai Population.

To identify novel host genetic variants that predispose to hepatitis B virus (HBV) persistence, we performed the first genome-wide association study in the Thai population involving 318 cases of chronic hepatitis B and 309 healthy controls after quality control measures. We detected the genome-wide significant association of the HLA class II region (HLA-DPA1/DPB1, rs7770370, p-value = 7.71 × 10-10, OR = 0.49) with HBV chronicity. Subsequent HLA allele imputation revealed HLA-DPA1*01:03 (Pc = 1.21 × 10-6, OR = 0.53), HLA-DPB1*02:01 (Pc = 2.17 × 10-3, OR = 0.50), and HLA-DQB1*06:09 (Pc = 2.17 × 10-2, OR = 0.07) as protective alleles, and HLA-DPA1*02:02 (Pc = 6.32 × 10-5, OR = 1.63), HLA-DPB1*05:01 (Pc = 1.13 × 10-4, OR = 1.72), HLA-DPB1*13:01 (Pc = 4.68 × 10-2, OR = 1.60), and HLA-DQB1*03:03 (Pc = 1.11 × 10-3, OR = 1.84) as risk alleles for HBV persistence. We also detected suggestive associations in the PLSCR1 (rs35766154), PDLIM5 (rs62321986), SGPL1 (rs144998273), and MGST1 (rs1828682) loci. Among single-nucleotide polymorphisms in the PLSCR1 locus, rs1061307 was identified as the primary functional variant by in silico/in vitro functional analysis. In addition to replicating the association of the HLA class II region, we detected novel candidate loci that provide new insights into the pathophysiology of chronic hepatitis B.

A single-nucleotide-polymorphism in the 5'-flanking region of MSX1 gene as a predictive marker candidate for platinum-based therapy of esophageal carcinoma.

BACKGROUND: Platinum derivatives are important treatment options for patients with esophageal carcinoma (EC), and a predictive marker for platinum-based therapy is needed for precision medicine. PATIENTS AND METHODS: This study contained two cohorts consisting of EC patients treated using platinum-based chemoradiation therapy (CRT) as the first-line and another external cohort of nationwide clinicogenomic data from the BioBank Japan (BBJ). RESULTS: Genome-wide association study (GWAS) of therapeutic outcomes, refractory disease or not, following platinum-based CRT as first-line in 94 patients in the first cohort suggested the association of 89 SNPs using p < 0.0001. The top 10 SNPs selected from each chromosomal region by odds ratio were evaluated for progression-free survival (PFS) and overall survival (OS) hazard ratios in the first cohort, resulting in four candidates (p < 0.0025). The four selected candidates were re-evaluated in another cohort of 24 EC patients, which included patients prospectively enrolled in this study to fulfill the sample size statistically suggested by the results of the first cohort, and of the four, only rs3815544 was replicated (p < 0.0125). Furthermore, this candidate genotype of rs3815544 proceeded to the re-evaluation study in an external cohort consisting of EC patients treated with platinum derivatives and/or by radiation therapy as the first-line treatment in BBJ, which confirmed that the alternative allele (G) of rs3815544 was statistically associated with non-response (SD or PD) to platinum-based therapy in EC patients (odds ratio = 1.801, p = 0.048). The methylation QTL database as well as online clinicogenomic databases suggested that the region including rs3815544 may regulate MSX1 expression through CpG methylation, and this down-regulation was statistically associated with poor prognosis after platinum-based therapies for EC. CONCLUSION: rs3815544 is a novel candidate predictive marker for platinum-based EC therapy.

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.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

Angiotensin-Converting Enzyme (ACE) 1 Gene Polymorphism and Phenotypic Expression of COVID-19 Symptoms.

The renin-angiotensin-aldosterone system (RAAS) appears to play an important role in SARS-CoV-2 infection. Polymorphisms within the genes that control this enzymatic system are candidates for elucidating the pathogenesis of COVID-19, since COVID-19 is not only a pulmonary disease but also affects many organs and systems throughout the body in multiple ways. Most striking is the fact that ACE2, one of the major components of the RAAS, is a prerequisite for SARS-COV-2 infection. Recently, we and other groups reported an association between a polymorphism of the ACE1 gene (a homolog of ACE2) and the phenotypic expression of COVID-19, particularly in its severity. The ethnic difference in ACE1 insertion (I)/deletion (D) polymorphism seems to explain the apparent difference in mortality between the West and East Asia. The purpose of this review was to further evaluate the evidence linking ACE1 polymorphisms to COVID-19. We searched the Medline database (2019-2021) for reference citations of relevant articles and selected studies on the clinical outcome of COVID-19 related to ACE1 I/D polymorphism. Although the numbers of patients are not large enough yet, most available evidence supports the notion that the DD genotype adversely influences COVID-19 symptoms. Surprisingly, small studies conducted in several countries yielded opposite results, suggesting that the ACE1 II genotype is a risk factor. This contradictory result may be the case in certain geographic areas, especially in subgroups of patients. It may also be due to interactions with other genes or to yet unexplained biochemical mechanisms. According to our hypothesis, such candidates are genes that are functionally involved in the pathophysiology of COVID-19, can act in concert with the ACE1 DD genotype, and that show differences in their frequency between the West and East Asia. For this, we conducted research focusing on Alu-related genes. The current study on the ACE1 genotype will provide potentially new clues to the pathogenesis, treatment, and diagnosis of SARS-CoV-2 infections.

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.

Exploring correlations in genetic and cultural variation across language families in northeast Asia.

Culture evolves in ways that are analogous to, but distinct from, genomes. Previous studies examined similarities between cultural variation and genetic variation (population history) at small scales within language families, but few studies have empirically investigated these parallels across language families using diverse cultural data. We report an analysis comparing culture and genomes from in and around northeast Asia spanning 11 language families. We extract and summarize the variation in language (grammar, phonology, lexicon), music (song structure, performance style), and genomes (genome-wide SNPs) and test for correlations. We find that grammatical structure correlates with population history (genetic history). Recent contact and shared descent fail to explain the signal, suggesting relationships that arose before the formation of current families. Our results suggest that grammar might be a cultural indicator of population history while also demonstrating differences among cultural and genetic relationships that highlight the complex nature of human history.

Genome-wide copy number variation analysis of hepatitis B infection in a Japanese population.

Genome-wide association studies have been performed to identify common genetic variants associated with hepatitis B (HB). However, little is known about copy number variations (CNVs) in HB. In this study, we performed a genome-wide CNV analysis between 1830 healthy controls and 1031 patients with HB infection after quality control. Using signal calling by the Axiom Analysis Suite and CNV detection by PennCNV software, we obtained a total of 4494 CNVs across all individuals. The genes with CNVs that were found only in the HB patients were associated with the immune system, such as antigen processing. A gene-level CNV association test revealed statistically significant CNVs in the contactin 6 (CNTN6) gene. Moreover, we also performed gene-level CNV association tests in disease subgroups, including hepatocellular carcinoma patients, liver cirrhosis patients, and HBV carriers, including asymptomatic carriers and patients with HBV-derived chronic hepatitis. Our findings from germline cells suggested that patient-specific CNVs may be inherent genetic risk factors for HB.