An improved method for the highly specific detection of transcription start sites.

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

ABA signaling converts stem cell fate by substantiating a tradeoff between cell polarity, growth and cell cycle progression and abiotic stress responses in the moss Physcomitrium patens.

Abscisic acid (ABA)-mediated abiotic stress tolerance causes plant growth inhibition. Under such stress conditions, some mosses generate de novo stress-resistant stem cells, also called brood cells or brachycytes, that do not exist under normal conditions. However, the cell physiological basis of the growth inhibition and the stem cell formation is not well understood. Here, we show that the ABA-induced growth inhibition of the moss Physcomitrium patens apical protonemal cells (protonemal stem cells) is mediated through a shift from asymmetric to symmetric cell division. This change of the cell division mode, and consequently change of stem cell activity, is substantiated by dampening cell polarity and cell proliferative activity through the altered distribution of cytoskeletal elements, the mitotic spindle and the vacuole, which results in the production of stress-resistant stem cells. Alteration of the cell physiological data is supported by the results of RNAseq analysis indicating rapid changes in both cell polarity and cell cycle regulation, while long-term treatments with ABA for 5 to 10 days impact mainly the transcriptional and translational regulation. The regulation of cell polarity and cell cycle genes suggests growth arrest mediated by small GTPases (ROPs) and their guanine exchange factors (ROPGEFs) and by cyclin and cyclin-dependent-kinase complex, respectively. Our data suggest that a tradeoff relationship between growth ability and abiotic stress response in the moss is substantiated by ABA signaling to suppress cell polarity and asymmetric cell growth and may play a pivotal role in stem cell fate conversion to newly produced stress-resistant stem cells.

FXYD3 functionally demarcates an ancestral breast cancer stem cell subpopulation with features of drug-tolerant persisters.

The heterogeneity of cancer stem cells (CSCs) within tumors presents a challenge in therapeutic targeting. To decipher the cellular plasticity that fuels phenotypic heterogeneity, we undertook single-cell transcriptomics analysis in triple-negative breast cancer (TNBC) to identify subpopulations in CSCs. We found a subpopulation of CSCs with ancestral features that is marked by FXYD domain-containing ion transport regulator 3 (FXYD3), a component of the Na+/K+ pump. Accordingly, FXYD3+ CSCs evolve and proliferate, while displaying traits of alveolar progenitors that are normally induced during pregnancy. Clinically, FXYD3+ CSCs were persistent during neoadjuvant chemotherapy, hence linking them to drug-tolerant persisters (DTPs) and identifying them as crucial therapeutic targets. Importantly, FXYD3+ CSCs were sensitive to senolytic Na+/K+ pump inhibitors, such as cardiac glycosides. Together, our data indicate that FXYD3+ CSCs with ancestral features are drivers of plasticity and chemoresistance in TNBC. Targeting the Na+/K+ pump could be an effective strategy to eliminate CSCs with ancestral and DTP features that could improve TNBC prognosis.

Multimodal single-cell analyses of peripheral blood mononuclear cells of COVID-19 patients in Japan.

SARS-CoV-2 continues to spread worldwide. Patients with COVID-19 show distinct clinical symptoms. Although many studies have reported various causes for the diversity of symptoms, the underlying mechanisms are not fully understood. Peripheral blood mononuclear cells from COVID-19 patients were collected longitudinally, and single-cell transcriptome and T cell receptor repertoire analysis was performed. Comparison of molecular features and patients' clinical information revealed that the proportions of cells present, and gene expression profiles differed significantly between mild and severe cases; although even among severe cases, substantial differences were observed among the patients. In one severely-infected elderly patient, an effective antibody response seemed to have failed, which may have caused prolonged viral clearance. Naïve T cell depletion, low T cell receptor repertoire diversity, and aberrant hyperactivation of most immune cell subsets were observed during the acute phase in this patient. Through this study, we provided a better understanding of the diversity of immune landscapes and responses. The information obtained from this study can help medical professionals develop personalized optimal clinical treatment strategies for COVID-19.

