A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk.

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Helicobacter pylori, Homologous-Recombination Genes, and Gastric Cancer.

BACKGROUND: Helicobacter pylori infection is a well-known risk factor for gastric cancer. However, the contribution of germline pathogenic variants in cancer-predisposing genes and their effect, when combined with H. pylori infection, on the risk of gastric cancer has not been widely evaluated. METHODS: We evaluated the association between germline pathogenic variants in 27 cancer-predisposing genes and the risk of gastric cancer in a sample of 10,426 patients with gastric cancer and 38,153 controls from BioBank Japan. We also assessed the combined effect of pathogenic variants and H. pylori infection status on the risk of gastric cancer and calculated the cumulative risk in 1433 patients with gastric cancer and 5997 controls from the Hospital-based Epidemiologic Research Program at Aichi Cancer Center (HERPACC). RESULTS: Germline pathogenic variants in nine genes (APC, ATM, BRCA1, BRCA2, CDH1, MLH1, MSH2, MSH6, and PALB2) were associated with the risk of gastric cancer. We found an interaction between H. pylori infection and pathogenic variants in homologous-recombination genes with respect to the risk of gastric cancer in the sample from HERPACC (relative excess risk due to the interaction, 16.01; 95% confidence interval [CI], 2.22 to 29.81; P = 0.02). At 85 years of age, persons with H. pylori infection and a pathogenic variant had a higher cumulative risk of gastric cancer than noncarriers infected with H. pylori (45.5% [95% CI, 20.7 to 62.6] vs. 14.4% [95% CI, 12.2 to 16.6]). CONCLUSIONS: H. pylori infection modified the risk of gastric cancer associated with germline pathogenic variants in homologous-recombination genes. (Funded by the Japan Agency for Medical Research and Development and others.).

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

Genotype imputation accuracy and the quality metrics of the minor ancestry in multi-ancestry reference panels.

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value.

Germ line DDX41 mutations define a unique subtype of myeloid neoplasms.

Germ line DDX41 variants have been implicated in late-onset myeloid neoplasms (MNs). Despite an increasing number of publications, many important features of DDX41-mutated MNs remain to be elucidated. Here we performed a comprehensive characterization of DDX41-mutated MNs, enrolling a total of 346 patients with DDX41 pathogenic/likely-pathogenic (P/LP) germ line variants and/or somatic mutations from 9082 MN patients, together with 525 first-degree relatives of DDX41-mutated and wild-type (WT) patients. P/LP DDX41 germ line variants explained ∼80% of known germ line predisposition to MNs in adults. These risk variants were 10-fold more enriched in Japanese MN cases (n = 4461) compared with the general population of Japan (n = 20 238). This enrichment of DDX41 risk alleles was much more prominent in male than female (20.7 vs 5.0). P/LP DDX41 variants conferred a large risk of developing MNs, which was negligible until 40 years of age but rapidly increased to 49% by 90 years of age. Patients with myelodysplastic syndromes (MDS) along with a DDX41-mutation rapidly progressed to acute myeloid leukemia (AML), which was however, confined to those having truncating variants. Comutation patterns at diagnosis and at progression to AML were substantially different between DDX41-mutated and WT cases, in which none of the comutations affected clinical outcomes. Even TP53 mutations made no exceptions and their dismal effect, including multihit allelic status, on survival was almost completely mitigated by the presence of DDX41 mutations. Finally, outcomes were not affected by the conventional risk stratifications including the revised/molecular International Prognostic Scoring System. Our findings establish that MDS with DDX41-mutation defines a unique subtype of MNs that is distinct from other MNs.

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

Different risk genes contribute to clear cell and non-clear cell renal cell carcinoma in 1532 Japanese patients and 5996 controls.

