Molecular subtypes of lung adenocarcinoma present distinct immune tumor microenvironments.

Overcoming resistance to immune checkpoint inhibitors is an important issue in patients with non-small-cell lung cancer (NSCLC). Transcriptome analysis shows that adenocarcinoma can be divided into three molecular subtypes: terminal respiratory unit (TRU), proximal proliferative (PP), and proximal inflammatory (PI), and squamous cell carcinoma (LUSQ) into four. However, the immunological characteristics of these subtypes are not fully understood. In this study, we investigated the immune landscape of NSCLC tissues in molecular subtypes using a multi-omics dataset, including tumor-infiltrating leukocytes (TILs) analyzed using flow cytometry, RNA sequences, whole exome sequences, metabolomic analysis, and clinicopathologic findings. In the PI subtype, the number of TILs increased and the immune response in the tumor microenvironment (TME) was activated, as indicated by high levels of tertiary lymphoid structures, and high cytotoxic marker levels. Patient prognosis was worse in the PP subtype than in other adenocarcinoma subtypes. Glucose transporter 1 (GLUT1) expression levels were upregulated and lactate accumulated in the TME of the PP subtype. This could lead to the formation of an immunosuppressive TME, including the inactivation of antigen-presenting cells. The TRU subtype had low biological malignancy and "cold" tumor-immune phenotypes. Squamous cell carcinoma (LUSQ) did not show distinct immunological characteristics in its respective subtypes. Elucidation of the immune characteristics of molecular subtypes could lead to the development of personalized immune therapy for lung cancer. Immune checkpoint inhibitors could be an effective treatment for the PI subtype. Glycolysis is a potential target for converting an immunosuppressive TME into an antitumorigenic TME in the PP subtype.

Genetic insights into ossification of the posterior longitudinal ligament of the spine.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.

Tumor-infiltrating Leukocyte Profiling Defines Three Immune Subtypes of NSCLC with Distinct Signaling Pathways and Genetic Alterations.

Resistance to immune checkpoint blockade remains challenging in patients with non-small cell lung cancer (NSCLC). Tumor-infiltrating leukocyte (TIL) quantity, composition, and activation status profoundly influence responsiveness to cancer immunotherapy. This study examined the immune landscape in the NSCLC tumor microenvironment by analyzing TIL profiles of 281 fresh resected NSCLC tissues. Unsupervised clustering based on numbers and percentages of 30 TIL types classified adenocarcinoma (LUAD) and squamous cell carcinoma (LUSQ) into the cold, myeloid cell-dominant, and CD8+ T cell-dominant subtypes. These were significantly correlated with patient prognosis; the myeloid cell subtype had worse outcomes than the others. Integrated genomic and transcriptomic analyses, including RNA sequencing, whole-exome sequencing, T-cell receptor repertoire, and metabolomics of tumor tissue, revealed that immune reaction-related signaling pathways were inactivated, while the glycolysis and K-ras signaling pathways activated in LUAD and LUSQ myeloid cell subtypes. Cases with ALK and ROS1 fusion genes were enriched in the LUAD myeloid subtype, and the frequency of TERT copy-number variations was higher in LUSQ myeloid subtype than in the others. These classifications of NSCLC based on TIL status may be useful for developing personalized immune therapies for NSCLC. Significance: The precise TIL profiling classified NSCLC into novel three immune subtypes that correlates with patient outcome, identifying subtype-specific molecular pathways and genomic alterations that should play important roles in constructing subtype-specific immune tumor microenvironments. These classifications of NSCLC based on TIL status are useful for developing personalized immune therapies for NSCLC.

Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases.

The overwhelming majority of participants in current genetic studies are of European ancestry. To elucidate disease biology in the East Asian population, we conducted a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 320 independent signals in 276 loci for 27 diseases, with 25 novel loci (P < 9.58 × 10-9). East Asian-specific missense variants were identified as candidate causal variants for three novel loci, and we successfully replicated two of them by analyzing independent Japanese cohorts; p.R220W of ATG16L2 (associated with coronary artery disease) and p.V326A of POT1 (associated with lung cancer). We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, and identified 378 significant enrichments across nine diseases (false discovery rate < 0.05) (for example, NKX3-1 for prostate cancer). This large-scale GWAS in a Japanese population provides insights into the etiology of complex diseases and highlights the importance of performing GWAS in non-European populations.

IMSindel: An accurate intermediate-size indel detection tool incorporating de novo assembly and gapped global-local alignment with split read analysis.

