Whole-genome sequencing reveals the molecular implications of the stepwise progression of lung adenocarcinoma.

The mechanism underlying the development of tumors, particularly at early stages, still remains mostly elusive. Here, we report whole-genome long and short read sequencing analysis of 76 lung cancers, focusing on very early-stage lung adenocarcinomas such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma. The obtained data is further integrated with bulk and spatial transcriptomic data and epigenomic data. These analyses reveal key events in lung carcinogenesis. Minimal somatic mutations in pivotal driver mutations and essential proliferative factors are the only detectable somatic mutations in the very early-stage of AIS. These initial events are followed by copy number changes and global DNA hypomethylation. Particularly, drastic changes are initiated at the later AIS stage, i.e., in Noguchi type B tumors, wherein cancer cells are exposed to the surrounding microenvironment. This study sheds light on the pathogenesis of lung adenocarcinoma from integrated pathological and molecular viewpoints.

Bioinformatic Identification of Potential RNA Alterations on the Atrial Fibrillation Remodeling from Human Pulmonary Veins.

Atrial fibrillation (AF) is the most frequent persistent arrhythmia. Many genes have been reported as a genetic background for AF. However, most transcriptome analyses of AF are limited to the atrial samples and have not been evaluated by multiple cardiac regions. In this study, we analyzed the expression levels of protein-coding and long noncoding RNAs (lncRNAs) in six cardiac regions by RNA-seq. Samples were donated from six subjects with or without persistent AF for left atria, left atrial appendages, right atria, sinoatrial nodes, left ventricles, right ventricles, and pulmonary veins (PVs), and additional four right atrial appendages samples were collected from patients undergoing mitral valve replacement. In total, 23 AF samples were compared to 23 non-AF samples. Surprisingly, the most influenced heart region in gene expression by AF was the PV, not the atria. The ion channel-related gene set was significantly enriched upon analysis of these significant genes. In addition, some significant genes are cancer-related lncRNAs in PV in AF. A co-expression network analysis could detect the functional gene clusters. In particular, the cancer-related lncRNA, such as SAMMSON and FOXCUT, belong to the gene network with the cancer-related transcription factor FOXC1. Thus, they may also play an aggravating role in the pathogenesis of AF, similar to carcinogenesis. In the least, this study suggests that (1) RNA alteration is most intense in PVs and (2) post-transcriptional gene regulation by lncRNA may contribute to the progression of AF. Through the screening analysis across the six cardiac regions, the possibility that the PV region can play a role other than paroxysmal triggering in the pathogenesis of AF was demonstrated for the first time. Future research with an increase in the number of PV samples will lead to a novel understanding of the pathophysiology of AF.

Comparative Study of Transcriptome in the Hearts Isolated from Mice, Rats, and Humans.

The heart is a significant organ in mammalian life, and the heartbeat mechanism has been an essential focus of science. However, few studies have focused on species differences. Accordingly, challenges remain in studying genes that have universal functions across species and genes that determine species differences. Here, we analyzed transcriptome data in mouse, rat, and human atria, ventricles, and sinoatrial nodes (SA) obtained from different platforms and compared them by calculating specificity measure (SPM) values in consideration of species differences. Among the three heart regions, the species differences in SA were the greatest, and we searched for genes that determined the essential characteristics of SA, which was SHOX2 in our criteria. The SPM value of SHOX2 was prominently high across species. Similarly, by calculating SPM values, we identified 3 atrial-specific, 11 ventricular-specific, and 17 SA-specific markers. Ontology analysis identified 70 cardiac region- and species-specific ontologies. These results suggest that reanalyzing existing data by calculating SPM values may identify novel tissue-specific genes and species-dependent gene expression. This study identified the importance of SHOX2 as an SA-specific transcription factor, a novel cardiac regional marker, and species-dependent ontologies.

Phasing analysis of lung cancer genomes using a long read sequencer.

Chromosomal backgrounds of cancerous mutations still remain elusive. Here, we conduct the phasing analysis of non-small cell lung cancer specimens of 20 Japanese patients. By the combinatory use of short and long read sequencing data, we obtain long phased blocks of 834 kb in N50 length with >99% concordance rate. By analyzing the obtained phasing information, we reveal that several cancer genomes harbor regions in which mutations are unevenly distributed to either of two haplotypes. Large-scale chromosomal rearrangement events, which resemble chromothripsis events but have smaller scales, occur on only one chromosome, and these events account for the observed biased distributions. Interestingly, the events are characteristic of EGFR mutation-positive lung adenocarcinomas. Further integration of long read epigenomic and transcriptomic data reveal that haploid chromosomes are not always at equivalent transcriptomic/epigenomic conditions. Distinct chromosomal backgrounds are responsible for later cancerous aberrations in a haplotype-specific manner.

Transcript Identification Through Long-Read Sequencing.

RNA-seq using long-read sequencing, such as nanopore and SMRT (Single Molecule, Real-Time) sequencing, enabled the identification of the full-length structure of RNA molecules. Several tools for long-read RNA-seq were developed recently. In this section, we introduce an analytical pipeline of long-read RNA-seq for isoform identification and the estimation of expression levels using minimap2, TranscriptClean, and TALON. We applied this pipeline to the public direct RNA-seq data of the HAP1 and HEK293 cell lines to identify transcript isoforms which can be detected only using long-read RNA-seq data.

