A Global Survey of Alternative Splicing of HBV Transcriptome Using Long-Read Sequencing.
Guiwen Guan,Jun Zou,Ting Zhang,Fengmin Lu,Xiangmei Chen
DOI: https://doi.org/10.1016/j.jhep.2021.07.023
IF: 25.7
2021-01-01
Journal of Hepatology
Abstract:Translatomic profiling reveals novel self-restricting virus-host interactions during HBV infectionJournal of HepatologyVol. 75Issue 1PreviewHBV remains a global threat to human health. It remains incompletely understood how HBV self-restricts in the host during most adult infections. Thus, we performed multi-omics analyses to systematically interrogate HBV-host interactions and the life cycle of HBV. Full-Text PDF Reply to: "A global survey of alternative splicing of HBV transcriptome using long-read sequencing"Journal of HepatologyVol. 76Issue 1PreviewThe unfolding multiomics landscapes of HBV-host interactions Full-Text PDF With great interest, we read the article "Translatomic profiling reveals novel self-restricting virus-host interactions during HBV infection" by Yuan et al. published in Journal of Hepatology.[1]Yuan S. Liao G. Zhang M. Zhu Y. Wang K. Xiao W. et al.Translatomic profiling reveals novel self-restricting virus-host interactions during HBV infection.J Hepatol. 2021; 75: 74-85https://doi.org/10.1016/j.jhep.2021.02.009Abstract Full Text Full Text PDF PubMed Scopus (8) Google Scholar By using RNA-seq and Ribo-seq, the authors found that alternative splicing and alternative translation initiation contributed to the translation of HpZ/P' and HBxZ, respectively, and thereby inhibited the transcription of covalently closed circular DNA. However, RNA-seq is deficient in determining transcript isoforms for short length reads and thus is not able to identify multispliced RNAs, as well as alternative transcription start sites (TSS) and transcription end sites (TES), especially for a complex genome like HBV with highly overlapping transcripts and open reading frames (ORFs).[2]Steijger T. Abril J.F. Engstrom P.G. Kokocinski F. Consortium R. Hubbard T.J. et al.Assessment of transcript reconstruction methods for RNA-seq.Nat Methods. 2013; 10: 1177-1184https://doi.org/10.1038/nmeth.2714Crossref PubMed Scopus (433) Google Scholar,[3]Boldogkoi Z. Moldovan N. Balazs Z. Snyder M. Tombacz D. Long-read sequencing - a powerful tool in viral transcriptome research.Trends Microbiol. 2019; 27: 578-592https://doi.org/10.1016/j.tim.2019.01.010Abstract Full Text Full Text PDF PubMed Scopus (35) Google Scholar These challenges can be overcome by long-read sequencing technology, such as Iso-seq, since it is able to provide full contig information regarding transcripts.[3]Boldogkoi Z. Moldovan N. Balazs Z. Snyder M. Tombacz D. Long-read sequencing - a powerful tool in viral transcriptome research.Trends Microbiol. 2019; 27: 578-592https://doi.org/10.1016/j.tim.2019.01.010Abstract Full Text Full Text PDF PubMed Scopus (35) Google Scholar Therefore, we aimed to explore the global transcriptome of HBV through Iso-seq to provide molecular evidence for the expression of HpZ/P' and HBxZ at the transcription level.We co-transfected the prcccDNA/pCMV-Cre plasmids (a gift from Professor Qiang Deng) into Huh-7 cells for 72 hours, at a similar timepoint to the previous work.[4]Qi Z. Li G. Hu H. Yang C. Zhang X. Leng Q. et al.Recombinant covalently closed circular hepatitis B virus DNA induces prolonged viral persistence in immunocompetent mice.J Virol. 2014; 88: 8045-8056https://doi.org/10.1128/JVI.01024-14Crossref PubMed Scopus (62) Google Scholar Then the total RNA was collected from the above transiently transfected cells, and sequenced by Iso-seq. A total of 1,059,497 full-length reads were obtained from 4.67 Gbp sequencing data, with lengths ranging from 531 to 26,451 bp (Fig. S1A,B). After mapping with the HBV genome (subtype ayw; GenBank accession no. V01460.1), we obtained 12,535 HBV-related full-length reads with lengths ranging between 1,140 to 12,931 bp (Fig. S1B). In addition to the classical 3.5, 2.4 and 2.1 kb HBV transcripts, we also detected abundant ultra-long HBV transcripts. Intriguingly, none of the 0.7 kb HBx transcript for X protein was detected (Fig. 1A, left lane). Among these ultra-long HBV transcripts, 3.9 kb and 6.5 kb transcripts were the major isoforms. Such ultra-long HBV transcripts were further verified by the northern blot assay (Fig. 1A, right lane). The 3.9 kb transcript contains 2 copies of HBx ORFs with TSS at 1,297-1,571 nt (Fig. S1C). It has been reported that the 0.7 kb HBx transcripts disappeared quickly after HBV infection,[5]Wu H.L. Chen P.J. Lin M.H. Chen D.S. Temporal aspects of major viral transcript expression in Hep G2 cells transfected with cloned hepatitis B virus DNA: with emphasis on the X transcript.Virology. 1991; 185: 644-651https://doi.org/10.1016/0042-6822(91)90535-jCrossref PubMed Scopus (0) Google Scholar and the 3.9 kb transcripts were able to translate the X protein.[6]Doitsh G. Shaul Y. A long HBV transcript encoding pX is inefficiently exported from the nucleus.Virology. 2003; 309: 339-349https://doi.org/10.1016/s0042-6822(03)00156-9Crossref PubMed Scopus (0) Google Scholar Therefore, we speculated that, in addition to the X protein, the HBxZ protein might also be mainly translated from the 3.9 kb HBx transcript. The 6.5 kb HBV transcript contains 2 complete copies of the HBV genome, and was named 6.5 kb PreC mRNA, since the TSS were concentrated around 1,785 nt, the classical TSS of PreC mRNA (Fig. S1D).Next, we analyzed the TSS and TES variants of the HBV transcripts obtained in this study. Four TSS peaks were detected and mapped to the promoter regions of CP/BCP, SP1, SP2, and XP, respectively (Fig. 1B). As for TES, 98.6% of these HBV transcripts ended at classical HBV TES (from 1,930 to 1,950 nt). While only 1.4% of them ended at non-classical HBV TES (Fig. 1B). Given the non-classical HBV TES are highly overlapped with HBV alternative splicing sites, we speculated that these non-classical TES of HBV might be caused by alternate termination, which is one type of alternative splicing event.[7]Wang E.T. Sandberg R. Luo S. Khrebtukova I. Zhang L. Mayr C. et al.Alternative isoform regulation in human tissue transcriptomes.Nature. 2008; 456: 470-476https://doi.org/10.1038/nature07509Crossref PubMed Scopus (3488) Google Scholar In addition, we also detected 3,591 transcripts with alternative splicing, which accounting for 28.6% total HBV RNA (Fig. 1B), consistent with the previously reported ratio of 30%.[8]Wu H.L. Chen P.J. Tu S.J. Lin M.H. Lai M.Y. Chen D.S. Characterization and genetic analysis of alternatively spliced transcripts of hepatitis B virus in infected human liver tissues and transfected HepG2 cells.J Virol. 1991; 65: 1680-1686https://doi.org/10.1128/JVI.65.4.1680-1686.1991Crossref PubMed Scopus (0) Google Scholar We found that 2,067 nt, 2,447 nt, 2,902 nt, 2,985 nt, 3,083 nt, 3,169 nt, 203 nt, 283 nt and 489 nt were the main splicing sites of HBV. The multiple combinations between these sites constituted the main splicing transcription isoforms of HBV (Fig. 1C). Similar to Yuan's study, we also found that 2,447/489 (the splicing editing between 2,447 nt and 489 nt) were the major HBV splicing variants, which accounting for 17.3% of spliced transcripts in this study (Fig. 1C). Of these 2,447/489 splicing variants, up to 97.8% were transcribed with TSS between 1,682 nt and 1,814 nt, prior to the start codon of the PreC protein. In addition to the translated initiation sites of PreC, Core and Polymerase protein, several more ATGs were located before the 2,446A2,447T489G formed by the 2,447/489 splicing variants. Therefore, the mechanism relevant to the translation of HpZ/P' remains to be clarified.In conclusion, our study explored the multispliced RNAs, alternative transcription initiation, and termination of HBV transcripts by Iso-seq, which provided molecular biology evidence for the expression of HpZ/P' and HBxZ suggested by Yuan and his colleagues.Financial supportThis work was supported by grants from Beijing Natural Science Foundation , China (No. 7212063 ), the Natural Science Foundation of China , China (Nos. 81572366 and 82072280 ), and the National S & T Major Project for Infectious Diseases of China , China (No. 2017ZX10202203 ).Authors' contributionsChen X, Lu F designed the research. Guan G, Zou J, Zhang T performed the experiments. All authors analyzed the data. Guan G wrote the paper. Chen X and Lu F revised the paper.Data availability statementThe data generated in this study have been deposited in Sequence Read Archive (SRA) repository (www.ncbi.nih.gov.sra) with accession number of PRJNA742195. With great interest, we read the article "Translatomic profiling reveals novel self-restricting virus-host interactions during HBV infection" by Yuan et al. published in Journal of Hepatology.[1]Yuan S. Liao G. Zhang M. Zhu Y. Wang K. Xiao W. et al.Translatomic profiling reveals novel self-restricting virus-host interactions during HBV infection.J Hepatol. 2021; 75: 74-85https://doi.org/10.1016/j.jhep.2021.02.009Abstract Full Text Full Text PDF PubMed Scopus (8) Google Scholar By using RNA-seq and Ribo-seq, the authors found that alternative splicing and alternative translation initiation contributed to the translation of HpZ/P' and HBxZ, respectively, and thereby inhibited the transcription of covalently closed circular DNA. However, RNA-seq is deficient in determining transcript isoforms for short length reads and thus is not able to identify multispliced RNAs, as well as alternative transcription start sites (TSS) and transcription end sites (TES), especially for a complex genome like HBV with highly overlapping transcripts and open reading frames (ORFs).[2]Steijger T. Abril J.F. Engstrom P.G. Kokocinski F. Consortium R. Hubbard T.J. et al.Assessment of transcript reconstruction methods for RNA-seq.Nat Methods. 2013; 10: 1177-1184https://doi.org/10.1038/nmeth.2714Crossref PubMed Scopus (433) Google Scholar,[3]Boldogkoi Z. Moldovan N. Balazs Z. Snyder M. Tombacz D. Long-read sequencing - a powerful tool in viral transcriptome research.Trends Microbiol. 2019; 27: 578-592https://doi.org/10.1016/j.tim.2019.01.010Abstract Full Text Full Text PDF PubMed Scopus (35) Google Scholar These challenges can be overcome by long-read sequencing technology, such as Iso-seq, since it is able to provide full contig information regarding transcripts.[3]Boldogkoi Z. Moldovan N. Balazs Z. Snyder M. Tombacz D. Long-read sequencing - a powerful tool in viral transcriptome research.Trends Microbiol. 2019; 27: 578-592https://doi.org/10.1016/j.tim.2019.01.010Abstract Full Text Full Text PDF PubMed Scopus (35) Google Scholar Therefore, we aimed to explore the global transcriptome of HBV through Iso-seq to provide molecular evidence for the expression of HpZ/P' and HBxZ at the transcription level. We co-transfected the prcccDNA/pCMV-Cre plasmids (a gift from Professor Qiang Deng) into Huh-7 cells for 72 hours, at a similar timepoint to the previous work.[4]Qi Z. Li G. Hu H. Yang C. Zhang X. Leng Q. et al.Recombinant covalently closed circular hepatitis B virus DNA induces prolonged viral persistence in immunocompetent mice.J Virol. 2014; 88: 8045-8056https://doi.org/10.1128/JVI.01024-14Crossref PubMed Scopus (62) Google Scholar Then the total RNA was collected from the above transiently transfected cells, and sequenced by Iso-seq. A total of 1,059,497 full-length reads were obtained from 4.67 Gbp sequencing data, with lengths ranging from 531 to 26,451 bp (Fig. S1A,B). After mapping with the HBV genome (subtype ayw; GenBank accession no. V01460.1), we obtained 12,535 HBV-related full-length reads with lengths ranging between 1,140 to 12,931 bp (Fig. S1B). In addition to the classical 3.5, 2.4 and 2.1 kb HBV transcripts, we also detected abundant ultra-long HBV transcripts. Intriguingly, none of the 0.7 kb HBx transcript for X protein was detected (Fig. 1A, left lane). Among these ultra-long HBV transcripts, 3.9 kb and 6.5 kb transcripts were the major isoforms. Such ultra-long HBV transcripts were further verified by the northern blot assay (Fig. 1A, right lane). The 3.9 kb transcript contains 2 copies of HBx ORFs with TSS at 1,297-1,571 nt (Fig. S1C). It has been reported that the 0.7 kb HBx transcripts disappeared quickly after HBV infection,[5]Wu H.L. Chen P.J. Lin M.H. Chen D.S. Temporal aspects of major viral transcript expression in Hep G2 cells transfected with cloned hepatitis B virus DNA: with emphasis on the X transcript.Virology. 1991; 185: 644-651https://doi.org/10.1016/0042-6822(91)90535-jCrossref PubMed Scopus (0) Google Scholar and the 3.9 kb transcripts were able to translate the X protein.[6]Doitsh G. Shaul Y. A long HBV transcript encoding pX is inefficiently exported from the nucleus.Virology. 2003; 309: 339-349https://doi.org/10.1016/s0042-6822(03)00156-9Crossref PubMed Scopus (0) Google Scholar Therefore, we speculated that, in addition to the X protein, the HBxZ protein might also be mainly translated from the 3.9 kb HBx transcript. The 6.5 kb HBV transcript contains 2 complete copies of the HBV genome, and was named 6.5 kb PreC mRNA, since the TSS were concentrated around 1,785 nt, the classical TSS of PreC mRNA (Fig. S1D). Next, we analyzed the TSS and TES variants of the HBV transcripts obtained in this study. Four TSS peaks were detected and mapped to the promoter regions of CP/BCP, SP1, SP2, and XP, respectively (Fig. 1B). As for TES, 98.6% of these HBV transcripts ended at classical HBV TES (from 1,930 to 1,950 nt). While only 1.4% of them ended at non-classical HBV TES (Fig. 1B). Given the non-classical HBV TES are highly overlapped with HBV alternative splicing sites, we speculated that these non-classical TES of HBV might be caused by alternate termination, which is one type of alternative splicing event.[7]Wang E.T. Sandberg R. Luo S. Khrebtukova I. Zhang L. Mayr C. et al.Alternative isoform regulation in human tissue transcriptomes.Nature. 2008; 456: 470-476https://doi.org/10.1038/nature07509Crossref PubMed Scopus (3488) Google Scholar In addition, we also detected 3,591 transcripts with alternative splicing, which accounting for 28.6% total HBV RNA (Fig. 1B), consistent with the previously reported ratio of 30%.[8]Wu H.L. Chen P.J. Tu S.J. Lin M.H. Lai M.Y. Chen D.S. Characterization and genetic analysis of alternatively spliced transcripts of hepatitis B virus in infected human liver tissues and transfected HepG2 cells.J Virol. 1991; 65: 1680-1686https://doi.org/10.1128/JVI.65.4.1680-1686.1991Crossref PubMed Scopus (0) Google Scholar We found that 2,067 nt, 2,447 nt, 2,902 nt, 2,985 nt, 3,083 nt, 3,169 nt, 203 nt, 283 nt and 489 nt were the main splicing sites of HBV. The multiple combinations between these sites constituted the main splicing transcription isoforms of HBV (Fig. 1C). Similar to Yuan's study, we also found that 2,447/489 (the splicing editing between 2,447 nt and 489 nt) were the major HBV splicing variants, which accounting for 17.3% of spliced transcripts in this study (Fig. 1C). Of these 2,447/489 splicing variants, up to 97.8% were transcribed with TSS between 1,682 nt and 1,814 nt, prior to the start codon of the PreC protein. In addition to the translated initiation sites of PreC, Core and Polymerase protein, several more ATGs were located before the 2,446A2,447T489G formed by the 2,447/489 splicing variants. Therefore, the mechanism relevant to the translation of HpZ/P' remains to be clarified. In conclusion, our study explored the multispliced RNAs, alternative transcription initiation, and termination of HBV transcripts by Iso-seq, which provided molecular biology evidence for the expression of HpZ/P' and HBxZ suggested by Yuan and his colleagues. Financial supportThis work was supported by grants from Beijing Natural Science Foundation , China (No. 7212063 ), the Natural Science Foundation of China , China (Nos. 81572366 and 82072280 ), and the National S & T Major Project for Infectious Diseases of China , China (No. 2017ZX10202203 ). This work was supported by grants from Beijing Natural Science Foundation , China (No. 7212063 ), the Natural Science Foundation of China , China (Nos. 81572366 and 82072280 ), and the National S & T Major Project for Infectious Diseases of China , China (No. 2017ZX10202203 ). Authors' contributionsChen X, Lu F designed the research. Guan G, Zou J, Zhang T performed the experiments. All authors analyzed the data. Guan G wrote the paper. Chen X and Lu F revised the paper. Chen X, Lu F designed the research. Guan G, Zou J, Zhang T performed the experiments. All authors analyzed the data. Guan G wrote the paper. Chen X and Lu F revised the paper. Data availability statementThe data generated in this study have been deposited in Sequence Read Archive (SRA) repository (www.ncbi.nih.gov.sra) with accession number of PRJNA742195. The data generated in this study have been deposited in Sequence Read Archive (SRA) repository (www.ncbi.nih.gov.sra) with accession number of PRJNA742195. The authors declare no conflicts of interest that pertain to this work. Please refer to the accompanying ICMJE disclosure forms for further details. Supplementary dataThe following is the supplementary data to this article: Download .pdf (.32 MB) Help with pdf files Multimedia component 1 Download .pdf (.19 MB) Help with pdf files Multimedia component 2 Download .pdf (5.56 MB) Help with pdf files Multimedia component 3 The following is the supplementary data to this article: Download .pdf (.32 MB) Help with pdf files Multimedia component 1 Download .pdf (.19 MB) Help with pdf files Multimedia component 2 Download .pdf (5.56 MB) Help with pdf files Multimedia component 3