Insertions and deletions (indels) have been implicated in dozens of human diseases through the radical alteration of gene function by short frameshift indels as well as long indels. However, the accurate detection of these indels from next-generation sequencing data is still challenging. This is particularly true for intermediate-size indels (≥50 bp), due to the short DNA sequencing reads. Here, we developed a new method that predicts intermediate-size indels using BWA soft-clipped fragments (unmatched fragments in partially mapped reads) and unmapped reads. We report the performance comparison of our method, GATK, PINDEL and ScanIndel, using whole exome sequencing data from the same samples. False positive and false negative counts were determined through Sanger sequencing of all predicted indels across these four methods. The harmonic mean of the recall and precision, F-measure, was used to measure the performance of each method. Our method achieved the highest F-measure of 0.84 in one sample, compared to 0.56 for GATK, 0.52 for PINDEL and 0.46 for ScanIndel. Similar results were obtained in additional samples, demonstrating that our method was superior to the other methods for detecting intermediate-size indels. We believe that this methodology will contribute to the discovery of intermediate-size indels associated with human disease.

A functional SNP in FLT1 increases risk of coronary artery disease in a Japanese population.

Coronary artery disease (CAD) including myocardial infarction is one of the leading causes of death in many countries. Similar to other common diseases, its pathogenesis is thought to result from complex interactions among multiple genetic and environmental factors. Recent large-scale genetic association analysis for CAD identified 15 new loci. We examined the reproducibility of these previous association findings with 7990 cases and 6582 controls in a Japanese population. We found a convincing association of rs9319428 in FLT1, encoding fms-related tyrosine kinase 1 (P=5.98 × 10(-8)). Fine mapping using tag single-nucleotide polymorphisms (SNPs) at FLT1 locus revealed that another SNP (rs74412485) showed more profound genetic effect for CAD (P=2.85 × 10(-12)). The SNP, located in intron 1 in FLT1, enhanced the transcriptional level of FLT1. RNA interference experiment against FLT1 showed that the suppression of FLT1 resulted in decreased expression of inflammatory adhesion molecules. Expression of FLT1 was observed in endothelial cells of human coronary artery. Our results indicate that the genetically coded increased expression of FLT1 by a functional SNP implicates activation in an inflammatory cascade that might eventually lead to CAD.

Molecular genetics of coronary artery disease.

Coronary artery disease (CAD) including myocardial infarction (MI) is a common disease and among the leading cause of death in the world. The onset of CAD depends on complex interactions of environmental and genetic factors. To clarify the genetic architecture of MI, we started a genome-wide association study (GWAS) using nearly 100 000 gene-based single-nucleotide polymorphisms (SNPs) from 2000, and identified LTA associated with the increased risk of MI in Japanese population. To our knowledge, this is the first study identified a genetic factor for common disease by GWAS in the worldwide. Through examining the LTA cascade by combination of molecular biological and genetic analyses, we have identified additional MI susceptible genes, LGALS2, PSMA6 and BRAP, so far. Nowadays a lot of large-scale GWAS have identified numerous genetic risk factors for common diseases. In CAD, 51 loci with GWAS significance (P<5 × 10(-8)) have collectively identified by recent large-scale GWAS mainly in Caucasian descent. In this review, we discuss recent advances in molecular genetics for CAD.

Decreased mortality associated with statin treatment in patients with acute myocardial infarction and lymphotoxin-alpha C804A polymorphism.

AIMS: We previously reported the association of single nucleotide polymorphisms in the lymphotoxin alpha (LTα) gene with susceptibility to acute myocardial infarction (AMI) and increased mortality after discharge. In the present study, we investigated whether the adverse effect of LTα C804A polymorphism on mortality could be pharmacologically modified by statin treatment after AMI. METHODS AND RESULTS: We conducted a multicenter study that included 3486 post-AMI patients between 1998 and 2008. During a median follow-up period of 1775 days, 247 deaths were recorded. The mortality rate was significantly higher in LTα 804A allele carriers compared to non-804A allele carriers (7.9% vs. 5.7%, p = 0.011). The LTα 804A allele was significantly associated with increased mortality for post-AMI patients not receiving statins (hazard ratio [HR]: 1.48, 95% confidence interval [CI]: 1.03-2.12, p = 0.034), but not for those receiving statins (HR: 1.22, 95% CI: 0.70-2.10, p = 0.486). In-vitro experimental analyses demonstrated that the LTα 804A polymorphic protein, 26Asn-LTα3, induced monocyte-endothelial interaction and endoplasmic reticulum (ER) stress in cardiomyocytes more strongly than the LTα3 804C polymorphic protein 26Thr-LTα3. However, the effects of both LTα3 proteins were decreased and became comparable by the pretreatment of cells with pravastatin. CONCLUSION: LTα C804A polymorphism was associated with an increased risk of mortality for AMI patients, although this effect was masked in patients treated with statins. This finding is supported by the observed attenuation of 26Asn-LTα3-mediated monocyte-endothelial interaction and ER stress in cardiomyocytes treated with pravastatin. LTα C804A polymorphism may have potential as a novel therapeutic target for secondary prevention after AMI.

