Recombination of repeat elements generates somatic complexity in human genomes.

Non-allelic recombination between homologous repetitive elements contributes to evolution and human genetic disorders. Here, we combine short- and long-DNA read sequencing of repeat elements with a new bioinformatics pipeline to show that somatic recombination of Alu and L1 elements is widespread in the human genome. Our analysis uncovers tissue-specific non-allelic homologous recombination hallmarks; moreover, we find that centromeres and cancer-associated genes are enriched for retroelements that may act as recombination hotspots. We compare recombination profiles in human-induced pluripotent stem cells and differentiated neurons and find that the neuron-specific recombination of repeat elements accompanies chromatin changes during cell-fate determination. Finally, we report that somatic recombination profiles are altered in Parkinson's and Alzheimer's disease, suggesting a link between retroelement recombination and genomic instability in neurodegeneration. This work highlights a significant contribution of the somatic recombination of repeat elements to genomic diversity in health and disease.

Embryonic LTR retrotransposons supply promoter modules to somatic tissues.

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

Clinical Utility of Inflammatory and Nutritious Index as Therapeutic Prediction of Nivolumab plus Ipilimumab in Advanced Non-Small Cell Lung Cancer.

INTRODUCTION: Biomarkers for predicting the outcome of ipilimumab plus nivolumab (Nivo-Ipi) treatment in cancer patients have not been identified. Herein, we investigated the prognostic significance of inflammatory and nutritional markers in patients with advanced non-small cell lung cancer (NSCLC) receiving Nivo-Ipi. METHODS: Our study retrospectively analyzed 101 patients with advanced NSCLC who received Nivo-Ipi at a single institution. Inflammatory and nutritional indices were correlated with patient outcomes and included the neutrophil to lymphocyte ratio (NLR), platelet to lymphocyte ratio (PLR), systemic immune-inflammation index (SII), prognostic nutritional index (PNI), advanced lung cancer inflammation index (ALI), and Glasgow prognostic score (GPS). RESULTS: The NLR significantly correlated with the PLR, SII, PNI, ALI, and GPS. Regarding therapeutic efficacy, the NLR, SII, and PNI predicted a partial response, and all indices predicted progressive disease. In subgroup analyses, the SII, PNI, and ALI predicted the outcome of patients with adenocarcinoma, whereas only the PNI predicted the outcome of patients with non-adenocarcinoma. The PNI and SII were the most useful indices in patients with a programmed death ligand-1 expression level of <1% and ≥1%, respectively. CONCLUSION: The NLR, PLR, SII, PNI, ALI, and GPS were significantly associated with the outcome of Nivo-Ipi treatment in patients with NSCLC. The PNI was the most suitable marker regardless of histological type. The SII and PNI were the most promising markers for patients with and without PD-L1 expression, respectively.

The Glasgow Prognostic Score as a Predictor of Survival after Chemoradiotherapy for Limited-Disease Small Cell Lung Cancer.

INTRODUCTION: Established biomarkers for predicting chemoradiotherapy efficacy for limited-disease small cell lung cancer (LD-SCLC) are lacking. The inflammation-based Glasgow Prognostic Score (GPS), comprising serum C-reactive protein (CRP) and albumin levels, can predict survival in advanced cancer. This study investigated whether metabolic and inflammatory markers, including the GPS, can predict the efficacy of chemoradiotherapy in patients with LD-SCLC. METHODS: We retrospectively analyzed 124 patients who underwent chemoradiotherapy for LD-SCLC at two institutions between April 2007 and June 2021, and assessed the prognostic significance of various metabolic and inflammatory markers. The GPS was calculated using the CRP and albumin concentrations, and categorized as follows: 0, CRP <1.0 mg/dL and albumin ≥3.5 mg/dL; 1, elevated CRP or decreased albumin; and 2, CRP ≥1.0 mg/dL and albumin<3.5 mg/dL. Differences in progression-free survival (PFS) and overall survival (OS) were examined using Kaplan-Meier curves and Cox proportional-hazard models. RESULTS: The overall response rate was 95.1% (95% confidence interval [CI]: 89.6-97.9%). The median PFS and OS from chemoradiotherapy initiation were 12.6 (95% CI: 9.9-15.4) and 29.0 (95% CI: 24.8-45.5) months, respectively. The GPS demonstrated independent predictive ability for the effectiveness of chemoradiotherapy, wherein favorable scores (GPS 0-1) were significantly correlated with superior PFS and OS compared to unfavorable scores (GPS 2: PFS: 14.8 vs. 6.7 months, p = 0.0001; OS: 35.4 vs. 11.0 months, p < 0.0001). CONCLUSION: This preliminary examination revealed that the GPS was significantly associated with PFS and OS in patients undergoing chemoradiotherapy for LD-SCLC, indicating its potential utility in assessing the therapeutic outcomes in LD-SCLC.

