严格去除基因组DNA（gDNA）后HCV感染Huh-7肝癌细胞的转录组分析：低丰度lncRNA检测及宿主基因表达重评估

《Viruses》：Transcriptome Analysis of Hepatitis C Virus-Infected Hepatoma Cells Following Stringent Removal of Genomic DNA

丙型肝炎病毒（HCV）是一种正链RNA病毒（positive-strand RNA virus），感染肝细胞并重编程细胞基因表达，进而促进慢性肝病和肝细胞癌（hepatocellular carcinoma, HCC）。HCV复制缓慢，对细胞转录组（transcriptome）诱导的变化有限。在此，研究人员改进了对低丰度转录本（low-abundance transcripts），包括长非编码RNA（long non-coding RNA, lncRNA）在内的分析。研究人员通过TRIzol、酸性酚、硅胶基质柱和重组DNase I的组合，从Huh-7肝癌细胞系总RNA中严格去除基因组DNA（genomic DNA, gDNA）。RNA测序（RNA sequencing, RNA-Seq）和靶向验证被用于表征Lnc-ITM2C-1剪接变体（splice variants）及其邻近基因G蛋白偶联受体55（G protein-coupled receptor 55, GPR55）的表达，以研究它们在HCV复制中的潜在作用。与研究人员先前的发现相反，在改进的实验条件下，Lnc-ITM2C-1和GPR55既不参与与HCV的正调控回路，也未诱导可检测的转录组变化。根据HCV感染状态对RNA-Seq数据的二次分析显示，HCV诱导21个宿主细胞基因出现统计学稳健的表达，涉及应激反应（stress response，尤其ATF3）、脂质代谢（lipid metabolism）、细胞骨架组织（cytoskeleton organization）和生长调节（growth regulation）等通路。相反，研究人员发现先前报道的HCV相关lncRNA中仅有限子集（例如EGOT）表达发生改变。这凸显了严格去除gDNA对于分析低丰度转录本如lncRNA的必要性。

论文主体内容总结：HCV是一种正链RNA病毒，感染肝细胞后可建立复制复合体并重编程细胞基因表达，约70%感染转为慢性，常导致肝硬化和肝癌。HCV复制缓慢，感染后约2天才达到病毒RNA复制高峰，细胞培养中每细胞正链RNA约10⁴分子，低于许多其他病毒；HCV感染引起的细胞差异表达基因（differentially expressed genes, DEGs）数量也有限，常为数十至数百个，取决于感染时间、实验设置及倍数变化和显著性阈值。研究人员先前报道Lnc-ITM2C-1及其邻近基因GPR55在HCV感染细胞中差异表达，并提出正调控回路：HCV刺激Lnc-ITM2C-1，Lnc-ITM2C-1刺激GPR55 mRNA，GPR55抑制干扰素刺激基因（interferon-stimulated gene, ISG）表达，从而提高HCV复制容许性。然而，低丰度转录本分析易受总RNA中残留gDNA污染影响，尤其qPCR中无逆转录酶对照（noRT）出现假阳性可干扰低丰度基因定量。因此，研究人员改进总RNA纯化和gDNA去除方案，重新评估该调控轴，鉴定Huh-7细胞中Lnc-ITM2C-1和GPR55剪接变体，并重新分析HCV感染对细胞转录组及既往报道lncRNA的影响。研究结论是：在严格去除gDNA的改进条件下，未证实Lnc-ITM2C-1—GPR55正调控回路；HCV感染主要稳健诱导有限数量的应激、脂质代谢、细胞骨架和生长调控相关基因；低丰度lncRNA分析必须严格去除gDNA并结合严格验证。其意义在于为低丰度转录本，尤其lncRNA的可靠检测提供技术警示和方法学依据。

