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.

Genome- and transcriptome-wide association studies of 386,000 Asian and European-ancestry women provide new insights into breast cancer genetics.

By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p < 5.0 × 10-8 and a Bonferroni-corrected p < 4.6 × 10-6, respectively. Of them, 32 loci and 15 genes showed a significantly different association between ER-positive and ER-negative breast cancer after Bonferroni correction. Significant ancestral differences in risk variant allele frequencies and their association strengths with breast cancer risk were identified. Of the significant associations identified in this study, 17 loci and 14 genes are located 1Mb away from any of the previously reported breast cancer risk variants. Pathways analyses including 221 putative risk genes identified multiple signaling pathways that may play a significant role in the development of breast cancer. Our study provides a comprehensive understanding of and new biological insights into the genetics of this common malignancy.

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.

Pharmacogenetic and clinical risk factors for bevacizumab-related gastrointestinal hemorrhage in prostate cancer patients treated on CALGB 90401 (Alliance).

The objective of this study was to discover clinical and pharmacogenetic factors associated with bevacizumab-related gastrointestinal hemorrhage in Cancer and Leukemia Group B (Alliance) 90401. Patients with metastatic castration-resistant prostate cancer received docetaxel and prednisone ± bevacizumab. Patients were genotyped using Illumina HumanHap610-Quad and assessed using cause-specific risk for association between single nucleotide polymorphisms (SNPs) and gastrointestinal hemorrhage. In 1008 patients, grade 2 or higher gastrointestinal hemorrhage occurred in 9.5% and 3.8% of bevacizumab (n = 503) and placebo (n = 505) treated patients, respectively. Bevacizumab (P < 0.001) and age (P = 0.002) were associated with gastrointestinal hemorrhage. In 616 genetically estimated Europeans (n = 314 bevacizumab and n = 302 placebo treated patients), grade 2 or higher gastrointestinal hemorrhage occurred in 9.6% and 2.0% of patients, respectively. One SNP (rs1478947; HR 6.26; 95% CI 3.19-12.28; P = 9.40 × 10-8) surpassed Bonferroni-corrected significance. Grade 2 or higher gastrointestinal hemorrhage rate was 33.3% and 6.2% in bevacizumab-treated patients with the AA/AG and GG genotypes, versus 2.9% and 1.9% in the placebo arm, respectively. Prospective validation of these findings and functional analyses are needed to better understand the genetic contribution to treatment-related gastrointestinal hemorrhage.

ABO O blood group as a risk factor for platelet reactivity in heparin-induced thrombocytopenia.

Heparin-induced thrombocytopenia (HIT) is an unpredictable, potentially catastrophic adverse effect resulting from an immune response to platelet factor 4 (PF4)/heparin complexes. We performed a genome-wide association study (GWAS) with positive functional assay as the outcome in a large discovery cohort of patients divided into 3 groups: (1) functional assay-positive cases (n = 1269), (2) antibody-positive (functional assay-negative) controls (n = 1131), and (3) antibody-negative controls (n = 1766). Significant associations (α = 5 × 10-8) were investigated in a replication cohort (α = 0.05) of functional assay-confirmed HIT cases (n = 177), antibody-positive (function assay-negative) controls (n = 258), and antibody-negative controls (n = 351). We observed a strong association for positive functional assay with increasing PF4/heparin immunoglobulin-G (IgG) level (odds ratio [OR], 16.53; 95% confidence interval [CI], 13.83-19.74; P = 1.51 × 10-209) and female sex (OR, 1.15; 95% CI, 1.01-1.32; P = .034). The rs8176719 C insertion variant in ABO was significantly associated with positive functional assay status in the discovery cohort (frequency = 0.41; OR, 0.751; 95% CI, 0.682-0.828; P = 7.80 × 10-9) and in the replication cohort (OR, 0.467; 95% CI, 0.228-0.954; P = .0367). The rs8176719 C insertion, which encodes all non-O blood group alleles, had a protective effect, indicating that the rs8176719 C deletion and the O blood group were risk factors for HIT (O blood group OR, 1.42; 95% CI, 1.26-1.61; P = 3.09 × 10-8). Meta-analyses indicated that the ABO association was independent of PF4/heparin IgG levels and was stronger when functional assay-positive cases were compared with antibody-positive (functional assay-negative) controls than with antibody-negative controls. Sequencing and fine-mapping of ABO demonstrated that rs8176719 was the causal single nucleotide polymorphism (SNP). Our results clarify the biology underlying HIT pathogenesis with ramifications for prediction and may have important implications for related conditions, such as vaccine-induced thrombotic thrombocytopenia.

Correction: Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

[This corrects the article DOI: 10.1371/journal.pgen.1009497.].

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

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.

Population-Based Impact of Smoking, Drinking, and Genetic Factors on HDL-cholesterol Levels in J-MICC Study Participants.

