Estimating gene-level false discovery probability improves eQTL statistical fine-mapping precision.

Statistical fine-mapping prioritizes putative causal variants from a large number of candidate variants, and is widely used in expression quantitative loci (eQTLs) studies. In eQTL fine-mapping, the existence of causal variants for gene expression is not guaranteed, since the genetic heritability of gene expression explained by nearby (cis-) variants is limited. Here we introduce a refined fine-mapping algorithm, named Knockoff-Finemap combination (KFc). KFc estimates the probability that the causal variant(s) exist in the cis-window of a gene through construction of knockoff genotypes (i.e. a set of synthetic genotypes that resembles the original genotypes), and uses it to adjust the posterior inclusion probabilities (PIPs). Utilizing simulated gene expression data, we show that KFc results in calibrated PIP distribution with improved precision. When applied to gene expression data of 465 genotyped samples from the Japan COVID-19 Task Force (JCTF), KFc resulted in significant enrichment of a functional score as well as reporter assay hits in the top PIP bins. When combined with functional priors derived from an external fine-mapping study (GTEx), KFc resulted in a significantly higher proportion of hematopoietic trait putative causal variants in the top PIP bins. Our work presents improvements in the precision of a major fine-mapping algorithm.

Subclonal accumulation of immune escape mechanisms in microsatellite instability-high colorectal cancers.

BACKGROUND: Intratumor heterogeneity (ITH) in microsatellite instability-high (MSI-H) colorectal cancer (CRC) has been poorly studied. We aimed to clarify how the ITH of MSI-H CRCs is generated in cancer evolution and how immune selective pressure affects ITH. METHODS: We reanalyzed public whole-exome sequencing data on 246 MSI-H CRCs. In addition, we performed a multi-region analysis from 6 MSI-H CRCs. To verify the process of subclonal immune escape accumulation, a novel computational model of cancer evolution under immune pressure was developed. RESULTS: Our analysis presented the enrichment of functional genomic alterations in antigen-presentation machinery (APM). Associative analysis of neoantigens indicated the generation of immune escape mechanisms via HLA alterations. Multiregion analysis revealed the clonal acquisition of driver mutations and subclonal accumulation of APM defects in MSI-H CRCs. Examination of variant allele frequencies demonstrated that subclonal mutations tend to be subjected to selective sweep. Computational simulations of tumour progression with the interaction of immune cells successfully verified the subclonal accumulation of immune escape mutations and suggested the efficacy of early initiation of an immune checkpoint inhibitor (ICI) -based treatment. CONCLUSIONS: Our results demonstrate the heterogeneous acquisition of immune escape mechanisms in MSI-H CRCs by Darwinian selection, providing novel insights into ICI-based treatment strategies.

Vector Autoregression for Forecasting the Number of COVID-19 Cases and Analyzing Behavioral Indicators in the Philippines: Ecologic Time-Trend Study.

BACKGROUND: Traditional surveillance systems rely on routine collection of data. The inherent delay in retrieval and analysis of data leads to reactionary rather than preventive measures. Forecasting and analysis of behavior-related data can supplement the information from traditional surveillance systems. OBJECTIVE: We assessed the use of behavioral indicators, such as the general public's interest in the risk of contracting SARS-CoV-2 and changes in their mobility, in building a vector autoregression model for forecasting and analysis of the relationships of these indicators with the number of COVID-19 cases in the National Capital Region. METHODS: An etiologic, time-trend, ecologic study design was used to forecast the daily number of cases in 3 periods during the resurgence of COVID-19. We determined the lag length by combining knowledge on the epidemiology of SARS-CoV-2 and information criteria measures. We fitted 2 models to the training data set and computed their out-of-sample forecasts. Model 1 contains changes in mobility and number of cases with a dummy variable for the day of the week, while model 2 also includes the general public's interest. The forecast accuracy of the models was compared using mean absolute percentage error. Granger causality test was performed to determine whether changes in mobility and public's interest improved the prediction of cases. We tested the assumptions of the model through the Augmented Dickey-Fuller test, Lagrange multiplier test, and assessment of the moduli of eigenvalues. RESULTS: A vector autoregression (8) model was fitted to the training data as the information criteria measures suggest the appropriateness of 8. Both models generated forecasts with similar trends to the actual number of cases during the forecast period of August 11-18 and September 15-22. However, the difference in the performance of the 2 models became substantial from January 28 to February 4, as the accuracy of model 2 remained within reasonable limits (mean absolute percentage error [MAPE]=21.4%) while model 1 became inaccurate (MAPE=74.2%). The results of the Granger causality test suggest that the relationship of public interest with number of cases changed over time. During the forecast period of August 11-18, only change in mobility (P=.002) improved the forecasting of cases, while public interest was also found to Granger-cause the number of cases during September 15-22 (P=.001) and January 28 to February 4 (P=.003). CONCLUSIONS: To the best of our knowledge, this is the first study that forecasted the number of COVID-19 cases and explored the relationship of behavioral indicators with the number of COVID-19 cases in the Philippines. The resemblance of the forecasts from model 2 with the actual data suggests its potential in providing information about future contingencies. Granger causality also implies the importance of examining changes in mobility and public interest for surveillance purposes.

