大刍草苗期转录组RNA—Seq数据的denovo拼接

合集下载

生物大数据技术在转录组组装与注释中的使用技巧

生物大数据技术在转录组组装与注释中的使用技巧

生物大数据技术在转录组组装与注释中的使用技巧转录组组装和注释是生物学研究中的关键步骤之一。

随着高通量测序技术的发展,生物大数据的应用越来越广泛,为转录组组装和注释提供了更多的工具和方法。

本文将介绍生物大数据技术在转录组组装与注释中的使用技巧。

转录组组装是将测序数据拼接成完整的转录本,以便分析其表达模式和功能。

传统的转录组组装方法通常基于参考基因组,但对于非模式物种或缺乏参考基因组的物种,这种方法可能有限。

在这种情况下,使用生物大数据技术可以克服这个问题。

首先,对于非模式物种,可以使用转录组组装的无参考方法。

这些方法利用图形理论和算法将测序数据组装成转录本。

例如,Trinity、SOAPdenovo-Trans、IDBA-Tran等是常用的无参考转录组组装工具。

它们能够识别转录本的跨越剪切位点、剪切异构体和基因家族重复等复杂结构,从而提高组装的准确性。

另一种生物大数据技术是参考转录组组装。

这种方法是基于已有的参考转录本序列进行组装。

常用的参考转录组组装工具包括Cufflinks、StringTie、Traph等。

这些工具使用的是启发式算法,能够根据已有的参考序列评估和组装新的转录本。

因此,相对于无参考组装方法,参考转录组组装更加快速和准确。

除了转录组组装外,注释转录组也是生物学研究的重要一步。

转录组注释包括基因功能注释和转录本表达注释。

生物大数据技术可以为转录组注释提供丰富的信息资源。

对于基因功能注释,生物大数据可以提供已知基因的功能和相互作用网络。

例如,Gene Ontology (GO) 数据库提供了基因功能、过程和组件的分类信息。

另外,KEGG数据库可以帮助研究人员了解基因在代谢途径和信号转导网络中的作用。

通过结合这些数据库,可以对转录组的基因功能进行更详细的注释。

转录本表达注释是了解转录组表达模式的重要方法。

在生物大数据中,储存了大量的转录组表达数据,例如RNA-Seq和Chip-Seq数据。

转录组Denovo手册(无答案)

转录组Denovo手册(无答案)

对于初级分析的项目,只需要给合作伙伴提供过滤后的数据即可,所以会对过滤后的数 据做dt1h LReenagdt2h
N3Z5123 FUNzPTEARAA 503060871000;536871 5725.71;52. 99.9971;99. 75 75
确定的?我们所说的插入片段长度是指 括了 read1 和 read2 本身的长度?
read1

read2
之间没有测到的那一段的长度还是包
1112..什解么释是Soilnedxeax测测序序中,几进个行关in键de的x技测术序:的边主合要成目边的测是序什(么S?BS),可逆阻断技术和桥式 。 PCR
2.信息分析流程:
软件(Conesa, A., S. Gotz, et al. (2005). "Blast2GO: a universal tool for annotation, visualization 得到 的 and analysis in functional genomics research." Bioinformatics 21(18): 3674-6.) Unigene
Unigene
相似根性的据蛋KE白G,G从注而释得信到息该我U们n能ige进ne一的步蛋得白到功U能ni注ge释ne信的息P。athway 注释。
统计我,如们将下图Un所ig示en:e 和 COG 数据库进行比对,预测 Unigene 可能的功能并对其做功能分类
根据nr注释信息我们能得到GO功能注释。我们根据nr注释信息,使用 ( ) Blast2GO 2.3.5
本节问题: 1.Q20 是什么意思? 2.BMS 系统上给出的 Q20%值是如何计算出来的? 3.转录组暂时执行数据质量标准是怎样的?你有什么更好的建议(拿出自己的测试数据)? 4.在统计数据信息时,read1 和 read2 长度相等吗? 5.read每个碱基测序错误率的分布如何?read测序长度增加有什么好处?为什么SOAP比对 的时候允许 3’端有更多的错配? 6.如何根据 BMS 上的碱基频率分布图查找建库或测序失败的问题?

华大基因转录组结题报告(de novo)

