Warning: mysqli::mysqli() [function.mysqli-mysqli]: (00000/1040): Too many connections in /home/kjbs/www/inc/db_con.inc on line 9

Warning: mysqli::query() [function.mysqli-query]: Couldn't fetch mysqli in /home/kjbs/www/inc/db_con.inc on line 10

Warning: main() [function.main]: Couldn't fetch mysqli in /home/kjbs/www/inc/db_con.inc on line 13
Assessment of the Candidate Genes of Expression Markers Associated with Drought Stress in Maize Seedlings
search for




 

Assessment of the Candidate Genes of Expression Markers Associated with Drought Stress in Maize Seedlings
옥수수 유묘기 한발 스트레스 평가를 위한 발현마커 후보군 탐색
Korean J Breed Sci 2018;50(3):224-235
Published online August 31, 2018
© 2018 Korean Society of Breeding Science.

Jun-Cheol Moon1, Seungho Shin2, Hyo Chul Kim2, Kitae Song2, Jae Yoon Kim3, Kyung-Hee Kim2, and Byung-Moo Lee2,*
문준철1, 신승호2, 김효철2, 송기태2, 김재윤3, 김경희2, 이병무2,*

1Agriculture and Life Sciences Research Institute, Kangwon National University, Chuncheon, 24341 Korea,
2Department of Life Science, Dongguk University-Seoul, Seoul, 04620, Korea,
3Department of Plant Resources, Kongju National University, Yesan, 32439, Korea
1강원대학교 농업생명과학연구원,
2동국대학교 생명과학과,
3공주대학교 식물자원학과
Correspondence to: (E-mail: bmlee@dongguk.edu, Tel: +82-31-961-5130, Fax: +82-31-961-5108)
Received August 8, 2018; Accepted August 14, 2018.
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract

Drought stress during the seedling stage has a disastrous effect on the growth of maize. The purpose of this study was to assess the developed expression markers that are related to drought stress in maize. For the selection of expressed genes by drought stress, co-expression analysis was carried out using published microarray data of drought stress in maize (Zea mays L.) seedlings. Six consensus modules were based on 4,770 stress responsive genes differentially expressed in drought stress, and the royal blue module was chosen. Thirty genes were selected according to different expression patterns between susceptible and tolerant types. Drought stress treatments were performed on both Ki3 and Ki11. Ki3 and Ki11 are widely known drought-susceptible and -tolerant types, respectively. At first, the 30 selected genes were compared to Ki3 and Ki11 using qRT-PCR. The gene expression values of eight genes (BU050895, BF728598, CK827168, CO524848, AF457983, CF037152, AJ606944, and BG836522) were significantly tolerant types rather than susceptible types in the roots. After applying the eight above-mentioned genes to nine cultivars, a different pattern was detected between susceptible and tolerant types. The results of the present study will show the possibility of developing novel expression markers and the application for various maize varieties.

Keywords : Co-expression, Drought stress, Maize seedling, Susceptible type, Tolerant type
서언

옥수수(Zea mays L.)는 세계적으로 3대 작물 중 하나로써, 식량뿐만 아니라 사료와 산업용 원자재에도 널리 이용되고 있으며, 생육기간 동안 다양한 생물학적 및 비생물학적 스트레스 요인들로 인해 생산 및 수량에 영향을 받는다(Betrán et al. 2003, Tollenaar & Lee 2002). 옥수수는 다른 작물들에 비해 키가 큰 편이고 잎이 넓기 때문에 특히 한발에 더욱 치명적이며, 생산성 손실에 가장 중요한 요인 중 하나이다. 또한 한발은 옥수수 유묘기, 개화기 직전, 개화기의 anthesis-silking interval(ASI), 잎의 노화등에 영향을 미칠 수 있다(Leach et al. 2011).

한편 많은 연구들에서 옥수수 전체 생육기간 중 개화기가 한발 스트레스에 가장 민감한 시기라고 생각하기 때문에 중점적으로 다양한 연구가 진행되고 있다(Buckler et al. 2009, Chardon et al. 2004, Thornsberry et al. 2001). 하지만 한발 스트레스는 계절을 불문하고 어떤 시기라도 영향을 줄 수 있어 개화기에 집중하는 연구뿐만 아니라 유묘기의 한발 연구도 필요하다. 옥수수의 유묘기는 식물의 생존 자체에 영향을 주기 때문에 중요성으로는 개화기와 동일하다고 볼 수 있다(Banziger et al. 2000, Leach et al. 2011). 특히 유묘기는 식물체의 구조나 식물집단의 원동력을 결정하는 중요한 생육단계 이기 때문에 유묘기에서 형태학적으로 평가를 하는 것은 매우 매력적인 접근이라고 볼 수 있다. 또한, 유묘기에서 평가를 수행하면 처리할 수 있는 재료도 매우 풍부하고 비용도 절감되며, 공간 절약 및 시간을 단축시키는 장점을 가지고 있다(De La Cruz et al. 2008, Meeks et al. 2013). 유묘기 스트레스에서 강하고 약한 품종들을 선발하는 또 다른 장점은 형태적 다양성으로 인해 발생되는 실험적 오류를 제어할 수 있다는 점이다. 배양실이나 실험실에서 실험을 진행할 때 유묘기의 식물들은 거의 유사한 외형적 형태를 가지고 있기 때문에 생육 초기의 발달 이상을 빠르게 잡아낼 수 있다(Wang et al. 2015).

