< 최종보고서 > 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 2017 년 10 월 31 일 국민대학교산학협력단 한국임업진흥원

< 최종보고서 > 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 2017 년 10 월 31 일 국민대학교산학협력단 한국임업진흥원

< 최종보고서 > 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 2017 년 10 월 31 일 국민대학교산학협력단 한국임업진흥원

제출문 한국임업진흥원장김남균귀하 이보고서를 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 과제의최종보고서로제출 합니다. 2017. 10. 31. 주관사업기관명 : 국민대학교산학협력단 주관사업책임자 : 참여연구원 : 국민대학교김태종 국민대학교이자민

요약문 Ⅰ. 제목 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 Ⅱ. 사업배경및목적 l 한국임업진흥원에서수행하고있는임업시험중수종감정은, 국내에서유일하게진행하고있는시험임 - 국내산목재공급실적증가에따른국산재이력관리에활용 - 목재이용량확대에따라목재수입량이매년증가하는추세로, 수입수종에대한정밀식별필요 l 현미경을통한수종감정은목재세포의종류와배열등을파악 - 같은속 ( 屬, genus) 에속하는목재수종중종 ( 種, species) 단위까지불가한경우가있음 l 목재문화재분야및불법벌채목재와관련하여수종감정의중요성이대두되고있음 - 문화재수리표준시방서 상에서목재문화재의복원과수리에사용되는목재는기존목재와동일한국산재수종을사용할것을원칙으로함 - 개발도상국목재의불법벌채가증가함에따라종및산지를식별할수있는기술개발이필요함 l 기존에개발된기술을활용하여목재유통시문제가되고있는수종에대 한종식별방법마련이시급함 * 목재수종감정및종식별을위한 DNA 분석방법타당성조사 (2016 년용역 ) 에서개발한 DNA 추출기술활용 수종감정고도화로의뢰인의불만해소및불법벌채목재의유통질 서확립 - 1 -

Ⅲ. 사업내용 l DNA 분석도입이필요한문제유통수종선별 -진흥원에의뢰되는수종중문제되는수종선별및 DNA 정보수집 * Eucalyptus diversicolor, Eucalyptus marginata, Pinus rigida, Pinus palustris, Chamaecyparis obtusa (5수종) l 문제유통수종의 DNA분석유전자 ( 분석마커 ) 도출및개발 -NCBI 정보및문헌조사를통한같은수종간에식별을위한분석유전자발굴 * NCBI (National Center for Biotechnology Information) 미국국립생물정보센터 l 도출된유전자 ( 분석마커 ) 를이용한목재 DNA 분석검증 - DNA 분석검증을통한재현성및실용성확인 - 2 -

Ⅳ. 사업결과 가. 문제유통수종의 DNA 분석유전자 ( 분석마커 ) 도출및개발 1) 문제유통수종선별및수종에대한 DNA 정보수집 Eucalyptus 속과 Fokienia hodgisii 수종은전체 chloroplast 유전자가등록. 그외의수종은엽록체유전자가 partial DNA로등록됨. P. resinosa 는 NCBI 에등록된 DNA의수가 13,497개이지만, 그에비하여등록된엽록체 DNA의개수는 partial DNA로 19개에불과. Acer 속과 Chamaecyparis 속, 그리고 Thuja 속또한 partial DNA 로등록. Pinus 속들은 partial DNA라고해도그크기가약 120,000bp, 이에비하여활엽수종들은평균약 2,000bp로매우작은크기만이 partial DNA로등록. Acer saccharinum 는 1,095개의 DNA 중 1,067개의엽록체 DNA가등록되었으나모두 partial DNA 등록. (113p~114p, 표 15) 2) NCBI 정보및문헌조사를통한수종식별을위한분석유전자발굴 32편의문헌검색을통하여 trnl intron, trnl-trnf intergenic spacer, trnk/matk intron 영역들은특정계통에대해서분류학적으로만족스러운정보를나타내며, 그외의그룹에대해서는식물분류학적으로낮은분해능이나타난다는점을확인. 또한특정한몇개의유전자부위에대한연구는그외의다양한유전자부위에대한데이터가부족함확인. (28-29p) 42편의문헌검색과 NCBI 데이터를바탕으로 P. densiflora와 P. resinosa 시료에대한 rpl20 - rps18 non-coding 유전자, P. palustris 와 P. rigida 시료에대한 trnd GUC - trny GUA non-coding 유전자,C. obtusa 와 T. occidentalis 시료에대한 petg trnp non-coding 유전자, E. diversicolor 와 E. marginata 시료에대한 trnw psaj non-coding 유전자, A. saccharum 과 A. saccharinum 시료에대한 trnh-psba non-coding 유전자발굴. - 3 -

나. 도출된유전자 ( 분석마커 ) 를이용한목재 DNA 분석방법검증 1) DNA 분석검증을통한재현성및실용성확인 가. 수종간차이를확인하기위하여본실험에서얻은 P. densiflora와 P. resinosa 시료에대한 rpl20 - rps18 유전자의염기서열을비교한결과 2개의부위에서차이가존재하는것을확인.( 그림52) 연구에사용된 P. densiflora와 P. resinosa 의목재시료로부터 DNA의안정적추출확인과동시에 DNA분석방법을이용하여 P. densiflora 와 P. resinosa 시료의구분가능확인. (119p ~ 120p) 나. 수종간의차이를확인하기위하여 P. palustris 와 P. rigida 시료에서얻은 trnd GUC - trny GUA 유전자의염기서열을비교한결과 1개의염기서열에서차이가있음을확인.( 그림 57) 본연구에사용된 P. palustris와 P. rigida 의목재시료로부터 DNA의안정적추출확인과동시에 DNA분석방법을이용하여 P. palustris 와 P. rigida 시료의구분가능확인.(128p) 다. 수종간의차이를확인하기위하여 C. obtusa 와 T. occidentalis 시료에서얻은 petg-trnp 유전자의염기서열을확인할수있음. 각수종의목재시료의안정적 DNA추출, 염기서열분석을위한유전자증폭재현성확인.( 그림 60, 61) 두수종의식별가능한지표유전자제시에따른향후수종식별기술확보. (136p) 라. 두수종간차이를확인하기위하여 E. diversicolor 와 E. marginata 시료에대한 trnw psaj PCR product sequence를비교한결과대부분의염기서열이불치함확인. ( 그림 67) 두목재시료로수종으로부터 DNA의추출, 유전자증폭, 염기서열분석, 시료구분가능.(145p) 마. 수종간의차이를확인하기위하여 A. saccharum 과 A. saccharinum 시료에대한 trnh-psba PCR 유전자의염기서열을비교한결과 4개의서열에서차이가있음을확인.( 그림 72) 두목재시료로수종으로부터 DNA의추출, 유전자증폭, 염기서열분석, 시료구분가능.(153p) - 4 -

다. 결과산출물 Ÿ 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 최종 보고서 Ÿ 특허출원 (1 건 ) - 사포를이용한목분의 DNA 분리방법, 대한민국특허출원번호 : 10-2017-0045388, 2017 년 04 월 07 일 Ÿ 포스터발표 (1건) - 분자생물학적분석을활용한수종감정에서 DNA 추출의효율성증대를위한목재의전처리, 이자민, 오정애, 손석규, 김태종, 2017 한국목재공학회춘계학술대회, 전북대학교, 전주, 대한민국 (2017년 4월 14 15일 ) Ÿ 홍보 (1 건 ) 2017 년 05 월 31 일 순번언론사제목 1 내일신문낡은목재도 DNA 로수종분석 2 이코노믹톡뉴스국내최초, 목재에서 DNA 추출기술개발해특허출원까지 3 동아일보국민대임산생명공학과, 국내최초로 ' 목재 DNA' 추출기술개발 - 5 -

< 목차 > < 표차례 > 10 < 그림차례 > 12 1. 사업의개요 17 가. 사업배경및필요성 17 나. 사업목적 19 2. 사업의범위및수행방법 20 가. 문헌정보분석을통한수종식별분석유전자선별 20 나. 수종식별을위한선별유전자의확보 21 다. 확보된유전자분석및평가 26 3. 연구결과 27 가. 목재의수종식별을위한국내외유전자분석사례 27 1) Chloroplast DNA 27 2) Non-coding region 28 3) 수종별분석유전자부위조사 30 나. 소나무 (Pinus densiflora) 를이용한수종식별을위한다양한유전자검증 34 1) trnh GUG psba 37 가 ) Primer 디자인 37 나 ) DNA 추출및유전자증폭 37 다 ) 소나무 (Pinus densiflora) 의 trnh GUG - psba 염기서열 40 라 ) 소나무 (Pinus densiflora) 의 trnh GUG - psba 염기서열분석 41 2) matk 44 가 ) Primer 디자인 44 나 ) DNA 추출및유전자증폭 44 3) trns GCU 5 trng 46 가 ) Primer 디자인 46 나 ) DNA 추출및유전자증폭 46 4) 5 trng 3 trng UUC 48 가 ) Primer 디자인 48-6 -

나 ) DNA 추출및유전자증폭 48 다 ) 소나무 (Pinus densiflora) 의 5 trng 3 trng UUC 염기서열 51 라 ) 소나무 (Pinus densiflora) 의 5 trng 3 trng UUC 염기서열분석 52 5) rpob trnc GCA 54 가 ) Primer 디자인 54 나 ) DNA 추출및유전자증폭 54 6) psbm trnd GUC 57 가 ) Primer 디자인 57 나 ) DNA 추출및유전자증폭 57 다 ) 소나무 (Pinus densiflora) 의 psbm trnd GUC 염기서열 60 라 ) 소나무 (Pinus densiflora) 의 psbm trnd GUC 염기서열분석 61 7) trnd GUC trne UUC 63 가 ) Primer 디자인 63 나 ) DNA 추출및유전자증폭 63 다 ) 소나무 (Pinus densiflora) 의 trnd GUC trne UUC 염기서열 66 라 ) 소나무 (Pinus densiflora) 의 trnd GUC trne UUC 염기서열분석 67 8) trnd GUC trny GUA 68 가 ) Primer 디자인 68 나 ) DNA 추출및유전자증폭 68 다 ) 소나무 (Pinus densiflora) 의 trnd GUC trny GUA 염기서열 71 라 ) 소나무 (Pinus densiflora) 의 trnd GUC trny GUA 염기서열분석 71 9) trnfm CAU trns UGA 73 가 ) Primer 디자인 73 나 ) DNA 추출및유전자증폭 73 다 ) 소나무 (Pinus densiflora) 의 trnfm CAU - trns UGA 염기서열 76 라 ) 소나무 (Pinus densiflora) 의 trnfm CAU - trns UGA 염기서열분석 77 10) 5 trnl UAA trnt UGU 79 가 ) Primer 디자인 79 나 ) DNA 추출및유전자증폭 79 다 ) 소나무 (Pinus densiflora) 의 5 trnl UAA trnt UGU 염기서열 82 라 ) 소나무 (Pinus densiflora) 의 5 trnl UAA - trnt UGU 염기서열분석 82 11) 3 trnl UAA - trnl5 - UAA 85 가 ) Primer 디자인 85-7 -

나 ) DNA 추출및유전자증폭 85 12) 5 rps12 rpl20 87 가 ) Primer 디자인 87 나 ) DNA 추출및유전자증폭 87 다 ) 소나무 (Pinus densiflora) 의 5 rps12 rpl20 염기서열 90 라 ) 소나무 (Pinus densiflora) 의 5 rps12 rpl20 염기서열분석 91 13) psbb psbh 93 가 ) Primer 디자인 93 나 ) DNA 추출및유전자증폭 93 14) rpl16 95 가 ) Primer 디자인 95 나 ) DNA 추출및유전자증폭 95 다 ) 소나무 (Pinus densiflora) 의 rpl16 염기서열 98 라 ) 소나무 (Pinus densiflora) 의 rpl16 염기서열분석 99 15) rpl20 rps18 102 가 ) Primer 디자인 102 나 ) DNA 추출및유전자증폭 102 다 ) 소나무 (Pinus densiflora) 의 rpl20 - rps18 염기서열 105 라 ) 소나무 (Pinus densiflora) 의 rpl20 - rps18 염기서열분석 105 16) 수종분석을위한다양한유전자분석 108 다. 다양한수종에서의 DNA 안정적추출검증및염기서열분석 111 1) Pinus densiflora와 Pinus resinosa의비교 115 가 ) Primer 디자인 115 나 ) DNA 추출및유전자증폭 115 다 ) rpl20 - rps18 염기서열 118 라 ) rpl20 - rps18 염기서열분석 119 2) Pinus palustris와 Pinus rigida의비교 124 가 ) Primer 디자인 124 나 ) DNA 추출및유전자증폭 124 다 ) trnd GUC - trny GUA 염기서열 127 라 ) trnd GUC - trny GUA 염기서열분석 128 3) Chamaecyparis obtusa와 Thuja plicata의비교 132 가 ) Primer 디자인 132-8 -

나 ) DNA 추출및유전자증폭 132 다 ) petg trnp 염기서열 135 라 ) petg trnp 염기서열분석 135 4) Eucalyptus diversicolor와 Eucalyptus marginata의비교 140 가 ) Primer 디자인 140 나 ) DNA 추출및유전자증폭 140 다 ) trnw psaj 염기서열 143 라 ) trnw psaj 염기서열분석 144 5) Acer saccharum와 Acer saccharinum의비교 149 가 ) Primer 디자인 149 나 ) DNA 추출및유전자증폭 149 다 ) trnh psba 염기서열 152 라 ) trnh psba 염기서열분석 153 6) 목재의유전자식별을위한다양한수종에서의 DNA 안정적추출검증 및염기서열분석결론 157 4. 과제수행달성도 159 가 ) 문제유통수종의 DNA 분석유전자 ( 분석마커 ) 도출및개발 159 나 ) 도출된유전자 ( 분석마커 ) 를이용한목재 DNA 분석방법검증 160 5. 정책제안 162 가 ) 산학연협약을통한유전자분석을통한목재수종식별서비스구축 162 나 ) 유전자분석을통한지표유전자발굴방안 164 참고문헌 165-9 -

< 표차례 > 표1. 다양한목본수종의계통발생및분류연구에사용되는유전자부위 31 표 2. primer 3 version 0.4.0을이용하여디자인된 P. densiflora의 non-coding 유전자부위 primer 목록 36 표 3. trnh-gug와 psba primer로증폭한유전자의염기서열을 NCBI data base 에검색한결과 41 표 4. 5 trng-2s와 3 trng-uuc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 52 표 5. psbmf와 trnd-guc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 61 표 6. trnd-guc와 trne-uuc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 66 표 7. trnd-guc와 trny-gua primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 71 표 8. trnfm-cau와 trns-uga primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 77 표 9. 5 trnl-uaa와 trnt-ugu primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 82 표 10. 5 rps12와 rpl20 primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 91 표 11. rpl16f와 rpl16r primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 99 표 12. rpl20-rps18f와 rpl20-rps18r primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 105 표 13. 각 primer에따른유전자증폭여부및염기서열분석여부, NCBI blast 검색결과요약 110 표 14. 해부학적식별의한계가있는유통수종선별 111 표 15. 수종별유전정보수집 114 표 16. P. densiflora와 P. resinosa 시료에대한 rpl20-rps18f, rpl20-rps18r primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 119-10 -

표 17. P. palustris와 P. rigida 시료에대한 trnd-guc, trny- GUA primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 128 표 18. C. obtusa와 T. plicata시료에대한 petg, trnp primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 136 표 19. E. diversicolor와 E. marginata 시료에대한 trnw, psaj primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 144 표 20. A. saccharum과 A. saccharinum 시료에대한 trnh, psba primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 153-11 -

< 그림차례 > 그림 1. 염기서열분석의뢰결과로분석된각각염기서열의 sequencing 정확도를확인하는데이터결과예 25 그림 2. 양방향으로읽힌염기서열을 align한결과의예 25 그림 3. 1995년에서 2002년사이에 American Botany, Systematic Botany, Molecular Phylogenetics and Evolution, Plant Systematics and Evolution에게재된 cpdna matk 및 non coding 영역을이용한 445건의계통발생연구추세 [3] 29 그림 4. 담배엽록체의 21가지 non-coding 유전자부위 [3] 35 그림 5. P. densiflora 의 trnh GUG -psba 유전자서열과디자인한 primer의위치 38 그림 6. trnhgug와 psba primer를이용하여소나무시료로부터얻은 DNA의 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 39 그림 7. NCBI에등록된 P. densiflora 염기서열과 trnh-gug와 psba primer로증폭한유전자의염기서열을 align한결과 42 그림 8. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이외에 NCBI data base에검색되는 P. densiflora 의 trnh GUG -psba 유전자들을 align) 43 그림 9. P. densiflora의 matk 유전자서열과디자인한 primer의위치 45 그림 10. P. densiflora의 trns GCU 5 trng 유전자서열과디자인한 primer 의위치 47 그림 11. P. densiflora의 5 trng 3 trng UUC 유전자서열과디자인한 primer의위치 49 그림 12. P. densiflora 시료로부터얻은 DNA를 template로하여 primer trng F와 trng-uuc로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 50 그림 13. NCBI에등록된 P. densiflora 염기서열과 trng F와 trng-uuc primer로증폭한유전자의염기서열을 align한결과 53 그림 14. P. densiflora의 rpob trnc GCA 유전자서열과디자인한 primer의위치 55 그림 15. P. densiflora 시료로부터얻은 DNA를 template로하여 primer rpob 와 trnc-gca로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 56-12 -

그림 16. P. densiflora의 psbm trnd GUC 유전자서열과디자인한 primer의 위치 58 그림 17. P. densiflora 시료로부터얻은 DNA를 template로하여 primer psbmf와 trnd-guc-r로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 59 그림 18. NCBI에등록된 P. densiflora 염기서열과 psbmf과 trnd-guc primer 로증폭한유전자의염기서열을 align한결과 62 그림 19. P. densiflora의 trnd GUC trne UUC 유전자서열과디자인한 primer의 위치 64 그림 20. P. densiflora 시료로부터얻은 DNA를 template로하여 primer trnd-guc-f와 trne-uuc로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 65 그림 21. NCBI에등록된 P. densiflora 염기서열과 trnd-guc와 trne-uuc primer로증폭한유전자의염기서열을 align한결과 67 그림 22. P. densiflora의 trnd GUC trny GUA 유전자염기서열과디자인한 primer의위치 69 그림 23. P. densiflora 시료로부터얻은 DNA를 template로하여 primer trnd-guc-f와 trny-gua로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 70 그림 24. NCBI에등록된 P. densiflora 염기서열과 trnd-guc와 trny-gua primer로증폭한유전자의염기서열을 align한결과 72 그림 25. P. densiflora의 trnfm CAU - trns UGA 유전자서열과디자인한 primer 의위치 74 그림 26. P. densiflora 시료로부터얻은 DNA를 template로하여 primer trnfm-cau와 trns-uga로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 75 그림 27. NCBI에등록된 P. densiflora 염기서열과 trnfm-cau와 trns-uga primer로얻은유전자의염기서열을 align한결과 78 그림 28. P. densiflora의 5 trnl UAA trnt UGU 유전자서열과디자인한 primer의위치 80 그림 29. P. densiflora 시료로부터얻은 DNA를 template로하여 primer 5 trnl-uaa와 trnt-ugu로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 81 그림 30. NCBI에등록된 P. densiflora 염기서열과 5 trnl-uaa와 trnt-ugu primer로증폭한유전자의염기서열을 align한결과 83 그림 31. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻 은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이 외에 NCBI data base에검색되는 P. densiflora의 5 trnl UAA trnt UGU 유전 자의 align) 84-13 -