Identifying causative genes via genetic testing is useful for screening, preventing and treating cancer. Several hereditary syndromes occur in patients with renal cell carcinoma (RCC). However, the evidence is from the European population; it remains unclear how the RCC-related genes and other cancer-predisposing genes contribute to RCC development in the Japanese population. A case-control study of 14 RCC-related genes and 26 cancer-predisposing genes was performed in 1563 Japanese patients with RCC and 6016 controls. The patients were stratified into clear cell RCC (ccRCC) or non-ccRCC (nccRCC). Gene-based analysis of germline pathogenic variants in patients with each subtype and cancer-free subjects was performed. Following quality control, 1532 patients with RCC and 5996 controls were analyzed. For ccRCC, 52 of 1283 (4.05%) patients carried pathogenic variants mainly in the cancer-predisposing genes such as TP53 (P = 1.73 × 10-4; OR, 5.8; 95% CI, 2.2-15.7). Approximately 80% of patients with pathogenic variants in TP53 had p.Ala189Val that was specific in East Asian population. For nccRCC, 14 of 249 (5.62%) patients carried pathogenic variants mainly in the RCC-related genes such as BAP1 and FH (P = 6.27 × 10-5; OR, Inf; 95% CI, 10.0-Inf). The patients with the pathogenic variants in the associated genes were diagnosed 15.8 years earlier and had a higher proportion of patients with a family history of RCC (OR, 20.0; 95% CI, 1.3-237.4) than the non-carriers. We showed different and population-specific contributions of risk genes between ccRCC and nccRCC in Japanese for improved personalized medicine.

Genomic variants associated with age at diagnosis of childhood-onset type 1 diabetes.

Age at diagnosis (AAD) of Type 1 diabetes (T1D) is determined by the age at onset of the autoimmune attack and by the rate of beta cell destruction that follows. Twin studies found that T1D AAD is strongly influenced by genetics, notably in young children. In young UK, Finnish, Sardinian patients AAD-associated genomic variants were previously identified, which may vary across populations and with time. In 1956 children of European ancestry born in mainland France in 1980-2008 who declared T1D before 15 years, we tested 94 T1D-associated SNPs for their association with AAD using nonparametric Kruskal-Wallis test. While high-risk HLA genotypes were not found to be associated with AAD, fourteen SNPs located in 12 non-HLA loci showed a strong association (2.9 × 10-12 < P < 1.4 × 10-3 after FDR correction). Four of these loci have been associated with AAD in previous cohorts (GSDMB, IL2, TNFAIP3, IL1), supporting a partially shared genetic influence on AAD of T1D in the studied European populations. In contrast, the association of 8 new loci CLEC16A, TYK2, ERBB3, CCR7, FCRL3, DNAH2, FGF3/4, and HPSE2 with AAD is novel. The 12 protein-coding genes located within these loci are involved in major immune pathways or in predisposition to other autoimmune diseases, which suggests a prominent role for these genes in the early immune mechanisms of beta cell destruction.

Identification of epistatic SNP combinations in rheumatoid arthritis using LAMPLINK and Japanese cohorts.

Genome-wide association studies have enabled the identification of important genetic factors in many trait studies. However, only a fraction of the heritability can be explained by known genetic factors, even in the most common diseases. Genetic loci combinations, or epistatic contributions expressed by combinations of single nucleotide polymorphisms (SNPs), have been argued to be one of the critical factors explaining some of the missing heritability, especially in oligogenic/polygenic diseases. Rheumatoid arthritis (RA) is a complex disease with more than 100 reported SNP associations, as well as various HLA haplotypes and amino acids; however, many associations between RA and inter-chromosomal SNP combinations are unknown. To discover novel associations of epistatic interactions with high odds ratios in RA, we applied the LAMPLINK method, a systematic enumerative procedure for identifying high-order SNP combinations, to a Japanese RA cohort (discovery cohort; 4024 patients with RA and 7731 controls). We validated the identified associations in a different Japanese cohort (validation cohort; 810 RA patients and 6303 controls). In this study, we identified 90 significant genetic associations in the discovery cohort. Among these, 74 (82.2%) associations were replicated in the validation cohort, and eight combinations were inter-chromosomal, all of which comprised rs7765379 or rs35265698 located in the HLA region. These two SNPs exhibited strong correlations with valine at amino acid position 11 in HLA-DRB1 (HLA-DRB1-11-Val). Finally, we discovered that rs9624 showed an association with RA through an epistatic interaction with HLA-DRB1-11-Val. Overall, LAMPLINK showed high reliability for identifying epistatic genetic contributions hidden in complex traits.