Whole-genome Methylation Analysis of APOBEC Enzyme-converted DNA (~5 kb) by Nanopore Sequencing.

In recent years, DNA methylation research has been accelerated by the advent of nanopore sequencers. However, read length has been limited by the constraints of base conversion using the bisulfite method, making analysis of chromatin content difficult. The read length of the previous method combining bisulfite conversion and long-read sequencing was ~1.5 kb, even using targeted PCR. In this study, we have improved read length (~5 kb), by converting unmethylated cytosines to uracils with APOBEC enzymes, to reduce DNA fragmentation. The converted DNA was then sequenced using a PromethION nanopore sequencer. We have also developed a new analysis pipeline that accounts for base conversions, which are not present in conventional nanopore sequencing, as well as errors produced by nanopore sequencing.

Intensive single-cell analysis reveals immune-cell diversity among healthy individuals.

Immune responses are different between individuals and personal health histories and unique environmental conditions should collectively determine the present state of immune cells. However, the molecular systems underlying such heterogeneity remain elusive. Here, we conducted a systematic time-lapse single-cell analysis, using 171 single-cell libraries and 30 mass cytometry datasets intensively for seven healthy individuals. We found substantial diversity in immune-cell profiles between different individuals. These patterns showed daily fluctuations even within the same individual. Similar diversities were also observed for the T-cell and B-cell receptor repertoires. Detailed immune-cell profiles at healthy statuses should give essential background information to understand their immune responses, when the individual is exposed to various environmental conditions. To demonstrate this idea, we conducted the similar analysis for the same individuals on the vaccination of influenza and SARS-CoV-2. In fact, we detected distinct responses to vaccines between individuals, although key responses are common. Single-cell immune-cell profile data should make fundamental data resource to understand variable immune responses, which are unique to each individual.

Extreme value theory as a framework for understanding mutation frequency distribution in cancer genomes.

Currently, the population dynamics of preclonal cancer cells before clonal expansion of tumors has not been sufficiently addressed thus far. By focusing on preclonal cancer cell population as a Darwinian evolutionary system, we formulated and analyzed the observed mutation frequency among tumors (MFaT) as a proxy for the hypothesized sequence read frequency and beneficial fitness effect of a cancer driver mutation. Analogous to intestinal crypts, we assumed that sample donor patients are separate culture tanks where proliferating cells follow certain population dynamics described by extreme value theory (EVT). To validate this, we analyzed three large-scale cancer genome datasets, each harboring > 10000 tumor samples and in total involving > 177898 observed mutation sites. We clarified the necessary premises for the application of EVT in the strong selection and weak mutation (SSWM) regime in relation to cancer genome sequences at scale. We also confirmed that the stochastic distribution of MFaT is likely of the Fréchet type, which challenges the well-known Gumbel hypothesis of beneficial fitness effects. Based on statistical data analysis, we demonstrated the potential of EVT as a population genetics framework to understand and explain the stochastic behavior of driver-mutation frequency in cancer genomes as well as its applicability in real cancer genome sequence data.

Identification of deleterious recessive haplotypes and candidate deleterious recessive mutations in Japanese Black cattle.