BACKGROUND: Environmental and genetic factors are suggested to exhibit factor-based association with HDL-cholesterol (HDL-C) levels. However, the population-based effects of environmental and genetic factors have not been compared clearly. We conducted a cross-sectional study using data from the Japan Multi-Institutional Collaborative Cohort (J-MICC) Study to evaluate the population-based impact of smoking, drinking, and genetic factors on low HDL-C. METHODS: Data from 11,498 men and women aged 35-69 years were collected for a genome-wide association study (GWAS). Sixty-five HDL-C-related SNPs with genome-wide significance (P < 5 × 10-8) were selected from the GWAS catalog, of which seven representative SNPs were defined, and the population-based impact was estimated using population attributable fraction (PAF). RESULTS: We found that smoking, drinking, daily activity, habitual exercise, egg intake, BMI, age, sex, and the SNPs CETP rs3764261, APOA5 rs662799, LIPC rs1800588, LPL rs328, ABCA1 rs2575876, LIPG rs3786247, and APOE rs429358 were associated with HDL-C levels. The gene-environmental interactions on smoking and drinking were not statistically significant. The PAF for low HDL-C was the highest in men (63.2%) and in rs3764261 (31.5%) of the genetic factors, and the PAFs of smoking and drinking were 23.1% and 41.8%, respectively. CONCLUSION: The present study showed that the population-based impact of genomic factor CETP rs3764261 for low HDL-C was higher than that of smoking and lower than that of drinking.

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.

Incorporating European GWAS findings improve polygenic risk prediction accuracy of breast cancer among East Asians.

Previous genome-wide association studies (GWASs) have been largely focused on European (EUR) populations. However, polygenic risk scores (PRSs) derived from EUR have been shown to perform worse in non-EURs compared with EURs. In this study, we aim to improve PRS prediction in East Asians (EASs). We introduce a rescaled meta-analysis framework to combine both EUR (N = 122,175) and EAS (N = 30,801) GWAS summary statistics. To improve PRS prediction in EASs, we use a scaling factor to up-weight the EAS data, such that the resulting effect size estimates are more relevant to EASs. We then derive PRSs for EAS from the rescaled meta-analysis results of EAS and EUR data. Evaluated in an independent EAS validation data set, this approach increases the prediction liability-adjusted Nagelkerke's pseudo R2 by 40%, 41%, and 5%, respectively, compared with PRSs derived from an EAS GWAS only, EUR GWAS only, and conventional fixed-effects meta-analysis of EAS and EUR data. The PRS derived from the rescaled meta-analysis approach achieved an area under the receiver operating characteristic curve (AUC) of 0.6059, higher than AUC = 0.5782, 0.5809, 0.6008 for EAS, EUR, and conventional meta-analysis of EAS and EUR. We further compare PRSs constructed by single-nucleotide polymorphisms that have different linkage disequilibrium (LD) scores and minor allele frequencies (MAFs) between EUR and EAS, and observe that lower LD scores or MAF in EAS correspond to poorer PRS performance (AUC = 0.5677, 0.5530, respectively) than higher LD scores or MAF (AUC = 0.589, 0.5993, respectively). We finally build a PRS stratified by LD score differences in EUR and EAS using rescaled meta-analysis, and obtain an AUC of 0.6096, with improvement over other strategies investigated.

Two polymorphic gene loci associated with treprostinil dose in pulmonary arterial hypertension.

OBJECTIVE: Prostacyclin infusion for pulmonary arterial hypertension (PAH) is an effective therapy with varied dosing requirements and clinical response. The major aim of this study was to determine new biologically-based predictors of prostacyclin treatment response heterogeneity. METHODS: Ninety-eight patients with hemodynamically defined PAH at two academic medical centers volunteered for registry studies. A stable dose of treprostinil was the quantitative phenotype for the genome-wide association study (GWAS). Candidate genes with the largest effect sizes and strongest statistical associations were further characterized with in silico and in-vitro assays to confirm mechanistic hypotheses. The clinical significance of these candidate predictors was assessed for mechanistically consistent physiologic effects in an independent cohort of patients. RESULTS: GWAS identified three loci for association with P < 10-6. All three loci had clinically significant effect sizes. Specific single-nucleotide polymorphisms (SNPs) at two of the loci: rs11078738 in phosphoribosylformylglycinamidine synthase and rs10023113 in CAMK2D encoded sequence changes with clear predicted consequences. Production of the primary mediator of prostacyclin-induced vasodilation, cyclic AMP, was reduced in human cell lines by the missense variant rs11078738 (p.L621P). Located in the promoter of CAMK2D, the allele of rs10023113 associated with a higher treprostinil dose has higher ventricular transcription of CAMK2δ. At initial diagnostic catheterization in a separate cohort of patients, the same allele of rs10023113 was associated with elevated right mean atrial and ventricular diastolic pressures. CONCLUSIONS: The quantitative phenotype of stable treprostinil dose identified two gene loci associated with pharmacodynamic response and right ventricular function in PAH worth further investigation.

Accuracy and applications of sequencing and genotyping approaches for CYP2A6 and homologous genes.