Clinical significance of antinuclear antibody as prognostic marker for first-line pembrolizumab in advanced non-small cell lung cancer.

BACKGROUND: The relationship between antinuclear antibody (ANA) and the efficacy of programmed death-1 (PD-1) blockade remains controversial. Here, we investigated the prognostic significance of ANA titer in patients with non-small cell lung cancer (NSCLC) receiving pembrolizumab monotherapy as the first-line treatment, compared with that of platinum-based chemotherapy with PD-1 blockade. METHODS: Our clinical data based on the ANA titer (1:80) were retrospectively reviewed for patients with advanced NSCLC, who were treated with first-line pembrolizumab monotherapy and platinum-based chemotherapy with PD-1 blockade. Immunohistochemical staining for tumor-infiltrating lymphocytes such as CD4, CD8 and Foxp3 was performed. RESULTS: Among 106 patients treated with pembrolizumab, 19 (17.9%) tested high for ANA. Progression-free survival (PFS) and overall survival (OS) were significantly better in patients with high ANA than in those with low ANA, and high ANA was identified as an independent prognostic predictor, particularly in the subgroup with programmed death ligand-1 (PD-L1) ≥ 50%. However, no statistically significant difference in PFS and OS based on the ANA titer was observed in 59 patients treated with combinational chemotherapy and immunotherapy. High numbers of intratumoral Foxp3 and stromal CD8 were significantly associated with low ANA. CONCLUSIONS: Assessment of preexisting ANA titers was useful to prognose PD-1 blockade as a first-line setting, particularly for the PD-L1 ≥ 50% subgroup, but not in the case of combined immunotherapy and chemotherapy.

Correction for Extrinsic Background in Raman Hyperspectral Images.

Raman hyperspectral microscopy is a valuable tool in biological and biomedical imaging. Because Raman scattering is often weak in comparison to other phenomena, prevalent spectral fluctuations and contaminations have brought advancements in analytical and chemometric methods for Raman spectra. These chemometric advances have been key contributors to the applicability of Raman imaging to biological systems. As studies increase in scale, spectral contamination from extrinsic background, intensity from sources such as the optical components that are extrinsic to the sample of interest, has become an emerging issue. Although existing baseline correction schemes often reduce intrinsic background such as autofluorescence originating from the sample of interest, extrinsic background is not explicitly considered, and these methods often fail to reduce its effects. Here, we show that extrinsic background can significantly affect a classification model using Raman images, yielding misleadingly high accuracies in the distinction of benign and malignant samples of follicular thyroid cell lines. To mitigate its effects, we develop extrinsic background correction (EBC) and demonstrate its use in combination with existing methods on Raman hyperspectral images. EBC isolates regions containing the smallest amounts of sample materials that retain extrinsic contributions that are specific to the device or environment. We perform classification both with and without the use of EBC, and we find that EBC retains biological characteristics in the spectra while significantly reducing extrinsic background. As the methodology used in EBC is not specific to Raman spectra, correction of extrinsic effects in other types of hyperspectral and grayscale images is also possible.