Intensive use of a few elite sires has increased the risk of the manifestation of deleterious recessive traits in cattle. Substantial genotyping data gathered using single-nucleotide polymorphism (SNP) arrays have identified the haplotypes with homozygous deficiency, which may compromise survival. We developed Japanese Black cattle haplotypes (JBHs) using SNP array data (4843 individuals) and identified deleterious recessive haplotypes using exome sequencing of 517 sires. We identified seven JBHs with homozygous deficiency. JBH_10 and JBH_17 were associated with the resuming of estrus after artificial insemination, indicating that these haplotypes carried deleterious mutations affecting embryonic survival. The exome data of 517 Japanese Black sires revealed that AC_000165.1:g.85341291C>G of IARS in JBH_8_2, AC_000174.1:g.74743512G>T of CDC45 in JBH_17, and a copy variation region (CNVR_27) of CLDN16 in JBH_1_1 and JBH_1_2 were the candidate mutations. A novel variant AC_000174.1:g.74743512G>T of CDC45 in JBH_17 was located in a splicing donor site at a distance of 5 bp, affecting pre-mRNA splicing. Mating between heterozygotes of JBH_17 indicated that homozygotes carrying the risk allele died around the blastocyst stage. Analysis of frequency of the CDC45 risk allele revealed that its carriers were widespread throughout the tested Japanese Black cattle population. Our approach can effectively manage the inheritance of recessive risk alleles in a breeding population.

Complete genome sequence and bioinformatics analysis of nine Egyptian females with clinical information from different geographic regions in Egypt.

Egyptians are at a crossroad between Africa and Eurasia, providing useful genomic resources for analyzing both genetic and environmental factors for future personalized medicine. Two personal Egyptian whole genomes have been published previously by us and here nine female whole genome sequences with clinical information have been added to expand the genomic resource of Egyptian personal genomes. Here we report the analysis of whole genomes of nine Egyptian females from different regions using Illumina short-read sequencers. At 30x sequencing coverage, we identified 12 SNPs that were shared in most of the subjects associated with obesity which are concordant with their clinical diagnosis. Also, we found mtDNA mutation A4282G is common in all the samples and this is associated with chronic progressive external ophthalmoplegia (CPEO). Haplogroup and Admixture analyses revealed that most Egyptian samples are close to the other north Mediterranean, Middle Eastern, and European, respectively, possibly reflecting the into-Africa influx of human migration. In conclusion, we present whole-genome sequences of nine Egyptian females with personal clinical information that cover the diverse regions of Egypt. Although limited in sample size, the whole genomes data provides possible geno-phenotype candidate markers that are relevant to the region's diseases.

MCM10 compensates for Myc-induced DNA replication stress in breast cancer stem-like cells.

Cancer stem-like cells (CSCs) induce drug resistance and recurrence of tumors when they experience DNA replication stress. However, the mechanisms underlying DNA replication stress in CSCs and its compensation remain unclear. Here, we demonstrate that upregulated c-Myc expression induces stronger DNA replication stress in patient-derived breast CSCs than in differentiated cancer cells. Our results suggest critical roles for mini-chromosome maintenance protein 10 (MCM10), a firing (activating) factor of DNA replication origins, to compensate for DNA replication stress in CSCs. MCM10 expression is upregulated in CSCs and is maintained by c-Myc. c-Myc-dependent collisions between RNA transcription and DNA replication machinery may occur in nuclei, thereby causing DNA replication stress. MCM10 may activate dormant replication origins close to these collisions to ensure the progression of replication. Moreover, patient-derived breast CSCs were found to be dependent on MCM10 for their maintenance, even after enrichment for CSCs that were resistant to paclitaxel, the standard chemotherapeutic agent. Further, MCM10 depletion decreased the growth of cancer cells, but not of normal cells. Therefore, MCM10 may robustly compensate for DNA replication stress and facilitate genome duplication in cancer cells in the S-phase, which is more pronounced in CSCs. Overall, we provide a preclinical rationale to target the c-Myc-MCM10 axis for preventing drug resistance and recurrence of tumors.

Corrigendum: Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

[This corrects the article DOI: 10.3389/fmicb.2020.01628.].

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera.

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister-taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault's rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.

Spindle assembly checkpoint competence in aneuploid canine malignant melanoma cell lines.