华大基因转录组结题报告(de novo)
华大基因转录组分析(de novo)结题报告
2009 年 月 日
华大基因转录组分析(de novo)结题报告
目录
一、 项目信息................................................................ 2 二、 工作流程说明............................................................ 3
Total Length N50 Length Mean Length
59,067,262 1,968 1,260
表 4 补洞后 gap 长度占 scaffold 长度的百分比统计:
Gap Length Percentage (to scaffold) Gap number Gap percentage (to total gaps)
3.2.5 Scaffold-gene表达差异分析
说明:我们利用软件 soap 将不同样本中得到的 Reads 比对到 scaffold-gene 上,获得 scaffold 上 reads 的数目,然后计算 scaffold 在不同样本之间表达差异的 P value,一般
6
华大基因转录组分析(de novo)结题报告
2.1 实验流程说明........................................................... 3 2.2 信息分析流程说明....................................................... 3 三、 项目结果报告............................................................ 4 3.1 数据处理和质控报告..................................................... 4

HiSeq+PacBio denovo拼接项目方案

HiSeq+PacBio denovo拼接项目方案

一DNA de novo测序服务介绍De novo 测序也叫从头测序,不需要任何基因序列信息即可对某个物种进行测序。

用生物信息学的分析方法对序列进行拼接、组装,从而获得该物种的基因组序列图谱。

目前广泛应用于从头解析未知物种的基因组序列、基因组成、进化特点等。

全基因组序列图谱完成后,可以构建该物种的基因组数据库,为该物种的后基因组学研究搭建一个高效的平台;为后续的基因挖掘、功能验证提供DNA序列信息。

上海伯豪利用新一代高通量测序技术进行全基因组从头测序,并通过生物信息学的手段,帮您完成物种基因组序列图谱。

二、样本的抽提与质检2.1 DNA的抽提与纯化本实验采用DNeasy Blood & Tissue Kit (Qiagen, 69506) 并且根据生产厂商提供的标准操作流程进行样本的DNA抽提及纯化。

2.2 DNA质检说明纯化后的DNA 经Qubit®2.0 Fluorometer和0.8%琼脂糖凝胶电泳进行质检,质检合格的DNA可进行后续的测序实验。

2.3样本DNA要求1) 样品纯度:OD 260/280值应在1.8~2.0 之间;RNA 应该去除干净。

2) 样品浓度:最低浓度不低于50ng/µl。

3) 样品总量:每个测序样本一般需要纯化后至少2微克(ug),pacbio测序样本一般需要纯化后至少10微克(ug)。

4) 样品溶剂:要溶解在H20或TE (pH 8.0)中。

5) 样品运输:DNA低温运输(-20℃);且在运输过程中请用parafilm将管口密封好,以防出现污染1 / 152 / 15三、HiSeq 测序实验3.1 实验流程3.2.1 测序文库构建1) 文库构建:按照实验说明书对基因组DNA 进行片段化,末端修复、3’末端加A 、连接接头、富集等步骤,完成测序样本文库构建。

2) 文库质检:所建文库使用 Qubit® 2.0 Fluorometer 检测浓度,Agilent2100检测文库的大小。

利用玉米高通量RNA-Seq数据预测长链非编码RNA(LncRNA)

利用玉米高通量RNA-Seq数据预测长链非编码RNA(LncRNA)

Ab s t r a c t : Nu me r o u s s t u d i e s h a v e s h o w n t h a t mo s t o f t h e p l a n t g e n o me i s t r a n s c i r b e d t o y i e l d n o t o n l y p r o t e i n — c o d i n g
a g r o u p c o n s i s t i n g o f n c RNA l o n g e r t h a n 2 0 0 n u c l e o t i d e s ,h a v e b e e n ou f n d t o p l a y a c i r t i c a l r o l e i n g e n e r e g u l a t i o n .T h u s , t h e d e s c i r p t i o n o f t h e L n c R NA i n t h e p l a n t g e n o me w i l l b e b e n e i f c i l a f o r t h e u n d e r s t a n d i n g f o g e n e e x p r e s s i o n v a ia r t i o n i n
Pr e d i c t i o n a n d c h a r a c t e r i z a t i o n o f l o ng n o n— c o d i n g RNAs u s i ng h i g h t h r o u g h pu t RNA- S e q d a t a i n ma i z e
吕远 大 , 李 坦 , 张晓林 , 赵 涵