옥수수의 유전체는 7,000만년 전부터 유전체 일부분의 복제가 이루어 졌다고 알려졌으며, 유전체 전체의 복제는 1,200만년 전부터 500만년 전에 일어났다(Blanc & Wolfe 2004, Paterson et al. 2004, Swigonova et al. 2004). 이로 인해 가깝게 연관되어 있던 작물인 수수로부터 구분되었다(Paterson et al. 2009). 이로 인해 옥수수 유전체 중 10개의 염색체는 구조적으로 다양해지고, 염색질 조성에서 역동적인 변화를 겪었다. 옥수수 연구에 어려움이 있다면 최근 300만 년간 교배를 통한 long terminal repeat retrotransposons (LTR retrotransposons)의 급격한 증가로 인해 유전체 크기가 타 작물보다 크다는 점이다(SanMiguel et al. 1998). 대표적인 모델 식물인 애기장대는 유전체의 크기가 125 Mbp데 반해 옥수수는 유전체의 크기는 20배나 더 큰 2,500 Mbp이다(Arumuganathan & Earle 1991). 이는 벼의 유전체 크기에 6배이기도 하다. 이러한 옥수수의 매우 큰 유전체를 연구하고 내병성 및 내재해성 품종을 개량하기 위한 다양한 마커들을 개발함에 있어 많은 양의 데이터베이스와 선별 마커들이 필요하다(Myles et al. 2009). 최근에 옥수수 연구의 어려움을 극복하기 위해 nested association mapping (NAM) population을 이용하여 분자 육종 연구가 진행되고 있다(Yu et al. 2008). NAM population은 B73을 공유하고 있는 교잡종 집단으로, B73의 염기서열 해독은 완료되었으며(Schnable et al. 2009), B73은 생물학적, 비생물학적 스트레스에 중간 또는 약한 성질을 띄고 있다. 이를 이용하여 많은 스트레스 내성 품종을 육종하기 위해 모본 또는 부본으로 사용되고 있다(Hao et al. 2010, Krill et al. 2010).

옥수수를 연구 및 이해함에 있어 지식적 한계를 극복하기 위해 방대한 양의 데이터베이스를 이용하여 생물정보학 기술 및 데이터베이스 기반 비교 유전체연구가 활발하게 진행되고 있다. 옥수수를 포함한 벼, 수수, 밀, 보리 등의 작물에 대한 유전체 비교분석은 전반적인 유전자 지도 수준에서 유전자의 정보와 그 염기서열이 보존되고 있다는 것을 보여준다(Gale & Devos 1998). 따라서 공개된 유전자 정보 및 염기서열을 가지고 많은 연구진이 여러 연구에 적용하고 있다. 옥수수의 경우 Long intergenic noncoding RNAs (lincRNAs)에 관련된 유전자 발현연구(Fan et al. 2015), 600k SNP genotyping array를 이용한 유전체 분석방법(Unterseer et al. 2014), 침수 스트레스에 관련된 유전체 전반적인 전사체 발현 연구(Thirunavukkarasu et al. 2013)등이 있다. 한발스트레스 관련 옥수수 연구에도 데이터베이스를 이용한 연구가 활발하게 진행되고 있다(Gulli et al. 2015, Zhang et al. 2014).

본 연구는 NCBI의 옥수수 유묘기 한발 관련 microarray 데이터를 분석하여 한발 내성 계통과 감수성 계통의 발현패턴 비교분석을 통해 유전자를 선발하고, 다양한 품종에 적용하여 한발에 대한 내성과 감수성을 선별하는 것을 목표로 한다. 궁극적으로 옥수수 한발 내성 품종 육성을 위하여 유묘기에 한발 내성 품종을 선발할 수 있는 유전자 기반 발현 마커 개발에 목적이 있다.

재료 및 방법

생육 조건 및 한발스트레스 처리

실험재료는 B73, CML103, CML228, CML277, CML322, Ki3, Ki11, M37W, NC358 총 9개의 옥수수 자식 계통을 USDA-ARS에서 분양 받아 사용하였다(Table 1). 옥수수 계통의 생육을 위해 지름 130 mm, 높이 180 mm, 바닥 80 mm의 포트에 토양은 질석(vermiculite)을 선택하여 종자를 파종하였다. 상토로 사용되는 질석은 공기 및 수분 친화력이 높기 때문에 뿌리샘플을 얻는 작업에서 뿌리에 붙은 토양을 제거하기가 용이하다. 각각의 포트마다 질석을 10 L씩 채웠으며, 복토의 깊이는 4 cm로 하였다. 수분공급은 각 포트마다 매일 1L씩 공급하였으며, 생육환경은 온도 25°C, 습도 40%가 항상 유지되며 16시간 광조건(LED 램프 및 삼파장 램프, 9500 lux), 8시간 암조건이 일정하게 반복되는 식물배양실 안에서 진행하였다.

List of 9 inbred lines of maize in study.

Inbred lineOrigin
B73USA, Iowa
CML 103Mexico, Federal District
CML 228Mexico, Federal District
CML 277Mexico, Federal District
CML 322Mexico, Federal District
Ki3Thailand
Ki11Thailand
M37WSouth Africa, KwaZulu-Nata
NC358USA, North Carolina

