# Genetic Dissection of Sympatric Populations of Brown Planthopper,Nilaparvata lugens (Stål), Using DALP-PCR Molecular Markers

Authors: M. A. Latif; M. Y. Rafii; M. S. Mazid; M. E. Ali; F. Ahmed; M. Y. Omar; S. G. Tan
Journal: The Scientific World Journal
 (2012)

---

## Abstract

Direct amplified length polymorphism (DALP) combines the advantages of a high-resolution fingerprint method and also characterizing the genetic polymorphisms. This molecular method was also found to be useful in brown planthopper,Nilaparvata lugens species complex for the analysis of genetic polymorphisms. A total of 11 populations of Nilaparvata spp. were collected from 6 locations from Malaysia. Two sympatric populations of brown planthopper, N. lugens, one from rice and the other from a weed grass (Leersia hexandra), were collected from each of five locations. N. bakeri was used as an out group. Three oligonucleotide primer pairs, DALP231/DALPR′5, DALP234/DALPR′5, and DALP235/DALPR′5 were applied in this study. The unweighted pair group method with arithmetic mean (UPGMA) dendrogram based on genetic distances for the 11 populations of Nilaparvata spp. revealed that populations belonging to the same species and the same host type clustered together irrespective of their geographical localities of capture. The populations of N. lugens formed into two distinct clusters, one was insects with high esterase activities usually captured from rice and the other was with low esterase activities usually captured from L. hexandra. N. bakeri, an out group, was the most isolated group. Analyses of principal components, molecular variance, and robustness also supported greatly to the findings of cluster analysis.


---

## Body

## 1. Introduction

The brown plant hopper,Nilaparvata lugens (Stål) (Homoptera: Delphacidae), is a major pest of rice, which is widely distributed from tropical to temperate areas of Asia and Australia. The insect is a phloem-feeder and is restricted to cultivated and wild rice as host plants. It causes “hopperburn” and complete wilting and drying of rice plants [1] and also transmits the grassy stunt and ragged stunt viral diseases [2]. Large-scale rice crop damage caused by the pest was reported in the 1970s in several South and Southeast Asian countries [1]. Another population of brown plant hopper was found to infest a weed grass, Leersia hexandra, which grows abundantly in canal near paddy fields in South East Asia [3, 4]. The weed infesting population ofN. lugens fails to survive on rice plants. Conversely, rice infesting population of N. lugens does not thrive on grass [5]. Based on nymphal survival, virulence, ovipositional preference, mate choice, and hybridization experiments, Claridge et al. [5] suggested that the rice and Leersiainfesting populations of brown planthopper (BPH) represented two distinct sympatric biological species. In recent studies, the analyses of isozymes and RAPD-PCR markers indicated that BPH with high esterase activity usually captured from rice plant, and those with low esterase activity usually captured fromL. hexandra in Malaysia represent two distinct closely related sibling species [6, 7].Direct amplification of length polymorphisms (DALPs) is a technique which uses arbitrarily primed PCR (AP-PCR) to produce genomic finger prints and to enable sequencing of DNA polymorphisms in any species. Oligonucleotide pairs were designed to produce a specific multibanded pattern for each individual of a population and between populations. This strategy combines the advantages of a high-resolution fingerprint technique and also characterizing the polymorphisms [8]. Higher number of polymorphic loci could be detected and isolated for sequencing in only one step. Therefore, this method is not simply another supplementary molecular fingerprinting technique but was designed from the very beginning to obtain nucleotide sequence information on DNA fragments from any genome with no need for a genomic library.Genetic polymorphic markers, such as isozymes, RAPD and SSR, and nuclear or organelle DNA polymorphism, have been developed for a variety of studies on genetic diversity, population structures, and subdivisions [6, 7, 9–11]. The present study was undertaken to analyze genetic diversity as well as to detect genetic structures between two sympatric populations of N. lugens, one from rice and the other from weed grass, L. hexandra. We hypothesized that the molecular method newly applied in rice brown planthopper could be able to detect structures among the populations of brown planthopper species complex.


## 2. Experimental Section

### 2.1. Collection of Insect Populations

A total of 11 populations were collected from 6 locations. Two sympatric populations ofN. lugens, one from rice and other fromL. hexandra, were collected from each of five locations. The locations were Universiti Putra Malaysia (UPM), Tanjung Karang (TK), Melaka (MK), Perak (PK), and Sabah (SB). An out group, N. bakeri was also collected from Cameron Highlands (CH), Malaysia. Locations, host type, date of collection, population code are shown in Table 1. Each population consisted of twenty insects. All collected insects were frozen at −70°C for further use.Table 1
Host types, sites of collection, and coding for 11 populations ofNilaparvata spp.
SpeciesLocationsHost plantsPopulation codeN. lugensUniversiti Putra Malaysia (UPM), Selangor, MalaysiaRiceUPM1N. lugensUPM, Selangor, MalaysiaL. hexandraUPM2N. lugensTanjung Karang (Tk), Selangor, MalaysiaRiceTK1N. lugensTanjung Karang (Tk), Selangor, MalaysiaL. hexandraTK2N. lugensMalim, Melaka (Mk), MalaysiaRiceMK1N. lugensMalim, Melaka(Mk), MalaysiaL. hexandraMK2N. lugensBander Seberang, Perak (Pk), MalaysiaRicePK1N. lugensBander, Seberang, Perak (Pk), MalaysiaL. hexandraPK2N. lugensTuaran, Sabah (SB), MalaysiaRiceSB1N. lugensTuaran, Sabah (SB), MalaysiaL. hexandraSB2N. bakeriCameron Highlands (CH), Pahang, MalaysiaL. hexandraCH


