排版相对好一些的原文链接如下（微信公众号）
https://mp.weixin.qq.com/s?__biz=MzI4MzU2Mjc1OA==&mid=2247484670&idx=1&sn=f4807aa5c5b4b79ec3c75416a4aeecf2&chksm=eb898297dcfe0b81a2113e2e22c6548eba7d3b5e0678541a8e8a7e161810d6b9a7660e9ff5fe&token=273967589&lang=zh_CN#rd
Identification of Diverse Mycovirusesthrough Metatranscriptomics Characterization of the Viromes of Five MajorFungal Plant Pathogens
分析真菌宏转录组寻找病毒
摘要
该研究对五种主要植物病原真菌（Colletotrichum truncatum（2株）, Diaporthe longicolla（3株）, Macrophomina phaseolina（48株）, Rhizoctonia solani（84株）和Sclerotinia sclerotiorum（138株））的275株菌株的total RNA进行深度测序。测序结果展示了植物病原真菌中病毒的多样性和分类，以及病毒与它们的寄主真菌的协同进化等信息。
知识背景
目前人们对真菌病毒的了解不到1%，通过病毒宏基因组（metagenomes，metatranscriptomic viromes）分析是发现真菌病毒多样性的主要渠道。单个真菌病毒对真菌和生态的作用我们知之甚少，因为一个真菌常常有多种病毒感染，很难去单个研究。
随着真菌转录组测序增多，人们对真菌和其病毒的起源和进化了解增多，一些植物内生病毒很有可能来源于真菌病毒，且两者之间遗传物质的交换可能产生新的有害的植物病毒。比如植物病毒Amalgaviridae科病毒通过一种dsRNA真菌病毒和负链RNA（nsRNA）病毒的重组产生。
研究发现一些病毒可以降低真菌致病力（capable of attenuating virulence），可作为防治真菌病害的策略。相同的策略已经被用在治疗人类细菌疾病中，人们利用噬菌体和噬菌体编码的蛋白质（比如内溶素）防治细菌病原菌。同样的，我们也可以考虑用真菌病毒编码的产物，干扰真菌RNA沉默，小RNA代谢和基因表达，从而防治真菌病害。
宏转录组揭示侵染真菌的病毒
1. 文库的构建如下图，使用 Illumina HiSeq2500进行测序。
2. 生信分析首先用trinity软件中k-mers 25和30拼装reads。接下来如下图：
3. qRT-PCR和RACE实验在分析了数据后进行，还用FirstChoice RLM-RACE kit和SMARTer RACE kit确定了tobamo-like和tombus-like病毒的两端，获得了完整序列。
4. 经过以上分析和实验，发现RNA-seq结果获得 0.7×10^8 到3.2×10^8个100nt的paired-end sequence reads。其中0.5%到15.1%显示是病毒源reads。C.truncatum (2 isolates)中2个病毒序列；M.phaseolina (48 isolates)中发现14个序列；D. longicolla中1个；R.solani (84 isolates)中27个；S.sclerotiorum (138 isolates)中28个，一共72个病毒序列，来自15个不同的科。depth为15到247,742-fold，平均20,687-fold，去除mitovirus的reads，平均为7,216-fold。
72个序列中，有五个与Benyviridae, Ophioviridae和Virgaviridae科亲缘关系近，而这三个科之前没有发现过真菌病毒。
检测到的病毒中，73%为ss（+）RNA病毒，15%为dsRNA病毒，12%为nsRNA病毒；49%为线粒体病毒。
针对72个病毒序列的分析和建树
1. Sequences related to members of the Mononegavirales
2. Fusariviridae- andHypoviridae-related sequences.
3. Endornaviridae-related sequences
4. Narnaviridae and Ourmiavirus-related sequences
5. Putative positive-sense single-stranded extramitochondrial RNA virus genomes
6. Putative double-stranded RNA virus genomes
文章中的概念解释
k-mer：是指将一条序列分成包含k个碱基的子字符串，如果reads长度为L，k-mer长度设为k，则产生的k-mers数目为：L-k+1，例如序列AACTGACT，设置k为3，则可以将其分割为AAC ACT CTG TGA GAC ACT共6个k-mers。其中k一定是奇数，如果是偶数遇到回文序列可能会产生完全相同的k-mers。
Mononegavirales：单分子负链RNA病毒目。
这里面很多系统发育分析，很基础很重要的分析内容。一个生物学家说过一句话，写在这篇推送最后吧！
“Nothing in Biology Makes Sense Except in the Light of Evolution.”
——Theodosius Dobzhansky
Marzano, S.-Y.L., Nelson, B.D., Ajayi-Oyetunde, O., Bradley, C.A., Hughes, T.J., Hartman, G.L. et al. (2016) Identification of diverse mycoviruses through metatranscriptomics characterization of the viromes of five major fungal plant pathogens. Journal of virology 90: 6846-6863.