Raman Multi-Omic Snapshot and Statistical Validation of Structural Differences between Herpes Simplex Type I and Epstein-Barr Viruses.

Raman spectroscopy was applied to study the structural differences between herpes simplex virus Type I (HSV-1) and Epstein-Barr virus (EBV). Raman spectra were first collected with statistical validity on clusters of the respective virions and analyzed according to principal component analysis (PCA). Then, average spectra were computed and a machine-learning approach applied to deconvolute them into sub-band components in order to perform comparative analyses. The Raman results revealed marked structural differences between the two viral strains, which could mainly be traced back to the massive presence of carbohydrates in the glycoproteins of EBV virions. Clear differences could also be recorded for selected tyrosine and tryptophan Raman bands sensitive to pH at the virion/environment interface. According to the observed spectral differences, Raman signatures of known biomolecules were interpreted to link structural differences with the viral functions of the two strains. The present study confirms the unique ability of Raman spectroscopy for answering structural questions at the molecular level in virology and, despite the structural complexity of viral structures, its capacity to readily and reliably differentiate between different virus types and strains.

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions.

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

Cold-induced vasodilation response in a Japanese cohort: insights from cold-water immersion and genome-wide association studies.

BACKGROUND: Cold-induced vasodilation (CIVD) occurs after blood vessels in the skin are constricted due to local cold exposure. Although many CIVD studies have been conducted, the underlying molecular mechanisms are yet to be clarified. Therefore, we explored genetic variants associated with CIVD response using the largest-scale dataset reported to date in a CIVD study involving wavelet analysis; thus, the findings improve our understanding of the molecular mechanisms that regulate the CIVD response. METHODS: We performed wavelet analysis of three skin blood flow signals [endothelial nitric oxide (eNO)-independent, eNO-dependent, and neurogenic activities] during finger cold-water immersion at 5 °C in 94 Japanese young adults. Additionally, we conducted genome-wide association studies of CIVD using saliva samples collected from the participants. RESULTS: We found that the mean wavelet amplitudes of eNO-independent and neurogenic activities significantly increased and decreased prior to CIVD, respectively. Our results also implied that as many as ~ 10% of the Japanese subjects did not show an apparent CIVD response. Our genome-wide association studies of CIVD using ~ 4,040,000 imputed data found no apparent CIVD-related genetic variants; however, we identified 10 genetic variants, including 2 functional genes (COL4A2 and PRLR) that are associated with notable blunted eNO-independent and neurogenic activity responses in individuals without CIVD response during local cold exposure. CONCLUSIONS: Our findings indicate that individuals without CIVD response differentiated by genotypes with COL4A2 and PRLR genetic variants exhibited notable blunted eNO-independent and neurogenic activity responses during local cold exposure.

OLIG2 is an in vivo bookmarking transcription factor in the developing neural tube in mouse.

Cells possess intrinsic features that are inheritable via epigenetic regulation, such as DNA methylation and histone modification. These inheritable features maintain a unique gene expression pattern, underlying cellular memory. Because of the degradation or displacement of mitotic chromosomes, most transcription factors do not contribute to cellular memory. However, accumulating in vitro evidence indicates that some transcription factors can be retained in mitotic chromosomes called as bookmarking. Such transcription factors may contribute to a novel third mechanism of cellular memory. Since most findings of transcription factor bookmarking have been reported in vitro, little is currently known in vivo. In the neural tube of mouse embryos, we discovered that OLIG2, a basic helix loop helix (bHLH) transcription factor that regulates proliferation of neural progenitors and the cell fate of motoneurons and oligodendrocytes, binds to chromatin through every cell cycle including M-phase. OLIG2 chromosomal localization coincides with mitotic cell features such as the phosphorylation of histone H3, KI67, and nuclear membrane breakdown. Chromosomal localization of OLIG2 is regulated by an N-terminus triple serine motif. Photobleaching analysis revealed slow OLIG2 mobility, suggesting a high affinity of OLIG2 to DNA. In Olig2 N-terminal deletion mutant mice, motoneurons and oligodendrocyte progenitor numbers are reduced in the neural tube, suggesting that the bookmarking regulatory domain is important for OLIG2 function. We conclude that OLIG2 is a de novo in vivo bookmarking transcription factor. Our results demonstrate the presence of in vivo bookmarking in a living organism and illustrate a novel function of transcription factors.