### 2.2. Esterase Activity Test

The individual insect used for DALP analysis was tested for esterase activity on a simple filter paper using the method reported by Pasteur and Georghiou [13].


### 2.3. DNA Extraction

DNA from individual insect was extracted by grinding single frozen adult insect with a glass rod in a 1.5 mL tube containing 20μL extraction buffer (0.1 M Nacl, 0.2 M Sucrose, 0.1 M Tris-HCL (pH 9.0) 0.05 M EDTA, 0.5% SDS). The glass rod was washed with an additional 40 μL of extraction buffer and the homogenate was incubated at 65°C for 40 min. An amount of 10 μL of 8 M potassium acetate was added and the tube was placed on ice for 40 min. The tube was spun at 14000 rpm for 20 min. The supernatant was transferred into a fresh 1.5 mL tube. One hundred microliters of chilled (−20°C) 100% ethanol was added and the DNA was allowed to precipitate at room temperature for 10 min. The tube was spun for 20 min and the ethanol was carefully removed with a pipette. The DNA pellet was washed with 100 μL of chilled 70% ethanol and spun for 10 min. The DNA pellets were dried by pouring off the ethanol. The tubes were kept for 10 min at room temperature. The dried DNA pellet was suspended in 50 μL TE (Tris EDTA, pH 8.0) and gently mixed for 10 min. The DNA concentration was measured using LKB-Ultrastep III UV/visible spectrophotometer at the absorbance of 260 nm and 280 nm. The DNA was considered pure if the ratio of OD260/OD280 was within the range of 1.6–1.9 [12].


### 2.4. Primers Used in This Study

A total of three forward sequencing primer denoted as DALP231, DALP234, DALP235 and a universal reverse primer, DALPR (5′TTTCACACAGGAAACAGCTATGAC-3′), were used for the PCR amplification (Table2). Primer DALPR was end labeled with Y33 PATP (10 m Cie/mL) [14].Table 2
Optimisation of DALP primers used for the PCR protocols.
PrimersSequenceMg+2 conc.Taq polymerase usedAnnealing temperature (°C)DALP231F5′-GTTTTCCCAGTCACGACAGC-3′1.8 mMPromega55DALP234F5′-GTTTTCCCAGTCACGACCAG-3′1.8 mMPromega55DALP235F5′-GTTTTCCCAGTCACGACCAC-3′1.8 mMPromega55Universal DALPR5′-TTTCACACAGGAAACAGCTATGAC-3′


### 2.5. PCR Protocols

The PCR reaction mixture contained 60 ng of insect DNA, 1.8 mM of Mg+, 0.15 μM of oligonucletide primer, 200 μM of each dNTP, 1 unit Taq DNA polymerase (Promega), and 1x PCR buffer in a total volume of 25 μL. Amplification reactions were carried out in a programmable thermal cycler (GeneAmp, PCR system 2400, Perkin Elmer) programmed as follows: predenaturation at 95°C for 2 min; followed by 30 cycles of denaturation at 91°C for 30 sec, annealing temperature at 55°C for 45 sec, and extension at 70°C for 30 sec. After the last cycle, final extension was at 70°C for 5 min. The protocols were modified from [8].


### 2.6. Electrophoresis of the Multilocus Amplification Products

Electrophoresis was performed on 6% denaturing polyacrylamide gels and run on a 50 cm long gel apparatus. The samples were mixed with 5μL 100% formamide loading dye and then heated for 10 min at 96°C before loading. The gel was run at 55 W for 3 hours.


### 2.7. Autoradiography

After electrophoresis, the gel was transferred to a Whatman paper and dried and developed after 5 days of exposures to X-ray film.


### 2.8. Statistical Analysis

