Deep whole-genome sequencing reveals recent selection signatures linked to evolution and disease risk of Japanese

A common deletion at BAK1 reduces enhancer activity and confers risk of intracranial germ cell tumors

Intracranial germ cell tumors (IGCTs) are rare brain neoplasms that mainly occur in children and adolescents with a particularly high incidence in East Asian populations. Here, we conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10⁻⁹, odds ratio = 2.46 [95% CI: 1.83–3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. In-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity, suggesting its contribution to IGCTs predisposition through altering BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10⁻⁴, Spearman’s ρ = 0.48). These results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.

Intracranial germ cell tumors (IGCTs) are rare brain tumors mainly diagnosed in children and young adults. Here, the authors conduct a genome-wide association study for IGCTs, identify a risk locus at BAK1, and characterize its functional consequences.

Genome-Wide Association Study of Intracranial Artery Stenosis Followed by Phenome-Wide Association Study

The genetic background of intracranial artery stenosis (ICAS), a major cause of ischemic stroke, remains elusive. We performed the world's first genome-wide association study (GWAS) of ICAS using DNA samples from Japanese subjects, to identify the genetic factors associated with ICAS and their correlation with clinical features. We also conducted a phenome-wide association study (PheWAS) of the top variant identified via GWAS to determine its association with systemic disease. The GWAS involved 408 patients with ICAS and 349 healthy controls and utilized an Asian Screening Array of venous blood samples. The PheWAS was performed using genotypic and phenotypic data of the Biobank Japan Project, which contained information on 46 diseases and 60 quantitative trait data from > 150,000 Japanese individuals. The GWAS revealed that the East Asian-specific functional variant of RNF213, rs112735431 (c.14429G > A, p.Arg4810Lys), was associated with ICAS (odds ratio, 12.3; 95% CI 5.5 to 27.5; P = 7.8 × 10^-10). Stratified analysis within ICAS cases demonstrated that clinical features of those with and without the risk allele were different. PheWAS indicated that high blood pressure and angina were significantly associated with RNF213 rs112735431. The first GWAS of ICAS, which stratifies subpopulations within the ICAS cases with distinct clinical features, revealed that RNF213 rs112735431 was the most significant variant associated with ICAS. Thus, RNF213 rs112735431 shows potential as an important clinical biomarker that characterizes pleiotropic risk in various vascular diseases, such as blood pressure and angina, thereby facilitating personalized medicine for systemic vascular diseases in East Asian populations.

ALDH2 variance in disease and populations

The ALDH2*2 missense variant that commonly causes alcohol flushing reactions is the single genetic polymorphism associated with the largest number of traits in humans. The dysfunctional ALDH2 variant affects nearly 8% of the world population and is highly concentrated among East Asians. Carriers of the ALDH2*2 variant commonly present alterations in a number of blood biomarkers, clinical measurements, biometrics, drug prescriptions, dietary habits and lifestyle behaviors, and they are also more susceptible to aldehyde-associated diseases, such as cancer and cardiovascular disease. However, the interaction between alcohol and ALDH2-related pathology is not clearly delineated. Furthermore, genetic evidence indicates that the ALDH2*2 variant has been favorably selected for in the past 2000-3000 years. It is therefore necessary to consider the disease risk and mechanism associated with ALDH2 deficiency, and to understand the possible beneficial or protective effect conferred by ALDH2 deficiency and whether the pleiotropic effects of ALDH2 variance are all mediated by alcohol use.

Genomic analyses of 10,376 individuals in the Westlake BioBank for Chinese (WBBC) pilot project

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

Crohn's Disease and Early Exposure to Thiopurines are Independent Risk Factors for Mosaic Chromosomal Alterations in Patients with Inflammatory Bowel Diseases

BACKGROUND AND AIMS

Mosaic chromosomal alterations [mCAs] increase the risk for haematopoietic malignancies and may be risk factors for several other diseases. Inflammatory bowel diseases [IBDs], including Crohn's disease [CD] and ulcerative colitis [UC], are associated with mCAs, and patients may be at risk for haematopoietic malignancy development and/or modification of IBD phenotypes. In the present study, we screened patients with IBD for the presence of mCAs and explored the possible pathophysiological and genetic risk factors for mCAs.

METHODS

We analysed mCAs in peripheral blood from 3339 patients with IBD and investigated the clinical and genetic risk factors for mCAs.

RESULTS

CD and exposure to thiopurines before the age of 20 years were identified as novel independent risk factors for mCAs [odds ratio = 2.15 and 5.68, p = 1.17e-2 and 1.60e-3, respectively]. In contrast, there were no significant associations of disease duration, anti-tumour necrosis factor alpha antibodies, or other clinical factors with mCAs. Gene ontology enrichment analysis revealed that genes specifically located in the mCAs in patients with CD were significantly associated with factors related to mucosal immune responses. A genome-wide association study revealed that ERBIN, CD96, and AC068672.2 were significantly associated with mCAs in patients with CD [p = 1.56e-8, 1.65e-8, and 4.92e-8, respectively].

CONCLUSIONS

The difference in mCAs between patients with CD and UC supports the higher incidence of haematopoietic malignancies in CD. Caution should be exercised when using thiopurines in young patients with IBD, particularly CD, in light of possible chromosomal alterations.

Association of drinking behaviors with TXNIP DNA methylation levels in leukocytes among the general Japanese population

Background: Thioredoxin-interacting protein (TXNIP) controls the cellular redox balance by binding to and inhibiting the expression and function of thioredoxin. DNA methylation of the TXNIP gene is involved in the regulation of TXNIP mRNA expression. Changes in TXNIP DNA methylation levels are associated with the development of various diseases such as type 2 diabetes mellitus (T2DM). However, few studies have focused on the influence of lifestyle factors such as alcohol intake on TXNIP DNA methylation.Objectives: This research examines the association of drinking behaviors with TXNIP DNA methylation levels in the general Japanese population.Methods: We conducted a cross-sectional study of 404 subjects (176 males and 228 females) who were divided into non-, moderate and heavy drinkers based on self-reported drinking behaviors. TXNIP DNA methylation levels in leukocytes were determined using a pyrosequencing assay.Results: The mean TXNIP DNA methylation level in heavy drinkers (74.2%) was significantly lower than that in non- and moderate drinkers (non: 77.7%, p < .001; moderate: 76.6%, p = .011). Multivariable linear regression analysis showed that log-transformed values of daily (b = -1.34; p < .001) and cumulative (b = -1.06; p = .001) alcohol consumption were associated with decreased TXNIP DNA methylation levels.Conclusion: TXNIP DNA methylation levels in heavy drinkers was lower than in non- and- moderate drinkers. Decreased TXNIP DNA methylation level increases the expression of TXNIP and elevates the risk of developing of diseases such as T2DM. Therefore, decreasing alcohol use in heavy drinkers may lessen the likelihood of some alcohol-related illnesses moderated through TXNIP DNA methylation.

Positive selection acts on regulatory genetic variants in populations of European ancestry that affect ALDH2 gene expression

