Time-series transcriptome analysis of peripheral blood mononuclear cells obtained from individuals who received the SARS-CoV-2 mRNA vaccine.

Messenger ribonucleic acid (mRNA) vaccination against coronavirus disease 2019 (COVID-19) is an effective prevention strategy, despite a limited understanding of the molecular mechanisms underlying the host immune system and individual heterogeneity of the variable effects of mRNA vaccination. We assessed the time-series changes in the comprehensive gene expression profiles of 200 vaccinated healthcare workers by performing bulk transcriptome and bioinformatics analyses, including dimensionality reduction utilizing the uniform manifold approximation and projection (UMAP) technique. For these analyses, blood samples, including peripheral blood mononuclear cells (PBMCs), were collected from 214 vaccine recipients before vaccination (T1) and on Days 22 (T2, after second dose), 90, 180 (T3, before a booster dose), and 360 (T4, after a booster dose) after receiving the first dose of BNT162b2 vaccine (UMIN000043851). UMAP successfully visualized the main cluster of gene expression at each time point in PBMC samples (T1-T4). Through differentially expressed gene (DEG) analysis, we identified genes that showed fluctuating expression levels and gradual increases in expression levels from T1 to T4, as well as genes with increased expression levels at T4 alone. We also succeeded in dividing these cases into five types based on the changes in gene expression levels. High-throughput and temporal bulk RNA-based transcriptome analysis is a useful approach for inclusive, diverse, and cost-effective large-scale clinical studies.

Polycomb group protein BMI1 protects neuroblastoma cells against DNA damage-induced apoptotic cell death.

The overexpression of BMI1, a polycomb protein, correlates with cancer development and aggressiveness. We previously reported that MYCN-induced BMI1 positively regulated neuroblastoma (NB) cell proliferation via the transcriptional inhibition of tumor suppressors in NB cells. To assess the potential of BMI1 as a new target for NB therapy, we examined the effects of reductions in BMI1 on NB cells. BMI1 knockdown (KD) in NB cells significantly induced their differentiation for up to 7 days. BMI1 depletion significantly induced apoptotic NB cell death for up to 14 days along with the activation of p53, increases in p73, and induction of p53 family downstream molecules and pathways, even in p53 mutant cells. BMI1 depletion in vivo markedly suppressed NB xenograft tumor growth. BMI1 reductions activated ATM and increased γ-H2AX in NB cells. These DNA damage signals and apoptotic cell death were not canceled by the transduction of the polycomb group molecules EZH2 and RING1B. Furthermore, EZH2 and RING1B KD did not induce apoptotic NB cell death to the same extent as BMI1 KD. Collectively, these results suggest the potential of BMI1 as a target of molecular therapy for NB and confirmed, for the first time, the shared role of PcG proteins in the DNA damage response of NB cells.

Comparative epigenomics by machine learning approach for neuroblastoma.

BACKGROUND: Neuroblastoma (NB) is the second most common pediatric solid tumor. Because the number of genetic mutations found in tumors are small, even in some patients with unfavorable NB, epigenetic variation is expected to play an important role in NB progression. DNA methylation is a major epigenetic mechanism, and its relationship with NB prognosis has been a concern. One limitation with the analysis of variation in DNA methylation is the lack of a suitable analytical model. Therefore, in this study, we performed a random forest (RF) analysis of the DNA methylome data of NB from multiple databases. RESULTS: RF is a popular machine learning model owing to its simplicity, intuitiveness, and computational cost. RF analysis identified novel intermediate-risk patient groups with characteristic DNA methylation patterns within the low-risk group. Feature selection analysis based on probe annotation revealed that enhancer-annotated regions had strong predictive power, particularly for MYCN-amplified NBs. We developed a gene-based analytical model to identify candidate genes related to disease progression, such as PRDM8 and FAM13A-AS1. RF analysis revealed sufficient predictive power compared to other machine learning models. CONCLUSIONS: RF is a useful tool for DNA methylome analysis in cancer epigenetic studies, and has potential to identify a novel cancer-related genes.

Polycomb EZH1 regulates cell cycle/5-fluorouracil sensitivity of neuroblastoma cells in concert with MYCN.