Band Scoring
DALP-PCR band profiles were scored visually for each DNA sample for each primer pair. The data was recorded according to the presence/absence criterion (1 = presence; 0 = absence of band).Cluster Analysis
The Dice algorithm was used for similarity index. The similarity index was calculated between two samples from within or between populations according to [15]:
(1)Sxy=2mxy(mx+my),
where mxy is the number of bands showed by sample x and sample y and mx and my are the number of bands in sample x and sample y, respectively. The value produced by this index ranges from 0 (representing no band sharing) to 1 (representing complete identity). The within or between population values are based on pairwise comparisons between individuals for a particular primer. The values obtained are then averaged over primers.
The between population similarity indices were also converted to distance values using the relationshipD=1-S [16, 17]. These distance matrices were used as the input matrix for the unweighted pair group method with arithmetic mean (UPGMA) tree [18] to find population relationships graphically using NTSYS-PC software (version 1.8; [19]).Test of Robustness
The test of robustness or bootstrapping was performed using the Phylogeny Inference Package (PHYLIP; version 3.5p) developed by Felsenstein [20]. The bootstrap values were obtained using gene frequencies option within the program PHYLIP. A consensus tree was produced based on the 1000 bootstrapped replicates as reported by Haymer et al. [21].Principle Component Analysis
A principal component analysis was performed based on the distance matrix among the populations using the NTSYS-PC software. The relationship among the populations is expressed in a three-dimensional graph based on the first three components.Analysis of Molecular Variance (AMOVA)
The distance between two samples was calculated according to the formula of Excoffier et al. [22]:
(2)D=N{1-(NxyN)},
where N is the total number of bands and Nxy is the number of bands shared by two samples. The resulting distance matrix was used in an AMOVA [22]. In the AMOVA, the sources of variation were divided into three nested levels: among the host types, among the populations within host types, and among individuals within populations. Mean square deviation was calculated by dividing sum of squared deviation by the degrees of freedom. The variance component was expressed as percentage. The significance of components of variance was tested by the random permutation.


## 2.1. Collection of Insect Populations

A total of 11 populations were collected from 6 locations. Two sympatric populations ofN. lugens, one from rice and other fromL. hexandra, were collected from each of five locations. The locations were Universiti Putra Malaysia (UPM), Tanjung Karang (TK), Melaka (MK), Perak (PK), and Sabah (SB). An out group, N. bakeri was also collected from Cameron Highlands (CH), Malaysia. Locations, host type, date of collection, population code are shown in Table 1. Each population consisted of twenty insects. All collected insects were frozen at −70°C for further use.Table 1
Host types, sites of collection, and coding for 11 populations ofNilaparvata spp.
SpeciesLocationsHost plantsPopulation codeN. lugensUniversiti Putra Malaysia (UPM), Selangor, MalaysiaRiceUPM1N. lugensUPM, Selangor, MalaysiaL. hexandraUPM2N. lugensTanjung Karang (Tk), Selangor, MalaysiaRiceTK1N. lugensTanjung Karang (Tk), Selangor, MalaysiaL. hexandraTK2N. lugensMalim, Melaka (Mk), MalaysiaRiceMK1N. lugensMalim, Melaka(Mk), MalaysiaL. hexandraMK2N. lugensBander Seberang, Perak (Pk), MalaysiaRicePK1N. lugensBander, Seberang, Perak (Pk), MalaysiaL. hexandraPK2N. lugensTuaran, Sabah (SB), MalaysiaRiceSB1N. lugensTuaran, Sabah (SB), MalaysiaL. hexandraSB2N. bakeriCameron Highlands (CH), Pahang, MalaysiaL. hexandraCH


## 2.2. Esterase Activity Test

The individual insect used for DALP analysis was tested for esterase activity on a simple filter paper using the method reported by Pasteur and Georghiou [13].


## 2.3. DNA Extraction

DNA from individual insect was extracted by grinding single frozen adult insect with a glass rod in a 1.5 mL tube containing 20μL extraction buffer (0.1 M Nacl, 0.2 M Sucrose, 0.1 M Tris-HCL (pH 9.0) 0.05 M EDTA, 0.5% SDS). The glass rod was washed with an additional 40 μL of extraction buffer and the homogenate was incubated at 65°C for 40 min. An amount of 10 μL of 8 M potassium acetate was added and the tube was placed on ice for 40 min. The tube was spun at 14000 rpm for 20 min. The supernatant was transferred into a fresh 1.5 mL tube. One hundred microliters of chilled (−20°C) 100% ethanol was added and the DNA was allowed to precipitate at room temperature for 10 min. The tube was spun for 20 min and the ethanol was carefully removed with a pipette. The DNA pellet was washed with 100 μL of chilled 70% ethanol and spun for 10 min. The DNA pellets were dried by pouring off the ethanol. The tubes were kept for 10 min at room temperature. The dried DNA pellet was suspended in 50 μL TE (Tris EDTA, pH 8.0) and gently mixed for 10 min. The DNA concentration was measured using LKB-Ultrastep III UV/visible spectrophotometer at the absorbance of 260 nm and 280 nm. The DNA was considered pure if the ratio of OD260/OD280 was within the range of 1.6–1.9 [12].


## 2.4. Primers Used in This Study

A total of three forward sequencing primer denoted as DALP231, DALP234, DALP235 and a universal reverse primer, DALPR (5′TTTCACACAGGAAACAGCTATGAC-3′), were used for the PCR amplification (Table2). Primer DALPR was end labeled with Y33 PATP (10 m Cie/mL) [14].Table 2
Optimisation of DALP primers used for the PCR protocols.
PrimersSequenceMg+2 conc.Taq polymerase usedAnnealing temperature (°C)DALP231F5′-GTTTTCCCAGTCACGACAGC-3′1.8 mMPromega55DALP234F5′-GTTTTCCCAGTCACGACCAG-3′1.8 mMPromega55DALP235F5′-GTTTTCCCAGTCACGACCAC-3′1.8 mMPromega55Universal DALPR5′-TTTCACACAGGAAACAGCTATGAC-3′