그림 32. P. densiflora의 3 trnl UAA - trnl5 UAA 유전자서열과디자인한 primer의위치 86 그림 33. P. densiflora의 5 rps12 - rpl20 유전자서열과디자인한 primer의 위치 88 그림 34. P. densiflora 시료로부터얻은 DNA를 template로하여 primer 5 rps12와 rpl20로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 89 그림 35. NCBI에등록된 P. densiflora 염기서열과 5 rps12와 rpl20 primer로 증폭한유전자의염기서열을 align한결과 92 그림 36. P. densiflora의 psbb psbh 유전자서열과디자인한 primer의위 치 94 그림 37. P. densiflora의 rpl16 염기서열과디자인한 primer의위치 96 그림 38. P. densiflora 시료로부터얻은 DNA를 template로하여 primer rpl16f와 rpl16r로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 97 그림 39. NCBI에등록된 P. densiflora 염기서열과 rpl16f와 rpl16r primer로 얻은유전자의염기서열을 align한결과 100 그림 40. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻 은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이 외에 NCBI data base에검색되는 P. densiflora의 rrpl16 non-coding 유전자 의 align) 101 그림 41. P. densiflora의 rpl20 - rps18 유전자서열과디자인한 primer의위 치 103 그림 42. P. densiflora 시료로부터얻은 DNA를 template로하여 primer rpl20-rps18f와 rpl20-rps18r로 1차 PCR한결과 (A) 와 2차 PCR한결과 (B) 104 그림 43. NCBI에 등록된 P. densiflora 염기서열과 rpl20-rps18f와 rpl20-rps18r primer로증폭한유전자의염기서열을 align한결과 106 그림 44. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻 은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이 외에 NCBI data base에검색되는 P. densiflora의 rpl20 - rps18 non-coding 유전자의 align) 107 그림 45. 세포의종류와구조가유사하여현미경상의식별이어려운수종의 현미경사진 ( 접선단면 ) (A) P. plaustris (B) P. rigida 112 그림 46. 세포의종류와구조가유사하여현미경상의식별이어려운수종의 현미경사진 ( 횡단면 ) (A) E. diversicolor (B) E. marginata 112-14 -

그림 47. 세포의종류와구조가유사하여현미경상의식별이어려운수종의현미경사진 ( 횡단면 ) (A) C. obtusa (B) T. plicata 113 그림 48. 목재시편에서얻은 DNA를 template로하여 rpl20 - rps18f와 rpl20-rps18r primer로증폭하여얻은유전자 (a) P. densiflora; (b) P. resinosa 116 그림 49. NCBI에등록되어있는 P. densiflora와 P. resinosa의 rpl20 rps18 non-coidng 유전자의염기서열비교및 primer 디자인 117 그림 50. NCBI에등록된 P. densiflora의 rpl20 - rps18 non-coding 염기서열과본실험에서 P. densiflora로부터얻은유전자염기서열의 align한결과 121 그림 51. NCBI에등록된 P. resinosa의 rpl20 - rps18 non-coding 염기서열과본실험에서 P. resinosa로부터얻은유전자염기서열의 align한결과 122 그림 52. 본실험에서 P. densiflora와 P. resinosa로부터얻은유전자염기서열의 align한결과 123 그림 53. 목재시편에서얻은 DNA를 template로하여 trnd-guc와 trny-gua primer로증폭하여얻은유전자 (a) P. palustris; (b) P. rigida 125 그림 54. NCBI에등록되어있는 P. palustris와 P. rigida의 trnd GUC - trny GUA non-coding 유전자의염기서열비교및 primer 디자인 126 그림 55. NCBI에등록된 P. palustris의 trnd GUC - trny GUA non-coding 염기서열과본실험에서 P. palustris로부터얻은유전자염기서열의 align한결과 129 그림 56. NCBI에등록된 P. rigida의 trnd GUC - trny GUA non-coding 염기서열과본실험에서 P. rigida로부터얻은유전자염기서열의 align한결과 130 그림 57. 본실험에서 P. palustris와 P. rigida로부터얻은유전자염기서열의 align한결과 131 그림 58. 목재시편에서얻은 DNA를 template로하여 petg와 trnp primer 로증폭하여얻은유전자 (a) C. obtusa (b) T. plicata 133 그림 59. NCBI에등록되어있는 C. obtusa와 T. occidentalis의 petg - trnp non-coding 유전자의염기서열비교및 primer 디자인 134 그림 60. NCBI에등록된 C. obtusa 의 petg - trnp non-coding 염기서열과본실험에서 C. obtusa로부터얻은유전자염기서열의 align한 137 그림 61. NCBI에등록된 T. occidentalis의 petg - trnp non-coding 염기서열과본실험에서 T. plicata로부터얻은유전자염기서열의 align한결과 138-15 -

그림 62. 본실험에서 C. obtusa와 T. plicata로부터얻은유전자염기서열의 align한결과시료에대한 139 그림 63. NCBI에등록되어있는 E. diversicolor와 E. marginata의 trnw psaj 유전자간의염기서열비교및 primer 디자인 141 그림 64. 목재시편에서얻은 DNA를 template로하여 trnw, psaj primer로증폭하여얻은유전자 (a) E. diversicolor (b) E. marginata 142 그림 65. NCBI에등록된 E. diversicolor의 trnw psaj non-coding 염기서열과본실험에서 E. diversicolor로부터얻은유전자염기서열의 align한결과 146 그림 66. NCBI에등록된 E. marginata의 trnw psaj non-coding 염기서열과본실험에서 E. marginata로부터얻은유전자염기서열의 align한결과 147 그림 67. 본실험에서 E. diversicolor와 E. marginata로부터얻은유전자염기서열의 align한결과 148 그림 68. NCBI에등록되어있는 A. saccharum과 A. saccharinum의 trnh - psba 유전자의염기서열비교및 primer 디자인 150 그림 69. 각수종으로부터얻은 DNA를 template로하여 trnh, psba primer로증폭하여얻은유전자 (a) A. saccharum (b) A.saccharinum 151 그림 70. NCBI에등록된 A. saccharum 의 trnh - psba non-coding 염기서열과본실험에서 A. saccharum로부터얻은유전자염기서열의 align한결과 154 그림 71. NCBI에등록된 A. saccharinum의 trnh - psba non-coding 염기서열과본실험에서 A. saccharinum로부터얻은유전자염기서열의 align한결과 155 그림 72. A. saccharum과 A. saccharinum로부터얻은유전자염기서열의 align한결과 156-16 -

1. 사업의개요 가. 사업배경및필요성 목재의수종식별은다음과같은목적에서필요하다. l l l l l l 목재활용을위한필수기초자료로이용이용목적에부합하는목재사용판별목조문화재의복원을위한동일한목재수종검증 ( 목재문화재는복원과수리를위해기존목재와동일한국내산수종으로대체해야함 ( 문화재수리표준시방서, 2014)) 법률적관리가필요한목재수종의관리목재한약재의감정목재의산지를식별 한국임업진흥원에서는목재의수종감정서비스를국내유일하게제공하고있다. 2014년에수종감정서비스가수요가 1,540건이었던것이 2015 년에 2,682건으로급격히 (74%) 증가하고있다. 현재의목재수종감정은현미경을이용하여해부학적으로목재가가지고있는세포조직의종류와배열등을관찰하고이를비교하여식별하는방법으로이루어지고있다. 이러한방법은가까운수종끼리식별하기어려운경우가일부발생하고있다. 이러한해부학적목재식별방법의한계를보완하기위한방안으로분자생물학적기술을이용한유전자분석을통하여수종을식별하는연구가주목받고있다. 그러나과거에는목재로부터 DNA를추출하고이를활용하여원하는유전자의염기서열을안정적으로얻을수없었다. 본연구에서는현재해부학적목재식별방법에의해감정이어려운 5쌍 (10종) 의수종을선정하여유전자를분석을통해주어진비교목재시료사이의식별이가능한지실증하고자한다. 이는선행연구 ( 용역과제계약번호 : 20170205951-00) 에서제안한목재로부터 DNA를추출하는새로운방법을다양한다른수종에적용하여실증함으로써제안된새로운 DNA 추출방법의가능성을확대하여실증한다. 비록많은수종의경우대표성을가지는충분한유전자가확보되지못해상 - 17 -

업적목재의수종식별을위한지표유전자가제안되지못하고있으나, 일부수종의경우지표유전자가제안되고있으며꾸준한연구를통해목재수종식별을위한기반자료가축적되고있다. 향후상업적수종식별을위한지표유전자가제안된다면본연구의결과를이용해다양한목재에서 DNA를추출하고유전자분석을통해목재의수종식별을위한가능성을제시하고자한다. 목재는이미죽은세포로다양한환경에오랜시간노출되어세포내의 DNA가절편화되어있으며추출하기에어려운상태이다. 이와같은이유로유전자분석을통하여수종을식별하는연구가주목받고있음에도불구하고상업적수종식별을위한방법의개발에진전이없었으며목재의식별을위한지표유전자의발굴에도실용성의한계때문에연구투자가미흡하였다. 본연구를통해다양한목재시료로부터안정적으로 DNA를추출하고, 이를통해원하는유전자의염기서열을확보될수있으며확보된염기서열을비교하여목재시료간의식별이가능함을실증하고자한다. 이는분자생물학적유전자분석에의한목재수종식별을위한방법의개발에있어가장큰문제였던목재로부터 DNA의추출방법을제시하여해결함으로써유전자분석에의한목재수종식별의실용성을제안한다. 본연구의실증을통해목재수종식별을위한지표유전자의발굴에많은연구가진행되고향후한국임업진흥원에서해부학적목재수종식별의보완방법으로유전자분석을통한목재수종식별서비스제공을위한계기가될수있을것이다. - 18 -

나. 사업목적 l 앞선용역과제 목재수종감정및종식별을위한 DNA 분석방법타당성 조사 의성공적인수행을통해유전자를이용한목재의수종식별타당성 확인 l 본연구에서는앞선용역과제의타당성결과를이용하여현재수종분석서비스가제공되고있으나, 해부학적수종식별에어려움이있는수종을중심으로 DNA에의한목재수종감정을실증하고서비스에적용하는기초가되고자함 - 19 -

2. 사업의범위및수행방법 가. 문헌정보분석을통한수종식별분석유전자선별 한국임업진흥원에의뢰되는여러수종중해부학적으로세포의종류와그배열이비슷하여현미경관찰을통한식별이어려운수종을 5쌍 (10종) 을선정하였다. 1991년부터 2014년도까지의목본식물의계통발생과식물분류에활용된유전자연구의참고문헌을조사하여현재계통분류학적으로활용되는유전자부위와그한계를확인하였다. Non-coding부위의 trnt, trnl-trnf, atpb-rbcl intergenic spacer, trnk/matk intron 부위의연구가많이이뤄졌으나특정수종에대해서만분류학적으로적절한정보를제공하며, 그외의여러수종에대해서는낮은분해능으로인하여추가적인유전자부위의연구가필요하다. 선정된두수종 (5쌍) Eucalyptus 속, Pinus 속, Acer 속, Thuja 속, Chamaecyparis 속관련하여분자생물학적인방법으로각쌍의두수종을식별하기위한유전자부위에대하여조사하였다. Pinus 속의경우식물계통, 분류학적연구가다른속에비하여비교적많은연구가되어있었지만그외의수종에대한연구는부족하였다. 또한각속수준에서특정수종에대해서만연구가주로되어있었고, 그외의종에대한연구는부족하였다. 이는유전자의염기서열을확인하여종수준의식별을하기위한분자생물학적연구가다양한종에이뤄지지않고있음을확인할수있고, 따라서보다나은분자생물학적방법을이용한수종감정을위해서는다양한수종에대한여러유전자부위의확인및그검증이필요하다. - 20 -

나. 수종식별을위한선별유전자의확보 많은식물의분류연구에상호비교에많은장점을가지고있는엽록체유전자를이용하고있으며, 이를통해많은엽록체유전자정보가축적되어가고있다. 목재식별을위한유전자를선정하는데이와같은유리한점이있는엽록체유전자를이용하는것이바람직하다. 또한, 목재는이미죽은세포로건조되어오랜시간동안다양한환경에노출된다. 이로인해목재안에존재하는 DNA는대부분손상되고절편화되어있다. 따라서수종식별을위한유전자는하나의세포에한쌍씩만존재하는 chromosomal DNA보다는많은수가존재하는엽록체의 DNA를확보하는것이보다안정적인목재식별을위해유리하다. 일반적으로수종식별을위한분석유전자는엽록체 DNA에서생물학적기능에관련되지않는염기서열인 non-coding 부위를이용한다. 그이유는생물학적기능과관련된유전자는염기서열변화 ( 돌연변이 ) 에의해기능이상실될가능성이있으며, 이로인해세포가사멸하게되어다양한염기서열의변화를관찰하기어려운문제점이발생하기때문이다. 그러나유전자를 non-coding하는부위는염기서열의변화를쉽게받아들일수있다. 하지만무엇보다식별이필요한두수종간의관계에따라 coding 부위와 non-coding 부위를적절히선택하여이용하는것이바람직하다. 목재의수종을식별하기위해서는시료의관계가중요하다. 비교하고자하는목재수종의관계에따라분석하는유전자를달리하여야하며, 동일한목재의시료를식별할때, 시료사이의분류학상차이가속 (genus) 간인지, 종 (species) 간인지, 아종간인지, 지역의집단간인지, 개체간인지에따라분석하여야하는유전자가달라져야한다. 따라서다양한수종식별을위한분석유전자를결정하기위해서는충분히많은대표성을가지는시료를확보하고이들의유전정보를체계적으로정리하여지표유전자를도출하는것이바람직하며이를위해많은연구가필요하다. 본연구에서는추후의연구에의해수종식별을위한적절한지표유전자가제시되면본연구그룹이새로이제안한특허방법에의해목재시료로부터유전자의확보가가능하다는것을확증하고자한다. 이에본실험에서는상업적으로이용가능한지표유전자를제시하지않지만문헌을통해제시된유전자를선정하고이를이용하여본연구에서주어진 5쌍의시료가유전자분석에의해수종식별이가능한지실증하였다. 이를위해다음과같은연구단계를수행하였다. - 21 -

l l l l 문헌분석을통해선별된유전자증폭을위한 primer 제작목재로부터 DNA 추출중합효소연쇄반응 (Polymerase Chain Reaction) 을이용한유전자증폭증폭된유전자의크기를이용한확보유전자의평가 구체적인실험방법은다음과같다. 1. 시료준비 1) 60grit 사포를이용하여멸균된시편의여섯면을사포질통하여표면을 두께약 0.1cm 제거하고, 사용한사포와시료를모두버린다. 2) 70% ethanol 을이용하여실험대의불순물을제거, 세척한다. 3) 멸균된새로운사포한장을꺼내어시편의횡단면방향으로갈아준다. 4) 분쇄한시료를 1.5ml 에펜튜브에 20mg 씩소분한다. 5) 수화를위해멸균된 3 차증류수 200 μl를첨가하여 vortex 한다. 6) 2 일동안 4 냉장고에서수화시키고, 24 시간마다 vortexing 하여시료 와용액이혼합되도록한다. 2. DNA 추출 _ Qiagen DNeasy Plant Mini kit (Cat. No. 69104) 1) 4 냉장고에서꺼내어 3 차증류수를처리한시료에 Ap1 buffer 600 μl, RNase A 6 μl (100mg/ml) 를첨가한다. 2) Vortexing 을강하게하여모든시료가충분히섞이도록한다. 3) 65 에서 10 분간반응시키면서 2 분간격으로꺼내어위아래로뒤집으 - 22 -

며섞는다. 4) P3 buffer 260 μl를넣고위아래로뒤집으며섞는다. 5) 충분히섞인시료는 5 분간얼음속에넣어정치시킨다. 6) 13,500rpm 으로 10 분간원심분리한다. 7) 끝을자른팁을이용하여상층액만취한후 QIAspin column 에옮겨담 은후 13,500rpm 으로 2 분간원심분리한다. 8) Collection tube 에모인용액을새로운에펜튜브 (#1) 에옮겨담는다. 9) 7) 번과정을반복하여, 또다른새로운에펜튜브 (#2) 에옮겨담는다. 10) 얻어진용액의 1.5 배의 AW1 buffer 를각에펜튜브에첨가한후바로 끝이잘린팁을이용하여피펫팅을한다. 11) 650 μl를따서 DNeasy Mini spin column 에옮겨담는다. 12) 8,000rpm 에서 1 분간원심분리를한다. 13) 흘러내려온용액은버리고, 나머지용액을 column에넣은후 8,000rpm 에서 1분간원심분리를한다. ( 에펜튜브 #1, #2에들어있는모든용액을하나의 DNeasy Mini spin column에 11~13번과정으로반복하여진행한다.) 14) Column 에새로운 collection tube 를장착한후 column 에 AW2 buffer 를 500 μl넣고 8000rpm 에서 1 분간원심분리한다. 15) 흘러내려온용액만버리고 collection tube 를재장착한후, AW2 buffer 500 μl를넣고 13,500rpm 에서 2 분간원심분리를한다. 16) 흘러내려온용액과 collection tube 를버리고새로운에펜튜브를장착 - 23 -

한후킴테크를덮어놓고상온에서 column 40 분동안건조시킨다. 17) AE buffer 50 μl를 column 중앙에직접적으로분지하고 5 분간정치한 다. 18) 8,000rpm 에서 1 분간원심분리한다. 19) -20 냉동보관한다. 3. Primer 디자인및유전자증폭 1) NCBI 데이터를바탕으로 primer3 (V.0.4.0) 를이용하여디자인한다. (http://bioinfo.ut.ee/primer3-0.4.0/) 2) DNA template, dntp, buffer, taq polymerase, primer forward, primer reverse 를혼합한다. 3) 일반적으로타겟으로하는유전자부위에 primer가붙을수있는적정온도는디자인된 primer 마다다르다. 이에맞는온도조건을조정하여 PCR을진행한다. ( 실험에서진행된온도조건은아래에구체적으로서술해놓았음.) 4) 전기영동 (electrophoresis) 을통하여타겟으로하는유전자가증폭되었 는지확인한다. 4. 염기서열분석의뢰 1) 증폭시킨유전자산물의염기서열분석을의뢰한다. 2) 결과데이터인염기서열 text file과 chromatogram file을비교를통하여염기서열분석시이뤄지는오류를확인한다. ( 일반적으로서열의초반 50bp까지는분석이깨끗하게이뤄지지않기때문에인위적으로삭제해주며, 이를보완하기위하여양방향으로염기서열분석을의뢰한다.) - 24 -

그림 1. 염기서열분석의뢰결과로분석된각각염기서열의 sequencing 정 확도를확인하는데이터결과예. 3) Vector NTI 프로그램을활용하여분석된양방향의염기서열을 align 한 다. 그림 2. 양방향으로읽힌염기서열을 align 한결과의예 4) Align 된염기서열을 assemble 하여하나의염기서열을만든다. - 25 -

다. 확보된유전자분석및평가 행하였다. 확보된염기서열의분석을위해다음과같은 4 단계의비교분석을수 Ÿ Ÿ Ÿ Ÿ 유전자의비교분석을위해선별된유전자의 Data base (NCBI 자료 ) 상의염기서열을두시료의등록된염기서열을비교분석하여원하는염기서열의차이가무엇인지제시하였다. 확보된비교수종 1번의유전자염기서열을이용하여기존의 Data base (NCBI 자료 ) 에서등록된유전자를탐색하여유사유전자들을도출하였으며이들의비교를통해수종을확인하였다. 확보된비교수종 2번의유전자염기서열을이용하여기존의 Data base (NCBI 자료 ) 에서등록된유전자를탐색하여유사유전자들을도출하였으며이들의비교를통해수종을확인하였다. 확보된비교수종 1번의유전자염기서열과확보된비교수종 2번의유전자염기서열을비교하여두수종간의유전자차이를확인하고앞서제시된기존의 Data base (NCBI 자료 ) 에서제시된차이점과비교하였다. 이를통해본연구에서제시된두수종간유전자염기서열을비교하 여수종식별가능성을평가하였다. - 26 -

3. 연구결과 가. 목재의수종식별을위한국내외유전자분석사례 1) Chloroplast DNA 대표적으로식물의식별을위한유전자로 chloroplast DNA (cpdna), mitochondria DNA, nuclear DNA의연구가이뤄지고있다.[30] 그중 cpdna 는식물이살아가는데필수적인유전자로써식물식별및식물종간의유전적관계를확립하기위한도구로광범위하게사용되고있다.[1,2] 서로다른엽록체의유전자좌 (locus) 는진화적으로가깝고, 먼거리에있는식물종을구별하기위해사용되어왔지만모든식물종을구별할수있는단일유전자좌는발견되지않은것으로알려져있다.[1] 이러한엽록체 DNA는다른 DNA에비하여단일세포당많은 copy수를가지고있기때문에상대적으로잔존되어있을확률이높은장점이있으며, 침엽수를제외한수종에서는모계유전이이뤄지고, 돌연변이가낮으며, 감수분열에의한교차가없기때문에과 (family), 강 (order) 수준에서고등식물에대한계통발생학적인정보를제공하기에적절한유전자라고알려져있다.[30] 초기의식물분자연구는이러한엽록체 DNA를이용한제한효소다형성연구를통하여식물분자수준의연구를진행하였으나 DNA sequencing 기술이가능해짐에따라 cpdna의유전자서열에대한비교연구가시작되었다.[7,8,9] 1993년, rbcl (ribulose-1,5-bisphosphate carboxylase) 의서열을이용하여속 (genus) 보다상위분류에대한연구가진행되었다.[10] 이유전자는광합성에중요한역할을하는유전자서열로이를이용한 DNA 종의식별에서식물의 과 와 강 을확인하는데사용되나대부분의식물에똑같은유전자염기서열이존재하기때문에보다자세한식물의식별이어렵다는단점이있었다.[30] - 27 -