Haplotypes with copy number and single nucleotide polymorphisms in CYP2A6 locus are associated with smoking quantity in a Japanese population.

Smoking is a major public health problem, but the genetic factors associated with smoking behaviors are not fully elucidated. Here, we have conducted an integrated genome-wide association study to identify common copy number polymorphisms (CNPs) and single nucleotide polymorphisms (SNPs) associated with the number of cigarettes smoked per day (CPD) in Japanese smokers (N = 17,158). Our analysis identified a common CNP with a strong effect on CPD (rs8102683; P=3.8 x 10(-42)) in the 19q13 region, encompassing the CYP2A6 locus. After adjustment for the associated CNP, we found an additional associated SNP (rs11878604; P=9.7 x 10(-30)) located 30 kb downstream of the CYP2A6 gene. Imputation of the CYP2A6 locus revealed that haplotypes underlying the CNP and the SNP corresponded to classical, functional alleles of CYP2A6 gene that regulate nicotine metabolism and explained 2% of the phenotypic variance of CPD (ANOVA F-test P=9.5 x 10(-52)). These haplotypes were also associated with smoking-related diseases, including lung cancer, chronic obstructive pulmonary disease and arteriosclerosis obliterans.

BRAP Activates Inflammatory Cascades and Increases the Risk for Carotid Atherosclerosis.

The BRCA-1 associated protein gene (BRAP) was recently identified as a susceptibility gene for myocardial infarction (MI). In the present study we aimed to decipher the association between the BRAP polymorphism and carotid atherosclerosis and the mechanism underlying its proatherogenic effect. A total of 1749 stroke/MI-free volunteers received carotid ultrasonic examinations for the measurement of intima-medial thickness (IMT) and plaque. The promoter polymorphism rs11066001 was selected because it affects the transcription of BRAP. We found that the GG genotype was associated with a 1.58-fold increased risk for having at least one plaque compared to carrying the A allele (P = 0.021). When subjects were divided by the cutoff value of IMT above the mean plus 1 standard deviation, there was an overrepresentation of the GG genotype in the subjects with thicker IMT (P = 0.004). The expression of BRAP increased significantly when human aortic smooth muscle cells (HASMCs) were treated with lipopolysaccharide (LPS). HASMCs were transfected with small interfering RNA against BRAP or scrambled sequences before treatment with LPS. Knockdown of BRAP led to attenuated HASMC proliferation and reduced secretion of monocyte chemoattractant protein-1 (MCP-1) and interleukin-8 (IL-8) in response to LPS. Downregulation of BRAP did not affect the protein levels of nuclear factor-κB (NF-κB), but prohibited its nuclear translocation. Coimmunoprecipitation experiments confirmed an interaction between BRAP and the two major components of the IKK signalosome, IκBß and IKKß. Collectively, BRAP conferred a risk for carotid plaque and IMT. Inflammatory stimuli upregulated BRAP expression, and BRAP activated inflammatory cascades by regulating NF-κB nuclear translocation.

SNPs in BRAP associated with risk of myocardial infarction in Asian populations.

Myocardial infarction is a common disease and among the leading causes of death in the world. We previously reported association of variants in LGALS2, encoding galectin-2, with myocardial infarction susceptibility in a case-control association study in a Japanese population. Here we identify BRAP (BRCA1-associated protein) as a galectin-2-binding protein. We report an association of SNPs in BRAP with myocardial infarction risk in a large Japanese cohort (P = 3.0 x 10(-18), OR = 1.48, 2,475 cases and 2,778 controls), with replication in additional Japanese and Taiwanese cohorts (P = 4.4 x 10(-6), 862 cases and 1,113 controls and P = 4.7 x 10(-3), 349 cases and 994 controls, respectively). BRAP expression was observed in smooth muscle cells (SMCs) and macrophages in human atherosclerotic lesions. BRAP knockdown by siRNA using cultured coronary endothelial cells suppressed activation of NF-kappaB, a central mediator of inflammation.