## 2.5. PCR Protocols

The PCR reaction mixture contained 60 ng of insect DNA, 1.8 mM of Mg+, 0.15 μM of oligonucletide primer, 200 μM of each dNTP, 1 unit Taq DNA polymerase (Promega), and 1x PCR buffer in a total volume of 25 μL. Amplification reactions were carried out in a programmable thermal cycler (GeneAmp, PCR system 2400, Perkin Elmer) programmed as follows: predenaturation at 95°C for 2 min; followed by 30 cycles of denaturation at 91°C for 30 sec, annealing temperature at 55°C for 45 sec, and extension at 70°C for 30 sec. After the last cycle, final extension was at 70°C for 5 min. The protocols were modified from [8].


## 2.6. Electrophoresis of the Multilocus Amplification Products

Electrophoresis was performed on 6% denaturing polyacrylamide gels and run on a 50 cm long gel apparatus. The samples were mixed with 5μL 100% formamide loading dye and then heated for 10 min at 96°C before loading. The gel was run at 55 W for 3 hours.


## 2.7. Autoradiography

After electrophoresis, the gel was transferred to a Whatman paper and dried and developed after 5 days of exposures to X-ray film.


## 2.8. Statistical Analysis

Band Scoring
DALP-PCR band profiles were scored visually for each DNA sample for each primer pair. The data was recorded according to the presence/absence criterion (1 = presence; 0 = absence of band).Cluster Analysis
The Dice algorithm was used for similarity index. The similarity index was calculated between two samples from within or between populations according to [15]:
(1)Sxy=2mxy(mx+my),
where mxy is the number of bands showed by sample x and sample y and mx and my are the number of bands in sample x and sample y, respectively. The value produced by this index ranges from 0 (representing no band sharing) to 1 (representing complete identity). The within or between population values are based on pairwise comparisons between individuals for a particular primer. The values obtained are then averaged over primers.
The between population similarity indices were also converted to distance values using the relationshipD=1-S [16, 17]. These distance matrices were used as the input matrix for the unweighted pair group method with arithmetic mean (UPGMA) tree [18] to find population relationships graphically using NTSYS-PC software (version 1.8; [19]).Test of Robustness
The test of robustness or bootstrapping was performed using the Phylogeny Inference Package (PHYLIP; version 3.5p) developed by Felsenstein [20]. The bootstrap values were obtained using gene frequencies option within the program PHYLIP. A consensus tree was produced based on the 1000 bootstrapped replicates as reported by Haymer et al. [21].Principle Component Analysis
A principal component analysis was performed based on the distance matrix among the populations using the NTSYS-PC software. The relationship among the populations is expressed in a three-dimensional graph based on the first three components.Analysis of Molecular Variance (AMOVA)
The distance between two samples was calculated according to the formula of Excoffier et al. [22]:
(2)D=N{1-(NxyN)},
where N is the total number of bands and Nxy is the number of bands shared by two samples. The resulting distance matrix was used in an AMOVA [22]. In the AMOVA, the sources of variation were divided into three nested levels: among the host types, among the populations within host types, and among individuals within populations. Mean square deviation was calculated by dividing sum of squared deviation by the degrees of freedom. The variance component was expressed as percentage. The significance of components of variance was tested by the random permutation.


## 3. Results