The spindle assembly checkpoint (SAC) is a surveillance mechanism that prevents unequal segregation of chromosomes during mitosis. Abnormalities in the SAC are associated with chromosome instability and resultant aneuploidy. This study was performed to evaluate the SAC competence in canine malignant melanoma (CMM) using four aneuploid cell lines (CMeC1, CMeC2, KMeC, and LMeC). After treatment with nocodazole, a microtubule disrupting agent, CMeC1, KMeC, and LMeC cells were arrested in M phase, whereas CMeC2 cells were not arrested, and progressed into the next cell cycle phase without cytokinesis. Chromosome spread analysis revealed a significantly increased rate of premature sister chromatid separation in CMeC2 cells. Expression of the phosphorylated form of the SAC regulator, monopolar spindle 1 (Mps1), was lower in CMeC2 cells than in the other CMM cell lines. These results indicate that the SAC is defective in CMeC2 cells, which may partially explain aneuploidy in CMM. Thus, CMeC2 cells may be useful for further studies of the SAC mechanism in CMM and in determining the relationship between SAC incompetence and aneuploidy.

The Origin and Composition of Korean Ethnicity Analyzed by Ancient and Present-Day Genome Sequences.

Koreans are thought to be an ethnic group of admixed northern and southern subgroups. However, the exact genetic origins of these two remain unclear. In addition, the past admixture is presumed to have taken place on the Korean peninsula, but there is no genomic scale analysis exploring the origin, composition, admixture, or the past migration of Koreans. Here, 88 Korean genomes compared with 91 other present-day populations showed two major genetic components of East Siberia and Southeast Asia. Additional paleogenomic analysis with 115 ancient genomes from Pleistocene hunter-gatherers to Iron Age farmers showed a gradual admixture of Tianyuan (40 ka) and Devil's gate (8 ka) ancestries throughout East Asia and East Siberia up until the Neolithic era. Afterward, the current genetic foundation of Koreans may have been established through a rapid admixture with ancient Southern Chinese populations associated with Iron Age Cambodians. We speculate that this admixing trend initially occurred mostly outside the Korean peninsula followed by continuous spread and localization in Korea, corresponding to the general admixture trend of East Asia. Over 70% of extant Korean genetic diversity is explained to be derived from such a recent population expansion and admixture from the South.

Characterization of cancer omics and drug perturbations in panels of lung cancer cells.

To better understand the disruptions of transcriptional regulations and gene expression in lung cancers, we constructed a multi-omics catalogue of the responses of lung cancer cells to a series of chemical compounds. We generated and analyzed 3,240 RNA-seq and 3,393 ATAC-seq libraries obtained from 23 cell lines treated with 95 well-annotated compounds. To demonstrate the power of the created multi-omics resource, we attempted to identify drugs that could induce the designated changes alone or in combination. The basal multi-omics information was first integrated into co-expression modules. Among these modules, we identified a stress response module that may be a promising drug intervention target, as new combinations of compounds that could be used to regulate this module and the consequent phenotypic appearance of cancer cells have been identified. We believe that the multi-omics profiles generated in this study and the strategy used to stratify them will lead to more rational and efficient development of anticancer drugs.

DNA aneuploidy and centrosome amplification in canine tumor cell lines.

DNA aneuploidy, the altered DNA content of a cell, is a common feature of canine tumors. However, it is unclear whether aneuploid DNA in canine tumor cells show centrosome amplification (CA), which contributes to numerical and structural chromosome aberrations that result in DNA aneuploidy. Here, we evaluated whether DNA aneuploidy and CA occur concurrently in canine tumor cell lines. Centrosome numbers were evaluated in 18 canine tumor cell lines by immunocytochemistry with anti-γ-tubulin antibody, and DNA content was evaluated by flow cytometry using propidium iodide. A total of 15 cell lines showed DNA aneuploidy, and CA was observed in 5 of these 15 cell lines. Together, our results suggest that DNA aneuploidy in canine tumor cells might be explained at least in part by CA. In addition, cell lines with CA may be useful tools to examine the detailed relationship between CA and DNA aneuploidy and the molecular mechanism of CA in canine tumor.