关键技术方法概括（不超过250字）：研究采用Huh-7.5细胞电转HCV J6/JFH1嵌合基因2a型（Jc1）全长RNA以产生感染性病毒，并以Huh-7-Lunet/CD81细胞进行HCV感染（感染复数（multiplicity of infection, MOI）0.3，感染后48 h（hpi）为主）及质粒过表达。RNA-Seq样本共24个：12个HCV感染与12个未感染样本，各含Lnc-ITM2C-1 SP8、SP12、GPR55 SP1和空载体4组、每组3个生物学重复；另设RT-qPCR时间点验证。关键技术包括：TRIzol联合酸性酚、硅胶柱和重组DNase I的严格gDNA去除；剪接变体特异性引物、noRT对照、熔解曲线及凝胶产物大小验证的逆转录定量实时PCR（reverse transcription quantitative real-time PCR, RT-qPCR）；rRNA去除RNA-Seq、唯一分子标识符（unique molecular identifier, UMI）去重、HISAT2比对和DESeq2差异分析；以及HCV NS3蛋白/RNA、焦点形成单位（focus-forming units, FFU）和过表达验证。

研究结果：
3.1 严格去除gDNA对低丰度基因可靠分析至关重要（Improved Removal of gDNA Is Essential for Reliable Analysis of Low-Abundance Genes）。通过比较标准TRIzol纯化与改进方案，并结合noRT对照、熔解曲线和凝胶产物大小验证，研究人员发现标准纯化可出现noRT扩增、宽熔解峰及错误大小产物，提示gDNA污染和非特异假阳性；改进后noRT阴性、熔解曲线单峰且凝胶显示正确大小产物。对低丰度靶标，仅采用noRT至40个循环无扩增的引物组。结论是低丰度转录本检测必须严格去除gDNA并验证qPCR产物特异性。
3.2 Huh-7肝癌细胞中Lnc-ITM2C-1和GPR55 mRNA剪接变体的表达（Expression of Lnc-ITM2C-1 and GPR55 mRNA Splice Variants in Huh-7 Hepatoma Cells）。通过剪接变体特异性引物RT-qPCR，研究人员明确鉴定Lnc-ITM2C-1的SP7和SP12在Huh-7细胞表达；SP8因缺乏独特序列，主要依据其NCBI RefSeq代表性和共享引物间接证据认为可能表达；GPR55的SP1和SP3被鉴定表达，SP3的5′端延长序列得到确认。结论是Huh-7细胞中Lnc-ITM2C-1和GPR55存在特定剪接变体表达。
3.3 HCV感染特异性上调的基因（Genes Specifically Upregulated by HCV Infection）。通过24个RNA-Seq样本分析，研究人员发现HCV感染后显著上调基因很少，ATF3最突出，CHAC1、TUBA1A及假基因ASNSP1、ROCK1P1也上调；Lnc-ITM2C-1或GPR55过表达未改变上述HCV应答。先前研究中14个显著上调基因本次也上调，包括PPP1R15A（GADD34）、BMP2、DAB2、DUSP1、JUN、GDF15、TRIB3（SINK）等。结论是HCV诱导的细胞基因应答有限且可重复，过表达Lnc-ITM2C-1或GPR55不改变该应答。
3.4 Lnc-ITM2C-1和GPR55不参与HCV感染期间的正调控回路（Lnc-ITM2C-1 and GPR55 Are Not Involved in a Positive Regulatory Circuit During HCV Infection）。通过RNA-Seq和RT-qPCR验证，研究人员未发现HCV感染显著刺激Lnc-ITM2C-1或GPR55表达；免疫和干扰素相关基因中仅CXCL6、IFNAR2轻微且低显著性变化，既往报道受HCV下调的免疫相关基因未下调。过表达Lnc-ITM2C-1或GPR55未促进HCV复制、FFU、NS3蛋白或HCV RNA；Lnc-ITM2C-1与GPR55互不刺激表达，也不影响邻近ITM2C。结论是不支持先前提出的HCV—Lnc-ITM2C-1—GPR55正调控回路。
3.5 Lnc-ITM2C-1和GPR55表达对细胞转录组的影响（Analysis of the Effect of Lnc-ITM2C-1 and GPR55 Expression on the Cellular Transcriptome）。通过RNA-Seq过表达组分析，研究人员发现除过表达基因外，仅Lnc-ITM2C-1 SP8过表达时RNU5B-1下调；其他处理在未感染和HCV感染细胞中均无显著DEG。结论是所测试的Lnc-ITM2C-1或GPR55 RNA对Huh-7-Lunet/CD81细胞转录组无明显影响。
3.6 HCV感染期间基因表达变化的严格再分析（Stringent Re-Analysis of Gene Expression Changes During HCV Infection）。通过合并12个未感染与12个HCV感染样本进行二次分析，统计功效提高；在padj < 10^-10且倍数变化>1的严格阈值下，发现21个已知上调基因而非仅ATF3。ATF3约3倍上调且显著性突出；新增INHBE、PLA2G4C、RasGEF1B、NEURL3、KSR2及TUBA1A、MYL7、FBN1、KRT80、FUT1、SOCS2等，涉及应激、脂质代谢、细胞骨架和生长信号；另有15个已知基因达padj < 10^-5，SLC6A12下调。主成分分析（principal component analysis, PCA）显示感染与未感染样本在top基因中分离，而过表达lncRNA或GPR55无影响。结论是HCV感染引起有限但稳健的宿主基因表达变化。
3.7 HCV感染期间其他lncRNA表达变化（Changes in Expression of Other lncRNAs During HCV Infection）。通过同一合并数据集分析18个既往报道与HCV复制相关的lncRNA，研究人员发现EGOT显著上调约1.75倍，LINC01419轻度但显著下调；HEIC、HOTAIR、HOTTIP、HULC、LUARIS、MALAT1、NORAD和UCA1未显著上调。结论是仅有限lncRNA表达变化可在改进条件下重复。