Deep learning-based predictions of clear and eosinophilic phenotypes in clear cell renal cell carcinoma.

We have recently shown that histological phenotypes focusing on clear and eosinophilic cytoplasm in clear cell renal cell carcinoma (ccRCC) correlated with prognosis and the response to angiogenesis inhibition and checkpoint blockade. This study aims to objectively show the diagnostic utility of clear or eosinophilic phenotypes of ccRCC by developing an artificial intelligence (AI) model using the TCGA-ccRCC dataset and to demonstrate if the clear or eosinophilic predicted phenotypes correlate with pathological factors and gene signatures associated with angiogenesis and cancer immunity. Before the development of the AI model, histological evaluation using hematoxylin and eosin whole-slide images of the TCGA-ccRCC cohort (n = 435) was performed by a urologic pathologist. The AI model was developed as follows. First, the highest-grade area on each whole slide image was captured for image processing. Second, the selected regions were cropped into tiles. Third, the AI model was trained using transfer learning on a deep convolutional neural network, and clear or eosinophilic predictions were scaled as AI scores. Next, we verified the AI model using a validation cohort (n = 95). Finally, we evaluated the accuracy of the prognostic predictions of the AI model and revealed that the AI model detected clear and eosinophilic phenotypes with high accuracy. The AI model stratified the patients' outcomes, and the predicted eosinophilic phenotypes correlated with adverse clinicopathological characteristics and high immune-related gene signatures. In conclusion, the AI-based histologic subclassification accurately predicted clear or eosinophilic phenotypes of ccRCC, allowing for consistently reproducible stratification for prognostic and therapeutic stratification.

Histologic-Based Tumor-Associated Immune Cells Status in Clear Cell Renal Cell Carcinoma Correlates with Gene Signatures Related to Cancer Immunity and Clinical Outcomes.

The three-tier immunophenotype (desert, excluded, and inflamed) and the four-tier immunophenotype (cold, immunosuppressed, excluded, and hot) have been linked to prognosis and immunotherapy response. This study aims to evaluate whether immunophenotypes of clear cell renal cell carcinoma, identified on hematoxylin and eosin-stained slides, correlate with gene expression signatures related to cancer immunity, and clinical outcomes. We evaluated tumor-associated immune cells (TAICs) status using three methodologies: three-tier immunophenotype based on the location of TAICs, four-tier immunophenotype considering both the location and degree of TAICs and inflammation score focusing only on the degree of TAICs, using a localized clear cell renal cell carcinoma cohort (n = 436) and The Cancer Genome Atlas (TCGA)-KIRC cohort (n = 162). We evaluated the association of the TAICs status assessed by three methodologies with CD8 and PD-L1 immunohistochemistry and immune gene expression signatures by TCGA RNA-sequencing data. All three methodologies correlated with immunohistochemical and immune gene expression signatures. The inflammation score and the four-tier immunophenotype showed similarly higher accuracy in predicting recurrence-free survival and overall survival compared to the three-tier immunophenotype. In conclusion, a simple histologic assessment of TIACs may predict clinical outcomes and immunotherapy responses.

PBRM1 Immunohistochemical Expression Profile Correlates with Histomorphological Features and Endothelial Expression of Tumor Vasculature for Clear Cell Renal Cell Carcinoma.