In the present study, we found that EZH1 depletion in MYCN-amplified neuroblastoma cells resulted in significant cell death as well as xenograft inhibition. EZH1 depletion decreased the level of H3K27me1; the interaction and protein stabilization of MYCN and EZH1 appear to play roles in epigenetic transcriptional regulation. Transcriptome analysis of EZH1-depleted cells resulted in downregulation of the cell cycle progression-related pathway. In particular, Gene Set Enrichment Analysis revealed downregulation of reactome E2F-mediated regulation of DNA replication along with key genes of this process, TYMS, POLA2, and CCNA1. TYMS and POLA2 were transcriptionally activated by MYCN and EZH1-related epigenetic modification. Treatment with the EZH1/2 inhibitor UNC1999 also induced cell death, decreased H3K27 methylation, and reduced the levels of TYMS in neuroblastoma cells. Previous reports indicated neuroblastoma cells are resistant to 5-fluorouracil (5-FU) and TYMS (encoding thymidylate synthetase) has been considered the primary site of action for folate analogues. Intriguingly, UNC1999 treatment significantly sensitized MYCN-amplified neuroblastoma cells to 5-FU treatment, suggesting that EZH inhibition could be an effective strategy for development of a new epigenetic treatment for neuroblastoma.

The PRC2 molecule EED is a target of epigenetic therapy for neuroblastoma.

Epigenetic modifications by polycomb repressive complex (PRC) molecules appear to play a role in the tumorigenesis and aggressiveness of neuroblastoma (NB). Embryonic ectoderm development (EED) is a member of the PRC2 complex that binds to the H3K27me3 mark deposited by EZH2 via propagation on adjacent nucleosomes. We herein investigated the molecular roles of EED in MYCN-amplified NB cells using EED-knockdown (KD) shRNAs, EED-knockout sgRNAs, and the EED small molecule inhibitor EED226. The suppression of EED markedly inhibited NB cell proliferation and flat and soft agar colony formation. A transcriptome analysis using microarrays of EED-KD NB cells indicated the de-repression of cell cycle-regulated and differentiation-related genes. The results of a GSEA analysis suggested that inhibitory cell cycle-regulated gene sets were markedly up-regulated. Furthermore, an epigenetic treatment with the EED inhibitor EED226 and the HDAC inhibitors valproic acid/SAHA effectively suppressed NB cell proliferation and colony formation. This combined epigenetic treatment up-regulated cell cycle-regulated and differentiation-related genes. The ChIP sequencing analysis of histone codes and PRC molecules suggested an epigenetic background for the de-repression of down-regulated genes in MYCN-amplified/PRC2 up-regulated NB.

Loss of p53 suppresses replication stress-induced DNA damage in ATRX-deficient neuroblastoma.

Genetic aberrations are present in the ATRX gene in older high-risk neuroblastoma (NB) patients with very poor clinical outcomes. Its loss-of-function (LoF) facilitates the alternative lengthening of telomeres (ALT) pathway in tumor cells and is strongly linked to replication stress (RS) and DNA damage through G-quadruplex (G4) DNA secondary structures. However, limited information is available on ATRX alteration-related NB tumorigenesis. We herein knocked out (KO) ATRX in MYCN-amplified (NGP) and MYCN single copy (SK-N-AS) NB cells with wild-type (wt) and truncated TP53 at the C terminus, respectively, using CRISPR/Cas9 technologies. The loss of ATRX increased DNA damage and G4 formation related to RS in TP53 wt isogenic ATRX KO NGP cells, but not in SK-N-AS clones. A gene set enrichment analysis (GSEA) showed that the gene sets related to DNA double-strand break repair, negative cell cycle regulation, the G2M checkpoint, and p53 pathway activation were enriched in NGP clones. The accumulation of DNA damage activated the ATM/CHK2/p53 pathway, leading to cell cycle arrest in NGP clones. Interestingly, ATRX loss did not induce RS related to DNA damage response (DDR) in TP53-truncated SK-N-AS cells. p53 inactivation abrogated cell cycle arrest and reduced G4 accumulation in NGP clones. The loss of p53 also induced G4 DNA helicases or Fanconi anemia group D2 protein (FANCD2) with ATRX deficiency, suggesting that ATRX maintained genome integrity and p53 deficiency attenuated RS-induced DNA damage in NB cells featuring inactivated ATRX by regulating DNA repair mechanisms and replication fork stability.