动植物Denovo测序知识大讲解

动植物Denovo测序知识大讲解

动植物Denovo测序知识⼤讲解⾼通量测序的技术开起我们探索动植物基因组奥秘的步伐,提到动植物基因组测序,这就不得不提⼀个概念——de novo测序。

那么什么是de nove测序呢,它与重测序有什么区别呢?De nove测序中Read、Contig和Scaffold等⼜代表什么呢?De nove测序中为什么要建不同⼤⼩⽚段的梯度⽂库?基因注释⼜是注释哪些内容?各位客官别急,且听⼩编给您细细讲来。

1De novo测序概念De novo是⼀个拉丁⽂,代表从头开始的意思,⽽de nove测序则是指在不需要任何参考序列的情况下对某⼀物种进⾏基因组测序,然后将测得的序列进⾏拼接、组装,从⽽绘制该物种的全基因组序列图谱。

由于⾼通量测序长度的限制,⽬前测序策略是先将基因组打断⼩的⽚段,然后再对测出序列⽚段进⾏拼接,最终得到物种的序列图谱如图1所⽰。

图1 ⾼通量测序模式图2De novo测序与重测序区别重测序概念:重测序是全基因组重新测序的简称,是指是对已知基因组序列的物种进⾏不同个体的基因组测序,并在此基础上对个体或群体进⾏差异性分析。

从概念上来看两者的区别在于de nove测序是对没有参考基因组的物种进⾏测序,⽽重测序是对已有基因组的物种进⾏测序,这只是它们区别很⼩的⼀部分。

从原理上来看de nove测序和重测序最根本的区别在于de nove测序需要对测序得到的Reads进⾏拼接组装,⽽重测序得到的数据则是没有组装的短的Reads序列。

值得注意的是,随着测序成本的降低以及组装算法的改进,de nove测序成本越来越低,⽬前来说de nove测序不只对于没有参考基因组物种进⾏测序,还可以对⼀些特有的亚种、品种以及变种等进⾏测序。

3Reads Conting Scaffold概念Reads:即我们通常说的读长的意思,它是指⾼通量测序平台直接产⽣的DNA序列。

Contig:是指Reads基于Overlap关系,拼接获得的长的序列;Scaffold:是指将获得的Contig根据⼤⽚段⽂库的Pair-end关系,将Contig进⼀步组装成更长的序列;关于三者之间的关系如图2所⽰,注意的是Contig是⽆Gap的连续的DNA序列,⽽Scaffold是存在Gap的DNA序列。

Dual RNA-seq技术及其在宿主-病原体相互作用中的研究进展

Dual RNA-seq技术及其在宿主-病原体相互作用中的研究进展邹宏;夏应菊;徐璐;赵俊杰;李玲;王震;刘业兵;王琴;宋振辉;张乾义【期刊名称】《中国兽药杂志》【年(卷),期】2022(56)9【摘要】利用转录组测序技术(RNA-sequencing,RNA-seq)进行转录组分析是了解病原体入侵宿主分子变化的重要工具,在同时分析病原体与宿主转录组时,RNA-seq技术需要分别构建病原体及宿主的cDNA文库,再将其各自映射到病原体及宿主参考基因组中,而互作转录组测序技术(Dual RNA-seq)无需分离两物种,只需构建一个转录组文库,便能同时对两个(或多个)研究对象进行测序和分析,可以直观地揭示病原体和宿主相互作用过程中转录组学动态变化,因此Dual RNA-seq技术被广泛应用到人类疾病和生物感染模型的相互作用研究中。

为了解Dual RNA-seq技术及其在宿主-病原体相互作用研究中的前景,就Dual RNA-seq技术概述以及近年来该技术在原核生物、真核生物以及病毒研究中的应用现状及发展前景进行综述。

Dual RNA-seq技术可为病原体与宿主相互作用的研究提供新视角,有助于更好地识别和理解感染过程中病原体和宿主的转录组学变化,从而揭示潜在的新靶点或生物标记物。