Lipid droplet accumulation and adipophilin expression in follicular thyroid carcinoma.

Follicular neoplasms of the thyroid include follicular thyroid carcinoma (FTC) and follicular thyroid adenoma (FTA). However, the differences in cytological findings between FTC and FTA remain undetermined. Here, we aimed to evaluate the accumulation of lipid droplets (LDs) and the expression of adipophilin (perilipin 2/ADRP/ADFP), a known LD marker, in cultured FTC cells. We also immunohistochemically compared adipophilin expression in the FTC and FTA of resected human thyroid tissues. Cultured FTC (FTC-133 and RO82W-1) possessed increased populations of LDs compared to thyroid follicular epithelial (Nthy-ori 3-1) cells. In vitro treatment with phosphatidylinositol-3-kinase (PI3K)/Akt/mammalian target of rapamycin (mTOR) signaling inhibitors (LY294002, MK2206, and rapamycin) in FTC-133 cells downregulated the PI3K/Akt/mTOR/sterol regulatory element-binding protein 1 (SREBP1) signaling pathway, resulting in a significant reduction in LD accumulation. SREBP1 is a master transcription factor that controls lipid metabolism. Fluorescence immunocytochemistry revealed adipophilin expression in the LDs of FTC-133 cells. Immunohistochemical analysis of surgically resected human thyroid tissues revealed significantly increased expression of adipophilin in FTC compared with FTA and adjacent non-tumorous thyroid epithelia. Taken together, LDs and adipophilin were abundant in cultured FTC; the evaluation of adipophilin expression can help distinguish FTC from FTA in surgical specimens.

Functional annotation of human long noncoding RNAs via molecular phenotyping.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

FANTOM enters 20th year: expansion of transcriptomic atlases and functional annotation of non-coding RNAs.

The Functional ANnoTation Of the Mammalian genome (FANTOM) Consortium has continued to provide extensive resources in the pursuit of understanding the transcriptome, and transcriptional regulation, of mammalian genomes for the last 20 years. To share these resources with the research community, the FANTOM web-interfaces and databases are being regularly updated, enhanced and expanded with new data types. In recent years, the FANTOM Consortium's efforts have been mainly focused on creating new non-coding RNA datasets and resources. The existing FANTOM5 human and mouse miRNA atlas was supplemented with rat, dog, and chicken datasets. The sixth (latest) edition of the FANTOM project was launched to assess the function of human long non-coding RNAs (lncRNAs). From its creation until 2020, FANTOM6 has contributed to the research community a large dataset generated from the knock-down of 285 lncRNAs in human dermal fibroblasts; this is followed with extensive expression profiling and cellular phenotyping. Other updates to the FANTOM resource includes the reprocessing of the miRNA and promoter atlases of human, mouse and chicken with the latest reference genome assemblies. To facilitate the use and accessibility of all above resources we further enhanced FANTOM data viewers and web interfaces. The updated FANTOM web resource is publicly available at https://fantom.gsc.riken.jp/.

Dynamics of cardiomyocyte transcriptome and chromatin landscape demarcates key events of heart development.

Organogenesis involves dynamic regulation of gene transcription and complex multipathway interactions. Despite our knowledge of key factors regulating various steps of heart morphogenesis, considerable challenges in understanding its mechanism still exist because little is known about their downstream targets and interactive regulatory network. To better understand transcriptional regulatory mechanism driving heart development and the consequences of its disruption in vivo, we performed time-series analyses of the transcriptome and genome-wide chromatin accessibility in isolated cardiomyocytes (CMs) from wild-type zebrafish embryos at developmental stages corresponding to heart tube morphogenesis, looping, and maturation. We identified genetic regulatory modules driving crucial events of heart development that contained key cardiac TFs and are associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. Loss of function of cardiac TFs Gata5, Tbx5a, and Hand2 affected the cardiac regulatory networks and caused global changes in chromatin accessibility profile, indicating their role in heart development. Among regions with differential chromatin accessibility in mutants were highly conserved noncoding elements that represent putative enhancers driving heart development. The most prominent gene expression changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred between heart tube morphogenesis and looping, and were associated with metabolic shift and hematopoietic/cardiac fate switch during CM maturation. Our results revealed the dynamic regulatory landscape throughout heart development and identified interactive molecular networks driving key events of heart morphogenesis.