Combined Genetic and Chromosomal Characterization of Wilms Tumors Identifies Chromosome 12 Gain as a Potential New Marker Predicting a Favorable Outcome.

To identify prognostic factors, array CGH (aCGH) patterns and mutations in WT1 and 9 other genes were analyzed in 128 unilateral Wilms tumors (WTs). Twenty patients had no aCGH aberrations, and 31 had WT1 alterations [silent and WT1 types: relapse-free survival (RFS), 95% and 83%, respectively]. Seventy-seven patients had aCGH changes without WT1 alterations (nonsilent/non-WT1 type) and were subtyped into those with or without +12, 11q-, 16q-, or HACE1 loss. RFS was better for those with than those without +12 (P = .010) and worse for those with than those without 11q-, 16q-, or HACE1 loss (P = .001, .025, or 1.2E-04, respectively). Silent and WT1 type and 8 subtype tumors were integrated and classified into 3 risk groups: low risk for the silent type and +12 subgroup; high risk for the no +12 plus 11q-, 16q-, or HACE1 loss subgroup; intermediate risk for the WT1 type and no +12 plus no 11q-, 16q-, or HACE1 loss subgroup. Among the 27 WTs examined, the expression of 146 genes on chromosome 12 was stronger in +12 tumors than in no +12 tumors, while that of 10 genes on 16q was weaker in 16q- tumors than in no 16q- tumors. Overexpression in 75 out of 146 upregulated genes and underexpression in 7 out of 10 downregulated genes correlated with better and worse overall survival, respectively, based on the public database. +12 was identified as a potential new marker predicting a favorable outcome, and chromosome abnormalities may be related to altered gene expression associated with these abnormalities.

EZH2 regulates neuroblastoma cell differentiation via NTRK1 promoter epigenetic modifications.

The polycomb repressor complex 2 molecule EZH2 is now known to play a role in essential cellular processes, namely, cell fate decisions, cell cycle regulation, senescence, cell differentiation, and cancer development/progression. EZH2 inhibitors have recently been developed; however, their effectiveness and underlying molecular mechanisms in many malignancies have not yet been elucidated in detail. Although the functional role of EZH2 in tumorigenesis in neuroblastoma (NB) has been investigated, mutations of EZH2 have not been reported. A Kaplan-Meier analysis on the event free survival and overall survival of NB patients indicated that the high expression of EZH2 correlated with an unfavorable prognosis. In order to elucidate the functional roles of EZH2 in NB tumorigenesis and its aggressiveness, we knocked down EZH2 in NB cell lines using lentivirus systems. The knockdown of EZH2 significantly induced NB cell differentiation, e.g., neurite extension, and the neuronal differentiation markers, NF68 and GAP43. EZH2 inhibitors also induced NB cell differentiation. We performed a comprehensive transcriptome analysis using Human Gene Expression Microarrays and found that NTRK1 (TrkA) is one of the EZH2-related suppression targets. The depletion of NTRK1 canceled EZH2 knockdown-induced NB cell differentiation. Our integrative methylome, transcriptome, and chromatin immunoprecipitation assays using NB cell lines and clinical samples clarified that the NTRK1 P1 and P2 promoter regions were regulated differently by DNA methylation and EZH2-related histone modifications. The NTRK1 transcript variants 1/2, which were regulated by EZH2-related H3K27me3 modifications at the P1 promoter region, were strongly expressed in favorable, but not unfavorable NB. The depletion and inhibition of EZH2 successfully induced NTRK1 transcripts and functional proteins. Collectively, these results indicate that EZH2 plays important roles in preventing the differentiation of NB cells and also that EZH2-related NTRK1 transcriptional regulation may be the key pathway for NB cell differentiation.

Sweepstake evolution revealed by population-genetic analysis of copy-number alterations in single genomes of breast cancer.