Loss of the polybromo-1 (PBRM1) protein has been expected as a possible biomarker for clear cell renal cell carcinoma (ccRCC). There is little knowledge about how PBRM1 immunohistochemical expression correlates with the histomorphological features of ccRCC and the endothelial expression of tumor vasculature. The present study evaluates the association of architectural patterns with the PBRM1 expression of cancer cells using a cohort of 425 patients with nonmetastatic ccRCC. Furthermore, we separately assessed the PBRM1 expression of the endothelial cells and evaluated the correlation between the expression of cancer cells and endothelial cells. PBRM1 loss in cancer cells was observed in 148 (34.8%) patients. In the correlation analysis between architectural patterns and PBRM1 expression, macrocyst/microcystic, tubular/acinar, and compact/small nested were positively correlated with PBRM1 expression, whereas alveolar/large nested, thick trabecular/insular, papillary/pseudopapillary, solid sheets, and sarcomatoid/rhabdoid were negatively correlated with PBRM1 expression. PBRM1 expression in vascular endothelial cells correlated with the expression of cancer cells (correlation coefficient = 0.834, p < 0.001). PBRM1 loss in both cancer and endothelial cells was associated with a lower recurrence-free survival rate (p < 0.001). Our PBRM1 expression profile indicated that PBRM1 expression in both cancer and endothelial cells may be regulated in an orchestrated manner.

Development and validation of a vascularity-based architectural classification for clear cell renal cell carcinoma: correlation with conventional pathological prognostic factors, gene expression patterns, and clinical outcomes.

The prognostic significance of an architectural grading system for clear cell renal cell carcinoma (ccRCC) has recently been demonstrated. The present study aimed to establish a vascularity-based architectural classification using the cohort of 436 patients with localized ccRCC who underwent extirpative surgery and correlated the findings with conventional pathologic factors, gene expression, and prognosis. First, we assessed architectural patterns in the highest-grade area on hematoxylin and eosin-stained slides, then separately evaluated our surrogate score for vascularity. We grouped nine architectural patterns into three categories based on the vascular network score. "Vascularity-based architectural classification" was defined: category 1: characterized by enrichment of the vascular network, including compact/small nested, macrocyst/microcystic, and tubular/acinar patterns; category 2: characterized by a widely spaced-out vascular network, including alveolar/large nested, thick trabecular/insular, papillary/pseudopapillary patterns; category 3: characterized by scattered vascularity without a vascular network, including solid sheets, rhabdoid and sarcomatoid patterns. Adverse pathological prognostic factors such as TNM stage, WHO/ISUP grade, and necrosis were significantly associated with category 3, followed by category 2 (all p < 0.001). We successfully validated the classification using The Cancer Genome Atlas (TCGA) cohort (n = 162), and RNA-sequencing data available from TCGA showed that the angiogenesis gene signature was significantly enriched in category 1 compared to categories 2 and 3, whereas the immune gene signature was significantly enriched in category 3 compared to categories 1 and 2. In univariate analysis, vascularity-based architectural classification showed the best accuracy in pathological prognostic factors for predicting recurrence-free survival (c-index = 0.786). The predictive accuracy of our model which integrated WHO/ISUP grade, necrosis, TNM stage, and vascularity-based architectural classification was greater than conventional risk models (c-index = 0.871 vs. 0.755-0.843). Our findings suggest that the vascularity-based architectural classification is prognostically useful and may help stratify patients appropriately for management based on their likelihood of post-surgical recurrence.

Upregulation of cathepsin L gene under mild cold conditions in young Japanese male adults.

BACKGROUND: Physiological thermoregulatory systems in humans have been a key factor for adaptation to local environments after their exodus from Africa, particularly, to cold environments outside Africa. Recent studies using high-throughput sequencing have identified various genes responsible for cold adaptation. However, the molecular mechanisms underlying initial thermoregulation in response to acute cold exposure remain unclear. Therefore, we investigated transcriptional profiles of six young Japanese male adults exposed to acute cold stress. METHODS: In a climatic chamber, the air temperature was maintained at 28°C for 65 min and was then gradually decreased to 19°C for 70 min. Saliva samples were obtained from the subjects at 28°C before and after 19°C cold exposure and were used for RNA sequencing. RESULTS: In the cold exposure experiment, expression levels of 14 genes were significantly changed [false discovery rate (FDR) < 0.05] although the degree of transcriptional changes was not high due to experimental conditions or blunted transcriptional reaction in saliva to cold stress. As a result, differential gene expression analyses detected the cathepsin L (CTSL) gene to be significantly upregulated, with FDR < 0.05 and log2 fold change value > 1; thus, this gene was identified as a differentially expressed gene. Given that the cathepsin L protein is related to invasion of the novel coronavirus (SARS-CoV-2), mild cold stress might alter the susceptibility to coronavirus disease-19 in humans. The gene ontology enrichment analysis for 14 genes with FDR < 0.05 suggested that immune-related molecules could be activated by mild cold stress. CONCLUSIONS: The results obtained from this study indicate that CTSL expression levels can be altered by acute mild cold stress.

