J Plant Biotechnol (2024) 51:143-151
Published online May 29, 2024
https://doi.org/10.5010/JPB.2024.51.015.143
© The Korean Society of Plant Biotechnology
조주련・박태호
대구대학교 원예학과
Correspondence to : e-mail: thzoo@daegu.ac.kr
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Solanum jamesii is a wild Solanum species with a native range from northwestern Mexico to southwestern United States. It is a resource for potato breeding due to its resistance to late blight. However, S. jamesii is diploid and has an endosperm balanced number (EBN) of one, which hampers the sexual hybridization between S. jamesii and cultivated potatoes. To overcome this problem, somatic fusion can be performed. After fusion, molecular markers are necessary to select fusion products efficiently. In this study, therefore, the chloroplast genome of S. jamesii was sequenced using next-generation sequencing (NGS) technology and was compared with those of eight other Solanaceae species to develop specific markers for S. jamesii. The total length of the chloroplast genome sequence of S. jamesii is 155,576 bp. The structural organization of the chloroplast genome is almost identical to those of other Solanum species. Phylogenetic analysis with 17 other Solanaceae species showed that S. jamesii is most closely related to S. cardiophyllum. Detailed analysis of the chloroplast genome sequence of nine Solanum species identified two S. jamesii-specific InDels and 11 S. jamesii-specific SNPs. Based on these results, six S. jamesii-specific PCR-based markers were developed. The results of this study can be used to identify S. jamesii from other Solanum species and will contribute to the development of new potato varieties using S. jamesii.
Keywords cpDNA, PCR-based marker, Potato, SNPs, Solanum jamesii
재배되고 있는 감자(Solanum tuberosum L.)의 야생종 중 하나인 S. jamesii는 멕시코의 북서부부터 미국의 남서부 지역까지 자생하며, 재배종 감자와 같이 괴경을 형성하는 반면 배수성이 2배체로 재배종 감자의 4배체와는 차이가 있다(Hanneman 1994; Hawkes 1990). S. jamesii는 풋마름병, 더뎅이병, 감자역병, 콜로라도감자잎벌레, 뿌리혹선충, 감자바이러스X, Y (PVX, PVY)와 함께 건조한 환경에도 강한 것으로 알려져 있어 신품종 감자를 육성하는데 중요한 재료로 활용될 수 있다(Bamberg et al. 2021; Kim-Lee et al. 2005; Kozlov et al. 2007; Machida-Hirano 2015; Srivastava et al. 2016; Watanabe et al. 2011). 저항성과 같은 유용 형질을 재배종 감자에 도입하여 감자 신품종을 육성하기 위해서는 전통적인 교배육종 방법을 이용할 수 있으나 S. jamesii의 EBN (Endosperm Balanced Number)은 1로 EBN이 4인 재배종 감자와의 차이가 있어 교배육종이 불가능하다(Cho et al. 1997; Hawkes 1990; Ortiz and Ehlenfeldt 1992; Spooner et al. 2014). 이러한 생리적 교잡 불화합성을 극복하기 위해서는 형질전환에 의한 유전자
도입(genetic modification, GM), 순차적 교배(bridge crosses), 체세포잡종(somatic hybridization) 등이 제시되고 있다(Cho and Park 2014). S. jamesii를 이용한 감자 신품종 육성을 위해 순차적 교배 방법이 시도된 바 있으나(Bamberg et al. 2021), S. jamesii와 S. tuberosum을 체세포융합 방법을 이용하여 체세포잡종 계통을 육성하고자 하는 시도는 아직 이루어진 바가 없어 연구가 진행중에 있으며, 향후 체세포잡종 계통을 확보하였을 때, 체세포잡종 계통에서의 엽록체 DNA의 구성을 확인하기 위한 목적으로 S. jamesii 전체 엽록체 유전체의 염기서열 정보를 구명하고 이를 이용하여 S. jamesii 특이적인 분자마커를 개발하였다.
속씨식물의 엽록체 유전체는 일반적으로 4분할의 원형 이중가닥 분자구조로 한 쌍의 inverted repeats (IRs) 영역과 각각 하나의 small single copy (SSC) 및 large single copy (LSC) 영역으로 구성되어 있다(Yurina and Odintosova 1998). 앞서 보고된 바와 같이 재배종 감자와 근연야생종의 엽록체 전체 염기서열과 그 구조는 서로 매우 유사하였다(Palmer 1991; Raubeson and Jansen 2005; Saski et al. 2005; Sugiura et al. 1998)(Table 1). 이러한 유사성에도 불구하고, 이들 엽록체 유전체에서는 유전자 역위나 재배열 등과 같은 구조적 변화에 의한 변이가 발생하여 엽록체 유전체 전체를 비교하면 식물체 종 특이적인 InDel이나 SNP를 발견할 수 있다(Calsa Junior et al. 2004; Cho et al. 2015; Jheng et al. 2012; Kim et al. 2005, 2018; Kim and Park 2019, 2020a, 2020b; Park 2021, 2022a, 2022b, 2023a, 2023b; Saski et al. 2005). 따라서, 본 연구에서는 S. jamesii의 엽록체 전장유전체를 구명하여 이에 대한 정보를 제공하고, 전장유전체의 염기서열 정보를 이용하여 개발한 PCR기반의 S. jamesii 특이적 분자마커에 대한 결과를 제공하고자 한다.
Table 1 . Results of comparative analysis of the chloroplast genome sequence of S. jamesii with those of 14 Solanaceae species
Species | Accession no. | Total length (bp) | GC content (%) | Total no. of genes | No. of tRNAs | No. of rRNAs | Reference |
---|---|---|---|---|---|---|---|
S. jamesii | MK690624 | 155,576 | 37.90 | 134 | 36 | 4 | This study |
S. hjertingii | MK690623 | 155,545 | 37.88 | 134 | 36 | 4 | Park (2023a) |
S. cardiophyllum | MK690622 | 155,570 | 37.90 | 134 | 36 | 4 | Park (2023b) |
S. acaule | MK036506 | 155,570 | 37.84 | 135 | 36 | 4 | Park (2022b) |
S. brevicaule | MK036507 | 155,531 | 37.87 | 135 | 36 | 4 | Park (2022a) |
S. demissum | MK036508 | 155,558 | 37.87 | 135 | 36 | 4 | Park (2021) |
S. hougasii | MF471372 | 155,549 | 37.87 | 135 | 36 | 4 | Kim and Park (2020b) |
S. stoloniferum | MF471373 | 155,567 | 37.87 | 135 | 36 | 4 | Kim and Park (2020a) |
S. chacoense | MF471371 | 155,532 | 37.89 | 136 | 36 | 4 | Kim and Park (2019) |
S. berthaultii | KY419708 | 155,533 | 37.88 | 137 | 39 | 4 | Kim et al. (2018) |
S. commersonii | KM489054 | 155,525 | 37.88 | 133 | 33 | 4 | Cho et al. (2016) |
S. nigrum | KM489055 | 155,432 | 37.90 | 139 | 39 | 4 | Cho and Park (2016) |
S. tuberosum | KM489056 | 155,312 | 37.88 | 130 | 30 | 4 | Cho et al. (2016) |
S. bulbocastanum | DQ347958 | 155,371 | 37.88 | 133 | 30 | 4 | Daniell et al. (2006) |
S. tuberosum | NC008096 | 155,296 | 37.88 | 131 | 36 | 4 | Gargano et al. (2005) |
*The data have been partially adopted from Park (2023b).
S. jamesii (PI578326) 엽록체의 전장유전체 분석에는 S. jamesii 계통 중 하나인 SJ-8이 이용되었으며, S. jamesii 특이적 분자마커 개발을 위해 S. jamesii의 SJ-6과 SJ-11 계통이 추가되었다. 이외에 3개의 재배종 감자 품종과 12종의 가지속(Solanum) 이용되었으며, 자체적으로 대구에서 수집한 S. nigrum (SN)을 제외한 모든 식물재료는 국립식량과학원 고령지농업연구소에서 분양 받아 기내식물체로 유지 및 증식하여 이용하였다. 실험에 이용된 재배종 감자의 품종은 ‘대지’(DJ), 대서’(DS), ‘수미’(SM)이고, S. nigrum 이외 11개의 야생종 계통은 S. acaule (PI310970, SA), S. brevicaule (PI205394, SB2), S. hougasii (PI161174, SH2), S. cardiophyllum (PI341233, SC1), S. commersonii (PI558050, SC2), S. hjertingii (PI186559, SH1), S. stoloniferum (PI160224, SS), S. kurtzianum (PI578236, SK), S. mochiquense (PI338616, SM1), S. verrucosum (PI160228, SV1), S. vernei (PI230468, SV2)이다.
DNA 분리는 Genomic DNA Extraction kit (Plants)(RBC, New Taipei City, Taiwan)를 이용하여 모든 식물체를 대상으로 계대배양 한 후 약 3주 이내 유식물체의 약 100 mg을 채취하여 진행하였다.
S. jamesii의 엽록체 전장유전체 분석은 S. jamesii 계통 SJ-8을 대상으로 진행하였다. 전체 염기서열 정보(raw data)는 Macrogen (Macrogen, Seoul, South Korea)의 MiSeq (Illumina, SanDiego, CA, USA) 플랫폼을 통해 확보하였으며, 염기서열의 분석과 조립은 파이젠(http://phyzen.com)에서 제공하는 de novo assembly가 이용되었다(Cho et al. 2015). 그 과정에서 전처리(trimming), high-quality 염기서열 추출(Phred score 20 이상), dnaLCW 방법으로 조립 등은 CLC quality trim, CLC de novo assemble 등의 프로그램이 이용되었으며, 이들 프로그램은 CLC genome assembler package version 4.2.1 (CLC Inc, Rarhus, Denmark)에서 제공된 것이다. 이후 조립된 염기서열 중, 엽록체 유전체 유래의 contig를 선발하고, gap-filling에 의한 확장과 통합을 통해 엽록체 유전체를 완성하였다. 전처리된 염기서열 전체를 완성된 엽록체 전장 유전체 서열에 mapping하여 각각의 영역별 read depth를 계산하였고, reference 종과 유전체의 구조 및 염기서열을 비교하였다. 또한, 엽록체 전장 유전체 서열을 이용하여 BLASTN을 수행하였으며, homology가 높은 종의 엽록체 전장 유전체 서열을 확인하고 BLASTZ (Schwartz et al. 2003) 분석을 실시하였다.
완성된 S. jamesii 엽록체 전장 유전체의 유전자는 Geseq 프로그램을 통해 확인하고, BLAST search를 기반으로 한 manual curation을 통해 최종적으로 유전자와 유전자의 위치를 결정하였다(Tillich et al. 2017). 이들은 OGDraw 프로그램을 통해 원형의 유전자지도 상에 표시되었다(Lohse et al. 2007).
S. jamesii와 유사종들과의 유전관계 분석을 위해 NCBI (National Center for Biotechnology Information)에 등록되어 있는 재배종 감자 S. tuberosum (NC008096 및 KM489056), 가지과의 까마중 S. nigrum (KM489055), 고추 Capsicum annuum (JX270811), 토마토 S. lycopersicum (NC007898), 감자의 근연야생종 12종 S. acaule (MK036506), S. brevicaule (MK036507), S. berthaultii (KY419708), S. bulbocastanum (DQ347958), S. demissum (MK036508), S. cardiophyllum (MK690622), S. stoloniferum (MF471373), S. chacoense (MF471371), S. commersonii (KM489054), S. hjertingii (MK690623), S. hougasii (MF471372), S. verrucosum (MK690625) 등 총 18개의 엽록체 전체 염기서열을 확보하여, 여기에 존재하는 단백질 코딩 서열을 MEGA11 프로그램을 이용하여 분석하였으며, 이 프로그램은 1,000번의 Bootstrap 옵션이 적용되었고 GTR (General Time Reversible)모델 기반의 Maximum Likelihood 방법이 이용되었다(Tamura et al. 2021).
분자마커의 개발은 S. jamesii의 엽록체 전체 염기서열과 NCBI로부터 얻은 S. verrucosum (MK690625), S. cardiophyllum (MK690622), S. bulbocastanum (DQ347958), S. berthaultii (KY419708), S. hjertingii (MK690623), S. tuberosum (KM489056), S. commersonii (KM489054), S. hougasii (MF471372)의 엽록체 전체 염기서열을 ClustalW2 (EMBL: http://www.ebi.ac.uk/Tools/msa/clustalw2)를 이용하여 다중정렬 하였고, 이를 통해 S. jamesii에 특이적으로 나타나는 InDel과 SNP를 확보할 수 있었다. InDel의 경우 InDel이 있는 영역에 primer를 제작하여 allele 특이적 분자마커를 개발하고자 하였으며, SNP의 경우 tetra-primer ARMS (Amplification Refractory Mutation System)-PCR 기술을 이용하여 primer를 제작하였다(Collins and Ke 2012; Ye et al. 2001). S. jamesii와 S. tuberosum을 포함한 총 18개의 계통, S. jamesii (SJ6, SJ8 및 SJ11), S. tuberosum (감자품종 ‘대지’(DJ), ‘수미’(SM), ‘대서’(DS)), S. acaule (SA), S. cardiophyllum (SC1), S. brevicaule (SB2), S. commersonii (SC2), S. hjertingii (SH1), S. hougasii (SH2), S. mochiquense (SM1), S. kurtzianum (SK), S. nigrum (SN), S. stoloniferum (SS), S. verrucosum (SV1), S. vernei (SV2)를 대상으로 PCR을 수행하였다. InDel 기반의 분자마커 개발은 위한 PCR의 경우, 총 20 ul (약 20 ng genomic DNA, 0.5 mM dNTPs, 10 pMol each primer, 1 U Taq polymerase (Genetbio, Daejeon, South Korea))으로, SNP 기반의 분자마커 개발을 위한 PCR의 경우, 총 20 ul (약 20 ng genomic DNA, 0.5 mM dNTPs, 10 pMol each inner primer, 1 pMol each outer primer, 1 U Taq polymerase (Genetbio, Daejeon, South Korea))으로 Thermocycler (Biometra, Göttingen, Germany)을 이용하였다. PCR의 결과는 핵산 염색용액인 StaySafe (RBC, New Taipei City, Taiwan)를 이용하여 1.5% agarose gel에서 확인하였다.
S. jamesii의 엽록체 전체 염기서열은 NGS 기술을 이용하여 완성되었다. NGS 기반으로 얻은 데이터는 평균 301.0 bp 길이의 총 염기서열 약 1.5 Gbp로, 이는 5,267,110 read를 통해 생산되었다. 이중 약 90.79%를 이용하여 전처리 및 Phred score 20 이상의 데이터 추출과정을 거쳐 high-quality 총 염기서열 정보 약 1.3 Gbp를 얻었으며, 평균 길이는 283.6 bp로 확인되었다. 이는 전체 raw data의 약 85.55%였으며, S. cardiophyllum 의 전체 엽록체 염기서열과 비교했을 때 사용된 총 염기서열이 약 0.3 Gbp 적었지만, 총 생산된 데이터 대비 추출 염기서열의 비율은 비슷한 수준이었다(Park 2023a, 2023b). 또한, de novo assembly와 엽록체 유래의 contig를 선발하고, gap-filling 과정을 통한 확장과 통합 등을 거쳐 총 염기서열은 155,576 bp의 S. jamesii 엽록체 전장 유전체가 완성되다. 최종적으로 aligned read의 개수는 435,578개, 전체 유전체에 대한 평균 coverage는 약 789.54x, 전체 trimmed read를 완성된 전체 유전체 서열에 mapping한 결과로 나타난 read depth는 최소 378 이상이었다. 구조적으로는 원형 이중가닥 분자 구조로 다른 식물의 엽록체 유전체와 동일한 형태였으며, 한 쌍의 25,596 bp의 IRa/IRb, 18,394 bp의 SSC, 85,990 bp의 LSC 영역으로 구성되어 있었다. S. bulbocastanum의 엽록체 전체 염기서열(DQ347958)과의 비교와 BLASTZ 분석의 결과로도 그 유사성을 확인할 수 있었다(Daniell et al. 2006; Schwartz et al. 2003) (Fig. 1). S. jamesii의 엽록체 전체 염기서열을 대상으로 NCBI의 BLASTN 분석 결과 S. pinnatisectum (MZ233590), S. cardiophyllum (MK690622)과 99.95%의 높은 염기서열 일치도를 보였으며, GC 비율 또한 37.90%로 매우 유사하였다(Table 1).
S. jamesii의 엽록체 전장 유전체에는 105개의 단백질 코딩, 45개의 tRNA, 8개의 rRNA를 포함는 총 158개의 유전자가 존재하고 있는데, 이중 11개의 단백질 코딩 유전자, 9개의 tRNA, 4개의 rRNA가 한 쌍의 Ira 및 IRb 영역에 역배열로 중복되어 있었다(Table 1, Fig. 2). 전체 염기서열 중 코딩 영역의 비율은 59.2%로 나타났으며, 158개 유전자의 평균 길이는 약 583.4 bp, 이중 각각의 평균은 단백질 코딩 유전자가 764.8 bp, tRNA가 62.4 bp, rRNA가 1,131.0 bp의 크기로 각각 51.6%, 1.8%, 5.8%의 비율로 분포되어 있으며, 이러한 결과는 다른 Solanum 종들과 매우 유사한 수준이었다(Table 1).
S. jamesii의 계통수 작성은 총 17종 가지과의 엽록체 전장 유전체 코딩 서열을 기반으로 실행되었다(Fig. 3). 그 결과, S. jamesii는 S. tuberosum 뿐만 아니라 다수의 가지속(Solanum) 종들과 같은 그룹에 존재함을 확인하였으며, 유연관계가 그중 S. cardiophyllum과 가장 가까웠으며, 상대적으로 까마중(S. nigrum), 고추(C. annuum), 토마토(S. lycopersicum)와는 다소 먼 것으로 확인되었다.
S. jamesii 엽록체의 전체 염기서열을 7종의 근연야생종과 재배종 감자의 엽록체 전체 염기서열과 ClustalW2 (EMBL; https://www.ebi.ac.uk/Tools/msa/clustalw2)를 이용하여 다중 정렬한 결과 S. jamesii 특이적인 2개의 InDel과 11개의 SNP를 구명하였다(Fig. 4). 상대적으로 S. jamesii 특이적 InDel과 SNP의 수가 매우 적은 편이었는데, 이는 다중 정렬에 S. jamesii의 엽록체 전체 염기서열과 유사도가 매우 높은 S. cardiophyllum이 포함되어 나타난 결과로 보인다. 하지민, 다양한 유전적 요인에 의한 변이 발생으로 인해 S. jamesii 특이적인 InDel과 SNP가 존재하며, 이를 이용하여 S. jamesii 특이적 분자마커를 개발할 수 있었다(Cho and Park 2016; Chung et al. 2006; Kim et al. 2018; Kim and Park 2019, 2020a, 2020b; Park 2021, 2022a, 2022b, 2023a, 2023b). 2개의 InDel은 각각 유전자 간 및 내 영역에 존재하고 있었으며, 11개의 SNP 중 8개는 유전자 간 영역에 3개는 유전자 내 영역에 위치하고 있었다.
S. jamesii를 포함한 총 9개 Solanum 종의 엽록체 전체 염기서열의 다중 정렬 결과로 확인된 S. jamesii 특이적 InDel과 SNP는 S. jamesii 특이적인 분자마커 개발에 활용되었다. InDel은 InDel 영역에서 특이적인 primer를 제작하여 PCR을 통해 다양성의 차이가 나타나는 S. jamesii 특이적 또는 비특이적인 3개의 SCAR (Sequence Characterized Amplified Region) 마커를 개발할 수 있었다(Cho et al. 2015; Garcia-Lor et al. 2013; Yamaki et al. 2013)(Table 2, Fig. 5A, 5B, 5C). SNP는 11개 중 3개의 영역을 대상으로 S. jamesii 특이적 SNP 위치를 기준으로 SNP가 포함된 DNA 단편과 SNP 위치을 대상으로 S. jamesii 특이적이거나 혹은 S. jamesii를 제외한 Solanum 종들에 특이적으로 DNA를 증폭할 수 있는 tetra-primer ARMS-PCR 방법을 적용하여 primer를 제작하였다(Collins and Ke 2012; Ye et al. 2001)(Table 2). 이러한 방법은 PCR로 염기서열이 매우 유사한 SNP를 대상으로 primer의 제작이 어려운 종 또는 SNP가 존재하는 계통 간의 차이를 비교적 손쉽게 확인할 수 있다는 장점이 있어 널리 활용되고 있으며, 본 연구에서도 Solanum 종간의 비교 결과로 나타난 SNP를 활용한 PCR 기반의 S. jamesii 특이적이거나 (Fig. 5D, 5E, 5F) 혹은 비특이적으로(Fig. 5D, 5E, 5F) 나타나는 3개의 분자마커를 개발하였다(Han et al. 2016; Kou et al. 2017; Mahmoudi et al. 2019; Zou et al. 2018). 최종적으로 본 연구에서는 InDel과 SNP 정보를 기반으로 한 S. jamesii 특이적 분자마커가 각각 3개씩 개발되었다. SJ8_InDel1 (Fig. 4, Fig. 5A, 5B)과 SJ8_InDel2 (Fig. 4, Fig. 5C)은 각각 trnL 유전자 내 및 petA와 psbJ 유전자 사이에 존재하는 SCAR 마커로 개발되었으며, SJ8_SNP1 (Fig. 4, Fig. 5D), SJ8_SNP3 (Fig. 4, Fig. 5E)과 SJ8_SNP5 (Fig. 4, Fig. 5F)는 각각 atpF 유전자 내, atpH와 atpI 유전자 사이 및 psaA 유전자 내에 존재하는 SNP 기반의 마커로 개발되었다.
Table 2 . Information of primers used to generate S. jamesii-specific and non-specific markers
InDel/SNP | Strand | Primer sequence (5’ → 3’) | Tm (°C) | Size |
---|---|---|---|---|
SJ8_InDel1(a) | Forward primer | AAACGTATATACATACTATATC | 56°C | 735 bp (delete sequence 11bp specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel1(b) | Forward primer | AAACGTATATACATACGTATTG | 58°C | 746 bp (delete sequence 11bp non-specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel2 | Forward primer | AGTATAGAAAGGATTAGATTAG | 60°C | 777 bp (insert sequence ‘ATTAGATTAGTATAGAAAGG’ specific to S. jamesii) |
Reverse primer | TATCAGCAATGCAGTTCATC | |||
SJ8_SNP1 | Forward inner primer | CCTATCTTCTCTTCTTTGTTCATATTTCC | 60°C | 396 bp (from two outer primers), 217 bp (C allele specific to S. jamesii), 238 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTGGTCTTACATTAGTTTCGATATCTTCTT | |||
Forward outer primer | TCAAGATCCTCTGTTTTCGATTATCTAA | |||
Reverse outer primer | GTAAATTGTAAAAACAGGGATTTTACCC | |||
SJ8_SNP3 | Forward inner primer | AATTCTCTTGAGCCACCCACCATATGCC | 56°C | 300 bp (from two outer primers), 178 bp (C allele specific to S. jamesii), 210 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTTTAGCTTAGCCCACAGCAATGTATAGAT | |||
Forward outer primer | TGAAGGAAAGGGGAATTTTAGGAAAAAGAT | |||
Reverse outer primer | GTTCCTATACCTGTCATGCTCCTTGGATT | |||
SJ8_SNP5 | Forward inner primer | CTGAATAGAAACATTAAACTAAAAGCCTAG | 58°C | 390 bp (from two outer primers), 243 bp (T allele specific to S. jamesii), 204 bp (G allele non-specific to S. jamesii) |
Reverse inner primer | CTTTTTTTCCTAGGTGCTCATTTTATA | |||
Forward outer primer | TACTGCAATAATTCTTGCTAAGAAGAAT | |||
Reverse outer primer | GTAAGTGATCAAGGGGTAGTAACTCATA |
감자의 신품종 육성을 위해 다양한 유용한 특성을 지니고 있는 근연야생종을 이용한 연구가 계속해서 수행되고 있으나, 교배육종을 통한 유용 형질의 도입은 다수의 근연야생종이 재배종 감자와의 EBN과 배수성 차이로 인해 어려움을 겪기 때문에, 체세포융합을 이용한 체세포잡종의 육성이 활용되고 있으며, 이에 따라 유전자형 확인이 필수적이다. 분자마커를 통한 유전자형의 확인은 단순히 염색체 조합만을 확인하는 것이 아니라 엽록체 및 미토콘드리아 유전체의 전달 확인도 필요하다. 특히 Solanum 속에 속하는 다양한 식물종의 경우, 체세포융합이나 재분화 과정에서 핵 내 및 미토콘드리아 유전체는 재조합 빈도가 높다는 것이 알려져 있으며, 엽록체 유전체의 경우에는 한 쪽 친의 엽록체 유전체가 무작위적으로 전달되는 것 외에도 양친의 유전체가 모두 전달되는 경우도 있었다(Chen et al. 2013; Cho et al. 2016; Lössl et al. 2000; Mohapatra et al. 1998; Smyda-Dajmund et al. 2016; Wang et al. 2011; Xiang et al. 2004). 이에 따라, 본 연구에서 개발된 분자마커는 체세포융합을 통한 체세포잡종의 육성 과정에서 특히 S. jamesii와 S. tuberosum 유래의 엽록체 유전체의 유전자형 구별에 이용이 될 것이며, 이러한 분자마커는 감자 육성 과정에서 신품종의 개발 및 진화학적 연구에 기여할 수 있을 것이다.
멕시코 북서부부터 미국 남서부까지 자생하는 감자 야생종 중 하나인 Solanum jamesii는 Phytophthora infestans에 대한 저항성으로 감자 육종에 훌륭한 재료이다. 그러나 S. jamesii는 이배체이며 EBN 값이 1로, 이는 S. jamesii와 재배용 감자 간 교배가 불가능하다. 이 문제를 극복하기 위해 체세포융합을 이용할 수 있다. 이후에는 분자마커를 사용하여 효율적으로 융합체를 선발하는 것이 필수적이다. 본 연구에서는 차세대 유전체 기술을 사용하여 S. jamesii의 엽록체 전장 유전체 서열을 확보하고, 이를 다른 8개의 가지과에 속한 종과 비교하여 S. jamesii 특이적 마커를 개발하였다. S. jamesii의 엽록체 유전체의 총 서열 길이는 155,576 bp이며, 그 구조는 다른 Solanum 종과 유사하다. 17개의 다른 Solanaceae 종과의 계통수 분석 결과, S. jamesii는 S. cardiophyllum과 가장 밀접한 유연관계를 가졌다. 8개의 다른 Solanum 종의 엽록체 전장 유전체 서열 비교에서 S. jamesii에 특이한 2개의 InDel 및 11개의 SNP가 확인되었다. 이러한 InDel 및 SNP를 기반으로 S. jamesii를 다른 Solanum 종과 구별하는 6개의 PCR 기반 분자마커가 개발되었다. 본 연구의 결과는 Solanum 종 중에서 S. Jamesii를 구별하고 S. Jamesii를 활용한 감자 신품종 개발에 기여할 것이다.
J Plant Biotechnol 2024; 51(1): 143-151
Published online May 29, 2024 https://doi.org/10.5010/JPB.2024.51.015.143
Copyright © The Korean Society of Plant Biotechnology.
조주련・박태호
대구대학교 원예학과
Ju-Ryeon Jo ・Tae-Ho Park
Department of Horticulture, Daegu University, Gyeongsan 38453, South Korea
Correspondence to:e-mail: thzoo@daegu.ac.kr
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Solanum jamesii is a wild Solanum species with a native range from northwestern Mexico to southwestern United States. It is a resource for potato breeding due to its resistance to late blight. However, S. jamesii is diploid and has an endosperm balanced number (EBN) of one, which hampers the sexual hybridization between S. jamesii and cultivated potatoes. To overcome this problem, somatic fusion can be performed. After fusion, molecular markers are necessary to select fusion products efficiently. In this study, therefore, the chloroplast genome of S. jamesii was sequenced using next-generation sequencing (NGS) technology and was compared with those of eight other Solanaceae species to develop specific markers for S. jamesii. The total length of the chloroplast genome sequence of S. jamesii is 155,576 bp. The structural organization of the chloroplast genome is almost identical to those of other Solanum species. Phylogenetic analysis with 17 other Solanaceae species showed that S. jamesii is most closely related to S. cardiophyllum. Detailed analysis of the chloroplast genome sequence of nine Solanum species identified two S. jamesii-specific InDels and 11 S. jamesii-specific SNPs. Based on these results, six S. jamesii-specific PCR-based markers were developed. The results of this study can be used to identify S. jamesii from other Solanum species and will contribute to the development of new potato varieties using S. jamesii.
Keywords: cpDNA, PCR-based marker, Potato, SNPs, Solanum jamesii
재배되고 있는 감자(Solanum tuberosum L.)의 야생종 중 하나인 S. jamesii는 멕시코의 북서부부터 미국의 남서부 지역까지 자생하며, 재배종 감자와 같이 괴경을 형성하는 반면 배수성이 2배체로 재배종 감자의 4배체와는 차이가 있다(Hanneman 1994; Hawkes 1990). S. jamesii는 풋마름병, 더뎅이병, 감자역병, 콜로라도감자잎벌레, 뿌리혹선충, 감자바이러스X, Y (PVX, PVY)와 함께 건조한 환경에도 강한 것으로 알려져 있어 신품종 감자를 육성하는데 중요한 재료로 활용될 수 있다(Bamberg et al. 2021; Kim-Lee et al. 2005; Kozlov et al. 2007; Machida-Hirano 2015; Srivastava et al. 2016; Watanabe et al. 2011). 저항성과 같은 유용 형질을 재배종 감자에 도입하여 감자 신품종을 육성하기 위해서는 전통적인 교배육종 방법을 이용할 수 있으나 S. jamesii의 EBN (Endosperm Balanced Number)은 1로 EBN이 4인 재배종 감자와의 차이가 있어 교배육종이 불가능하다(Cho et al. 1997; Hawkes 1990; Ortiz and Ehlenfeldt 1992; Spooner et al. 2014). 이러한 생리적 교잡 불화합성을 극복하기 위해서는 형질전환에 의한 유전자
도입(genetic modification, GM), 순차적 교배(bridge crosses), 체세포잡종(somatic hybridization) 등이 제시되고 있다(Cho and Park 2014). S. jamesii를 이용한 감자 신품종 육성을 위해 순차적 교배 방법이 시도된 바 있으나(Bamberg et al. 2021), S. jamesii와 S. tuberosum을 체세포융합 방법을 이용하여 체세포잡종 계통을 육성하고자 하는 시도는 아직 이루어진 바가 없어 연구가 진행중에 있으며, 향후 체세포잡종 계통을 확보하였을 때, 체세포잡종 계통에서의 엽록체 DNA의 구성을 확인하기 위한 목적으로 S. jamesii 전체 엽록체 유전체의 염기서열 정보를 구명하고 이를 이용하여 S. jamesii 특이적인 분자마커를 개발하였다.
속씨식물의 엽록체 유전체는 일반적으로 4분할의 원형 이중가닥 분자구조로 한 쌍의 inverted repeats (IRs) 영역과 각각 하나의 small single copy (SSC) 및 large single copy (LSC) 영역으로 구성되어 있다(Yurina and Odintosova 1998). 앞서 보고된 바와 같이 재배종 감자와 근연야생종의 엽록체 전체 염기서열과 그 구조는 서로 매우 유사하였다(Palmer 1991; Raubeson and Jansen 2005; Saski et al. 2005; Sugiura et al. 1998)(Table 1). 이러한 유사성에도 불구하고, 이들 엽록체 유전체에서는 유전자 역위나 재배열 등과 같은 구조적 변화에 의한 변이가 발생하여 엽록체 유전체 전체를 비교하면 식물체 종 특이적인 InDel이나 SNP를 발견할 수 있다(Calsa Junior et al. 2004; Cho et al. 2015; Jheng et al. 2012; Kim et al. 2005, 2018; Kim and Park 2019, 2020a, 2020b; Park 2021, 2022a, 2022b, 2023a, 2023b; Saski et al. 2005). 따라서, 본 연구에서는 S. jamesii의 엽록체 전장유전체를 구명하여 이에 대한 정보를 제공하고, 전장유전체의 염기서열 정보를 이용하여 개발한 PCR기반의 S. jamesii 특이적 분자마커에 대한 결과를 제공하고자 한다.
Table 1 . Results of comparative analysis of the chloroplast genome sequence of S. jamesii with those of 14 Solanaceae species.
Species | Accession no. | Total length (bp) | GC content (%) | Total no. of genes | No. of tRNAs | No. of rRNAs | Reference |
---|---|---|---|---|---|---|---|
S. jamesii | MK690624 | 155,576 | 37.90 | 134 | 36 | 4 | This study |
S. hjertingii | MK690623 | 155,545 | 37.88 | 134 | 36 | 4 | Park (2023a) |
S. cardiophyllum | MK690622 | 155,570 | 37.90 | 134 | 36 | 4 | Park (2023b) |
S. acaule | MK036506 | 155,570 | 37.84 | 135 | 36 | 4 | Park (2022b) |
S. brevicaule | MK036507 | 155,531 | 37.87 | 135 | 36 | 4 | Park (2022a) |
S. demissum | MK036508 | 155,558 | 37.87 | 135 | 36 | 4 | Park (2021) |
S. hougasii | MF471372 | 155,549 | 37.87 | 135 | 36 | 4 | Kim and Park (2020b) |
S. stoloniferum | MF471373 | 155,567 | 37.87 | 135 | 36 | 4 | Kim and Park (2020a) |
S. chacoense | MF471371 | 155,532 | 37.89 | 136 | 36 | 4 | Kim and Park (2019) |
S. berthaultii | KY419708 | 155,533 | 37.88 | 137 | 39 | 4 | Kim et al. (2018) |
S. commersonii | KM489054 | 155,525 | 37.88 | 133 | 33 | 4 | Cho et al. (2016) |
S. nigrum | KM489055 | 155,432 | 37.90 | 139 | 39 | 4 | Cho and Park (2016) |
S. tuberosum | KM489056 | 155,312 | 37.88 | 130 | 30 | 4 | Cho et al. (2016) |
S. bulbocastanum | DQ347958 | 155,371 | 37.88 | 133 | 30 | 4 | Daniell et al. (2006) |
S. tuberosum | NC008096 | 155,296 | 37.88 | 131 | 36 | 4 | Gargano et al. (2005) |
*The data have been partially adopted from Park (2023b)..
S. jamesii (PI578326) 엽록체의 전장유전체 분석에는 S. jamesii 계통 중 하나인 SJ-8이 이용되었으며, S. jamesii 특이적 분자마커 개발을 위해 S. jamesii의 SJ-6과 SJ-11 계통이 추가되었다. 이외에 3개의 재배종 감자 품종과 12종의 가지속(Solanum) 이용되었으며, 자체적으로 대구에서 수집한 S. nigrum (SN)을 제외한 모든 식물재료는 국립식량과학원 고령지농업연구소에서 분양 받아 기내식물체로 유지 및 증식하여 이용하였다. 실험에 이용된 재배종 감자의 품종은 ‘대지’(DJ), 대서’(DS), ‘수미’(SM)이고, S. nigrum 이외 11개의 야생종 계통은 S. acaule (PI310970, SA), S. brevicaule (PI205394, SB2), S. hougasii (PI161174, SH2), S. cardiophyllum (PI341233, SC1), S. commersonii (PI558050, SC2), S. hjertingii (PI186559, SH1), S. stoloniferum (PI160224, SS), S. kurtzianum (PI578236, SK), S. mochiquense (PI338616, SM1), S. verrucosum (PI160228, SV1), S. vernei (PI230468, SV2)이다.
DNA 분리는 Genomic DNA Extraction kit (Plants)(RBC, New Taipei City, Taiwan)를 이용하여 모든 식물체를 대상으로 계대배양 한 후 약 3주 이내 유식물체의 약 100 mg을 채취하여 진행하였다.
S. jamesii의 엽록체 전장유전체 분석은 S. jamesii 계통 SJ-8을 대상으로 진행하였다. 전체 염기서열 정보(raw data)는 Macrogen (Macrogen, Seoul, South Korea)의 MiSeq (Illumina, SanDiego, CA, USA) 플랫폼을 통해 확보하였으며, 염기서열의 분석과 조립은 파이젠(http://phyzen.com)에서 제공하는 de novo assembly가 이용되었다(Cho et al. 2015). 그 과정에서 전처리(trimming), high-quality 염기서열 추출(Phred score 20 이상), dnaLCW 방법으로 조립 등은 CLC quality trim, CLC de novo assemble 등의 프로그램이 이용되었으며, 이들 프로그램은 CLC genome assembler package version 4.2.1 (CLC Inc, Rarhus, Denmark)에서 제공된 것이다. 이후 조립된 염기서열 중, 엽록체 유전체 유래의 contig를 선발하고, gap-filling에 의한 확장과 통합을 통해 엽록체 유전체를 완성하였다. 전처리된 염기서열 전체를 완성된 엽록체 전장 유전체 서열에 mapping하여 각각의 영역별 read depth를 계산하였고, reference 종과 유전체의 구조 및 염기서열을 비교하였다. 또한, 엽록체 전장 유전체 서열을 이용하여 BLASTN을 수행하였으며, homology가 높은 종의 엽록체 전장 유전체 서열을 확인하고 BLASTZ (Schwartz et al. 2003) 분석을 실시하였다.
완성된 S. jamesii 엽록체 전장 유전체의 유전자는 Geseq 프로그램을 통해 확인하고, BLAST search를 기반으로 한 manual curation을 통해 최종적으로 유전자와 유전자의 위치를 결정하였다(Tillich et al. 2017). 이들은 OGDraw 프로그램을 통해 원형의 유전자지도 상에 표시되었다(Lohse et al. 2007).
S. jamesii와 유사종들과의 유전관계 분석을 위해 NCBI (National Center for Biotechnology Information)에 등록되어 있는 재배종 감자 S. tuberosum (NC008096 및 KM489056), 가지과의 까마중 S. nigrum (KM489055), 고추 Capsicum annuum (JX270811), 토마토 S. lycopersicum (NC007898), 감자의 근연야생종 12종 S. acaule (MK036506), S. brevicaule (MK036507), S. berthaultii (KY419708), S. bulbocastanum (DQ347958), S. demissum (MK036508), S. cardiophyllum (MK690622), S. stoloniferum (MF471373), S. chacoense (MF471371), S. commersonii (KM489054), S. hjertingii (MK690623), S. hougasii (MF471372), S. verrucosum (MK690625) 등 총 18개의 엽록체 전체 염기서열을 확보하여, 여기에 존재하는 단백질 코딩 서열을 MEGA11 프로그램을 이용하여 분석하였으며, 이 프로그램은 1,000번의 Bootstrap 옵션이 적용되었고 GTR (General Time Reversible)모델 기반의 Maximum Likelihood 방법이 이용되었다(Tamura et al. 2021).
분자마커의 개발은 S. jamesii의 엽록체 전체 염기서열과 NCBI로부터 얻은 S. verrucosum (MK690625), S. cardiophyllum (MK690622), S. bulbocastanum (DQ347958), S. berthaultii (KY419708), S. hjertingii (MK690623), S. tuberosum (KM489056), S. commersonii (KM489054), S. hougasii (MF471372)의 엽록체 전체 염기서열을 ClustalW2 (EMBL: http://www.ebi.ac.uk/Tools/msa/clustalw2)를 이용하여 다중정렬 하였고, 이를 통해 S. jamesii에 특이적으로 나타나는 InDel과 SNP를 확보할 수 있었다. InDel의 경우 InDel이 있는 영역에 primer를 제작하여 allele 특이적 분자마커를 개발하고자 하였으며, SNP의 경우 tetra-primer ARMS (Amplification Refractory Mutation System)-PCR 기술을 이용하여 primer를 제작하였다(Collins and Ke 2012; Ye et al. 2001). S. jamesii와 S. tuberosum을 포함한 총 18개의 계통, S. jamesii (SJ6, SJ8 및 SJ11), S. tuberosum (감자품종 ‘대지’(DJ), ‘수미’(SM), ‘대서’(DS)), S. acaule (SA), S. cardiophyllum (SC1), S. brevicaule (SB2), S. commersonii (SC2), S. hjertingii (SH1), S. hougasii (SH2), S. mochiquense (SM1), S. kurtzianum (SK), S. nigrum (SN), S. stoloniferum (SS), S. verrucosum (SV1), S. vernei (SV2)를 대상으로 PCR을 수행하였다. InDel 기반의 분자마커 개발은 위한 PCR의 경우, 총 20 ul (약 20 ng genomic DNA, 0.5 mM dNTPs, 10 pMol each primer, 1 U Taq polymerase (Genetbio, Daejeon, South Korea))으로, SNP 기반의 분자마커 개발을 위한 PCR의 경우, 총 20 ul (약 20 ng genomic DNA, 0.5 mM dNTPs, 10 pMol each inner primer, 1 pMol each outer primer, 1 U Taq polymerase (Genetbio, Daejeon, South Korea))으로 Thermocycler (Biometra, Göttingen, Germany)을 이용하였다. PCR의 결과는 핵산 염색용액인 StaySafe (RBC, New Taipei City, Taiwan)를 이용하여 1.5% agarose gel에서 확인하였다.
S. jamesii의 엽록체 전체 염기서열은 NGS 기술을 이용하여 완성되었다. NGS 기반으로 얻은 데이터는 평균 301.0 bp 길이의 총 염기서열 약 1.5 Gbp로, 이는 5,267,110 read를 통해 생산되었다. 이중 약 90.79%를 이용하여 전처리 및 Phred score 20 이상의 데이터 추출과정을 거쳐 high-quality 총 염기서열 정보 약 1.3 Gbp를 얻었으며, 평균 길이는 283.6 bp로 확인되었다. 이는 전체 raw data의 약 85.55%였으며, S. cardiophyllum 의 전체 엽록체 염기서열과 비교했을 때 사용된 총 염기서열이 약 0.3 Gbp 적었지만, 총 생산된 데이터 대비 추출 염기서열의 비율은 비슷한 수준이었다(Park 2023a, 2023b). 또한, de novo assembly와 엽록체 유래의 contig를 선발하고, gap-filling 과정을 통한 확장과 통합 등을 거쳐 총 염기서열은 155,576 bp의 S. jamesii 엽록체 전장 유전체가 완성되다. 최종적으로 aligned read의 개수는 435,578개, 전체 유전체에 대한 평균 coverage는 약 789.54x, 전체 trimmed read를 완성된 전체 유전체 서열에 mapping한 결과로 나타난 read depth는 최소 378 이상이었다. 구조적으로는 원형 이중가닥 분자 구조로 다른 식물의 엽록체 유전체와 동일한 형태였으며, 한 쌍의 25,596 bp의 IRa/IRb, 18,394 bp의 SSC, 85,990 bp의 LSC 영역으로 구성되어 있었다. S. bulbocastanum의 엽록체 전체 염기서열(DQ347958)과의 비교와 BLASTZ 분석의 결과로도 그 유사성을 확인할 수 있었다(Daniell et al. 2006; Schwartz et al. 2003) (Fig. 1). S. jamesii의 엽록체 전체 염기서열을 대상으로 NCBI의 BLASTN 분석 결과 S. pinnatisectum (MZ233590), S. cardiophyllum (MK690622)과 99.95%의 높은 염기서열 일치도를 보였으며, GC 비율 또한 37.90%로 매우 유사하였다(Table 1).
S. jamesii의 엽록체 전장 유전체에는 105개의 단백질 코딩, 45개의 tRNA, 8개의 rRNA를 포함는 총 158개의 유전자가 존재하고 있는데, 이중 11개의 단백질 코딩 유전자, 9개의 tRNA, 4개의 rRNA가 한 쌍의 Ira 및 IRb 영역에 역배열로 중복되어 있었다(Table 1, Fig. 2). 전체 염기서열 중 코딩 영역의 비율은 59.2%로 나타났으며, 158개 유전자의 평균 길이는 약 583.4 bp, 이중 각각의 평균은 단백질 코딩 유전자가 764.8 bp, tRNA가 62.4 bp, rRNA가 1,131.0 bp의 크기로 각각 51.6%, 1.8%, 5.8%의 비율로 분포되어 있으며, 이러한 결과는 다른 Solanum 종들과 매우 유사한 수준이었다(Table 1).
S. jamesii의 계통수 작성은 총 17종 가지과의 엽록체 전장 유전체 코딩 서열을 기반으로 실행되었다(Fig. 3). 그 결과, S. jamesii는 S. tuberosum 뿐만 아니라 다수의 가지속(Solanum) 종들과 같은 그룹에 존재함을 확인하였으며, 유연관계가 그중 S. cardiophyllum과 가장 가까웠으며, 상대적으로 까마중(S. nigrum), 고추(C. annuum), 토마토(S. lycopersicum)와는 다소 먼 것으로 확인되었다.
S. jamesii 엽록체의 전체 염기서열을 7종의 근연야생종과 재배종 감자의 엽록체 전체 염기서열과 ClustalW2 (EMBL; https://www.ebi.ac.uk/Tools/msa/clustalw2)를 이용하여 다중 정렬한 결과 S. jamesii 특이적인 2개의 InDel과 11개의 SNP를 구명하였다(Fig. 4). 상대적으로 S. jamesii 특이적 InDel과 SNP의 수가 매우 적은 편이었는데, 이는 다중 정렬에 S. jamesii의 엽록체 전체 염기서열과 유사도가 매우 높은 S. cardiophyllum이 포함되어 나타난 결과로 보인다. 하지민, 다양한 유전적 요인에 의한 변이 발생으로 인해 S. jamesii 특이적인 InDel과 SNP가 존재하며, 이를 이용하여 S. jamesii 특이적 분자마커를 개발할 수 있었다(Cho and Park 2016; Chung et al. 2006; Kim et al. 2018; Kim and Park 2019, 2020a, 2020b; Park 2021, 2022a, 2022b, 2023a, 2023b). 2개의 InDel은 각각 유전자 간 및 내 영역에 존재하고 있었으며, 11개의 SNP 중 8개는 유전자 간 영역에 3개는 유전자 내 영역에 위치하고 있었다.
S. jamesii를 포함한 총 9개 Solanum 종의 엽록체 전체 염기서열의 다중 정렬 결과로 확인된 S. jamesii 특이적 InDel과 SNP는 S. jamesii 특이적인 분자마커 개발에 활용되었다. InDel은 InDel 영역에서 특이적인 primer를 제작하여 PCR을 통해 다양성의 차이가 나타나는 S. jamesii 특이적 또는 비특이적인 3개의 SCAR (Sequence Characterized Amplified Region) 마커를 개발할 수 있었다(Cho et al. 2015; Garcia-Lor et al. 2013; Yamaki et al. 2013)(Table 2, Fig. 5A, 5B, 5C). SNP는 11개 중 3개의 영역을 대상으로 S. jamesii 특이적 SNP 위치를 기준으로 SNP가 포함된 DNA 단편과 SNP 위치을 대상으로 S. jamesii 특이적이거나 혹은 S. jamesii를 제외한 Solanum 종들에 특이적으로 DNA를 증폭할 수 있는 tetra-primer ARMS-PCR 방법을 적용하여 primer를 제작하였다(Collins and Ke 2012; Ye et al. 2001)(Table 2). 이러한 방법은 PCR로 염기서열이 매우 유사한 SNP를 대상으로 primer의 제작이 어려운 종 또는 SNP가 존재하는 계통 간의 차이를 비교적 손쉽게 확인할 수 있다는 장점이 있어 널리 활용되고 있으며, 본 연구에서도 Solanum 종간의 비교 결과로 나타난 SNP를 활용한 PCR 기반의 S. jamesii 특이적이거나 (Fig. 5D, 5E, 5F) 혹은 비특이적으로(Fig. 5D, 5E, 5F) 나타나는 3개의 분자마커를 개발하였다(Han et al. 2016; Kou et al. 2017; Mahmoudi et al. 2019; Zou et al. 2018). 최종적으로 본 연구에서는 InDel과 SNP 정보를 기반으로 한 S. jamesii 특이적 분자마커가 각각 3개씩 개발되었다. SJ8_InDel1 (Fig. 4, Fig. 5A, 5B)과 SJ8_InDel2 (Fig. 4, Fig. 5C)은 각각 trnL 유전자 내 및 petA와 psbJ 유전자 사이에 존재하는 SCAR 마커로 개발되었으며, SJ8_SNP1 (Fig. 4, Fig. 5D), SJ8_SNP3 (Fig. 4, Fig. 5E)과 SJ8_SNP5 (Fig. 4, Fig. 5F)는 각각 atpF 유전자 내, atpH와 atpI 유전자 사이 및 psaA 유전자 내에 존재하는 SNP 기반의 마커로 개발되었다.
Table 2 . Information of primers used to generate S. jamesii-specific and non-specific markers.
InDel/SNP | Strand | Primer sequence (5’ → 3’) | Tm (°C) | Size |
---|---|---|---|---|
SJ8_InDel1(a) | Forward primer | AAACGTATATACATACTATATC | 56°C | 735 bp (delete sequence 11bp specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel1(b) | Forward primer | AAACGTATATACATACGTATTG | 58°C | 746 bp (delete sequence 11bp non-specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel2 | Forward primer | AGTATAGAAAGGATTAGATTAG | 60°C | 777 bp (insert sequence ‘ATTAGATTAGTATAGAAAGG’ specific to S. jamesii) |
Reverse primer | TATCAGCAATGCAGTTCATC | |||
SJ8_SNP1 | Forward inner primer | CCTATCTTCTCTTCTTTGTTCATATTTCC | 60°C | 396 bp (from two outer primers), 217 bp (C allele specific to S. jamesii), 238 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTGGTCTTACATTAGTTTCGATATCTTCTT | |||
Forward outer primer | TCAAGATCCTCTGTTTTCGATTATCTAA | |||
Reverse outer primer | GTAAATTGTAAAAACAGGGATTTTACCC | |||
SJ8_SNP3 | Forward inner primer | AATTCTCTTGAGCCACCCACCATATGCC | 56°C | 300 bp (from two outer primers), 178 bp (C allele specific to S. jamesii), 210 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTTTAGCTTAGCCCACAGCAATGTATAGAT | |||
Forward outer primer | TGAAGGAAAGGGGAATTTTAGGAAAAAGAT | |||
Reverse outer primer | GTTCCTATACCTGTCATGCTCCTTGGATT | |||
SJ8_SNP5 | Forward inner primer | CTGAATAGAAACATTAAACTAAAAGCCTAG | 58°C | 390 bp (from two outer primers), 243 bp (T allele specific to S. jamesii), 204 bp (G allele non-specific to S. jamesii) |
Reverse inner primer | CTTTTTTTCCTAGGTGCTCATTTTATA | |||
Forward outer primer | TACTGCAATAATTCTTGCTAAGAAGAAT | |||
Reverse outer primer | GTAAGTGATCAAGGGGTAGTAACTCATA |
감자의 신품종 육성을 위해 다양한 유용한 특성을 지니고 있는 근연야생종을 이용한 연구가 계속해서 수행되고 있으나, 교배육종을 통한 유용 형질의 도입은 다수의 근연야생종이 재배종 감자와의 EBN과 배수성 차이로 인해 어려움을 겪기 때문에, 체세포융합을 이용한 체세포잡종의 육성이 활용되고 있으며, 이에 따라 유전자형 확인이 필수적이다. 분자마커를 통한 유전자형의 확인은 단순히 염색체 조합만을 확인하는 것이 아니라 엽록체 및 미토콘드리아 유전체의 전달 확인도 필요하다. 특히 Solanum 속에 속하는 다양한 식물종의 경우, 체세포융합이나 재분화 과정에서 핵 내 및 미토콘드리아 유전체는 재조합 빈도가 높다는 것이 알려져 있으며, 엽록체 유전체의 경우에는 한 쪽 친의 엽록체 유전체가 무작위적으로 전달되는 것 외에도 양친의 유전체가 모두 전달되는 경우도 있었다(Chen et al. 2013; Cho et al. 2016; Lössl et al. 2000; Mohapatra et al. 1998; Smyda-Dajmund et al. 2016; Wang et al. 2011; Xiang et al. 2004). 이에 따라, 본 연구에서 개발된 분자마커는 체세포융합을 통한 체세포잡종의 육성 과정에서 특히 S. jamesii와 S. tuberosum 유래의 엽록체 유전체의 유전자형 구별에 이용이 될 것이며, 이러한 분자마커는 감자 육성 과정에서 신품종의 개발 및 진화학적 연구에 기여할 수 있을 것이다.
멕시코 북서부부터 미국 남서부까지 자생하는 감자 야생종 중 하나인 Solanum jamesii는 Phytophthora infestans에 대한 저항성으로 감자 육종에 훌륭한 재료이다. 그러나 S. jamesii는 이배체이며 EBN 값이 1로, 이는 S. jamesii와 재배용 감자 간 교배가 불가능하다. 이 문제를 극복하기 위해 체세포융합을 이용할 수 있다. 이후에는 분자마커를 사용하여 효율적으로 융합체를 선발하는 것이 필수적이다. 본 연구에서는 차세대 유전체 기술을 사용하여 S. jamesii의 엽록체 전장 유전체 서열을 확보하고, 이를 다른 8개의 가지과에 속한 종과 비교하여 S. jamesii 특이적 마커를 개발하였다. S. jamesii의 엽록체 유전체의 총 서열 길이는 155,576 bp이며, 그 구조는 다른 Solanum 종과 유사하다. 17개의 다른 Solanaceae 종과의 계통수 분석 결과, S. jamesii는 S. cardiophyllum과 가장 밀접한 유연관계를 가졌다. 8개의 다른 Solanum 종의 엽록체 전장 유전체 서열 비교에서 S. jamesii에 특이한 2개의 InDel 및 11개의 SNP가 확인되었다. 이러한 InDel 및 SNP를 기반으로 S. jamesii를 다른 Solanum 종과 구별하는 6개의 PCR 기반 분자마커가 개발되었다. 본 연구의 결과는 Solanum 종 중에서 S. Jamesii를 구별하고 S. Jamesii를 활용한 감자 신품종 개발에 기여할 것이다.
Table 1 . Results of comparative analysis of the chloroplast genome sequence of S. jamesii with those of 14 Solanaceae species.
Species | Accession no. | Total length (bp) | GC content (%) | Total no. of genes | No. of tRNAs | No. of rRNAs | Reference |
---|---|---|---|---|---|---|---|
S. jamesii | MK690624 | 155,576 | 37.90 | 134 | 36 | 4 | This study |
S. hjertingii | MK690623 | 155,545 | 37.88 | 134 | 36 | 4 | Park (2023a) |
S. cardiophyllum | MK690622 | 155,570 | 37.90 | 134 | 36 | 4 | Park (2023b) |
S. acaule | MK036506 | 155,570 | 37.84 | 135 | 36 | 4 | Park (2022b) |
S. brevicaule | MK036507 | 155,531 | 37.87 | 135 | 36 | 4 | Park (2022a) |
S. demissum | MK036508 | 155,558 | 37.87 | 135 | 36 | 4 | Park (2021) |
S. hougasii | MF471372 | 155,549 | 37.87 | 135 | 36 | 4 | Kim and Park (2020b) |
S. stoloniferum | MF471373 | 155,567 | 37.87 | 135 | 36 | 4 | Kim and Park (2020a) |
S. chacoense | MF471371 | 155,532 | 37.89 | 136 | 36 | 4 | Kim and Park (2019) |
S. berthaultii | KY419708 | 155,533 | 37.88 | 137 | 39 | 4 | Kim et al. (2018) |
S. commersonii | KM489054 | 155,525 | 37.88 | 133 | 33 | 4 | Cho et al. (2016) |
S. nigrum | KM489055 | 155,432 | 37.90 | 139 | 39 | 4 | Cho and Park (2016) |
S. tuberosum | KM489056 | 155,312 | 37.88 | 130 | 30 | 4 | Cho et al. (2016) |
S. bulbocastanum | DQ347958 | 155,371 | 37.88 | 133 | 30 | 4 | Daniell et al. (2006) |
S. tuberosum | NC008096 | 155,296 | 37.88 | 131 | 36 | 4 | Gargano et al. (2005) |
*The data have been partially adopted from Park (2023b)..
Table 2 . Information of primers used to generate S. jamesii-specific and non-specific markers.
InDel/SNP | Strand | Primer sequence (5’ → 3’) | Tm (°C) | Size |
---|---|---|---|---|
SJ8_InDel1(a) | Forward primer | AAACGTATATACATACTATATC | 56°C | 735 bp (delete sequence 11bp specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel1(b) | Forward primer | AAACGTATATACATACGTATTG | 58°C | 746 bp (delete sequence 11bp non-specific to S. jamesii) |
Reverse primer | TTTCAGTCCTCTGCTCTAC | |||
SJ8_InDel2 | Forward primer | AGTATAGAAAGGATTAGATTAG | 60°C | 777 bp (insert sequence ‘ATTAGATTAGTATAGAAAGG’ specific to S. jamesii) |
Reverse primer | TATCAGCAATGCAGTTCATC | |||
SJ8_SNP1 | Forward inner primer | CCTATCTTCTCTTCTTTGTTCATATTTCC | 60°C | 396 bp (from two outer primers), 217 bp (C allele specific to S. jamesii), 238 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTGGTCTTACATTAGTTTCGATATCTTCTT | |||
Forward outer primer | TCAAGATCCTCTGTTTTCGATTATCTAA | |||
Reverse outer primer | GTAAATTGTAAAAACAGGGATTTTACCC | |||
SJ8_SNP3 | Forward inner primer | AATTCTCTTGAGCCACCCACCATATGCC | 56°C | 300 bp (from two outer primers), 178 bp (C allele specific to S. jamesii), 210 bp (A allele non-specific to S. jamesii) |
Reverse inner primer | TTTTAGCTTAGCCCACAGCAATGTATAGAT | |||
Forward outer primer | TGAAGGAAAGGGGAATTTTAGGAAAAAGAT | |||
Reverse outer primer | GTTCCTATACCTGTCATGCTCCTTGGATT | |||
SJ8_SNP5 | Forward inner primer | CTGAATAGAAACATTAAACTAAAAGCCTAG | 58°C | 390 bp (from two outer primers), 243 bp (T allele specific to S. jamesii), 204 bp (G allele non-specific to S. jamesii) |
Reverse inner primer | CTTTTTTTCCTAGGTGCTCATTTTATA | |||
Forward outer primer | TACTGCAATAATTCTTGCTAAGAAGAAT | |||
Reverse outer primer | GTAAGTGATCAAGGGGTAGTAACTCATA |
Ju-Ryeon Jo ・Tae-Ho Park
J Plant Biotechnol 2024; 51(1): 158-166Seoyeon Son ・Tae-Ho Park
J Plant Biotechnol 2024; 51(1): 121-128Tae-Ho Park
J Plant Biotechnol 2023; 50(1): 45-55
Journal of
Plant Biotechnology