Single-cell sequencing is a promising technology that can address cancer cell evolution by identifying genetic alterations in individual cells. In a recent study, genome-wide DNA copy numbers of single cells were accurately quantified by single-cell sequencing in breast cancers. Phylogenetic-tree analysis revealed genetically distinct populations, each consisting of homogeneous cells. Bioinformatics methods based on population genetics should be further developed to quantitatively analyse the single-cell sequencing data. We developed a bioinformatics framework that was combined with molecular-evolution theories to analyse copy-number losses. This analysis revealed that most deletions in the breast cancers at the single-cell level were generated by simple stochastic processes. A non-standard type of coalescent theory, the multiple-merger coalescent model, aided by approximate Bayesian computation fit well with the data, allowing us to estimate the population-genetic parameters in addition to false-positive and false-negative rates. The estimated parameters suggest that the cancer cells underwent sweepstake evolution, where only one or very few parental cells produced a descendent cell population. We conclude that breast cancer cells successively substitute in a tumour mass, and the high reproduction of only a portion of cancer cells may confer high adaptability to this cancer.

CFC1 is a cancer stemness-regulating factor in neuroblastoma.

BACKGROUND: Despite the use of aggressive therapy, survival rates among high-risk neuroblastoma (NB) patients remain poor. Cancer stem cells (CSCs) are considered to be critically involved in the recurrence and metastasis of NB and are isolated as NB spheres. METHODS: The gene expression profiling of adherent (control) and sphere-forming primary NB cells was conducted using a gene expression microarray. CFC1, which functions in the development of embryos and decides the left-right axis, was strongly expressed in sphere-forming cells only and was related to the unfavorable prognosis of NB patients. The knockdown and overexpression of CFC1 were performed using a lentiviral system in NB cell lines. Sphere formation, cell proliferation, colony formation in soft agar, and xenograft tumor formation were analyzed. RESULTS: The overexpression of CFC1 increased sphere formation, cell growth, and colony formation. These phenotypes, particularly sphere formation, and xenograft tumor formation were significantly suppressed by the knockdown of CFC1. CFC1 inhibited Activin A-induced NB cell differentiation and Smad2 phosphorylation in NB cell lines, indicating its involvement in tumorigenesis related to EGF-CFC co-receptor family molecule pathways. Collectively, these results indicate that CFC1 is a candidate molecule for the development of CSC-targeted therapy for NB.

BRCA1 alterations with additional defects in DNA damage response genes may confer chemoresistance to BRCA-like breast cancers treated with neoadjuvant chemotherapy.

The BRCA-like phenotype is a feature that some sporadic breast cancers share with those occurring in BRCA1 or BRCA2 mutation carriers. As tumors with the phenotype have defects in the DNA damage response pathway, which may increase sensitivity to drugs such as DNA cross-linking agents and PARP inhibitors, a method to identify this phenotype is important. The prediction of chemoresistance, which frequently develops in these tumors, is also crucial for improving therapy. We examined genomic aberrations and BRCA1 promoter methylation in tumors of 73 breast cancer (20 HR-/HER2- and 53 HR+/HER2-) patients, who received neoadjuvant chemotherapy with anthracycline, cyclophosphamide, and taxane, using SNP array CGH and quantitative PCR. The methylation and/or loss or uniparental disomy (UPD) of BRCA1 (BRCA1 alterations) and the loss or UPD of BRCA2 (BRCA2 alterations) were detected in 27 (37%) and 21 (29%), respectively, of the 73 tumors. Tumors with BRCA1 or BRCA2 alterations were associated with a higher number of genomic aberrations (P < 0.001 and P < 0.001) and higher percentage of TP53 alterations (P < 0.001 and P < 0.001) than those without. Overall survival (OS) rates were similar between patients with or without BRCA1 or BRCA2 alterations. However, when 27 patients with BRCA1-altered tumors were classified into those with or without the loss or UPD of PALB2, PAGR1, RAD51B, FANCM, MLL4, or ERCC1/2 in tumors, patients with additional defects in DNA damage response genes had worse OS (P = 0.037, 0.045, 0.038, 0.044, 0.041, or 0.019) than those without. These defects may confer chemoresistance and predict poor outcomes in patients with BRCA1-altered breast cancer.

Population genomics of the fission yeast Schizosaccharomyces pombe.

The fission yeast Schizosaccharomyces pombe has been widely used as a model eukaryote to study a diverse range of biological processes. However, population genetic studies of this species have been limited to date, and we know very little about the evolutionary processes and selective pressures that are shaping its genome. Here, we sequenced the genomes of 32 worldwide S. pombe strains and examined the pattern of polymorphisms across their genomes. In addition to introns and untranslated regions (UTRs), intergenic regions also exhibited lower levels of nucleotide diversity than synonymous sites, suggesting that a considerable amount of noncoding DNA is under selective constraint and thus likely to be functional. A number of genomic regions showed a reduction of nucleotide diversity probably caused by selective sweeps. We also identified a region close to the end of chromosome 3 where an extremely high level of divergence was observed between 5 of the 32 strains and the remain 27, possibly due to introgression, strong positive selection, or that region being responsible for reproductive isolation. Our study should serve as an important starting point in using a population genomics approach to further elucidate the biology of this important model organism.