Establishing a whole blood CD4+ T cell immunity measurement to predict response to anti-PD-1.

BACKGROUND: Biomarkers that can accurately predict the efficacy of immune checkpoint inhibitors (ICIs) against programmed death 1 (PD-1) ligand in cancer immunotherapy are urgently needed. We have previously reported a novel formula that predicts the response to treatment with second-line nivolumab with high sensitivity and specificity in patients with non-small cell lung cancer (NSCLC) previously treated with chemotherapy. The formula was based on the percentages of CD62LlowCD4+ T cells (effector T cells; %Teff) and CD4+CD25+FOXP3+ T cells (regulatory T cells; %Treg) in the peripheral blood before treatment estimated using the peripheral blood mononuclear cell (PBMC) method. Here, we investigated the applicability of the formula (K-index) to predict the response to treatment with another ICI to expand its clinical applicability. Furthermore, we developed a simpler assay method based on whole blood (WB) samples to overcome the limitations of the PBMC method, such as technical difficulties, in obtaining the K-index. METHODS: The K-index was evaluated using the PBMC method in 59 patients with NSCLC who received first-line pembrolizumab treatment. We also assessed the K-index using the WB method and estimated the correlation between the measurements obtained using both methods in 76 patients with lung cancer. RESULTS: This formula consistently predicted the response to first-line pembrolizumab therapy in patients with NSCLC. The WB method correlated well with the PBMC method to obtain %Teff, %Treg, and the formula value. The WB method showed high repeatability (coefficient of variation, < 10%). The data obtained using WB samples collected in tubes containing either heparin or EDTA-2K and stored at room temperature (18-24 °C) for one day after blood sampling did not differ. Additionally, the performance of the WB method was consistent in different flow cytometry instruments. CONCLUSIONS: The K-index successfully predicted the response to first-line therapy with pembrolizumab, as reported earlier for the second-line therapy with nivolumab in patients with NSCLC. The WB method established in this study can replace the cumbersome PBMC method in obtaining the K-index. Overall, this study suggests that the K-index can predict the response to anti-PD-1 therapy in various cancers, including NSCLC.

Single-cell transcriptomics reveals expansion of cytotoxic CD4 T cells in supercentenarians.

Supercentenarians, people who have reached 110 y of age, are a great model of healthy aging. Their characteristics of delayed onset of age-related diseases and compression of morbidity imply that their immune system remains functional. Here we performed single-cell transcriptome analysis of 61,202 peripheral blood mononuclear cells (PBMCs), derived from 7 supercentenarians and 5 younger controls. We identified a marked increase of cytotoxic CD4 T cells (CD4 cytotoxic T lymphocytes [CTLs]) as a signature of supercentenarians. Furthermore, single-cell T cell receptor sequencing of 2 supercentenarians revealed that CD4 CTLs had accumulated through massive clonal expansion, with the most frequent clonotypes accounting for 15 to 35% of the entire CD4 T cell population. The CD4 CTLs exhibited substantial heterogeneity in their degree of cytotoxicity as well as a nearly identical transcriptome to that of CD8 CTLs. This indicates that CD4 CTLs utilize the transcriptional program of the CD8 lineage while retaining CD4 expression. Indeed, CD4 CTLs extracted from supercentenarians produced IFN-γ and TNF-α upon ex vivo stimulation. Our study reveals that supercentenarians have unique characteristics in their circulating lymphocytes, which may represent an essential adaptation to achieve exceptional longevity by sustaining immune responses to infections and diseases.