옥수수 유묘기(V2 stage, 발아 후 2주)에 우선적으로 대조군 식물체의 잎과 뿌리에서 각각 RNA추출을 진행하였다. 한발 스트레스 처리군은 뿌리가 손상되지 않도록 물을 이용하여 질석을 제거한 뒤 6시간 및 12시간 동안 페이퍼 타올 위에서 한발 처리를 진행하였다. 한발 처리가 끝난 유묘기 옥수수의 잎과 뿌리에서 각각 RNA를 추출하였다. 물을 재공급한 처리군은 6시간, 12시간 동안 한발 처리를 진행한 유묘기 옥수수에 각각 24시간 동안 뿌리가 잠길 정도로 물을 채운 뒤 뿌리 부분이 빛에 노출되지 않도록 페이퍼 타올과 호일을 덮어 수분을 유지하였다. 물 재공급 처리가 끝난 유묘기 옥수수에서도 잎과 뿌리에서 RNA를 추출하였다.

한발 스트레스 처리는 세번 수행되었으며, RNA은 세번의 처리에서 각각 추출되어 cDNA합성 후 qRT-PCR분석에 사용하였다.

RNA 추출 및 cDNA 합성

RNA 추출은 TRIzol Reagent (Life technologies)를 이용하여 추출하였다. 대조군, 한발 처리 및 물 재공급 처리군의 샘플을 액체질소에 동결시킨 후 막자와 막자사발을 이용하여 곱게 분쇄하였다. 분쇄된 뿌리 조직을 2.0 ml tube에 각각 75 mg씩 넣고 TRIzol Reagent를 1 ml을 첨가하여 균질하게 섞어주었으며, 이후 과정은 제조사의 매뉴얼에 따랄 수행되었다. 추출한 RNA는 nano drop (MAESTROGEN)을 이용하여 정량하였다. 정량된 RNA는 PrimeScript™ 1st strand cDNA Synthesis Kit (TaKaRa)를 이용하여 메뉴얼에 따라 cDNA를 합성하였다.

In silico 분석

