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.

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.

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.

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.

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 nine novel loci related to hematological traits in a Japanese population.

Recent genome-wide association studies have identified various genetic variants associated with hematological traits. Although it is possible that quantitative data of hematological traits are varied among the years examined, conventional genome-wide association studies have been conducted in a cross-sectional manner that measures traits at a single point in time. To address this issue, we have traced blood profiles in 4,884 Japanese individuals who underwent annual health check-ups for several years. In the present study, longitudinal exome-wide association studies were conducted to identify genetic variants related to 13 hematological phenotypes. The generalized estimating equation model showed that a total of 67 single nucleotide polymorphisms (SNPs) were significantly [false discovery rate (FDR) of <0.01] associated with hematological phenotypes. Of the 67 SNPs, nine SNPs were identified as novel hematological markers: rs4686683 of SENP2 for red blood cell count (FDR = 0.008, P = 5.5 × 10-6); rs3917688 of SELP for mean corpuscular volume (FDR = 0.005, P = 2.4 × 10-6); rs3133745 of C8orf37-AS1 for white blood cell count (FDR = 0.003, P = 1.3 × 10-6); rs13121954 at 4q31.2 for basophil count (FDR = 0.007, P = 3.1 × 10-5); rs7584099 at 2q22.3 (FDR = 2.6 × 10-5, P = 8.8 × 10-8), rs1579219 of HCG17 (FDR = 0.003, P = 2.0 × 10-5), and rs10757049 of DENND4C (FDR = 0.008, P = 5.6 × 10-5) for eosinophil count; rs12338 of CTSB for neutrophil count (FDR = 0.007, P = 2.9 × 10-5); and rs395967 of OSMR-AS1 for monocyte count (FDR = 0.008, P = 3.2 × 10-5).

Identification of three genetic variants as novel susceptibility loci for body mass index in a Japanese population.

Recent genome-wide association studies have identified various obesity or metabolic syndrome (MetS) susceptibility loci. However, most studies were conducted in a cross-sectional manner. To address this gap, we performed a longitudinal exome-wide association study to identify susceptibility loci for obesity and MetS in a Japanese population. We traced clinical data of 6,022 Japanese subjects who had annual health check-ups for several years (mean follow-up period, 5 yr) and genotyped ~244,000 genetic variants. The association of single nucleotide polymorphisms (SNPs) with body mass index (BMI) or the prevalence of obesity and MetS was examined in a generalized estimating equation model. Our longitudinal exome-wide association studies detected 21 BMI- and five MetS-associated SNPs (false discovery rate, FDR <0.01). Among these SNPs, 16 have not been previously implicated as determinants of BMI or MetS. Cross-sectional data for obesity- and MetS-related phenotypes in 7,285 Japanese subjects were examined in a replication study. Among the 16 SNPs, three ( rs9491140 , rs145848316 , and rs7863248 ) were related to BMI in the replication cohort ( P < 0.05). In conclusion, three SNPs [ rs9491140 of NKAIN2 (FDR = 0.003, P = 1.9 × 10-5), rs145848316 of KMT2C (FDR = 0.007, P = 4.5 × 10-5), and rs7863248 of AGTPBP1 (FDR = 0.006, P = 4.2 × 10-5)] were newly identified as susceptibility loci for BMI.

Identification of CDC42BPG as a novel susceptibility locus for hyperuricemia in a Japanese population.

Chronic kidney disease and hyperuricemia are serious global health problems. Recent genome-wide association studies have identified various genetic variants related to these disorders. However, most studies have been conducted in a cross-sectional manner. To identify novel susceptibility loci for chronic kidney disease or hyperuricemia, we performed longitudinal exome-wide association studies (EWASs), using ~ 244,000 genetic variants and clinical data of Japanese individuals who had undergone annual health checkups for several years. After establishing quality controls, the association of renal function-related traits in 5648 subjects (excluding patients with dialysis and population outliers) with 24,579 single nucleotide variants (SNVs) for three genetic models (P < 3.39 × 10- 7) was tested using generalized estimating equation models. The longitudinal EWASs revealed novel relations of five SNVs to renal function-related traits. Cross-sectional data for renal function-related traits in 7699 Japanese subjects were examined in a replication study. Among the five SNVs, rs55975541 in CDC42BPG was significantly (P < 4.90 × 10- 4) related to the serum concentration of uric acid in the replication cohort. We also examined the SNVs detected in our longitudinal EWASs with the information on P values in GKDGEN meta-analysis data. Four SNVs in SLC15A2 were significantly associated with the estimated glomerular filtration rate in European ancestry populations, although these SNVs were related to the serum concentration of uric acid with borderline significance in our longitudinal EWASs. Our findings indicate that CDC42BPG may be a novel susceptibility locus for hyperuricemia.

Elucidating the origin of HLA-B*73 allelic lineage: Did modern humans benefit by archaic introgression?

A previous study reported that some of the human leukocyte antigen (HLA) alleles and haplotypes in present-day humans were acquired by admixture with archaic humans; specifically, an exceptionally diverged HLA-B*73 allele was proposed to be transmitted from Denisovans, although the DNA sequence of HLA-B*73 has not been detected in the Denisovan genome. Here, we argue against the hypothesis that HLA-B*73 introgressed from Denisovans into early modern humans. A phylogenetic analysis revealed that HLA-B*73:01 formed a monophyletic group with a chimpanzee MHC-B allele, strongly suggesting that the HLA-B*73 allelic lineage has been maintained in humans as well as in chimpanzees since the divergence of humans and chimpanzees. The global distribution of HLA-B*73 allele showed that the population frequency of HLA-B*73 in west Asia (0.24 %)-a possible site of admixture with Denisovans-is lower than that in Europe (0.72 %) and in south Asia (0.69 %). Furthermore, HLA-B*73 is not observed in Melanesia even though the Melanesian genome contains the highest proportion of Denisovan ancestry in present-day human populations. Single nucleotide polymorphisms in HLA-A*11-HLA-C*12:02 or HLA-A*11-C*15 haplotypes, one of which was assumed to be transmitted together with HLA-B*73 from Denisovans by the study of Abi-Rached and colleagues, were not differentiated from those in other HLA-A-C haplotypes in modern humans. These results do not support the introgression hypothesis. Thus, we conclude that it is highly likely that HLA-B*73 allelic lineage has been maintained in the direct ancestors of modern humans.