Efficacy and Safety of Anti-Programed Death-1 Blockade in Previously Treated Large-Cell Neuroendocrine Carcinoma.

BACKGROUND: Large-cell neuroendocrine carcinoma (LCNEC) of the lung is a rare tumor with an aggressive clinical course. However, there is limited knowledge of its treatment strategy. This retrospective study aimed to assess the efficacy and safety of anti-programed death-1 (PD-1) blockade monotherapy in previously treated advanced LCNEC. METHODS: Eleven patients with previously treated advanced LCNEC who received immune checkpoint inhibitor monotherapy between January 2015 and November 2020 were retrospectively analyzed for efficacy and safety. RESULTS: Of a total of 11 patients (median [range] age, 66 [37-79] years; 8 men [73%] and 3 women [27%]), 8 patients had performance status (PS) 0-1 [73%] and 3 patients had PS 2 [27%]; 9 patients received 1 prior chemotherapy [82%] and 2 patients received 2 prior chemotherapies [18%]. The median follow-up duration was 4.6 months. Although PD-1 blockade was administered at median cycles of 3 (range, 1-12), overall response rate, median progression-free survival, and median overall survival were 9.1%, 2.7 months, and 4.6 months, respectively. Any adverse events were observed in 9 patients (82%), including 1 patient with grade 3 pneumonitis as a serious adverse event. CONCLUSION: Anti-PD-1 blockade monotherapy as a subsequent line for previously treated advanced LCNEC exhibited usefulness and tolerability and was identified as a valid treatment option.

Prevention of hepatocellular carcinoma by targeting MYCN-positive liver cancer stem cells with acyclic retinoid.

Hepatocellular carcinoma (HCC) is a highly lethal cancer that has a high rate of recurrence, in part because of cancer stem cell (CSC)-dependent field cancerization. Acyclic retinoid (ACR) is a synthetic vitamin A-like compound capable of preventing the recurrence of HCC. Here, we performed a genome-wide transcriptome screen and showed that ACR selectively suppressed the expression of MYCN, a member of the MYC family of basic helix-loop-helix-zipper transcription factors, in HCC cell cultures, animal models, and liver biopsies obtained from HCC patients. MYCN expression in human HCC was correlated positively with both CSC and Wnt/ß-catenin signaling markers but negatively with mature hepatocyte markers. Functional analysis showed repressed cell-cycle progression, proliferation, and colony formation, activated caspase-8, and induced cell death in HCC cells following silencing of MYCN expression. High-content single-cell imaging analysis and flow cytometric analysis identified a MYCN+ CSC subpopulation in the heterogeneous HCC cell cultures and showed that these cells were selectively killed by ACR. Particularly, EpCAM+ cells isolated using a cell-sorting system showed increased MYCN expression and sensitivity to ACR compared with EpCAM- cells. In a long-term (>10 y) follow-up study of 102 patients with HCC, MYCN was expressed at higher levels in the HCC tumor region than in nontumor regions, and there was a positive correlation between MYCN expression and recurrence of de novo HCC but not metastatic HCC after curative treatment. In summary, these results suggest that MYCN serves as a prognostic biomarker and therapeutic target of ACR for liver CSCs in de novo HCC.

Raman Study on Lipid Droplets in Hepatic Cells Co-Cultured with Fatty Acids.

The purpose of the present study was to investigate molecular compositions of lipid droplets changing in live hepatic cells stimulated with major fatty acids in the human body, i.e., palmitic, stearic, oleic, and linoleic acids. HepG2 cells were used as the model hepatic cells. Morphological changes of lipid droplets were observed by optical microscopy and transmission electron microscopy (TEM) during co-cultivation with fatty acids up to 5 days. The compositional changes in the fatty chains included in the lipid droplets were analyzed via Raman spectroscopy and chemometrics. The growth curves of the cells indicated that palmitic, stearic, and linoleic acids induced cell death in HepG2 cells, but oleic acid did not. Microscopic observations suggested that the rates of fat accumulation were high for oleic and linoleic acids, but low for palmitic and stearic acids. Raman analysis indicated that linoleic fatty chains taken into the cells are modified into oleic fatty chains. These results suggest that the signaling pathway of cell death is independent of fat stimulations. Moreover, these results suggest that hepatic cells have a high affinity for linoleic acid, but linoleic acid induces cell death in these cells. This may be one of the causes of inflammation in nonalcoholic fatty liver disease (NAFLD).