Aberrant splicing isoforms detected by full-length transcriptome sequencing as transcripts of potential neoantigens in non-small cell lung cancer.

BACKGROUND: Long-read sequencing of full-length cDNAs enables the detection of structures of aberrant splicing isoforms in cancer cells. These isoforms are occasionally translated, presented by HLA molecules, and recognized as neoantigens. This study used a long-read sequencer (MinION) to construct a comprehensive catalog of aberrant splicing isoforms in non-small-cell lung cancers, by which novel isoforms and potential neoantigens are identified. RESULTS: Full-length cDNA sequencing is performed using 22 cell lines, and a total of 2021 novel splicing isoforms are identified. The protein expression of some of these isoforms is then validated by proteome analysis. Ablations of a nonsense-mediated mRNA decay (NMD) factor, UPF1, and a splicing factor, SF3B1, are found to increase the proportion of aberrant transcripts. NetMHC evaluation of the binding affinities to each type of HLA molecule reveals that some of the isoforms potentially generate neoantigen candidates. We also identify aberrant splicing isoforms in seven non-small-cell lung cancer specimens. An enzyme-linked immune absorbent spot assay indicates that approximately half the peptide candidates have the potential to activate T cell responses through their interaction with HLA molecules. Finally, we estimate the number of isoforms in The Cancer Genome Atlas (TCGA) datasets by referring to the constructed catalog and found that disruption of NMD factors is significantly correlated with the number of splicing isoforms found in the TCGA-Lung Adenocarcinoma data collection. CONCLUSIONS: Our results indicate that long-read sequencing of full-length cDNAs is essential for the precise identification of aberrant transcript structures in cancer cells.

Loss of the small GTPase Arl8b results in abnormal development of the roof plate in mouse embryos.

Lysosomes are acidic organelles responsible for degrading both exogenous and endogenous materials. The small GTPase Arl8 localizes primarily to lysosomes and is involved in lysosomal function. In the present study, using Arl8b gene-trapped mutant (Arl8b-/- ) mice, we show that Arl8b is required for the development of dorsal structures of the neural tube, including the thalamus and hippocampus. In embryonic day (E) 10.5 Arl8b-/- embryos, Sox1 (a neuroepithelium marker) was ectopically expressed in the roof plate, whereas the expression of Gdf7 and Msx1 (roof plate markers) was reduced in the dorsal midline of the midbrain. Ectopic expression of Sox1 in Arl8b-/- embryos was detected also at E9.0 in the neural fold, which gives rise to the roof plate. In addition, the levels of Bmp receptor IA and phosphorylated Smad 1/5/8 (downstream of BMP signaling) were increased in the neural fold of E9.0 Arl8b-/- embryos. These results suggest that Arl8b is involved in the development of the neural fold and the subsequently formed roof plate, possibly via control of BMP signaling.

Arl8b is required for lysosomal degradation of maternal proteins in the visceral yolk sac endoderm of mouse embryos.

The small GTPase Arl8b localizes primarily to lysosomes and is involved in lysosomal motility and fusion. Here, we show that Arl8b is required for lysosomal degradation of maternal proteins in the visceral yolk sac endoderm (VYSE), an apical cell layer of the visceral yolk sac, of mouse embryos. The VYSE actively takes up maternal materials from uterine fluid and degrades them in lysosomes to provide breakdown products to the embryo. Arl8b gene-trap mice (Arl8b-/- ) displayed decreased early embryo body size. The Arl8b-/-  VYSE exhibited defective endocytic trafficking to the lysosome and accumulation of maternal proteins such as albumin and immunoglobulin G in late endocytic organelles. Furthermore, Transthyretin-Cre;Arl8bflox/flox mice in which Arl8b was ablated specifically in the VYSE also showed decreased embryo body size, defects in trafficking to the lysosome and reduction of the free amino acid level in the embryos. Taken together, these results suggest that Arl8b mediates lysosomal degradation of maternal proteins in the VYSE, thereby contributing to mouse embryonic development.

Pochonicine, a polyhydroxylated pyrrolizidine alkaloid from fungus Pochonia suchlasporia var. suchlasporia TAMA 87 as a potent beta-N-acetylglucosaminidase inhibitor.

A new polyhydroxylated pyrrolizidine alkaloid designated as pochonicine (1) was isolated from a solid fermentation culture of the fungal strain Pochonia suchlasporia var. suchlasporia TAMA 87. The structure of 1 was determined using NMR and MS techniques as (1R*, 3S*, 5S*, 6S*, 7R*, 7a S*)-5-acetamidomethyl-3-hydroxymethyl-1,6,7-trihydroxypyrrolizidine. Pochonicine (1) showed potent inhibition against beta-N-acetylglucosaminidases (GlcNAcases) of various organisms including insects, fungi, mammals, and a plant but no inhibition against beta-glucosidase of almond, alpha-glucosidase of yeast, or chitinase of Bacillus sp. The GlcNAcase inhibitory activity of pochonicine (1) was comparable to nagstatin, a potent GlcNAcase inhibitor of natural origin.