Band or marker frequency was calculated for each marker pair for each population for DALP primers. Figure1 shows the banding patterns obtained from rice and Leersia infesting populations of N. lugens using primer DALP235/DALPR. A hundred percent marker frequency represented monomorphism while 0% showed complete absence of the particular marker. The data showed a range of 28.3–42.9% polymorphic markers for rice infesting populations of N. lugens while Leersia infesting populations and an out group,N. bakeri showed 31.9–45.5% and 17.1% polymorphic markers, respectively. The overall data for 10 populations of N. lugens based on three primers showed 29 (42.6%) polymorphic markers. Frequency of DALP markers, total number of markers, number of polymorphic markers, % polymorphic markers for each population are shown in Table 3.Table 3
Frequency of presence of bands (%) markers inNilaparvata spp. obtained from three DALP primers.
PrimerMarker (M. Wt.)UPMTanjung KarangMelakaPerakSabahCameron High landsN.  lugens (Rice)N.  lugens  (Leersia)N.  lugens (Rice)N.  lugens  (Leersia)N.  lugens (Rice)N.lugens (Leersia)N.  lugens (Rice)N.lugens (Leersia)N. lugens (Rice)N.lugens (Leersia)N.bakeri (Leersia)DALP231.1 (2.15 kb)101000102010045100DALP231.2 (2.0 kb)025004030150501030DALP231.3 (1.9 kb)250004035151025100DALP231.4 (1.8 kb)100651009010090100807010080DALP231.5 (1.7 kb)10010010010010010010010010010090DALP231.6 (1.6 kb)1001001001001001001001001001000DALP231.7 (1.4 kb)601001008050100100607010090DALP231.8 (1.3 kb)100100100100100100100100100100100DALP231.9 (1.2 kb)1001001001001001001001001001000DALP231.10 (1.0 kb)060050040060040100DALP231DALP231.11 (0.9 kb)010100020000300DALP231.12 (0.85 kb)3040255030302030403580DALP231.13 (0.75 kb)0000000000100DALP231.14 (0.70 kb)0000000000100DALP231.15 (0.60 kb)0000000030300DALP231.16 (0.50 kb)0000000000100DALP231.17 (0.40 kb)100100100100100100100100100100100DALP231.18 (0.30 kb)1001001001001001001001001001000DALP231.19 (0.25 kb)100100100100100100100100100100100DALP231.20 (0.20 kb)80601001001001001001001001000DALP231.21 (0.10 kb)100100100100100100100100100100100DALP231.22 (0.09 kb)0000000000100DALP234.1  (2.0 kb)1001001001001001001001001001000DALP234.2  (1.9 kb)1001001001001001001001001001000DALP234.3  (1.8 kb)205025502050203030500DALP234.4  (1.75 kb)0000000000100DALP234.5 (1.70 kb)1001001001001001001001001001000DALP234.6 (1.65 kb)070040040060040100DALP234DALP234.7  (1.6 kb)010300020000300DALP234.8  (1.5 kb)40406050303020304035100DALP234.9  (1.45 kb)0000000000100DALP234.10  (1.3 kb)100100100100100100100100100100100DALP234.11 (1.25 kb)0000000030300DALP234.12  (1.10 kb)0000000000100DALP234.13  (1.0 kb)100100100100100100100100100100100DALP234.14  (0.9 kb)1001001001001001001001001001000DALP234.15 (0.85 kb)100100100100100100100100100100100DALP234.16 (0.7 kb)1001001001001001001001001001000DALP234.17 (0.65 kb)90901001001001009010010090100DALP234.18  (0.5 kb)0000000000100DALP234.19  (0.35 kb)40501001001001001001008090100DALP234.20  (0.30 kb)80208070604050907080100DALP234.21 (0.15 kb)405090100801001001008090100DALP235.1 (2.2 kb)51010801092040350DALP235.2 (2.1 kb)100100100100100100100100100100100DALP235.3  (2.0 kb)1001001001001001001001001001000DALP235.4  (1.9 kb)100100100100100100100100100100100DALP235.5  (1.8 kb)104015502035163010500DALP235.6  (1.75 kb)0000000000100DALP235.7 (1.70 kb)1001001001001001001001001001000DALP235.8 (1.65 kb)060050040060040100DALP235.9  (1.6 kb)010100020000300DALP235.10  (1.5 kb)30402550303020304035100DALP235.11  (1.4 kb)0000000000100DALP235.12  (1.3 kb)0000000000100DALP235.13 (1.2 kb)0000000030300DALP235.14 (1.15 kb)0000000000100DALP235.15  (1.0 kb)100100100100100100100100100100100DALP235DALP235.16  (0.9 kb)100100100100100100100100100100100DALP235.17 (0.8 kb)100100100100100100100100100100100DALP235.18 (0.7 kb)1001001001001001001001001001000DALP235.19 (0.6 kb)100100100100100100100100100100100DALP235.20  (0.5 kb)0000000000100DALP235.21  (0.35 kb)3050100100100100100100809070DALP235.22  (0.30 kb)90208070604050907080100DALP235.23  (0.20 kb)609050308010090604075100DALP235.24 (0.15 kb)100100801001001001001009010090DALP235.25 (0.10 kb)100100100100100100100100100100100Total Number of markers4551464645524647495541Number of polymorphic markers172314141318131521257% Polymorphic markers37.7745.0930.4330.4328.8834.6128.2631.9142.8545.4517.08Overall number of markers based on the populations ofN. lugens68Overall number of polymorphic markers based on the populations ofN. lugens29Overall % poymorphism based on the populations ofN. lugens42.64Frequency based upon 20 individuals per population.Figure 1
DALP-PCR amplicons obtained from rice andLeersia infesting populations of N. lugens using primer DALP235/DALPR (lanes 1–3 = rice infesting populations, lanes 5–7 = Leersia infesting populations; lane 4 = An out group, N. bakeri). Polymorphic markers showed in arrow sign.


### 3.1. Cluster Analysis