Genome-Wide Association Study with Three Control Cohorts of Japanese Patients with Esotropia and Exotropia of Comitant Strabismus and Idiopathic Superior Oblique Muscle Palsy.

Esotropia and exotropia in the entity of comitant strabismus are multifactorial diseases with both genetic and environmental backgrounds. Idiopathic superior oblique muscle palsy, as the predominant entity of non-comitant (paralytic) strabismus, also has a genetic background, as evidenced by varying degrees of muscle hypoplasia. A genome-wide association study (GWAS) was conducted of 711 Japanese patients with esotropia (n= 253), exotropia (n = 356), and idiopathic superior oblique muscle palsy (n = 102). The genotypes of single nucleotide polymorphisms (SNPs) were determined by Infinium Asian Screening Array. Three control cohorts from the Japanese population were used: two cohorts from BioBank Japan (BBJ) and the Nagahama Cohort. BBJ (180K) was genotyped by a different array, Illumina Infinium OmniExpressExome or HumanOmniExpress, while BBJ (ASA) and the Nagahama Cohort were genotyped by the same Asian array. After quality control of SNPs and individuals, common SNPs between the case cohort and the control cohort were chosen in the condition of genotyping by different arrays, while all SNPs genotyped by the same array were used for SNP imputation. The SNPs imputed with R-square values ≥ 0.3 were used to compare the case cohort of each entity or the combined entity with the control cohort. In comparison with BBJ (180K), the esotropia group and the exotropia group showed CDCA7 and HLA-F, respectively, as candidate genes at a significant level of p < 5 × 10-8, while the idiopathic superior oblique muscle palsy group showed DAB1 as a candidate gene which is involved in neuronal migration. DAB1 was also detected as a candidate in comparison with BBJ (ASA) and the Nagahama Cohort at a weak level of significance of p < 1 × 10-6. In comparison with BBJ (180K), RARB (retinoic acid receptor-ß) was detected as a candidate at a significant level of p < 5 × 10-8 in the combined group of esotropia, exotropia, and idiopathic superior oblique muscle palsy. In conclusion, a series of GWASs with three different control cohorts would be an effective method with which to search for candidate genes for multifactorial diseases such as strabismus.

Genetic Risk Score for Intracranial Aneurysms: Prediction of Subarachnoid Hemorrhage and Role in Clinical Heterogeneity.

BACKGROUND: Recently, common genetic risk factors for intracranial aneurysm (IA) and aneurysmal subarachnoid hemorrhage (ASAH) were found to explain a large amount of disease heritability and therefore have potential to be used for genetic risk prediction. We constructed a genetic risk score to (1) predict ASAH incidence and IA presence (combined set of unruptured IA and ASAH) and (2) assess its association with patient characteristics. METHODS: A genetic risk score incorporating genetic association data for IA and 17 traits related to IA (so-called metaGRS) was created using 1161 IA cases and 407 392 controls from the UK Biobank population study. The metaGRS was validated in combination with risk factors blood pressure, sex, and smoking in 828 IA cases and 68 568 controls from the Nordic HUNT population study. Furthermore, we assessed association between the metaGRS and patient characteristics in a cohort of 5560 IA patients. RESULTS: Per SD increase of metaGRS, the hazard ratio for ASAH incidence was 1.34 (95% CI, 1.20-1.51) and the odds ratio for IA presence 1.09 (95% CI, 1.01-1.18). Upon including the metaGRS on top of clinical risk factors, the concordance index to predict ASAH hazard increased from 0.63 (95% CI, 0.59-0.67) to 0.65 (95% CI, 0.62-0.69), while prediction of IA presence did not improve. The metaGRS was statistically significantly associated with age at ASAH (ß=-4.82×10-3 per year [95% CI, -6.49×10-3 to -3.14×10-3]; P=1.82×10-8), and location of IA at the internal carotid artery (odds ratio=0.92 [95% CI, 0.86-0.98]; P=0.0041). CONCLUSIONS: The metaGRS was predictive of ASAH incidence, although with limited added value over clinical risk factors. The metaGRS was not predictive of IA presence. Therefore, we do not recommend using this metaGRS in daily clinical care. Genetic risk does partly explain the clinical heterogeneity of IA warranting prioritization of clinical heterogeneity in future genetic prediction studies of IA and ASAH.

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.