Co-expression network 분석을 위하여 NCBI GEO에서 약 1,000개의 microarray에 대한 CEL 파일들을 다운로드하였다. 수집된 CEL 파일을 R package affy를 통하여 RMA 방식으로 규격화하여 통합하였으며, 1,000개의 CEL 파일들을 통하여 구축된 microarray 데이터베이스는 R package the weighted correlation network analysis (WGCNA)를 통하여 co-expression network 분석하였다(Langfelder & Horvath 2008). 내성 타입과 감수성 타입과의 비교 분석을 통해 옥수수에 한발 스트레스 관련 유전자들을 선발하고, 내성 타입과 감수성 타입에서 유전자들의 발현랑을 MeV program ( http://www.tm4.org/mev.html)을 사용하여 분석하였다. 분석된 모듈은 DAVID ( http://david.abcc.ncifcrf.gov/home.jsp)를 이용하여 functional enrichment analysis을 실시하여 발현양상을 heatmap과 유전자의 수를 venn diagram으로 표시하였다. 클러스터링을 통해 내성 타입과 감수성 타입 사이의 consensus한 모듈을 구축하였다.

Primer 제작 및 qRT-PCR 분석

Primer는 royal blue 모듈에 포함된 유전자들 중 물 재공급 처리군에서 한발 내성타입과 한발 감수성타입 간 발현 차이를 보이는 유전자 상위 30개를 선발하였다(Table 2). Reference gene으로는 다양한 환경에서도 가장 안정적으로 발현되는 housekeeping gene인 18s rRNA를 선택하였다(Manoli et al. 2012). 18s rRNA의 forward primer는 5’-CCATCCCTCCGTAGTTAGCTTCT-3’이며, reverse primer는 5’-CCTGTCGGCCAAGGCTATATAC-3’이다. 한발 처리의 확인하기 위해 internal control로 Dehydration Responsive Element Binding proteins(DREB) 유전자를 이용하였다(Liu et al. 2013). 본 실험에서는 ZmDREB1A의 subgroup인 ZmDREB1.7을 이용하였다. ZmDREB1.7의 forward primer는 5’-CAGCGACTGCTGGCAAAGC-3’이며, reverse primer는 5’-CAGTAGCTCCAGAGCGTCATG-3’이다.

Information of primer sequences.

NCBI Accession no. Genechip ID Forward primer (5’ → 3’)Annealing Temp.(°C)

 Reverse primer (5’ → 3’)
BU050895Zm.12616.1GCTAATTCCGCAGGAGAGAAA56
TGCCAGCCCTTCTTCTTAAC
CO527263Zm.43.1CTTGTTCGCGAGTCCAGTAT56
CAGCAACGTACCTGCAAATC
AI834513Zm.2073.2CGTGGGTTTGGCTTTGTTAC56
CATTCAGAAGCCTCCCACTAAT
BF728598Zm.1769.1AACCAAGATCGAGCGAATCC56
CATGCGGTTTATTGGCATACAA
CK827168Zm.1987.2ACGGCATGGACACAACTAAA56
ACGGTCCAACTCCTTTGATATG
CD978285Zm.2321.4GCGAAAGTAGTGGCAGAGAA56
CATGCAGGTGCAGATGATAGA
AY104757Zm.14392.1CTCAAGGGCCCACTGTATTATG60
CTCCCAACAACATCGGCTATAC
CO524848Zm.18537.1GAGAGGTCGCGTAGATGATTAG56
TAGGTATTAGCTGCCGATGC
AF457983Zm.13389.1GTTCTTTCTTGCTGTTGCTTAGT57
GATTCCACACTCTCTTCCCTTC
CF627469Zm.2139.1CGTGGGTTTGGCTTTGTTAC56
CATTCAGAAGCCTCCCACTAAT
BG841205Zm.2171.1CGTGGGTTTGGCTTTGTTAC56
CATTCAGAAGCCTCCCACTAAT
BM072816Zm.8068.1GAGAGTCAAAGGGAAGGTGAAG56
ACTTGTTGCCTGACGAGAAATA
CF631845Zm.17625.1ACGAAGCTCGTCGAAGAATATG56
TTGCTCACAGGCGAACTAAA
CF037152Zm.15086.1GGGTAACCAGGATGTGTCTTT56
ACCCTGGAGATGGATTCTTATTG
AY108162Zm.4489.1ATGTCACTCATGACGACTTCTAC56
TCCGTCCTTCCTGCATTAAC
AJ606944Zm.16971.1CCTGTCTTTGCTTCCCTCTTAC53
AACGTGCCCAGCTTCTTT
CO520926Zm.19277.1TCTACTATCTGCCCGTCTACAA58
GAAGTAGTGGAAGCACCAGTAG
BM340313Zm.9828.1GATGAACGAGCGAGGATTGA56
CATAGCAAGTATACCGGAGACC
AW927712Zm.2258.1CGTGGGTTTGGCTTTGTTAC56
CATTCAGAAGCCTCCCACTAAT
BM381236Zm.11849.1CGTGTGCTAGAAAGAACCTGTA56
ACAGAACGAAGACGGAACAA
BG836522Zm.1967.1ACGGCATGGACACAACTAAA56
ACGGTCCAACTCCTTTGATATG
CA403577Zm.4270.2CGAGATGCAAGGTACAAGAAGA57
CCTGAACAGGCAGACCATATT
CB351077Zm.14663.1GTCATCTGGTCGAGTGGAAATAA59
GTCACAGATCGAAGAACCCTAAA
CD219246Zm.16973.1CCTGTCTTTGCTTCCCTCTTAC53
AACGTGCCCAGCTTCTTT
BM078934Zm.15119.1GTGCTGTGGAGCAAGGTAAA56
CCTCATCACTGCCAATCTCTTC
AY103786Zm.281.1GCCACCACCGAAGAAGATAAA57
GAAACGTAAGCGGAGGATCTG
CK145306Zm.17181.1GAGCAGAGAGCAGCACATTA56
CATCTGTAGCACCGCAAATTC
AY110813Zm.3423.1GGGAATTATCTCTCAGCTTCCC57
GCTTGGTCTCTGCAGTCTTAT
BM335820Zm.2227.1CGTGGGTTTGGCTTTGTTAC56
CATTCAGAAGCCTCCCACTAAT
AF329371Zm.445.1TTCTGCGATATGGTCCACAAG56
TTCCCAGCGAGTCAGATAGA

qRT-PCR은 CFX Connect Real-Time PCR Detection System (Bio-Rad)을 이용하였다. Real-Time용 pcr tube에 1/4로 희석한 cDNA 1μl, primer 10pmol, RNase free water 8μl, SYBG mix (Enzynomics) 10μl를 넣었다. qRT-PCR의 조건은 95°C에서 10분간 유지시킨 다음, 45 cycle (95°C 10초, 56°C 10초, 72°C 30초)을 진행한 후 melting curve를 65°C에서 95°C까지 진행하여 상대적 발현 값을 측정하였다. Bio-Rad CFX Manager 프로그램을 이용하여 melting curve 분석을 통해 primer dimer생성 여부를 확인하였다. 또한 qRT-PCR 결과 Cq값 및 Cq평균값을 산출했으며, ΔΔCt method[2^ (Ct sample − Ct ctrl.)]를 이용하여 각각의 상대적 발현값을 확인하였다. 세번의 반복을 통해 산출한 데이터를 바탕으로 발현 패턴을 분석하였다.

결과 및 고찰

In silico 분석

옥수수에서 한발 스트레스에 반응하는 4,771개의 유전자를 바탕으로 6개의 모듈(yellow, green yellow, turquoise, royal blue, brown4, plum1)을 구축하였다. 내성 타입에서는 turquoise, brown5 모듈이, 감수성 타입에서는 green yellow, turquoise, royal blue, brown4 모듈이 drought stress에서 대응관계의 모듈 구조를 보였다(Fig. 1). 클러스터링을 통해 내성 타입과 감수성 타입 간의 consensus한 모듈을 구축하였으며 이를 토대로 클러스터링 트리를 확인하였다. Heatmap을 분석한 결과 6개의 모듈에서 각각 유사한 발현 패턴을 보였고, 모듈별로 yellow 318개, green yellow 2,433개, turquoise 375개, royal blue 183개, brown4 1,405개, 그리고 plum1 56개의 유전자들이 포함되었다. 6개의 모듈 중 royalblue 모듈에서 한발에서는 대조군보다 증가하는 패턴을 보였지만 물 재공급 과정에서 내성 타입에서는 발현이 감소하고 감수성 타입에서는 발현이 증가하는 경향을 보였다. Royal blue 모듈에 포함된 유전자 183개 중에서 110개의 유전자가 내성 및 감수성 타입 모두 물 재공급 처리군에서 발현이 증가하였으며, 45개의 유전자는 감소하는 패턴을 보였다. 또한 4개의 유전자는 내성 타입에서 증가하였고, 24개의 유전자 감수성 타입에서 증가하여 내성 타입과는 반대의 경향을 보였다(Fig. 2, Table 3).

Fig. 1.

Clustering dendrogram of genes and consensus modules found in drought stress-related genes of maize. The correspondence between consensus modules and modules found individually in tolerant and susceptible types based on the expression values of the common genes. Each row corresponds to a drought stress specific module (labeled by color as well as text along with the number of genes in the module), and each column corresponds to one consensus module. Numbers indicate gene counts in the intersection of the corresponding modules. The stronger red color indicates a more significant overlap (left: tolerance type, right: susceptible type).


Number of genes differentially expressed in 6 modules.

ModuleTolerance typeSusceptible typeOverlap

MDSDRWMDSDRWMDSDRW
Yellowup11288949190105858784
down206230224227228213200227203
Green-yellowup9579741,2699569421,228880885907
down1,4761,4591,1641,4771,4911,2051,4001,402843
Turquoiseup191172155198170155130107148
down184203220177205220116140213
Royalblueup131131114134134134128128110
down525269494949464645
Brown4up874877625877874698856862356
down532529781529532708511517439
Plum1  up272724283030262720
down292932282626272622

MD: moderate drought, SD: severe drought, RW: rewatering


Fig. 2.

Heatmap of a royalblue module under drought stress in maize. The x-axis represents microarray samples grouped into tolerant and susceptible types and the y-axis represents genes found in the module. Below the heatmap, the corresponding module eigengene expression values are shown.


qRT-PCR을 이용한 발현양상 분석

한발 내성타입과 감수성타입 간 발현양상 비교전에 한발 처리가 제대로 수행되었는지 확인하기 위해 ZmDREB1.7의 발현량을 확인하였다. 식물에서 DREB 유전자는 식물에서 수분 스트레스에 관련된 유전자 발현을 조절한다. 옥수수에서 DREB 유전자인 ZmDREB1A와 ZmDREB2A는 수분 스트레스를 받았을 때 발현하게 된다(Liu et al. 2013). 한발 처리군의 뿌리 부분에서 Ki3와 Ki11은 한발 처리시간이 증가함에 따라 발현량이 증가하였으며, 물 재공급 처리군에서 Ki11은 시간이 증가함에 따라 발현량이 증가한 것을 확인할 수 있었다. 이는 한발 처리된 두 품종 모두 한발 스트레스가 제대로 가해졌음을 확인하였다(Fig. 3a). 이후 30개 유전자들의 발현 양상을 분석한 결과, BU050895, BF728598, CK827168, CO524848, AF457983, CF037152, AJ606944, BG836522에서 공통적으로 한발 처리군의 뿌리 부분에서 발현 패턴은 처리시간이 증가할수록 Ki3와 Ki11의 발현 패턴이 상반되는 것을 확인할 수 있었다(Fig. 3b). BF728598과 CO524848의 경우 물 재공급 처리군에서 Ki3와 Ki11에서 다른 패턴을 보이는 것을 확인하였다. 분석 결과를 통해 Ki3와 Ki11의 발현 패턴이 차이를 보이는 유전자 8개(BU050895, BF728598, CK827168, CO524848, AF457983, CF037152, AJ606944, BG836522)를 확인하였고, 그 차이는 주로 뿌리 부분에서 이루어지고 있음을 유추할 수 있었다. 따라서 유전자 30개 중에서 Ki3와 Ki11의 발현 패턴이 차이가 나는 유전자 BU050895, BF728598, CK827168, CO524848, AF457983, CF037152, AJ606944, BG836522를 한발스트레스 관련 발현마커 후보로 선발 하였다.

Fig. 3.

Expression profiles of Ki3 (susceptible type) and Ki11 (tolerant type) with selected genes under drought stress in maize seedling. A) Internal control (ZmDREB1.7), B) candidate eight genes from selected 30 genes. Data represent the means ±SD of three biological replicates. Asterisks (*, ** and ***) represent significant differences (p<0.05, p<0.01 and p<0.001, respectively) according to the t-test.