All data from three pairs of DALP primers were incorporated for cluster analysis. In this analysis, pairwise genetic distances were calculated between all individuals in order to make comparison. The distances within rice infesting populations ofN. lugens ranged from 0.112487 to 0.285200 (average 0.2245843) while distances within Leersia infesting populations ranged from 0.152379 to 0.235396 (average 0.2078274). The genetic distances between two sympatric populations of N. lugens, one from rice and the other from grass, ranged from 0.24019 to 0.390182 (average 0.31672).In addition to that genetic distances between rice infesting population ofN. lugens and Leersiainfesting populations of N. bakeri (out group) ranged from 0.555389 to 0.564963 (average 0.540276) but it was ranged from 0.499403 to 0.578171 (average 0.539168) between the populations of Leersia infesting N. lugens and N. bakeri (Table 4).Table 4
Genetic distance matrix of the 11 populations ofNilaparvata spp. based on Nei and Li’s similarity index.
UPM 1TK1MK1PK1SB1UPM2TK2MK2PK2SB2CHUPM10TK10.1124870MK10.1691360.1580450PK10.2816820.2852000.2440570SB10.2428070.2507720.2280620.2735950UPM20.2902650.2927170.2995150.3722070.3278290TK20.3002140.2852500.3131830.3781350.3433460.1523790MK20.2401930.2534260.2737730.3317070.2840360.1672850.1831290PK20.3060630.3163610.3112060.3384660.3265830.2089540.2228980.1971680SB20.3288760.3327740.346280.3901820.3356170.2353960.248590.2285970.2338780CH0.5553890.5090430.4973370.5746490.5649630.4994030.5259740.5546520.5376380.5781710UPM1, TK1, MK1, PK1, and SB1 are rice-infesting populations; UPM2, TK2, MK2, PK2, and SB2 areLeersia-infesting populations of N. lugens. CH is N. bakeri.UPGMA dendrogram revealed the genetic relationships among the 11 populations ofNilaparvata species. The cluster analysis divided the individuals into three main clusters. Among the three clusters, one was the most distinct and distant and the other two were closely related. All rice infesting populations like UPM1, MK1, TK1, PK1, and SB1 were included in one cluster, likewise the Leersia infesting populations such as UPM2, TK2, PK2, MK2, and SB2 were separated into another group. A common branch was shared by both groups. The isolated population CH (N. bakeri) was far away from either rice or Leersia infesting populations of N. lugens (Figure 2).Figure 2
UPGMA dendrogram of the 11 populations ofNilaparvata spp. based on genetic distance from Dice's index for DALP markers (Rice infesting population of N. lugens = UPM1, TK1, MK1, PK1, and SB1; Leersia-infesting population of N. lugens = UPM2, TK2, MK2, PK2, and SB2; an out group, N. bakeri = CH); bootstrap values from 1000 bootstraps are given at each fork.


### 3.2. Test of Robustness

The UPGMA tree was subjected to numerical resampling by bootstrapping [23] and the resultant bootstrap values were shown at the tree branch points. Each value represents the number of times that the represented groupings occurred in the resamplings. The consensus tree showed 100% confidence levels between rice (MK1, UPM1, TK1, PK1, and SB1) and Leersia infesting (MK2, TK2, UPM2, SB2, and PK2) population. Within rice and Leersia infesting populations, confidence level ranged from 37–56% to 35–51%, respectively (Figure 2). The confidence level between N. lugens and N. bakeri was also 100%.


### 3.3. Principal Component Analysis (PCA)

A principal component analysis was performed based on the distance matrix among the populations using the NTSYS-PC software. The relationship among the populations was expressed in a three-dimensional graph. In PCA graph, 11 populations were clustered into 3 groups. The cluster I consisted of rice infesting population ofN. lugens while Leersia infesting populations showed another group. The population of N. bakari showed an out group. The first three principal components accounted for 78.31% of the total variation among the 11 populations of Nilaparvata spp. and these 3 components, PC1, PC2, and PC3 showed 41.19, 27.42, 9.70% variation, respectively (Figure 3).Figure 3
Patterns of relationships of the 11 populations revealed by the principal component analysis based on short primer DALP data. Proportion of the total variance explained by the first three principal components (PCs) is 78.31%: PC1 = 41.19%; PC2 = 27.42%; PC3 = 9.70% (Rice infesting populations ofN.  lugens: 1 = UPM1, 2 = TK1, 3 = MK1, 4 = PK1 and 5 = SB1; Grass infesting populations of N.  lugens: 6 = UPM2, 7 = TK2, 8 = MK2, 9 = PK2, and 10 = SB2; An out group, N.  bakeri,  11 = CH).


### 3.4. Analysis of Molecular Variance (AMOVA)

Three level nested structures for each pair of primer of DALP are shown in Table5. All primers showed variance among the host types, among the populations, and among the individuals in a population. Out of three primers, DALP 235 determined the highest variance among the groups (rice versus Leersia) (26.90%), followed by DALP231 (10.10%), and DALP234 (9.68%). The percentage of the variance component among groups (rice versus Leersia) was greater than the percentage of the variance component among the populations detected by the three primers DALP235, DALP234 and DALP231. The results of AMOVA as well as dendrogram confirmed that genetic variation exists between the brown plant hopper of rice and Leersia.Table 5
Analysis of molecular variance (AMOVA) of 10 populations ofN. lugens based on three DALP primers.
PrimerSource of variationDegree of freedomSum of squared deviationMean squared deviationVariance component% TotalP
*DALP235/DALPRRice versusLeersia (Among groups)119.5016.810.15926.90<0.001Populations within group822.342.870.17226.29<0.001Individuals within population19049.500.290.28052.03<0.001DALP231/DALPRRice versusLeersia(Among groups)19.587.580.07510.10<0.001Populations within group89.992.000.0416.98<0.001Individuals within population19075.430.480.47485.06<0.001DALP234/DALPRRice versusLeersia (Among groups)13.043.090.0349.680.031Populations within group87.070.820.035.780.006Individuals within population19066.900.350.3895.02<0.001*After 1000 random permutations; P = Probability level.