Genome-Wide Association Study of Age-Related Macular Degeneration Reveals 2 New Loci Implying Shared Genetic Components with Central Serous Chorioretinopathy.

PURPOSE: To investigate the genetic architecture of age-related macular degeneration (AMD) in a Japanese population. DESIGN: Genome-wide association study (GWAS). PARTICIPANTS: Three thousand seven hundred seventy-two patients with AMD and 16 770 control participants from the Japanese population were enrolled in the association analyses. METHODS: We conducted a meta-analysis of 2 independent GWASs that included a total of 2663 patients with AMD and 9471 control participants using the imputation reference panel for genotype imputation specified for the Japanese population (n = 3541). A replication study was performed using an independent set of 1109 patients with AMD and 7299 control participants. MAIN OUTCOME MEASURES: Associations of genetic variants with AMD. RESULTS: A meta-analysis of the 2 GWASs identified 6 loci significantly associated with AMD (P < 5.0 × 10-8). Of these loci, 4 were known to be associated with AMD (CFH, C2/FB, TNFRSF10A, and ARMS2), and 2 were novel (rs4147157 near WBP1L and rs76228488 near GATA5). The newly identified associations were confirmed in a replication study (P < 0.01). After the meta-analysis of all datasets, we observed strong associations in these loci (P = 1.88 × 10-12 and P = 1.35 × 10-9 for meta-analysis for rs4147157 and rs76228488, respectively). When we looked up the associations in the reported central serous chorioretinopathy (CSC) GWAS conducted in the Japanese population, both loci were associated significantly with CSC (P = 4.86 × 10-3 and P = 4.28 × 10-3 for rs4147157 and rs76228488, respectively). We performed a genetic colocalization analysis for these loci and estimated that the posterior probabilities of shared causal variants between AMD and CSC were 0.39 and 0.60 for WBP1L and GATA5, respectively. Genetic correlation analysis focusing on the epidemiologically suggested clinical risk factors implicated shared polygenic architecture between AMD and smoking cessation (rg [the measure of genetic correlation] = -0.33; P = 0.01; false discovery rate, 0.099). CONCLUSIONS: Our findings imply shared genetic components conferring the risk of both AMD and CSC. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

Genetic risk score of cerebral infarction in atrial fibrillation genome-wide association study.

INTRODUCTION: Stroke is a leading cause of death and the primary cause of adult-acquired disability. Patients with cardiogenic embolic stroke also have higher mortality and recurrence rates than patients with other stroke subtypes. Atrial fibrillation (AF) is a major risk factor for cerebral infarction (CI). The large-scale study identified 32 loci in the MEGASTROKE study. However, few studies have attempted to identify novel stroke risk variants in patients with a history of AF. Our overall aim was to identify novel CI risk variants in AF cases and explore whether their associations with the CI risk were affected by the CHADS2 and CHA2DS2-VASc scores. METHODS: We performed association study with CI using 8181 AF cases in previous genome-wide association study (GWAS) and imputation data without controls. We classified AF cases into those with or without past history of CI, and the genetic associations with the CI risk were examined. RESULTS: GWAS identified eight associated loci. The generated genetic risk score (GRS) for the eight loci was significantly associated with CI in patients with AF (1.46 × 10-8 ). We estimated bivariate logistic regression model which contained GRS and CHADS2 score (GRS: p-Value = 7.41 × 10-9 , CHADS2 score: p-Value <2.0 × 10-16 ) or CHA2DS2-VASc scores (GRS: p-Value = 2.52 × 10-10 , CHA2DS2-VASc score: p-Value <2.0 × 10-16 ). CONCLUSION: We identified eight genetic variants that were potentially associated with the risk of CI of AF cases and the significant GRS, whose associations were independent of the CHADS2 or CHA2DS2-VASc score.

Variants encoding a restricted carboxy-terminal domain of SLC12A2 cause hereditary hearing loss in humans.