【总页数】8页(P71-78)【作者】邹宏;夏应菊;徐璐;赵俊杰;李玲;王震;刘业兵;王琴;宋振辉;张乾义【作者单位】西南大学动物医学院;中国兽医药品监察所【正文语种】中文【中图分类】S855【相关文献】1.泛素系统在结核分枝杆菌与宿主相互作用中的调控机制研究进展2.酵母双杂交系统在痘病毒与宿主相互作用中的研究进展3.伯氏疏螺旋体传播过程中病原体、媒介、宿主间的相互作用(英文)4.转录组测序——RNA-Seq技术在动物疾病检测中的应用研究进展5.坏死性凋亡在细菌与宿主相互作用中的机制研究进展因版权原因,仅展示原文概要,查看原文内容请购买。

大刍草苗期转录组RNA-Seq数据的de novo拼接

大刍草苗期转录组RNA-Seq数据的de novo拼接肖之夏;郑用琏【期刊名称】《华中农业大学学报》【年(卷),期】2014(33)2【摘要】采用solexa测序技术,对大刍草苗期全株各组织转录组进行了RNA-Seq 测序、de novo拼接和信息比对研究。

结果表明:转录组测序共得到了46.4 GB的原始数据,归并整理后获得长76 bp的序列有175 101 250条,经质量控制和de novo拼接后,共获得了58 147条大刍草转录本,其平均长度为1 335 bp。

比对分析发现其中94.3%的转录本和玉米B73自交系的cDNA序列有较好的匹配,与水稻匹配的有84.1%,高粱84.6%,短柄草83.9%,共56 036条转录本。

【总页数】7页(P15-21)【关键词】大刍草;转录组;RNA-Seq;de;novo拼接【作者】肖之夏;郑用琏【作者单位】华中农业大学作物遗传改良国家重点实验室【正文语种】中文【中图分类】S519.24【相关文献】1.基于转录组de novo拼接鉴定辣椒(Capsicum fruteseens)中辣椒素生物合成相关基因 [J], Liu Shaoqun;Li Wanshun;Wu Yimin;Chen Changming;Lei Jianjun;高凌云;孔秋生2.貉外周血单核细胞转录组RNA-Seq数据的de novo拼接和信息比对研究 [J], 易立;仝明薇;程悦宁;程世鹏3.基于RNA-Seq数据集的转录组从头拼接算法 [J], 武思文;李静;张少强4.基因组De novo组装结合比较RNA-Seq策略在VBNC状态乳杆菌转录组学研究中的应用 [J], 徐振波;侯玉超;刘君彦;邓阳;张霞;李冰;李琳5.盐胁迫下秋茄根系转录组的De novo测序及分析 [J], 邢建宏;陈伟因版权原因,仅展示原文概要,查看原文内容请购买。

转录组测序在牧草抗逆基因挖掘方面的应用前景

10.16863/ki.1003-6377.2018.03.002转录组测序在牧草抗逆基因挖掘方面的应用前景朱昊12,张博\张荟荟2袁王玉祥\柯梅2袁李学森2*(1.新疆农业大学草业与环境科学学院,乌鲁木齐830052;2.新疆畜牧科学院草业研究所,乌鲁木齐830011)摘要:本文系统阐述了转录组测序技术(RNA-seq)的原理和技术流程。

查阅近五年国内外的文献,梳理 了植物在干旱、盐胁迫及温度胁迫等逆境下的转录组表达方面的研究情况,发现有关牧草方面的转录组 研究相对较少且不够深入。

RNA-seq在牧草优良种质资源抗逆基因挖掘方面的应用前景十分广阔。

关键词:转录组测序;牧草;基因挖掘;应用前景中图分类号:S540.34 文献标识码:A文章编号:1003-6377(2018)03-0009-04我国的牧草种质资源极为丰富,目前的研究主要集中在以全球变化与农业多样性、农业系统固碳减 排、草地农业决策系统及林草生态等宏观方面[1]曰同时以分子生物学为基础的微观方面也有许多研究,主 要包括DNA水平上的遗传多样性和分子标记辅助选择育种等方面[2-3]。

利用新一代高通量测序技术进行 牧草种质资源创新及抗逆基因挖掘方面的研究相对较少[4]遥1RNA-se q技术原理RNA-seq即转录组测序技术。

克里克(FrancisHanyComptonCrick,1916-2004)于1958年提出“中心法则”,并于1970年在Nature上的一篇文章中进行了重申。

