免费文献传递   相关文献

Study on quality evaluation of sequence and SSR information in transcriptome of Astragalus membranacus

膜荚黄芪的转录组测序质量评估及其SSR位点信息分析的研究



全 文 :2016
y

z Ó
41
âÓ

}
Vol41,No.8 April,2016

st)(
] 20151110

†‡ˆR
] 
JÖ×fñŸdˆ¥
(81274010);
*H’ºY+j̟dˆ¥
(JC201101);
*H’Sehƒñ

,ԄE-å&.£¤

ˆ¥
(2012001);
/01o,Ôñf23-Ÿdˆ¥
(2014RFXXJ122);
*H’ºä´Èfñ†ijˆ¥
(12541743)

uq?v
] 
¢-

ij‰
,Tel:(0451)82193430,Email:mawei@hljucmnet;


ij‰
,Tel:(0451)82196224,Email:1447111706@
qqcom

?v‹Œ
] 
™š

Õij¦

Š‹ŒSh…4©Ehij
,Email:459820677@qqcom
ïðñòÄâã´©óôõö÷­®
SSR
øùvíµîą†
™š

›œ

$ž

Ÿ- 


,
¢-


*H’Sehƒñ hñg

*H’ /01
150040)

wx
] 
ã½.“°
Roche454GSFLX
©óÚâñòÄâã´=¼

ü°
454SequencingSystemSoftware
µî†p‘nâ
ã´q2T«

Ü°
MISA
g:=€óñòâã´©óâÄ
9893
Ï
unigenes,
¶®
SSR
øùví‘nóµî

ä+