ALDH2 is a key enzyme in alcohol metabolism that protects cells from acetaldehyde toxicity. Using iHS, iSAFE and FST statistics, we identified regulatory acting variants affecting ALDH2 gene expression under positive selection in populations of European ancestry. Several SNPs (rs3184504, rs4766578, rs10774625, rs597808, rs653178, rs847892, rs2013002) that function as eQTLs for ALDH2 in various tissues showed evidence of strong positive selection. Very large pairwise FST values indicated high genetic differentiation at these loci between populations of European ancestry and populations of other global ancestries. Estimating the timing of positive selection on the beneficial alleles suggests that these variants were recently adapted approximately 3000-3700 years ago. The derived beneficial alleles are in complete linkage disequilibrium with the derived ALDH2 promoter variant rs886205, which is associated with higher transcriptional activity. The SNPs rs4766578 and rs847892 are located in binding sequences for the transcription factor HNF4A, which is an important regulatory element of ALDH2 gene expression. In contrast to the missense variant ALDH2 rs671 (ALDH2*2), which is common only in East Asian populations and is associated with greatly reduced enzyme activity and alcohol intolerance, the beneficial alleles of the regulatory variants identified in this study are associated with increased expression of ALDH2. This suggests adaptation of Europeans to higher alcohol consumption.

Decoding the diversity of killer immunoglobulin-like receptors by deep sequencing and a high-resolution imputation method

The killer cell immunoglobulin-like receptor (KIR) recognizes human leukocyte antigen (HLA) class I molecules and modulates the function of natural killer cells. Despite its role in immunity, the complex genomic structure has limited a deep understanding of the KIR genomic landscape. Here we conduct deep sequencing of 16 KIR genes in 1,173 individuals. We devise a bioinformatics pipeline incorporating copy number estimation and insertion or deletion (indel) calling for high-resolution KIR genotyping. We define 118 alleles in 13 genes and demonstrate a linkage disequilibrium structure within and across KIR centromeric and telomeric regions. We construct a KIR imputation reference panel (n_reference = 689, imputation accuracy = 99.7%), apply it to biobank genotype (n_total = 169,907), and perform phenome-wide association studies of 85 traits. We observe a dearth of genome-wide significant associations, even in immune traits implicated previously to be associated with KIR (the smallest p = 1.5 × 10^-4). Our pipeline presents a broadly applicable framework to evaluate innate immunity in large-scale datasets.

Whole exome analysis of patients in Japan with hearing loss reveals high heterogeneity among responsible and novel candidate genes

Background

Heterogeneous genetic loci contribute to hereditary hearing loss; more than 100 deafness genes have been identified, and the number is increasing. To detect pathogenic variants in multiple deafness genes, in addition to novel candidate genes associated with hearing loss, whole exome sequencing (WES), followed by analysis prioritizing genes categorized in four tiers, were applied.

Results

Trios from families with non-syndromic or syndromic hearing loss (n = 72) were subjected to WES. After segregation analysis and interpretation according to American College of Medical Genetics and Genomics guidelines, candidate pathogenic variants in 11 previously reported deafness genes (STRC, MYO15A, CDH23, PDZD7, PTPN11, SOX10, EYA1, MYO6, OTOF, OTOG, and ZNF335) were identified in 21 families. Discrepancy between pedigree inheritance and genetic inheritance was present in one family. In addition, eight genes (SLC12A2, BAIAP2L2, HKDC1, SVEP1, CACNG1, GTPBP4, PCNX2, and TBC1D8) were screened as single candidate genes in 10 families.

Conclusions

Our findings demonstrate that four-tier assessment of WES data is efficient and can detect novel candidate genes associated with hearing loss, in addition to pathogenic variants of known deafness genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02262-4.

Qatar Genome: Insights on Genomics from the Middle East

Despite recent biomedical breakthroughs and large genomic studies growing momentum, the Middle Eastern population, home to over 400 million people, is under-represented in the human genome variation databases. Here we describe insights from phase 1 of the Qatar Genome Program with whole genome sequenced 6,047 individuals from Qatar. We identified more than 88 million variants of which 24 million are novel and 23 million are singletons. Consistent with the high consanguinity and founder effects in the region, we found that several rare deleterious variants were more common in the Qatari population while others seem to provide protection against diseases and have shaped the genetic architecture of adaptive phenotypes. These results highlight the value of our data as a resource to advance genetic studies in the Arab and neighbouring Middle Eastern populations and will significantly boost the current efforts to improve our understanding of global patterns of human variations, human history and genetic contributions to health and diseases in diverse populations. This article is protected by copyright. All rights reserved.

Genetic Background of Mesalamine-induced Fever and Diarrhea in Japanese Patients with Inflammatory Bowel Disease

BACKGROUND

Some patients with inflammatory bowel disease (IBD) who were under mesalamine treatment develop adverse reactions called "mesalamine allergy," which includes high fever and worsening diarrhea. Currently, there is no method to predict mesalamine allergy. Pharmacogenomic approaches may help identify these patients. Here we analyzed the genetic background of mesalamine intolerance in the first genome-wide association study of Japanese patients with IBD.

METHODS

Two independent pharmacogenetic IBD cohorts were analyzed: the MENDEL (n = 1523; as a discovery set) and the Tohoku (n = 788; as a replication set) cohorts. Genome-wide association studies were performed in each population, followed by a meta-analysis. In addition, we constructed a polygenic risk score model and combined genetic and clinical factors to model mesalamine intolerance.

RESULTS

In the combined cohort, mesalamine-induced fever and/or diarrhea was significantly more frequent in ulcerative colitis vs Crohn's disease. The genome-wide association studies and meta-analysis identified one significant association between rs144384547 (upstream of RGS17) and mesalamine-induced fever and diarrhea (P = 7.21e-09; odds ratio = 11.2). The estimated heritability of mesalamine allergy was 25.4%, suggesting a significant correlation with the genetic background. Furthermore, a polygenic risk score model was built to predict mesalamine allergy (P = 2.95e-2). The combined genetic/clinical prediction model yielded a higher area under the curve than did the polygenic risk score or clinical model alone (area under the curve, 0.89; sensitivity, 71.4%; specificity, 90.8%).

CONCLUSIONS

Mesalamine allergy was more common in ulcerative colitis than in Crohn's disease. We identified a novel genetic association with and developed a combined clinical/genetic model for this adverse event.

Analyzing the Korean reference genome with meta-imputation increased the imputation accuracy and spectrum of rare variants in the Korean population

Genotype imputation is essential for enhancing the power of association-mapping and discovering rare and indels that are missed by most genotyping arrays. Imputation analysis can be more accurate with a population-specific reference panel or a multi-ethnic reference panel with numerous samples. The National Institute of Health, Republic of Korea, initiated the Korean Reference Genome (KRG) project to identify variants in whole-genome sequences of ∼20,000 Korean participants. In the pilot phase, we analyzed the data from 1,490 participants. The genetic characteristics and imputation performance of the KRG were compared with those of the 1,000 Genomes Project Phase 3, GenomeAsia 100K Project, ChinaMAP, NARD, and TOPMed reference panels. For comparison analysis, genotype panels were artificially generated using whole-genome sequencing data from combinations of four different ancestries (Korean, Japanese, Chinese, and European) and two population-specific optimized microarrays (Korea Biobank Array and UK Biobank Array). The KRG reference panel performed best for the Korean population (R ² = 0.78-0.84, percentage of well-imputed is 91.9% for allele frequency >5%), although the other reference panels comprised a larger number of samples with genetically different background. By comparing multiple reference panels and multi-ethnic genotype panels, optimal imputation was obtained using reference panels from genetically related populations and a population-optimized microarray. Indeed, the reference panels of KRG and TOPMed showed the best performance when applied to the genotype panels of KBA (R ² = 0.84) and UKB (R ² = 0.87), respectively. Using a meta-imputation approach to merge imputation results from different reference panels increased the imputation accuracy for rare variants (∼7%) and provided additional well-imputed variants (∼20%) with comparable imputation accuracy to that of the KRG. Our results demonstrate the importance of using a population-specific reference panel and meta-imputation to assess a substantial number of accurately imputed rare variants.