中心法则指出遗传信息从DNA传递给 RNA,再从RNA传递给蛋白质,即完成遗传信息的转录和翻译的过程。

依据中心法则,随着二代测序技术 的兴起和成熟,转录组学研究得到了蓬勃发展。

简而言之,转录组学(transcriptomics)就是从RNA水平研 究基因表达的情况,是一门在整体水平上研究细胞中基因转录的情况及转录调控规律的学科。

转录组(transcriptome)广义上指某一生理条件下,细胞内所有转录产物的集合,包括mRNA、rRNA、tRNA及 Non-codingRNA;狭义上指所有 mRNA 的集合。

高通量转录组测序的数据分析与基因发掘

高通量转录组测序的数据分析与基因发掘周华;张新;刘腾云;余发新【摘要】高通量转录组测序(RNA-seq)是在转录组水平上进行深度测序的一项技术,为真核生物转录组学的研究开创了新平台,但同时测序所得到的海量数据的生物信息学分析成为科研工作者的一大挑战。

对转录组测序技术进行了阐述,重点介绍了转录组测序后的数据分析,以及在真核生物尤其是非模式物种中的基因发掘方法。

%High-throughput transcriptome sequencing (RNA-seq) is a recently developed approach to transcriptorne profiling that uses deep-sequencing technologies. It provided a novel platform for eu- karyotie transeripome researches, but the bioinformatics analysis of sequencing data became the chal- lenge of scientific worker. In this review, we described the researches process of high-throughput transcriptome sequencing technology, focusing on sequencing data analysis and gene discovery of dif- ferent species, especially of non-model species.【期刊名称】《江西科学》【年(卷),期】2012(030)005【总页数】5页(P607-611)【关键词】转录组测序;数据分析;基因发掘【作者】周华;张新;刘腾云;余发新【作者单位】江西省科学院生物资源研究所,江西省观赏植物遗传改良重点实验室,江西南昌330029;南京市林业站,江苏南京210036;江西省科学院生物资源研究所,江西省观赏植物遗传改良重点实验室,江西南昌330029;江西省科学院生物资源研究所,江西省观赏植物遗传改良重点实验室,江西南昌330029【正文语种】中文【中图分类】Q987转录组研究是一个发掘功能基因的重要途径,是基因功能及结构研究的基础和出发点。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

玉米 属 于 禾本 科 ( P o a c e a e ) 玉蜀黍属 ( Z e a ) , 玉 义上 包 括 了 mRNA、 n o n - c o d i n g RNA 等 ) , 进 而 推
蜀黍 属 又 被 划 分 成 5个 种 , 分别是 Z . d i p l o p e r e n — 断完整的基因结构, 确定选择性剪切事件 , 研究在不
比对研 究 。结 果 表 明 : 转 录组 测 序 共 得 到 了 4 6 . 4 G B的原始数据, 归并 整理后 获得长 7 6 b p的序 列 有
1 7 5 1 0 1 2 5 0 条, 经 质量控制 和 d e 舢l U o拼接后 , 共获得 了 5 8 1 4 7条大刍 草转 录本 , 其平均长 度为 1 3 3 5 b p 。比对 分 析发现其中 9 4 . 3 %的转录本和玉米 B 7 3自交系 的 c D NA序 列有较好 的匹配 , 与水 稻 匹配 的有 8 4 . 1 , 高粱
Z . ma y s s s p . me x i c a n a和 Z. ma y s s s p . ma y s聚在 本 ; 而 基 于 高 通 量 测 序 的 RNA— S e q技 术 则 可 以较