2) Non-coding region 대부분의식물의분류학적으로낮은수준의종간관계를해결하기위한 coding (genic) DNA의한계로인하여엽록체의 non-coding 시퀀스의비교분석을통한식물분류학연구에대한관심이증가하고있다.[3,6] ( 그림3) 이러한 non-coding 부위는염기서열의돌연변이발생빈도가높고, coding (genic) 부위에비해진화의속도가빠르다는특징을가지고있기때문에, 진화연구와종간유전자마커개발에많이활용되고있고, 실제로상당수 non-coding 영역의계통발생학연구에서잠재적인유용성을입증하였다.[4,5] 이러한 non-coding 영역은 trnt, trnl-trnf, atpb-rbcl intergenic spacer, trnk/matk intron 영역이우선적으로연구가되었고, 그림3 에서확인할수있듯이, 445건의연구중 342건 (77%) 이 trnk/matk intron영역또는 trnl-trnf intergenic space영역을이용하여식물분류학연구를진행하고있다. 이는엽록체 DNA를이용한계통발생연구가매년증가함에따라특정유전자부위에대한의존도역시증가하고있음을의미한다. 다시말해, 특정유전자에대한연구가편중되어있고, 이에비해다른유전자부위는느린속도로진행되고있음을의미한다. 문제는 trnl intron, trnl-trnf intergenic spacer, trnk/matk intron 영역들은특정계통에대해서분류학적으로만족스러운정보를나타내며 [11], 그외의그룹에대해서는식물분류학적으로낮은분해능이나타난다는점이다 [39]. 또한특정한몇개의유전자부위에대한연구는그외의다양한유전자부위에대한데이터가부족하다는것을의미한다. - 28 -

그림 3. 1995년에서 2002년사이에 American Botany, Systematic Botany, Molecular Phylogenetics and Evolution, Plant Systematics and Evolution에게재된 cpdna matk 및 non coding영역을이용한 445건의계통발생연구추세 [3] - 29 -

3) 수종별분석유전자부위조사 1991년부터 2014년도까지의목본식물의계통발생과식물분류에활용된유전자연구의참고문헌을조사한결과, 32편의문헌중활엽수를대상으로연구된결과가주를이루었으나연구된수종이다양하지않음을확인할수있었다. 침엽수를대상으로한유전자연구는 12편에불과하였다. 주로연구가되고있는 non-coding 영역인 trnt, trnl-trnf, atpb-rbcl intergenic spacer, trnk/matk intron 부위에관련한연구논문은 32편의중 22편에서확인할수있었고, 그외의유전자부위에대한연구는미비한것을확인할수있었다. 이는각수종의다양한유전자에대한연구가아직부족함을의미한다. 이번과제에서연구될수종에따라살펴보면, Eucalyptus 속, Pinus 속, Acer 속, Thuja 속, Chamaecyparis 속중, Pinus 속과 Acer 속에관한유전자연구가가장많았으나본연구에서사용될수종에대한구체적논문은 Acer 속에서찾아보기어려웠다. 이는속수준의연구는비교적많이이뤄지고있으나다양한수종에대해서는연구가미비함을의미한다. Eucalpytus 속은대부분의연구된수종이 Eucalyptus globulus 로매우제한적이었다. 또한 Thuja 속, Chamaecyparis 속에관한연구논문이가장미비한것으로나타났다. 결론적으로같은속에대해서다양한수종에대한연구가미비할뿐만아니라한수종에대한식별가능한여러유전자부위의연구또한부족한것을확인할수있다. - 30 -

표 1. 다양한목본수종의계통발생및분류연구에사용되는유전자부위 수종유전자참고문헌 Pinus nigra ( 유럽곰솔 ) trnl gene, trnt trnl 5' exon, trnl 3' exon, trnf 5 Cunninghamia Lanceolata ( 넓은잎삼나무 ) psbc and trns intergenic spacer, ycf3 intron, trnn - trnr intergenic spacer Aquilaria malaccensis ( 침향속 ) trnl-trnf intergenic spacer 26 Pinus sp. ( 소나무속 ) rbcl, matk, trnh -psba 27 Cryptomeria japonica ( 삼나무 ) rbcl 29 24 침엽수 C. pisifera, C. obtusa, C. lawsoniana, C. formosensis, C. taiwanensis, C. thyoides, C. nootkatensis, Thuja occidentalis, Platycladus orientalis, Cupressus cashmeriana,calocedrus formosana, Juniperus chinensis trnv intron, petg - trnp intergenic spacer Thuja plicata ( 서약측백나무 ) microsatellit 34 P. densiflora, P. densiflora for. multicaulis, P. sylvestris, P. rigida, P. rigitaeda, P. koraiensis, P. bungeana Pinus densiflora, Pinus sylvestris, Thuja plicata, Acer carpinifolium. Juniperus ( 노간주나무속 ), Picea ( 가문비나무 ), Pinus sp. ( 피너스속 ) trnt(ugu) - trnl(uaa) intergenic spacer 35 trns gene and the adjacent psbc gene 38 matk 39 33 Mexican pine species, Picea species (Pinus nelsonii, P. pinceana, P. hartwegii ) trnl gene, trnt trnl 5' exon, trnl 3' exon, trnf 40 Strobus,Cathaya argyrophylla, Larix decidua, Pseudolarix amabilis, Keteleeria rbcl, matk, trnv intron, rpl20-rps18 41-31 -

davidiana, Abiesnumidica, Pinaceae Ginkgo biloba ( 은행나무 ), Salix babylonica ( 수양버들 ), Robinia pseudacacia ( 아까시나무 ), A. pseudoplatanus ( 개버즘단풍 ), A. platanoides ( 노르웨이단풍 ) trnl gene, trnt trnl 5' exon, trnl 3' exon, trnf 5 활엽수 Cyclobalanopsis species ( 가시나무아속 ) trnl-trnt intergenic spacer 13 Quercus petraea ( 페트라참나무 ) trnd - trnt, 5S ribosomal DNA, mitochondrial: cox2-1/2 nad4-3/4 14 Quercus ( 상수리나무 ) trnd - trnt intergenic spacer, trnc - trnd intergenic spacer, trnk intron 15 Quercus ( 상수리나무 ) trnl intron trnd - trnt intergenic spacer 16 48 French oak populations (comprising Quercus robur trnd - trnt intergenic regions, tnrc and Quercus petraea) ( 참나무 ) -trnd intergenic regions, trnk intron Shorea selanica, Shorea leprosula, Shorea parvifolia trnk intron, 5 trns, trng intron, ORF77 ( 사라수속 ), Dryobalanobs, Hopea plagata, Dipterocarpus - ORF82 intergenic, rpl2-rps19 kerrii intergenic, rpl2-rps19 intergenic spacer Quercus ( 상수리 ) trnt - trnf intergenic spacer, trnd - trnt intergenic spacer Olive wood trnt-trnl spacer, ITS1 20 Quercus robur ( 상수리 ) trnd-trnt intergenic spacer, trnc - trnd intergenic spacer, psa-trns intergenic spacer 17 18 19 21 Cannabis sativa ( 대마, 뽕나무과의 1 년생식물 ) Nuclear and chloroplast DNA: Inter Transcribed spacer(its) region of the ribosomal DNA ITS4, ITS5, trnl - trnf intergenic spacer 22-32 -

Dipterocarpaceae ( 이엽시과식물 ), dipterocarp wood, non-dipterocarp wood trnl - trnf intergenic spacer 23 Catalpa bungei ( 오동나무 ), Fraxinus chinesis ( 물푸레나무 ) psbc and trns intergenic spacer, ycf3 intron, trnn - trnr intergenic spacer 24 Shorea species ( 이엽시과 ) trnl intron, trnl - trnf intergenic spacer, trnh - trnk intergenic spacer, trnk - psbc intergenic spacer Quercus petraea rbcl, matk, trnh -psba 27 Sandal wood ( 단향속 ) rbcl, matk, trnh-psba 28 Eucalyptus globulus microsatellites 31 25 Eucalyptus globulus, E. dalrympleana, E. archeri, E. morrisbyi, E. vernicosa Aceraceae (Acer palmatum,a. robustum, A. miaotaiense, A. negundo, A. truncatum ) rpl2 - trnh, rpl2 - rpl22 32 matk, rbcl, ITS, trns-trng 36 Blank mape, sugar maple ndha intron 37 Acer and Dipteronia ( 칠엽수아과 ) psbm-trnd and trnd-trnt intergenic spacer 42 Acer pycnanthum ( 꽃단풍 ) trnl - trnf intergenic spacer, atpb-rbcl intergenic spacer, trnh-psba intergenic spacer 43-33 -

나. 소나무 (Pinus densiflora) 를이용한수종식별을위한다양한유 전자검증 선행연구 ( 용역과제계약번호 : 20170205951-00) 에서새로운방법에의해엽록체유전자가추출되고확보된소나무 (P. densiflora) 목재로부터 DNA 를추출하는새로운방법을제안하고이를활용하여 3가지수종 ( 소나무, 구주적송, 해송 ) 을성공적으로식별할수있음을보여주었다. 2005년발표된논문에서담배로부터 21개의다양한엽록체 non-coding 부위를유전자로분석하여계통학적분석에이용한연구가보고되었다.[3] 이에본연구에서는제안된 21개의유전자부위를 National Center for Biotechnology Information의 data base에등록되어있는소나무 (P. densiflora) 엽록체의전체유전자 (NCBI 등록번호 : JN854210) 와그림 4의 noncoding 부위를비교하여 15개의유전자부위를확인하였다. 이 15개소나무유전자를 P. densiflora 목재시료로부터확보할수있는지확인하기위하여제안된새로운방법을이용하여 DNA를추출하고이를 template로이용하여중합효소연쇄반응을수행하여유전자를확보하였다. 확인된 15개유전자부위를증폭하기위해소나무 (P. densiflora) 엽록체의전체유전자 (NCBI 등록번호 : JN854210) 를바탕으로 primer 3 version 0.4.0을이용하여 primer를디자인하였다.( 표 2) 증폭하고자하는유전자의크기는 1,000개의염기쌍이하이었으며 15개유전자의평균염기서열크기는 630개의염기쌍이었다. Primer 의 melting 온도 (Tm) 는 60 부근이었다. GC 함량 (GC%) 은 36 63% 로분석하고자하는유전자의부위에따라다양하였으나대부분은 45 55% 범위의 GC 함량을나타냈다. - 34 -

그림 4. 담배엽록체의 21 가지 non-coding 유전자부위 [3] - 35 -

표 2. primer 3 version 0.4.0 을이용하여디자인된 P. densiflora 의 non-coding 유전자부위 primer 목록 연번 유전자부위 방향 Primer 이름 염기서열 (5 3 ) T m GC% 크기 (bp) 1 trnh GUG - psba F trnh-gug ttggtccacttggctacgtc 61.1 55 R psba atgcacgaacgtaatgctca 60.29 45 653 2 matk F matk cagaaatttgatccgggttg 60.3 45 R matk ttcctggcaaaagatacagga 59.69 42.86 274 3 trns GCU - 5'trnG F trns-gcu gtccgctcagccatctctc 61.09 63.16 R trng R ctatacccgccacgatctga 61.01 55 873 4 5'trnG - 3'trnG UUC F trng F cagatcgtggcgggtatagt 59.98 55 R trng-uuc cgttagcttggaaggctagg 59.13 55 797 5 rpob - trnc GCA F rpob ggattgtaaacgctccctca 60 50 R trnc-gca gaactggggatggaggattt 60 50 510 6 psbm - trnd GUC F psbm actggccgttttgacgtaag 60.17 50 R trnd-guc R agagtaccgccctgtcaaga 59.87 55 673 7 trnd GUC - trne UUC F trnd-guc F cagggcggtactctaaccaa 60.12 55 R trne-uuc tctctttcaaggaggcaacg 60.51 50 475 8 trnd GUC - trny GUA F trnd-guc cagggcggtactctaaccaa 60.12 55 R trny-gua atgcccgagtggttaatgag 59.96 50 354 9 trnfm CAU - trns UGA F trnfm-cau acgggttcaaatcctgtctc 58.99 50 R trns-uga cattaaccactcggccatct 59.96 50 956 10 5 trnl UAA - trnt UGU F 5 trnl-uaa ccgtagcgtctaccaattcc 59.59 55 R trnt-ugu tgcgatgctctaacctctga 59.7 50 515 11 3 trnl UAA - trnl5' UAA F 3 trnl-uaa gacttgaaccctcacggtct 59.15 55 R trnl5'-uaa ggaattggtagacgctacgg 59.59 55 552 12 5 rps12 - rpl20 F rps12 ctcctgctcttcgaggatgt 59.55 55 R rpl20 ttccgacgttttcgagctat 59.85 45 805 13 psbb - psbh F psbb tccaaaaactgggagatcca 60.43 45 R psbh tggcgtggttttggaagtat 60.37 45 564 14 rpl16 F rpl16f tgcttagtgtgcgactcgtt 59.66 50 R rpl16r tcctttcattcttcccctatg 58.12 42.86 887 15 rpl20 - rps18 F rpl20-rps18f ccaaatcgtaaaacacttgcac 59.63 36.36 R rpl20-rps18r ttatccggccgagtgaatag 60.05 50 560-36 -

1) trnh GUG - psba 가 ) Primer 디자인 이두유전자사이에존재하는 non-coding region은비교적평균유전자길이가짧기때문에대부분의계통에서유전자증폭과염기서열분석이다고보고되어있다.[3] 이에본연구에서는그림 5와같이 P. densiflora의 NCBI data base를바탕으로, trnh GUG 와 psba coding region에서 primer를디자인하였으며, 이를이용하여두유전자사이의 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnh-gug와 psba를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의용량으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 52 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림6 (A) 에서확인할수있듯이약 600bp에서증폭하고자하는유전자부위사이즈의밴드를확인할수있었으나그농도가매우약하여염기서열분석이어렵다고판단되었다. 이에그림6 (A) 의화살표로표시된증폭된유전자부위를분리정제하고, 이를 template로하여 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 27회로변경하여진행하였다. 그림6 (B) 에서증폭하고자하는유전자크기인약 600bp에서농도가충분히확보되어염기서열을분석하였다. - 37 -

그림 5. P. densiflora 의 trnh GUG -psba 유전자서열과디자인한 primer 의위치 - 38 -

그림 6. trnh-gug 와 psba primer 를이용하여소나무시료로부터얻은 DNA 의 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 39 -

다 ) 소나무 (Pinus densiflora) 의 trnh GUG -psba 염기서열 실험에서 trnh-gug와 psba primer를이용하여증폭된유전자를양방향으로각각염기서열을분석한결과를 vector NTI program을통하여 align한후에 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TTTGGTCCACTTGGCTACGTCCGCCCCGGAAAAAACCAATTTATCTATCTACAG TCATTTCTTCCAAGAATAAAAAATGAGCTAACGTTTGTTCTTATGTGGGTCGGT GACCTCTGATAGGGGATCAAAGCAAGATCGATGACTAACTCTCAACTCTCGAA CAAGTTGTATGGCTTATTATCAAATTAATTCAATGAAATGTTATTTGAATGTCT ATTGAGAATATTTGACATGAATCAAGTATGAGAGAAAGAATGTAAATGGGTCC CAATCCCGAACTACCCAAATTTAGATCTGAGTCTATACTCATATTATAGAATGA ATATGTTGATGATCTTTATCCAGCATCCATGGCTGAATGGTTAAAGCGCCCAAC TCATAATTGGCGAACTCGCGGGTTCAATTCCTGCTGGATGCAGCGGAACTGCA CATATCCATTATTGATGGAATGGGAAGAAGTATGAAGAATAACACCAAACAAT TGAAGCATACCAAGCCTTTCAATAAAATGAATGAAAGGCTTGGTATGCTTCAAT TAAGCAATACACGATAACAATCCATCAGAGTATTATCCACCTATTGAAATAGAT TCAACAGCGGCTAGATCCAGAGGAAAGTTGTGAGCATTACTTTTCGTGCAT - 40 -

라 ) 소나무 (Pinus densiflora) 의 trnh GUG -psba 염기서열분석 얻어진 trnh GUG -psba의염기서열을 NCBI의 nucleotide blast에서검색한결과, 표 3과같이 Pinus taeda가 query coverage 99%, identities 100% 로검색되었다. 한편연구에사용된목재인 P. densiflora는 query coverage 98%, identities 97% 로검색되었다. 표 3. trnh-gug와 psba primer로증폭한유전자의염기서열을 NCBI data base 에검색한결과 Description Query coverage Identities Pinus taeda 99% 100% Pinus taeda, psba-trnh 96% 100% Pinus rigida, psba gene 97% 99% Pinus latteri, psba gene 98% 99% Pinus latteri, psba gene 98% 99%... Pinus densiflora, psba gene 98% 97% NCBI에등록되어있는 P. densiflora의 trnh GUG -psba non-coding 염기서열 (NCBI 등록번호 : JN854210) 과 PCR을통해증폭된 PCR product의염기서열을 align 하여비교하였다. 그림 7에서확인할수있듯이붉은색화살표로표시한 11개의염기서열부위에서서로불일치하는것을확인할수있었다. 이렇게불일치하는서열의개체간차이를확인하기위하여 primer를디자인한 base data외에 NCBI에서추가로검색되는다른 P. densiflora 의 trnh GUG -psba 유전자 data base와 align 하여비교하였다. 그결과앞서 NCBI 에등록되어있는 trnh GUG -psba 유전자 base는모두일치하나분석된시료의 PCR product 유전자와는동일하게 11개의염기서열부위에서불일치함을확인할수있었다.( 그림8) 이와같은불일치부위의존재는본연구에서사용한소나무시료가 NCBI에등록되어있는소나무와분류학상어느수준에서의차이인지는 ( 종간의차이인지개체간의차이인지 ) 밝히기위해서는추가적인연구가필요하다. 본연구에서 non-coding 부위에서제작한 primer를이용하여유전자증폭을통해소나무의 trnh GUG -psba 유전자를확보할수있음을보여주었다. - 41 -

그림 7. NCBI 에등록된 P. densiflora 염기서열과 trnh-gug 와 psba primer 로 증폭한유전자의염기서열을 align 한결과 - 42 -

그림 8. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이외에 NCBI data base에검색되는 P. densiflora 의 trnh GUG -psba 유전자들을 align) - 43 -

2) matk 가 ) Primer 디자인 MaturaseK (matk) 유전자는식물색소체유전자로써 intron maturase 를암호화하는유전자부위이다. 따라서식물계통학적으로매우잘보존되어있다고알려져있다.[44] 이에본연구에서는그림 9와같이 P. densiflora 의 NCBI data base를바탕으로, matk primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서제안된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 matkf와 makr을각각 0.8 μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의용량으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 51 ~59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였지만 51 ~ 59 의모든온도조건에서유전자가증폭되지않았다. 따라서유전자염기서열분석이불가능하였다. - 44 -

그림 9. P. densiflora 의 matk 유전자서열과디자인한 primer 의위치 - 45 -

3) trns GCU - 5 trng 가 ) Primer 디자인 그림 10과같이 P. densiflora 의 NCBI data base를바탕으로, trns GCU - 5 trng coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trns-gcu와 5 trng를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 51 ~59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였지만 51 ~59 모든온도조건에서유전자가증폭되지않았다. 따라서유전자염기서열분석이불가능하였다. - 46 -

그림 10. P. densiflora 의 trns GCU - 5 trng 유전자서열과디자인한 primer 의위치 - 47 -

