Page 4 - Fiorentino_alii_2008
P. 4

812   V. FIORENTINO ET AL.

             stored in 95% ethanol for morphological examination  treated as unknown and were estimated during the
             and kept in the F. Giusti Collection (Dipartimento di  run. The separate genes were treated as ‘unlinked’, so
             Scienze Ambientali, Siena, Italy).                that separate parameter estimates were obtained for
                                                               each gene. The default value of four Markov chains
                                                               was used (one cold and three heated) and each chain
               DNA EXTRACTION, POLYMERASE CHAIN REACTION
                                                               was started from a random tree. The ‘temperature’
                          (PCR) AND SEQUENCING                 parameter was set at 0.2. The Markov Chain Monte
             Total genomic DNA was extracted from foot muscle of  Carlo length was 4 000 000 generations and trees
             fresh specimens using C-TAB buffer (0.1 M Tris-HCl,  were sampled every 100 generations.
             pH 8.0, 1.4 M NaCl, 0.02 M ethylenediaminetetracetic  Parsimony analysis was performed with PAUP*,
             acid, 2% cetyl trimethylammonium bromide, 0.2%    version 4.0 (Swofford, 2001) using a heuristic search
             2-mercaptoethanol) and standard phenol–chloroform/  under equal weighting of all characters (1000 random
             ethanol extraction (Hillis et al., 1996).         stepwise addition and tree bisection-reconnection
               A fragment of CO I was amplified by PCR using    branch swapping). To evaluate clade support, non-
             the specifically-designed primer pair (5′-CCTATTATA  parametric bootstrap resampling was used with 1000
             ATTGGGGGTTTTGG-3′) and (5′-GTATCGGCTCTAA          replicate datasets and 100 random-addition sequ-
             AATAAGC-3′). A fragment of the mitochondrial gene  ences per dataset.
             coding for the large ribosomal subunit (16S rRNA)
             was also amplified, using the primer combination
             LR-J-12887 [5′-CGATTTGAACTCAGATCA-3′ (Simon                         RESULTS
             et al., 1994) and 5′-GTGCAAAGGTAGCATAATCA-3′
                                                                          MORPHOLOGICAL ANALYSIS
             (Gantenbein et al., 1999)].
               PCR reactions were run in a total volume of 50 mL  Regarding the DFA of shell, the first discriminant
             under the following conditions: 95 °C for 20 s, 48 °C  function accounted for most morphological variance
             for 30 s and 72 °C for 30 s (repeated for 25-s cycles),  (74%) and the second for only 12%. Considering the
             plus a final extension step at 72 °C for 5 min. Ampli-  canonical and structural coefficients making up the
             fied products were examined on 1.5% agarose gels.  first function, the highest loadings (Table 1) were for
             DNA fragments were then excised and purified with a  maximum diameter (MaxD) (1.42) and shell height
             ‘Nucleospin extract’ (Genenco) column kit, following  (ShH) (-1.29) respectively. A sequential chi-square
             the manufacturer’s instructions.                  test showed that the first function contributed to
               Detecting the appropriate amplification products  population discrimination to a very large extent
                                                                 2
             was straightforward for both genes because no sec-  (c = 3506, d.f. = 319, P < 0.001) whereas the contribu-
                                                                                2
             ondary bands due to nonspecific priming occurred.  tion of the second (c = 1981, d.f. = 280, P < 0.001) and
             Approximately 30 mg of the purified products was   subsequent functions was less significant. Indeed, the
             then used as templates for direct sequencing of both  first function discriminated scores from negative to
             strands, using the same amplification primers and
             a CEQ dye terminator cycle sequencing kit. DNA
             sequences then underwent electrophoresis on a CEQ  Table 1. Canonical coefficients of discriminant analysis of
             8000XL (Beckman Coulter). All sequences generated  shell measurements
             in this study were deposited in GenBank (see
             Appendix).                                                              Canonical variables
                                                               Variable              First              Second
                          GENETIC DATA ANALYSIS
                                                               AH                    -0.27              -0.06
             For both mitochondrial regions, sequences were easily
                                                               AW                     0.45              0.44
             aligned unambiguously by eye using the computer   MaxD                   1.42              0.98
             program BIOEDIT, version 7.0 (Hall, 1999).
                                                               MinD                   0.21              0.10
               Phylogenetic relationships were estimated by Baye-  HMaxD             -0.61              -0.19
             sian inference (BI) and maximum parsimony (MP).
                                                               LWW                   -0.26              -0.15
             Prior to Bayesian inference, we used MRMODELT-    ShH                   -1.29              -0.56
             EST, version 2.2 (Nylander, 2004) to determine the  SpH                  0.30              0.23
             best model of nucleotide evolution for each dataset.  MaxLWH            -0.76              0.12
             The selected models were then used in a partitioned  a                  -0.07              -0.17
             Bayesian analysis carried out with MRBAYES,       b                     -0.11              0.12
             version 3.1 (Huelsenbeck & Ronquist, 2001; Ronquist
             & Huelsenbeck, 2003). Model parameter values were  For abbreviations, see Fig. 1.

                                © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 94, 809–823


  Downloaded from https://academic.oup.com/biolinnean/article-abstract/94/4/809/2701256
  by guest
  on 18 November 2017
   1   2   3   4   5   6   7   8   9