The A118G single-nucleotide polymorphism in OPRM1 is a risk factor for asthma severity

BACKGROUND

Although population studies have implicated emotional burden in asthma severity, the underlying genetic risk factors are not completely understood. We aimed to evaluate the genetic influence of a functional single-nucleotide polymorphism (SNP) in the stress-related μ-opioid receptor gene (OPRM1; A118G SNP, rs1799971) on asthma severity.

METHODS

We initially assessed disease severity in asthmatic outpatients carrying A118G. Using an ovalbumin-induced experimental asthma rodent model harboring the functionally equivalent SNP, we investigated the mechanism by which this SNP influences the allergic immune response.

RESULTS

Among 292 outpatients, 168 underwent airway hyperresponsiveness (AHR) to methacholine testing. Compared with patients carrying the AA and AG genotypes, those carrying the GG genotype exhibited enhanced AHR. The stress levels were presumed to be moderate among patients and were comparable among genotypes. Compared with Oprm1 AA mice, GG mice demonstrated aggravated asthma-related features and increased pulmonary interleukin-4⁺CD4⁺ effector and effector memory T cells under everyday life stress conditions. Intraperitoneal naloxone methiodide injection reduced effector CD4⁺ T cell elevation associated with increased eosinophil numbers in bronchoalveolar lavage fluid of GG mice to the levels in AA mice, suggesting that elevated Th2 cell generation in the bronchial lymph node (BLN) of GG mice induces enhanced eosinophilic inflammation.

CONCLUSIONS

Without forced stress exposure, patients with asthma carrying the OPRM1 GG genotype exhibit enhanced AHR, attributable to enhanced Th2 cell differentiation in the regional lymph node. Further research is necessary to elucidate the role of the OPRM1 A118G genotype in the Th2 cell differentiation pathway in the BLN.

Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population

Understanding the genetic effects on non-coding RNA (ncRNA) expression facilitates functional characterization of disease-associated genetic loci. Among several classes of ncRNAs, microRNAs (miRNAs) are key post-transcriptional gene regulators. Despite its biological importance, previous studies on the genetic architecture of miRNA expression focused mostly on the European individuals, underrepresented in other populations. Here, we mapped miRNA expression quantitative trait loci (miRNA-eQTL) for 343 miRNAs in 141 Japanese using small RNA sequencing (sRNA-seq) and whole-genome sequencing (WGS), identifying 1275 cis-miRNA-eQTL variants for 40 miRNAs (false discovery rate < 0.2). Of these, 25 miRNAs having eQTL were unreported in the European studies, including 5 miRNAs with their lead variant monomorphic in the European populations, which demonstrates the value of miRNA-eQTL analysis in diverse ancestral populations. MiRNAs with eQTL effect showed allele-specific expression (ASE) (e.g. miR-146a-3p), and ASE analysis further detected cis-regulatory variants not captured by the conventional miRNA-eQTL mapping (e.g. miR-933). We identified a copy number variation (CNV) associated with miRNA expression (e.g. miR-570-3p, P = 7.2 × 10-6), which contributes to a more comprehensive landscape of miRNA-eQTLs. To elucidate a post-transcriptional modification in miRNAs, we created a catalog of miRNA-editing sites, including ten canonical and six non-canonical sites. Finally, by integrating the miRNA-eQTLs and Japanese genome-wide association studies of 25 complex traits (mean n = 192 833), we conducted a transcriptome-wide association study (TWAS), identifying miR-1908-5p as a potential mediator for adult height, colorectal cancer, and type 2 diabetes (P < 9.1 × 10-5). Our study broadens the population diversity in ncRNA-eQTL studies and contributes to functional annotation of disease-associated loci found in non-European populations.

Envisioning the next human genome reference

Summary: We provide an Editorial perspective on approaches to improve ethnic representation in the human genome reference sequence, enabling its widespread use in genomic studies and precision medicine to benefit all peoples.

Comparison of sequencing data processing pipelines and application to underrepresented African human populations

Background

Population genetic studies of humans make increasing use of high-throughput sequencing in order to capture diversity in an unbiased way. There is an abundance of sequencing technologies, bioinformatic tools and the available genomes are increasing in number. Studies have evaluated and compared some of these technologies and tools, such as the Genome Analysis Toolkit (GATK) and its “Best Practices” bioinformatic pipelines. However, studies often focus on a few genomes of Eurasian origin in order to detect technical issues. We instead surveyed the use of the GATK tools and established a pipeline for processing high coverage full genomes from a diverse set of populations, including Sub-Saharan African groups, in order to reveal challenges from human diversity and stratification.

Results

We surveyed 29 studies using high-throughput sequencing data, and compared their strategies for data pre-processing and variant calling. We found that processing of data is very variable across studies and that the GATK “Best Practices” are seldom followed strictly. We then compared three versions of a GATK pipeline, differing in the inclusion of an indel realignment step and with a modification of the base quality score recalibration step. We applied the pipelines on a diverse set of 28 individuals. We compared the pipelines in terms of count of called variants and overlap of the callsets. We found that the pipelines resulted in similar callsets, in particular after callset filtering. We also ran one of the pipelines on a larger dataset of 179 individuals. We noted that including more individuals at the joint genotyping step resulted in different counts of variants. At the individual level, we observed that the average genome coverage was correlated to the number of variants called.

Conclusions

We conclude that applying the GATK “Best Practices” pipeline, including their recommended reference datasets, to underrepresented populations does not lead to a decrease in the number of called variants compared to alternative pipelines. We recommend to aim for coverage of > 30X if identifying most variants is important, and to work with large sample sizes at the variant calling stage, also for underrepresented individuals and populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04407-x.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

A cross-population atlas of genetic associations for 220 human phenotypes

Current genome-wide association studies do not yet capture sufficient diversity in populations and scope of phenotypes. To expand an atlas of genetic associations in non-European populations, we conducted 220 deep-phenotype genome-wide association studies (diseases, biomarkers and medication usage) in BioBank Japan (n = 179,000), by incorporating past medical history and text-mining of electronic medical records. Meta-analyses with the UK Biobank and FinnGen (n_total = 628,000) identified ~5,000 new loci, which improved the resolution of the genomic map of human traits. This atlas elucidated the landscape of pleiotropy as represented by the major histocompatibility complex locus, where we conducted HLA fine-mapping. Finally, we performed statistical decomposition of matrices of phenome-wide summary statistics, and identified latent genetic components, which pinpointed responsible variants and biological mechanisms underlying current disease classifications across populations. The decomposed components enabled genetically informed subtyping of similar diseases (for example, allergic diseases). Our study suggests a potential avenue for hypothesis-free re-investigation of human diseases through genetics.

Unhealthy lifestyles and brain condition: Examining the relations of BMI, living alone, alcohol intake, short sleep, smoking, and lack of exercise with gray matter volume