Monoallelic, antisense and total RNA transcription in an in vitro neural differentiation system based on F1 hybrid mice.

We developed an in vitro system to differentiate embryonic stem cells (ESCs) derived from reciprocally crossed F1 hybrid mice into neurons, and used it to investigate poly(A)+ and total RNA transcription at different stages of cell differentiation. By comparing expression profiles of transcripts assembled from 20 RNA sequencing datasets [2 alleles×(2 cell lines×4 time-points+2 mouse brains)], the relative influence of strain, cell and parent specificities to overall expression could be assessed. Divergent expression profiles of ESCs converged tightly at neural progenitor stage. Patterns of temporal variation of monoallelically expressed transcripts and antisense transcripts were quantified. Comparison of sense and antisense transcript pairs within the poly(A)+ sample, within the total RNA sample, and across poly(A)+ and total RNA samples revealed distinct rates of pairs showing anti-correlated expression variation. Unique patterns of sharing of poly(A)+ and poly(A)- transcription were identified in distinct RNA species. Regulation and functionality of monoallelic expression, antisense transcripts and poly(A)- transcription remain elusive. We demonstrated the effectiveness of our approach to capture these transcriptional activities, and provided new resources to elucidate the mammalian developmental transcriptome.

Quantifying local malignant adaptation in tissue-specific evolutionary trajectories by harnessing cancer's repeatability at the genetic level.

Cancer is a potentially lethal disease, in which patients with nearly identical genetic backgrounds can develop a similar pathology through distinct combinations of genetic alterations. We aimed to reconstruct the evolutionary process underlying tumour initiation, using the combination of convergence and discrepancies observed across 2,742 cancer genomes from nine tumour types. We developed a framework using the repeatability of cancer development to score the local malignant adaptation (LMA) of genetic clones, as their potential to malignantly progress and invade their environment of origin. Using this framework, we found that premalignant skin and colorectal lesions appeared specifically adapted to their local environment, yet insufficiently for full cancerous transformation. We found that metastatic clones were more adapted to the site of origin than to the invaded tissue, suggesting that genetics may be more important for local progression than for the invasion of distant organs. In addition, we used network analyses to investigate evolutionary properties at the system-level, highlighting that different dynamics of malignant progression can be modelled by such a framework in tumour-type-specific fashion. We find that occurrence-based methods can be used to specifically recapitulate the process of cancer initiation and progression, as well as to evaluate the adaptation of genetic clones to given environments. The repeatability observed in the evolution of most tumour types could therefore be harnessed to better predict the trajectories likely to be taken by tumours and preneoplastic lesions in the future.

Evaluation and application of RNA-Seq by MinION.

The current RNA-Seq method analyses fragments of mRNAs, from which it is occasionally difficult to reconstruct the entire transcript structure. Here, we performed and evaluated the recent procedure for full-length cDNA sequencing using the Nanopore sequencer MinION. We applied MinION RNA-Seq for various applications, which would not always be easy using the usual RNA-Seq by Illumina. First, we examined and found that even though the sequencing accuracy was still limited to 92.3%, practically useful RNA-Seq analysis is possible. Particularly, taking advantage of the long-read nature of MinION, we demonstrate the identification of splicing patterns and their combinations as a form of full-length cDNAs without losing precise information concerning their expression levels. Transcripts of fusion genes in cancer cells can also be identified and characterized. Furthermore, the full-length cDNA information can be used for phasing of the SNPs detected by WES on the transcripts, providing essential information to identify allele-specific transcriptional events. We constructed a catalogue of full-length cDNAs in seven major organs for two particular individuals and identified allele-specific transcription and splicing. Finally, we demonstrate that single-cell sequencing is also possible. RNA-Seq on the MinION platform should provide a novel approach that is complementary to the current RNA-Seq.