Multiancestry genomic and transcriptomic analysis of gastric cancer.

Gastric cancer is among the most common malignancies worldwide, characterized by geographical, epidemiological and histological heterogeneity. Here, we report an extensive, multiancestral landscape of driver events in gastric cancer, involving 1,335 cases. Seventy-seven significantly mutated genes (SMGs) were identified, including ARHGAP5 and TRIM49C. We also identified subtype-specific drivers, including PIGR and SOX9, which were enriched in the diffuse subtype of the disease. SMGs also varied according to Epstein-Barr virus infection status and ancestry. Non-protein-truncating CDH1 mutations, which are characterized by in-frame splicing alterations, targeted localized extracellular domains and uniquely occurred in sporadic diffuse-type cases. In patients with gastric cancer with East Asian ancestry, our data suggested a link between alcohol consumption or metabolism and the development of RHOA mutations. Moreover, mutations with potential roles in immune evasion were identified. Overall, these data provide comprehensive insights into the molecular landscape of gastric cancer across various subtypes and ancestries.

Single-cell DNA and RNA sequencing reveals the dynamics of intra-tumor heterogeneity in a colorectal cancer model.

BACKGROUND: Intra-tumor heterogeneity (ITH) encompasses cellular differences in tumors and is related to clinical outcomes such as drug resistance. However, little is known about the dynamics of ITH, owing to the lack of time-series analysis at the single-cell level. Mouse models that recapitulate cancer development are useful for controlled serial time sampling. RESULTS: We performed single-cell exome and transcriptome sequencing of 200 cells to investigate how ITH is generated in a mouse colorectal cancer model. In the model, a single normal intestinal cell is grown into organoids that mimic the intestinal crypt structure. Upon RNAi-mediated downregulation of a tumor suppressor gene APC, the transduced organoids were serially transplanted into mice to allow exposure to in vivo microenvironments, which play relevant roles in cancer development. The ITH of the transcriptome increased after the transplantation, while that of the exome decreased. Mutations generated during organoid culture did not greatly change at the bulk-cell level upon the transplantation. The RNA ITH increase was due to the emergence of new transcriptional subpopulations. In contrast to the initial cells expressing mesenchymal-marker genes, new subpopulations repressed these genes after the transplantation. Analyses of colorectal cancer data from The Cancer Genome Atlas revealed a high proportion of metastatic cases in human subjects with expression patterns similar to the new cell subpopulations in mouse. These results suggest that the birth of transcriptional subpopulations may be a key for adaptation to drastic micro-environmental changes when cancer cells have sufficient genetic alterations at later tumor stages. CONCLUSIONS: This study revealed an evolutionary dynamics of single-cell RNA and DNA heterogeneity in tumor progression, giving insights into the mesenchymal-epithelial transformation of tumor cells at metastasis in colorectal cancer.

Comprehensive analysis of DNA adducts (DNA adductome analysis) in the liver of rats treated with 1,4-dioxane.

1,4-Dioxane is a genotoxic carcinogen, and its mutagenic properties were recently observed in the liver of guanine phosphoribosyl transferase (gpt) delta transgenic rats. However, the mechanisms of its genotoxicity remain unclear. We analyzed DNA adduct formation in rat livers following 1,4-dioxane treatment. After administering 1,4-dioxane in drinking water at doses of 0, 20, 200, and 5,000 ppm, liver adduct formation was analyzed by DNA adductome analysis. Adducts in treated rat livers were dose-dependently increased compared with those in the control group. Principal component analysis-discriminant analysis (PCA-DA) clearly revealed two clusters of DNA adducts, associated with 0 ppm and low-dose (20 ppm) 1,4-dioxane-treatment versus middle- and high-dose (200, 5,000 ppm)-treated rats. After confirming the intensity of each adduct, three adducts were screened as characteristic of 1,4-dioxane treatment. Two of the three candidates contained thymine or cytidine/uracil moieties. Another candidate was identified as 8-oxo-dG based on mass fragmentation together with high-resolution accurate-mass (HRAM) mass spectrometry data. Oxidative stress responses may partly explain the mechanisms of increased mutations in the liver of gpt delta rats following 1,4-dioxane treatment.

Removal of redundant contigs from de novo RNA-Seq assemblies via homology search improves accurate detection of differentially expressed genes.