후보 유전자의 발현양상 확인

선발된 8개의 유전자를 데이터베이스를 활용하여 기능을 유추하였다(Table 4). BU050895는 ferredoxin-6으로 엽록체에서 전자전달을 수행하며, 몇몇 연구에서 한발스트레스에서 광합성계의 전자 전달에 영향을 주는 것으로 알려졌다(Lehtimäki et al. 2010, Mulo 2011). CK827168은 homeobox-leucine zipper HOX24-like로써 세포핵에서 전사인자의 활성과 시퀀스 특이적인 DNA 결합 등을 수행한다. 옥수수에서는 ABA-responsive element, low temperature responsive element 그리고 dehydration responsive element로 알려져 있으며, 한발 스트레스에서 발현이 감소한다(Zhao et al. 2011). AF457983은 종자 성숙에 역할을 하며, AJ606944는 cytokinin oxidase 3으로 세포 외 공간에서 탈수소 효소의 활성 등을 수행하며, 옥수수의 뿌리에서 다른 조직에 비해 발현이 높고, 사니토키닌의 발현을 조절하고 다양한 비생물학적 스트레스에 반응을 하는 것으로 보고되었다(Brugière et al. 2003). BG836522는 pathogenesis-related 1-like로 생물학적 자극에 반응하는 방어 기작을 수행한다. 유묘기 한발에 강한 계통 4개(CML103, CML228, CML322, Ki11), 유묘기 한발에 약한 계통 2개(Ki3, M37w), 성숙기 한발에 약한 계통 2개(CML277, NC358), reference 품종인 B73에 최종 선발된 유전자들을 적용하여 추가적인 발현 패턴 비교분석을 수행하였다.