Lymphotoxin-alpha3 mediates monocyte-endothelial interaction by TNFR I/NF-kappaB signaling.

We recently reported that the single nucleotide polymorphisms of the lymphotoxin-(LT)alpha gene, a member of the tumor necrosis factor (TNF) family, are closely related to acute myocardial infarction; however, the precise mechanism of LTalpha signaling in atherogenesis remains unclear. We investigated the role of LTalpha3, a secreted homotrimer of LTalpha, in monocyte-endothelial cell adhesion using cultured human umbilical vein endothelial cells (HUVEC). We found that LTalpha3 induced cell adhesion molecules and activated NF-kappaB p50 and p65. LTalpha3 also induced phosphorylation of Akt, phosphorylation and degradation of IkappaB, nuclear translocation of p65, and increased adhesion of THP1 monocytes to HUVEC. These effects were mediated by TNF receptor (TNFR) I and attenuated by the phosphatidylinositol triphosphate-kinase (PI3K) inhibitors LY294002 and Wortmannin. Thus, LTalpha3 mediates the monocyte-endothelial interaction via the classical NF-kappaB pathway following TNFR I/PI3K activation, indicating it may play a role in the development of coronary artery disease.

A functional polymorphism in THBS2 that affects alternative splicing and MMP binding is associated with lumbar-disc herniation.

Lumbar-disc herniation (LDH), one of the most common musculoskeletal diseases, has strong genetic determinants. Recently, several genes that encode extracellular matrix (ECM) proteins in the intervertebral disc have been reported to associate with LDH. Thrombospondins (THBSs) 1 and 2 are good candidates for the LDH susceptibility gene: They are intervertebral disc ECM proteins that regulate the effective levels of matrix metalloproteinases (MMPs) 2 and 9, which are key effectors of ECM remodeling. Here, we report that THBS2 is associated with LDH in Japanese populations. An intronic SNP in THBS2 (IVS10-8C --> T; rs9406328) showed significant association (p = 0.0000028) with LDH in two independent Japanese populations. This SNP, located in a polypyrimidine tract upstream of the 3' splice site of intron 10, exerts allelic differences on exon 11 skipping rates in vivo, with the susceptibility allele showing increased skipping. Skipping of exon 11 results in decreased THBS2 interaction with MMP2 and MMP9. Further, a missense SNP in MMP9 (Q279R; rs17576) is also strongly associated with LDH in the Japanese population (p = 0.00049) and shows a combinatorial effect with THBS2 (odds ratio 3.03, 95% confidence interval 1.58-5.77). Thus, a splicing-affecting SNP in THBS2 and a missense SNP in MMP9 are associated with susceptibility to LDH. Our data indicate that regulation of intervertebral disc ECM metabolism by the THBS2-MMP system plays an essential role in the etiology and pathogenesis of LDH.

Genome-wide association study to identify single-nucleotide polymorphisms conferring risk of myocardial infarction.

Myocardial infarction (MI) is characterized by abrupt occlusion of coronary artery resulting in irreversible damage to cardiac muscle. This disease might result from the interactions of multiple genetic and environmental factors, none of which can cause disease solely by themselves. To reveal the genetic bases of MI, we performed a large-scale, case-control association study using 92,788 gene-based single-nucleotide polymorphism (SNP) markers. We have identified functional SNPs within the lymphotoxin-alpha gene located on chromosome 6p21 conferred susceptibility to MI. This chapter describes a detailed protocol for performing a genome-wide association study as used in our MI study.

Inflammation as a risk factor for myocardial infarction.