Unhealthy lifestyles are damaging to the brain. Previous studies have indicated that body mass index (BMI), alcohol intake, short sleep, smoking, and lack of exercise are negatively associated with gray matter volume (GMV). Living alone has also been found to be related to GMV through lowered subjective happiness. However, to our knowledge, no GMV study has dealt with these unhealthy lifestyles simultaneously. By our analyses based on 142 healthy Japanese participants, BMI, alcohol intake, living alone, and short sleep were negatively associated with the gray-matter brain healthcare quotient (GM-BHQ), an MRI-based normalized GMV, after controlling for age, sex, and facility, not only individually but also when they were entered into a single regression model. Moreover, there were small but significant differences in the proportion of the variance for GM-BHQ explained by variables in a regression model (measured by R squared) between when these unhealthy variables were entered in an equation at the same time and when they were entered separately, with the former larger than the latter. However, smoking and lack of exercise were not significantly associated with GM-BHQ. Results indicate that some kinds of unhealthy lifestyles are somewhat harmful on their own, but may become more noxious to brain condition if practiced simultaneously, although its difference may not be large. To our knowledge, this study is the first to show that overlapping unhealthy lifestyles affects the brains of healthy adults.

An exploratory assessment of the applicability of direct-to-consumer genetic testing to translational research in Japan

Objective

In order to assess the applicability of a direct-to-consumer (DTC) genetic testing to translational research for obtaining new knowledge on relationships between drug target genes and diseases, we examined possibility of these data by associating SNPs and disease related phenotype information collected from healthy individuals.

Results

A total of 12,598 saliva samples were collected from the customers of commercial service for SNPs analysis and web survey were conducted to collect phenotype information. The collected dataset revealed similarity to the Japanese data but distinguished differences to other populations of all dataset of the 1000 Genomes Project. After confirmation of a well-known relationship between ALDH2 and alcohol-sensitivity, Phenome-Wide Association Study (PheWAS) was performed to find association between pre-selected drug target genes and all the phenotypes. Association was found between GRIN2B and multiple phenotypes related to depression, which is considered reliable based on previous reports on the biological function of GRIN2B protein and its relationship with depression. These results suggest possibility of using SNPs and phenotype information collected from healthy individuals as a translational research tool for drug discovery to find relationship between a gene and a disease if it is possible to extract individuals in pre-disease states by properly designed questionnaire.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05696-4.

Fine-Scale Genetic Structure and Demographic History in the Miyako Islands of the Ryukyu Archipelago

The Ryukyu Archipelago is located in the southwest of the Japanese islands and is composed of dozens of islands, grouped into the Miyako Islands, Yaeyama Islands, and Okinawa Islands. Based on the results of principal component analysis on genome-wide single-nucleotide polymorphisms, genetic differentiation was observed among the island groups of the Ryukyu Archipelago. However, a detailed population structure analysis of the Ryukyu Archipelago has not yet been completed. We obtained genomic DNA samples from 1,240 individuals living in the Miyako Islands, and we genotyped 665,326 single-nucleotide polymorphisms to infer population history within the Miyako Islands, including Miyakojima, Irabu, and Ikema islands. The haplotype-based analysis showed that populations in the Miyako Islands were divided into three subpopulations located on Miyakojima northeast, Miyakojima southwest, and Irabu/Ikema. The results of haplotype sharing and the D statistics analyses showed that the Irabu/Ikema subpopulation received gene flows different from those of the Miyakojima subpopulations, which may be related with the historically attested immigration during the Gusuku period (900 − 500 BP). A coalescent-based demographic inference suggests that the Irabu/Ikema population firstly split away from the ancestral Ryukyu population about 41 generations ago, followed by a split of the Miyako southwest population from the ancestral Ryukyu population (about 16 generations ago), and the differentiation of the ancestral Ryukyu population into two populations (Miyako northeast and Okinawajima populations) about seven generations ago. Such genetic information is useful for explaining the population history of modern Miyako people and must be taken into account when performing disease association studies.

Deep sequencing of 1320 genes reveals the landscape of protein-truncating variants and their contribution to psoriasis in 19,973 Chinese individuals

Protein-truncating variants (PTVs) have important impacts on phenotype diversity and disease. However, their population genetics characteristics in more globally diverse populations are not well defined. Here, we describe patterns of PTVs in 1320 genes sequenced in 10,539 healthy controls and 9434 patients with psoriasis, all of Han Chinese ancestry. We identify 8720 PTVs, of which 77% are novel, and estimate 88% of all PTVs are deleterious and subject to purifying selection. Furthermore, we show that individuals with psoriasis have a significantly higher burden of PTVs compared to controls (P = 0.02). Finally, we identified 18 PTVs in 14 genes with unusually high levels of population differentiation, consistent with the action of local adaptation. Our study provides insights into patterns and consequences of PTVs.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

Whole genome sequencing of 45 Japanese patients with intellectual disability

Intellectual disability (ID) is characterized by significant limitations in both intellectual functioning and adaptive behaviors, originating before the age of 18 years. However, the genetic etiologies of ID are still incompletely elucidated due to the wide range of clinical and genetic heterogeneity. Whole genome sequencing (WGS) has been applied as a single-step clinical diagnostic tool for ID because it detects genetic variations with a wide range of resolution from single nucleotide variants (SNVs) to structural variants (SVs). To explore the causative genes for ID, we employed WGS in 45 patients from 44 unrelated Japanese families and performed a stepwise screening approach focusing on the coding variants in the genes. Here, we report 12 pathogenic and likely pathogenic variants: seven heterozygous variants of ADNP, SATB2, ANKRD11, PTEN, TCF4, SPAST, and KCNA2, three hemizygous variants of SMS, SLC6A8, and IQSEC2, and one homozygous variant in AGTPBP1. Of these, four were considered novel. Furthermore, a novel 76 kb deletion containing exons 1 and 2 in DYRK1A was identified. We confirmed the clinical and genetic heterogeneity and high frequency of de novo causative variants (8/12, 66.7%). This is the first report of WGS analysis in Japanese patients with ID. Our results would provide insight into the correlation between novel variants and expanded phenotypes of the disease.

Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases

Variants of transcription factor binding sites (TFBSs) constitute an important part of the human genome. Current evidence demonstrates close links between nucleotides within TFBSs and gene expression. There are multiple pathways through which genomic sequences located in TFBSs regulate gene expression, and recent genome-wide association studies have shown the biological significance of TFBS variation in human phenotypes. However, numerous challenges remain in the study of TFBS polymorphisms. This article aims to cover the current state of understanding as regards the genomic features of TFBSs and TFBS variants; the mechanisms through which TFBS variants regulate gene expression; the approaches to studying the effects of nucleotide changes that create or disrupt TFBSs; the challenges faced in studies of TFBS sequence variations; the effects of natural selection on collections of TFBSs; in addition to the insights gained from the study of TFBS alleles related to gout, its associated comorbidities (increased body mass index, chronic kidney disease, diabetes, dyslipidemia, coronary artery disease, ischemic heart disease, hypertension, hyperuricemia, osteoporosis, and prostate cancer), and the treatment responses of patients.

Progress in Defining the Genetic Contribution to Type 2 Diabetes in Individuals of East Asian Ancestry

PURPOSE OF REVIEW