Analysis of splice variants of the human protein disulfide isomerase (P4HB) gene.

BACKGROUND: Protein Disulfide Isomerases are thiol oxidoreductase chaperones from thioredoxin superfamily with crucial roles in endoplasmic reticulum proteostasis, implicated in many diseases. The family prototype PDIA1 is also involved in vascular redox cell signaling. PDIA1 is coded by the P4HB gene. While forced changes in P4HB gene expression promote physiological effects, little is known about endogenous P4HB gene regulation and, in particular, gene modulation by alternative splicing. This study addressed the P4HB splice variant landscape. RESULTS: Ten protein coding sequences (Ensembl) of the P4HB gene originating from alternative splicing were characterized. Structural features suggest that except for P4HB-021, other splice variants are unlikely to exert thiol isomerase activity at the endoplasmic reticulum. Extensive analyses using FANTOM5, ENCODE Consortium and GTEx project databases as RNA-seq data sources were performed. These indicated widespread expression but significant variability in the degree of isoform expression among distinct tissues and even among distinct locations of the same cell, e.g., vascular smooth muscle cells from different origins. P4HB-02, P4HB-027 and P4HB-021 were relatively more expressed across each database, the latter particularly in vascular smooth muscle. Expression of such variants was validated by qRT-PCR in some cell types. The most consistently expressed splice variant was P4HB-021 in human mammary artery vascular smooth muscle which, together with canonical P4HB gene, had its expression enhanced by serum starvation. CONCLUSIONS: Our study details the splice variant landscape of the P4HB gene, indicating their potential role to diversify the functional reach of this crucial gene. P4HB-021 splice variant deserves further investigation in vascular smooth muscle cells.

CAGE profiling of ncRNAs in hepatocellular carcinoma reveals widespread activation of retroviral LTR promoters in virus-induced tumors.

An increasing number of noncoding RNAs (ncRNAs) have been implicated in various human diseases including cancer; however, the ncRNA transcriptome of hepatocellular carcinoma (HCC) is largely unexplored. We used CAGE to map transcription start sites across various types of human and mouse HCCs with emphasis on ncRNAs distant from protein-coding genes. Here, we report that retroviral LTR promoters, expressed in healthy tissues such as testis and placenta but not liver, are widely activated in liver tumors. Despite HCC heterogeneity, a subset of LTR-derived ncRNAs were more than 10-fold up-regulated in the vast majority of samples. HCCs with a high LTR activity mostly had a viral etiology, were less differentiated, and showed higher risk of recurrence. ChIP-seq data show that MYC and MAX are associated with ncRNA deregulation. Globally, CAGE enabled us to build a mammalian promoter map for HCC, which uncovers a new layer of complexity in HCC genomics.

Discrimination analysis of excitatory and inhibitory neurons using Raman spectroscopy.

We have succeeded in discriminating between intact excitatory and inhibitory neuronal cells with Raman analysis. Excitatory and inhibitory neurons have several differences in their electric activities, but it can be difficult to determine their types based only on visual appearances. As Raman spectroscopy does not require any staining or labeling, its use in live neuronal cells is possible. In the present study, we used primary neurons obtained from rat cerebral cortexes, which we cultured on a glial feeder layered culturing dish for 15 days. The Raman spectra of the intact neurons on the dish were obtained; the neurons were then immunostained and their types determined. Partial least squares regression-discriminant analysis (PLSR-DA) was employed for classification of the excitatory and inhibitory neurons. The results demonstrated a high feasibility for use of Raman spectroscopy for discrimination analysis of inhibitory and excitatory neurons in a nondestructive manner.