QTL map meets population genomics: an application to rice.

Genes involved in the transition from wild to cultivated crop species should be of great agronomic importance. Population genomic approaches utilizing genome resequencing data have been recently applied for this purpose, although it only reports a large list of candidate genes with no biological information. Here, by resequencing more than 30 genomes altogether of wild rice Oryza rufipogon and cultivated rice O. sativa, we identified a number of regions with clear footprints of selection during the domestication process. We then focused on identifying candidate domestication genes in these regions by utilizing the wealth of QTL information in rice. We were able to identify a number of interesting candidates such as transcription factors that should control key domestication traits such as shattering, awn length, and seed dormancy. Other candidates include those that might have been related to the improvement of grain quality and those that might have been involved in the local adaptation to dry conditions and colder environments. Our study shows that population genomic approaches and QTL mapping information can be used together to identify genes that might be of agronomic importance.

The linkage method: a novel approach for SNP detection and haplotype reconstruction from a single diploid individual using next-generation sequence data.

When we sequence a diploid individual, the output actually comprises two genomes: one from the paternal parent and the other from the maternal parent. In this study, we introduce a novel heuristic algorithm for distinguishing single-nucleotide polymorphisms (SNPs) from the two parents and phasing them into haplotypes. The algorithm is unique because it simultaneously performs SNP calling and haplotype phasing. This approach can exploit the linkage information of nearby SNPs, which facilitates the efficient removal of haplotypes that originate from incorrectly mapped short reads. Using simulated data we demonstrated that our approach increased the accuracy of SNP calls. The haplotype reconstruction performance depended largely on the density of SNPs. Using current next-generation sequence technology with a relatively short read length, reasonable performance is expected when this approach is applied to species with an average of five heterozygous sites per 1 kb. The algorithm was implemented as the program "linkSNPs."

Natural selection on gene order in the genome reorganization process after whole-genome duplication of yeast.

A genome must locate its coding genes on the chromosomes in a meaningful manner with the help of natural selection, but the mechanism of gene order evolution is poorly understood. To explore the role of selection in shaping the current order of coding genes and their cis-regulatory elements, a comparative genomic approach was applied to the baker's yeast Saccharomyces cerevisiae and its close relatives. S. cerevisiae have experienced a whole-genome duplication followed by an extensive reorganization process of gene order, during which a number of new adjacent gene pairs appeared. We found that the proportion of new adjacent gene pairs in divergent orientation is significantly reduced, suggesting that such new divergent gene pairs may be disfavored most likely because their coregulation may be deleterious. It is also found that such new divergent gene pairs have particularly long intergenic regions. These observations suggest that selection specifically worked against deletions in intergenic regions of new divergent gene pairs, perhaps because they should be physically kept away so that they are not coregulated. It is indicated that gene regulation would be one of the major factors to determine the order of coding genes.

Population genomics in bacteria: a case study of Staphylococcus aureus.

We analyzed the genome-wide pattern of single nucleotide polymorphisms (SNPs) in a sample with 12 strains of Staphylococcus aureus. Population structure of S. aureus seems to be complex, and the 12 strains were divided into five groups, named A, B, C, D, and E. We conducted a detailed analysis of the topologies of gene genealogies across the genomes and observed a high rate and frequency of tree-shape switching, indicating extensive homologous recombination. Most of the detected recombination occurred in the ancestral population of A, B, and C, whereas there are a number of small regions that exhibit evidence for homologous recombination with a distinct related species. As such regions would contain a number of novel mutations, it is suggested that homologous recombination would play a crucial role to maintain genetic variation within species. In the A-B-C ancestral population, we found multiple lines of evidence that the coalescent pattern is very similar to what is expected in a panmictic population, suggesting that this population is suitable to apply the standard population genetic theories. Our analysis showed that homologous recombination caused a dramatic decay in linkage disequilibrium (LD) and there is almost no LD between SNPs with distance more than 10 kb. Coalescent simulations demonstrated that a high rate of homologous recombination-a relative rate of 0.6 to the mutation rate with an average tract length of about 10 kb-is required to produce patterns similar to those observed in the S. aureus genomes. Our results call for more research into the evolutionary role of homologous recombination in bacterial populations.

