Identification of RNA editing sites in the SNP database

The relationship between human inherited genomic variations and phenotypic differences has been the focus of much research effort in recent years. These studies benefit from millions of single-nucleotide polymorphism (SNP) records available in public databases, such as dbSNP. The importance of identifying false dbSNP records increases with the growing role played by SNPs in linkage analysis for disease traits. In particular, the emerging understanding of the abundance of DNA and RNA editing calls for a careful distinction between inherited SNPs and somatic DNA and RNA modifications. In order to demonstrate that some of the SNP database records are actually somatic modification, we focus on one type of these modifications, namely A-to-I RNA editing, and present evidence for hundreds of dbSNP records that are actually editing sites. We provide a list of 102 RNA editing sites previously annotated in dbSNP database as SNPs, and experimentally validate seven of these. Interestingly, we show how dbSNP can serve as a starting point to look for new editing sites. Our results, for this particular type of RNA editing, demonstrate the need for a careful analysis of SNP databases in light of the increasing recognition of the significance of somatic sequence modifications.     

Influence of gag and RRE Sequences on HIV-1 RNA Packaging Signal Structure and Function

The packaging signal (Ψ) and Rev-responsive element (RRE) enable unspliced HIV-1 RNAs' export from the nucleus and packaging into virions. For some retroviruses, engrafting Ψ onto a heterologous RNA is sufficient to direct encapsidation. In contrast, HIV-1 RNA packaging requires 5′ leader Ψ elements plus poorly defined additional features. We previously defined minimal 5′ leader sequences competitive with intact Ψ for HIV-1 packaging, and here examined the potential roles of additional downstream elements. The findings confirmed that together, HIV-1 5′ leader Ψ sequences plus a nuclear export element are sufficient to specify packaging. However, RNAs trafficked using a heterologous export element did not compete well with RNAs using HIV-1's RRE. Furthermore, some RNA additions to well-packaged minimal vectors rendered them packaging-defective. These defects were rescued by extending gag sequences in their native context. To understand these packaging defects' causes, in vitro dimerization properties of RNAs containing minimal packaging elements were compared to RNAs with sequence extensions that were or were not compatible with packaging. In vitro dimerization was found to correlate with packaging phenotypes, suggesting that HIV-1 evolved to prevent 5′ leader residues' base pairing with downstream residues and misfolding of the packaging signal. Our findings explain why gag sequences have been implicated in packaging and show that RRE's packaging contributions appear more specific than nuclear export alone. Paired with recent work showing that sequences upstream of Ψ can dictate RNA folds, the current work explains how genetic context of minimal packaging elements contributes to HIV-1 RNA fate determination.     

基于全基因组重测序技术的‘红叶’杜仲SNP位点开发

为了挖掘与‘红叶’杜仲（Eucommia ulmoides ‘Hongye’）红叶性状紧密联系的SNP位点,进一步揭示红叶性状的遗传基础和分子机理。以‘红叶’杜仲和普通绿叶杜仲‘小叶’杜仲（Eucommia ulmoides ‘Xiaoye’）为研究材料,进行覆盖深度约为10x的全基因组重测序。使用SnpEff软件预测变异位点对蛋白编码的影响,结合花色苷的代谢通路和关键酶基因,筛选与‘红叶’杜仲叶色形成相关的差异位点。利用Sanger测序二代测序筛选的SNP位点,分子标记验证群体是‘红叶’杜仲和‘小叶’杜仲。结果表明,‘红叶’杜仲测序产生Clean data为14.16 Gb,‘小叶’杜仲产生Clean data为14.29 Gb。在‘红叶’杜仲中注释到严重影响蛋白质功能的有1 516个SNP,中度影响的41 328个SNP,在‘小叶’杜仲中存在严重影响蛋白质功能的SNP为1 640个,中度影响功能的SNP为47 192个。测得26 722条基因中有228条基因是与花色苷或类黄酮合成相关的酶基因。经过筛选,确定了12个特异性的SNP位点,均属于外显子区域的错义突变。利用一代测序验证,根据SNP位置设计了7对引物,SNP准确率达到100%。     

Identification and remediation of biases in the activity of RNA ligases in small-RNA deep sequencing

Deep sequencing of small RNAs (sRNA-seq) is now the gold standard for small RNA profiling and discovery. Biases in sRNA-seq have been reported, but their etiology remains unidentified. Through a comprehensive series of sRNA-seq experiments, we establish that the predominant cause of the bias is the RNA ligases. We further demonstrate that RNA ligases have strong sequence-specific biases which distort the small RNA profiles considerably. We have devised a pooled adapter strategy to overcome this bias, and validated the method through data derived from microarray and qPCR. In light of our findings, published small RNA profiles, as well as barcoding strategies using adapter-end modifications, may need to be revisited. Importantly, by providing a wide spectrum of substrate for the ligase, the pooled-adapter strategy developed here provides a means to overcome issues of bias, and generate more accurate small RNA profiles.     

Compositional complementarity between genomic RNA and coat proteins in positive-sense single-stranded RNA viruses