‘n©óVÄ
reads
Ä͑(ré
413bp,
š
86%
Ä
reads
ÐÔóT«

T«Ä
N50
(ré
1205bp,
V©Ä
unigene

×ugõów1âã´ví

ñòâã´pbT
1729
Y
SSR
øù
,SSR
Ämː˜é
924%,SSR
¦ñòÌYâã´$å×
Иé
1342%,SSR
Ä͑cªé
797kb。
[amכS;uó¯
127
,

ۗ¶Äæ؛™çq$Ä
TG/AC
q

å×
ИÛu
SSR
øùÄ
425%。
ñòâã´Ä©óäÙHóñòâã´ÄÌm+XÆÇ

ñToõñòâã´
unigene
ó¯

ñSñòâã´
SSR
øùåא˜O

GqØ&

Ø¿¹ãöO


yz{
] 
ïðñò

âã´©ó
;SSR;
øùví
StudyonqualityevaluationofsequenceandSSRinformationin
transcriptomeofAstragalusmembranacus
CHANGYue,YANSong,LIUZhenpeng,RENWeichao,LIUJie,MAWei
(ColegeofPharmaceuticalSciences,HeilongjiangUniversityofChineseMedicine,Harbin150040,China)
[Abstract] Inthisstudy,454/RocheGSFLXsequencingtechnologywasusedtoobtainthedataoftheAstragalusmembranaceus
FourhundredandfiftyfourSequencingSystemSoftwarewasappliedtocaryoutthetranscriptionofthegroupfromscratchUsingMI
SAtools,9893unigeneswereselectedforthesequenceofthegenomeofAmembranaceus,andtheinformationofSSRlocuswasana
lyzedAccordingtotheresult,theaveragelengthofreadswas413bp,about86% ofthereadswasinvolvedinthesplicing,thelength
oftheN50was1205bp,thenumberofunigeneswasmeasuredbythewholetranscript1729SSRlociintheAmembranaceustran
scriptomeweresearched,theoccurencefrequencyofSSRwas924%,thefrequencyofSSRinthewholetranscriptomewas1342%,
theaveragelengthofSSRwas797kbOnehundredandtwentysevenkindsofcorerepeatsequenceswerefound,thedominanttype
wasTG/ACtypeofdinucleotide,itappearedtoaccountfor425% ofthetotalSSRlocusTheresultsofthesequenceofthetranscrip
tionoftheAmembranaceustranscriptomerevealedtheoveralexpression,andalargenumberofunigenessequencewasobtained,and
theSSRlocusinthegenomeoftheAmembranaceusishigh,andthetypeisdiverse,andthepolymorphismofthegeneishigh
[Keywords] Astragalusmembranaceus;transcriptomesequencing;SSR;informationofloci
doi:10.4268/cjcmm20160810
·0341·
© 


ïðñòÄâã´©óôõö÷­®
SSR
øùvíµîą†
  
$,Xñòˆ

$%,.
》(2010


!

æ— ƒ ñ ò
Astragalusmembranaceus(Fisch)
BgevarmongholicusHsiao Bge
 ï ð ñ ò
Amembranaceus(Fisch)Bge
ċŒ¥


‰

:B>ïp

jÜ4

ðºðû

–—˦Ö
ʬ
[1]。
âã´
(transcriptome)
æ`Æ!56¦Y
[Êö¾¿Àw1+XÄ×Øu»

P+ó|[Y
×ØÄXJ»+X‚Í

âã´©óöw9ÄFÙ
HËÉYm¦Æ!´7»Æ!O}Äwo×ØÄ+
XáÎ

úî;uó¯
(simplesequencerepeats,
SSR),
[—ç
1~6
YŸ×é›Só¯

:BOrØ
¿¹

øùµÜ°
SSR
õýbódØÉ,Än
Rm@Wq¥

ñZúûê°Þ×Ø!ø­A˜©
Ùµî

[,d!»&ÈÉÓ,֌`
[2]。
u½.
™¦sÊñòâã´©óñµî®
SSR
øùví

éñò‘[¸¦@Wq¥ýb

,ôd!

@WØ&
¹µîÖ9ą†D!ó×ë

1 
XLÔ#
11 
ÈÉ
uB.V°ïðñò,-

UŒ¡p$A,o
B.$Sǎ…†Gd!

¦Œ¡p$A,o
,°ÈÉI{¤‰þ8fË>??É

€ <š
40

Äñò

 X¥Ä¥

4

é

“Žž°PëM

-80℃
B

12 
ÈÉu
RNA
ğ 
ñòu
RNA
Ÿ ÐÜ0Ä
CTAB
#
[3],
:
m¸1IÀ

 õÈÉXLޅ›$‘iPë…
œ“‡¢¾

-‡¢âê“
20mL
ªS*$

‘i
1mL
f‚¯Œ

?µò͞

¶rƒž
10min,4℃
12000r·min-1
ªS
2~5min;
:Sörh

¦ª
S*$‘i
900μLCTABŸ P,òÍO=ž,65

‚ã
5min,
ƒž¶r

‘iÖmËȞ
,4℃
12000r·min-1
ªS
5~10min;
:S“ rh“ß
ĪS*$

‘i
1/2
mËf‚¯Œ­
08


mol·L-1LiCl,
O=

¶rƒž
10min,4℃ 12000r
·min-1
ªS
10min,
rh
;70%
¯ŒwߍI

Š

þ‰ƒž
5min
‹Œ
RNA,
‘i
50μLË­‚N
œ
RNA,-80℃
ì°

13 mRNA
&[‚ì
131 
u
RNA
Äôõ¨©Ô!õ
 
ü°á£µ¤
¤rD

¶u
RNA
ēr

uõ­
A260/280» A260/230
‘nµî

“°
AgilentRNA6000PicoKit;Agilent,
50671513,
¶u
RNA
Äs̹

7r­›œr‘n
©!

132 mRNA
Ä7_Ô!õ
 
“°
OligotexmRNA
MiniKit;Qiagen,70022;AgilentRNA6000PicoKit;
Agilent,50671513
B„¼

‘n
mRNA
7_Ô!õ

14 cDNA
^V±ì
mRNA
€Ü_“ªC¦
450~1200nt;
ç
1∶1∶1
OH¥

4

éu
RNA;
“°
DNASynthesissystem,
Roche,11117831001
»
AgencourtAMPureXPBeads,
Beckman,A63881
B„¼

‘nÝ_
cDNA
HÊÔ7
_

“°
RocheRapidLibraryPreparationKit,Roche,
5608228001
B„¼

‘n€Ü¢Vtu

ˆ««2

“°
AgencourtAMPureXPBeads,Beckman,A63881
B„¼

ù­:€Ü

ü°ÚÃ
QuantifluorSTfluo
rometer,Promega,E6090,
‘n
cDNA
^V!õ

ü°
ÚÃ
Agilent2100Bioanalyzer,Agilent,2100
»
Agi
lentHighSensitivityDNAKit,Agilent,50674626
B
„¼

‘n
cDNA
^Vôõ¨©

‚ìg¯P

15 emPCR
Äc
@•ŸP !©ó !$!
emPCR
Äc
$VXÄ
DNA
^VÄõ

‚ì
emPCR
xB„

Ÿ

、MockAmplificationMix
»ÑŸP
、LiveAm
plificationMix;
ü°
DNACaptureBeads
-
DNA
^V¡â

Ÿ_

Äc

ɔ
DNACaptureBeads;
x
DNA
^VÄ
DNACaptureBeads
įŽ

©ó
=ÉhÔ

16 Roche454GSFLX+
©óÚrf©ó
âã´©óg¯F‘rH¢ÊºËÉ&(BC
DEsÊ

17 
2=¼Ìè

••­ôõö÷
Àf=¼’•
454Newbler
ù«2

ù±ôõŸ
×±è

ü°
454SequencingSystemSoftware(ht
tp://www454com/)
µî†p‘nâã´q2

[45]。
18 
©ó£»rµî
©ó£»ræ12“&à Ä©óõ
(sampled
readsnumber)
Äc‘

í¨©TÄ
unigene
=øā
_áÎ

“&à Ä©óõT[!–r

í¨©

unigene
ËõÒc‘Ãwkc‘

ô©óý
»

¶Ã°æ©óõÒ÷NBXT£»

¦uŠ©
ó=¼$

ç
15000
¯éŠr

¶=¼×Š1f€
&

S€å(Ä78
reads
µe¨©TØk
unigene。
·1341·
2016
y

z Ó
41
âÓ

}
Vol41,No.8 April,2016
nž|“&à Ä©óõMÔ¤»¨©TÄ
uni
gene
Ä=õM¥¤¦[Yâ§

/S7쉤1
2“&à =¼õÄc‘
unigene
æ¶:B
£»¹

19 
âã´
SSR
øùÄvíµî
-âã´=¼°
MISAMIcroSAteliteidentifica
tiontool(htp://pgrcipkgaterslebende/misa/)
‘
n
SSR
µî

–óÿžé
:1/10,2/6,3/5,4/5,5/5,
6/5(
ST;uěSó¯

]:;uîø=ø
);
u
Hq
SSR
$JƒST›Só¯Ä½;uŸ×=]
oÒN•
100[68]。
2 
äÔµî
21 
u
RNA
ôõµî
ôõ“r≥250mg·L-1,uõ≥50μg;A260/280
18~22,A260/230ê≥20。UN¨©28S∶18S“k

15,
lq
1;RIN≥80;ñ$RNAf›œ,f
nn

l+
1。
1,2
(>Þñò¥
;3,4
(>Þñò4
;5,6
(>Þñòé

q
1 RNA
 ~UN¨©
Fig1 ElectrophoresisimageoftotalRNA
+
1 
ñòu
RNA
Äᣨ©
Table1 SpectrophotometryoftotalRNA
No
(>
A260/280 A260/230 ôõ“r/mg·L-1 uõ/μg

¥
2068 2331 705 282

¥
2079 2353 655 262

4
2097 2351 606 2424

4
2071 2325 698 2792

é
2074 2337 959 3836

é
2087 2365 653 2612
  
g

Ÿ {}é
2012
y

z
27
{

‚õ
40μL。
qçr=¼ÕçSå

u
RNA
ôõHÐ

Õç
‘[¸‘nâã´©óg¯

22 
©óõ=D
©2=¼õš
513Mb,
@•••TO
HT««vÄB¬=¼
306Mb。
=¼u?ž

=
¼õé
306805437,
×ØÏ=é
742721
Ï

Í
‘(ré
413kb;
×Ø=Dž

×Øé
742721
Ï

à(oÞ
20bp
ó¯=é
742602,
2=¼
é
513365449,
••ž=¼é
306805437,

Ü°˜é
5976%。
’•âã´q2T«
,8604%
Ä
reads
ÐÔT
«

aD
639061
Ï
reads。
ó¯T«Ä
N50
(ré
1205bp,
(r
>N50
é
3823bp,
͑(ré
10695bp。
T«T
contig15167
Y

78
contig
íT«Êé
12851
Y
isotig。
f#‘[¸T«Ê
é
isotig
Ä
contig
Ô78
isotig

,´
Ê
10742
Y
isogroup,
uaT
12880
Ï
Isotig(
ýþf#‘[¸
T«Êé
isotig
Ä
contig),Unigene
=õé
9893,
Unigene
͑(ré
109712898bp。
23 
©ó£»rµî
“&à ©óõXT
2×102kb
O
,unigene
Ä=øµµÞ£»

“&à =¼õXT
3×102
kbunigene
Ä=¼µ’sw£»

7e¶ÞuŠ
ñò¥

4

éÄ֕»O&©ó

©óõµ’XT
£»

©óõ¡¢B.«v

lq
2。
q
2 
©ó£»rµî
Fig2 Sequencingsaturationanalysis
24 SSR
øùÄ=õÔµ@
¶ñòâã´Ä
12880
Ï
isotig
‘n
SSR
Äp
b

T
SSR
ÄuŸ×=æ
23834bp,
aTT
SSR
øù
1729
Y

uHq
SSR
øù
435
Y

xB
SSR
ø
ùÄ
isotig
Ï=é
1190
Ï

ýx[Yçr
SSR
øù
Ä
isotig
ÄÏ=æ
252
Ï
。SSR
Ämː˜

xB
SSR
Ä
unigene
=øÔu
unigene
=øĕ^

924%;SSR
¦ñòÌYâã´$å×И
(SSR
øùY=»u
unigene
=øĕ^
)1342%;
ñò
âã´$
SSR
Ä͑cª

u
unigene
Ä(rÔ
SSR
=øĕ^
)797kb,
l+
2。
·2341·
© 


ïðñòÄâã´©óôõö÷­®
SSR
øùvíµîą†
+
2 SSR
¦ñòâã´$å×И
Table2 OccurenceSSRsinAstragalusmembranaceustranscriptome
;uGq =ø xGq•»
/%
˜
/%
͑cª
/kb
u(r
/bp
͑(r
/bp
ç
653 3777 507 2109 1919 12
؛™ç
282 1631 219 4885 4472 18
L›™ç
730 4222 567 1887 11433 18
–›™ç
32 185 027 39357 732 24
/›™ç
14 081 009 114790 300 27
ԛ™ç
18 104 013 81028 486 32
  
q+

ÕçSåñòâã´
SSR
,G§¯

q
çTԛ™çÄx,›ç;uGqö÷S
T

fxÿQÚĕ¯–o
。SSR
Gq؎$¦î
›™ç

؛™ç

L›™çr

Û
963%,
®xG
qQ¶ÄÛBõw±

ñòâã´
SSR
øùó¯u(rXT
23834
bp,
®$L›™ç›S;uó¯Äu(r](

é
11433bp;
®Šæ؛™ç

é
4472bp;
«À(æ
ç

æ
1919bp。SSR
øùÄ͑(ræ
24
bp,
xGq
SSR
øùÄ͑(rµeæ
12,18,18,
24,27,32bp。
25 SSR
Äƹ
¦ñòâã´Ä
1729
Y
SSR
øù$amכ
S;uó¯
127
,

®$çq

,

؛™ç
q
11
,

L›™çq
57
,

–›™çq
28
,


™çq
11
,

ԛ™çq
16
,

7

,;uó¯
GqÄ;uŠ=o؎$¦
5~10
Š

XT
1173
Š

ÛÌY;uŠ=Ä
6784%;
®Šæ
10~15
Š

;uŠ=æ
356
Š

ÛÌY;uŠ=Ä
2059;
ÓL
Äæ
16~20
Š

;uŠ=æ
149
Š

ÛÌY;uŠ

862%。
A°æe

›Só¯;uŠ=¦
5~
20
ŠÄ
SSR
Äå׊=é
1678
Š

ÛÌY;uŠ

9705%,
›Só¯;uŠ=N•
20
ŠÄ
SSR
å×Ċ=•»oÒT
3%,
lq
3。
qå×Ä:m
Gqµî

ۗ¶Äæ؛™çq$Ä
TG/AC
q

å×ИÛu
SSR
øùÄ
425%;
®Šæ
GA/CT
q

å×ИÛu
SSR
øùÄ
394%;
«À(æ
GAA/CTT
q

˜æ
356%;CCT/GGA
q

˜æ
309%,
l+
3。
26 SSR
ÄÕ°¹öÏ
š»
SSR
Õ°¹Ä;«×¼æ
SSR
µ-õ
ÄØ¿¹

[— é
SSR
Ä(ræÂîؿ¹O
±Ä;«ØF


SSR≥20bpO,Ø¿¹–O;
   
q
3 SSR
›Só¯;uŠ=µ@
Fig3 Distributionoftherepeatsnumber
+
3 
ñò×~GqÄ;u˜
Table3 Distributionofrepeatmotifaccordingtorepetitionfre
quencyforAstragalusmembranaceus
;uGq ;u˜
/%
؛™ç
TG/AC
q
425
؛™ç
GA/CT
q
394
L›™ç
GAA/CTT
q
356
L›™ç
CCT/GGA
q
309

12bp≤SSR≤20bpO,Ø¿¹$Ö; SSR≤12
bp
O

Ø¿¹;±

ñòâã´=¼øH
SSR
Ä(
r؎$¦
12~45bp,
®$
12bp≤SSR≤20bp:
B$rØ¿¹Ä
SSR
aB
861
Y

Ûu
SSR
Ä
4980%;SSR≥20bp:BOØ¿¹Ä SSRa 422
Y

Ûu
SSR
Ä
2441%。
UOB…†+

O/×
~
SSR
ÄØ¿¹•±/×~ÄØ¿¹Ü¨ö±

¦ñòâã´=¼=Dµî$m×

(roÞ
20bp
Ä
SSR
oØÇÞ±/×~

aýx
SSR
øù
282
Y

ÕçÑD71µØ¿¹ãöOÄ
SSR
¦ñò
×ØÊö…†r:B–OÄÜ°Ï^

3 
ºÒÔäÒ
31 
ñòâã´©óôõö÷
ñòâã´©óV°Ä¥

4

韠Äu
RNA
¨©äHÐ

‘n©óVÄ
reads
Ä͑
·3341·
2016
y

z Ó
41
âÓ

}
Vol41,No.8 April,2016
(ré
413bp,
š
86%
Ä
reads
ÐÔóT«

T«Ä
N50
(ré
1205bp。
©ó£»rä+


Ä
unigene
=õ×ugõów1âã´ví

q
çräÕàuâã´Ä©óäX½Õv

Õç
¶=¼‘n‘[¸µî

32 
ñòâã´
SSR
øùvíµî
SSR
úûµ@Þx,X›ËÉÄ×Ø´$

o
š|ƒ
10~50kb
°B¦

Y
SSR。
¦ÈÉ$

͑
233kb
°B

Y
SSR;
Ý-éÈÉ$Ä
SSR
=õo
Þî-éÈÉ

ùÃ

Y
SSR
›JÄ͑Jcé
212kb,
žÃé
646kb。
@•¶ñòâã´=¼
Ä
SSR
µî

ñòâã´$
SSR
Ä͑cªæ
797
kb,
qçGqTԛ™çGq‘:ì

›S;
uó¯
127
,

7+ñò×Ø´8:B–O§r
Ä
SSR。
o1µÈÉÄ
SSR
;u×~ª«çØ

L›™
çqéª

fÒUÉ,›JĪ{
SSR
;uGqB
VÚÛ

u…†m×ïðñòâã´
SSR
;uG
qª«çL›™çéª

Ûw1
SSR
Ä
4222%,
[
›™çVە»A–O

Ûw1
SSR
Ä
3777%。
7Ôo{

‚˜

og

||

‚GÖª«’©¯ÉÄ
…†äQU

78ÈÉAæçL›™ç;uGq
éª
[910]。
¦ïðñòÄ
SSR
$Ä؛™ç×~$
TG/AC
Gq]Ø

L›™ç×~$
GAA/CTT
Gq


78;uGqÔwÐ
[11]、
ÏÐ
[12]、
X¢˜
[13]
ÖÈÉ$Ø~Ÿ×ç
AG,TC,CT
骫GqÒU

¦LŸ×;u×~$ç
AAG/TTC,GAA/TCC
éª
ÔñòQU

¼_Ì»7,;u×~ÄÚeÕöÔ
É,ÚÛB©T

q
SSR
ÄGq(S

ñòâã´
SSR
GqØ&

78GqØ&

=õ§¯Ä
SSR
é
ñò‘[¸¦@Wq¥ýb

,ôd!

@WØ&¹
µî

õjs€†
(MAS,markerassistantseletion,
markeraidedseletion)、
×Ø!ø

=õ¹¾×ØT
(QTL)
µî

Ù¥µî»A>©Ùd!Ö9ą
†D!ó×ë


~=€

[1] 
$%,.

[1
[S].2010:11.
[2] 
:š

×ÞËÉví#qª»«…+Xó¯y$=
€
SSR
õÄ·¸…†
[D].
§

K22<&(o

2007.
[3] 


hϓ

ßm¬

[,LÜB¬ÄŸ ×D
RNA
#
[J].
•£<`&(
,2008(5):25
[4] MarguliesM,EgholmM,AltmanWE,etal.Genomesequen
cinginmicrofabricatedhighdensitypicolitrereactors[J].Na
ture,2005,437(7057):376.
[5] KumarS,BlaxterML.Comparingdenovoassemblersfor454
transcriptomedata[J].BMCGenomics,2010,11(2):237.
[6] ArgoutX,FouetO,WinckerP,etal.Towardstheunderstand
ingofthecocoatranscriptome:productionandanalysisofanex
haustivedatasetofESTsofTheobromacacaoL.generatedfrom
varioustissuesandundervariousconditions[J].BMCGenom
ics,2008,9(11):512.
[7] LuroFL,CostantinoG,TerolJ,etal.Transferabilityofthe
ESTSSRsdevelopedonNulesclementine(CitrusclementinaHort
exTan)tootherCitrusspeciesandtheirefectivenessforgenetic
mapping[J].BMCGenomics,2008,9(12):287.
[8] SimbaquebaJ,SanchezP,SanchezE,etal.Developmentand
characterizationofmicrosatelitemarkersforthecapegoosebery
physalisperuviana[J].PLoSONE,2011,6(10):e26719.
[9] CardleL,RamsayL,MilbourneD,etal.Computationalandex
perimentalcharacterizationofphysicalyclusteredsimplese
quencerepeatsinplants[J].Genetics.2000,156(2):847.
[10] VarshneyRK,GranerA,SorelsME.Genicmicrosatelite
markersinplants:featuresandapplications[J].TrendsBiotech
nol,2005,23(1):48.
[11] LiC,ZhuY,GuoX,etal.Transcriptomeanalysisrevealsgin
senosidesbiosyntheticgenes,microRNAsandsimplesequencere
peatsinPanaxginsengC.A.Meyer[J].BMCGenomics,
2013,14:245.
[12] 
ž

!lf



Ö

ÏÐßÄ
ESTSSR
µ@ž­µ
-õÄbc
[J].
$%$,-.
,2011,36(3):289.
[13] 
:%<

X¢˜
EST
~>Ä
SSR
víµî
[J].
úKÈÉ

2011,31(1):43.

‚-0
 
‘’

·4341·