Diversidad y control biológico de insectos - CyberTesis UACh ...

Diversidad y control biológico de insectos - CyberTesis UACh ... Diversidad y control biológico de insectos - CyberTesis UACh ...

cybertesis.uach.cl
from cybertesis.uach.cl More from this publisher
31.08.2013 Views

products were precipitated adding ethanol. The sequencing reactions were suspended in deionized formamide, heat denatured and run on an ABI 3100 Genetic Analyzer (Applied Systems, Foster City, CA). The resulting chromatographs were imported into Sequencher 4.1 (Gene Codes Corp., Ann Arbor, Michigan) for visual inspection and editing. Multiple sequence alignments were constructed with the MegAlign module of DNASTAR 5 (Lasergene, Madison, WI). Data analysis. Phylogenetic approaches. The ninety seven sequences were aligned using the Clustal W module of MEGA version 3.1 (Kumar et al., 2004) and collapsed into haplotypes using the DNASP 4.10 software (Rozas et al., 2003). The phylogenetic relationships among the whole sequences and the extracted haplotypes were evaluated by using neighbour-joining (NJ), maximum parsimony (MP) and minimum evolution (ME) methods. All analyses were conducted using the Phylogeny module of MEGA version 3.1 (Kumar et al., 2004). The NJ tree was constructed from the matrix of pairwise p distances. This distance measure was chosen because it has low variance and then it is suitable to deal with low variation sequences (Nei and Kumar, 2000). Both MP and ME trees were constructed using unweighted data. Bootstrapping values (600 pseudoreplicates) were used to determine nodal support. The trees were rooted using a sequence obtained from the rice blast fungus Magnaporthe rosea, retrieved from Genbank (XM_368948, http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=Nucleotide&list_uids=3 9975114&dopt=GenBank). We considered gapped positions as unreliable characters and excluded them from further analysis (Swofford et al, 1996). To further analyze the genetic relatedness of haplotypes, we constructed a haplotype network using the median-joining network method as implemented in Network software (Bandelt et al., 1999; www.fluxus- engineering.com). Analyses of demographic history. The mismatch analysis (the distribution of the observed number of differences between pairs of haplotypes) was adopted to examine the demographic history of the species 23

(Schneider and Excoffier, 1999; Althoff and Pellmyr, 2002). The resulting distribution was compared against a simulated distribution assuming the constant growth model as implemented in DNASP (Rozas et., 2003). A unimodal distribution is expected if the lineages have undergone a recent bottleneck or population expansion, while a multimodal or ragged distribution is expected for a lineage whose populations are at demographic equilibrium (Althoff and Pellmyr, 2002). In addition, we used Tajima’s D to examine further the historical demography of B. bassiana sensu lato (Althoff and Pellmyr, 2002). This test for selective neutrality can be used to examine demography: a significant negative value indicate a deviation from the expectations of the mutation-drift equilibrium and can indicate population expansion, while a positive value is expected under population subdivision. Under neutrality the number of nucleotides differences between sequences from a random sample should be equal to the number of differences between the polymorphic sites only, but population expansions can cause negative departures of Tajima’s D from zero. The molecular diversity of the data set was calculated using the Arlequin 3.01 software (Schneider et al., 2000), including the average number of nucleotides per sites (nucleotide diversity π; Nei, 1987) and haplotype diversity. These parameters were calculated because centers of origin should be more diverse than more recently founded populations (Althoff and Pellmyr, 2002). Analyses of population structure. Hierarchical structuring of genetic variation was determined using analysis of molecular variance (AMOVA, Excoffier et al., 1992) as implemented in the Arlequin 3.01 software (Schneider et al., 2000). This analysis computes fixation indices φST analogous to Wright’s (1951) F statistics to divide the variance into the different components. The significance of the a priori geographical groupings was tested by bootstrapping (1023 permutations). The molecular distances between pairs of sequences necessary to run AMOVA were computed under the Kimura 2P model. 24

(Schnei<strong>de</strong>r and Excoffier, 1999; Althoff and Pellmyr, 2002). The resulting distribution was<br />

compared against a simulated distribution assuming the constant growth mo<strong>de</strong>l as<br />

implemented in DNASP (Rozas et., 2003). A unimodal distribution is expected if the<br />

lineages have un<strong>de</strong>rgone a recent bottleneck or population expansion, while a multimodal<br />

or ragged distribution is expected for a lineage whose populations are at <strong>de</strong>mographic<br />

equilibrium (Althoff and Pellmyr, 2002).<br />

In addition, we used Tajima’s D to examine further the historical <strong>de</strong>mography of B.<br />

bassiana sensu lato (Althoff and Pellmyr, 2002). This test for selective neutrality can be<br />

used to examine <strong>de</strong>mography: a significant negative value indicate a <strong>de</strong>viation from the<br />

expectations of the mutation-drift equilibrium and can indicate population expansion, while<br />

a positive value is expected un<strong>de</strong>r population subdivision. Un<strong>de</strong>r neutrality the number of<br />

nucleoti<strong>de</strong>s differences between sequences from a random sample should be equal to the<br />

number of differences between the polymorphic sites only, but population expansions can<br />

cause negative <strong>de</strong>partures of Tajima’s D from zero.<br />

The molecular diversity of the data set was calculated using the Arlequin 3.01 software<br />

(Schnei<strong>de</strong>r et al., 2000), including the average number of nucleoti<strong>de</strong>s per sites (nucleoti<strong>de</strong><br />

diversity π; Nei, 1987) and haplotype diversity. These parameters were calculated because<br />

centers of origin should be more diverse than more recently foun<strong>de</strong>d populations (Althoff<br />

and Pellmyr, 2002).<br />

Analyses of population structure.<br />

Hierarchical structuring of genetic variation was <strong>de</strong>termined using analysis of molecular<br />

variance (AMOVA, Excoffier et al., 1992) as implemented in the Arlequin 3.01 software<br />

(Schnei<strong>de</strong>r et al., 2000). This analysis computes fixation indices φST analogous to Wright’s<br />

(1951) F statistics to divi<strong>de</strong> the variance into the different components. The significance of<br />

the a priori geographical groupings was tested by bootstrapping (1023 permutations). The<br />

molecular distances between pairs of sequences necessary to run AMOVA were computed<br />

un<strong>de</strong>r the Kimura 2P mo<strong>de</strong>l.<br />

24

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!