Description of candidate 8 genes from selected 30 genes.

Gene ID (on B73 RefGen_v3)NCBI Accession no. Description GO Terms
GRMZM2G063126BU050895ferredoxin-6, chloroplastic-likechloroplast;electron carrier activity;metal ion binding;2 iron,
-BF728598   -   -
GRMZM2G407223CK827168TPA: hypothetical protein ZEAMMB73_259306   -
GRMZM2G117164CO524848homeobox-leucine zipper HOX24-likechloroplast;electron carrier activity;metal ion binding;2 iron, sequence-specific DNA binding;transcription, DNA-templated;regulation of transcription, DNA-templated
GRMZM2G448627AF457983seed maturationembryo development
-CF037152   -   -
GRMZM2G167220AJ606944cytokinin oxidase 3extracellular space;cytoplasmic membrane-bounded vesicle; UDP-N-acetylmuramate dehydrogenase activity; cytokinin dehydrogenase activity;flavin adenine dinucleotide binding;
GRMZM2G112538BG836522pathogenesis-related 1-likedefense response;response to biotic stimulus

첫 번째로 한발 처리군에서 한발 감수성 타입인 Ki3와 한발 내성 타입인 Ki11간 발현 패턴을 비교한 결과, BF728598, CO524848, CF037152에서 공통적으로 증가하는 패턴의 발현 양상을 보였다(Fig. 4). Ki3는 CK827168, AF457983, BG836522에서 6시간 한발 스트레스 처리군에서 증가하였으나, 이후 12시간 한발 스트레스 처리군에서 감소하는 패턴의 발현 양상이 나타났다. 반면, Ki11은 12시간 한발 스트레스 처리군에서 증가하는 패턴의 발현 양상을 확인함으로써 패턴의 차이가 있는 것을 확인하였다. 유묘기 한발에 강한 품종은 Ki11과 유사하게 BF728598, CO524848, CF037152, CK827168, AF457983, BG836522에서 시간이 지날수록 증가하는 패턴의 발현 양상이 나타났고, 유묘기 한발에 약한 품종에서는 Ki3와 유사한 발현 양상을 보였다. B73은 Ki3와 Ki11의 두 발현 패턴이 모두 나타났으며, 성숙기 한발에 약한 품종은 Ki3와 Ki11과는 다르게 BF728598이 감소하거나 소폭 증가하는 것을 확인하였다.

Fig. 4.

Expression profiles of four inbred lines with selected candidate eight genes under drought stress in maize root. Ki11 (tolerant type, upper left), CML322 (tolerant type, upper right), Ki3 (susceptible type, lower left), and M37w (susceptible type, lower right) were represented, respectively. Data represent the means ±SD of three biological replicates. Asterisks (*, ** and ***) represent significant differences (p<0.05, p<0.01 and p<0.001, respectively) according to the t-test.


두 번째로 물 재공급 처리군에서 한발 감수성 타입인 Ki3와 한발 내성 타입인 Ki11간 발현 패턴을 비교한 결과, 공통적으로 CK827168, AF457983, CF037152, BG836522는 증가하는 패턴을, BU050895, CO524848, AJ606944는 감소하는 패턴을 확인하였다(Fig. 5). BF728598에서 Ki3는 증가 이후 유지하는 패턴의 발현 양상을 보였고 Ki11은 시간이 지날수록 감소하는 패턴의 발현 양상을 확인하였다. 유묘기 한발에 강한 품종에서는 Ki11과 유사하지만 BF728598에서 소폭 증가하거나 감소하는 패턴을 보였고, 유묘기 한발에 약한 품종은 Ki3와 유사한 발현 양상을 보였다. B73과 성숙기 한발에 약한 품종에서는 Ki3와 Ki11과는 다른 새로운 패턴의 발현 양상을 확인하였다.

결론적으로, 최종 선발된 8개의 유전자를 이용하여 9개 품종에 적용한 결과 한발 감수성 타입과 내성 타입의 유전자 발현 양상이 다른 것을 확인하였다. 추후에 더 많은 품종들에 적용한다면, 이 8개의 유전자를 이용하여 한발에 감수성인지 내성인지를 유묘기에 판단할 수 있는 한발 관련 발현 마커로써 개발이 가능할 것이다.