OBJECTIVES: We evaluated multiple genotyping/sequencing approaches in a homologous region of chromosome 19, and investigated associations of two common 3'-UTR CYP2A6 variants with activity in vivo. METHODS: Individuals (n = 1704) of European and African ancestry were phenotyped for the nicotine metabolite ratio (NMR), an index of CYP2A6 activity, and genotyped/sequenced using deep amplicon exon sequencing, SNP array, genotype imputation and targeted capture sequencing. Amplicon exon sequencing was the gold standard to which other methods were compared within-individual for CYP2A6, CYP2A7, CYP2A13, and CYP2B6 exons to identify highly discordant positions. Linear regression models evaluated the association of CYP2A6*1B and rs8192733 genotypes (coded additively) with logNMR. RESULTS: All approaches were ≤2.6% discordant with the gold standard; discordant calls were concentrated at few positions. Fifteen positions were discordant in >10% of individuals, with 12 appearing in regions of high identity between homologous genes (e.g. CYP2A6 and CYP2A7). For six, allele frequencies in our study and online databases were discrepant, suggesting errors in online sources. In the European-ancestry group (n = 935), CYP2A6*1B and rs8192733 were associated with logNMR (P < 0.001). A combined model found main effects of both variants on increasing logNMR. Similar trends were found in those of African ancestry (n = 506). CONCLUSION: Multiple genotyping/sequencing approaches used in this chromosome 19 region contain genotyping/sequencing errors, as do online databases. Gene-specific primers and SNP array probes must consider gene homology; short-read sequencing of related genes in a single reaction should be avoided. Using improved sequencing approaches, we characterized two gain-of-function 3'-UTR variants, including the relatively understudied rs8192733.

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.

Assessing the Relationship Between High-sensitivity C-reactive Protein and Kidney Function Employing Mendelian Randomization in the Japanese Community-based J-MICC Study.

BACKGROUND: Inflammation is thought to be a risk factor for kidney disease. However, whether inflammatory status is either a cause or an outcome of chronic kidney disease remains controversial. We aimed to investigate the causal relationship between high-sensitivity C-reactive protein (hs-CRP) and estimated glomerular filtration rate (eGFR) using Mendelian randomization (MR) approaches. METHODS: A total of 10,521 participants of the Japan Multi-institutional Collaborative Cohort Study was analyzed in this study. We used two-sample MR approaches (the inverse-variance weighted (IVW), the weighted median (WM), and the MR-Egger method) to estimate the effect of genetically determined hs-CRP on kidney function. We selected four and three hs-CRP associated single nucleotide polymorphisms (SNPs) as two instrumental variables (IV): IVCRP and IVAsian, based on SNPs previously identified in European and Asian populations. IVCRP and IVAsian explained 3.4% and 3.9% of the variation in hs-CRP, respectively. RESULTS: Using the IVCRP, genetically determined hs-CRP was not significantly associated with eGFR in the IVW and the WM methods (estimate per 1 unit increase in ln(hs-CRP), 0.000; 95% confidence interval [CI], -0.019 to 0.020 and -0.003; 95% CI, -0.019 to 0.014, respectively). For IVAsian, we found similar results using the IVW and the WM methods (estimate, 0.005; 95% CI, -0.020 to 0.010 and -0.004; 95% CI, -0.020 to 0.012, respectively). The MR-Egger method also showed no causal relationships between hs-CRP and eGFR (IVCRP: -0.008; 95% CI, -0.058 to 0.042; IVAsian: 0.001; 95% CI, -0.036 to 0.036). CONCLUSION: Our two-sample MR analyses with different IVs did not support a causal effect of hs-CRP on eGFR.

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 (PASMC = 2.2 × 10-36), followed by relatively strong selection at the FAM96A (15q22), MYOF (10q23), 13q21, GRIA2 (4q32), and ASAP2 (2p25) loci (PASMC < 1.0 × 10-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.

Genome-wide association study of cervical cancer suggests a role for ARRDC3 gene in human papillomavirus infection.

The development of cervical cancer is initiated by human papillomavirus (HPV) infection and involves both viral and host genetic factors. Genome-wide association studies (GWAS) of cervical cancer have identified associations in the HLA locus and two loci outside HLA, but the principal genes that control infection and pathogenesis have not been identified. In the present study, we performed GWAS of cervical cancer in East Asian populations, involving 2609 cases and 4712 controls in the discovery stage and 1461 cases and 3295 controls in the follow-up stage. We identified novel-significant associations at 5q14 with the lead single nucleotide polymorphism (SNP) rs59661306 (P = 2.4 × 10-11) and at 7p11 with the lead SNP rs7457728 (P = 1.2 × 10-8). In 5q14, the chromatin region of the GWAS-significant SNPs was found to be in contact with the promoter of the ARRDC3 (arrestin domain-containing 3) gene. In our functional studies, ARRDC3 knockdown in HeLa cells caused significant reductions in both cell growth and susceptibility to HPV16 pseudovirion infection, suggesting that ARRDC3 is involved in the infectious entry of HPV into the cell. Our study advances the understanding of host genes that are responsible for cervical cancer susceptibility and guides future research on HPV infection and cancer development.