4) 5 trng 3 trng UUC 가 ) Primer 디자인 그림 11과같이 P. densiflora 의 NCBI data base를바탕으로, 5 trng 3 trng UUC coding region에 primer를디자인하려했으나 5 trng coding region에적절한 Tm값또는 GC content를갖는 primer가검색되지않았다. 따라서 5 trng exon 부근에서 non-coding region을포함하는 primer를디자인하였다이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trng F와 trng-uuc를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 57 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 12(A) 에서약 700bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 1차 PCR product를 template로하여 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 27회로변경하여진행하였다. 그림 12(B) 에서약 700bp에서증폭하고자하는유전자를확인할수있었고, 이를염기서열분석하였다. - 48 -

그림 11. P. densiflora의 5 trng 3 trng UUC 유전자서열과디자인한 primer의위치 - 49 -

그림 12. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer trng F 와 trng-uuc 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 50 -

다 ) 소나무 (Pinus densiflora) 의 5 trng 3 trng UUC 염기서열 증폭된 PCR product 를 trng F 와 trng-uuc primer 를이용하여양방 향염기서열분석을하였고, 그결과를 vector NTI 를통하여 align 한후에 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. ACGTCAGATCGTGGGCGGGTATAGTTTAGTGGTAAAAGTGTGATTCGTTACAT TATCAACTCGAATAGTTGAAGGATCTTTTACATGGATCTTTGATCCATCTAAAG CTCTATTTCCAGATAGGTTCAAAACCTTCCCTATGAATACCCAGGAATAGCATA CCTTCTCGTAATCTGGGTCTGAAATGATGAAAGAAAGATAAACACATTTTTGGT TCCAAGATAGCGACACAGAATGTTTGAAAGGTTAGATAAATTACAGATCTATG GGGAAGATATTTTTTCATCGTGACTAAACCTTCAAATTATGGATGAAAGGACCC CAGGTTGAAGCCGAATAGCTATAGAACTATGACTTTTGTTTACTTGGGTTATCG GCATTTCGTTCGAGCTCGGTGGAATTTGTTCTCCTAGGGATCAGAATCAGTAG AAATAGGGAACGAAGTAACTAGAAAGATTTGTGATTATTTGAAACTTTTCAAAT TGATCTTTCTATCAGGATCATCTAGAAAGTGAGTAAGTATTCTACGATTTTGGG CCCTTCACTAAAAAAAGAGAATAAACTGACGTAGATGCCATGGGTACGATAAG AATTTAGGGTTTTATTAGAAAAAAAAAAAAGGAAAAAGAAAGGAATTTTAAAT TCCTTTCAAAAAAATTTGATATAAATACTCCGCTTCTTCCCTCATACCGCTCTCT ATCTCCCATGAAGGAGCCAAATGACATGAAATTCTCAGGTTCGGTTCTGAATTA GAGACATTGAATAATCAATGTCAACTATAACCCCTAGCCTTCCAAAGCTAAGGA - 51 -

라 ) 소나무 (Pinus densiflora) 의 5 trng 3 trng UUC 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표4와같이 Pinus taeda가 query coverage 99%, identities 99% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 99%, identities 96% 로검색되었다. 표 4. trng F와 trng-uuc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 99% 99% Pinus clausa isolate CLAU02 99% 99% Pinus rigida isolate RIGI01 99% 99% Pinus pungens isolate PUNG01 99% 99% Pinus sabiniana isolate SABI04 99% 99%. Pinus densiflora isolate DENS02 99% 96% NCBI에등록되어있는 P. densiflora의 5 trng 3 trng UUC 유전자간 non-coding 염기서열 (NCBI 등록번호 :JN854210) 과앞서증폭된 PCR product의염기서열을 align하였다. 그림13에서확인할수있듯이화살표로표시한염기서열에서 NCBI data와 22개의부위에서염기서열이불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한유전자이외에 NCBI에추가로검색되는 P. densiflora의 5 trng 3 trng UUC 유전자 base data를검색하였으나 P. densiflora의다른개체의 5 trng 3 trng UUC 유전자부위는등록되지않은것을확인할수있었다. 따라서본연구를통해소나무 (P. densiflora) 의 5 trng 3 trng UUC 유전자가확보될수있음을보여주고이는소나무의분류학상차이 ( 종간, 개체간 ) 연구의하나의자료가될수있음을보여준다. - 52 -

그림 13. NCBI 에등록된 P. densiflora 염기서열과 trng F 와 trng-uuc primer 로증폭한유전자의염기서열을 align 한결과 - 53 -

5) rpob trnc GCA 가 ) Primer 디자인 그림14와같이 P. densiflora 의 NCBI data base(ncbi 등록번호 : JN854210) 를바탕으로, rpob trnc GCA coding region 에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하고염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora 로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 rpob와 trnc-gca를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 15(A) 에서약 500bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 27회로변경하여진행하였다. 그림 15(B) 에서약 500bp의크기의유전자를확인할수있었고, 밴드의농도가충분히높아염기서열분석을수행하였으나, 염기서열분석결과를얻지못하였다. 이는염기서열분석시염기서열상의 A와 T 염기가반복되는서열을가지고있어서 slippage현상이일어나염기서열을읽는데어려움이있는것으로판단된다. - 54 -

그림 14. P. densiflora의 rpob trnc GCA 유전자서열과디자인한 primer의위치 - 55 -

그림 15. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer rpob 와 trnc-gca 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 56 -

6) psbm trnd GUC 가 ) Primer 디자인 그림 16과같이 P. densiflora 의 NCBI data base를바탕으로, psbm trnd GUC coding region에대한 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 psbmf와 trnd-guc-r를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 17(A) 에서약 670bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 1000배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 33회로변경하여진행하였다. 그림 17(B) 에서약 670bp에서증폭하고자하는유전자를확인할수있었고, 밴드의농도가충분히진하기때문에염기서열을분석하였다. - 57 -

그림 16. P. densiflora 의 psbm trnd GUC 유전자서열과디자인한 primer 의위치 - 58 -

그림 17. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer psbmf 와 trnd-guc-r 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 59 -

다 ) 소나무 (Pinus densiflora) 의 psbm trnd GUC 염기서열 증폭된 PCR product 를 psbmf 와 trnd-guc-r primer 를이용하여양 방향으로염기서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. ACTGGCCGTTTTGACGTAAGGGATAAGTAGAAAAGCAGTAGGAACTGCAAGGA ACAGTAGAACAGCAATAAATGCAAGGGTATTTACTTCCATGATCTCCTTATTTT AACTTTGTTCAAAAATAGCTCGAGATTTGATCTCATTTTGATCTCATTTTGATC TCATAGAGATGAATGAATCTTTCTTTCAAGATCCATGAAATAGATGTATGTCAT TGATAGAATTGGTTCGAGTAACGGAATCTAACTAATGGATTGAATGAATTCTT CGATGGAAATAGTATAACCTAATACGATCACTTTTGATAAGACTAGTAGTAAAA ACTTACGTGCATCAGAAAACGTTCCGATCAACACATGTCTGGTATAATTTCGAA AGAATAGGATTCGTACGAATAAGAACCTTATGCTCCTCCAGAAGACGTTCGTC AGACATGAGAGATATTTTGCTTGGATCTCCACGATTTAGATGGAATTGTGAAT CTATCCCGCTTCGGTGATGATATAGAAAGAAGTATCCCATATCGAATTTATCCA GAGGCCCACCCATGACCCTGGAATGATAAGGACTAGTCCGTCCAGATCCTGAC ATGATCATTCACAGTTCACTTACTATTGGATCGGGACTGACGGGACTCGAACCC GCAACTTCCGTCTTGACAGGGCGGTACTCT - 60 -

라 ) 소나무 (Pinus densiflora) 의 psbm trnd GUC 염기서열분석 얻어진염기서열을 NCBI 의 nucleotide blast 검색을하였고, 표 5 와같 이 P. densiflora 가 query coverage 100%, identities 100% 로검색되었다. 표 5. psbmf와 trnd-guc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus sylvestris, complete genome 100% 100% Pinus taiwanensis, complete genome 100% 100% Pinus massoniana, complete genome 100% 100% Pinus densiflora isolate DENS02 100% 100% Pinus massoniana isolate MASS02 100% 100% NCBI에등록되어있는 P. densiflora의 psbm trnd GUC non-coding 염기서열과증폭된유전자의염기서열을 align하였다. 그림 18에서확인할수있듯이모든서열이일치하였다. 따라서본연구를통해소나무 (P. densiflora) 의 psbm trnd GUC non-coding 유전자가확보될수있음을보여주고, 소나무시료에대한종식별후보유전자가될수있음을확인하였다. - 61 -

그림 18. NCBI 에등록된 P. densiflora 염기서열과 psbmf 과 trnd-guc primer 로증폭한유전자의염기서열을 align 한결과 - 62 -

7) trnd GUC trne UUC 가 ) Primer 디자인 그림 19와같이 P. densiflora 의 NCBI data base를바탕으로, trnd GUC trne UUC coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하고염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnd-guc, trne-uuc를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 20(A) 에서약 470bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 10,000배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 32회로변경하여진행하였다. 그림 20(B) 에서약 470bp에서증폭하고자하는유전자를확인할수있었고, 이유전자의염기서열을분석하였다. - 63 -

그림 19. P. densiflora 의 trnd GUC trne UUC 유전자서열과디자인한 primer 의위치 - 64 -

그림 20. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer trnd-guc-f 와 trne-uuc 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 65 -

다 ) 소나무 (Pinus densiflora) 의 trnd GUC trne UUC 염기서열 PCR 로증폭된유전자를 trnd-guc 와 trne-uuc primer 를이용하여 양방향으로염기서열을분석하였다. 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TTCTCTTTCAAGGAGGCAACGGGGATTCAACTTCCCCTGGGGGTACCATGATC CATTCGTTAGGTATCCTAAATAATTTTCAATTACTGTCTTGTTCCTGGGTCGAT GCCCGAGTGGTTAATGAGGACGGACTGTAAATCCGTTGGCAATATGCCTACGC TGGTTCAAATCCAGCTCGACCCAACCAATATCATCAATACCAATCTATCGGTAT GGTTATATTCATGCTAATCATGAATGAAAAGATAATTGGTTATTGAACCCCGGC ATCGATATCTATTCATATCGTAGATGCCCGATTCCGTATGAATCGATACATTTC AAAAATGAAAGAGAAGATTAAGATATTTAATCGACCTAGCACTCGCATAATCTA TATGACCTATCTAGCACCTATCATGGAATAAATATTATCCAATAGTAGGTGAAC TGTACTGTATTGGGATTGTAGCTCAATTGGTTAGAGTACCGCCCT 라 ) 소나무 (Pinus densiflora) 의 trnd GUC trne UUC 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표6과같이 Pinus taeda가 query coverage 99%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 99%, identities 99% 로검색되었다 표 6. trnd-guc와 trne-uuc primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 99% 100% Pinus cubensis isolate CUBE01 99% 100% Pinus echinata isolate ECHI01 99% 100% Pinus elliottii isolate ELLI01 99% 100% Pinus densiflora isolate DENS02 99% 99% - 66 -

NCBI에등록되어있는 P. densiflora 의 trnd GUC trne UUC non-coding 염기서열 (NCBI 등록번호 : JN854210) 과증폭된유전자의염기서열을 align하였다. 그림 21에서확인할수있듯이화살표로표시한염기서열에서 NCBI data 5개의부위에서염기서열이불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한 base data외에 NCBI에추가로검색되는다른 P. densiflora의 trnd GUC - trne UUC 유전자 data base를검색하였으나다른개체의 trnd GUC - trne UUC 유전자부위는등록되지않은것을확인할수있었다. 따라서본연구를통해소나무 (P. densiflora) 의 trnd GUC - trne UUC 유전자가확보될수있음을보여주고이는소나무의분류학상차이 ( 종간, 개체간 ) 연구의하나의 base data 가될수있음을보여준다. 그림 21. NCBI 에등록된 P. densiflora 염기서열과 trnd-guc 와 trne-uuc primer 로증폭한유전자의염기서열을 align 한결과 - 67 -

8) trnd GUC trny GUA 가 ) Primer 디자인 그림 22과같이 P. densiflora의 NCBI data base를바탕으로, trnd GUC trny GUA coding region에서 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하고그염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnd-guc와 trny-gua를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1 μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 57 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 23(A) 에서약 350bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 37회로변경하여진행하였다. 그림 23(B) 에서약 350bp에서증폭하고자하는유전자를확인할수있었고, 밴드의농도가충분히진하기때문에염기서열을분석하였다. - 68 -

그림 22. P. densiflora 의 trnd GUC trny GUA 유전자염기서열과디자인한 primer 의위치 - 69 -

그림 23. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer trnd-guc-f 와 trny-gua 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 70 -

다 ) 소나무 (Pinus densiflora) 의 trnd GUC trny GUA 염기서열 증폭된유전자를 trnd-guc 와 trny-gua primer 를이용하여양방향으 로염기서열분석을하였다. 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TCAGGGCGGTACTCTAACCAATTGAGCTACAATCCCAATACAGTACAGTTCACC TACTATTGGATAATATTTATTCCATGATAGGTGCTAGATAGGTCATATAGATTA TGCGAGTGCTAGGTCGATTAAATATCTTAATCTTCTCTTTCATTTTTGAAATGT ATCGATTCATACGGAATCGGGCATCTACGATATGAATAGATATCGATGCCGGG GTTCAATAACCAATTATCTTTTCATTCATGATTAGCATGAATATAACCATACCG ATAGATTGGTATTGATGATATTGGTTGGGTCGAGCTGGATTTGAACCAGCGTA GGCATATTGCCAACGGATTTACAGTCCGTCCTCATTAACCACTCGGGCATATG 라 ) 소나무 (Pinus densiflora) 의 trnd GUC trny GUA 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표 7과같이 Pinus taeda가 query coverage 98%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 98%, identities 98% 로검색되었다. 표 7. trnd-guc와 trny-gua primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 98% 100% Pinus cubensis isolate CUBE01 98% 100% Pinus echinata isolate ECHI01 98% 100% Pinus elliottii isolate ELLI01 98% 100% Pinus greggii isolate GREG02 98% 100% Pinus densiflora isolate DENS02 98% 98% - 71 -

NCBI에등록되어있는 P. densiflora의 trnd GUC trny GUA non-coding 염기서열 (NCBI 등록번호 : JN854210) 과앞서증폭된유전자의염기서열을 align하였다. 그림 24에서확인할수있듯이화살표로표시한 4개의염기서열부위에서불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한 base data외에 NCBI에추가로검색되는다른 P. densiflora 의 ttrnd GUC trny GUA 유전자 base data 를검색하였으나 P. densiflora의다른개체의 trnd GUC trny GUA 유전자부위는등록되지않은것을확인할수있었다. 따라서본연구를통해소나무 (P. densiflora) 의 trnd GUC trny GUA 유전자가확보될수있음을보여주고이는소나무의분류학상차이 ( 종간, 개체간 ) 연구의하나의 base data가될수있음을보여준다. 그림 24. NCBI 에등록된 P. densiflora 염기서열과 trnd-guc 와 trny-gua primer 로증폭한유전자의염기서열을 align 한결과 - 72 -

9) trnfm CAU - trns UGA 가 ) Primer 디자인 그림 25과같이 P. densiflora의 NCBI data base를바탕으로, trnfm CAU - trns UGA coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하고염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnfm-cau와 trns-uga를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1 μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 26(A) 에서약 950bp에서증폭하고자하는유전자의위치를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 27회로변경하여진행하였다. 그림 26(B) 에서약 950bp에서증폭하고자하는유전자를확인할수있었고, 유전자의농도가충분히진해서염기서열을분석하였다. - 73 -

그림 25. P. densiflora의 trnfm CAU - trns UGA 유전자서열과디자인한 primer의위치 - 74 -

그림 26. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer trnfm-cau 와 trns-uga 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 75 -

다 ) 소나무 (Pinus densiflora) 의 trnfm CAU - trns UGA 염기서열 증폭된유전자를 trnfm-cau 와 trns-uga primer 를이용하여양방향 으로염기서열을분석하였다. 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TCAATTTATAATTGAATTATAAGTGAGGAAAAGAAACGATCAATACATGAACTA TTCTTTTTCATATTCTATGTGAAGGGCGGGTAGCGGGGATCGAACCCGCATCT TCTCCTTGGCAAAGAGAAATTTTACCATTCAACCATACCCGCATTTCAATTATA TGAATATATATATATTCATATAATTGAAATGGAAAATACATATATGTTCAATCC CATTCGTAAGTAGATACCATTTTTGAATGGTCCAATTATTAATTAAATTACGGA TATGGAAAACCCCTCCTCCTTGGGCCTGAAGTAAAAGGCCCTTCAGAAAAGCT ATAAAGCCCTCCGGGAGGAGGGGGAGGGAGTGTGAACCTAAAACATCTAGAAC ATATCTAAGATATGAAAGAATTAAGGATACCTACAAAAAGGACTGATCCGATC CATAATGATACACCCGAAAATACAACATTTTTGCTACTTGACCAACCATCGGGA GAGGCAAGTGCAACAGGTACACCAATCACTAGTAAAAATGAAATAGCAATTAA TGCAAACACAGCCGATTGGAAAGCAATAGTCATGGTTGTGATCTTACAAGCTC CCAACAGATTGAATAGACTATACCATCCGATCCCCAATTTTCAATTCTGATCAA ATGCACTATGGAAAGGAAAGGATTTGACCAGAAATTGATCGAGCTTTCAAACT TTTTCAATTTGATATTTGAGTGTGAGAGGAGAGGGAAATTCGTTTCTTTTTTTT TCCATTCCAGATGATGAAAAATCATCATATATATGTCACAGGTATAGTTGGTAT CGACCCGA - 76 -

라 ) 소나무 (Pinus densiflora) 의 trnfm CAU - trns UGA 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표 8와같이 Pinus taeda가 query coverage 100%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 100%, identities 98% 로검색되었다 표 8. trnfm-cau와 trns-uga primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 100% 100% Pinus pungens isolate PUNG01 100% 100% Pinus serotina isolate SERO01 100% 100% Pinus cubensis isolate CUBE01 100% 99% Pinus densiflora isolate DENS02 100% 98% NCBI에등록되어있는 P. densiflora의 trnfm CAU - trns UGA non-coding 염기서열 (NCBI 등록번호 : JN854210) 과증폭된유전자의염기서열을 align하였다. 그림 27에서확인할수있듯이화살표로표시한 8개의염기서열부위에서 NCBI에등록되어있는유전자와불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한유전자이외에 NCBI에추가로검색되는 P. densiflora의 trnfm CAU - trns UGA 유전를검색하였으나, P. densiflora의다른개체의 trnfm CAU - trns UGA 유전자부위는등록되지않은것을확인할수있었다. 따라서본연구를통해소나무 (P. densiflora) 의 trnfm CAU- trns UGA 유전자가확보될수있음을보여주고이는소나무의분류학상차이 ( 종간, 개체간 ) 연구의하나의 base data가될수있음을보여준다. - 77 -

그림 27. NCBI 에등록된 P. densiflora 염기서열과 trnfm-cau 와 trns-uga primer 로증폭한유전자의염기서열을 align 한결과 - 78 -

10) 5 trnl UAA - trnt UGU 가 ) Primer 디자인 5 trnl UAA - trnt UGU 부위는목재의엽록체 DNA 중에서보편적으로수종을식별하는데이용되는부위로보고되어있다. 그림 28과같이 P. densiflora 의 NCBI data base(ncbi 등록번호 : JN854210) 를바탕으로, 5 trnl UAA - trnt UGU coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 5 trnl-uaa와 trnt-ugu 를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 57 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 29(A) 에서약 500bp에서유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된유전자를 1,000배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 27회로변경하여진행하였다. 그림 29(B) 에서약 500bp에서증폭하고자하는유전자를확인할수있었고, 이유전자의염기서열을분석하였다. - 79 -

그림 28. P. densiflora 의 5 trnl UAA\ - trnt UGU 유전자서열과디자인한 primer 의위치 - 80 -

그림 29. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer 5 trnl-uaa 와 trnt-ugu 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 81 -