类, 后续 对 叶绿体 、 核糖 体 的研 究也 得 到 了相 似结 为全 面地 、 对 几乎 全部 的 R NA 转 录本 进行分 析 。 根据 所 研 究 的 物 种 是 否 有 参 考 基 因 组 信 息 , 的序 列拼 接 和 d e n O V O序 列 拼接 , 或 当基 因组 信 息
n i s 、 Z . p e r e n n i s 、 Z . 1 u x u r i a n s 、Z . n i c a r a g u e n s i s 和 同组 织 、 不 同发育 阶段 、 不 同实验 处 理 中的相关 基 因
Z . ma y s , 其中 Z . ma y s有 4个 亚 种 , 分别是 Z . 在转 录水 平上 的表 达 水 平 。近 年 来 , 随着 转 录 组学
a y s , 前 3 个亚种被统称为大刍草 , 第 4个 是 经 人 法” 主要是 利 用高 密度 的商 业化 芯 片进行 研究 , 但 限
类长 期栽 培 驯化 的玉米 。果 穗化 石证 据 的缺失 以及 制这 一方 法广 泛应 用 的瓶 颈是 必须 预先得 到试 材对
玉米 和大 刍草果 穗 上 的 巨大 差异 , 一 度 让科 学 家无 象 的全基 因组序列 信息 , 检 测杂 交信 号误差 也较 法确 定 玉 米 的 直 系 祖 先 。L o n g l e y [ 1 ] 通过 对 Z . 大[ 8 ] ; “ 测序法” 包 括 基 于 标 签 序 列 代 表 基 因 的
大 刍 草 苗期 转 录 组 RN A — S e q数 据 的 d e n o v o拼接
肖之夏 郑用琏
华 中农业 大学作物遗传 改 良国家重点 实验 室, 武汉 4 3 0 0 7 0
摘 要 采用 s o l e x a 测序技术 , 对大刍草苗期全 株各组 织转 录组进 行 了 R NA - S e q 测序、 d e ; ' 1 0 " 0 0拼 接和信 息
8 4 . 6 , 短柄草 8 3 . 9 , 共5 6 0 3 6条转录本 。 关键词 大刍草 ; 转 录组 ; R NA- S e q ;d e ? I O W O拼接 S 5 1 9 . 2 4 文献标识码 A 文章 编号 1 0 0 0 - 2 4 2 1 ( 2 0 1 4 ) 0 2 — 0 0 1 5 — 0 7 中图分类号
me x i c a n a可 能 就是 玉 米 的祖 先 , Ka t o等 _ 2 通 过 对 录组 的部 分序 列进 行 分 析 , 无 法 得 出完 整 的基 因结
染色体 形 态 的研 究 则将 Z . ma y s s s p . p a r v i g l u mi s 、 构信 息 , 也无 法分 辨 出选 择 性剪 切 产 生 的不 同转 录
果_ 4 ] , 而 D o e b l e y等 ¨ 6 ] 基 于 同 工 酶 以 及 Ma t s u o k a ma y s s s p . p a r v i g l u mi s 是玉 米 的祖先 。
等 ] 基于 S S R 标记 的遗 传 多 态 性 研 究 , 则认 为 Z e a RNA- S e q序 列拼 接 的策 略 也 分 为 基 于 参 考 基 因组 转 录组 是 一个 细 胞 中的所 有 转 录 本信 息 , 包 括 不完整时 , 将二者结合起来进行序列拼接。基于参 转 录本 的数 量 、 特定 发育 阶段 的表 达动 态 、 转 录后 的 考 基 因组 的拼 接方 法是 将测 序得 到 的序列 匹配 到基 修 饰 以及 n o n — c o d i n g RNA 的调 控 表达 情 况 等 。转 因组 上去 , 根据 r e a d s的 重 叠 情 况 和 是 否 跨 越 剪 切
ma y s s s p . me x i c a n a染 色体 臂 的长 度 、 中心 粒 的 位 S AGE 、 C AGE、 MP S S 以 及 基 于 高 通 量 测 序 的
置、 结 节 大 小 和 位 置 的研 究 指 出, Z . ma y s s s p . R NA— S e q技 术[ 9 ] 。杂交 法 和 标 签 序 列 法 都 只对 转
m a y s s s p .m e o c i c a n a、 Z. ma y s s s p .p a r v i gl u m i s 、
研究 的迅 速展 开 , 在 研 究 方 法 上也 有 较 多 的发 展 与
Z . ma y s s s p . “ P
P 盘 g 5 s和 Z. ma y s s s p . 创新 , 主要 包 括基 于杂 交 和基于 测序 的方 法 。“ 杂交
第 3 3卷 第 2期
2 0 1 4 年




Vo L 3 3 No . 2
Ma r .2 O1 4, 1 5~ 21
J o u r n a l o f Hu a z h o n g Ag r i c u l t u r a l Un i v e r s i t y
相关文档
最新文档