The current status of emergency departments in secondary emergency medical institutions in Japan: a questionnaire survey.

BACKGROUND: While emergency medicine (ER)-based emergency care is prevalent in many countries, in Japan, the "department-specific emergency care model" and the "emergency center model" are mainstream. We hypothesized that many secondary emergency medical institutions in Japan have inadequate systems. Using a questionnaire, we investigated the status of and problems in the emergency medical services system in secondary emergency medical institutions in Japan. Until date, there has not been an exhaustive survey of emergency facilities on a countrywide scale. The main objective of this study was to investigate problems in the Japanese emergency medical services system and thereby improve optimal care for emergency patients. RESULTS: A nationwide questionnaire survey involving 4063 facilities (all government-approved emergency medical facilities certified by prefectural governors) in Japan was conducted. Of the facilities that responded, all secondary emergency facilities were included in the analysis. Responses from 1289 facilities without a tertiary emergency medical care center were analyzed. Among them, 61% (792/1289) had ≤ 199 beds, and 8% were emergency department specialty training program core facilities. Moreover, 42% had an annual patient acceptance number of ≤ 500, 19% did not calculate the number of acceptances, 29% had an acceptance rate of ≥ 81%, and 25% had an acceptance rate of 61-80%. Pregnant women (63%) and children (56%) were the major types of patients that affected the acceptance rate. Factors affecting facilities with a response rate of 81% or higher were "hospitals designated for residency training" and "facilities making some efforts to improve the quality of emergency care and the emergency medical system" (logistic analysis, P < .001). CONCLUSION: Relevant authorities and core regional facilities should consider and implement specific measures for regions and hospitals with a shortage of emergency medicine specialists and physicians (e.g., development of ER-based emergency medicine and provision of education). This study may lead to further improvement in the optimal care of emergency patients through the nationwide establishment of the proposed measures as well as through grouping and integrating the structures and systems in emergency and other medical facilities.

Convergent genomic diversity and novel BCAA metabolism in intrahepatic cholangiocarcinoma.

BACKGROUND: Driver alterations may represent novel candidates for driver gene-guided therapy; however, intrahepatic cholangiocarcinoma (ICC) with multiple genomic aberrations makes them intractable. Therefore, the pathogenesis and metabolic changes of ICC need to be understood to develop new treatment strategies. We aimed to unravel the evolution of ICC and identify ICC-specific metabolic characteristics to investigate the metabolic pathway associated with ICC development using multiregional sampling to encompass the intra- and inter-tumoral heterogeneity. METHODS: We performed the genomic, transcriptomic, proteomic and metabolomic analysis of 39-77 ICC tumour samples and eleven normal samples. Further, we analysed their cell proliferation and viability. RESULTS: We demonstrated that intra-tumoral heterogeneity of ICCs with distinct driver genes per case exhibited neutral evolution, regardless of their tumour stage. Upregulation of BCAT1 and BCAT2 indicated the involvement of 'Val Leu Ile degradation pathway'. ICCs exhibit the accumulation of ubiquitous metabolites, such as branched-chain amino acids including valine, leucine, and isoleucine, to negatively affect cancer prognosis. We revealed that this metabolic pathway was almost ubiquitously altered in all cases with genomic diversity and might play important roles in tumour progression and overall survival. CONCLUSIONS: We propose a novel ICC onco-metabolic pathway that could enable the development of new therapeutic interventions.

The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force.

Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection.