Prevalence of type 2 diabetes (T2D) and progression of complications differ between worldwide populations. While obesity is a major contributing risk factor, variations in physiological manifestations, e.g., developing T2D at lower body mass index in some populations, suggest other contributing factors. Early T2D genetic associations were mostly discovered in European ancestry populations. This review describes the progression of genetic discoveries associated with T2D in individuals of East Asian ancestry in the last 10 years and highlights the shared genetic susceptibility between the population groups and additional insights into genetic contributions to T2D.

RECENT FINDINGS

Through increased sample size and power, new genetic associations with T2D were discovered in East Asian ancestry populations, often with higher allele frequencies than European ancestry populations. As we continue to generate maps of T2D-associated variants across diverse populations, there will be a critical need to expand and diversify other omics resources to enable integration for clinical translation.

Genome-Wide Natural Selection Signatures Are Linked to Genetic Risk of Modern Phenotypes in the Japanese Population

Elucidation of natural selection signatures and relationships with phenotype spectra is important to understand adaptive evolution of modern humans. Here, we conducted a genome-wide scan of selection signatures of the Japanese population by estimating locus-specific time to the most recent common ancestor using the ascertained sequentially Markovian coalescent (ASMC), from the biobank-based large-scale genome-wide association study data of 170,882 subjects. We identified 29 genetic loci with selection signatures satisfying the genome-wide significance. The signatures were most evident at the alcohol dehydrogenase (ADH) gene cluster locus at 4q23 (P_ASMC = 2.2 × 10⁻³⁶), followed by relatively strong selection at the FAM96A (15q22), MYOF (10q23), 13q21, GRIA2 (4q32), and ASAP2 (2p25) loci (P_ASMC < 1.0 × 10⁻¹⁰). The additional analysis interrogating extended haplotypes (integrated haplotype score) showed robust concordance of the detected signatures, contributing to fine-mapping of the genes, and provided allelic directional insights into selection pressure (e.g., positive selection for ADH1B-Arg48His and HLA-DPB1*04:01). The phenome-wide selection enrichment analysis with the trait-associated variants identified a variety of the modern human phenotypes involved in the adaptation of Japanese. We observed population-specific evidence of enrichment with the alcohol-related phenotypes, anthropometric and biochemical clinical measurements, and immune-related diseases, differently from the findings in Europeans using the UK Biobank resource. Our study demonstrated population-specific features of the selection signatures in Japanese, highlighting a value of the natural selection study using the nation-wide biobank-scale genome and phenotype data.

A genome-wide association study on confection consumption in a Japanese population: the Japan Multi-Institutional Collaborative Cohort Study

Differences in individual eating habits may be influenced by genetic factors, in addition to cultural, social or environmental factors. Previous studies suggested that genetic variants within sweet taste receptor genes family were associated with sweet taste perception and the intake of sweet foods. The aim of this study was to conduct a genome-wide association study (GWAS) to find genetic variations that affect confection consumption in a Japanese population. We analysed GWAS data on confection consumption using 14 073 participants from the Japan Multi-Institutional Collaborative Cohort study. We used a semi-quantitative FFQ to estimate food intake that was validated previously. Association of the imputed variants with confection consumption was performed by linear regression analysis with adjustments for age, sex, total energy intake and principal component analysis components 1-3. Furthermore, the analysis was repeated adjusting for alcohol intake (g/d) in addition to the above-described variables. We found 418 SNP located in 12q24 that were associated with confection consumption. SNP with the ten lowest P-values were located on nine genes including at the BRAP, ACAD10 and aldehyde dehydrogenase 2 regions on 12q24.12-13. After adjustment for alcohol intake, no variant was associated with confections intake with genome-wide significance. In conclusion, we found a significant number of SNP located on 12q24 genes that were associated with confections intake before adjustment for alcohol intake. However, all of them lost statistical significance after adjustment for alcohol intake.

Association of Rare CYP39A1 Variants With Exfoliation Syndrome Involving the Anterior Chamber of the Eye

IMPORTANCE

Exfoliation syndrome is a systemic disorder characterized by progressive accumulation of abnormal fibrillar protein aggregates manifesting clinically in the anterior chamber of the eye. This disorder is the most commonly known cause of glaucoma and a major cause of irreversible blindness.

OBJECTIVE

To determine if exfoliation syndrome is associated with rare, protein-changing variants predicted to impair protein function.

DESIGN, SETTING, AND PARTICIPANTS

A 2-stage, case-control, whole-exome sequencing association study with a discovery cohort and 2 independently ascertained validation cohorts. Study participants from 14 countries were enrolled between February 1999 and December 2019. The date of last clinical follow-up was December 2019. Affected individuals had exfoliation material on anterior segment structures of at least 1 eye as visualized by slit lamp examination. Unaffected individuals had no signs of exfoliation syndrome.

EXPOSURES

Rare, coding-sequence genetic variants predicted to be damaging by bioinformatic algorithms trained to recognize alterations that impair protein function.

MAIN OUTCOMES AND MEASURES

The primary outcome was the presence of exfoliation syndrome. Exome-wide significance for detected variants was defined as P < 2.5 × 10-6. The secondary outcomes included biochemical enzymatic assays and gene expression analyses.

RESULTS

The discovery cohort included 4028 participants with exfoliation syndrome (median age, 78 years [interquartile range, 73-83 years]; 2377 [59.0%] women) and 5638 participants without exfoliation syndrome (median age, 72 years [interquartile range, 65-78 years]; 3159 [56.0%] women). In the discovery cohort, persons with exfoliation syndrome, compared with those without exfoliation syndrome, were significantly more likely to carry damaging CYP39A1 variants (1.3% vs 0.30%, respectively; odds ratio, 3.55 [95% CI, 2.07-6.10]; P = 6.1 × 10-7). This outcome was validated in 2 independent cohorts. The first validation cohort included 2337 individuals with exfoliation syndrome (median age, 74 years; 1132 women; n = 1934 with demographic data) and 2813 individuals without exfoliation syndrome (median age, 72 years; 1287 women; n = 2421 with demographic data). The second validation cohort included 1663 individuals with exfoliation syndrome (median age, 75 years; 587 women; n = 1064 with demographic data) and 3962 individuals without exfoliation syndrome (median age, 74 years; 951 women; n = 1555 with demographic data). Of the individuals from both validation cohorts, 5.2% with exfoliation syndrome carried CYP39A1 damaging alleles vs 3.1% without exfoliation syndrome (odds ratio, 1.82 [95% CI, 1.47-2.26]; P < .001). Biochemical assays classified 34 of 42 damaging CYP39A1 alleles as functionally deficient (median reduction in enzymatic activity compared with wild-type CYP39A1, 94.4% [interquartile range, 78.7%-98.2%] for the 34 deficient variants). CYP39A1 transcript expression was 47% lower (95% CI, 30%-64% lower; P < .001) in ciliary body tissues from individuals with exfoliation syndrome compared with individuals without exfoliation syndrome.

CONCLUSIONS AND RELEVANCE

In this whole-exome sequencing case-control study, presence of exfoliation syndrome was significantly associated with carriage of functionally deficient CYP39A1 sequence variants. Further research is needed to understand the clinical implications of these findings.

Genetic determinants of risk in autoimmune pulmonary alveolar proteinosis