Hereditary hearing loss is challenging to diagnose because of the heterogeneity of the causative genes. Further, some genes involved in hereditary hearing loss have yet to be identified. Using whole-exome analysis of three families with congenital, severe-to-profound hearing loss, we identified a missense variant of SLC12A2 in five affected members of one family showing a dominant inheritance mode, along with de novo splice-site and missense variants of SLC12A2 in two sporadic cases, as promising candidates associated with hearing loss. Furthermore, we detected another de novo missense variant of SLC12A2 in a sporadic case. SLC12A2 encodes Na+, K+, 2Cl- cotransporter (NKCC) 1 and plays critical roles in the homeostasis of K+-enriched endolymph. Slc12a2-deficient mice have congenital, profound deafness; however, no human variant of SLC12A2 has been reported as associated with hearing loss. All identified SLC12A2 variants mapped to exon 21 or its 3'-splice site. In vitro analysis indicated that the splice-site variant generates an exon 21-skipped SLC12A2 mRNA transcript expressed at much lower levels than the exon 21-included transcript in the cochlea, suggesting a tissue-specific role for the exon 21-encoded region in the carboy-terminal domain. In vitro functional analysis demonstrated that Cl- influx was significantly decreased in all SLC12A2 variants studied. Immunohistochemistry revealed that SLC12A2 is located on the plasma membrane of several types of cells in the cochlea, including the strial marginal cells, which are critical for endolymph homeostasis. Overall, this study suggests that variants affecting exon 21 of the SLC12A2 transcript are responsible for hereditary hearing loss in humans.

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.

A polygenic risk score improves risk stratification of coronary artery disease: a large-scale prospective Chinese cohort study.

AIMS: To construct a polygenic risk score (PRS) for coronary artery disease (CAD) and comprehensively evaluate its potential in clinical utility for primary prevention in Chinese populations. METHODS AND RESULTS: Using meta-analytic approach and large genome-wide association results for CAD and CAD-related traits in East Asians, a PRS comprising 540 genetic variants was developed in a training set of 2800 patients with CAD and 2055 controls, and was further assessed for risk stratification for CAD integrating with the guideline-recommended clinical risk score in large prospective cohorts comprising 41 271 individuals. During a mean follow-up of 13.0 years, 1303 incident CAD cases were identified. Individuals with high PRS (the highest 20%) had about three-fold higher risk of CAD than the lowest 20% (hazard ratio 2.91, 95% confidence interval 2.43-3.49), with the lifetime risk of 15.9 and 5.8%, respectively. The addition of PRS to the clinical risk score yielded a modest yet significant improvement in C-statistic (1%) and net reclassification improvement (3.5%). We observed significant gradients in both 10-year and lifetime risk of CAD according to the PRS within each clinical risk strata. Particularly, when integrating high PRS, intermediate clinical risk individuals with uncertain clinical decision for intervention would reach the risk levels (10-year of 4.6 vs. 4.8%, lifetime of 17.9 vs. 16.6%) of high clinical risk individuals with intermediate (20-80%) PRS. CONCLUSION: The PRS could stratify individuals into different trajectories of CAD risk, and further refine risk stratification for CAD within each clinical risk strata, demonstrating a great potential to identify high-risk individuals for targeted intervention in clinical utility.

Association between germline pathogenic variants in cancer-predisposing genes and lymphoma risk.

The application of advanced molecular technology has significantly expanded lymphoma classification, allowing risk stratification and treatment optimization. Limited evidence suggests the presence of a genetic predisposition in lymphoma, indicating the potential for better individualized clinical management based on a novel lymphoma classification. Herein, we examined the impact of germline pathogenic variants in 27 cancer-predisposing genes with lymphoma risk and explored the clinical characteristics of pathogenic variant carriers. This study included 2,066 lymphoma patients and 38,153 cancer-free controls from the Japanese population. Following quality control of sequencing data, samples from 1,982 lymphoma patients and 37,592 controls were further analyzed. We identified 309 pathogenic variants among 4,850 variants in the 27 cancer-predisposing genes. Pathogenic variants in the following four cancer-predisposing genes were associated with a high risk of lymphoma: ATM (odds ratio [OR], 2.63; 95% confidence interval [CI], 1.25-5.51; p = 1.06 × 10-2 ), BRCA1 (OR, 5.88; 95% CI, 2.65-13.02; p = 1.27 × 10-5 ), BRCA2 (OR, 2.94; 95% CI, 1.60-5.42; p = 5.25 × 10-4 ), and TP53 (OR, 5.22; 95% CI, 1.43-19.02; p = 1.23 × 10-2 ). The proportion of carriers of these genes was 1.6% of lymphoma patients. Furthermore, pathogenic variants in these genes were especially associated with a higher risk of mantle cell lymphoma (OR, 21.57; 95% CI, 7.59-61.26; p = 8.07 × 10-9 ). These results provide novel insights concerning monogenic form into lymphoma classification. Some lymphoma patients may benefit from surveillance and targeted treatment, such as other neoplasms.

