콩명나방(Maruca vitrata)(나비목: 포충나방과)은 주로 콩 과(Fabaceae) 작물의 꽃과 꼬투리, 줄기에 피해를 주는 해충인 데, 동부(Vigna unguiculata subsp. unguiculata)와 녹두(V. radiata), 팥(V. angularis) 등의 Vigna속 작물과 콩(Glycine max)이 주요 피해작물로 알려져 있다(Sharma, 1998;Jung et al., 2009;Margam et al., 2011;Grigolli et al., 2015;Srinivasan et al., 2021). 작물 피해는 아프리카와 동남아시아, 중남미아메 리카를 포함하는 주로 열대와 아열대 지역에서 보고되었으나 (Jackai and Daoust, 1986;Sharma, 1998;Grigolli et al., 2015), 온대지역인 우리나라에서도 팥의 생식생장기 중 같은 시기에 발생한 팥나방(Matsumuraeses phaseoli)(잎말이나방과)과 콩 줄기명나방(Ostrinia scapulalis)(포충나방과) 보다 우점하면 서(조사식물체의 약 83%에 분포) 팥에 심한 피해를 준 사례가 보고되었다(Jung et al., 2009;2012).
국내에서의 연중 발생 시기와 세대는 현재까지 뚜렷하게 밝혀 져 있지 않다. 겨울철 온도가 영하 이하로 떨어지는 국내 수원 지 방에서 콩명나방의 야외 월동이 불가능했었던 결과를 근거로 하 여 콩명나방은 적어도 국내 중부지역에서는 월동할 수 없을 것으 로 추정되었다(Jung et al., 2016). 대신 남쪽 지역이나 해외로부 터 비래하는 집단에 의해 매년 새로운 개체군이 형성된 것으로 보고있다. 콩명나방 성충의 이동성에 대해서는, 서아프리카에서 연중 성충 발생양상과 초위성체 유전자 분석으로 기주식물을 찾 아 장거리 이동을 하는 것으로 보고되었으며(Bottenberg et al., 1997;Agunbiade et al., 2012), 일본 서남쪽 해상에서 6월 중 콩 명나방 성충이 기후관측 선박에 설치된 유인등에 포획된 것으로 중국으로부터의 비행 발원지를 추정한 결과가 있으며(Kawazu et al., 2008), 중국 동쪽 Jiangpu지역에서 공중에 설치된 트랩에 성충이 포획된 보고가 있어, 콩명나방 성충이 장거리 이동 능력 이 있을 것으로 짐작되고 있다(Riley et al., 1995). 국내에서는 약 30년동안(1971-2000년) 채집된 성충 표본 기록(Bae, 2001) 분 석을 통해 콩명나방은 연중 4월과 11월 사이에 발생한 것으로 확 인되어(Jung et al., 2016), 첫 세대 발생 이후 다음 세대들이 진전 할 것으로 짐작된다. 따라서 늦여름과 가을에 걸쳐 생식생장기 를 겪는 팥 재배기간 중에는 적어도 그 해 1세대 이상을 거친 개 체군이 팥 포장으로 이입되는 것으로 가정된다.
곤충은 변온동물로서 발육과 생존, 이동 등의 생명현상들이 온도에 크게 영향을 받는다(Hoffmann, 1985;Ratte, 1985). 온 도를 매개변수로 하여 온도와 곤충 생명현상과의 관계를 설명 할 수 있는 수리함수들이 개발되어 왔고, 이를 통해 곤충 개체 군의 밀도 변동 과정을 이해하고 해충 방제에 응용하여 왔다 (Campbell et al., 1974;Logan et al., 1976;Sharpe and DeMichele, 1977;Schoolfield et al., 1981;Lactin et al., 1995;Briere et al., 1999;Kim and Lee, 2003;Ikemoto 2005;Kim et al., 2017;Ahn et al., 2019a, b;Jung et al., 2021). 콩명나방의 발육 에 미치는 온도 영향의 결과는 대만(Huang and Peng, 2002), 나 이지리아(Adati et al., 2004), 일본(Chi et al., 2006) 개체군에 대해서 일부 알려져 있다. 본 연구에서는 국내에서 채집한 콩명 나방 개체군을 대상으로 인공사료를 이용하여 미성숙태 발육 과 성충의 생식 특성에 대한 온도 영향을 조사하였고 이들 자료 를 분석하여 발육 및 산란 시기 예측가능 모델 개발에 필요한 다양한 수리모형을 추정하였다. 조사 및 분석 결과는 향후 콩명 나방의 개체군 밀도 변동 시기를 추정하는데 응용하여 콩명나 방 피해에 대비하고자 하였다.
재료 및 방법
실험곤충
콩명나방은 2004년 경기도 수원시에 소재한 국립식량과학 원 팥포장에서 팥 생식기관을 가해하는 유충 상태로 채집하여, Jung et al. (2016)의 보고와 동일한 방법으로 유충을 인공사료 로 사육하면서 실험실 집단으로 유지하였다.
온도별 발육 및 생식 실험
실험실에서 유지한 콩명나방 집단에서 알, 유충, 성충을 각 각 선택하여 13~34°C 온도범위에서 3°C 간격으로 8개 항온조 건에서 발육실험을 진행하였다. 16°C 이하의 온도는 2연실 저 온배양기(Dasol, Hwaseong, Korea)를, 19°C 이상의 온도는 4 연실 배양기(Dasol, Hwaseong, Korea)를 이용하였다. 일일 광 주기를 15:9(명:암) 시간으로, 습도를 60±5%로 설정하였다.
알 발육실험은 실내 집단의 성충 교미상자에서 낱개로 산란 된 알이 부착된 키친타월을 산란 당일 수거하여 미리 온도를 설 정해 놓은 배양기에 옮겨 실험을 실시하였다. 알이 30개 붙은 키친타월 조각을 플라스틱 페트리접시(50 mm in diameter, 10 mm in height)에 넣고 부화할 때까지 매일 관찰하면서 부화한 유충 수를 세고 제거하였다. 부화하지 않은 알은 온도별로 1 주~2주 더 두어 부화 여부를 확인하였다. 각 온도별로 페트리 접시 4개를 반복으로 처리하였다. 부화율은 페트리접시 수를 반복으로, 발육기간은 각 온도에서 부화한 모든 알 수를 반복으 로 하여 평균값을 산출하였다.
유충과 번데기 기간 동안의 발육 특성은 실내 집단에서 부화 후 12시간이 경과하지 않은 유충을 채집하여 페트리접시(50 mm in diameter, 10 mm in height)에 한 마리씩 인공사료와 함께 넣어 사육하면서 조사하였다. 성충으로 우화할 때까지 매일 조사하면 서 탈피와 사망 여부를 기록하였는데, 유충 영기와 번데기 시기는 허물을 벗은 날을 기준으로 구분하였다. 각 온도별로 100마리씩 을 처리하였고, 발육태별로 발육기간과 생존율을 산출하였다.
성충 생식 특성은 실내 집단으로부터 성별이 구분된 번데기 를 개체별로 놓고, 우화 당일 암수 1쌍을 플라스틱 상자[72 (width) × 72(depth) × 100(height) mm](SPL Life Sciences, Pocheon, Korea)에 넣어 사육하면서 상자 안쪽 면에 산란하도 록 하였다. 먹이로 증류수와 10% 설탕물을 유리바이알에 넣어 공급하였다. 암컷 성충이 사망할 때까지 매일 산란과 사망 여부 를 관찰하면서 산란한 교미쌍은 새 상자에 옮겼다. 알이 부착된 상자는 25°C에 옮겨 부화한 유충 수를 조사하였다. 온도별로 30쌍을 처리하였다. 산란전 기간과 산란기간, 산란후 기간 및 성충 수명, 암컷당 총산란수를 산출하였다.
자료분석
콩명나방 각 발육단계별 발육기간, 생존율, 성충의 수명, 산 란수는 SAS에 있는 일반선형모형(PROC GLM)을 이용하여 통계적 유의성을 분석하였다(SAS Institute, 2004). 실험에 이 용된 온도별 각 발육단계의 발육기간, 성충의 수명, 산란수가 통계적 차이를 보이면 Tukey의 평균간 다중비교를 이용하여 처리온도간의 평균값을 비교하였다. 선형 및 비선형함수의 매 개변수들은 TableCurve2D (Jandel Scientific, 1994) 혹은 R statistics (2015)를 이용하여 추정하였다.
발육모형 매개변수 추정
온도와 발육단계별 발육율(발육기간의 역수)과의 관계는 선 형과 비선형함수를 이용하여 분석하였다. 선형함수는 Campbell et al. (1974)이 제시한 방법을 이용하여 각 발육단계별 선형함 수(y = aT+b, y = 발육율, T = 실험온도)에서 얻어진 결과를 바 탕으로 발육영점온도()와 유효적산온일도()을 계산하였다. 비선형함수는 Ratkowsky and Reddy (2017)와 Shi et al. (2017) 의 평가를 기초로 Lobry-Rosso-Flandrois (LRF)함수(식 1)과 Sharpe-Schoolfield-Ikemoto (SSI)함수(식 2)를 이용하였다.
이 식에서 r(T)는 실험온도 T(절대온도)에서 발육율, μopt 는 T opt에서 발육율, T opt는 발육최적온도, T max는 발육이 진 행되는 최고온도, T min는 발육이 진행되는 최저온도이다.
이 식에서 r(T)는 절대온도 oK에서 발육율, R은 기체상수 (1.987 cal degree-1 mole-1), ρϕ는 효소 활성의 저해가 없다고 가 정한 기준 온도, Tϕ는 공통고유최적온도로서 활성화상태에서 속도조절효소의 활성이 최고조에 이른 온도, ΔHA는 발육에 관여하는 속도조절효소가 촉매하는 반응의 활성화 엔탈피, ΔHL은 속도조절효소의 활성이 저온영역에서 저해되는 것과 관련된 엔탈피의 변화(HL), TL은 속도조절효소가 1/2 불활성 화되어 1/2 활성화 되어 있는 저온 영역의 절대온도(TL), ΔHH 는 속도조절효소의 활성이 고온영역에서 저해되는 것과 관련 된 엔탈피의 변화(HH), TH는 속도조절효소가 1/2 불활성화되 어 1/2 활성화 되어 있는 고온 영역의 절대온도(TH)이다(Sharpe and DeMichele, 1977;Schoolfield et al., 1981;Ikemoto, 2005). Shi et al. (2017)가 개발한 R 프로그램을 이용하여 SSI모형의 매개변수를 계산하였다.
발육단계별 다음 발육단계 출현 분포에 대한 누적비율과 정 규화한 발육시간(발육단계별 발육기간/평균발육기간)의 관계 (발육완료분포모형)는 3-매개변수 Weibull함수(식 3)를 이용 하였다(Weibull, 1951;Wagner et al., 1984).
F(t)는 정규화된 시간(각 개체의 발육기간/평균발육기간), γ, α, β는 Weibull함수의 매개변수들이다. 정온조건에서 콩명 나방 성충이 출현하는 빈도를 비선형함수와 발육완료분포모형 을 이용하여 추정하였다(식 4).
F(x, T)는 시간 x와 온도 T에서 콩명나방 성충이 출현하 는 빈도로서 r(T)는 온도와 발육율간의 관계를 설명하는 비선 형함수 γ, α, β는 발육완료분포모형에서 얻은 매개변수들이 다. F(x, T)의 결과는 xt+1과 xt 결과의 차이를 나타낸다.
산란모형 매개변수 추정
성충산란모형은 성충노화율함수, 연령별 생존율함수, 연령 별 누적산란율함수, 온도별 총산란수함수로 구성되어 있다.
성충노화율함수
성충노화율(1/성충수명)과 온도와의 관계는 비선형함수(식 5)를 이용하였다.
r(T)는 온도(T)에서 성충의 노화율, a, b는 비선형함수의 매 개변수이다.
연령별 생존율함수
콩명나방의 연령별 생존율은 2-매개변수 Weibull함수(식 6) 를 이용하였다(Pinder III et al., 1978).
S(t)는 동일 연령집단의 개체들에서 시간 t (정규화한 시간 = 각 개체의 성충연령/평균연령)에 살아있는 개체수의 비율을 의미하며, α, β는 함수식의 매개변수들이다.
연령별 누적산란율함수
성충의 연령변화에 따른 누적산란율은 2-매개변수 Weibull 함수(식 7)를 이용하였다.
F(t)는 정규화된 시간(각 성충암컷의 산란시기/평균성충수 명), c, d는 Weibull함수의 매개변수들이다.
온도별 총산란수함수
온도에 따른 총산란수 함수는 gaussian 함수(식 8)를 이용하 였다.
f(t)는 온도 T에서 암컷 성충의 총산란수, f, g, h는 함수와 관련된 매개변수이다.
동일 연령집단의 암컷으로부터 생산될 수 있는 암컷당 일일 산란수(시간 ti와 ti+1의 차이)는 다음과 같이 추정하였다(식 9).
f(t), F(t), s(t)는 위 식에서 정의하였다.
결 과
발육단계별 온도발육
콩명나방은 알 발육단계에서 온도조건 13~34°C에서 부화 율은 23.3~100.0%였고 발육기간은 2.3~15.3일로 온도처리에 따라 유의한 차이를 나타내었다 (Table 1). 유충에서 성충으로 발육은 13°C, 34°C를 제외한 항온조건(16~31°C)에서 성충으 로 발육이 가능하였다(Table 2). 유충부터 성충이 출현하는데 필요한 발육기간은 16.9일(31°C)~78.4일(16°C)로 실험 온도 가 올라갈수록 짧아졌다 (Table 2). 유충과 번데기시기의 생존 율은 41%~94%, 18%~92%로 고온조건 (31°C)에서 가장 낮은 생존율을 보였다. 콩명나방의 발육단계별 발육율과 온도와의 관계는 선형과 비선형함수를 이용하여 나타내었다(Fig. 1). 알, 유충, 번데기, 부화부터 성충출현까지 추정된 발육영점온도는 9.8°C, 12.5°C, 13.3°C, 12.8°C였다(Table 3). 알, 유충, 번데기, 유충부터 성충출현까지 유효적산온일도는 49.14, 166.09, 95.16, 280.834DD였다.
SSI함수를 이용하여 분석한 콩명나방 알 발육기간과 부화부 터 성충출현까지 TL은 14.8°C, 14.3°C, TΦ는 29.7°C, 27.5°C, TH는 34.7°C, 31.9°C였다. 1령부터 5령유충까지 최적공통온 도는 25.5°C, 26.2°C, 25.9°C, 35.7°C, 27.1°C였다. 콩명나방의 알 발육단계 및 유충에서 성충출현까지 발육적정온도 즉 TL과 TH의 차이는 19.9°C, 17.7°C로 알의 발육적정온도범위가 넓었 다(Table 4).
3-매개변수 Weibull함수는 콩명나방 발육단계별 발육완료 분포와 정규화된 시간과의 관계를 잘 설명해주었다(Fig. 2) (알: F2, 11 = 1753.59, P < 0.0001, r2 = 0.99, 전체유충기간: F2, 42 = 437.62, P < 0.0001, r2 = 0.95, 번데기: F2, 29 = 2 81.70, P < 0.0001, r2 = 0.95, 부화부터 성충출현: F2, 43 = 685.69, P < 0.0001, r2 = 0.97). 12~36°C의 항온조건에서 콩명나방 알 부화 부터 성충출현기간의 온도발육함수과 발육완료분포함수의 매 개변수들을 이용하여 콩명나방 성충출현빈도를 추정한 결과 30°C (16일)에서 성충 출현빈도가 가장 높았다(Fig. 3).
온도별 성충의 수명, 산란수 및 산란모형
콩명나방 성충은 13~34°C 범위에서 생존과 산란이 가능하 였다(Table 5). 암, 수컷의 평균수명은 저온에서 온도가 올라갈 수록 줄어드는 경향이었고 산란수는 25°C (511개)에서 가장 많았다. 콩명나방 성충의 산란전기간, 산란기간, 산란후기간, 수명, 산란수는 온도의 영향을 많이 받았다. Fig. 4는 온도별 콩 명나방 성충의 일별 생존율 및 산란수의 변화를 보여주고 있다. 생존율은 13°C에서 34일 이후, 16°C에서 26일 이후, 19°C에서 18일 이후, 22°C에서 7일 이후, 25°C에서 9일 이후, 28°C에서 6일 이후, 31°C에서 12일 이후, 34°C에서 6일 이후에 급속히 감소하였다.
콩명나방 산란모형과 관련된 여러 함수들의 매개변수는 Table 6에 제시하였다. 산란이 가능한 추정온도범위는 11~37°C 이며 최대산란수를 생산하는 추정온도는 24.5°C이다(Fig. 5). 암컷 성충의 연령과 항온 조건을 고려하여 추정한 일별 산란수 는 Fig. 6에 나타내었다.
고 찰
콩명나방은 주로 열대지방에 서식하며 개체군 변동 및 작물 피해에 대한 연구들이 이들 지역을 중심으로 많이 이루어져 왔다 (Agunbiade et al., 2014, Sreekanth et al., 2015, Taggar et al., 2019, Shejulpatil et al., 2020). 본 연구는 수원지역에서 채집하 여 인공사료를 이용하여 누대사육하고 있던 콩명나방을 이용하 여 발육단계별 발육기간, 생존율, 성충수명과 산란을 13~34°C 범 위 8개 항온조건에서 조사하였다. 콩명나방 알 발육단계는 13~ 34°C 범위에서 부화할 수 있었고 유충은 16~31°C범위에서 성충 으로 발육할 수 있었다(Table 1 and 2). Jung et al. (2007)은 13°C 에서 알이 부화하지 못하였다고 기술하였는데 본 연구에서는 23.3% 부화하여 실험에 이용된 개체에 차이가 있는 것으로 사료 된다. 알, 유충, 번데기의 발육기간과 성충의 수명은 온도가 증가 할수록 짧아지는 경향을 보였다. Jung et al. (2007)이 25°C에서 조사한 발육단계별 발육기간은 본 연구와 유사하였다.
Adati et al. (2004)는 나이지리아 Ibadan (12°09´N, 8°40´E) 에 위치한 열대농업국제연구소(International Institute of Tropical Agriculture)에서 얻은 콩명나방을 이용하여 인공사료를 이용 하여 실험을 진행하였다. 본 연구에서 얻은 알의 부화율(23.3~ 100%)이 Adati et al. (2004)의 결과보다 높았고(7.0~97.2%), 발육기간은 유사하였다. Chi et al. (2006)은 10°C 에서 알이 부 화하지 못하였다고 보고하였다. 유충생존율은 19°C를 제외한 온도범위에서 본 연구의 결과가 Adati et al. (2004)의 결과보다 높았고 발육기간은 28°C이상 조건에서 본 연구의 결과가 짧았 다(Table 7). Chi et al. (2006)은 20, 25, 30°C(21.3, 12.0, 10.0 일)에서 본 연구보다 짧은 발육기간을 보여주었다. 본 연구와 Adati et al. (2004)은 34°C조건에서 콩명나방 유충이 번데기가 되지 못하였다. 번데기 생존율은 Adati et al. (2004)의 결과가 모 든 온도범위에서 본 연구의 결과보다 높았고 발육기간은 31°C 조건에서만 본 연구의 결과가 짧았다. Chi et al. (2006)은 20, 25°C(13.7, 6.7일)에서 본 연구보다 짧은 발육기간을 보여주었 다. Adati et al. (2004)는 발육영점온도와 유효적산온일도를 계 산하는 방법에 대해 기술하지 않고 제시한 알, 유충, 번데기의 발육영점온도는 10.5°C, 9.8°C, 10.9°C였으며, Chi et al. (2006) 의 경우는 11.03°C, 10.50°C, 11.05°C로 보고하였다. 그들의 결과를 본 연구의 결과 9.8°C, 12.5°C, 13.3°C와 비교할 때 알 발육단계의 발육영점온도는 그들의 연구 결과가 높았고 유충 과 번데기의 발육영점온도는 본 연구 결과가 상대적으로 높게 나타났다. Huang and Peng (2002)의 경우는 7.1°C, 8.2°C, 9.1°C로 본 연구의 결과보다 모든 발육단계에서 낮은 발육영점 온도를 나타내었다. SSI함수를 이용하여 얻은 본 연구의 알, 유 충 번데기 발육단계의 TL은 14.8°C, 14.0°C, 15.3°C, TΦ는 29.7°C, 27.4°C, 27.6°C, TH는 34.7°C, 31.9°C, 32.5°C였고 Adati et al. (2004)의 결과를 재분석한 결과 알, 유충 번데기 발 육단계의 TL은 11.2°C, 16.8°C, 8.0°C, TΦ는 24.5°C, 28.8°C, 17.6°C, TH는 33.4°C, 35.2°C, 31.8°C였다. 콩명나방의 알 발 육단계 및 유충에서 성충출현까지 발육적정온도 즉 TL과 TH의 차이는 본 연구에서 19.9°C, 17.7°C, Adati et al. (2004)은 22.3°C, 20.7°C였다. 본 연구에 이용된 콩명나방보다 나이지리 아 계통 개체군이 넓은 온도범위에 적응된 것으로 추정되었다. Chi et al. (2005)가 조사한 성충의 산란전기간(20~27°C, 2.9~ 5.2일)은 20°C를 제외하고 본 연구결과보다 길었다. 암컷당 산 란수는 20~27°C 범위에서 Chi et al. (2005)의 결과(644~856 개)가 본 연구보다 많았고 성충수명(20.7~26.5일)도 길었다.
Jung et al. (2009)은 2005년 수원지역 팥 포장에서 나비목 해충의 발생과 피해양상 조사를 통해 콩명나방 1령 유충이 전 체 기간 중 8월 29일 조사된 것이 가장 높은 비율을 차지하는 것을 이용하여 Adati et al. (2004)이 발표한 알 발육영점온도 (10.5°C)와 유효적산온일도(51.1DD), 성충 우화 후 최대산란 일(20~27°C에서 약 7~9일)을 고려하여 (Chi et al., 2005) 콩명 나방 성충의 우화시점을 유충발견 11~13일 이전인 8월 17일 로 추정하였다. 본 연구를 통해 얻어진 알 발육영점온도 (9.8°C)와 유효적산온일도(49.1DD), 산란전기간(22~25°C에서 약 2~4일), 성충 우화 후 최대산란일(22~25°C에서 약 6~7일) 과 실제 2005년 수원지역 평균기온값을 이용하여 계산을 한 결과 콩명나방 성충의 우화시점은 유충발견 12~14일 이전인 8 월 16일로 추정되어 Jung et al. (2009)이 추정한 결과와 유사 한 결과를 얻을 수 있었다. 현재 콩과 팥에 발생하는 콩명나방 을 방제하기 위해 등록되어있는 약제는 65개 이며 (http://www. nonsaro.go.kr)(Accessed on 15 Oct. 2022) 성충 발생 최성기 3~7일 이후에 약제를 살포하는 것이 방제적기로 사료된다.
본 연구에서는 인공사료를 이용하여 콩명나방 알, 유충의 온 도발육과 생존율, 성충의 수명과 산란을 조사하였으며 이를 바 탕으로 온도발육모형, 발육완료모형, 산란모형을 작성하였다. 본 연구에서 제시한 생물적 자료들과 수학적 모델들은 콩과작 물에 피해를 주는 콩명나방의 생물계절학을 이해하고 콩명나 방을 포함한 해충들의 종합적인 개체군관리를 수행하기 위한 기본적인 자료로 이용될 수 있을 것이라 사료된다.