다 ) 소나무 (Pinus densiflora) 의 5 trnl UAA trnt UGU 염기서열 증폭된유전자를 trnl-uaa 와 trnt-ugu primer 를이용하여양방향으 로염기서열을분석하였고, 그결과를 vector NTI 를통하여 align 한후에 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TTGCGATGCTCTAACCTCTGAGCTAAGCAGGCTCGATGGAATATCATTTCCCAT TTCATTGGGGTGAGATCCATAGATTCCTTTGGAATCCTAACGATTATAGCGCG AATCAGATTCAATGACGCAATCCAGATTACGATTACAACGGAACATTGTCTGAA TGTAGATTCTTTCAAGGGAAAGAAAGGAGAAGTAAGGGTCAATTGATATCGAT CGTTTTACTCACTCTTCCAAATCGACTAGGGGAGGATAATAACATTGCATTTCA AATGGATAAATAATATAATGATTCCAGTCATATTCGATTGGGGTAGAGATAGA GAAGGGGAGAGATGGGTAGAGAAGGGGAGAGAAGGGGAGAGAAGGGGAGTA GGGGAGGATATTTCTCCCACCCAATATTGAACAAATATCCAATGAATAATGCTG ATGG 라 ) 소나무 (Pinus densiflora) 의 5 trnl UAA trnt UGU 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표 9와같이 Pinus pungens가 query coverage 99%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 99%, identities 97% 로검색되었다. 표 9. 5 trnl-uaa와 trnt-ugu primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Description Query coverage Identities Pinus pungens isolate PUNG01 99% 100% Pinus rigida isolate RIGI01 99% 100% Pinus serotina isolate SERO01 99% 100% Pinus taeda, complete genome 99% 99% Pinus densiflora isolate DENS02 99% 97% - 82 -

NCBI에등록되어있는 P. densiflora의 5 trnl UAA - trnt UGU non-coding 염기서열 (NCBI 등록번호 : JN854210) 과앞서증폭된유전자의염기서열을 align하였다. 그림 30에서확인할수있듯이화살표로표시한 8개의염기서열부위에서 NCBI에등록되어있는유전자와불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한 data base이외에 NCBI에추가로검색되는 P. densiflora의 5 trnl UAA - trnt UGU 유전자 (JX504689.1, JX504684.1) 와 align하였다. 그결과앞서 NCBI에등록되어있는 5 trnl UAA trnt UGU 유전자 base는모두일치하나분석된시료의유전자는동일하게 8개의염기서열부위에서불일치함을확인할수있었다.( 그림31) 이와같은불일치부위의존재는시료가 NCBI 에등록되어있는소나무와분류학상어느수준에서의차이인지는 ( 종간의차이인지개체간의차이인지 ) 추가적인연구를통해밝혀야될것이다. 본연구를통해소나무 (P. densiflora) 의 5 trnl UAA - trnt UGU 유전자가확보될수있음을보여주었다. 그림 30. NCBI 에등록된 P. densiflora 염기서열과 5 trnl-uaa 와 trnt-ugu primer 로증폭한유전자의염기서열을 align 한결과 - 83 -

그림 31. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이외에 NCBI data base에검색되는 P. densiflora의 5 trnl UAA trnt UGU 유전자의 align) - 84 -

11) 3 trnl UAA - trnl5 UAA 가 ) Primer 디자인 그림 32와같이 P. densiflora의 NCBI data base를바탕으로, 3 trnl UAA trnl5 UAA coding region 에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 3 trnl UAA 와 trnl5 UAA 를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1 μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 51 ~59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였지만 51 ~59 모든온도조건에서유전자가증폭되지않았다. 따라서유전자염기서열분석이불가능하였다. - 85 -

그림 32. P. densiflora의 3 trnl UAA - trnl5 UAA 유전자서열과디자인한 primer의위치 - 86 -

12) 5 rps12 - rpl20 가 ) Primer 디자인 그림 33과같이 P. densiflora의 NCBI data base를바탕으로, 5 rps12 - rpl20와 coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하고유전자의염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 rps12와 rpl20를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 34 (A) 에서약 800bp에서유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 37회로변경하여진행하였다. 그림 34 (B) 에서약 800bp에서증폭하고자하는유전자를확인할수있었고, 밴드의농도가진하기때문에유전자의염기서열을분석하였다. - 87 -

그림 33. P. densiflora 의 5 rps12 - rpl20 유전자서열과디자인한 primer 의위치 - 88 -

그림 34. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer rps12 와 rpl20 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 89 -

다 ) 소나무 (Pinus densiflora) 의 5 rps12 rpl20 염기서열 증폭된유전자를 rps12 와 rpl20 primer 를이용하여양방향으로염기 서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TCTCCTGCTCTTCGAGGATGTCCTCAGCGTAGAGGAGTATGTGCTAGGGTGTA TGTGCGACTCGTTTGGATCAGAAACTTAAACAGACTGGGAAATTCTTACAACG GAATAACAGAAGAATTTTCCCGCTGTAATCAAGATCCATCATCTAACCTTACCG GGGATGAAATTTATGGTTTCCATTGGTGCAAATCCAATCACCTTGATGTGGGA TGAAAAGCAATTCCTCATTGGTAGCGAATGGTCATCAATCCATTGAGCGGGGA AATCATGCAAATTGAAGAAGCATAAAGTTTATGCTCATATTGCCGCGAGAACG GAACAACAGGGTCAGCTACCTGGCCAACCCCAGAATTACATGTCGTTACTGTAT TAATAGATCTTGTAATGAGAGTATTTATTACTTAATGATTCAAGAGTAGAGCTG GAGATGGTAAACCATACAAGTTACCGATAACATACTCATCTCATATTTCGGAAA TGAATAAAAGGCTCCGGCGTATAGAGAGGACCTTACCGTTGGAGAAAGAACCA TAGAAACGATGAAACCCATGATTTTTTCCCATTCTCAGTACAAGGTCCCAGTCC TTGCCTACCAGGAAGATGCCTATCCAAAGGGAATTATGGTTGGGAAGGTTGCA GTAGCAAAAGCCATTGGAACTTTTATTTTCTAAATAAGAAGAATCAGTGTTATT CTCAATGGACCAAACAAAACAAATGACTAACCGAGAAAAATATTAGCGATTCAT GAGAGTAAACTTCAATGACCAGAGTGAAACGTGGATATATAGCTCGAAAACGT CGGAAAA - 90 -

라 ) 소나무 (Pinus densiflora) 의 5 rps12 rpl20 염기서열분석 얻어진염기서열을이용하여 NCBI의 nucleotide blast 검색을하였고, 표 10와같이 Pinus taeda가 query coverage 100%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 100%, identities 98% 로검색되었다. 표 10. 5 rps12 와 rpl20 primer 로증폭한유전자의염기서열을 NCBI data base 에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 100% 100% Pinus cubensis isolate CUBE01 100% 100% Pinus echinata isolate ECHI01 100% 100% Pinus elliottii isolate ELLI01 100% 100% Pinus pringlei isolate PRIN02 100% 100% Pinus densiflora isolate DENS02 100% 98% NCBI에등록되어있는 P. densiflora의 5 rps12 - rpl20 non-coding 염기서열 (NCBI 등록번호 : JN854210) 과앞서증폭된유전자의염기서열을 align하였다. 그림 35에서확인할수있듯이 6개의화살표로표시한염기서열부위에서 NCBI에등록된염기서열과불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한 base data외에 NCBI에추가로검색되는 P. densiflora의 5 rps12 - rpl20 유전자 base data를검색하였으나 P. densiflora의다른개체의 5 rps12 - rpl20 유전자부위는등록되지않은것을확인할수있었다. 따라서본연구를통해소나무 (P. densiflora) 의 5 rps12 - rpl20 유전자가확보될수있음을보여주고이는소나무의분류학상차이 ( 종간, 개체간 ) 연구의하나의자료가될수있음을보여준다. - 91 -

그림 35. NCBI 에등록된 P. densiflora 염기서열과 5 rps12 와 rpl20 primer 로 증폭한유전자의염기서열을 align 한결과 - 92 -

13) psbb - psbh 가 ) Primer 디자인 본연구에서는그림 36과같이 P. densiflora의 NCBI data base를바탕으로, psbb psbh coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 psbb와 psbh를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 51 ~59 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였지만 51 ~59 모든온도조건에서유전자가증폭되지않았다. 따라서유전자염기서열분석이불가능하였다. - 93 -

그림 36. P. densiflora 의 psbb psbh 유전자서열과디자인한 primer 의위치 - 94 -

14) rpl16 가 ) Primer 디자인 본연구에서는그림 37과같이 P. densiflora의 NCBI data base를바탕으로, rpl16 유전자 coding 부위에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 rpl16f와 rpl16r을각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 38(A) 에서약 900bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 100배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 25회로변경하여진행하였다. 그림 38(B) 에서약 900bp에서증폭하고자하는유전자를확인할수있었고, 밴드의농도가충분히진해염기서열을분석하였다. - 95 -

그림 37. P. densiflora 의 rpl16 염기서열과디자인한 primer 의위치 - 96 -

그림 38. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer rpl16f 와 rpl16r 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 97 -

다 ) 소나무 (Pinus densiflora) 의 rpl16 염기서열 증폭된유전자를 rpl16f 와 rpl16r primer 를이용하여양방향으로염 기서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였 다. 이를통해얻어진유전자의염기서열은다음과같다. GATTGGTCTCTATAAATAAACTAGAGACAACTCATTGCTTCATATTGCCAAGAT CCAATAACTATGGGTAGATACTATCACCTCGTAAATACTTTCTTCCACGAGATA AAAAATGATATGATATTTTGATCTCTCGTGAAGCGAGAAAGGTGTTTGGTAAT GCGGAAAAAGAAGAAAAAGGAGAAATATTCTCCCAATATTGTATGGTTCATTA ACTACCCTCCAAAAAGCAGTGTGATAAAGCATAAGAATCATCCTGATTCATTCA AGCAAGATTGAAGGAATTCAGTTTATGAAGGAGAAAGAGCTTCGGGTCGATGG AAACTAAGGAAATTGATTCCGTAACAGATGAAAAGAAAGCTCCATTGCGGAGG GAAGACCTGGGCATTTAGATAGAAGCTATGGAACGATGGAACCTATGACTGCA TAAGATCATATAGGAATCTTAATCATTCATTGGATAGGATGGCGAAATAAACCA AAAACCAATTAATATCATTGGAGTCTATAGGTAAGTCGACCCCATGGAGCAAA TACTGGAAGAGTCAATATTCGCCTGTGAAATTCTTTATTAAGACATTGGATGGA AATGTAAGAGTTATTCGTCCTGTAAATTATATTCTATATACAATATGTGTAATG CATAGATATAGACTCATTTATCTACACATTGATTATTCACGAGGAGCCGGATGA GAAGAAACTTTCATGTCCGGTTCTGGATTAGAGATGGGATCAAAATACTATTA ACCATCGACTATAACCCAAAAAGAACAAAA - 98 -

라 ) 소나무 (Pinus densiflora) 의 rpl16 염기서열분석 얻어진염기서열을 NCBI의 nucleotide blast 검색을하였고, 표 11과같이 Pinus taeda가 query coverage 100%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 100%, identities 98% 로검색되었다. 표 11. rpl16f 와 rpl16r primer 로증폭한유전자의염기서열을 NCBI data base 에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 100% 100% Pinus coulteri isolate COUL03 100% 100% Pinus cubensis isolate CUBE01 100% 100% Pinus devoniana isolate DEVO02 100% 100% Pinus hartwegii isolate DONN02 100% 100% Pinus densiflora isolate DENS02 100% 98% NCBI에등록되어있는 P. densiflora의 rpl16 염기서열 (NCBI 등록번호 : JN854210) 과앞서증폭된유전자의염기서열을 align하였다. 그림 39에서확인할수있듯이화살표로표시한 8개의염기서열부위에서 NCBI에등록된유전자와불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한 base data외에 NCBI에추가로검색되는 P. densiflora의 rpl16 유전자 (JF701594.1) 와 align하였다. 그결과앞서 NCBI에등록되어있는 rpl16 유전자 base는모두일치하나분석된시료의유전자는동일하게 8개의염기서열부위에서불일치함을확인할수있었다.( 그림40) 이와같은불일치부위의존재는시료가 NCBI에등록되어있는소나무와분류학상어느수준에서의차이인지는 ( 종간의차이인지개체간의차이인지 ) 추가적인연구를통해밝혀야될것이다. 본연구를통해소나무 (P. densiflora) 의 rpl16 유전자가확보될수있음을보여주었다. - 99 -

그림 39. NCBI 에등록된 P. densiflora 염기서열과 rpl16f 와 rpl16r primer 로 얻은유전자의염기서열을 align 한결과 - 100 -

그림 40. NCBI 에등록된여러개의 P. densiflora 염기서열과본실험에서얻 은유전자의염기서열을 align 한결과 (primer 를디자인한보고된유전자이 외에 NCBI data base 에검색되는 P. densiflora 의 rrpl16 유전자의 align) - 101 -

15) rpl20 - rps18 가 ) Primer 디자인 본연구에서는그림 41과같이 P. densiflora의 NCBI data base를바탕으로, rpl20 - rps18 coding region에 primer를디자인하였으며, 이를이용하여 non-coding region을증폭하여염기서열을분석하였다. primer에대한정보는표 2에정리하였다. 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열분석을하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2 μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 rpl20-rps18f와 rpl20-rps18r를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 35회반복하였다. 그후 7 2 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 42(A) 에서약 500bp에서증폭하고자하는유전자를확인할수있었으나그농도가매우약하기때문에유전자서열분석이어렵다고판단되어 2차 PCR을진행하였다. 1차 PCR 조건과동일하되증폭된 PCR product 를 10,000배희석하여 template로사용하고 denaturation 단계부터 extension 단계를 35회에서 32회로변경하여진행하였다. 그림 42(B) 에서약 500bp에서증폭하고자하는유전자를확인할수있었고, 밴드의농도가충분히진해염기서열을분석하였다. - 102 -

그림 41. P. densiflora 의 rpl20 - rps18 유전자서열과디자인한 primer 의위치 - 103 -

그림 42. P. densiflora 시료로부터얻은 DNA 를 template 로하여 primer rpl20-rps18f 와 rpl20-rps18r 로 1 차 PCR 한결과 (A) 와 2 차 PCR 한결과 (B) - 104 -

다 ) 소나무 (Pinus densiflora) 의 rpl20 - rps18 염기서열 증폭된유전자를 rpl20-rps18f 와 rpl20-rps18r primer 를이용하여양 방향으로연기서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. TTTATCCGGCCGAGTGAATAGATTAACTTCAAAACAACAACGTCTAATGACTAA CGCTATTAAACGAGCTCGTATTCTATCTTTGTTACCTTTTCTTTATAATGAAAA CTAATTTGATCCATGGGAAATAAGCCTACTCCTTAATCAAATCGGAGCTGGAAT ATCTGTTCGATAAAGATGCAGATCTTGAGTCTATTGTTGAAAAAGTCAACAGG ATTTCATTTTTGGAAAAGAATTGGATTGAATTCTTTTCGTTCATTTCCAATCCA ATATGAAAAAACTTTTCAATATTTCGACCTTCCCGGAGGGAATTCTCCGGGAGA ATCTGTTCTATCCATTCCATCTTTACCTATTTATTTCATCTATACCTAACCTATT TCATCATACGATAAGTTGTTCAAGATGGTGGAAAAGCAATTACTATCTAATATA GCTATTTGTGCAAGTGTTTTACGATTTGG 라 ) 소나무 (Pinus densiflora) 의 rpl20 - rps18 염기서열분석 얻어진염기서열을이용하여 NCBI의 nucleotide blast 검색을하였고, 표 12와같이 Pinus taeda가 query coverage 99%, identities 100% 로검색되었다. 한편연구에사용된샘플인 P. densiflora는 query coverage 99%, identities 99% 로검색되었다. 표 12. rpl20-rps18f 와 rpl20-rps18r primer 로증폭한유전자의염기서열을 NCBI data base 에검색한결과 Description Query coverage Identities Pinus taeda, complete genome 99% 100% Pinus arizonica isolate COOP01 99% 100% Pinus coulteri isolate COUL03 99% 100% Pinus cubensis isolate CUBE01 99% 100% Pinus devoniana isolate DEVO02 99% 100% Pinus densiflora isolate DENS02 99% 99% - 105 -

NCBI에등록되어있는 P. densiflora의 rpl20 - rps18 non-coding 염기서열 (NCBI 등록번호 : JN854210) 과앞서증폭된유전자의염기서열을 align하였다. 그림 43에서확인할수있듯이화살표로표시한 4개의염기서열부위에서 NCBI에등록된유전자와불일치되는것을확인할수있었다. 이렇게불일치되는서열의개체간차이를확인하기위하여 primer를디자인한유전자이외에 NCBI에추가로검색되는 P. densiflora의 rpl20 - rps18 유전자 (AB019925.1) 와 align하였다. 그결과앞서 NCBI에등록되어있는 rpl20 - rps18 non-coding 유전자 base는모두일치하나분석된시료의유전자는동일하게 4개의염기서열부위에서불일치함을확인할수있었다.( 그림 44) 이와같은불일치부위의존재는시료가 NCBI에등록되어있는소나무와분류학상어느수준에서의차이인지는 ( 종간의차이인지개체간의차이인지 ) 추가적인연구를통해밝혀야될것이다. 본연구를통해소나무 (P. densiflora) 의 rpl20 - rps18 non-coding 유전자가확보될수있음을보여주었다. 그림 43. NCBI 에등록된 P. densiflora 염기서열과 rpl20-rps18f 와 rpl20-rps18r primer 로얻은유전자의염기서열을 align 한결과 - 106 -

그림 44. NCBI에등록된여러개의 P. densiflora 염기서열과본실험에서얻은유전자의염기서열을 align한결과 (primer를디자인한보고된유전자이외에 NCBI data base에검색되는 P. densiflora의 rpl20 - rps18 non-coding 유전자의 align) - 107 -

16) 수종분석을위한다양한유전자분석 소나무 (P. densiflora) 를이용하여다양한유전자부위를목재시료에서확보할수있는지실증하였다. NCBI에등록되어있는 P. densiflora의엽록체전체 DNA 정보를바탕으로참고문헌에서제안된 non-coding 부위의증폭을위한 15개의 primer 쌍을제작하였다. 이를이용하여중합효소연쇄반응을통한유전자증폭으로총 11개의유전자를확보하였다. 확보한분석유전자부위는특정 non-coding 부위에국한되어증폭된것이아니라다양한 non-coding 부위에서증폭되었다. 이는효과적으로다양한엽록체유전자를증폭하는데이용될수있는 DNA 추출이이루어지고있다고평가할수있다. 증폭된유전자부위의길이는최소 354bp에서최대 956bp로 1kb 길이의유전자에서는길이에따른증폭가능여부에큰영향을주지않는것으로판단된다. 이는 1kb 정도크기의 DNA들이본연구과제에서새로이제안하는 DNA 추출방법에의해안정적으로추출되고있음을보여준다. 또한선행연구에서확립한건조된목재로부터 DNA를추출하는기본매뉴얼의재현성을높음을확인하였다. 증폭된 11개의분석유전자의염기서열을분석한결과, 10개의유전자염기서열이분석가능하였다. 염기서열분석이불가능한 rpob-trnc GCA 부위는 A와 T 염기가반복되는서열을가지고있어서염기서열분석시에 slippage현상, primer binding 으로인한문제인것으로판단된다. 분석이이루어진 10개의염기서열을 NCBI blast 검색을통하여수종식별을진행하였다. 수종분석결과 8개는 Pinus taeda로검색되었고, 1개는 Pinus pungens, psbm-trnd GUC primer를이용한나머지 1개에대해서는 P. densiflora로검색되었다. 선행과제 ( 용역과제계약번호 : 20170205951-00) 에서 P. densiflora의수종분석을위한 primer로써 pdest와 cp1을이용하였으며, 추가적으로가능성이있는후보유전자부위를확보할수있었다. 그러나후보가목재식별을위한지표유전자가되기위해서는많은대표개체에대한추가실험이필요하다고사료된다. Pinus taeda와 Pinus pungens는현미경을이용하여해부학적분석시 P. densiflora와비교하여세포의구조가구분이잘되는수종이다. 이러한결과는현재까지 NCBI 에등록된수종별엽록체염기서열의데이터베이스가아직충분히구축되지않았기때문에, 어떠한수종인지판단하기에는그정확도가매우낮으며신뢰도가떨어진다고판단된다. 또한본연구에서진행한분석된목재시료의수가제한적이기때문에, 분석유전자부위의염기서열분석결과가개체간의차이인지또는종 - 108 -

간의차이인지확인하기어려운것으로판단된다. 본연구에서는엽록체의다양한부위에서유전의염기서열확보가가능하며향후도출되는지표유전자의선정에특정엽록체부위에대한제약이없음을증명하였다. - 109 -