DOCK2 is involved in the host genetics and biology of severe COVID-19.

Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge1-5. Here we conducted a genome-wide association study (GWAS) involving 2,393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3,289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target.

Impact of salivary and pancreatic amylase gene copy numbers on diabetes, obesity, and functional profiles of microbiome in Northern Japanese population.

Amylase genes reside in a structurally complex locus, and their copy numbers vary greatly, and several studies have reported their association with obesity. The mechanism of this effect was partially explained by changes in the oral and gut microbiome compositions; however, a detailed mechanism has been unclarified. In this study, we showed their association with diabetes in addition to obesity, and further discovered a plausible mechanism of this association based on the function of commensal bacteria. First, we confirmed that the amylase copy number in the population tends to be larger than that reported in other studies and that there is a positive association between obesity and diabetes (p = 1.89E-2 and 8.63E-3). Second, we identified that relative abundance of some genus level microbiome, Capnocytophaga, Dialister, and previously reported bacteria, were significantly associated with amylase copy numbers. Finally, through functional gene-set analysis using shotgun sequencing, we observed that the abundance of genes in the Acarbose pathway in the gut microbiome was significantly decreased with an increase in the amylase copy number (p-value = 5.80E-4). Our findings can partly explain the mechanism underlying obesity and diabetes in populations with high amylase copy numbers.

Protease-digested egg-white products induce oral tolerance in mice but elicit little IgE production upon epicutaneous exposure.

BACKGROUND: Early food introduction induces tolerance, but epicutaneous exposure, especially via eczema lesions, promotes IgE sensitization. Aiming for safe and effective primary prevention of egg allergy, we examined several protease-digested egg-white (EW) products for three properties: 1) induction of oral tolerance that prevents IgE sensitization, 2) weak IgE binding that can prevent allergic reactions even in IgE-sensitized mice, and 3) minimal epicutaneous IgE sensitization even when in contact with inflamed skin. METHODS: Heated EW was digested with several proteases under optimal conditions. First, three-week-old BALB/c female mice were intragastrically administered EW or each protease-digested EW product, followed by intraperitoneal ovalbumin (OVA) or ovomucoid (OVM) injection with alum. Serum OVA- and OVM-specific IgE titers were measured. Second, six-week-old mice were sensitized with OVA/OVM, and the rectal temperature was measured after intraperitoneal administration of EW or each protease-digested EW. Third, EW or each protease-digested EW product was applied to the tape-stripped skin for 3 days/week for 3 weeks. Serum OVA- and OVM-specific IgE titers were measured. RESULTS: Orally administered pepsin-digested EW product (PDEW) and Thermoase PC10F-digested EW product (TDEW) significantly suppressed OVA-/OVM-specific IgE production. Neither product elicited a body temperature decline (anaphylaxis) in OVA-/OVM-sensitized mice. Serum OVA-/OVM-specific IgE levels were significantly lower in mice epicutaneously exposed to PDEW or TDEW than in EW-exposed mice. CONCLUSIONS: Two protease-digested EWs showed potential as optimal EW products for early introduction for primary prevention of egg allergy.

Immunogenomic pan-cancer landscape reveals immune escape mechanisms and immunoediting histories.

Immune reactions in the tumor microenvironment are an important hallmark of cancer, and emerging immune therapies have been proven effective against several types of cancers. To investigate cancer genome-immune interactions and the role of immunoediting or immune escape mechanisms in cancer development, we analyzed 2834 whole genome and RNA sequencing datasets across 31 distinct tumor types with respect to key immunogenomic aspects and provided comprehensive immunogenomic profiles of pan-cancers. We found that selective copy number changes in immune-related genes may contribute to immune escape. Furthermore, we developed an index of the immunoediting history of each tumor sample based on the information of mutations in exonic regions and pseudogenes and evaluated the immunoediting history of each tumor. Our immuno-genomic analyses of pan-cancers have the potential to identify a subset of tumors with immunogenicity and diverse backgrounds or intrinsic pathways associated with their immune status and immunoediting history.

The Evolving Genomic Landscape of Esophageal Squamous Cell Carcinoma Under Chemoradiotherapy.

