Immuno-genomic analysis reveals eosinophilic feature and favorable prognosis of female non-smoking esophageal squamous cell carcinomas.

Most of esophageal squamous cell carcinoma (ESCC) develop in smoking males in Japan, but the genomic etiology and immunological characteristics of rare non-smoking female ECSS remain unclear. To elucidate the genomic and immunological features of ESCC in non-smoking females, we analyzed whole-genome or transcriptome sequencing data from 94 ESCCs, including 20 rare non-smoking female cases. In addition, 31,611 immune cells were extracted from four ESCC tissues and subject to single-cell RNA-seq. We compared their immuno-genomic and microbiome profiles between non-smoking female and smoking ESCCs. Non-smoking females showed much better prognosis. Whole-genome sequencing analysis showed no significant differences in driver genes or copy number alterations depending on smoking status. The mutational signatures specifically observed in non-smoking females ESCC could be attributed to aging. Immune profiling from RNA-seq revealed that ESCC in non-smoking females had high tumor microenvironment signatures and a high abundance of eosinophils with a favorable prognosis. Single-cell RNA-sequencing of intratumor immune cells revealed gender differences of eosinophils and their activation in female cases. ESCCs in non-smoking females have age-related mutational signatures and gender-specific tumor immune environment with eosinophils, which is likely to contribute to their favorable prognosis.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma.

BACKGROUND: Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated. METHODS: To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours. FINDINGS: RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients. INTERPRETATION: These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology. FUNDING: This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Genotype-phenotype mapping of a patient-derived lung cancer organoid biobank identifies NKX2-1-defined Wnt dependency in lung adenocarcinoma.

Human lung cancer is a constellation of tumors with various histological and molecular properties. To build a preclinical platform that covers this broad disease spectrum, we obtained lung cancer specimens from multiple sources, including sputum and circulating tumor cells, and generated a living biobank consisting of 43 lines of patient-derived lung cancer organoids. The organoids recapitulated the histological and molecular hallmarks of the original tumors. Phenotypic screening of niche factor dependency revealed that EGFR mutations in lung adenocarcinoma are associated with the independence from Wnt ligands. Gene engineering of alveolar organoids reveals that constitutive activation of EGFR-RAS signaling provides Wnt independence. Loss of the alveolar identity gene NKX2-1 confers Wnt dependency, regardless of EGFR signal mutation. Sensitivity to Wnt-targeting therapy can be stratified by the expression status of NKX2-1. Our results highlight the potential of phenotype-driven organoid screening and engineering for the fabrication of therapeutic strategies to combat cancer.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Hereditary cancer variants and homologous recombination deficiency in biliary tract cancer.

BACKGROUND & AIMS: The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs. METHODS: We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status. RESULTS: Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p <0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants without a second hit or variants of other homologous recombination-related genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes. CONCLUSIONS: In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs. IMPACT AND IMPLICATIONS: We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Immuno-genomic profiling of biopsy specimens predicts neoadjuvant chemotherapy response in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

Association of an IGHV3-66 gene variant with Kawasaki disease.

In a meta-analysis of three GWAS for susceptibility to Kawasaki disease (KD) conducted in Japan, Korea, and Taiwan and follow-up studies with a total of 11,265 subjects (3428 cases and 7837 controls), a significantly associated SNV in the immunoglobulin heavy variable gene (IGHV) cluster in 14q33.32 was identified (rs4774175; OR = 1.20, P = 6.0 × 10-9). Investigation of nonsynonymous SNVs of the IGHV cluster in 9335 Japanese subjects identified the C allele of rs6423677, located in IGHV3-66, as the most significant reproducible association (OR = 1.25, P = 6.8 × 10-10 in 3603 cases and 5731 controls). We observed highly skewed allelic usage of IGHV3-66, wherein the rs6423677 A allele was nearly abolished in the transcripts in peripheral blood mononuclear cells of both KD patients and healthy adults. Association of the high-expression allele with KD strongly indicates some active roles of B-cells or endogenous immunoglobulins in the disease pathogenesis. Considering that significant association of SNVs in the IGHV region with disease susceptibility was previously known only for rheumatic heart disease (RHD), a complication of acute rheumatic fever (ARF), these observations suggest that common B-cell related mechanisms may mediate the symptomology of KD and ARF as well as RHD.

Improved Parkinsons disease motor score in a single-arm open-label trial of febuxostat and inosine.