BACKGROUND: For plant species with unsequenced genomes, cDNA contigs created by de novo assembly of RNA-Seq reads are used as reference sequences for comparative analysis of RNA-Seq datasets and the detection of differentially expressed genes (DEGs). Redundancies in such contigs are evident in previous RNA-Seq studies, and such redundancies can lead to difficulties in subsequent analysis. Nevertheless, the effects of removing redundancy from contig assemblies on comparative RNA-Seq analysis have not been evaluated. RESULTS: Here we describe a method for removing redundancy from raw contigs that were primarily created by de novo assembly of Arabidopsis thaliana RNA-Seq reads. Specifically, the contigs with the highest bit scores were selected from raw contigs by a homology search against the gene dataset in the TAIR10 database. The two existing methods for removal of redundancy based on contig length or clustering analysis used to eliminate redundancies from raw contigs. Contig number was reduced most effectively with the method based on homology search. In a comparative analysis of RNA-Seq datasets, DEGs detected in contigs that underwent redundancy removal via the homology search method showed the highest identity to the DEGs detected when the TAIR10 gene dataset was used as an exact reference. Redundancy in raw contigs could also be removed by a homology search against integrated protein datasets from several plant species other than A. thaliana. DEGs detected using contigs that underwent such redundancy-removed also showed high homology to DEGs detected using the TAIR10 gene dataset. CONCLUSION: Here we describe a method for removing redundant contigs within raw contigs; this method involves a homology search against a gene or protein database. In principal, this method can be used with unsequenced plant genomes that lack a well-developed gene database. Redundant contigs were not removed adequately via either of two existing methods, but our method allowed for removal of all redundant contigs. To our knowledge, this is the first reported improvement in accurate detection of DEGs via comparative RNA-Seq analysis that involved preparation of a non-redundant reference sequence. This method could be used to rapidly and cost-effectively detect useful genes in unsequenced plants.

Accumulation of long-lived mRNAs associated with germination in embryos during seed development of rice.

Mature dry seeds contain translatable mRNAs called long-lived mRNAs. Early studies have shown that protein synthesis during the initial phase of seed germination occurs from long-lived mRNAs, without de novo transcription. However, the gene expression systems that generate long-lived mRNAs in seeds are not well understood. To examine the accumulation of long-lived mRNAs in developing rice embryos, germination tests using the transcriptional inhibitor actinomycin D (Act D) were performed with the Japonica rice cultivar Nipponbare. Although over 70% of embryos at 10 days after flowering (DAF) germinated in the absence of the inhibitor, germination was remarkably impaired in embryos treated with Act D. In contrast, more than 70% of embryos at 20, 25, 30 and 40 DAF germinated in the presence of Act D. The same results were obtained when another cultivar, Koshihikari, was used, indicating that the long-lived mRNAs required for germination predominantly accumulate in embryos between 10 and 20 DAF during seed development. RNA-Seq identified 529 long-lived mRNA candidates, encoding proteins such as ABA, calcium ion and phospholipid signalling-related proteins, and HSP DNA J, increased from 10 to 20 DAF and were highly abundant in 40 DAF embryos of Nipponbare and Koshihikari. We also revealed that these long-lived mRNA candidates are clearly up-regulated in 10 DAF germinating embryos after imbibition, suggesting that the accumulation of these mRNAs in embryos is indispensable for the induction of germination. The findings presented here may facilitate in overcoming irregular seed germination or producing more vigorous seedlings.

Length of the dark period affects flower opening and the expression of circadian-clock associated genes as well as xyloglucan endotransglucosylase/hydrolase genes in petals of morning glory (Ipomoea nil).

KEY MESSAGE: We isolated differentially expressed and dark-responsive genes during flower development and opening in petals of morning glory. Flower opening usually depends on petal expansion and is regulated by both genetic and environmental factors. Flower opening in morning glory (Ipomoea nil) is controlled by the dark/light regime just prior to opening. Opening was normal after 8- or 12-h dark periods but progressed very slowly after a 4-h dark period or in continuous light. Four genes (InXTH1-InXTH4) encoding xyloglucan endotransglucosylase/hydrolases (XTHs) and three genes (InEXPA1-InEXPA3) encoding alpha-expansins (EXPAs) were isolated. The expression patterns of InXTH2, InXTH3, and InXTH4 in petals were closely correlated with the rate of flower opening controlled by the length of the dark period prior to opening, but those of the EXPA genes were not. The expression pattern of InXTH1 gene was closely correlated with petal elongation. Suppression subtractive hybridization was used to isolate dark-responsive genes accompanying flower opening. The expressions of ten isolated genes were associated with the length of the dark period prior to flower opening. One gene was highly homologous to Arabidopsis pseudo-response regulator7, which is associated with the circadian clock and phytochrome signaling; another to Arabidopsis REVEILLE1, which affects the output of the circadian clock. Other genes were related to light responses, plant hormone effects and signal transduction. The possible roles of these genes in regulation of flower opening are discussed.

Zigzag motions of the myosin-coated beads actively sliding along actin filaments suspended between immobilized beads.

The motions of myosin filaments actively sliding along suspended actin filaments were studied. By manipulating a double-beam laser tweezers, single actin filaments were suspended between immobilized microbeads. When another beads coated with myosin filaments were dragged to suspended actin filaments, the beads instantly and unidirectionally slid along the actin filaments. The video image analysis showed that the beads slid at a velocity of ca. 3-5 microm/s accompanied with zigzag motions. When beads were densely coated with myosin filaments, the sliding motions became straight and smooth. The obtained results indicate that (1) during the sliding motions, the interaction between myosin heads and actin filaments is weak and susceptible to random thermal agitations, (2) the effects of thermal agitations to the sliding motions of myofilaments are readily suppressed by mechanical constraints imposed to the filaments, and (3) the active sliding force is produced almost in parallel to the filaments axis.