Fig. 5.

Expression profiles of four inbred lines with selected candidate eight genes under rewatering stress in maize root. Ki11 (tolerant type, upper left), CML322 (tolerant type, upper right), Ki3 (susceptible type, lower left), and M37w (susceptible type, lower right) were represented, respectively. Data represent the means ±SD of three biological replicates. Asterisks (*, ** and ***) represent significant differences (p<0.05, p<0.01 and p<0.001, respectively) according to the t-test.


적요

본 연구는 유묘기 옥수수에서 한발 관련 유전자를 이용하여 발현 마커의 탐색 및 개발을 위해 NCBI GEO dataset에서 옥수수 microarray CEL 파일을 이용하여 한발에 관련된 유전자들의 발현을 분석하였다. 유전자들의 co-expression 분석을 통해 6개의 모듈을 구성하였으며, 이 중 한발 내성 타입과 감수성 타입에서 발현 차이를 보이는 royal blue모듈 내의 유전자들 중 상위 30개의 유전자들을 선발하였다. Ki11(한발 내성)과 Ki3(한발 감수성)의 유묘기에 한발 처리 후 qRT-PCR을 이용하여 선발된 유전자들의 발현 패턴을 비교 분석을 통하여 8개의 유전자(BU050895, BF728598, CK827168, CO524848, AF457983, CF037152, AJ606944, BG836522)가 유묘기 옥수수의 뿌리에서 발현 차이를 보이는 것을 확인하였다. 최종 선발된 유전자들을 9개의 품종에 한발 처리 후 뿌리 부분을 타겟으로 qRT-PCR을 진행한 결과, 한발 처리군에서 CK827168과 BG836522가 한발 내성 타입과 감수성 타입 간 발현 양상의 차이를 확인할 수 있었다. 이 결과를 통해 스트레스에 관련된 유전자를 이용하여 발현 마커의 개발 및 적용 가능성을 보여준다.

사사

본 연구는 농림축산식품부, 해양수산부, 농촌진흥청, 산림청 Golden Seed 프로젝트 사업(식량종자사업단, 과제번호: 213009-05-2-SB710)과 농촌진흥청 공동연구사업(과제번호: PJ012649022018)의 지원에 의해 이루어진 것임.