讨论部分总结：讨论指出，HCV复制缓慢并可持续数十年；感染后48 hpi病毒基因组读段仅约占样本总读段2%，约80个细胞基因显著上调或下调超过2倍。上调基因可归为四类功能：内质网（endoplasmic reticulum, ER）应激、未折叠蛋白反应（unfolded protein response, UPR）和氧化还原平衡相关基因（ATF3、CHAC1、DDIT3/CHOP、TRIB3、PPP1R15A、INHBE）；脂质代谢调控基因（PLA2G4C、INHBE、GDF15）；细胞骨架组织相关基因（TUBA1A、MYL7、KRT80、RASGEF1B）；生长调控和信号相关基因（JUN、KSR2、NEURL3、SOCS2）。原文认为这种平衡应答可能有助于HCV低水平复制和持续存在，而不导致严重细胞死亡或完全病毒清除。讨论强调，SOCS2、PLA2G4C、ASNS、NEURL3和KSR2等显著变化基因属低丰度表达，log₂读段计数约2至7，RPKM（每千碱基转录本每百万比对读段，reads per kilobase of transcript per million mapped reads）约0.01至0.75，若仅分析RPKM>1基因将被忽略。技术层面，三步额外gDNA去除使gDNA污染降低至少9或10个log₂，约500至1000倍；结合随机RT引物、合适退火条件和经验证引物，可可靠检测Lnc-ITM2C-1和GPR55剪接变体。与Liu等报道一致，低gDNA污染样本呈双峰读段分布并保留低表达基因群；该研究样本呈双峰分布，Lnc-ITM2C-1平均仅2个读段（log₂读段计数=1）。旧2021年RNA-Seq数据中低丰度lncRNA在部分或全部重复中检测不到，而当前数据中图11所列lncRNA在所有重复中均有读段，Lnc-ITM2C-1仅一个重复为0。原文指出，先前回路未被证实可能与旧方案TRIzol后直接DNase I处理有关，酚—水混合水相约含7%酚，残留酚可能抑制DNase I，使gDNA污染波动影响结果；新方案在DNase处理前用硅胶柱去除残留酚，并采用三步纯化。低表达范围内IGF2-AS、HOTAIR、LNCATV、HOTTIP和LUARIS在24个生物学样本中一致检测，提示这些低表达值反映真实RNA表达而非gDNA污染导致的技术波动。局限为MOI 0.3可能降低效应强度；但采用温和阈值仍检出477个DEGs，与其他MOI研究数量和通路可比。研究结论翻译：总之，对于通过RNA-Seq或RT-qPCR研究低水平表达基因，研究人员强烈建议严格去除gDNA（如采用改进的RNA纯化方案），并结合严格验证（noRT对照和检查正确产物大小的凝胶）。这对低丰度表达长非编码RNA和其他基因的研究可能特别重要。