Esophageal squamous cell carcinoma (ESCC) often recurs after chemoradiotherapy, and the prognosis of ESCC after chemoradiotherapy has not improved over the past few decades. The mutation process in chemoradiotherapy-resistant clones and the functional relevance of genetic alterations remain unclear. To address these problems, we performed whole-exome sequencing of 52 tumor samples from 33 patients with ESCC who received radiotherapy combined with 5-fluorouracil/platinum. In multiregion analyses of pretreatment and locally recurrent lesions from five cases, most driver gene-altered clones remained under chemoradiotherapy selection pressure, while few driver gene alterations were acquired at recurrence. The mutation signatures of recurrent ESCC, including increased deletion frequency and platinum dose-dependent base substitution signatures, were substantially different from those of primary ESCC and reflected the iatrogenic impacts of chemoradiotherapy. Single-region analysis of 28 pretreatment tumors indicated that focal copy-number gain at the MYC locus was significantly associated with poor progression-free survival and overall survival after chemoradiotherapy. MYC gain remained throughout the chemoradiotherapy course and potentially contributes to intrinsic resistance to chemoradiotherapy. Consistent with these findings, MYC copy number and mRNA and protein levels in ESCC cell lines correlated positively with resistance to radiotherapy, and MYC knockdown improved sensitivity to radiotherapy. Overall, these data characterize the clonal evolution process induced by chemoradiotherapy and clinically relevant associations for genetic alterations in ESCC. These findings increase our understanding of therapeutic resistance and support the rationale for precision chemoradiotherapy. SIGNIFICANCE: Whole-exome sequencing reveals the genetic evolution of ESCC during chemoradiotherapy, highlighting MYC gain in pretreatment tumors as a potential marker of therapy resistance.

Deficiency of TMEM53 causes a previously unknown sclerosing bone disorder by dysregulation of BMP-SMAD signaling.

Bone formation represents a heritable trait regulated by many signals and complex mechanisms. Its abnormalities manifest themselves in various diseases, including sclerosing bone disorder (SBD). Exploration of genes that cause SBD has significantly improved our understanding of the mechanisms that regulate bone formation. Here, we discover a previously unknown type of SBD in four independent families caused by bi-allelic loss-of-function pathogenic variants in TMEM53, which encodes a nuclear envelope transmembrane protein. Tmem53-/- mice recapitulate the human skeletal phenotypes. Analyses of the molecular pathophysiology using the primary cells from the Tmem53-/- mice and the TMEM53 knock-out cell lines indicates that TMEM53 inhibits BMP signaling in osteoblast lineage cells by blocking cytoplasm-nucleus translocation of BMP2-activated Smad proteins. Pathogenic variants in the patients impair the TMEM53-mediated blocking effect, thus leading to overactivated BMP signaling that promotes bone formation and contributes to the SBD phenotype. Our results establish a previously unreported SBD entity (craniotubular dysplasia, Ikegawa type) and contribute to a better understanding of the regulation of BMP signaling and bone formation.

Impaired tumor immune response in metastatic tumors is a selective pressure for neutral evolution in CRC cases.

A Darwinian evolutionary shift occurs early in the neutral evolution of advanced colorectal carcinoma (CRC), and copy number aberrations (CNA) are essential in the transition from adenoma to carcinoma. In light of this primary evolution, we investigated the evolutionary principles of the genome that foster postoperative recurrence of CRC. CNA and neoantigens (NAG) were compared between early primary tumors with recurrence (CRCR) and early primary tumors without recurrence (precancerous and early; PCRC). We compared CNA, single nucleotide variance (SNV), RNA sequences, and T-cell receptor (TCR) repertoire between 9 primary and 10 metastatic sites from 10 CRCR cases. We found that NAG in primary sites were fewer in CRCR than in PCRC, while the arm level CNA were significantly higher in primary sites in CRCR than in PCRC. Further, a comparison of genomic aberrations of primary and metastatic conditions revealed no significant differences in CNA. The driver mutations in recurrence were the trunk of the evolutionary phylogenic tree from primary sites to recurrence sites. Notably, PD-1 and TIM3, T cell exhaustion-related molecules of the tumor immune response, were abundantly expressed in metastatic sites compared to primary sites along with the increased number of CD8 expressing cells. The postoperative recurrence-free survival period was only significantly associated with the NAG levels and TCR repertoire diversity in metastatic sites. Therefore, CNA with diminished NAG and diverse TCR repertoire in pre-metastatic sites may determine postoperative recurrence of CRC.

Metagenomic profiling of gut microbiome in early chronic kidney disease.