Effect of EGLN1 Genetic Polymorphisms on Hemoglobin Concentration in Andean Highlanders.

The physiological characteristics of Andean natives living at high altitudes have been investigated extensively, with many studies reporting that Andean highlanders have a higher hemoglobin (Hb) concentration than other highlander populations. It has previously been reported that positive natural selection has acted independently on the egl-9 family hypoxia inducible factor 1 (EGLN1) gene in Tibetan and Andean highlanders and is related to Hb concentration in Tibetans. However, no study has yet revealed the genetic determinants of Hb concentration in Andeans even though several single-nucleotide polymorphisms (SNPs) in EGLN1 have previously been examined. Therefore, we explored the relationship between hematological measurements and tag SNPs designed to cover the whole EGLN1 genomic region in Andean highlanders living in Bolivia. Our findings indicated that haplotype frequencies estimated from the EGLN1 SNPs were significantly correlated with Hb concentration in the Bolivian highlanders. Moreover, we found that an Andean-dominant haplotype related to high Hb level may have expanded rapidly in ancestral Andean highlander populations. Analysis of genotype data in an ~436.3 kb genomic region containing EGLN1 using public databases indicated that the population structure based on EGLN1 genetic markers in Andean highlanders was largely different from that in other human populations. This finding may be related to an intrinsic or adaptive physiological characteristic of Andean highlanders. In conclusion, the high Hb concentrations in Andean highlanders can be partly characterized by EGLN1 genetic variants.

Individual variations and sex differences in hemodynamics with percutaneous arterial oxygen saturation (SpO2) in young Andean highlanders in Bolivia.

BACKGROUND: Many studies have reported specific adaptations to high altitude, but few studies have focused on physiological variations in high-altitude adaptation in Andean highlanders. This study aimed to investigate the relationships between SpO2 and related factors, including individual variations and sex differences, in Andean highlanders. METHODS: The participants were community-dwelling people in La Paz, Bolivia, aged 20 years and over (age range 20-34 years). A total of 50 men and 50 women participated in this study. Height, weight, SpO2, hemoglobin concentration, finger temperature, heart rate, and blood pressure were measured. Information about lifestyle was also obtained by interview. RESULTS: There were individual variations of SpO2 both in men (mean 89.9%, range 84.0-95.0%) and women (mean 91.0%, range 84.0-96.0%). On Student's t test, men had significantly lower heart rate (p = 0.046) and SpO2 (p = 0.030) than women. On the other hand, men had significantly higher SBP (p < 0.001), hemoglobin (p < 0.001), and finger temperature (p = 0.004). In men, multiple stepwise regression analysis showed that a higher SpO2 was correlated with a lower heart rate (ß = - 0.089, p = 0.007) and a higher finger temperature (ß = 0.308, p = 0.030) (r2 for model = 0.18). In women, a higher SpO2 was significantly correlated with a higher finger temperature (ß = 0.391, p = 0.015) (r2 for model = 0.12). A higher SpO2 was related to a higher finger temperature (ß = 0.286, p = 0.014) and a lower heart rate (ß = - 0.052, p = 0.029) in all participants (r2 for model = 0.21). Residual analysis showed that individual SpO2 values were randomly plotted. CONCLUSION: Random plots of SpO2 on residual analysis indicated that these variations were random error, such as biological variation. A higher SpO2 was related to a lower heart rate and finger temperature in men, but a higher SpO2 was related to finger temperature in women. These results suggest that there are individual variations and sex differences in the hemodynamic responses of high-altitude adaptation in Andean highlanders.

Transcriptomic Changes in Young Japanese Males After Exposure to Acute Hypobaric Hypoxia.