Pulmonary alveolar proteinosis (PAP) is a devastating lung disease caused by abnormal surfactant homeostasis, with a prevalence of 6-7 cases per million population worldwide. While mutations causing hereditary PAP have been reported, the genetic basis contributing to autoimmune PAP (aPAP) has not been thoroughly investigated. Here, we conducted a genome-wide association study of aPAP in 198 patients and 395 control participants of Japanese ancestry. The common genetic variant, rs138024423 at 6p21, in the major-histocompatibility-complex (MHC) region was significantly associated with disease risk (Odds ratio [OR] = 5.2; P = 2.4 × 10-12). HLA fine-mapping revealed that the common HLA class II allele, HLA-DRB1*08:03, strongly drove this signal (OR = 4.8; P = 4.8 × 10-12), followed by an additional independent risk allele at HLA-DPβ1 amino acid position 8 (OR = 0.28; P = 3.4 × 10-7). HLA-DRB1*08:03 was also associated with an increased level of anti-GM-CSF antibody, a key driver of the disease (β = 0.32; P = 0.035). Our study demonstrated a heritable component of aPAP, suggesting an underlying genetic predisposition toward an abnormal antibody production.

Detection of selection signatures for response to Aleutian mink disease virus infection in American mink

Aleutian disease (AD) is the most significant health issue for farmed American mink. The objective of this study was to identify the genomic regions subjected to selection for response to infection with Aleutian mink disease virus (AMDV) in American mink using genotyping by sequencing (GBS) data. A total of 225 black mink were inoculated with AMDV and genotyped using a GBS assay based on the sequencing of ApeKI-digested libraries. Five AD-characterized phenotypes were used to assign animals to pairwise groups. Signatures of selection were detected using integrated measurement of fixation index (F_ST) and nucleotide diversity (θπ), that were validated by haplotype-based (hap-FLK) test. The total of 99 putatively selected regions harbouring 63 genes were detected in different groups. The gene ontology revealed numerous genes related to immune response (e.g. TRAF3IP2, WDR7, SWAP70, CBFB, and GPR65), liver development (e.g. SULF2, SRSF5) and reproduction process (e.g. FBXO5, CatSperβ, CATSPER4, and IGF2R). The hapFLK test supported two strongly selected regions that contained five candidate genes related to immune response, virus–host interaction, reproduction and liver regeneration. This study provided the first map of putative selection signals of response to AMDV infection in American mink, bringing new insights into genomic regions controlling the AD phenotypes.

Genome-wide SNP data of Izumo and Makurazaki populations support inner-dual structure model for origin of Yamato people

The “Dual Structure” model on the formation of the modern Japanese population assumes that the indigenous hunter-gathering population (symbolized as Jomon people) admixed with rice-farming population (symbolized as Yayoi people) who migrated from the Asian continent after the Yayoi period started. The Jomon component remained high both in Ainu and Okinawa people who mainly reside in northern and southern Japan, respectively, while the Yayoi component is higher in the mainland Japanese (Yamato people). The model has been well supported by genetic data, but the Yamato population was mostly represented by people from Tokyo area. We generated new genome-wide SNP data using Japonica Array for 45 individuals in Izumo City of Shimane Prefecture and for 72 individuals in Makurazaki City of Kagoshima Prefecture in Southern Kyushu, and compared these data with those of other human populations in East Asia, including BioBank Japan data. Using principal component analysis, phylogenetic network, and f4 tests, we found that Izumo, Makurazaki, and Tohoku populations are slightly differentiated from Kanto (including Tokyo), Tokai, and Kinki regions. These results suggest the substructure within Mainland Japanese maybe caused by multiple migration events from the Asian continent following the Jomon period, and we propose a modified version of “Dual Structure” model called the “Inner-Dual Structure” model.

Novel candidates of pathogenic variants of the BRCA1 and BRCA2 genes from a dataset of 3,552 Japanese whole genomes (3.5KJPNv2)

Identification of the population frequencies of definitely pathogenic germline variants in two major hereditary breast and ovarian cancer syndrome (HBOC) genes, BRCA1/2, is essential to estimate the number of HBOC patients. In addition, the identification of moderately penetrant HBOC gene variants that contribute to increasing the risk of breast and ovarian cancers in a population is critical to establish personalized health care. A prospective cohort subjected to genome analysis can provide both sets of information. Computational scoring and prospective cohort studies may help to identify such likely pathogenic variants in the general population. We annotated the variants in the BRCA1 and BRCA2 genes from a dataset of 3,552 whole-genome sequences obtained from members of a prospective cohorts with genome data in the Tohoku Medical Megabank Project (TMM) with InterVar software. Computational impact scores (CADD_phred and Eigen_raw) and minor allele frequencies (MAFs) of pathogenic (P) and likely pathogenic (LP) variants in ClinVar were used for filtration criteria. Familial predispositions to cancers among the 35,000 TMM genome cohort participants were analyzed to verify the identified pathogenicity. Seven potentially pathogenic variants were newly identified. The sisters of carriers of these moderately deleterious variants and definite P and LP variants among members of the TMM prospective cohort showed a statistically significant preponderance for cancer onset, from the self-reported cancer history. Filtering by computational scoring and MAF is useful to identify potentially pathogenic variants in BRCA genes in the Japanese population. These results should help to follow up the carriers of variants of uncertain significance in the HBOC genes in the longitudinal prospective cohort study.

Whole genome sequence analyses of eGFR in 23,732 people representing multiple ancestries in the NHLBI trans-omics for precision medicine (TOPMed) consortium

BACKGROUND

Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants.

METHODS

We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity.

FINDINGS

When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10-11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10-9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10-9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10-9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10-9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants.

INTERPRETATION

This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

A Mendelian randomization study identified obesity as a causal risk factor of uterine endometrial cancer in Japanese

Causal inference is one of the challenges in epidemiologic studies. Gynecologic diseases have been reported to have association with obesity, however the causality remained controversial except for uterine endometrial cancer. We conducted two-sample Mendelian randomization (MR) analysis using the large-scale genome-wide association study (GWAS) results of gynecologic diseases and body mass index (BMI) in the Japanese population to assess causal effect of BMI on gynecologic diseases. We first conducted GWAS of ovarian cancer, uterine endometrial cancer, uterine cervical cancer, endometriosis, and uterine fibroid (n = 647, 909, 538, 5236, and 645 cases, respectively, and 39 556 shared female controls), and BMI (81 610 males and non-overlapping 23 924 females). We then applied two-sample MR using 74 BMI-associated variants as instrumental variables. We observed significant causal effect of increased BMI on uterine endometrial cancer (β = 0.735, P = .0010 in inverse variance-weighted analysis), which is concordant with results of European studies. Causal effect of obesity was not apparent in the other gynecologic diseases tested. Our MR analyses provided strong evidence of the causal role of obesity in gynecologic diseases etiology, and suggested a possible preventive effect of intervention for obesity.

Causal Inference between Rheumatoid Arthritis and Breast Cancer in East Asian and European Population: A Two-Sample Mendelian Randomization

Simple Summary

Rheumatoid arthritis (RA) is one of the chronic autoimmune diseases that affects about 0.5 to 1.0% of the general population worldwide. The main symptom of RA is the destruction of the synovial joint, leading to a reduced quality of life and increased mortality. RA may be accompanied by several comorbidities, on which several studies have been conducted on the association between RA and breast cancer. However, the association between RA and breast cancer has shown different directions and has not been clearly established. In this study, we tried to determine whether RA had a causal effect on breast cancer using Mendelian randomization (MR) analysis, but causal evidence was not found. Therefore, additional studies are needed to determine whether RA patients are at high risk of breast cancer, based on large-scale cohorts to validate these results.