References
  1. Arumuganathan K, and Earle ED. Nuclear DNA content of some important plant species. Plant Mol Biol Rep 1991;9:208-218.
    CrossRef
  2. Banziger M, Edmeades GO, Beck DL, and Bellon MR. Breeding for drought and nitrogen stress tolerance in maize:from theory to practice. Mexico CIMMYT 2000.
  3. Betrán FJ, Beck D, Bänziger M, and Edmeades GO. Genetic analysis of inbred and hybrid grain yield under stress and nonstress environments in tropical maize. Crop Sci 2003;43:807-817.
    CrossRef
  4. Blanc G, and Wolfe KH. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell 2004;16:1667-1678.
    Pubmed KoreaMed CrossRef
  5. Brugière N, Jiao S, Hantke S, Zinselmeier C, Roessler JA, Niu X, Jones RJ, and Habben JE. Cytokinin Oxidase Gene Expression in Maize Is Localized to the Vasculature, and Is Induced by Cytokinins, Abscisic Acid, and Abiotic Stress. Plant Physiol 2003;132:1228-1240.
    CrossRef
  6. Buckler ES, Holland JB, Bradbury PJ, Acharya CB, Brown PJ, Browne C, Ersoz E, Flint-Garcia S, Garcia A, and Glaubitz JC et al. The genetic architecture of maize flowering time. Science 2009;325:714-718.
    Pubmed CrossRef
  7. Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L, Murigneux A, and Charcosset A. Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics 2004;168:2169-2185.
    Pubmed KoreaMed CrossRef
  8. De La Cruz M, Romao RL, Escudero A, and Maestre FT. Where do seedlings go? A spatio-temporal analysis of seedling mortality in a semi-arid gypsophyte. Ecography 2008;31:720-730.
    CrossRef
  9. Fan CY, Hao ZQ, Yan JH, and Li GL. Genome-wide identification and functional analysis of lincRNAs acting as miRNA targets or decoys in maize. BMC Genomics 2015;16.
    CrossRef
  10. Gale MD, and Devos KM. Comparative genetics in the grasses. Proc Natl Acad Sci U S A 1998;95:1971-1974.
    Pubmed CrossRef
  11. Gulli M, Salvatori E, Fusaro L, Pellacani C, Manes F, and Marmiroli N. Comparison of drought stress response and gene expression between a GM maize variety and a near-isogenic non-GM variety. PLoS One 2015;10.
    CrossRef
  12. Hao ZF, Li XH, Liu XL, Xie CX, Li MS, Zhang DG, and Zhang SH. Meta-analysis of constitutive and adaptive QTL for drought tolerance in maize. Euphytica 2010;174:165-177.
    CrossRef
  13. Krill AM, Kirst M, Kochian LV, Buckler ES, and Hoekenga OA. Association and linkage analysis of aluminum tolerance genes in maize. PLoS One 2010;5.
    CrossRef
  14. Langfelder P, and Horvath S. WGCNA:an R package for weighted correlation network analysis. BMC Bioinformatics 2008;9.
    CrossRef
  15. Leach KA, Hejlek LG, Hearne LB, Nguyen HT, Sharp RE, and Davis GL. Primary root elongation rate and abscisic acid levels of maize in response to water stress. Crop Sci 2011;51:157-172.
    CrossRef
  16. Lehtimäki N, Lintala M, Allahverdiyeva Y, Aro E-M, and Mulo P. Drought stress-induced upregulation of components involved in ferredoxin-dependent cyclic electron transfer. J Plant Physiol 2010;167:1018-1022.
    Pubmed CrossRef
  17. Liu SX, Wang XL, Wang HW, Xin HB, Yang XH, Yan JB, Li JS, Tran LSP, Shinozaki K, Yamaguchi-Shinozaki K, and Qin F. Genome-wide analysis of ZmDREB genes and their association with natural variation in drought tolerance at seedling stage ofZea maysL. PLoS Genet 2013;9.
    CrossRef
  18. Manoli A, Sturaro A, Trevisan S, Quaggiotti S, and Nonis A. Evaluation of candidate reference genes for qPCR in maize. J Plant Physiol 2012;169:807-815.
    Pubmed CrossRef
  19. Meeks M, Murray SC, Hague S, and Hays D. Measuring maize seedling drought response in search of tolerant germplasm. Agronomy-Basel 2013;3:135-147.
    CrossRef
  20. Mulo P. Chloroplast-targeted ferredoxin-NADP+oxidoreductase (FNR):Structure, function and location. Biochim Biophys Acta 2011;1807:927-934.
    Pubmed CrossRef
  21. Myles S, Peiffer J, Brown PJ, Ersoz ES, Zhang Z, Costich DE, and Buckler ES. Association mapping:critical considerations shift from genotyping to experimental design. Plant Cell 2009;21:2194-2202.
    Pubmed KoreaMed CrossRef
  22. Paterson AH, Bowers JE, and Chapman BA. Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics. Proc Natl Acad Sci U S A 2004;101:9903-9908.
    Pubmed KoreaMed CrossRef
  23. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, and Poliakov A et al. The Sorghum bicolor genome and the diversification of grasses. Nature 2009;457:551-556.
    Pubmed CrossRef
  24. SanMiguel P, Gaut BS, Tikhonov A, Nakajima Y, and Bennetzen JL. The paleontology of intergene retrotransposons of maize. Nat Genet 1998;20:43-45.
    Pubmed CrossRef
  25. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, and Graves TA et al. The B73 maize genome:complexity, diversity, and dynamics. Science 2009;326:1112-1115.
    CrossRef
  26. Swigonova Z, Lai JS, Ma JX, Ramakrishna W, Llaca V, Bennetzen JL, and Messing J. Close split of sorghum and maize genome progenitors. Genome Res 2004;14:1916-1923.
    Pubmed KoreaMed CrossRef
  27. Thirunavukkarasu N, Hossain F, Mohan S, Shiriga K, Mittal S, Sharma R, Singh RK, and Gupta HS. Genome-wide expression of transcriptomes and their co-expression pattern in subtropical maize (Zea maysL.) under waterlogging stress. PLoS One 2013;8.
    CrossRef
  28. Thornsberry JM, Goodman MM, Doebley J, Kresovich S, Nielsen D, and Buckler ESt. Dwarf8 polymorphisms associate with variation in flowering time. Nat Genet 2001;28:286-289.
    Pubmed CrossRef
  29. Tollenaar M, and Lee EA. Yield potential, yield stability and stress tolerance in maize. Field Crops Res 2002;75:161-169.
    CrossRef
  30. Unterseer S, Bauer E, Haberer G, Seidel M, Knaak C, Ouzunova M, Meitinger T, Strom TM, Fries R, and Pausch H et al. A powerful tool for genome analysis in maize:development and evaluation of the high density 600 k SNP genotyping array. BMC Genomics 2014;15.
    CrossRef
  31. Wang L, Liu K, Mao SS, Li ZY, Lu YL, Wang JR, Liu YX, Wei YM, and Zheng YL. Large-scale screening for Aegilops tauschii tolerant genotypes to phosphorus deficiency at seedling stage. Euphytica 2015;204:571-586.
    CrossRef
  32. Yu J, Holland JB, McMullen MD, and Buckler ES. Genetic design and statistical power of nested association mapping in maize. Genetics 2008;178:539-551.
    Pubmed KoreaMed CrossRef
  33. Zhang W, Han ZX, Guo QL, Liu Y, Zheng YX, Wu FL, and Jin WB. Identification of maize long non-coding RNAs responsive to drought stress. PLoS One 2014;9.
    CrossRef
  34. Zhao Y, Zhou Y, Jiang H, Li X, Gan D, Peng X, Zhu S, and Cheng B. Systematic Analysis of Sequences and Expression Patterns of Drought-Responsive Members of the HD-Zip Gene Family in Maize. PLoS One 2011;6:e28488.
    Pubmed KoreaMed CrossRef


December 2018, 50 (4)
Full Text(PDF) Free

Social Network Service
Services

Cited By Articles
  • CrossRef (0)

Funding Information

Warning: mysqli::close() [function.mysqli-close]: Couldn't fetch mysqli in /home/kjbs/www/inc/bottom.inc on line 49