After the genomic era, the development of high-throughput sequencing technologies has allowed us to advance our understanding of genetic variants responsible for adaptation to high altitude in humans. However, transcriptomic characteristics associated with phenotypic plasticity conferring tolerance to acute hypobaric hypoxic stress remain unclear. To elucidate the effects of hypobaric hypoxic stress on transcriptional variability, we aimed to describe transcriptomic profiles in response to acute hypobaric hypoxia in humans. In a hypobaric hypoxic chamber, young Japanese males were exposed to a barometric pressure of 493 mmHg (hypobaric hypoxia) for 75 min after resting for 30 min at the pressure of 760 mmHg (normobaric normoxia) at 28°C. Saliva samples of the subjects were collected before and after hypobaric hypoxia exposure, to be used for RNA sequencing. Differential gene expression analysis identified 30 significantly upregulated genes and some of these genes may be involved in biological processes influencing hematological or immunological responses to hypobaric hypoxic stress. We also confirmed the absence of any significant transcriptional fluctuations in the analysis of basal transcriptomic profiles under no-stimulus conditions, suggesting that the 30 genes were actually upregulated by hypobaric hypoxia exposure. In conclusion, our findings showed that the transcriptional profiles of Japanese individuals can be rapidly changed as a result of acute hypobaric hypoxia, and this change may influence the phenotypic plasticity of lowland individuals for acclimatization to a hypobaric hypoxic environment. Therefore, the results obtained in this study shed light on the transcriptional mechanisms underlying high-altitude acclimatization in humans.

Effects of acute hypobaric hypoxia on thermoregulatory and circulatory responses during cold air exposure.

BACKGROUND: The thermoregulatory responses during simultaneous exposure to hypoxia and cold are not well understood owing to the opposite reactions of vasomotor tone in these two environments. Therefore, the purpose of this study was to investigate the influences of hypobaric hypoxia on various thermoregulatory responses, including skin blood flow (SkBF) during cold exposure. METHODS: Ten subjects participated in two experimental conditions: normobaric normoxia with cold (NC, barometric pressure (PB) = 760 mmHg) and hypobaric hypoxia with cold (HC, PB = 493 mmHg). The air temperature was maintained at 28 °C for 65 min and gradually decreased to 19 °C for both conditions. The total duration of the experiment was 135 min. RESULTS: The saturation of percutaneous oxygen (SpO2) was maintained at 98-99% in NC condition, but decreased to around 84% in HC condition. The rectal and mean skin temperatures showed no significant differences between the conditions; however, the forehead temperature was higher in HC condition than in NC condition. The pulse rate increased in HC condition, and there was a strong negative relationship between SpO2 and pulse rate (r = - 0.860, p = 0.013). SkBF and blood pressure showed no significant differences between the two conditions. CONCLUSION: These results suggest that hypobaric hypoxia during cold exposure did not alter the overall thermoregulatory responses. However, hypobaric hypoxia did affect pulse rate regardless of cold exposure.

Development of a novel monoclonal antibody that binds to most HLA-A allomorphs in a conformation-dependent yet peptide-promiscuous fashion.

Specificity analyses of peptide binding to human leukocyte antigen (HLA)-A molecules have been hampered due to a lack of proper monoclonal antibodies (mAbs) for certain allomorphs, such as the prevalent HLA-A1 for Caucasians and HLA-A11 for Asians. We developed a mAb that recognizes a conformational epitope common to most HLA-A allomorphs. The mAb, named A-1, does not discriminate peptides by amino acid sequences, making it suitable for measuring peptide binding. A stabilization assay using TAP-deficient cell lines and A-1 was developed to investigate the specificity of peptide binding to HLA-A molecules. Regarding the evolution of HLA-A genes, the A-1 epitope has been conserved among most HLA-A allomorphs but was lost when the HLA-A gene diversified into the HLA-A*32, HLA-A*31, and HLA-A*33 lineages together with HLA-A*29 after bifurcating from the HLA-A*25 and HLA-A*26 branchs. The establishment of A-1 is expected to help researchers investigate the peptide repertoire and develop computational tools to identify cognate peptides. Since no HLA-A locus-specific mAb has been available, A-1 will also be useful for analyzing the locus-specific regulation of the HLA gene expression.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies.