Abstract

Previous studies have been reported that the association between rheumatoid arthritis (RA) and breast cancer remains inconclusive. A two-sample Mendelian randomization (MR) analysis can reveal the potential causal association between exposure and outcome. A two-sample MR analysis using the penalized robust inverse variance weighted (PRIVW) method was performed to analyze the association between RA and breast cancer risk based on the summary statistics of six genome-wide association studies (GWAS) targeting RA in an East Asian population along with summary statistics of the BioBank Japan (BBJ), Breast Cancer Association Consortium (BCAC), and Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) targeting breast cancer. We found that the direction of the effect of RA on breast cancer varied among GWAS-summary data from BBJ, BCAC, and CIMBA. Significant horizontal pleiotropy based on a penalized robust MR-Egger regression was observed only for BBJ and CIMBA BRCA2 carriers. As the results of the two-sample MR analyses were inconsistent, the causal association between RA and breast cancer was inconclusive. The biological mechanisms explaining the relationship between RA and breast cancer were unclear in Asian as well as in Caucasians. Further studies using large-scale patient cohorts are required for the validation of these results.

Prevalence and Spectrum of Pathogenic Germline Variants in Japanese Patients With Early-Onset Colorectal, Breast, and Prostate Cancer

PUPOSE

We investigated the prevalence and spectrum of pathogenic germline variants in patients with early-onset colorectal cancer (CRC), breast cancer (BC), and prostate cancer (PCA) in the Japanese population. We also identified pathogenic variants in other cancer risk genes, giving consideration to future multigene testing panels for this population.

METHODS

We performed whole-genome sequencing for 1,037 Japanese individuals, including patients with early-onset CRC (n = 196), BC (n = 237), and PCA (n = 215) and controls (n = 389). We screened for pathogenic variants, including single nucleotide variants and copy number variants, among well-established first-tier cancer genes for each cancer type and examined an expended second-tier panel including cancer-predisposing genes from the Cancer Gene Census.

RESULTS

Proportions of patients with germline pathogenic variants differed by cancer subgroup, with the highest in BC (14.8%), followed by CRC (9.2%), and PCA (3.7%). In contrast, 2 of 389 control subjects (0.5%) carried a germline pathogenic variant. In comparison with controls, the proportion of patients with pathogenic variants in the second-tier panel was increased significantly for PCA (3.7% to 11.6%, P = 2.96 × 10^-4), but not for CRC or BC, after multitesting adjustment. In patients with PCA, DNA repair pathway genes in the extended panel often contained pathogenic variants (P = .011).

CONCLUSION

Our analyses support the clinical usefulness of established cancer gene panels in the Japanese population for 3 major cancer types. Additional genes, especially those involved in DNA repair, might be considered for developing multipanel testing in Japanese patients with early-onset PCA.

Self-Peptidome Variation Shapes Individual Immune Responses

The relationship between human genetic variation and disease has not been fully elucidated. According to the present view on infectious diseases pathogen resistance is linked to human leukocyte antigen (HLA) class I/II variants and their individual capacity to present pathogen-derived peptides. Yet, T cell education in the thymus occurs through negative and positive selection, and both processes are controlled by a combination of HLA class I/II variants and peptides from the self. Therefore, the capacity of given HLA class I/II variants to bind pathogen-derived peptides is only one part of the selective process to generate effective immune responses. We thus propose that peptidome variation contributes to shaping T cell receptor (TCR) repertoires and hence individual immune responses, and that this variation represents inherent modulator epitopes.

TCR repertoires emerge in the thymus in each individual as T cells undergo positive and negative selection.

T cell education is controlled by the combination of HLA class I/II molecules and their peptide pools (peptidome).

HLA class I/II molecules are highly plastic in human populations but the peptidome is also a source of variation. Hence combined diversity of HLA class I/II molecules and of self-peptides shapes individual immune responses.

Self-peptide variants that affect T cell repertoires represent inherent modulator epitopes.

Association of the RPA3-UMAD1 locus with interstitial lung diseases complicated with rheumatoid arthritis in Japanese

OBJECTIVES

The genetic background of rheumatoid arthritis-interstitial lung disease (RA-ILD) has been evaluated in Europeans, but little knowledge has been obtained in non-Europeans. This study aimed to elucidate genome-wide risk of RA-ILD in non-Europeans.

METHODS

We performed an initial genome-wide association study (GWAS) of RA-ILD in the Japanese population. By conducting the meta-analysis of the three GWAS datasets of the RA cohorts and biobank of Japanese, our study included 358 RA-ILD cases and 4550 RA subjects without ILD. We then conducted the stratified analysis of the effect of the GWAS risk allele in each CT image pattern.

RESULTS

We identified one novel RA-ILD risk locus at 7p21 that satisfied the genome-wide significance threshold (rs12702634 at RPA3-UMAD1, OR=2.04, 95% CI 1.59 to 2.60, p=1.5×10-8). Subsequent stratified analysis based on the CT image patterns demonstrated that the effect size of the RA-ILD risk allele (rs12702634-C) was large with the UIP pattern (OR=1.86, 95% CI 0.97 to 3.58, p=0.062) and the probable UIP pattern (OR=2.26, 95% CI 1.36 to 3.73, p=0.0015).

CONCLUSION

We revealed one novel genetic association with RA-ILD in Japanese. The RA-ILD risk of the identified variant at RPA3-UMAD1 was relatively high in the CT image patterns related to fibrosis. Our study should contribute to elucidation of the complicated aetiology of RA-ILD.

A genome-wide association study on fish consumption in a Japanese population-the Japan Multi-Institutional Collaborative Cohort study

BACKGROUND/OBJECTIVE

Although benefits of fish consumption for health are well known, a significant percentage of individuals dislike eating fish. Fish consumption may be influenced by genetic factors in addition to environmental factors. We conducted a genome-wide association study (GWAS) to find genetic variations that affect fish consumption in a Japanese population.

METHODS

We performed a two-stage GWAS on fish consumption using 13,739 discovery samples from the Japan Multi-Institutional Collaborative Cohort study, and 2845 replication samples from the other population. We used a semi-quantitative food frequency questionnaire to estimate food intake. Association of the imputed variants with fish consumption was analyzed by separate linear regression models per variant, with adjustments for age, sex, energy intake, principal component analysis components 1-10, and alcohol intake (g/day). We also performed conditional analysis.

RESULTS

We found 27 single nucleotide polymorphisms (SNPs) located in 12q24 and 14q32.12 that were associated with fish consumption. The 19 SNPs were located at 11 genes including six lead SNPs at the BRAP, ACAD10, ALDH2, NAA25, and HECTD4 regions on 12q24.12-13, and CCDC197 region on 14q32.12. In replication samples, all five SNPs located on chromosome 12 were replicated successfully, but the one on chromosome 14 was not. Conditional analyses revealed that the five lead variants in chromosome 12 were in fact the same signal.

CONCLUSION

We found that new SNPs in the 12q24 locus were related to fish intake in two Japanese populations. The associations between SNPs on chromosome 12 and fish intake were strongly confounded by drinking status.