## 3.1. Cluster Analysis

All data from three pairs of DALP primers were incorporated for cluster analysis. In this analysis, pairwise genetic distances were calculated between all individuals in order to make comparison. The distances within rice infesting populations ofN. lugens ranged from 0.112487 to 0.285200 (average 0.2245843) while distances within Leersia infesting populations ranged from 0.152379 to 0.235396 (average 0.2078274). The genetic distances between two sympatric populations of N. lugens, one from rice and the other from grass, ranged from 0.24019 to 0.390182 (average 0.31672).In addition to that genetic distances between rice infesting population ofN. lugens and Leersiainfesting populations of N. bakeri (out group) ranged from 0.555389 to 0.564963 (average 0.540276) but it was ranged from 0.499403 to 0.578171 (average 0.539168) between the populations of Leersia infesting N. lugens and N. bakeri (Table 4).Table 4
Genetic distance matrix of the 11 populations ofNilaparvata spp. based on Nei and Li’s similarity index.
UPM 1TK1MK1PK1SB1UPM2TK2MK2PK2SB2CHUPM10TK10.1124870MK10.1691360.1580450PK10.2816820.2852000.2440570SB10.2428070.2507720.2280620.2735950UPM20.2902650.2927170.2995150.3722070.3278290TK20.3002140.2852500.3131830.3781350.3433460.1523790MK20.2401930.2534260.2737730.3317070.2840360.1672850.1831290PK20.3060630.3163610.3112060.3384660.3265830.2089540.2228980.1971680SB20.3288760.3327740.346280.3901820.3356170.2353960.248590.2285970.2338780CH0.5553890.5090430.4973370.5746490.5649630.4994030.5259740.5546520.5376380.5781710UPM1, TK1, MK1, PK1, and SB1 are rice-infesting populations; UPM2, TK2, MK2, PK2, and SB2 areLeersia-infesting populations of N. lugens. CH is N. bakeri.UPGMA dendrogram revealed the genetic relationships among the 11 populations ofNilaparvata species. The cluster analysis divided the individuals into three main clusters. Among the three clusters, one was the most distinct and distant and the other two were closely related. All rice infesting populations like UPM1, MK1, TK1, PK1, and SB1 were included in one cluster, likewise the Leersia infesting populations such as UPM2, TK2, PK2, MK2, and SB2 were separated into another group. A common branch was shared by both groups. The isolated population CH (N. bakeri) was far away from either rice or Leersia infesting populations of N. lugens (Figure 2).Figure 2
UPGMA dendrogram of the 11 populations ofNilaparvata spp. based on genetic distance from Dice's index for DALP markers (Rice infesting population of N. lugens = UPM1, TK1, MK1, PK1, and SB1; Leersia-infesting population of N. lugens = UPM2, TK2, MK2, PK2, and SB2; an out group, N. bakeri = CH); bootstrap values from 1000 bootstraps are given at each fork.


## 3.2. Test of Robustness

The UPGMA tree was subjected to numerical resampling by bootstrapping [23] and the resultant bootstrap values were shown at the tree branch points. Each value represents the number of times that the represented groupings occurred in the resamplings. The consensus tree showed 100% confidence levels between rice (MK1, UPM1, TK1, PK1, and SB1) and Leersia infesting (MK2, TK2, UPM2, SB2, and PK2) population. Within rice and Leersia infesting populations, confidence level ranged from 37–56% to 35–51%, respectively (Figure 2). The confidence level between N. lugens and N. bakeri was also 100%.


## 3.3. Principal Component Analysis (PCA)

A principal component analysis was performed based on the distance matrix among the populations using the NTSYS-PC software. The relationship among the populations was expressed in a three-dimensional graph. In PCA graph, 11 populations were clustered into 3 groups. The cluster I consisted of rice infesting population ofN. lugens while Leersia infesting populations showed another group. The population of N. bakari showed an out group. The first three principal components accounted for 78.31% of the total variation among the 11 populations of Nilaparvata spp. and these 3 components, PC1, PC2, and PC3 showed 41.19, 27.42, 9.70% variation, respectively (Figure 3).Figure 3
Patterns of relationships of the 11 populations revealed by the principal component analysis based on short primer DALP data. Proportion of the total variance explained by the first three principal components (PCs) is 78.31%: PC1 = 41.19%; PC2 = 27.42%; PC3 = 9.70% (Rice infesting populations ofN.  lugens: 1 = UPM1, 2 = TK1, 3 = MK1, 4 = PK1 and 5 = SB1; Grass infesting populations of N.  lugens: 6 = UPM2, 7 = TK2, 8 = MK2, 9 = PK2, and 10 = SB2; An out group, N.  bakeri,  11 = CH).