BACKGROUND: Our longitudinal exome-wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single-nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow-up period, 5 years) in 4884-6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease-susceptibility loci may help us uncover the pathogenesis of the related complex disorders. METHODS: In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease-susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project. RESULTS: Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out-of-Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non-Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease-associated loci may be affected by demographic events. CONCLUSION: Our findings indicate that derived allele frequencies of nine disease-associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome-wide association studies largely increased in non-Africans after Out-of-Africa.

Two novel susceptibility loci for type 2 diabetes mellitus identified by longitudinal exome-wide association studies in a Japanese population.

Recent genome-wide association studies identified genetic variants that confer susceptibility to type 2 diabetes mellitus (T2DM). However, few longitudinal genome-wide association studies of this metabolic disorder have been reported to date. Therefore, we performed a longitudinal exome-wide association study of T2DM, using 24,579 single nucleotide polymorphisms (SNPs) and repeated measurements from 6022 Japanese individuals. The generalized estimating equation model was applied to test relations of SNPs to three T2DM-related parameters: prevalence of T2DM, fasting plasma glucose level, and blood glycosylated hemoglobin content. Three SNPs that passed quality control were significantly (P<2.26×10-7) associated with two of the three T2DM-related parameters in additive and recessive models. Of the three SNPs, rs6414624 in EVC and rs78338345 in GGA3 were novel susceptibility loci for T2DM. In the present study, the SNP of GGA3 was predicted to be a genetic variant whose minor allele frequency has recently increased in East Asia.

Identification of six novel susceptibility loci for dyslipidemia using longitudinal exome-wide association studies in a Japanese population.

Recent genome-wide association studies have identified various dyslipidemia-related genetic variants. However, most studies were conducted in a cross-sectional manner. We thus performed longitudinal exome-wide association studies of dyslipidemia in a Japanese population. We used ~244,000 genetic variants and clinical data of 6022 Japanese individuals who had undergone annual health checkups for several years. After quality control, the association of dyslipidemia-related phenotypes with 24,691 single nucleotide polymorphisms (SNPs) was tested using the generalized estimating equation model. In total, 82 SNPs were significantly (P < 2.03 × 10-6) associated with dyslipidemia phenotypes. Of these SNPs, four (rs74416240 of TCHP, rs925368 of GIT2, rs7969300 of ATXN2, and rs12231744 of NAA25) and two (rs34902660 of SLC17A3 and rs1042127 of CDSN) were identified as novel genetic determinants of hypo-HDL- and hyper-LDL-cholesterolemia, respectively. A replication study using the cross-sectional data of 8310 Japanese individuals showed the association of the six identified SNPs with dyslipidemia-related traits.

Identification of 12 novel loci that confer susceptibility to early-onset dyslipidemia.