BACKGROUND: Cellular energetics play an important role in Parkinsons disease etiology, but no treatments directly address this deficiency. Our past research showed that treatment with febuxostat and inosine increased blood hypoxanthine and ATP in healthy adults, and a preliminary trial in 3 Parkinson's disease patients suggested some symptomatic improvements with no adverse effects. METHODS: To examine the efficacy on symptoms and safety in a larger group of Parkinsons disease patients, we conducted a single-arm, open-label trial at 5 Japanese neurology clinics and enrolled thirty patients (nmales = 11; nfemales = 19); 26 patients completed the study (nmales = 10; nfemales = 16). Each patient was administered febuxostat 20 mg and inosine 500 mg twice-per-day (after breakfast and dinner) for 8 weeks. The primary endpoint was the difference of MDS-UPDRS Part III score immediately before and after 57 days of treatment. RESULTS: Serum hypoxanthine concentrations were raised significantly after treatment (Pre = 11.4 µM; Post = 38.1 µM; P < .0001). MDS-UPDRS Part III score was significantly lower after treatment (Pre = 28.1 ±â€Š9.3; Post = 24.7 ±â€Š10.8; mean ±â€ŠSD; P = .0146). Sixteen adverse events occurred in 13/29 (44.8%) patients, including 1 serious adverse event (fracture of the second lumbar vertebra) that was considered not related to the treatment. CONCLUSIONS: The results of this study suggest that co-administration of febuxostat and inosine is relatively safe and effective for improving symptoms of Parkinsons disease patients. Further controlled trials need to be performed to confirm the symptomatic improvement and to examine the disease-modifying effect in long-term trials.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Shortage of Cellular ATP as a Cause of Diseases and Strategies to Enhance ATP.

Germline mutations in cellular-energy associated genes have been shown to lead to various monogenic disorders. Notably, mitochondrial disorders often impact skeletal muscle, brain, liver, heart, and kidneys, which are the body's top energy-consuming organs. However, energy-related dysfunctions have not been widely seen as causes of common diseases, although evidence points to such a link for certain disorders. During acute energy consumption, like extreme exercise, cells increase the favorability of the adenylate kinase reaction 2-ADP -> ATP+AMP by AMP deaminase degrading AMP to IMP, which further degrades to inosine and then to purines hypoxanthine -> xanthine -> urate. Thus, increased blood urate levels may act as a barometer of extreme energy consumption. AMP deaminase deficient subjects experience some negative effects like decreased muscle power output, but also positive effects such as decreased diabetes and improved prognosis for chronic heart failure patients. That may reflect decreased energy consumption from maintaining the pool of IMP for salvage to AMP and then ATP, since de novo IMP synthesis requires burning seven ATPs. Similarly, beneficial effects have been seen in heart, skeletal muscle, or brain after treatment with allopurinol or febuxostat to inhibit xanthine oxidoreductase, which catalyzes hypoxanthine -> xanthine and xanthine -> urate reactions. Some disorders of those organs may reflect dysfunction in energy-consumption/production, and the observed beneficial effects related to reinforcement of ATP re-synthesis due to increased hypoxanthine levels in the blood and tissues. Recent clinical studies indicated that treatment with xanthine oxidoreductase inhibitors plus inosine had the strongest impact for increasing the pool of salvageable purines and leading to increased ATP levels in humans, thereby suggesting that this combination is more beneficial than a xanthine oxidoreductase inhibitor alone to treat disorders with ATP deficiency.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Genome-wide association study in Japanese females identifies fifteen novel skin-related trait associations.

Skin trait variation impacts quality-of-life, especially for females from the viewpoint of beauty. To investigate genetic variation related to these traits, we conducted a GWAS of various skin phenotypes in 11,311 Japanese women and identified associations for age-spots, freckles, double eyelids, straight/curly hair, eyebrow thickness, hairiness, and sweating. In silico annotation with RoadMap Epigenomics epigenetic state maps and colocalization analysis of GWAS and GTEx Project eQTL signals provided information about tissue specificity, candidate causal variants, and functional target genes. Novel signals for skin-spot traits neighboured AKAP1/MSI2 (rs17833789; P = 2.2 × 10-9), BNC2 (rs10810635; P = 2.1 × 10-22), HSPA12A (rs12259842; P = 7.1 × 10-11), PPARGC1B (rs251468; P = 1.3 × 10-21), and RAB11FIP2 (rs10444039; P = 5.6 × 10-21). HSPA12A SNPs were the only protein-coding gene eQTLs identified across skin-spot loci. Double edged eyelid analysis identified that a signal around EMX2 (rs12570134; P = 8.2 × 10-15) was also associated with expression of EMX2 and the antisense-RNA gene EMX2OS in brain putamen basal ganglia tissue. A known hair morphology signal in EDAR was associated with both eyebrow thickness (rs3827760; P = 1.7 × 10-9) and straight/curly hair (rs260643; P = 1.6 × 10-103). Excessive hairiness signals' top SNPs were also eQTLs for TBX15 (rs984225; P = 1.6 × 10-8), BCL2 (rs7226979; P = 7.3 × 10-11), and GCC2 and LIMS1 (rs6542772; P = 2.2 × 10-9). For excessive sweating, top variants in two signals in chr2:28.82-29.05 Mb (rs56089836; P = 1.7 × 10-11) were eQTLs for either PPP1CB or PLB1, while a top chr16:48.26-48.45 Mb locus SNP was a known ABCC11 missense variant (rs6500380; P = 6.8 × 10-10). In total, we identified twelve loci containing sixteen association signals, of which fifteen were novel. These findings will help dermatologic researchers better understand the genetic underpinnings of skin-related phenotypic variation in human populations.