Population Genomics of sub-saharan Drosophila melanogaster: African diversity and non-African admixture.

Drosophila melanogaster has played a pivotal role in the development of modern population genetics. However, many basic questions regarding the demographic and adaptive history of this species remain unresolved. We report the genome sequencing of 139 wild-derived strains of D. melanogaster, representing 22 population samples from the sub-Saharan ancestral range of this species, along with one European population. Most genomes were sequenced above 25X depth from haploid embryos. Results indicated a pervasive influence of non-African admixture in many African populations, motivating the development and application of a novel admixture detection method. Admixture proportions varied among populations, with greater admixture in urban locations. Admixture levels also varied across the genome, with localized peaks and valleys suggestive of a non-neutral introgression process. Genomes from the same location differed starkly in ancestry, suggesting that isolation mechanisms may exist within African populations. After removing putatively admixed genomic segments, the greatest genetic diversity was observed in southern Africa (e.g. Zambia), while diversity in other populations was largely consistent with a geographic expansion from this potentially ancestral region. The European population showed different levels of diversity reduction on each chromosome arm, and some African populations displayed chromosome arm-specific diversity reductions. Inversions in the European sample were associated with strong elevations in diversity across chromosome arms. Genomic scans were conducted to identify loci that may represent targets of positive selection within an African population, between African populations, and between European and African populations. A disproportionate number of candidate selective sweep regions were located near genes with varied roles in gene regulation. Outliers for Europe-Africa F(ST) were found to be enriched in genomic regions of locally elevated cosmopolitan admixture, possibly reflecting a role for some of these loci in driving the introgression of non-African alleles into African populations.

FancJ/Brip1 helicase protects against genomic losses and gains in vertebrate cells.

Defects in the FANCJ/BRIP1 helicase gene are associated with genome instability disorders such as familial breast cancer or Fanconi anemia (FA). Although FANCJ has an in vitro activity to resolve G-quadruplex (G4) structures, and FANCJ ortholog in C. elegans prevents G4-associated deletions during replication, how FANCJ loss affects genome integrity in higher organisms remains unclear. Here, we report that FANCJ, but not other FA genes FANCD2 or FANCC, protected against large-scale genomic deletion that occurred frequently at the rearranged immunoglobulin heavy chain (IgH) locus in chicken DT40 cell line, suggesting that FancJ protects the genome independently of the FA ubiquitination pathway. In a more unbiased approach using array-comparative genomic hybridization, we identified de novo deletions as well as amplifications in fancj cells kept in culture for 2 months. A cluster of G4 sequence motifs was found near the breakpoint of one amplified region, but G4 sequence motifs were not detected at the breakpoints of two deleted regions. These results collectively suggest that, unlike in C. elegans, actions of vertebrate FANCJ to promote genome stability may not be limited to protection against the G4-mediated gene deletions.

Selection for more of the same product as a force to enhance concerted evolution of duplicated genes.

The duration of concerted evolution after gene duplication is highly variable across genes. To identify the cause of the variation, we analyzed of duplicated genes in yeast that originate from a whole genome duplication event. There appears to be a strong positive correlation between the duration of concerted evolution and the gene expression level. This observation can be explained by selection favoring more of the same product, which could enhance concerted evolution in dosage-sensitive genes.

Estimating the time to the whole-genome duplication and the duration of concerted evolution via gene conversion in yeast.

A maximum-likelihood (ML) method is developed to estimate the duration of concerted evolution and the time to the whole-genome duplication (WGD) event in baker's yeast (Saccharomyces cerevisiae). The models with concerted evolution fit the data significantly better than the molecular clock model, indicating a crucial role of concerted evolution via gene conversion after gene duplication in yeast. Our ML estimate of the time to the WGD is nearly identical to the time to the speciation event between S. cerevisiae and Kluyveromyces waltii, suggesting that the WGD occurred in very early stages after speciation or the WGD might have been involved in the speciation event.