Population-based Screening for Hereditary Colorectal Cancer Variants in Japan.

BACKGROUND & AIMS: Colorectal cancer (CRC) is one of the most common cancers in the world. A small proportion of CRCs can be attributed to recognizable hereditary germline variants of known CRC susceptibility genes. To better understand cancer risk, it is necessary to explore the prevalence of hereditary CRC and pathogenic variants of multiple cancer-predisposing genes in non-European populations. METHODS: We analyzed the coding regions of 27 cancer-predisposing genes in 12,503 unselected Japanese CRC patients and 23,705 controls by target sequencing and genome-wide SNP chip. Their clinical significance was assessed using ClinVar and the guidelines by ACMG/AMP. RESULTS: We identified 4,804 variants in the 27 genes and annotated them as pathogenic in 397 and benign variants in 941, of which 43.6% were novel. In total, 3.3% of the unselected CRC patients and 1.5% of the controls had a pathogenic variant. The pathogenic variants of MSH2 (odds ratio (OR) = 18.1), MLH1 (OR = 8.6), MSH6 (OR = 4.9), APC (OR = 49.4), BRIP1 (OR=3.6), BRCA1 (OR = 2.6), BRCA2 (OR = 1.9), and TP53 (OR = 1.7) were significantly associated with CRC development in the Japanese population (P-values<0.01, FDR<0.05). These pathogenic variants were significantly associated with diagnosis age and personal/family history of cancer. In total, at least 3.5% of the Japanese CRC population had a pathogenic variant or CNV of the 27 cancer-predisposing genes, indicating hereditary cancers. CONCLUSIONS: This largest study of CRC heredity in Asia can contribute to the development of guidelines for genetic testing and variant interpretation for heritable CRCs.

Multi-trait and cross-population genome-wide association studies across autoimmune and allergic diseases identify shared and distinct genetic component.

OBJECTIVES: Autoimmune and allergic diseases are outcomes of the dysregulation of the immune system. Our study aimed to elucidate differences or shared components in genetic backgrounds between autoimmune and allergic diseases. METHODS: We estimated genetic correlation and performed multi-trait and cross-population genome-wide association study (GWAS) meta-analysis of six immune-related diseases: rheumatoid arthritis, Graves' disease, type 1 diabetes for autoimmune diseases and asthma, atopic dermatitis and pollinosis for allergic diseases. By integrating large-scale biobank resources (Biobank Japan and UK biobank), our study included 105 721 cases and 433 663 controls. Newly identified variants were evaluated in 21 778 cases and 712 767 controls for two additional autoimmune diseases: psoriasis and systemic lupus erythematosus. We performed enrichment analyses of cell types and biological pathways to highlight shared and distinct perspectives. RESULTS: Autoimmune and allergic diseases were not only mutually classified based on genetic backgrounds but also they had multiple positive genetic correlations beyond the classifications. Multi-trait GWAS meta-analysis newly identified six allergic disease-associated loci. We identified four loci shared between the six autoimmune and allergic diseases (rs10803431 at PRDM2, OR=1.07, p=2.3×10-8, rs2053062 at G3BP1, OR=0.90, p=2.9×10-8, rs2210366 at HBS1L, OR=1.07, p=2.5×10-8 in Japanese and rs4529910 at POU2AF1, OR=0.96, p=1.9×10-10 across ancestries). Associations of rs10803431 and rs4529910 were confirmed at the two additional autoimmune diseases. Enrichment analysis demonstrated link to T cells, natural killer cells and various cytokine signals, including innate immune pathways. CONCLUSION: Our multi-trait and cross-population study should elucidate complex pathogenesis shared components across autoimmune and allergic diseases.