Single-Nucleotide Resolution Mapping of Hepatitis B Virus Promoters in Infected Human Livers and Hepatocellular Carcinoma.

Hepatitis B virus (HBV) is a major cause of liver diseases, including hepatocellular carcinoma (HCC), and more than 650,000 people die annually due to HBV-associated liver failure. Extensive studies of individual promoters have revealed that heterogeneous RNA 5' ends contribute to the complexity of HBV transcriptome and proteome. Here, we provide a comprehensive map of HBV transcription start sites (TSSs) in human liver, HCC, and blood, as well as several experimental replication systems, at a single-nucleotide resolution. Using CAGE (cap analysis of gene expression) analysis of 16 HCC/nontumor liver pairs, we identify 17 robust TSSs, including a novel promoter for the X gene located in the middle of the gene body, which potentially produces a shorter X protein translated from the conserved second start codon, and two minor antisense transcripts that might represent viral noncoding RNAs. Interestingly, transcription profiles were similar in HCC and nontumor livers, although quantitative analysis revealed highly variable patterns of TSS usage among clinical samples, reflecting precise regulation of HBV transcription initiation at each promoter. Unlike the variety of TSSs found in liver and HCC, the vast majority of transcripts detected in HBV-positive blood samples are pregenomic RNA, most likely generated and released from liver. Our quantitative TSS mapping using the CAGE technology will allow better understanding of HBV transcriptional responses in further studies aimed at eradicating HBV in chronic carriers. IMPORTANCE: Despite the availability of a safe and effective vaccine, HBV infection remains a global health problem, and current antiviral protocols are not able to eliminate the virus in chronic carriers. Previous studies of the regulation of HBV transcription have described four major promoters and two enhancers, but little is known about their activity in human livers and HCC. We deeply sequenced the HBV RNA 5' ends in clinical human samples and experimental models by using a new, sensitive and quantitative method termed cap analysis of gene expression (CAGE). Our data provide the first comprehensive map of global TSS distribution over the entire HBV genome in the human liver, validating already known promoters and identifying novel locations. Better knowledge of HBV transcriptional activity in the clinical setting has critical implications in the evaluation of therapeutic approaches that target HBV replication.

Deficiency of multidrug resistance 2 contributes to cell transformation through oxidative stress.

Multidrug resistance 2 (Mdr2), also called adenosine triphosphate-binding cassette B4 (ABCB4), is the transporter of phosphatidylcholine (PC) at the canalicular membrane of mouse hepatocytes, which plays an essential role for bile formation. Mutations in human homologue MDR3 are associated with several liver diseases. Knockout of Mdr2 results in hepatic inflammation, liver fibrosis and hepatocellular carcinoma (HCC). Whereas the pathogenesis in Mdr2 (-/-) mice has been largely attributed to the toxicity of bile acids due to the absence of PC in the bile, the question of whether Mdr2 deficiency per se perturbs biological functions in the cell has been poorly addressed. As Mdr2 is expressed in many cell types, we used mouse embryonic fibroblasts (MEF) derived from Mdr2 (-/-) embryos to show that deficiency of Mdr2 increases reactive oxygen species accumulation, lipid peroxidation and DNA damage. We found that Mdr2 (-/-) MEFs undergo spontaneous transformation and that Mdr2 (-/-) mice are more susceptible to chemical carcinogen-induced intestinal tumorigenesis. Microarray analysis in Mdr2-/- MEFs and cap analysis of gene expression in Mdr2 (-/-) HCCs revealed extensively deregulated genes involved in oxidation reduction, fatty acid metabolism and lipid biosynthesis. Our findings imply a close link between Mdr2 (-/-) -associated tumorigenesis and perturbation of these biological processes and suggest potential extrahepatic functions of Mdr2/MDR3.