표 13. 각 primer에따른유전자증폭여부및염기서열분석여부, NCBI blast 검색결과요약 연번 분석유전자 유전자유전자 Database search Query 증폭여부분석여부결과 coverage Identities 1 trnh GUG - psba O O Pinus taeda 99% 100% 2 matk X - - 3 trns GCU - 5'trnG X - - 4 5'trnG - 3'trnG UUC O O Pinus taeda 99% 99% 5 rpob - trnc GCA O X - 6 psbm - trnd GUC O O P. densiflora 100% 100% 7 trnd GUC - trne UUC O O Pinus taeda 99% 100% 8 trnd GUC - trny GUA O O Pinus taeda 98% 100% 9 trnfm CAU - trns UGA O O Pinus taeda 100% 100% 10 5'trnL UAA (TabB) - trnt UGU (TabA) O O Pinus pungens 99% 100% 11 3'trnL UAA (TabD) - trnl5' UAA (TabC) X - - 12 5'rpS12 - rpl20 O O Pinus taeda 100% 100% 13 psbb - psbh X - - 14 rpl16 O O Pinus taeda 100% 100% 15 rpl20 - rps18 O O Pinus taeda 99% 100% 선행연구 pdest cp1 O O P. densiflora 99% 100% - 110 -

다. 다양한수종에서의 DNA 안정적추출검증및염기서열분석 선행연구에서여러가지 non-coding 부위의분석유전자를이용하여 P. densiflora의특정유전자를증폭시켜 NCBI에등록된염기서열과비교, 분석한결과, 증폭되어염기서열을확보하여도 NCBI에등록된자료가미흡하여정확한비교, 수종분석에어려움이있었다. 그러나다양한 non-coding 부위를분석유전자롤사용하여도전체유전자로부터 PCR을통한선택적증폭이가능한것을확인하였다. 이에본연구에서는국내유통되는여러수종중, 해부학적분석만으로는정확한식별이불가능한대표적인수종을다음표 14와같이선별하였고, 그림에해부학적식별이어려운수종의현미경사진을비교해놓았다.( 한국임업진흥원제공 ) 표 14. 해부학적식별의한계가있는유통수종선별 수종 ( ) 침엽수 활엽수 Pinus densiflora Pinus palustris Chamaecyparis obtusa Eucalyptus diversicolor Acer saccharum Pinus resinosa Pinus rigida Thuja plicta Eucalyptus marginata Acer saccharinum - 111 -

그림 45. 세포의종류와구조가유사하여현미경상의식별이어려운수종의 현미경사진 ( 방사단면 ) (A) P. palustris (B) P. rigida 그림 46. 세포의종류와구조가유사하여현미경상의식별이어려운수종의 현미경사진 ( 횡단면 ) (A) E. diversicolor (B) E. marginata - 112 -

그림 47. 세포의종류와구조가유사하여현미경상의식별이어려운수종의 현미경사진 ( 방사단면 ) (A) C. obtusa (B) T. plicata 각수종별유전정보를 NCBI로부터수집한결과 Eucalyptus 속과 Fokienia hodgisii 수종은전체 chloroplast 유전자가등록되어있다. 그외의수종은엽록체유전자가 partial DNA로만등록되어있는것으로확인되었다. 특히 P. resinosa는 NCBI 에등록된 DNA의수가 13,497개이지만, 그에비하여등록된엽록체 DNA 의개수는 19개에불과하며이또한 partial DNA일뿐이었다. Acer 속과 Chamaecyparis 속, 그리고 Thuja 속또한 partial DNA 로등록되어있었다. 그러나 Pinus 속들은 partial DNA라고해도그크기가약 120,000bp였으나, 이에비하여활엽수종들은평균약 2,000bp로매우작은크기만이 partial DNA 로등록되어있었다. 또한, Acer saccharinum는 1,095 개의 DNA 중 1,067개의엽록체 DNA 가등록되었으나모두 partial DNA로등록되어있었다.( 표 15) 이러한사실로미루어보아현재 NCBI에등록된 DNA 데이터베이스의정보가매우부족한실정이며, 추후분자생물학적기술을이용한수종식별을위해서는다양한수종에대한유전자데이터베이스의구축이필요하다고사료된다. - 113 -

표 15. 수종별유전정보수집 연번수종일반명 NCBI DNA 등록개수 등록최대 DNA 완성도 크기 (bp) NCBI 등록번호 1 Pinus densiflora 소나무 1,210(103) Partial 119,799 JN854210 2 Pinus resinosa Red pine 13,497(19) Partial 120,179 FJ899556 3 Pinus palustris Long leaf pine 121(34) Partial 120,397 JN854176 4 Pinus rigida pitch pine 165 (54) Partial 120,844 JN854163 5 Eucalyptus diversicolor Karri 52(9) Complete 160,214 KC180795 6 Eucalyptus marginata Jarrah 194(30) Complete 160,076 NC_022390 7 Acer saccharum Sugarmaple 809(96) Partial 1,440 L13181 8 Acer saccharinum Silver maple 1,095(1067) Partial 1,404 DQ659815 9 Chamaecyparis obtusa 편백 563(58) Partial 2,439 AB232479 10 Thuja plicata Red cedar 260(36) Partial 3,575 AY664823 11 Fokienia hodginsii Fukien cypress 327(30) Complete 128,789 KX832623-114 -

1) Pinus densiflora 와 Pinus resinosa 의비교 가 ) Primer 디자인 앞서진행한나. 소나무 (Pinus densiflora) 를이용한수종식별을위한후보유전자확보검증을통해 P. densiflora의 NCBI 데이터를기준으로총 15개의분석유전자부위를 primer로제작하였고, 이중유전자증폭하여염기서열이가능한 primer는총 10개로선별되었다 (primer 4, 6, 7, 8, 9, 10, 12, 14, 15 번 ). 이 10개의분석유전자부위를 NCBI에등록된또다른 Pinus 속인 P. resinosa의유전자서열과비교한결과, primer 9, 10, 12, 15번이상 4개의유전자부위가존재하는것을확인하였다. 이중 primer 15번에의한분석유전자부위가 (rpl20-rps18) 다른부위에비하여염기서열의길이가 560bp로비교적짧기때문에쉽게증폭이가능할것이라예상되어선택하였다. 본연구에서는 P. resinosa로부터 DNA 추출가능성을확인하고자하였다. 또한, 현재 P. densiflora와 P. resinosa는현미경적분석으로구별이어려운수종으로알려져있다. NCBI에등록된데이터상으로는 P. densiflora와 P. resinosa의이유전자부위염기서열간차이가있는것으로확인되었다.( 그림49) 이에 P. densiflora와 P. resionosa에대하여 15번 rpl20-rps18 분석유전자부위를각각증폭, sequencing 결과를비교분석하여, 서로다른유전자부위를확인함으로써, 분자생물학적기술을이용한두수종간의식별가능성을확인하고자하였다. 나 ) DNA 추출및유전자증폭 염기서열을분석하기위하여선행과제 ( 용역과제계약번호 : 20170205951-00) 에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. densiflora와 P. resinosa로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 rpl20-rps18f와 rpl20-rps18r를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의용량으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, - 115 -

denaturation 단계는 94 에서 30초, annealing 단계는 53 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 50회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 2% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 48. 에서약 500bp에서증폭하고자하는유전자를확인할수있었고, 얻은유전자의염기서열을분석하였다. 그림 48. 목재시편에서얻은 DNA 를 template 로하여 rrpl20-rps18f 와 rpl20-rps18r primer 로증폭하여얻은유전자. (a) P. densiflora; (b) P. resinosa - 116 -

그림 49. NCBI 에등록되어있는 P. densiflora 와 P. resinosa 의 rpl20 - rps18 유전자의염기서열비교및 primer 디자인 - 117 -

다 ) rpl20 - rps18 염기서열 증폭된유전자를 rrpl20-rps18f 와 rpl20-rps18r primer 를이용하여양 방향으로염기서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. P. densiflora CCAAATCGTAAAACACTTGCACAAATAGCTGTATTAGATAGTAATTGCTTTTCC ACCATCTTTAACAACTTATCGTATGATGAAATAGGTTAGGTATAGATGAAATAA ATAGGTAAAGATGGAATGGATAGAACAGATTCTCCCGGAGAATTCCCTCCGGG AAGGTCGAAATATTGAAAAGTTTTTTCATATTGGATTGGAAATGAACGAAAAG AATTCAATCCAATTCTTTTCCAAAAATGAAATCCTGTTGACTTTTTCAACAATA GACTCAAGATCTGCATCTTTATCGAACATATATTCCAGCTCCGATTTGATTAAG GAGTAGGCTTATTTCCCATGGATCAAATTAGTTTTCATTATAAAGAAAAGGTAA CAAAGATAGAATACGAGCTCGTTTAATAGCGTTAGTCATTAGACGTTGTTGTTT TGAAGTTAATCTATTCACTCGGCCGGATAAA P. resinosa CCAAATCGTAAAACACTTGCACAAATAGCTGTATTAGATAGTAATTGCTTTTCC ACCATCTTTAACAACTTATCGTATGATGAAATAGGTTAGGTATAGATGAAATAA ATAGGTAAAGATGGAATGGATAGAACAGATTCTCCCGGAGAATTCCCTCCGGG AAGGTCGAAATATTGAAAAGTTTTTTCATATTGGATTGGAAATGAACGAAAAG AATTCAATCCAATTCTTTTCCAAAAATGAAATCCTGTTGACTTTTTCAACAATA GACTCAAGATCTGCATCTTTATCGAACATATATTCCAGCTCCGATTTGATTAAG GAGTAGGCTTATTTCCCATGGATCAAATTAGTTTTCATTATAAAGAAAAGGTAA CAAAGATAGAATACGAGCTCGTTTAATAGCGTTAGTCATTAGACGTTGTTGTTT TGAAGTTAATCTATTCACTCGGCCGGATAAA - 118 -

라 ) rpl20 - rps18 염기서열분석 얻어진 rpl20 - rps18의염기서열을 NCBI의 nucleotide blast 검색을하였다. 표 16와같이 P.densiflora는 query coverage 98%, identities 100% 로 P. taeda 검색되었고, P. resinosa는 query coverage 99%, identities 100% 로 P. tabuliformis가검색되었다. 표 16. P. densiflora와 P. resinosa 시료에대한 rpl20-rps18f와 rpl20-rps18r primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 수종 NCBI검색결과 Query cover Identity Pinus taeda, complete genome 98% 100% Pinus cubensis isolate CUBE01 98% 100% Pinus densiflora Pinus echinata isolate ECHI01 98% 100% Pinus elliottii isolate ELLI01 98% 100% Pinus greggii isolate GREG02 98% 100% Pinus tabuliformis 99% 100% Pinus taivanensis 99% 100% Pinus resinosa Pinus yunnanensis var. tenuifolia 99% 100% pinus fragilissima isolate FRAG01 99% 100% Pinus massoniana chloroplast 99% 100% NCBI에등록되어있는 P. densiflora와 P. resinosa의 rpl20 - rps18 non-coding 염기서열과앞서증폭된 PCR product의염기서열을각각 align하였다. 본실험에서얻은 P. densiflora의유전자염기서열은 NCBI에등록되어있는유전자와 4개의부위에서차이가있는것을확인할수있었고 ( 그림 50), 실험에서얻은 P. resinosa의유전자염기서열은 NCBI에등록되어있는유전자와 1개의부위에서차이가있는것을확인할수있었다.( 그림 51) 수종간차이를확인하기위하여본실험에서얻은 P. densiflora와 P. resinosa 시료에대한 rpl20 - rps18 유전자의염기서열을비교한결과두개의부위에서차이가존재하는것을확인할수있었다.( 그림 52) 이와같은차이는시료가 NCBI에등록되어있는소나무와분류학상종간의차이인지개체간의차이인지는추가적인연구를통해밝혀야될것이며, 본연구를통 - 119 -

해사용된 P. densiflora와 P. resinosa의목재시료로부터 DNA가안정적으로추출됨을확인하였다. 각수종에대한 rpl20 - rps18 유전자의데이터베이스가구축되고, 이를활용하여판단한다면본실험에서제시하는두부분의차이가 P. densiflora와 P. resinosa의수종식별에이용될수있는지판단될수있을것이다. - 120 -

그림 50. NCBI 에등록된 P. densiflora 의 rpl20 - rps18 non-coding 염기서열과본실험에서 P. densiflora 로부터얻은유전 자염기서열의 align 한결과 - 121 -

그림 51. NCBI 에등록된 P. resinosa 의 rpl20 - rps18 non-coding 염기서열과본실험에서 P. resinosa 로부터얻은유전자 염기서열의 align 한결과 - 122 -

그림 52. 본실험에서 P. densiflora 와 P. resinosa 로부터얻은유전자염기서열의 align 한결과 - 123 -

2) Pinus palustris 와 Pinus rigida 의비교 가 ) Primer 디자인 앞서진행한나. 소나무 (Pinus densiflora) 를이용한수종식별을위한후보유전자확보검증을통해 P. densiflora의 NCBI 데이터를기준으로총 15개의분석유전자부위를 primer 로제작하였고, 이중유전자증폭하여염기서열이가능한 primer 는총 10개로선별되었다 (primer 4, 6, 7, 8, 9, 10, 12, 14, 15 번 ). 이 10개의분석유전자부위를 NCBI에등록된또다른 Pinus 속인 P. palustris, P. rigida의유전자서열과비교한결과, 두수종에모두 primer 서열이존재하며, 그길이가짧은 8번유전자 trnd GUC trny GUA 가증폭이쉽게이뤄질것이라고예상되어선택하였다. NCBI에등록되어있는 P. palustris와 P. rigida의 trnd GUC - trny GUA 유전자의염기서열비교한결과그림 54와같이두개의부위에서두수종간에는염기서열의차이가없는것으로확인되었다. 비록두수종을구별하기에 trnd GUC - trny GUA 유전자는적당하지못하나수종식별을위한지표유전자가제시되지못한상황에서본연구의목적인유전자의안정적인추출의가능여부는확인할수있으므로 trnd GUC - trny GUA 유전자의분석을진행하였다. 사용한 primer 정보는표 2에나타내었다. 나 ) DNA 추출및유전자증폭 trnd GUC - trny GUA 유전자 primer를이용하여유전자를증폭하고염기서열을분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 P. palustris와 P. rigida로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnd-guc와 trny-gua를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 55 에서 30 초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 50회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과 - 124 -

는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 53에서약 350bp에서증폭하고자하는유전자를확인할수있었고, 염기서열을분석하였다. 그림 53. 목재시편에서얻은 DNA 를 template 로하여 trnd-guc 와 trny-gua primer 로증폭하여얻은유전자 (a) P. palustris; (b) P. rigida - 125 -

그림 54. NCBI 에등록되어있는 P. palustris 와 P. rigida 의 trnd GUC - trny GUA non-coding 유전자의염기서열비교및 primer 디자인 - 126 -

다 ) trnd GUC - trny GUA 염기서열 증폭된유전자를 trnd-guc 와 trny-gua primer 를이용하여양방향 유전자분석을의뢰하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. P. palustris TTTTCAGGGCGGTACTCTAACCAATTGAGCTACAATCCCAATACAGTACAGTTC ACCTACTATTGGATAATATTTATTCCATGATAGGTGCTAGATAGGTCATATAGA TTATGCGAGTGCTAGGTCGATTAAATATCTTAATCTTCTCTTTCATTTTTGAAA TGTATCGATTCATACGGAATCGGGCATCTACGATATGAATAGATATCGATGCC GGGGTTCAATAACCAATTATCTTTTCATTCATGATTAGCATGAATATAACCATA CCGATAGATTGGTATTGATGATATTGGTTGGGTCGAGCTGGATTTGAACCAGC GTAGGCATATTGCCAACGGATTTACAGTCCGTCCTCATTAACCACTCGGGCAT P. rigida TTTCAGGGCGGTACTCTAACCAATTGAGCTACAATCCCAATACAGTACAGTTCA CCTACTATTGGATAATATTTATTCCATGATAGGTGCTAGATAGGTCATATAGAT TATGCGAGTGCTAGGTCGATTAAATATCTTAATCTTCTCTTTCATTTTTGAAAT GTATCGATTCATACGGAATCGGGCATCTACGATATGAATAGATATCGATGCCG GGGTTCAATAACCAATTATCTTTTCATTCATGATTAGCATGAATATAACCATAC CGATAGACTGGTATTGATGATATTGGTTGGGTCGAGCTGGATTTGAACCAGCG TAGGCATATTGCCAACGGATTTACAGTCCGTCCTCATTAACCACTCGGGCAT - 127 -

라 ) trnd GUC - trny GUA 염기서열분석 얻어진 trnd GUC - trny GUA 의염기서열을 NCBI의 nucleotide blast 검색을하였다. 표 17와같이 P. palustris는 query coverage 99%, identities 100% 로 P. palustris isolate PALU02로검색되었고, P. rigida는 query coverage 98%, identities 99% 로 P.rigida가검색되었다. 표 17. P. palustris와 P. rigida시료에대한 trnd-guc, trny- GUA primer로 증폭한유전자의염기서열을 NCBI data base에검색한결과 수종 NCBI검색결과 Query cover Identity Pinus palustris isolate PALU02 98% 100% Pinus patula isolate PATU01 98% 100% P. palustrispinus pungens isolate PUNG01 98% 100% Pinus radiata isolate RADI01 98% 100% Pinus rigida isolate RIGI01 98% 100% Pinus rigida isolate RIGI01 98% 99% Pinus serotina isolate SERO01 98% 99% P. rigida Pinus caribaea isolate CARI01 98% 99% Pinus oocarpa voucher DSG711 plastid, 98% 99% Pinus clausa isolate CLAU02 98% 99% NCBI에등록되어있는 P. palustris와 P. rigida의 trnd GUC - trny GUA non-coding 염기서열과본실험에서얻은유전자의연기서열을각각해당하는유전자와 align하였다. P. palustris에서얻은유전자의염기서열은 NCBI에등록되어있는유전자염기서열과모두일치함을확인할수있었고 ( 그림 55), P. rigida에서얻은유전자의염기서열은하나의염기서열에서 NCBI에등록되어있는유전자와불일치되는것을확인할수있었다.( 그림 56) 분석된수종간의차이를확인하기위하여 P. palustris와 P. rigida 시료에서얻은 trnd GUC - trny GUA 유전자의염기서열을비교한결과 1개의염기서열에서차이가있음을확인할수있었다.( 그림 57) 이와같은결과는본실험에사용한 P. palustris와 P. rigida 시료에대하여 DNA의안정적추출로수종식별은가능하였으나, 하나의개체에대하여이뤄진연구이기때문에개체의차이, 종간의차이, 지역적차이인지에대한추가적인연구가필요하다. - 128 -

그림 55. NCBI에등록된 P. palustris의 trnd GUC - trny GUA non-coding 염기서열과본실험에서 P. palustris로부터얻은유 전자염기서열의 align 한결과 - 129 -

그림 56. NCBI에등록된 P. rigida의 trnd GUC - trny GUA non-coding 염기서열과본실험에서 P. rigida로부터얻은유전자 염기서열의 align 한결과 - 130 -

그림 57. 본실험에서 P. palustris 와 P. rigida 로부터얻은유전자염기서열의 align 한결과 - 131 -

3) Chamaecyparis obtusa 와 Thuja plicata 의비교 가 ) Primer 디자인 Chamaecyparis obtusa와 Thuja plicata는속수준의분류학적차이에도불구하고해부학적분석방법으로식별이어려운것으로알려져있는수종이다. 본연구에서는 Chamaecyparis 속의계통발생학적관계에관련한논문으로참고문헌으로하였고이에제시된 primer는본연구에서선정한 C. obtusa와 T. plicata와같은속인 T. occidentalis를증폭할수있었기에이와같은 primer를이용하였다.[33] 그유전자부위는 petg - trnp로 T. plicata에대하여이유전자부위는 NCBI에등록되지않았기때문에같은속인 Thuja occidentalis의 NCBI data를활용하여다음그림 58과같이비교하였고, 두수종간의차이는그림에서처럼 18개의부위에서차이가있는것을확인할수있었다.( 그림 58) 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열을분석하기위하여선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 C. obtusa와 T. plicata로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 petg-f와 trnp-r를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1 μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 3분, denaturation 단계는 94 에서 1분, annealing 단계는 58 에서 1분, extension 단계는 72 에서 1분 30초로하여, denaturation 단계부터 extension 단계를 45회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 59에서약 500bp에서증폭하고자하는유전자를확인할수있었고, 염기서열을분석하였다. - 132 -