Association of single nucleotide polymorphisms in the NRF2 promoter with vascular stiffness with aging

Purpose

Pulse wave velocity (PWV), an indicator of vascular stiffness, increases with age and is increasingly recognized as an independent risk factor for cardiovascular disease (CVD). Although many mechanical and chemical factors underlie the stiffness of the elastic artery, genetic risk factors related to age-dependent increases in PWV in apparently healthy people are largely unknown. The transcription factor nuclear factor E2 (NF-E2)-related factor 2 (Nrf2), which is activated by unidirectional vascular pulsatile shear stress or oxidative stress, regulates vascular redox homeostasis. Previous reports have shown that a SNP in the NRF2 gene regulatory region (−617C>A; hereafter called SNP−617) affects NRF2 gene expression such that the minor A allele confers lower gene expression compared to the C allele, and it is associated with various diseases, including CVD. We aimed to investigate whether SNP−617 affects vascular stiffness with aging in apparently healthy people.

Methods

Analyzing wide-ranging data obtained from a public health survey performed in Japan, we evaluated whether SNP−617 affected brachial-ankle PWV (baPWV) in never-smoking healthy subjects (n = 642). We also evaluated the effects of SNP−617 on other cardiovascular and blood test measurements.

Results

We have shown that not only AA carriers (n = 55) but also CA carriers (n = 247) show arterial stiffness compared to CC carriers (n = 340). Furthermore, SNP−617 also affected blood pressure indexes such as systolic blood pressure and mean arterial pressure but not the ankle brachial pressure index, an indicator of atherosclerosis. Multivariate analysis showed that SNP−617 accelerates the incremental ratio of baPWV with age.

Conclusions

This study is the first to show that SNP−617 affects the age-dependent increase in vascular stiffness. Our results indicate that low NRF2 activity induces premature vascular aging and could be targeted for the prevention of cardiovascular diseases associated with aging.

Endogenization and excision of human herpesvirus 6 in human genomes

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.

Differences in local population history at the finest level: the case of the Estonian population

Several recent studies detected fine-scale genetic structure in human populations. Hence, groups conventionally treated as single populations harbour significant variation in terms of allele frequencies and patterns of haplotype sharing. It has been shown that these findings should be considered when performing studies of genetic associations and natural selection, especially when dealing with polygenic phenotypes. However, there is little understanding of the practical effects of such genetic structure on demography reconstructions and selection scans when focusing on recent population history. Here we tested the impact of population structure on such inferences using high-coverage (~30×) genome sequences of 2305 Estonians. We show that different regions of Estonia differ in both effective population size dynamics and signatures of natural selection. By analyzing identity-by-descent segments we also reveal that some Estonian regions exhibit evidence of a bottleneck 10-15 generations ago reflecting sequential episodes of wars, plague and famine, although this signal is virtually undetected when treating Estonia as a single population. Besides that, we provide a framework for relating effective population size estimated from genetic data to actual census size and validate it on the Estonian population. This approach may be widely used both to cross-check estimates based on historical sources as well as to get insight into times and/or regions with no other information available. Our results suggest that the history of human populations within the last few millennia can be highly region specific and cannot be properly studied without taking local genetic structure into account.

PLA2G6 variants associated with the number of affected alleles in Parkinson's disease in Japan

This study aimed to evaluate genotype-phenotype correlations of Parkinson's disease (PD) patients with phospholipase A2 group V (PLA2G6) variants. We analyzed the DNA of 798 patients with PD, including 78 PD patients reported previously, and 336 in-house controls. We screened the exons and exon-intron boundaries of PLA2G6 using the Ion Torrent system and Sanger method. We identified 21 patients with 18 rare variants, such that 1, 9, and 11 patients were homozygous, heterozygous, and compound heterozygous, respectively, with respect to PLA2G6 variants. The allele frequency was approximately equal between patients with familial PD and those with sporadic PD. The PLA2G6 variants detected frequently were identified in the early-onset sporadic PD group. Patients who were homozygous for a variant showed more severe symptoms than those who were heterozygous for the variant. The most common variant was p.R635Q in our cohort, which was considered a risk variant for PD. Thus, the variants of PLA2G6 may play a role in familial PD and early-onset sporadic PD.

Evidence of Polygenic Adaptation in Sardinia at Height-Associated Loci Ascertained from the Biobank Japan

Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

The identified clinical features of Parkinson's disease in homo-, heterozygous and digenic variants of PINK1

To investigate the prevalence and genotype-phenotype correlations of phosphatase and tensin homolog induced putative kinase 1 (PINK1) variants in Parkinson's disease (PD) patients, we analyzed 1700 patients (842 familial PD and 858 sporadic PD patients from Japanese origin). We screened the entire exon and exon-intron boundaries of PINK1 using Sanger sequencing and target sequencing by Ion torrent system. We identified 30 patients with heterozygous variants, 3 with homozygous variants, and 3 with digenic variants of PINK1-PRKN. Patients with homozygous variants presented a significantly younger age at onset than those with heterozygous variants. The allele frequency of heterozygous variants in patients with age at onset at 50 years and younger with familial PD and sporadic PD showed no differences. [¹²³I]meta-iodobenzylguanidine (MIBG) myocardial scintigraphy indicated that half of patients harboring PINK1 heterozygous variants showed a decreased heart to mediastinum ratio (12/23). Our findings emphasize the importance of PINK1 variants for the onset of PD in patients with age at onset at 50 years and younger and the broad spectrum of clinical symptoms in patients with PINK1 variants.

Chromosomal alterations among age-related haematopoietic clones in Japan

The extent to which the biology of oncogenesis and ageing are shaped by factors that distinguish human populations is unknown. Haematopoietic clones with acquired mutations become common with advancing age and can lead to blood cancers^1-10. Here we describe shared and population-specific patterns of genomic mutations and clonal selection in haematopoietic cells on the basis of 33,250 autosomal mosaic chromosomal alterations that we detected in 179,417 Japanese participants in the BioBank Japan cohort and compared with analogous data from the UK Biobank. In this long-lived Japanese population, mosaic chromosomal alterations were detected in more than 35.0% (s.e.m., 1.4%) of individuals older than 90 years, which suggests that such clones trend towards inevitability with advancing age. Japanese and European individuals exhibited key differences in the genomic locations of mutations in their respective haematopoietic clones; these differences predicted the relative rates of chronic lymphocytic leukaemia (which is more common among European individuals) and T cell leukaemia (which is more common among Japanese individuals) in these populations. Three different mutational precursors of chronic lymphocytic leukaemia (including trisomy 12, loss of chromosomes 13q and 13q, and copy-neutral loss of heterozygosity) were between two and six times less common among Japanese individuals, which suggests that the Japanese and European populations differ in selective pressures on clones long before the development of clinically apparent chronic lymphocytic leukaemia. Japanese and British populations also exhibited very different rates of clones that arose from B and T cell lineages, which predicted the relative rates of B and T cell cancers in these populations. We identified six previously undescribed loci at which inherited variants predispose to mosaic chromosomal alterations that duplicate or remove the inherited risk alleles, including large-effect rare variants at NBN, MRE11 and CTU2 (odds ratio, 28-91). We suggest that selective pressures on clones are modulated by factors that are specific to human populations. Further genomic characterization of clonal selection and cancer in populations from around the world is therefore warranted.