The circulating concentrations of triglycerides, high density lipoprotein (HDL)­cholesterol, and low density lipoprotein (LDL)­cholesterol have a substantial genetic component, and the heritability of early­onset dyslipidemia might be expected to be higher compared with late­onset forms. In the present study, exome­wide association studies (EWASs) were performed for early­onset hypertriglyceridemia, hypo­HDL­cholesterolemia, and hyper­LDL­cholesterolemia, with the aim to identify genetic variants that confer susceptibility to these conditions in the Japanese population. A total of 8,073 individuals aged ≤65 years were enrolled in the study. The EWASs for hypertriglyceridemia (2,664 cases and 5,294 controls), hypo­HDL­cholesterolemia (974 cases and 7,085 controls), and hyper­LDL­cholesterolemia (2,911 cases and 5,111 controls) were performed with Illumina Human Exome­12 v1.2 DNA Analysis BeadChip or Infinium Exome­24 v1.0 BeadChip arrays. The association of allele frequencies for 31,198, 31,133, or 31,175 single nucleotide polymorphisms (SNPs) to hypertriglyceridemia, hypo­HDL­cholesterolemia, or hyper­LDL­cholesterolemia, respectively, was examined with Fisher's exact test. To compensate for multiple comparisons of genotypes with each of the three conditions, Bonferroni's correction was applied for statistical significance of association. The results demonstrated that 25, 28 and 65 SNPs were significantly associated with hypertriglyceridemia, hypo­HDL­cholesterolemia and hyper­LDL­cholesterolemia, respectively. Multivariable logistic regression analysis with adjustment for age and sex revealed that all 25, 28 and 65 of these SNPs were significantly associated with hypertriglyceridemia, hypo­HDL­cholesterolemia and hyper­LDL­cholesterolemia, respectively. Following examination of the association of the identified SNPs to serum concentrations of triglycerides, HDL­cholesterol, or LDL­cholesterol, linkage disequilibrium of the SNPs, and results of previous genome­wide association studies, we newly identified chromosomal region 19p12 as a susceptibility locus for hypertriglyceridemia, eight loci (MOB3C­TMOD4, LPGAT1, EHD3, COL6A3, ZNF860­CACNA1D, COL6A5, DCLRE1C, ZNF77) for hypo­HDL­cholesterolemia, and three loci (KIAA0319­FAM65B, UBD, LOC105375015) for hyper­LDL­cholesterolemia. The present study thus identified 12 novel loci that may confer susceptibility to early­onset dyslipidemia. Determination of genotypes for the SNPs at these loci may prove informative for assessment of genetic risk for hypertriglyceridemia, hypo­HDL­cholesterolemia, or hyper­LDL­cholesterolemia in the Japanese population.

Identification of 26 novel loci that confer susceptibility to early-onset coronary artery disease in a Japanese population.

Early-onset coronary artery disease (CAD) has a strong genetic component. Although genome-wide association studies have identified various genes and loci significantly associated with CAD mainly in European populations, genetic variants that contribute toward susceptibility to this condition in Japanese patients remain to be definitively identified. In the present study, exome-wide association studies (EWASs) were performed to identify genetic variants that confer susceptibility to early-onset CAD in Japanese. A total of 7,256 individuals aged ≤65 years were enrolled in the present study. EWAS were conducted on 1,482 patients with CAD and 5,774 healthy controls. Genotyping of single nucleotide polymorphisms (SNPs) was performed using Illumina Human Exome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The association between allele frequencies for 31,465 SNPs that passed quality control and CAD was examined using Fisher's exact test. To compensate for multiple comparisons of allele frequencies with CAD, a false discovery rate (FDR) of <0.05 was applied for statistically significant associations. The association between allele frequencies for 31,465 SNPs and CAD, as determined by Fisher's exact test, demonstrated that 170 SNPs were significantly (FDR <0.05) associated with CAD. Multivariable logistic regression analysis with adjustment for age, sex, and the prevalence of hypertension, diabetes mellitus and dyslipidemia revealed that 162 SNPs were significantly (P<0.05) associated with CAD. A stepwise forward selection procedure was performed to examine the effects of genotypes for the 162 SNPs on CAD. The 54 SNPs were significant (P<0.05) and independent [coefficient of determination (R2), 0.0008 to 0.0297] determinants of CAD. These SNPs together accounted for 15.5% of the cause of CAD. Following examination of results from previous genome-wide association studies and linkage disequilibrium of the identified SNPs, 21 genes (RNF2, YEATS2, USP45, ITGB8, TNS3, FAM170B-AS1, PRKG1, BTRC, MKI67, STIM1, OR52E4, KIAA1551, MON2, PLUT, LINC00354, TRPM1, ADAT1, KRT27, LIPE, GFY and EIF3L) and five chromosomal regions (2p13, 4q31.2, 5q12, 13q34 and 20q13.2) that were significantly associated with CAD were newly identified in the present study. Gene ontology analysis demonstrated that various biological functions were predicted in the 18 genes identified in the present study. The network analysis revealed that the 18 genes had potential direct or indirect interactions with the 30 genes previously revealed to be associated with CAD or with the 228 genes identified in previous genome-wide association studies. The present study newly identified 26 loci that confer susceptibility to CAD. Determination of genotypes for the SNPs at these loci may prove informative for assessment of the genetic risk for CAD in Japanese patients.