全 文 :2016
y
4
z Ó
41
âÓ
8
}
Vol41,No.8 April,2016
[
st)(
] 20151110
[
R
]
JÖ×fñd¥
(81274010);
*HºY+jÌd¥
(JC201101);
*HSehñ
“
,ÔE-å&.£¤
”
¥
(2012001);
/01o,Ôñf23-d¥
(2014RFXXJ122);
*Hºä´Èfñij¥
(12541743)
[
uq?v
]
¢-
,
ij
,Tel:(0451)82193430,Email:mawei@hljucmnet;
$¡
,
ij
,Tel:(0451)82196224,Email:1447111706@
qqcom
[
?v
]
,
Õij¦
,
Sh
4©Ehij
,Email:459820677@qqcom
ïðñòÄâã´©óôõö÷®
SSR
øùvíµîÄ
,
,
$
,
-
,
$¡
,
¢-
(
*HSehñ hñg
,
*H /01
150040)
[
wx
]
ã½.°
Roche454GSFLX
©óÚâñòÄâã´=¼
,
ü°
454SequencingSystemSoftware
µîpnâ
ã´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Ëé
924%,SSR
¦ñòÌYâã´$å×
Äé
1342%,SSR
ÄÍcªé
797kb。
[am×S;uó¯
127
,
,
Û¶ÄæØçq$Ä
TG/AC
q
,
å×
ÄÛu
SSR
øùÄ
425%。
ñòâã´Ä©óäÙHóñòâã´ÄÌm+XÆÇ
,
ñToõñòâã´
unigene
ó¯
,
ñSñòâã´
SSR
øùå×O
,
GqØ&
,
Ø¿¹ãöO
。
[
yz{
]
ïðñò
;
âã´©ó
;SSR;
øùví
StudyonqualityevaluationofsequenceandSSRinformationin
transcriptomeofAstragalusmembranacus
CHANGYue,YANSong,LIUZhenpeng,RENWeichao,LIUJie,MAWei
(ColegeofPharmaceuticalSciences,HeilongjiangUniversityofChineseMedicine,Harbin150040,China)
[Abstract] Inthisstudy,454/RocheGSFLXsequencingtechnologywasusedtoobtainthedataoftheAstragalusmembranaceus
FourhundredandfiftyfourSequencingSystemSoftwarewasappliedtocaryoutthetranscriptionofthegroupfromscratchUsingMI
SAtools,9893unigeneswereselectedforthesequenceofthegenomeofAmembranaceus,andtheinformationofSSRlocuswasana
lyzedAccordingtotheresult,theaveragelengthofreadswas413bp,about86% ofthereadswasinvolvedinthesplicing,thelength
oftheN50was1205bp,thenumberofunigeneswasmeasuredbythewholetranscript1729SSRlociintheAmembranaceustran
scriptomeweresearched,theoccurencefrequencyofSSRwas924%,thefrequencyofSSRinthewholetranscriptomewas1342%,
theaveragelengthofSSRwas797kbOnehundredandtwentysevenkindsofcorerepeatsequenceswerefound,thedominanttype
wasTG/ACtypeofdinucleotide,itappearedtoaccountfor425% ofthetotalSSRlocusTheresultsofthesequenceofthetranscrip
tionoftheAmembranaceustranscriptomerevealedtheoveralexpression,andalargenumberofunigenessequencewasobtained,and
theSSRlocusinthegenomeoftheAmembranaceusishigh,andthetypeisdiverse,andthepolymorphismofthegeneishigh
[Keywords] Astragalusmembranaceus;transcriptomesequencing;SSR;informationofloci
doi:10.4268/cjcmm20160810
·0341·
©
vÖ
:
ïðñòÄâã´©óôõö÷®
SSR
øùvíµîÄ
$,Xñò
《
$%,.
》(2010
yÒ
)
!
,
æ ñ ò
Astragalusmembranaceus(Fisch)
BgevarmongholicusHsiao Bge
ï ð ñ ò
Amembranaceus(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Ô#
11
ÈÉ
uB.V°ïðñò,-
,
U¡p$A,o
B.$SÇ
Gd!
。
¦¡p$A,o
,°ÈÉI{¤þ8fË>??É
,
<
40
d
Äñò
。
X¥Ä¥
、
4
、
é
,
°PëM
&
-80℃
B
。
12
ÈÉu
RNA
Ä
ñòu
RNA
ÐÜ0Ä
CTAB
#
[3],
:
m¸1IÀ
:
õÈÉXLÞ
$iPë
¢¾
;
-¢âê
20mL
ª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¯
08
mË
5
mol·L-1LiCl,
O=
,
¶r
10min,4℃ 12000r
·min-1
ªS
10min,
rh
;70%
¯wßI
2
,
þ
5min
RNA,
i
50μLËN
RNA,-80℃
ì°
。
13 mRNA
&[ì
131
u
RNA
Äôõ¨©Ô!õ
ü°á£µ¤
¤rD
,
¶u
RNA
Är
、
uõ
A260/280» A260/230
nµî
。
°
AgilentRNA6000PicoKit;Agilent,
50671513,
¶u
RNA
Äs̹
、
7rrn
©!
。
132 mRNA
Ä7_Ô!õ
°
OligotexmRNA
MiniKit;Qiagen,70022;AgilentRNA6000PicoKit;
Agilent,50671513
B¼
,
n
mRNA
7_Ô!õ
。
14 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¼
,
ù:Ü
;
ü°ÚÃ
QuantifluorSTfluo
rometer,Promega,E6090,
n
cDNA
^V!õ
;
ü°
ÚÃ
Agilent2100Bioanalyzer,Agilent,2100
»
Agi
lentHighSensitivityDNAKit,Agilent,50674626
B
¼
,
n
cDNA
^Vôõ¨©
;
ìg¯P
。
15 emPCR
Äc
@P !©ó !$!
emPCR
Äc
$VXÄ
DNA
^VÄõ
;
ì
emPCR
xB
、
_
、MockAmplificationMix
ȄP
、LiveAm
plificationMix;
ü°
DNACaptureBeads
-
DNA
^V¡â
;
_
;
Äc
;
É
DNACaptureBeads;
x
DNA
^VÄ
DNACaptureBeads
į
;
©ó
=ÉhÔ
。
16 Roche454GSFLX+
©óÚrf©ó
âã´©óg¯FrH¢ÊºËÉ&(BC
DEsÊ
。
17
2=¼Ìè
、
ôõö÷
Àf=¼
454Newbler
ù«2
、
ù±ôõ
×±è
。
ü°
454SequencingSystemSoftware(ht
tp://www454com/)
µîpnâã´q2
T«
[45]。
18
©ó£»rµî
©ó£»ræ12&à Ä©óõ
(sampled
readsnumber)
Äc
,
í¨©TÄ
unigene
=øÄ
_áÎ
。
&à Ä©óõT[!r
,
í¨©
TÄ
unigene
ËõÒcÃwkc
,
ô©óý
»
,
¶Ã°æ©óõÒ÷NBXT£»
。
¦u©
ó=¼$
,
ç
15000
¯ér
,
¶=¼×1f
&
,
Så(Ä78
reads
µe¨©TØk
unigene。
·1341·
2016
y
4
z Ó
41
âÓ
8
}
Vol41,No.8 April,2016
n|&à Ä©óõMÔ¤»¨©TÄ
uni
gene
Ä=õM¥¤¦[Yâ§
,
/S7쉤1
2&à =¼õÄc
unigene
æ¶:B
£»¹
。
19
âã´
SSR
øùÄvíµî
-âã´=¼°
MISAMIcroSAteliteidentifica
tiontool(htp://pgrcipkgaterslebende/misa/)
n
SSR
µî
。
óÿé
:1/10,2/6,3/5,4/5,5/5,
6/5(
ST;uÄSó¯
/
]:;uîø=ø
);
u
Hq
SSR
$JSTSó¯Ä½;u×=]
oÒN
100[68]。
2
äÔµî
21
u
RNA
ôõµî
ôõr≥250mg·L-1,uõ≥50μg;A260/280
18~22,A260/230ê≥20。UN¨©28S∶18Sk
oÞ
15,
lq
1;RIN≥80;ñ$RNAf,f
nn
,
l+
1。
1,2
(>Þñò¥
;3,4
(>Þñò4
;5,6
(>Þñòé
。
q
1 RNA
~UN¨©
Fig1 ElectrophoresisimageoftotalRNA
+
1
ñòu
RNA
Äᣨ©
Table1 SpectrophotometryoftotalRNA
No
(>
A260/280 A260/230 ôõr/mg·L-1 uõ/μg
1
¥
2068 2331 705 282
2
¥
2079 2353 655 262
3
4
2097 2351 606 2424
4
4
2071 2325 698 2792
5
é
2074 2337 959 3836
6
é
2087 2365 653 2612
g
:
{}é
2012
y
8
z
27
{
,
õ
40μL。
qçr=¼ÕçSå
,
u
RNA
ôõHÐ
,
Õç
[¸nâã´©óg¯
。
22
©óõ=D
©2=¼õ
513Mb,
@TO
HT««vÄB¬=¼
306Mb。
=¼u?
,
=
¼õé
306805437,
×ØÏ=é
742721
Ï
,
Í
(ré
413kb;
×Ø=D
,
×Øé
742721
Ï
,
à(oÞ
20bp
ó¯=é
742602,
2=¼
é
513365449,
=¼é
306805437,
T«
Ü°é
5976%。
âã´q2T«
,8604%
Ä
reads
ÐÔT
«
,
aD
639061
Ï
reads。
ó¯T«Ä
N50
(ré
1205bp,
(r
>N50
é
3823bp,
Í(ré
10695bp。
T«T
contig15167
Y
,
78
contig
íT«Êé
12851
Y
isotig。
f#[¸T«Ê
é
isotig
Ä
contig
Ô78
isotig
[Ö
,´
Ê
10742
Y
isogroup,
uaT
12880
Ï
Isotig(
ýþf#[¸
T«Êé
isotig
Ä
contig),Unigene
=õé
9893,
Unigene
Í(ré
109712898bp。
23
©ó£»rµî
&à ©óõXT
2×102kb
O
,unigene
Ä=øµµÞ£»
,
&à =¼õXT
3×102
kbunigene
Ä=¼µsw£»
。
7e¶Þu
ñò¥
、
4
、
éÄÖ»O&©ó
,
©óõµXT
£»
,
©óõ¡¢B.«v
,
lq
2。
q
2
©ó£»rµî
Fig2 Sequencingsaturationanalysis
24 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
=øÄ^
)
924%;SSR
¦ñòÌYâã´$å×Ä
(SSR
øùY=»u
unigene
=øÄ^
)1342%;
ñò
âã´$
SSR
ÄÍcª
(
u
unigene
Ä(rÔ
SSR
=øÄ^
)797kb,
l+
2。
·2341·
©
vÖ
:
ïðñòÄâã´©óôõö÷®
SSR
øùvíµîÄ
+
2 SSR
¦ñòâã´$å×Ä
Table2 OccurenceSSRsinAstragalusmembranaceustranscriptome
;uGq =ø xGq»
/%
/%
Ícª
/kb
u(r
/bp
Í(r
/bp
îç
653 3777 507 2109 1919 12
Øç
282 1631 219 4885 4472 18
Lç
730 4222 567 1887 11433 18
ç
32 185 027 39357 732 24
/ç
14 081 009 114790 300 27
Ôç
18 104 013 81028 486 32
q+
2
ÕçSåñòâã´
SSR
,G§¯
,
q
îçTÔçÄx,ç;uGqö÷S
T
,
fxÿQÚįo
。SSR
GqØ$¦î
ç
、
Øç
、
Lçr
,
Û
963%,
®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。
25 SSR
Äƹ
¦ñòâã´Ä
1729
Y
SSR
øù$am×
S;uó¯
127
,
,
®$îçq
4
,
;
Øç
q
11
,
;
Lçq
57
,
;
çq
28
,
;
/
çq
11
,
;
Ôçq
16
,
。
7
6
,;uó¯
GqÄ;u=oØ$¦
5~10
,
XT
1173
,
ÛÌY;u=Ä
6784%;
®æ
10~15
,
;u=æ
356
,
ÛÌY;u=Ä
2059;
ÓL
Äæ
16~20
,
;u=æ
149
,
ÛÌY;u
=Ä
862%。
A°æe
,
Só¯;u=¦
5~
20
Ä
SSR
Äå×=é
1678
,
ÛÌY;u
=Ä
9705%,
Só¯;u=N
20
Ä
SSR
å×Ä=»oÒT
3%,
lq
3。
qå×Ä:m
Gqµî
,
Û¶ÄæØçq$Ä
TG/AC
q
,
å×ÄÛu
SSR
øùÄ
425%;
®æ
GA/CT
q
,
å×ÄÛu
SSR
øùÄ
394%;
«À(æ
GAA/CTT
q
,
æ
356%;CCT/GGA
q
,
æ
309%,
l+
3。
26 SSR
ÄÕ°¹öÏ
»
SSR
Õ°¹Ä;«×¼æ
SSR
µ-õ
ÄØ¿¹
,
[ é
SSR
Ä(ræÂîؿ¹O
±Ä;«ØF
,
SSR≥20bpO,Ø¿¹O;
q
3 SSR
Só¯;u=µ@
Fig3 Distributionoftherepeatsnumber
+
3
ñò×~GqÄ;u
Table3 Distributionofrepeatmotifaccordingtorepetitionfre
quencyforAstragalusmembranaceus
;uGq ;u
/%
Øç
TG/AC
q
425
Øç
GA/CT
q
394
Lç
GAA/CTT
q
356
Lç
CCT/GGA
q
309
12bp≤SSR≤20bpO,Ø¿¹$Ö; SSR≤12
bp
O
,
Ø¿¹;±
。
ñòâã´=¼øH
SSR
Ä(
rØ$¦
12~45bp,
®$
12bp≤SSR≤20bp:
B$rØ¿¹Ä
SSR
aB
861
Y
,
Ûu
SSR
Ä
4980%;SSR≥20bp:BOØ¿¹Ä SSRa 422
Y
,
Ûu
SSR
Ä
2441%。
UOB
+
,
O/×
~
SSR
ÄØ¿¹±/×~ÄØ¿¹Ü¨ö±
。
¦ñòâã´=¼=Dµî$m×
,
(roÞ
20bp
Ä
SSR
oØÇÞ±/×~
,
aýx
SSR
øù
282
Y
。
ÕçÑD71µØ¿¹ãöOÄ
SSR
¦ñò
×ØÊö
r:BOÄÜ°Ï^
。
3
ºÒÔäÒ
31
ñòâã´©óôõö÷
ñòâã´©óV°Ä¥
、
4
、
é Äu
RNA
¨©äHÐ
。
n©óVÄ
reads
ÄÍ
·3341·
2016
y
4
z Ó
41
âÓ
8
}
Vol41,No.8 April,2016
(ré
413bp,
86%
Ä
reads
ÐÔóT«
,
T«Ä
N50
(ré
1205bp。
©ó£»rä+
,
V©
Ä
unigene
=õ×ugõów1âã´ví
。
q
çräÕàuâã´Ä©óäX½Õv
,
Õç
¶=¼n[¸µî
。
32
ñòâã´
SSR
øùvíµî
SSR
úûµ@Þx,XËÉÄ×Ø´$
,
o
|
10~50kb
°B¦
1
Y
SSR。
¦ÈÉ$
,
Í
233kb
°B
1
Y
SSR;
Ý-éÈÉ$Ä
SSR
=õo
Þî-éÈÉ
,
ùÃ
2
Y
SSR
JÄÍJcé
212kb,
Ãé
646kb。
@¶ñòâã´=¼
Ä
SSR
µî
,
ñòâã´$
SSR
ÄÍcªæ
797
kb,
qîçGqTÔçGq:ì
,
S;
uó¯
127
,
,
7+ñò×Ø´8:BO§r
Ä
SSR。
o1µÈÉÄ
SSR
;u×~ª«çØ
、
L
çqéª
,
fÒUÉ,JĪ{
SSR
;uGqB
VÚÛ
。
u
m×ïðñòâã´
SSR
;uG
qª«çLçéª
,
Ûw1
SSR
Ä
4222%,
[
çVÛ»AO
,
Ûw1
SSR
Ä
3777%。
7Ôo{
、
、
og
、
||
、
GÖª«©¯ÉÄ
äQU
,
78ÈÉAæçLç;uGq
éª
[910]。
¦ïðñòÄ
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,markerassistantseletion,
markeraidedseletion)、
×Ø!ø
、
=õ¹¾×Ø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
cinginmicrofabricatedhighdensitypicolitrereactors[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
ESTSSRsdevelopedonNulesclementine(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
,
O
,
Ö
.
ÏÐßÄ
ESTSSR
µ@µ
-õÄbc
[J].
$%$,-.
,2011,36(3):289.
[13]
:%<
.
X¢
EST
~>Ä
SSR
víµî
[J].
úKÈÉ
,
2011,31(1):43.
[
-0
]
·4341·