Myocardial infarction (MI) is one of the common diseases whose pathogenesis includes genetic factors. To reveal genetic backgrounds of this clinically heterogeneous disorder, we started our case-control association study by examining large-scale gene-based single nucleotide polymorphism (SNP) sets for approximately 1,000 patients and controls. As a core genotyping method, a combination of multiplex PCR and Invader method was used, and after genotyping approximately 65,000 SNPs, we found two SNPs located within lymphotoxin-alpha (LTA) gene showing significant association with MI. LTA is one of the cytokines produced in the early stages of vascular inflammatory processes. These SNPs seem to be involved in inflammation by both qualitatively and quantitatively modifying the function of LTA protein, thereby conferring a risk of MI. The genetic association was further confirmed by other researchers using white European trios. To further understand the roles of LTA protein in the pathogenesis of MI, we searched for proteins that interact directly with LTA protein and identified galectin-2 protein as a binding partner of LTA protein. It is a member of galactose-binding lectin family whose function has not been well characterized. Genetic association study again revealed that an SNP in LGALS2 encoding galectin-2 was also associated with susceptibility to MI. This genetic substitution seemed to affect the transcriptional level of galectin-2, which led to altered secretion of LTA, thereby affecting the degree of inflammation. Thus, our findings indicate the importance of inflammation, especially the LTA cascade, in the pathogenesis of MI. Also, combined strategy of genetic and molecular-cellular biological approaches may be useful for clarification of the pathogenesis of common diseases in general.

Impact of atherosclerosis-related gene polymorphisms on mortality and recurrent events after myocardial infarction.

Although previous epidemiologic studies have suggested an association between the onset of myocardial infarction (MI) and some genetic variations, the impact of these variants on recurrent cardiovascular events after MI has not been fully elucidated. We genotyped 87 polymorphisms of 73 atherosclerosis-related genes in consecutive acute MI patients registered in the Osaka Acute Coronary Insufficiency Study and compared the incidence of death and major adverse cardiac events (MACE) among the polymorphisms of each gene. After initial screening in 507 patients, we selected nine polymorphisms for screening in all 1586 patients. Multivariate Cox regression analysis revealed that G allele carriers at the position 252 of the lymphotoxin alpha (LTA) gene were independently associated with an increased risk of death (hazard ratio [HR]: 2.46; 95% CI: 1.24-4.86). In conclusion, a 252G allele of LTA is associated with an increased risk of death after AMI and may be a useful genetic predictor.

Functional variation in LGALS2 confers risk of myocardial infarction and regulates lymphotoxin-alpha secretion in vitro.

Myocardial infarction (MI) has become one of the leading causes of death in the world. Its pathogenesis includes chronic formation of plaque inside the vessel wall of the coronary artery and acute rupture of the artery, implicating a number of inflammation-mediating molecules, such as the cytokine lymphotoxin-alpha (LTA). Functional variations in LTA are associated with susceptibility to MI. Here we show that LTA protein binds to galectin-2, a member of the galactose-binding lectin family. Our case-control association study in a Japanese population showed that a single nucleotide polymorphism in LGALS2 encoding galectin-2 is significantly associated with susceptibility to MI. This genetic substitution affects the transcriptional level of galectin-2 in vitro, potentially leading to altered secretion of LTA, which would then affect the degree of inflammation; however, its relevance to other populations remains to be clarified. Smooth muscle cells and macrophages in the human atherosclerotic lesions expressed both galectin-2 and LTA. Our findings thus suggest a link between the LTA cascade and the pathogenesis of MI.

Functional SNPs in the lymphotoxin-alpha gene that are associated with susceptibility to myocardial infarction.

By means of a large-scale, case-control association study using 92,788 gene-based single-nucleotide polymorphism (SNP) markers, we identified a candidate locus on chromosome 6p21 associated with susceptibility to myocardial infarction. Subsequent linkage-disequilibrium (LD) mapping and analyses of haplotype structure showed significant associations between myocardial infarction and a single 50 kb halpotype comprised of five SNPs in LTA (encoding lymphotoxin-alpha), NFKBIL1 (encoding nuclear factor of kappa light polypeptide gene enhancer in B cells, inhibitor-like 1) and BAT1 (encoding HLA-B associated transcript 1). Homozygosity with respect to each of the two SNPs in LTA was significantly associated with increased risk for myocardial infarction (odds ratio = 1.78, chi(2) = 21.6, P = 0.00000033; 1,133 affected individuals versus 1,006 controls). In vitro functional analyses indicated that one SNP in the coding region of LTA, which changed an amino-acid residue from threonine to asparagine (Thr26Asn), effected a twofold increase in induction of several cell-adhesion molecules, including VCAM1, in vascular smooth-muscle cells of human coronary artery. Moreover, the SNP, in intron 1 of LTA, enhanced the transcriptional level of LTA. These results indicate that variants in the LTA are risk factors for myocardial infraction and implicate LTA in the pathogenesis of the disorder.