Genetic evidence for contribution of human dispersal to the genetic diversity of EBA-175 in Plasmodium falciparum.

BACKGROUND: The 175-kDa erythrocyte binding antigen (EBA-175) of Plasmodium falciparum plays a crucial role in merozoite invasion into human erythrocytes. EBA-175 is believed to have been under diversifying selection; however, there have been no studies investigating the effect of dispersal of humans out of Africa on the genetic variation of EBA-175 in P. falciparum. METHODS: The PCR-direct sequencing was performed for a part of the eba-175 gene (regions II and III) using DNA samples obtained from Thai patients infected with P. falciparum. The divergence times for the P. falciparum eba-175 alleles were estimated assuming that P. falciparum/Plasmodium reichenowi divergence occurred 6 million years ago (MYA). To examine the possibility of diversifying selection, nonsynonymous and synonymous substitution rates for Plasmodium species were also estimated. RESULTS: A total of 32 eba-175 alleles were identified from 131 Thai P. falciparum isolates. Their estimated divergence time was 0.13-0.14 MYA, before the exodus of humans from Africa. A phylogenetic tree for a large sequence dataset of P. falciparum eba-175 alleles from across the world showed the presence of a basal Asian-specific cluster for all P. falciparum sequences. A markedly more nonsynonymous substitutions than synonymous substitutions in region II in P. falciparum was also detected, but not within Plasmodium species parasitizing African apes, suggesting that diversifying selection has acted specifically on P. falciparum eba-175. CONCLUSIONS: Plasmodium falciparum eba-175 genetic diversity appeared to increase following the exodus of Asian ancestors from Africa. Diversifying selection may have played an important role in the diversification of eba-175 allelic lineages. The present results suggest that the dispersals of humans out of Africa influenced significantly the molecular evolution of P. falciparum EBA-175.

Current perspectives on the intensity of natural selection of MHC loci.

Polymorphism of genes in the major histocompatibility complex (MHC) is believed to be maintained by balancing selection. However, direct evidence of selection has proven difficult to demonstrate. In 1994, Satta and colleagues estimated the selection intensity of the human MHC (human leukocyte antigen (HLA)) loci; however, at that time the number of HLA sequences was limited. By comparing five different methods, this study demonstrated the best way to calculate the selection coefficient, through a computer simulation study. Since the study, many HLA nucleotide sequences have been made available. Our new analysis takes advantage of these newly available sequences and compares new estimates with those of the previous study. Generally, our new results are consistent with those of the 1994 study. Our results show that, even after 20 years of exhaustive sequencing of human HLA, the number of dominant HLA alleles, on which our original estimate of selection intensity depended, appears to be conserved. Indeed, according to the frequency distribution for each HLA allele, most sequences in the database were minor or private alleles; therefore, we conclude that the selection intensities of HLA loci are at most 4.4 % even though the HLA is the prominent example on which the natural selection has been operating.

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.

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.

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.

MHC class II DQB diversity in the Japanese black bear, Ursus thibetanus japonicus.

BACKGROUND: The major histocompatibility complex (MHC) genes are one of the most important genetic systems in the vertebrate immune response. The diversity of MHC genes may directly influence the survival of individuals against infectious disease. However, there has been no investigation of MHC diversity in the Asiatic black bear (Ursus thibetanus). Here, we analyzed 270-bp nucleotide sequences of the entire exon 2 region of the MHC DQB gene by using 188 samples from the Japanese black bear (Ursus thibetanus japonicus) from 12 local populations. RESULTS: Among 185 of 188 samples, we identified 44 MHC variants that encoded 31 different amino acid sequences (allotypes) and one putative pseudogene. The phylogenetic analysis suggests that MHC variants detected from the Japanese black bear are derived from the DQB locus. One of the 31 DQB allotypes, Urth-DQB*01, was found to be common to all local populations. Moreover, this allotype was shared between the black bear on the Asian continent and the Japanese black bear, suggesting that Urth-DQB*01 might have been maintained in the ancestral black bear population for at least 300,000 years. Our findings, from calculating the ratio of non-synonymous to synonymous substitutions, indicate that balancing selection has maintained genetic variation of peptide-binding residues at the DQB locus of the Japanese black bear. From examination of genotype frequencies among local populations, we observed a considerably lower level of observed heterozygosity than expected. CONCLUSIONS: The low level of observed heterozygosity suggests that genetic drift reduced DQB diversity in the Japanese black bear due to a bottleneck event at the population or species level. The decline of DQB diversity might have been accelerated by the loss of rare variants that have been maintained by negative frequency-dependent selection. Nevertheless, DQB diversity of the black bear appears to be relatively high compared with some other endangered mammalian species. This result suggests that the Japanese black bears may also retain more potential resistance against pathogens than other endangered mammalian species. To prevent further decline of potential resistance against pathogens, a conservation policy for the Japanese black bear should be designed to maintain MHC rare variants in each local population.

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.

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.