## 3.4. Analysis of Molecular Variance (AMOVA)

Three level nested structures for each pair of primer of DALP are shown in Table5. All primers showed variance among the host types, among the populations, and among the individuals in a population. Out of three primers, DALP 235 determined the highest variance among the groups (rice versus Leersia) (26.90%), followed by DALP231 (10.10%), and DALP234 (9.68%). The percentage of the variance component among groups (rice versus Leersia) was greater than the percentage of the variance component among the populations detected by the three primers DALP235, DALP234 and DALP231. The results of AMOVA as well as dendrogram confirmed that genetic variation exists between the brown plant hopper of rice and Leersia.Table 5
Analysis of molecular variance (AMOVA) of 10 populations ofN. lugens based on three DALP primers.
PrimerSource of variationDegree of freedomSum of squared deviationMean squared deviationVariance component% TotalP
*DALP235/DALPRRice versusLeersia (Among groups)119.5016.810.15926.90<0.001Populations within group822.342.870.17226.29<0.001Individuals within population19049.500.290.28052.03<0.001DALP231/DALPRRice versusLeersia(Among groups)19.587.580.07510.10<0.001Populations within group89.992.000.0416.98<0.001Individuals within population19075.430.480.47485.06<0.001DALP234/DALPRRice versusLeersia (Among groups)13.043.090.0349.680.031Populations within group87.070.820.035.780.006Individuals within population19066.900.350.3895.02<0.001*After 1000 random permutations; P = Probability level.


## 4. Discussion

Cluster and principal component analyses revealed the genetic relationships among the different populations ofNilaparvata species. Three major clusters were observed in the dendrogram as well as in the graph. The results showed that population of N. bakeri formed the most isolated cluster from populations of either rice or Leersia infesting populations ofN. lugens. The rice infesting populations of UPM, Tanjung Karang, Melaka, Perak, and Sabah, Malaysia clustered together as a group. On the other hand, Leersiainfesting populations of the same localities formed another distinct cluster.Leersia infesting populations with low esterase activities seem to be formed, a different structure from rice infesting populations of brown planthopper,N. lugens. These results were also confirmed by bootstrapping analysis as described by Felsenstein [23] and Latif et al. [6]. Bootstrapping was initially used to evaluate the accuracy of a tree obtained by the parsimony method and could increase the confidence level of the results obtained from the DALP assay. The results showed 100% confidence level for the separate clusterings between the rice and Leersia infesting populations of N. lugens and also for the genetically isolated group, N. bakeri.In DALP fingerprinting method, we did not get any diagnostic markers between two sympatric populations ofN. lugens. Saxena and Barrion [24] reported that karyoytpe, idiogram, nuclear organelles, chromosomes with nucleolus organizing region (site of RNA synthesis) showed clear differences between rice and Leersia infesting populations. Despite their morphological similarities, a distinct cytological incongruity and a certain degree of genetic isolation between the two populations were inferred. Species differentiation in early stages of a species formation may not be associated with substantial genetic change [25–27]. Many ecologists have accepted that the evolutionary processes are common in animals with specialized food habits [28, 29]. There was no distinct electrophoretic differentiation between Lethe eurydice L. and Lethe appalacia L., although the two species were found to be good species [27, 30]. Latif et al. [6] reported that the closely related sibling species in the N. lugens complex might have developed through insecticide exposures that were heavier in rice-infesting populations than in grass populations, through RAPD-PCR analysis.The genetic distance indicates the magnitude of genetic variation between populations. Genetic distance commonly ranged from nearly 0.01 for populations within species, 0.1 for different subspecies, and 1.0 for different species [31]. So, the genetic distances (average 0.31672) between rice infesting populations (high esterase activities) and Leersia infesting populations (low esterase activities) of brown planthoppers indicated that these sympatric populations represented two distinct but closely related biological species.The results of AMOVA in single primer yielded highly significant variance among group (rice versusLeersia) and among population components. The total genetic variation, an average 15.56% was attributable to group divergence (Rice versus Leersia), 13.01% to population differences and 77.37% to individual differences within a population. The percentage of variance component among groups (rice versus Leersia) was larger than the percentage of variance component among populations for bands detected by three DALP primers and these were tested by random permutation. These results revealed that there was genetic differentiation between the brown planthopper of rice versus Leersia (two sympatric populations of N. lugens). AMOVA was performed and was confirmed the differentiation into two groups of Aphid gossypii [32], five groups of Acorus gramineus [33], and two groups of natural populations of the wild rice, Oryza rufipogon [34]. Therefore, our molecular data of DALP-PCR indicated that brown plant hopper (BPH) with high esterase activity usually captured from rice plant and those with low esterase activity, usually captured fromL. hexandra in Malaysia, represent two distinct closely related species and supported previous results as reported by Latif et al. [6, 7, 35]. Although DALP molecular method is not new, but so far to our knowledge this study is the first to detect genetic polymorphism in rice brown planthopper complex using this method.



---

*Source: 586831-2012-04-19.xml*