Japanese GWAS identifies variants for bust-size, dysmenorrhea, and menstrual fever that are eQTLs for relevant protein-coding or long non-coding RNAs.

Traits related to primary and secondary sexual characteristics greatly impact females during puberty and day-to-day adult life. Therefore, we performed a GWAS analysis of 11,348 Japanese female volunteers and 22 gynecology-related phenotypic variables, and identified significant associations for bust-size, menstrual pain (dysmenorrhea) severity, and menstrual fever. Bust-size analysis identified significant association signals in CCDC170-ESR1 (rs6557160; P = 1.7 × 10-16) and KCNU1-ZNF703 (rs146992477; P = 6.2 × 10-9) and found that one-third of known European-ancestry associations were also present in Japanese. eQTL data points to CCDC170 and ZNF703 as those signals' functional targets. For menstrual fever, we identified a novel association in OPRM1 (rs17181171; P = 2.0 × 10-8), for which top variants were eQTLs in multiple tissues. A known dysmenorrhea signal near NGF replicated in our data (rs12030576; P = 1.1 × 10-19) and was associated with RP4-663N10.1 expression, a putative lncRNA enhancer of NGF, while a novel dysmenorrhea signal in the IL1 locus (rs80111889; P = 1.9 × 10-16) contained SNPs previously associated with endometriosis, and GWAS SNPs were most significantly associated with IL1A expression. By combining regional imputation with colocalization analysis of GWAS/eQTL signals along with integrated annotation with epigenomic data, this study further refines the sets of candidate causal variants and target genes for these known and novel gynecology-related trait loci.

Genome-wide association study of self-reported food reactions in Japanese identifies shrimp and peach specific loci in the HLA-DR/DQ gene region.

Food allergy is an increasingly important health problem in the world. Several genome-wide association studies (GWAS) focused on European ancestry samples have identified food allergy-specific loci in the HLA class II region. We conducted GWAS of self-reported reactivity with common foods using the data from 11011 Japanese women and identified shrimp and peach allergy-specific loci in the HLA-DR/DQ gene region tagged by rs74995702 (P = 6.30 × 10-17, OR = 1.91) and rs28359884 (P = 2.3 × 10-12, OR = 1.80), respectively. After HLA imputation using a Japanese population-specific reference, the most strongly associated haplotype was HLA-DRB1*04:05-HLA-DQB1*04:01 for shrimp allergy (P = 3.92 × 10-19, OR = 1.99) and HLA-DRB1*09:01-HLA-DQB1*03:03 for peach allergy (P = 1.15 × 10-7, OR = 1.68). Additionally, both allergies' associated variants were eQTLs for several HLA genes, with HLA-DQA2 the single eQTL gene shared between the two traits. Our study suggests that allergy to certain foods may be related to genetic differences that tag both HLA alleles having particular epitope binding specificities as well as variants modulating expression of particular HLA genes. Investigating this further could increase our understanding of food allergy aetiology and potentially lead to better therapeutic strategies for allergen immunotherapies.

A genome-wide association analysis identifies NMNAT2 and HCP5 as susceptibility loci for Kawasaki disease.

Kawasaki disease (KD), a systemic vasculitis of infants and children, manifests as fever and mucocutaneous inflammation. Although its etiology is largely unknown, the epidemiological data suggest that genetic factors are important in KD susceptibility. To identify genetic variants influencing KD susceptibility, we performed a genome-wide association study (GWAS) and replication study using a total of 915 children with KD and 4553 controls in the Korean population. Six single-nucleotide polymorphisms (SNPs) in three loci were associated significantly with KD susceptibility (P<1.0 × 10-5), including the previously reported BLK locus (rs6993775, odds ratio (OR)=1.52, P=2.52 × 10-11). The other two loci were newly identified: NMNAT2 on chromosome 1q25.3 (rs2078087, OR=1.33, P=1.15 × 10-6) and the human leukocyte antigen (HLA) region on chromosome 6p21.3 (HLA-C, HLA-B, MICA and HCP5) (rs9380242, rs9378199, rs9266669 and rs6938467; OR=1.33-1.51, P=8.93 × 10-6 to 5.24 × 10-8). Additionally, SNP rs17280682 in NLRP14 was associated significantly with KD with a family history (18 cases vs 4553 controls, OR=6.76, P=5.46 × 10-6). These results provide new insights into the pathogenesis and pathophysiology of KD.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.