BACKGROUND: The relationship between chronic kidney disease (CKD) and the gut microbiome, which interact through chronic inflammation, uraemic toxin production and immune response regulation, has gained interest in the development of CKD therapies. However, reports using shotgun metagenomic analysis of the gut microbiome are scarce, especially for early CKD. Here we characterized gut microbiome differences between non-CKD participants and ones with early CKD using metagenomic sequencing. METHODS: In total, 74 non-CKD participants and 37 participants with early CKD were included based on propensity score matching, controlling for various factors including dietary intake. Stool samples were collected from participants and subjected to shotgun sequencing. Bacterial and pathway abundances were profiled at the species level with MetaPhlAn2 and HUMAnN2, respectively, and overall microbiome differences were determined using Bray-Curtis dissimilarities. Diabetic and non-diabetic populations were analysed separately. RESULTS: For diabetic and non-diabetic participants, the mean estimated glomerular filtration rates of the CKD group were 53.71 [standard deviation (SD) 3.87] and 53.72 (SD 4.44), whereas those of the non-CKD group were 72.63 (SD 7.72) and 76.10 (SD 9.84), respectively. Alpha and beta diversities were not significantly different between groups. Based on taxonomic analysis, butyrate-producing species Roseburia inulinivorans, Ruminococcus torques and Ruminococcus lactaris were more abundant in the non-CKD group, whereas Bacteroides caccae and Bacteroides coprocora were more abundant in the non-diabetic CKD group. CONCLUSIONS: Although gut microbiome changes in individuals with early CKD were subtle, the results suggest that changes related to producing short-chain fatty acids can already be observed in early CKD.

Data science and precision health care.

Precision health care plays a crucial role in an elderly society by providing personalized health care plans for improving an individual's health conditions and preventing disease. To realize precision health care, data science is key; it allows for analyses of health-related big data. In this article, an actual analysis of time-series health check-up data is presented and as is a discussion of how personalized simulation models of health conditions are constructed and used to modify individual behavior. Future directions for precision health care based on the integration of genetic variations and the microbiome are also discussed.

Neoantimon: a multifunctional R package for identification of tumor-specific neoantigens.

SUMMARY: It is known that some mutant peptides, such as those resulting from missense mutations and frameshift insertions, can bind to the major histocompatibility complex and be presented to antitumor T cells on the surface of a tumor cell. These peptides are termed neoantigen, and it is important to understand this process for cancer immunotherapy. Here, we introduce an R package termed Neoantimon that can predict a list of potential neoantigens from a variety of mutations, which include not only somatic point mutations but insertions, deletions and structural variants. Beyond the existing applications, Neoantimon is capable of attaching and reflecting several additional information, e.g. wild-type binding capability, allele specific RNA expression levels, single nucleotide polymorphism information and combinations of mutations to filter out infeasible peptides as neoantigen. AVAILABILITY AND IMPLEMENTATION: The R package is available at http://github/hase62/Neoantimon.

Association of single nucleotide polymorphisms in the NRF2 promoter with vascular stiffness with aging.

PURPOSE: Pulse wave velocity (PWV), an indicator of vascular stiffness, increases with age and is increasingly recognized as an independent risk factor for cardiovascular disease (CVD). Although many mechanical and chemical factors underlie the stiffness of the elastic artery, genetic risk factors related to age-dependent increases in PWV in apparently healthy people are largely unknown. The transcription factor nuclear factor E2 (NF-E2)-related factor 2 (Nrf2), which is activated by unidirectional vascular pulsatile shear stress or oxidative stress, regulates vascular redox homeostasis. Previous reports have shown that a SNP in the NRF2 gene regulatory region (-617C>A; hereafter called SNP-617) affects NRF2 gene expression such that the minor A allele confers lower gene expression compared to the C allele, and it is associated with various diseases, including CVD. We aimed to investigate whether SNP-617 affects vascular stiffness with aging in apparently healthy people. METHODS: Analyzing wide-ranging data obtained from a public health survey performed in Japan, we evaluated whether SNP-617 affected brachial-ankle PWV (baPWV) in never-smoking healthy subjects (n = 642). We also evaluated the effects of SNP-617 on other cardiovascular and blood test measurements. RESULTS: We have shown that not only AA carriers (n = 55) but also CA carriers (n = 247) show arterial stiffness compared to CC carriers (n = 340). Furthermore, SNP-617 also affected blood pressure indexes such as systolic blood pressure and mean arterial pressure but not the ankle brachial pressure index, an indicator of atherosclerosis. Multivariate analysis showed that SNP-617 accelerates the incremental ratio of baPWV with age. CONCLUSIONS: This study is the first to show that SNP-617 affects the age-dependent increase in vascular stiffness. Our results indicate that low NRF2 activity induces premature vascular aging and could be targeted for the prevention of cardiovascular diseases associated with aging.