그림 58. 목재시편에서얻은 DNA 를 template 로하여 petg 와 trnp primer 로증폭하여얻은유전자 (a) C. obtusa; (b) T. plicata - 133 -

그림 59. NCBI 에등록되어있는 C. obtusa 와 T. occidentalis 의 petg - trnp 유전자의염기서열비교및 primer 디자인 - 134 -

다 ) petg trnp 염기서열 증폭된유전자를 petg 와 trnp primer 를이용하여양방향으로염기서 열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. C. obtusa GGCTGGATTATTTGTAACTGCATATTTACAATACAGGCGCGGCGATCAGTTGG ACCTTTGATTAATTAACATCTCTTTTTTTGATTGACCTCCTCCTTTCTTGAATAT CCAGGAGGTCAAATTCAGATTAATGTTCAAGTTAGTGTTTCAGCCAAATTCGAT TGAAGAAGATCCGAATCACGCTCTGTAGGATTTGAACCTACGACATCGGGTTT TGGAGACCCACGTTCTACCGAACTGAACTAAGAGCGCTTTCTTATCATAAAAGA TAAGACTGTAAAGAAAAGGATTGGCCCCAATACATCTTGTATGCATATACTATA TAGTATCATAAGAATGAAAGATTCTATATGATATGTCCAATGTGAATTGATCTC AAAAAATACCTCGTTACTGCTCCTTTGAGCAGTAATAGGTAGGGATGACAGGA TTTGAACCCGTGACATTTTGTACCCAAAACAAACGCGCTACC T. plicata TGGATTATTTGTAACTGCATATTTACAATACAGGCGCGGCGATCAGTTGGACC TTTGATTAATTAACATCTCTTTTTTTGATTGACCTCCTCCTTTCTTGAATATCCA GGAGGTCAAATTCAGATTAATGTTCAAGTTAGTGTTTCAGCCAAATTCGATTGA AGAAGATCCGAATCACGCTCTGTAGGATTTGAACCTACGACATCGGGTTTTGG AGACCCACGTTCTACCGAACTGAACTAAGAGCGCTTTCTTATCATAAAAGATAA GACTGTAAAGAAAAGGATTGGCCCCAATACATCTTGTATGCATATACTATATAG TATCATAAGAATGAAAGATTCTATATGATATGTCCAATGTGAATTGATCTCAAA AAATACCTCGTTACTGCTCCTTTGAGCAGTAATAGGTAGGGATGACAGGATTT GAACCCGTGACATTTTGTACCCAAAACAAACGCGCTA 라 ) petg trnp 염기서열분석 얻어진 petg trnp 유전자의염기서열을 NCBI 의 nucleotide blast 검 색을하였다. 표 18 과같이두수종에서모두 query coverage 100%, identities 100% 로 Carthamus tinctorius chloroplast ( 잇꽃 ) 로검색되었다. - 135 -

표 18. C. obtusa와 T. plicata 시료에대한 petg와 trnp primer로증폭한 유전자의염기서열을 NCBI data base에검색한결과 수종 NCBI검색결과 Query cover Identity Carthamus tinctorius chloroplast 100% 100% Carthamus tinctorius chloroplast 100% 100% C. obtusa Centaurea diffusa chloroplast 100% 99% Cynara syriaca chloroplast 100% 98% Cynara cardunculus var. sylvestris 100% 98% Carthamus tinctorius chloroplast 100% 100% Carthamus tinctorius chloroplast 100% 100% T. plicata Centaurea diffusa chloroplast 100% 99% Cynara syriaca chloroplast 100% 98% Cynara cardunculus var. sylvestris 100% 98% NCBI에등록되어있는 C. obtusa와 T. occidentalis의 petg - trnp non-coding 염기서열과앞서증폭된 PCR product의염기서열을각각 align하였다. C. obtusa의 PCR product 염기서열은 NCBI의 data base와차이가있는것을확인할수있었고 ( 그림 60), T. plicata의 PCR product도 NCBI의 data base와차이가있는것을확인할수있었다 ( 그림 61) 두수종간차이를확인하기위하여 C. obtusa와 T. plicata 시료에대한 petg-trnp PCR product sequence를비교한결과모든시퀀스가일치함을확인할수있었다. ( 그림 62) 두수종으로부터 DNA의추출, 염기서열분석을위한유전자증폭재현성을확인하였다. 그러나 NCBI 데이터검색결과 C. obtusa와 T. plicata가아닌 Carthamus tinctorius ( 잇꽃 ) 이검색되었으며이는각수종에대한데이터베이스구축의필요성을나타낸다. 또한두시료샘플의 align결과차이가없지만이는한개체에대한결과이기때문에많은개체의 petg-trnp 유전자증폭및염기서열비교연구를통한데이터베이스구축후수종식별유전자에대한가능성을판단할필요가있다. - 136 -

그림 60. NCBI 에등록된 C. obtusa 의 petg-trnp non-coding 염기서열과본실험에서 C. obtusa 로부터얻은유전자염기서 열의 align 한 - 137 -

그림 61. NCBI 에등록된 T. occidentalis 의 petg-trnp non-coding 염기서열과본실험에서 T. plicata 로부터얻은유전자 염기서열의 align 한결과 - 138 -

그림 62. 본실험에서 C. obtusa 와 T. plicata 로부터얻은유전자염기서열의 align 한결과시료에대한 - 139 -

4) Eucalyptus diversicolor 와 Eucalyptus marginata 의비교 가 ) Primer 디자인앞서진행한나. 소나무 (Pinus densiflora) 를이용한수종식별을위한후보유전자확보검증단계에서이용한 15가지 non-coding 유전자부위를 NCBI에등록되어있는 E. diversicolor와 E. marginata의데이터베이스를기준으로하였고, 두수종사이의차이가존재하여식별이가능하다고판단되는부위인 trnw psaj 의유전자간부위의 primer를디자인하였다.( 그림 63) 나 ) DNA 추출및유전자증폭 디자인한 primer를이용하여유전자를증폭후염기서열을분석하기위하여 DNA를추출하였고, 선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 E. diversicolor와 E. marginata으로부터 DNA를추출하였다. 이를 template로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnw, psaj 를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 E. diversicolor 55 E. marginata 54 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 45회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 64에서확인할수있듯이약 270bp에서증폭하고자하는유전자크기의밴드를확인할수있었고, sequencing을통한염기서열분석이가능하다고판단되어 코스모진텍에염기서열분석을의뢰하였다. - 140 -

그림 63. NCBI 에등록되어있는 E. diversicolor 와 E. marginata 의 trnw psaj 유전자간의염기서열비교및 primer 디 자인 - 141 -

그림 64. 목재시편에서얻은 DNA 를 template 로하여 trnw, psaj primer 로 증폭하여얻은유전자 (a) E. diversicolor (b) E. marginata - 142 -

다 ) trnw psaj 염기서열 앞의실험에서증폭된 PCR product 를 trnw, psaj primer 를이용하여 양방향유전자분석을의뢰하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였다. 이를통해얻어진유전자의염기서열은다음과같다. E. diversicolor TTGGAGACCCACGTTCTACCGAACTGAACTAAGAGCGCTTTATTATCACAATAA ATATTTTGTGACCTAACTCATCTTGGATATATATACTATCATAAATAAGAAAAT TAAAGATTGATTATGTATATCCAATGTTAATCAATCTCAATTGACCCCTCGTTA CTGTTCATAAGATAAGTAATAGGTAGGGATGACAGGATTTGAACCCGTGACAT TTTGTACCCAAAACAAACGCGCTACCGAGCTGCGCTACATCCCTTTCAATTAGC CTACAGTGTCATTGTAGAGAATCCTTGTCTTGTTTTCCACATCTTTATTTCTTT CAGTGATATACCAAATTATCTTGTCATTTCTTAATTTTTTTTAGTCTCATATCAT TATAACATATAATAATAGACTTATATTATATACGTACTAAAGAAATATGAAACC AAAAAAATGAATATTTTTTGGGAATTCTCAGACAAAAAGGGACGTATTTTTTTT TGTTTTACAAAAAGGAATATCTACTGAATCGTTGTACATTTCAATTTGCATTAG GAATTCTGCGTATAAATCCAAGTGTCAGTCATTAACTACATATACACATCTGGT CATATATGTATTACCATTATGTATTACAAATTGTAATAAAATTACAATAACAAA TAACAAAGAAGGAGGATTTTAAATGCGAAATCTAAAAACATACCTTTCCGTGG CACCGGTAGTAAGTACTCTATGGTTTGGGGCTTTAGCAA E. marginata ACCCCTCGTTACTGTTCATAAGATAAGTAATAGGTAGGGATGACAGGATTTGA ACCCGTGACATTTTGTACCCAAAACAAACGCGCTACCGAGCTGCGCTACATCCC TTTCAATTAGCCTACAGTGTCATTGTAGAGAATCCTTGTCTTGTTTTCCACATC TTTATTTCTTTCAGTGATATACCAAATTATCTTGTCATTTCTTAATTTTTTTTAG TCTCATATCATTATAACATATAATAATAGACTTATATTATATACGTACTAAAGA AATATGAAACCAAAAAAATGAATATTTTTTGGGAATTCTCAGACAAAAAGGGA CGTATTTTTTTTTGTTTTACGAAAAGGAATATCTACTGAATCGTTGtACATTTCA ATTTGCATTAGGAATTCTGCGTATAAATCCAAGTGTCAGtCATTAACTACATAT ACACATCTGGtCATATATGTATTACCATTATGTATTACAAATTGTAATAAAATTA CAATACCAATTAACAAAGAAGGAGGA - 143 -

라 ) trnw psaj 염기서열분석 얻어진 trnw psaj 유전자의염기서열을 NCBI blast 검색결과 E. diversicolor는 Morus australis plastid ( 뽕나무 ) 가 query cover 99%, identity 99% 로검색되었으며, E. marginata 또한 Morus australis plastid ( 뽕나무 ) query cover 100%, identity 99% 로검색되었다.( 표 19) 표 19. E. diversicolor와 E. marginata 시료에대한 trnw psaj primer로 증폭한유전자의염기서열을 NCBI data base에검색한결과 수종 NCBI검색결과 Query cover Identity Morus australis plastid 99% 99% Morus alba var. multicaulis plastid 99% 99% E. diversicolormorus alba var. multicaulis chloroplst 99% 99% Morus mongolica chloroplast 99% 99% Morus cathayana plastid 99% 99% Morus australis plastid 100% 99% Morus alba var. multicaulis plastid 100% 99% E. marginata Morus alba var. multicaulis chloroplst 100% 99% Morus mongolica chloroplast 100% 99% Morus cathayana plastid 100% 99% - 144 -

NCBI에등록되어있는 E. diversicolor와 E. marginata의 trnw psaj non-coding 염기서열과앞서증폭된 PCR product의염기서열을각각 align하였다. E. diversicolor의 PCR product 염기서열은 NCBI의 data base와차이가있는것을확인할수있었고 ( 그림 65), E. marginata의 PCR product도 NCBI 의 data base와차이가있는것을확인할수있었다 ( 그림 66) 두수종간차이를확인하기위하여 E. diversicolor와 E. marginata 시료에대한 trnw psaj PCR product sequence를비교한결과대부분의염기서열이불치함을확인할수있었다. ( 그림 67) 두수종으로부터 DNA의추출, 유전자증폭, 염기서열분석은이루어졌으나 NCBI 데이터검색결과 E. diversicolor와 E. marginata가아닌 Morus australis가각각검색되었다. 이러한결과는목재로부터 DNA의추출은안정적으로이루어지며두수종간의구분은가능하나, 분석된유전자에대한데이터베이스의부족으로인하여수종의식별이불가능한것으로판단된다. - 145 -

그림 65. NCBI 에등록된 E. diversicolor 의 trnw psaj non-coding 염기서열과본실험에서 E. diversicolor 로부터얻은유 전자염기서열의 align 한결과 - 146 -

그림 66. NCBI 에등록된 E. marginata 의 trnw psaj non-coding 염기서열과본실험에서 E. marginata 로부터얻은유전 자염기서열의 align 한결과 - 147 -

그림 67. 본실험에서 E. diversicolor 와 E. marginata 로부터얻은유전자염기서열의 align 한결과 - 148 -

5) Acer saccharum 와 Acer saccharinum 의비교 가 ) Primer 디자인 두수종은모두전체시퀀스가 NCBI에등록되지않았기때문에, NCBI에등록되어있는 A. saccharum과 A. saccharinum 엽록체데이터중두수종모두존재하는유전자간 intergenic 부위 (non-coding) 를비교하였다. 그결과두수종간차이가존재하며, 비교적등록된데이터의수가많은유전자간부위인 trnh - psba 유전자간부위의 primer를디자인하였다.( 그림 68) 나 ) DNA 추출및유전자증폭 디자인한 primer 를이용하여유전자를증폭후염기서열을분석하기위하여 DNA를추출하였고, 선행연구에서도출된 목재 DNA 추출기본매뉴얼 을토대로 A. saccharum과 A. saccharinum으로부터 DNA를추출하였다. 이를 template 로하여유전자증폭을실시하였다. 멸균된 PCR tube에 3차증류수 13.9μl, 10x Taq reaction buffer 2μl, 10mM dntp mix 0.4μl, primer 10pmole/ μl의 trnh와 psba를각각 0.8μl씩첨가해주었다. 마지막에 Taq polymerase (5U/ μl BioFACT ) 0.1μl, DNA template 2μl순으로첨가하여총 20μl의볼륨으로 Thermal cycler (GenePro Thermal Cycler TC-E-48D) 를이용하여 PCR을진행하였다. Pre-denaturation 단계는 94 에서 5분, denaturation 단계는 94 에서 30초, annealing 단계는 A. saccharum은 52, Acer saccharium은 51 에서 30초, extension 단계는 72 에서 30초로하여, denaturation 단계부터 extension 단계를 45회반복하였다. 그후 72 에서 10분을 1회추가하였다. 실험결과는 1.5% agarose gel에 marker와 PCR product 모두 2μl씩 loading하여전기영동하여확인하였다. 그림 69에서확인할수있듯이약 480bp에서증폭하고자하는유전자크기의밴드를확인할수있었고, sequencing을통한염기서열분석이가능하다고판단되어 코스모진텍에염기서열분석을의뢰하였다. - 149 -

그림 68. NCBI 에등록되어있는 A. saccharum 과 A. saccharinum 의 trnh - psba 유전자의염기서열비교및 primer 디 자인 - 150 -

그림 69. 각수종으로부터얻은 DNA 를 template 로하여 trnh, psba primer 로 증폭하여얻은유전자 (a) A. saccharum (b) A.saccharinum - 151 -

다 ) trnh psba 염기서열 증폭된유전자를 trnh-f 와 psba-r primer 를이용하여양방향으로염 기서열을분석하였고, 그결과를 vector NTI 를통하여 align 후 assemble 하였 다. 이를통해얻어진유전자의염기서열은다음과같다. A. saccharum TGCCTTGATCCACTTGGCTACCAGTTACCAAGTATGAAGAGTGGGAGGATGAC TTTGACACCTCAGTTCCACCAGAGAAAGGGTTCAACTTCAAGAGTCCACTTCCA AAGTTACTCCAAGGGAAGAAGGAGCCACTGGTGCAACCTAAGCCAAAGTATGA CCCACTTACTGAGACTGACCCTGACATTATTCGCTCAATGCCAGACAGTGACTT CAAGCTTGCTATTGGTAACATTAAGGACACAGCCCTGAAGCTTGCCATGATTAA GGAACGAGCTGAACAAATGCTTGTGAGAGATAATGGCTTGAAAGAGGTTCGTG TTGCTCCAATGTTTATGTGGTATGGGTTCAACTCAAAGGAACTTACTGGCTCAT TTGTAG A. saccharinum TGCCTTGATTTCTCTTGGCTACCAGTTACCAAGTATGAAGAGTGGGAGGATGA CTTTGACACCTCAGTTCCACCAGAGAAAGGGTTCAACTTCAAGAGTCCACTTCC AAAGTTACTCCAAGGGAAGAAGGAGCCACTGGTGCAACCTAAGCCAAAGTATG ACCCACTTACTGAGACTGACCCTGACATTATTCGCTCAATGCCAGACAGTGACT TCAAGCTTGCTATTGGTAACATTAAGGACACAGCCCTGAAGCTTGCCATGATTA AGGAACGAGCTGAACAAATGCTTGTGAGAGATAATGGCTTGAAAGAGGTTCGT GTTGCTCCAATGTTTATGTGGTATGGGTTCAACTCAGAGGAACTTACTGGCTC ATTTGAGATTTATGACGAGCTTGTGACCTTCTACATGACGCAAGAGACCAGAG CAGACGTTGAGGAGCTTGTGAGACAACTTCAAGGAGCTTCAAAGC - 152 -

라 ) trnh psba 염기서열분석 얻어진 trnh psba 유전자의염기서열을 NCBI blast 검색결과 A. saccharum 과 A. saccharinum 두수종모두특정한수종이검색되지않았다. 표 20. A. saccharum과 A. saccharinum 시료에대한 trnh, psba primer로증폭한유전자의염기서열을 NCBI data base에검색한결과 Query 수종 NCBI검색결과 Identity cover A. saccharum 검색결과없음 A. saccharinum 검색결과없음 NCBI에등록되어있는 A. saccharum과 A. saccharinum의 trnh-psba 의 non-coding 염기서열과앞서증폭된 PCR product의염기서열을각각 align 하였다. 두가지수종모두 NCBI에등록되어있는데이터와대부분일치하지않는것으로확인되었다.( 그림 70, 71) 수종간의차이를확인하기위하여 A. saccharum과 A. saccharinum 시료에대한 trnh-psba PCR product sequence를비교한결과 4개의서열에서차이가있음을확인할수있었다.( 그림 72) 두수종으로부터 DNA의추출, 유전자증폭, 염기서열분석은이루어졌으나 NCBI 데이터검색결과어떠한수종도검색되지않았다. 이는 Acer 속수종에대한데이터베이스구축의필요성을나타낸다. 또한이는 NCBI에등록되지않은새로운개체의 Acer 속수종일가능성이있다고판단된다. 두시료샘플의 align결과, 4개의염기서열에서차이가있는것으로보아증폭시킨 trnh-psba 유전자간부위는 A. saccharum과 A. saccharinum 간수종의식별이가능한후보유전자가될가능성이있다고판단된다. 그러나이는하나의개체에대한결과이기때문에많은개체의 trnh-psba 유전자간부위증폭및염기서열비교연구를통한데이터베이스구축후수종식별유전자에대한가능성을판단할필요가있다. - 153 -

그림 70. NCBI 에등록된 A. saccharum 의 trnh - psba non-coding 염기서열과 A. saccharum 로부터얻은유전자염기서열 의 align 한결과 - 154 -

그림 71. NCBI 에등록된 A. saccharinum 의 trnh - psba non-coding 염기서열과 A. saccharinum 로부터얻은유전자염기 서열의 align 한결과 - 155 -

그림 72. A. saccharum 과 A. saccharinum 로부터얻은유전자염기서열의 align 한결과 - 156 -

6) 목재의유전자식별을위한다양한수종에서의 DNA 안정적추출검증및 염기서열분석결론 본연구는국내유통되는여러수종중, 해부학적분석만으로는정확한식별이불가능한대표적인수종을선정하여목재의유전자식별을위해다양한수종으로부터 DNA가안정적으로추출되고, 그결과주어진비교시료간의식별이가능한지확인하였다. 현미경상으로종수준의식별이어려운국내주요유통수종으로 P. densiflora와 P. resinosa, P. palustris와 P. rigida, E. diversicolor와 E. marginata, A. saccharum과 A. saccharium를이용하였으며, 속수준의식별이어려운수종으로는 C. obtusa와 T. plicata를이용하였다. 이를대상으로유전자증폭및염기서열분석, NCBI 데이터검색을진행하였으며, 그결과 DNA 분석방법으로 P. densiflora와 P. resinosa, P. palustris와 P. rigida는염기서열의비교를통하여두수종간의구분은가능하였으나, NCBI 데이터검색결과같은 Pinus 속수준까지검색이되어식별에는한계가존재하였다. NCBI 데이터와분석시료의염기서열을비교해보면그차이는 10개미만으로존재하였다. 이러한결과는지역적, 지리적, 종간, 개체간등의차이에따라존재할수도있다. 따라서 DNA분석을통한목재의수종식별이이뤄지기위해서는각수종의한유전자간 non-coding부위에대하여많은개체에대한분석이필요하다고사료된다. 따라서본연구를통하여데이터베이스가구축이된후이를이용하기에 DNA가안정적으로추출이되었다고판단된다. 다음으로 A. saccharum과 A. saccharium으로부터 DNA 추출후유전자증폭, 염기서열분석, NCBI 데이터검색결과, 어떠한수종도검색되지않았다. 이러한결과는데이터베이스에분석에필요한충분한자료가결여된것으로판단된다. 두시료간의염기서열비교결과차이가존재하는것으로보아 A. saccharum과 A. saccharium 두수종에대한 trnh-psba non-coding 유전자는많은개체에대하여검증이이뤄져야한다고판단된다. 따라서두수종으로부터 DNA의안정적추출이이뤄졌고, 추후많은개체에대한식별가능유무의유전자검증이이뤄진다면종간의식별이이뤄질수있다고판단된다. 속수준에서해부학적수종식별이어려운 C. obtusa와 T. plicata에대하여 DNA 추출후유전자증폭하고염기서열을확보하여 NCBI data base 검색결과, 두수종의증폭된유전자의염기서열이일치하며, 모두 Carthamus tinctorius 잇꽃이검색되었다. NCBI 데이터베이스와두시료의염 - 157 -