During packaging in positive-sense single-stranded RNA (+ssRNA) viruses, coat proteins (CPs) interact directly with multiple regions in genomic RNA (gRNA), but the underlying physicochemical principles remain unclear. Here we analyze the high-resolution cryo-EM structure of bacteriophage MS2 and show that the gRNA/CP binding sites, including the known packaging signal, overlap significantly with regions where gRNA nucleobase-density profiles match the corresponding CP nucleobase-affinity profiles. Moreover, we show that the MS2 packaging signal corresponds to the global minimum in gRNA/CP interaction energy in the unstructured state as derived using a linearly additive model and knowledge-based nucleobase/amino-acid affinities. Motivated by this, we predict gRNA/CP interaction sites for a comprehensive set of 1082 +ssRNA viruses. We validate our predictions by comparing them with site-resolved information on gRNA/CP interactions derived in SELEX and CLIP experiments for 10 different viruses. Finally, we show that in experimentally studied systems CPs frequently interact with autologous coding regions in gRNA, in agreement with both predicted interaction energies and a recent proposal that proteins in general tend to interact with own mRNAs, if unstructured. Our results define a self-consistent framework for understanding packaging in +ssRNA viruses and implicate interactions between unstructured gRNA and CPs in the process.     

Analysis of context-dependent errors for illumina sequencing

The new generation of short-read sequencing technologies requires reliable measures of data quality. Such measures are especially important for variant calling. However, in the particular case of SNP calling, a great number of false-positive SNPs may be obtained. One needs to distinguish putative SNPs from sequencing or other errors. We found that not only the probability of sequencing errors (i.e. the quality value) is important to distinguish an FP-SNP but also the conditional probability of "correcting" this error (the "second best call" probability, conditional on that of the first call). Surprisingly, around 80% of mismatches can be "corrected" with this second call. Another way to reduce the rate of FP-SNPs is to retrieve DNA motifs that seem to be prone to sequencing errors, and to attach a corresponding conditional quality value to these motifs. We have developed several measures to distinguish between sequence errors and candidate SNPs, based on a base call's nucleotide context and its mismatch type. In addition, we suggested a simple method to correct the majority of mismatches, based on conditional probability of their "second" best intensity call. We attach a corresponding second call confidence (quality value) of being corrected to each mismatch.     

基于转录组测序数据识别黑猩猩RNA编辑位点

使用转录组测序(RNA-Seq)数据识别黑猩猩RNA编辑位点,探索了RNA编辑的识别机制以及潜在的功能影响．基于黑猩猩RNA-Seq数据与基因组序列的比对信息发现RNA-DNA错配位点,并构建编辑位点候选集．从中滤除基因组或转录组测序质量低的位点,其他的过滤条件包括3′端测不准、覆盖度、SNP位点以及估算的编辑水平．构建二项分布统计模型和Bonferroni多重检验滤除候选集中的随机错误,得到RNA编辑位点．选取落在已知基因上的编辑位点进行功能分析,并用Two Sample Logo软件分析编辑位点上下游序列的特征．识别出黑猩猩12种碱基替换型RNA编辑位点8 334个,其中有41个编辑位点改变原有的氨基酸,另有3个编辑位点落在microRNA(miRNA)潜在靶基因的种子结合区．统计学分析表明,分别有640和872个RNA编辑位点存在组织和性别差异．上下游碱基频率分析表明,多种类型的编辑位点紧邻碱基具有显著偏好．结果显示, RNA编辑在黑猩猩体内大量存在,且潜在具有重要的生物学功能,为进一步深入研究灵长类RNA编辑的机制奠定了基础．     

Comparative structural effects of HIV-1 Gag and nucleocapsid proteins in binding to and unwinding of the viral RNA packaging signal

The major RNA binding region of the HIV-1 Gag polyprotein is the nucleocapsid (NC) domain, which is responsible for the specific capture of the genomic RNA genome during viral assembly. The Gag polyprotein has other RNA chaperone functions, which are mirrored by the isolated NC protein after physiological cleavage from Gag. Gag, however, is suggested to have superior nucleic acid chaperone activity. Here we investigate the interaction of Gag and NC with the core RNA structure of the HIV-1 packaging signal (Ψ), using 2-aminopurine substitution to create a series of modified RNAs based on the Ψ helix loop structure. The effects of 2-aminopurine substitution on the physical and structural properties of the viral Ψ were characterized. The fluorescence properties of the 2-aminopurine substitutions showed features consistent with the native GNAR tetraloop. Dissociation constants (K(d)) of the two viral proteins, measured by fluorescence polarization (FP), were similar, and both NC and Gag affected the 2-aminopurine fluorescence of bases close to the loop binding region in a similar fashion. However, the influence of Gag on the fluorescence of the 2-aminopurine nucleotides at the base of the helix implied a much more potent helix destabilizing action on the RNA stem loop (SL) versus that seen with NC. This was further supported when the viral Ψ SL was tagged with a 5' fluorophore and 3' quencher. In the absence of any viral protein, minimal fluorescence was detected; addition of NC yielded a slight increase in fluorescence, while addition of the Gag protein yielded a large change in fluorescence, further suggesting that, compared to NC, the Gag protein has a greater propensity to affect RNA structure and that Ψ helix unwinding may be an intrinsic step in RNA encapsidation.