Tet DNA demethylase is required for plasma cell differentiation by controlling expression levels of IRF4.

Antibodies produced by plasma cells are critical for protection from infection. It has been demonstrated that global epigenetic modification, such as changes in DNA methylation, occurs during differentiation of plasma cells from B cells. However, the precise mechanisms by which DNA methylation controls plasma cell differentiation are not fully understood. We examined the effect of deficiency of DNA demethylases, Tet2 and Tet3, on B-cell activation and plasma cell differentiation, by generating conditional Tet2/3 double-KO (Tet dKO) B cells. We found that Tet dKO B cells failed to differentiate into plasma cells upon immunization with antigens. Tet dKO B cells proliferated normally and were capable of generating cells with IRF4int, but not with IRF4hi, the majority of which were CD138+ plasma cells. IRF4 overexpression rescued the defect of Tet dKO B cells in plasma cell differentiation, suggesting that Tet2/3-dependent high IRF4 expression is required for plasma cell differentiation. We identified CpG sites in the Irf4 locus that were demethylated specifically in plasma cells and in a Tet2/3-dependent manner. Our results suggest that Tet2/3-dependent demethylation of these CpG sites is dispensable for initial IRF4 expression but is essential for high IRF4 expression which is prerequisite for plasma cell differentiation.

The relationship between cigarette smoking and the tongue microbiome in an East Asian population.

Background: The oral microbiome, which consists of various habitats, has been shown to be influenced by smoking. However, differences in the tongue microbiomes of current and former smokers, as well as their resultant functional consequences, have rarely been investigated in East Asian populations. Methods: We used 16S rRNA amplicon sequencing of tongue-coating samples obtained from East Asian subjects who were current, former, or never smokers to identify differences in their tongue microbiomes and related metagenomic functions. Two sets of participants from 2016 to 2017 (n = 657 and n = 187, respectively) were analyzed separately. Results: We found significant differences between the overall microbiome compositions of current versus never smokers (p = 0.0015), but not between former versus never smokers (p = 0.43) based on the weighted UniFrac distance. Twenty-nine of 43 investigated genera showed significantly different expression levels in current versus never smokers. Neisseria and Capnocytophaga were less abundant, and Streptococcus and Megasphaera were more abundant in current smokers. Moreover, the abundances of metagenomic pathways, including those related to nitrate reduction and the tricarboxylic acid cycle, were significantly different between current and never smokers. Conclusions: The tongue microbiomes and related metagenomic pathways of current smokers differ from those of never smokers among East Asians.

A unified simulation model for understanding the diversity of cancer evolution.

Because cancer evolution underlies the therapeutic difficulties of cancer, it is clinically important to understand the evolutionary dynamics of cancer. Thus far, a number of evolutionary processes have been proposed to be working in cancer evolution. However, there exists no simulation model that can describe the different evolutionary processes in a unified manner. In this study, we constructed a unified simulation model for describing the different evolutionary processes and performed sensitivity analysis on the model to determine the conditions in which cancer growth is driven by each of the different evolutionary processes. Our sensitivity analysis has successfully provided a series of novel insights into the evolutionary dynamics of cancer. For example, we found that, while a high neutral mutation rate shapes neutral intratumor heterogeneity (ITH) characterized by a fractal-like pattern, a stem cell hierarchy can also contribute to shaping neutral ITH by apparently increasing the mutation rate. Although It has been reported that the evolutionary principle shaping ITH shifts from selection to accumulation of neutral mutations during colorectal tumorigenesis, our simulation revealed the possibility that this evolutionary shift is triggered by drastic evolutionary events that occur in a short time and confer a marked fitness increase on one or a few cells. This result helps us understand that each process works not separately but simultaneously and continuously as a series of phases of cancer evolution. Collectively, this study serves as a basis to understand in greater depth the diversity of cancer evolution.