기서열을비교한결과대부분의염기서열에차이가있는것을확인할수있었다. 이는 NCBI에등록되어있는해당수종의엽록체데이터개수는 C. obtusa 58개, T. plicata 36개로그수는 P. densiflora 103개에비하여적은수준이다. 또한등록되어있는데이터의종류도다양하지않다. 따라서다양한부위의유전자에대한개체간의데이터베이스가구축된다면종식별이가능할것으로사료된다. E. diversicolor와 E. marginata 또한 DNA추출, 유전자증폭염기서열을각각 NCBI 데이터와비교한결과완전히다른수종이검색되었다. 게다가유전자증폭부위끼리비교한결과도염기서열간에많은차이가존재했다. 이는두수종으로부터 DNA의추출, 유전자증폭, 염기서열분석은이루어졌으나, 이유전자부위는분석하고자하는종의식별은물론서로다른두수종의구분도불가능하다는것을의미하므로새로운유전자부위의탐색이필요하다고사료된다. 이러한결과를종합해봤을때, 목재로부터 DNA의추출은안정적으로이루어지며, 염기서열분석을위한유전자증폭또한가능하다는것을확인하였다. 그러나증폭된유전자의염기서열분석을통한종식별및수종구분을하기위해서는앞선결론에서지속적으로강조한것처럼, 보다더많은개체의 NCBI 데이터베이스구축, 수종식별이가능한분석유전자부위탐색및검증에대한연구가선행되어야한다고제안하는바이다. - 158 -

4. 과제수행달성도 가. 문제유통수종의 DNA 분석유전자 ( 분석마커 ) 도출및개발 1) 문제유통수종선별및수종에대한 DNA 정보수집 Eucalyptus 속과 Fokienia hodgisii 수종은전체 chloroplast 유전자가등록. 그외의수종은엽록체유전자가 partial DNA로등록됨. P. resinosa 는 NCBI 에등록된 DNA의수가 13,497개이지만, 그에비하여등록된엽록체 DNA의개수는 partial DNA로 19개에불과. Acer 속과 Chamaecyparis 속, 그리고 Thuja 속또한 partial DNA 로등록. Pinus 속들은 partial DNA라고해도그크기가약 120,000bp, 이에비하여활엽수종들은평균약 2,000bp로매우작은크기만이 partial DNA로등록. Acer saccharinum 는 1,095개의 DNA 중 1,067개의엽록체 DNA가등록되었으나모두 partial DNA 등록. (113p~114p, 표 15) 2) NCBI 정보및문헌조사를통한수종식별을위한분석유전자발굴 32편의문헌검색을통하여 trnl intron, trnl-trnf intergenic spacer, trnk/matk intron 영역들은특정계통에대해서분류학적으로만족스러운정보를나타내며, 그외의그룹에대해서는식물분류학적으로낮은분해능이나타난다는점을확인. 또한특정한몇개의유전자부위에대한연구는그외의다양한유전자부위에대한데이터가부족함확인. (28-29p) 42편의문헌검색과 NCBI 데이터를바탕으로 P. densiflora와 P. resinosa 시료에대한 rpl20 - rps18 non-coding 유전자, P. palustris 와 P. rigida 시료에대한 trnd GUC - trny GUA non-coding 유전자,C. obtusa 와 T. occidentalis 시료에대한 petg trnp non-coding 유전자, E. diversicolor 와 E. marginata 시료에대한 trnw psaj non-coding 유전자, A. saccharum 과 A. saccharinum 시료에대한 trnh-psba non-coding 유전자발굴. - 159 -

나. 도출된유전자 ( 분석마커 ) 를이용한목재 DNA 분석방법검증 1) DNA 분석검증을통한재현성및실용성확인 가. 수종간차이를확인하기위하여본실험에서얻은 P. densiflora와 P. resinosa 시료에대한 rpl20 - rps18 유전자의염기서열을비교한결과 2개의부위에서차이가존재하는것을확인.( 그림52) 연구에사용된 P. densiflora와 P. resinosa 의목재시료로부터 DNA의안정적추출확인과동시에 DNA분석방법을이용하여 P. densiflora 와 P. resinosa 시료의구분가능확인. (119p ~ 120p) 나. 수종간의차이를확인하기위하여 P. palustris 와 P. rigida 시료에서얻은 trnd GUC - trny GUA 유전자의염기서열을비교한결과 1개의염기서열에서차이가있음을확인.( 그림 57) 본연구에사용된 P. palustris와 P. rigida 의목재시료로부터 DNA의안정적추출확인과동시에 DNA분석방법을이용하여 P. palustris 와 P. rigida 시료의구분가능확인.(128p) 다. 수종간의차이를확인하기위하여 C. obtusa 와 T. occidentalis 시료에서얻은 petg-trnp 유전자의염기서열을확인할수있음. 각수종의목재시료의안정적 DNA추출, 염기서열분석을위한유전자증폭재현성확인.( 그림 60, 61) 두수종의식별가능한지표유전자제시에따른향후수종식별기술확보. (136p) 라. 두수종간차이를확인하기위하여 E. diversicolor 와 E. marginata 시료에대한 trnw psaj PCR product sequence를비교한결과대부분의염기서열이불치함확인. ( 그림 67) 두목재시료로수종으로부터 DNA의추출, 유전자증폭, 염기서열분석, 시료구분가능.(145p) 마. 수종간의차이를확인하기위하여 A. saccharum 과 A. saccharinum 시료에대한 trnh-psba PCR 유전자의염기서열을비교한결과 4개의서열에서차이가있음을확인.( 그림 72) 두목재시료로수종으로부터 DNA의추출, 유전자증폭, 염기서열분석, 시료구분가능.(153p) - 160 -

다. 결과산출물 Ÿ 목재종식별을위한주요유통수종의 DNA 분석방안마련및검증 최종 보고서 Ÿ 특허출원 (1 건 ) - 사포를이용한목분의 DNA 분리방법, 대한민국특허출원번호 : 10-2017-0045388, 2017 년 04 월 07 일 Ÿ 포스터발표 (1건) - 분자생물학적분석을활용한수종감정에서 DNA 추출의효율성증대를위한목재의전처리, 이자민, 오정애, 손석규, 김태종, 2017 한국목재공학회춘계학술대회, 전북대학교, 전주, 대한민국 (2017년 4월 14 15일 ) Ÿ 홍보 (1 건 ) 2017 년 05 월 31 일 순번언론사제목 1 내일신문 낡은목재도 DNA로수종분석 2 이코노믹톡뉴스 국내최초, 목재에서 DNA추출기술개발해특허출원까지 3 동아일보 국민대임산생명공학과, 국내최초로 ' 목재 DNA' 추출기술개발 - 161 -

5. 정책제안 가. 산학연협약을통한유전자분석을통한목재수종식별서비 스구축 목재의식별은현미경관찰에의한해부학적차이에의해이루어지고있다. 그러나일부수종의경우, 종 (Species, 種 ) 간의식별이나특히종 ( 種 ) 이하의분류학적차이는해부학적관찰로구별하기가어려운경우가종종발생한다. 이를해결하기위한보조수단으로유전자의염기서열을분석하여비교하는방법이제안되고있는데많은제약이있다. 본연구에서는목재로부터안정적으로 DNA를추출하는방법을제안하여하나의큰기술적문제를해결하고있다. 또한유전자를선별하여비교시료간의종의식별이가능할수있음을제안하고있다. 그러나유전자분석을통한목재수종식별서비스를제공하기위해서는여러가지해결해야할문제점이있다. 첫번째문제점으로목재수종식별을위한지표유전자의확보가필요하다. 이를위해서는관련되는많은목재의유전자분석자료를확보하고비교평가하여지표유전자를도출하여야하는데많은시간과노력이필요하다. 국립산림과학원에서는목재의수종식별을위한지표유전자확보를위해연구를수행하고있는데이미소나무 (Pinus densiflora) 와구주적송 (Pinus sylvestris) 을구별하기위한지표유전자를선별하여특허를제출하였다. 또한 2018년부터 3년간의연구개발로국내에서유통되는목재의식별을위한지표유전자를발굴하고자계획하고있다. 두번째문제점인목재로부터안정적인유전자의추출은본연구과제를통 해해결하였다. 세번째문제점으로지표유전자의선택적인증폭이필요하다. 유전자를증폭은중합효소연쇄반응 (PCR) 을통해이루어지는데이때사용하는 DNA polymerase가목재의추출물성분에의해활성저해를받는다. 이를극복하기위해서는효소활성저해성분에의해활성이영향을받지않는효소의개발이필요하다. 최근에국내에서유전자분석의주요기업인 SolGent에서는효소의활성저해물질에영향을받지않고안정적으로유 - 162 -

전자를증폭하는 DNA polymerase를개발하여시판을앞두고있다. 이를활용하면안정적으로원하는유전자를목재의추출물 DNA를이용하여선택적으로증폭할수있으며본보고서에는나타내지않았으나본연구진의연구를통해서도우수한선택적유전자증폭을확인하였다. 이에유전자분석을통한목재수종식별을서비스를위한대부분의문제를 해결했다고판단되며다음과같은 유전자분석을통한목재수종식별을 서비스시스템 을제안한다. 다음의각기술을양해각서를통해확보하고 Ÿ 지표유전자의발굴 국립산림과학원산림유전자원부 Ÿ 목재로부터 DNA의추출 국민대학교임산생명공학과 Ÿ 유전자의증폭기술 SolGent 주식회사유전자분석을통한목재수종식별을서비스를한국임업진흥원에서수행한다. 현재국립산림과학원산림유전자원부에서지표유전자를제안한수종에대해서만목재의수종식별서비스를해부학적수종식별의보완서비스로이용한다. 향후국립산림과학원산림유전자원부의연구에따라지표유전자가제공되는대로추가적인서비스수종을확대한다. 유전자분석을통한목재수종식별을보완서비스가필요한경우한국임업진흥원에서국민대학교에의뢰하여 DNA를추출하고이를 SolGent에의뢰하여유전자를증폭하고염기서열을분석하여수종의식별한결과를한국임업진흥원에통보한다. ( 또는 SolGent에서 DNA추출기술을이전받아한국임업진흥원으로부터직접서비스를의뢰받아유전자의염기서열분석을통한수종을식별하고이를한국임업진흥원에통보한다.) 비용은추가유전자분석이필요한수종의경우서비스의뢰인에게추가비용을계산한다. 이와같은방법으로현재소나무 (Pinus densiflora) 와구주적송 (Pinus sylvestris) 을구별하기위한서비스를곧바로시작할수있으며, 추후국립산림과학원산림유전자원부의지표유전자제안에따라서비스를확대해갈수있다. - 163 -

나. 유전자분석을통한지표유전자발굴방안 유전자분석을통한목재수종식별을서비스를수행하기위해서는다른모든문제가해소된현재의상황에서지표유전자의도출이매우중요한다. 국립산림과학원산림유전자원부에서는 3개년의연구과제수행을통해목재수종식별을위한지표유전자를제시하고자한다. 그러나많은유전자가등록되어있고이들을통합적으로분석한다면보다정확하고목적에맞는지표유전자를제안할수있을것이다. 이를위해서는많은목재수종의유전자등록이필요하나이를위한비용은막대하다. 이러한유전자를수집하고확보하기위해다음과같은제도를제안한다. 국제기구를통해목재의생산량이부피나금액단위로일정량이상되면제재업자는벌목현장에서의무적으로나뭇잎등의시료를정해진수만큼채취하여공신력이인정된기관에제공하고유전자분석을의뢰한다. 분석된엽록체전체유전자는모든사람들이접근할수있는 Data base에등록한다. 추후에목재를거래하거나이용하여제품을생산하면유전자분석을통해등록되어있는자료와비교할수있게하는방법이다. 또한국제간의거래시에이와같은유전자등록이되어있지않은목재는수출입을엄격히금한다. 이때유전자분석과 data base등록은목재생산업자가부담하며소량생산업자는경제적부담이클수있으므로일정규모이상의목재생산자에게의무를지우며일정규모이하의목재생산자에게는현재의인허가방식을이용한다. 이와같은새로운제도를다음과같은효과를얻을수있다. 1. 유전자 data base가부족하여도유통수종의유전자가등록되어있으므로직접적으로비교가능하며수종식별과원산지확인을위한보조자료로할용가능 2. 목재수종의유전자 data base가지속적으로확충되며수종식별의목적에정확한지표유전자의제시가쉬워짐. - 164 -

참고문헌 1. Chloroplast genome sequences from total DNA for plant identification. 2011. 9:328-333 2. The Use of Chloroplast DNA to Resolve Plant Phylogenies: Noncoding versus rbcl Sequences. Molecular Biology and Evolution. 1994. 1:769-777 3. The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. American Journal of Botany. 2005. 92:142-166 4. The evolution of non coding chloroplast dna and its application in plant systematics. Annals of the Missouri Botanical Garden. 2000. 87:482-498 5. Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Molecular Biology. 1991. 17:1105-1109 6. Hairpins create minute inversions in non-coding regions of chloroplast DNA. Curr Genet. 1996. 30:259 262 7. Phylogenetic utility of the internal transcribed spacers of nuclear ribosomal DNA in plants: an example from the Compositae. Molecular phylogenetics and evolution. 1992. 1:3-16 8. The ITS region of nuclear ribosomal DNA: a valuable source of evidence on angiosperm phylogeny. Annals of the Missouri Botanical Garden. 1995. 82: 247 277. 9. The phylogeny of the Alstroemeriaceae, based on morphology, rps16 intron, and rbcl sequence data. Systematic Botany. 2003. 28:47 69-165 -

10. Phylogenetics of seed plants: an analysis of nucleotide sequences from the plastid gene rbcl. Annals of the Missouri Botanical Garden. 1993. 80: 528-580 11. Phylogenetic relationships in disa based on non-coding trnl-trnf chloroplast sequences: evidence of numerous repeat regions. American Journal of Botany. 2001. 88:2088 2100. 12. Power and limitation of the chloroplast trnl(uaa) intron for plant DNA barcoding. Nucleic acids research. 2007. 35:e14 13. Wood identification of Japanese cyclobalanopsis species (Fagaceae) based on DNA polymorphism of the intergenic spacer between trnt and trnl 5' exon. Journal of Wood Science. 2001. 47:81-86 14. Novel perspectives in wood certification and forensics: dry wood as a source of DNA. Proceeding of the Royal Society of London B. 2002. 269:1039-1046 15. A set of primers for the amplification of chloroplast microsatellites in Quercus. Molecular Ecology Notes. 2003. 3:24-27 16. Checking the geographical origin of oak wood: molecular and statistical tools. Molecular Ecology Notes. 2003. 12:1629-1636. 17. Use of chloroplast microsatellites to differentiate oak populations. Annals of Forest Science. 2004. 61:825-830 18. Extraction, amplification and characterization of wood DNA from Dipterocarpaceae. Plant Molecular Biology. 2006. 24:45-55 19. Genetic analysis of archaeological wood remains: first results and prospects. Journal of Archaeological Science.2006. 33:1216-1227 - 166 -

20. Ancient olive DNA in pits: Preservation, amplification and sequence analysis. Journal of Archaeological Science. 2006. 33:77-88 21. Chloroplast DNA analysis in oak stands (Quercus robur L.) in North Rhine-Westphalia with presumably Slavonian orirgin: Is there an association between geographic origin and bud phenology? Journal of Applied Botany and Food Quality. 2007. 81:1165-171. 22. Results of molecular analysis of an archaeological hemp (Cannabis sativa L.) DNA sample from North West China. Genetic Resources and Crop Evolution. 2008. 55:481-485. 23. DNA from processed and unprocessed wood: Factors influencing the isolation success. Forensic Science International Genetics. 2009. 3:185-192. 24. Wood identification with PCR targeting noncoding chloroplast DNA. Plant Molecular Biology. 2011. 77:609-617. 25. Molecular database for classifying Shorea species (Dipterocarpaceae) and techniques for checking the legitimacy of timber and wood products. Journal of Plant Research. 2011. 124:35-48. 26. A real-time PCR method for the detection of trnl-trnf sequence in agarwood and products from Aquilaria (Thymelaeaceae). Conservation Genetics Resources. 2012. 4:803-806. 27. Identification of plant remains in underwater archaeological areas by morphological analysis and DNA barcoding. Advances in Anthropology. 2013. 3:240-248. 28. Identification of market adulterants in East Indian sandalwood using DNA barcoding. Annals of Forest Science. 2014. 72:517-522. - 167 -

29. Quantitative evaluation of properties of residual DNA in Cryptomeria japonica wood. Journal of Wood Science. 2015. 61:1-9. 30. Review: Ancient plant DNA in archaeobotany. Vegetation History and Archaeobotany. 2008. 17:233-244. 31. Development and Characterisation of Microsatellite Loci in Eucalyptus globulus (Myrtaceae). Silvae Genetica. 2001. 50:89-91. 32. A chloroplast DNA hypervariable region in eucalypts. Theoretical and Applied Genetics. 2000 101:473 477. 33. Historical biogeography and phylogenetic relationships of the genus Chamaecyparis (Cupressaceae) inferred from chloroplast DNA polymorphism. Plant Systematics and Evolution. 2003. 241:13-28. 34. Genetic diversity and differentiation of core vs. peripheral populations of eastern white cedar, Thuja occidentalis(cupressaceae). American Journal of Botany. 2012. 99:690 699. 35. Phylogenetic Analysis of Pines Based on Chloroplast trnt-trnl Intergenic Spacer DNA Sequences. Journal of Forest and Environmental Science. 2014. 30: 307-313. 36. Efficient Identification of the Forest Tree Species in Aceraceae Using DNA Barcodes. Frontiers in plant science. 1016. 7: 1707. 37. Variation in TaqI-digested DNA of Sugar and Black Maples Is Independent of Taxon and Plant Origin HortScience. 2001. 36:1327-1328. 38. Variation in the psbc gene region of gymnosperms and angiosperms as detected by a single restriction site polymorphism. Theoretical and Applied Genetics. 1997 94:1065-1071. - 168 -

39. Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well. PLoS One. 2008. 30:2802. 40. Noncoding chloroplast DNA variation in Mexican pines. Theoretical and Applied Genetics. 1995. 91:1101-1106. 41. Phylogenetic relationships of Eurasian pines (Pinus, Pinaceae) based on chloroplast rbcl, MATK, RPL20-RPS18spacer, and TRNV intron sequences. American Journal of Botany. 1999. 12:1742 1753. 42. Phylogenetics of Acer (Aceroideae, Sapindaceae) based on nucleotide sequences of two chloroplast non-coding regions. Harvard Papers in Botany. 2006. 11:101-115. 43. Chloroplast DNA phylogeography of the endangered Japanese red maple (Acer pycnanthum): the spatial configuration of wetlands shapes genetic diversity. Diversity and Distributions. 2009. 15:917 927 44. An organellar maturase associates with multiple group II introns. Proceedings of the National Academy of Sciences. 107:3245 3250-169 -

부록 목재시료사진 (scale bar 1cm) Pinus densiflora Pinus resinosa Pinus palustris Pinus rigida - 170 -

Chamaecyparis obtusa Thuja plicata Eucalyptus diversicolor Eucalyptus marginata - 171 -

Acer saccharum Acer saccharinum - 172 -