Comparison among Methods and Statistical Software Packages to Analyze Germplasm Genetic Diversity by Means of Codominant Markers

Co-dominant markers’ data are often analysed as if they were dominant markers, an over-simplification that may be misleading. Addressing this, the present paper aims to provide a practical guide to the analysis of co-dominant data and selection of suitable software. An overview is provided of the computational methods and basic principles necessary for statistical analyses of co-dominant molecular markers to determine genetic diversity and molecular characterization of germplasm collections. The Hardy–Weinberg principle is at the base of statistical methods to determine genetic distance, genetic diversity, and its distribution among and within populations. Six statistical software packages named GenAlEx, GDA, Power Marker, Cervus, Arlequin, and Structure are compared and contrasted. The different software packages were selected based on: (i) The ability to analyze co-dominant data, (ii) open access software, (iii) ease of downloading, and (iv) ease of running using a Microsoft Window interface. The software packages are compared analyzing the same dataset. Differences among parameters are discussed together with the comments on some of the software outputs.


Introduction
Genetic diversity of germplasm is assessed by collecting key information, especially: (i) Allele number per locus; (ii) genotype number per locus; (iii) gene diversity; (iv) PIC (polymorphism information content) values; (v) observed and expected heterozygosity; (vi) partition of the diversity into its components within and between populations; and (vii) the genetic distance among the analyzed populations.The analyses are usually performed using a variety of molecular markers grouped into two categories: Co-dominant markers, such as SSR (single sequence repeat) and SNP (single nucleotide polymorphism), which are able to identify the allelic situation at each locus, and dominant markers, such as ISSR (inter simple sequence repeats), RAPD (random amplified polymorphic DNA), and AFLP (amplified fragment length polymorphism), which usually have a multi-band pattern and are unable to recognize allelic variants [1].The latter produce a series of bands with unknown relationships (i.e., could be allelic variants of the same genes or mark different genome regions).Hence, without knowing the allelic situation, each band is recorded as a locus with two possible alleles' band presence (scored as 1) or band absence (scored as 0) and the relative 0/1 matrix is used in statistical analyses.The papers reviewed here comprise data based on co-dominant markers that were often wrongly recorded as the presence/absence of possible bands, leading to a loss of information on allelic variance and the presence of heterozygosity (observed heterozygosity, Ho).
The present paper offers a short and simple guide to the principles that form the base of the most common analyses.It focuses on some of the most widely-used computer programs in population genetics, run under Windows, to highlight the advantages and disadvantages of the various software packages, thus facilitating appropriate selection and use.

Hardy-Weinberg Principle
Most of the statistical computations use parameters based on the Hardy-Weinberg principle [2,3].Here, the basis of the principle and its applications are highlighted.As it is widely known, the Hardy-Weinberg principle considers the genetic and genotype frequency for a single locus in a population and states: "allele and genotype frequencies in a population will remain constant from generation to generation in the absence of other evolutionary influences".These potential evolutionary forces include: (i) Migration, (ii) mutation, (iii) selection, (iv) population size sufficient to avoid drift, and (v) random mating.Unfortunately, this definition of the Hardy-Weinberg does not sufficiently focus on other important consequences of the principle such as: "if a population is in equilibrium it is possible to compute the allele frequencies knowing the genotype frequencies and vice-versa by the formula of binomial square development i.e., (p + q) 2 = p 2 + q 2 + 2pq = 1", where p 2 is the frequency of the AA genotype, q 2 indicates the aa genotype frequency, 2pq the Aa genotype frequency, p the A allele frequency, and q the a allele frequency.This equation is true only for a population in the Hardy-Weinberg equilibrium where it is possible to compute allele frequencies from knowing the genotype frequencies and vice versa.The above is if only two alleles, A and a, are possible for that locus.If, instead, three alleles may occur at a locus, the formula would be a trinomial square development ((p + q + r) 2 = p 2 + q 2 + r 2 + 2pq + 2pr + 2qr = 1) and so on for higher numbers of alleles.It should be noted that the square terms (i.e., p 2 + q 2 + r 2 , etc.) are homozygote frequencies while the others (i.e., 2pq + 2pr + 2qr, etc.) are heterozygotes.Considering several alleles, I, with a frequency, p i , the homozygote frequency is Σp i 2 and heterozygote frequency can be calculated as the complementary difference from the homozygote frequency (i.e., 2pq = 1 − (p 2 + q 2 ) or 1 − Σp i 2 ).

Genetic Diversity
The gene diversity index is calculated for each locus and population according to Nei [4], utilizing the Hardy-Weinberg formula, He = 1 − ∑ n i=1 p 2 i , hereafter simplified as He = 1 − Σp i 2 , which is the heterozygosity expected if the population is in Hardy-Weinberg equilibrium.In analogy, the genetic identity (J) is Σp i 2 (homozygotes).However, since He could be computed for all populations, including non-random mating systems (e.g., autogamus, which, by definition, will not in Hardy-Weinberg equilibrium being a pure line with homozygosity for all loci), the terminology for He is thus gene diversity, rather than expected heterozygosity.
In a small population, the alleles per locus can be skewed, especially when compared to large populations [5].Unbiased heterozygosity is as for the above-mentioned heterozygosity multiplied by the factor, 2n/(2n − 1) [6].As a result, the larger the population, the lower are the differences between the biased and unbiased expected heterozygosity.This detail is often not sufficiently elaborated upon in the literature, as many papers do not mention whether unbiased or biased He is used.
The variability between and within populations can be calculated according to Nei [4] by taking into account different allele frequencies in whole populations or only in subpopulations.The nomenclature used is: H T for total observed diversity; H S for within-population diversity; and D ST for the between-population diversity, with H T = H S + D ST .
Similarly, the Wright's fixation indices, F IS , F ST , and F IT [7], are often used, also the F-statistics are based on the expected level of heterozygosity.The measures describe the different levels of population structures, such as variance of allele frequencies within populations (F IS ), variance of allele frequencies between populations (F ST ), and an inbreeding coefficient of an individual relative to the total population (F IT ), all of which are related to heterozygosity at various levels of population structure.The terms mentioned above are represented by the formula, 1 − F IT = 1 − F IS + 1 − F ST , where I is the individual, S the subpopulation, and T the total population.F IT thus refers to the individual in comparison with the total, F IS is the individual in comparison with the subpopulation, and F ST is the subpopulation in comparison with the total.As shown in Figure 1, total F, indicated by F IT , can be partitioned into F IS (or f ) and F ST (or θ).
comparison with the total, FIS is the individual in comparison with the subpopulation, and FST is the subpopulation in comparison with the total.As shown in Figure 1, total F, indicated by FIT, can be partitioned into FIS (or f) and FST (or θ).FST can be calculated using the formula: FST = (HT − HS)/HT, where HT is the proportion of the heterozygotes in the total population and HS the average proportion of heterozygotes in subpopulations.
In a series of loci, l, in n populations and using the complementary sum of allele frequency (1 − Ʃpi 2 ), different figures can be obtained.In particular: 1.For each locus and each population, He = (1 − Ʃpi(lg) 2 ), where pi(lg) is the ith allele frequency of the lth locus in the gth population.2. The average of the above He over populations gives the genetic diversity within a population for each locus, while the average of all the loci within a population diversity gives HS.The formula can thus be written as: HS = (Ʃl(Ʃg(1 − Ʃpi(lg) 2 )/g)/l), where (1 − Ʃpi(lg) 2 ) indicates the expected heterozygosity for each locus in each population, g indicates the number of populations, and l the loci number.3. The total genetic diversity, HT, is calculated using the allele frequency, pi(l), for each locus over all populations and calculating the mean over loci: HT = Ʃ(1 − Ʃpi(l) 2 )/l).4. The between population component of diversity is calculated using the formula: DST = HT − HS. 5.The between population component may also be expressed in relation to the total genetic diversity (for each locus and overall loci) as GST = HT/DST [4].Table 1 shows an example extracted from Turpeinen et al. [8], where different parameters for three populations were analyzed using two markers.The HT for each locus corresponds to the polymorphic information content (PIC) of that locus, which in other words, consists in the capacity of that locus (or better a marker) to assess polymorphism and diversity.Botstein et al. [9] proposed an adjustment of this value as: where pi and pj are the population frequency of the ith and jth alleles.The PIC proposed by Botstein and colleagues [9] subtracts from the He value an additional probability (ƩƩ2pi 2 pj 2 ) due to the fact that linked individuals do not add information to the overall variation.F ST can be calculated using the formula: F ST = (H T − H S )/H T , where H T is the proportion of the heterozygotes in the total population and H S the average proportion of heterozygotes in subpopulations.
In a series of loci, l, in n populations and using the complementary sum of allele frequency (1 − Σp i 2 ), different figures can be obtained.In particular: 1.
For each locus and each population, He = (1 − Σp i(lg) 2 ), where p i(lg) is the ith allele frequency of the lth locus in the gth population.

2.
The average of the above He over populations gives the genetic diversity within a population for each locus, while the average of all the loci within a population diversity gives H S .The formula can thus be written as: 2 )/ g )/ l ), where (1 − Σp i(lg) 2 ) indicates the expected heterozygosity for each locus in each population, g indicates the number of populations, and l the loci number.

3.
The total genetic diversity, H T , is calculated using the allele frequency, p i(l) , for each locus over all populations and calculating the mean over loci: 2 )/ l ).

4.
The between population component of diversity is calculated using the formula: The between population component may also be expressed in relation to the total genetic diversity (for each locus and overall loci) as G ST = H T /D ST [4].Table 1 shows an example extracted from Turpeinen et al. [8], where different parameters for three populations were analyzed using two markers.The H T for each locus corresponds to the polymorphic information content (PIC) of that locus, which in other words, consists in the capacity of that locus (or better a marker) to assess polymorphism and diversity.Botstein et al. [9] proposed an adjustment of this value as: where p i and p j are the population frequency of the ith and jth alleles.The PIC proposed by Botstein and colleagues [9] subtracts from the He value an additional probability (ΣΣ2p i 2 p j 2 ) due to the fact that linked individuals do not add information to the overall variation.

Genetic Distance
Genetic diversity (He) and genetic identity (J or Ho) are also used to estimate the genetic distance within and between populations, since two populations with high identity in their genes are closer than two with high diversity.If J x = Σp xi 2 is the probability of identity in population x with p xi the frequency of the i-th allele and J y = Σp yi 2 is the probability of identity in population y, the probability of identity in both populations is J xy = Σp xi p yi as described by Nei [10,11].The probability of identity in population x for all normalized loci is I = J xy / √ (J x J y ) and, in turn, the genetic distance is D = −LnI = −Ln (J xy / √ (J x J y )).In a small sample set with many loci, any biases can be corrected using , where G x and G y are (2n x J x − 1)/(2n x − 1) and (2n y J y − 1)/(2n y − 1) over the l loci studied, respectively, and G xy = J xy [12].In this case, Ď could be negative, due to sampling errors, and hence considered as zero.

Data Input
One of the first issues is the data format required as the various software packages use different data-file formats.A difference of a single comma or space can make the data unreadable or misclassified.As a result, it often takes more time to organize the data into the correct format than to run the analysis itself.Some programs may offer the possibility of importing/exporting data from/to other formats, thereby avoiding reformatting data manually and making it easier and faster to analyze a given data set with different programs.This is of particular importance where the data set may require the use of more than one application and/or analyses offered by different software packages.
The amplicons generated from markers are distinguished by submarine gel electrophoresis or a capillary in a sequencer; in the later cases, the results, as alleles call, can be exported from the sequencer into a Microsoft Excel file.Excel seems to be the easiest and most universal way to insert data.As such, GenAlEx [13], which is an Excel macro rather than a full software package, is first to be considered.GenAlEx software, as its 6.5 version [21], can be downloaded from http://biology.anu.edu.au/GenAlEx/Download.html and has a template function for co-dominant, binary, and haploid data, creating a framework on which data insertion can be easily carried out starting from the cell C4.After the data are inserted, they can be analyzed directly by GenAlEx or alternatively be exported to other formats specific to other commonly used statistical software.The present example entails seven loci, 90 samples, nine populations, and three regions, which are indicated in the template (Figure 2a).
The results are stored in an Excel sheet where the loci and the populations are indicated with consecutive numbering; it is possible, however, to change these to the correct locus and population names.Being co-dominant data, each locus will have two columns for the two alleles (Figure 2b).GenAlEx can also be used to import or export data from or to other software packages, although it is very important to pay attention to the codes used by the different software to indicate missing data.For example, the alleles can be easily named with their molecular weight in bp, however, the null allele (which is not missing data) could be named as zero, but zero is considered missing for some software, such as GenAlEx, when co-dominance is the option selected.In these cases, it is important to rename the null allele, for example, by substituting zero with 1.

Data Input
One of the first issues is the data format required as the various software packages use different data-file formats.A difference of a single comma or space can make the data unreadable or misclassified.As a result, it often takes more time to organize the data into the correct format than to run the analysis itself.Some programs may offer the possibility of importing/exporting data from/to other formats, thereby avoiding reformatting data manually and making it easier and faster to analyze a given data set with different programs.This is of particular importance where the data set may require the use of more than one application and/or analyses offered by different software packages.
The amplicons generated from markers are distinguished by submarine gel electrophoresis or a capillary in a sequencer; in the later cases, the results, as alleles call, can be exported from the sequencer into a Microsoft Excel file.Excel seems to be the easiest and most universal way to insert data.As such, GenAlEx [13], which is an Excel macro rather than a full software package, is first to be considered.GenAlEx software, as its 6.5 version [21], can be downloaded from http://biology.anu.edu.au/GenAlEx/Download.html and has a template function for co-dominant, binary, and haploid data, creating a framework on which data insertion can be easily carried out starting from the cell C4.After the data are inserted, they can be analyzed directly by GenAlEx or alternatively be exported to other formats specific to other commonly used statistical software.The present example entails seven loci, 90 samples, nine populations, and three regions, which are indicated in the template (Figure 2a).
The results are stored in an Excel sheet where the loci and the populations are indicated with consecutive numbering; it is possible, however, to change these to the correct locus and population names.Being co-dominant data, each locus will have two columns for the two alleles (Figure 2b).GenAlEx can also be used to import or export data from or to other software packages, although it is very important to pay attention to the codes used by the different software to indicate missing data.For example, the alleles can be easily named with their molecular weight in bp, however, the null allele (which is not missing data) could be named as zero, but zero is considered missing for some software, such as GenAlEx, when co-dominance is the option selected.In these cases, it is important to rename the null allele, for example, by substituting zero with 1.

Data Analysis
The same data was then analyzed using various software packages and the various outputs compared and reported here.

GenAlEx
GenAlEx is available at http://biology-assets.anu.edu.au/GenAlEx/Welcome.html, as mentioned above.It is an Excel macro used for statistical genetic analysis, so the user should be registered for an Office package which is not open source.By using the "Frequency…" option, it is possible to compute allele frequency, heterozygosity, F-stat, and polymorphism by population and by locus, some genetic distances (i.e., Nei distance, Nei unbiased distance, pairwise FST) together with some graphic options (Figure 3).

Data Analysis
The same data was then analyzed using various software packages and the various outputs compared and reported here.

GenAlEx
GenAlEx is available at http://biology-assets.anu.edu.au/GenAlEx/Welcome.html, as mentioned above.It is an Excel macro used for statistical genetic analysis, so the user should be registered for an Office package which is not open source.By using the "Frequency . . ." option, it is possible to compute allele frequency, heterozygosity, F-stat, and polymorphism by population and by locus, some genetic distances (i.e., Nei distance, Nei unbiased distance, pairwise F ST ) together with some graphic options (Figure 3).
One of the positive aspects of GenAlEx is that the different output-sheets display the base of the statistic used.There are also options for graphics (i.e., Allele Frequencies by Population with Graph over Loci or Graphs by Population and Locus) that provide a quick overview of allele distribution among populations.The most important outputs are in the sheets "HFP" and "HFL", where the different statistical parameters by locus (Table 2) and/or by populations (Table 3) are provided.The parameters are: One of the positive aspects of GenAlEx is that the different output-sheets display the base of the statistic used.There are also options for graphics (i.e., Allele Frequencies by Population with Graph over Loci or Graphs by Population and Locus) that provide a quick overview of allele distribution among populations.The most important outputs are in the sheets "HFP" and "HFL", where the different statistical parameters by locus (Table 2) and/or by populations (Table 3) are provided.The parameters are: Where tpi is the frequency of the ith allele for the total and Σtpi 2 is the sum of the squared total allele frequencies.Where tp i is the frequency of the ith allele for the total and Σtp i 2 is the sum of the squared total allele frequencies.The three levels of the fixation indexes (F IS , F IT , F ST ) are computed per locus and not per population as in other programs, such as Arlequin (see below).
The output of different genetic distances, such as the Nei's distance, Nei's unbiased distance, and pairwise F ST , are reported in Table 4. GenAlEx can calculate the molecular analysis of variance (AMOVA), which partitions genetic variability into different components (Table 5), including, or not, the individual level.

GDA
GDA can be downloaded at http://en.bio-soft.net/dna/gda.htmlor now at https://phylogeny.uconn.edu/software/.Data can be exported from GenAlEx to GDA, but it is necessary to manually change the file extension.A useful tool of GDA is the possibility of easily re-running the analysis excluding/including loci and/or populations.
The descriptive statistics offered by GDA are: (i) Number of alleles per population (A), (Na in GenAlEx); (ii) polymorphic alleles per locus, (not in GenAlex); (iii) expected (He); and (iv) observed (Ho) heterozygosity.Observed heterozygosity is in line with GenAlEx output, while the He is here the unbiased expected heterozygosity (uHe in GenAlEx).GDA outputs per population and per locus are reported in Table 6.Table 7 shows the private alleles, another useful option present in GDA.In Table 8, genetic distances computed in agreement with Nei (1972) [10] and Nei (1978) [12] are shown; the first is the unbiased genetic distance of GenAlEx, while the second is equal to the genetic distance reported in GenAlEx.The graphic output is as a text file.To improve the options for the quality of graphs, it is necessary to use other software, such as TreeView [22].The graphic quality and options are not considered here as it is the ability of the statistic software to export the dendrogram codes to be then used in the graphical software that is of prime importance.

Popgene
Popgene offers two versions for either 32 or 16 bit Windows operating environments, and can be downloaded at https://sites.ualberta.ca/~fyeh/popgene.html.It immediately divided the analysis depending on whether it deals with dominant or codominant markers.For diploid data, it performs a genotypic frequency, HW test (not commonly found in other packages), fixation index, allele The graphic output is as a text file.To improve the options for the quality of graphs, it is necessary to use other software, such as TreeView [22].The graphic quality and options are not considered here as it is the ability of the statistic software to export the dendrogram codes to be then used in the graphical software that is of prime importance.Nei (1978) [12] distance; below diagonal Nei (1972) [10]

Popgene
Popgene offers two versions for either 32 or 16 bit Windows operating environments, and can be downloaded at https://sites.ualberta.ca/~{}fyeh/popgene.html.It immediately divided the analysis depending on whether it deals with dominant or codominant markers.For diploid data, it performs a genotypic frequency, HW test (not commonly found in other packages), fixation index, allele frequency, allele number, effective allele number, polymorphic loci, observed and expected homozygosity and heterozygosity, Shannon index, homogeneity test, F-statistics (F IT , F ST , F IS ), gene flow, and genetic distance (following Nei 1972 [10] and Nei 1978 [6]).It also produces a dendrogram using UPGMA of the Nei's distance, neutrality test, and the linkage disequilibrium (LD) between two loci.In the cases of several alleles per locus, the required input is not straightforward, based on the Mendelian convention (Figure 5), i.e., providing a letter for each allele, but it is possible to export the Popgene format from GenAlEx.However, a significant disadvantage is that it assigns the same letter to alleles from different loci, as if they were the same allele.This creates confusions and errors especially when reading the tables of "Allele Frequency".

Power Marker
Power Marker, like GDA, was developed at the North Carolina State University and uses as a reference the Genetic Data Analysis by Weir [23].The original download source for Power Markers, http://www.powermarker.net/[16], seems to be expired, however, the program and the manual can be found at http://statgen.ncsu.edu/powermarker/index.html.
Data input is very easy, entering the allelic phase separated by space, tab, and/or commas.It is possible to indicate up to three category levels.In this example, we used: Genotype, populations,

Power Marker
Power Marker, like GDA, was developed at the North Carolina State University and uses as a reference the Genetic Data Analysis by Weir [23].The original download source for Power Markers, http://www.powermarker.net/[16], seems to be expired, however, the program and the manual can be found at http://statgen.ncsu.edu/powermarker/index.html.
Data input is very easy, entering the allelic phase separated by space, tab, and/or commas.It is possible to indicate up to three category levels.In this example, we used: Genotype, populations, and regions.The program is suitable for microsatellite data; however, it also works with haplotypes.The data can be reduced by a sub-selection of genotypes or markers based on particular parameters, such as the level of missing data, heterozygosity, or diversity.Outputs have their own format, which can be easily converted into Excel files.A very useful tool is the internal link with the TreeView [24] graphic program (http://taxonomy.zoology.gla.ac.uk/rod/treeview.html)used to display genotype relationships (Trees) with good graphical resolution.However, to use this function, the user must also install the TreeView program.
The summary table (Figure 6) illustrates information, such as: (i) Allele frequency, (ii) genotype number, (iii) number of observations, (iv) allele, (v) gene diversity, (vi) heterozygosity, and (vii) PIC.Number of observations, allele, gene diversity, heterozygosity, and PIC are equivalent to the values reported in GDA, respectively, as n, A, He, and Ho.In Power Marker, the expected heterozygosity (which is not unbiased expected heterozygosity as in GDA) is named "gene diversity".It should be noted that the PIC values are here computed according to Botstein et al. [9].The main disadvantage of Power Marker is that outputs always refer to the markers rather than to the population as per GDA.To show values per population, it is necessary to create a subset of data where only one population is considered each time.Another disadvantage is that the output does not report the options chosen, so naming the folders with self-explaining labels is an imperative.After the user has computed the allele frequency by using the "phylogeny" option, it is possible to calculate the frequency based distance utilizing several methods.The only equivalent method to the other software in this paper is Nei's genetic distance 1972 [10].In addition, Power Maker can compute the pairwise linkage disequilibrium, where the output is displayed for each marker in the order they were inserted in the data file (Figure 7).Therefore, it is crucial that the marker results be entered in the "right" order, which is important only if a genetic map with marker positions along the chromosomes is available.After the user has computed the allele frequency by using the "phylogeny" option, it is possible to calculate the frequency based distance utilizing several methods.The only equivalent method to the other software in this paper is Nei's genetic distance 1972 [10].In addition, Power Maker can compute the pairwise linkage disequilibrium, where the output is displayed for each marker in the order they were inserted in the data file (Figure 7).Therefore, it is crucial that the marker results be entered in the "right" order, which is important only if a genetic map with marker positions along the chromosomes is available.ter the user has computed the allele frequency by using the "phylogeny" option, it is late the frequency based distance utilizing several methods.The only equivalent m er software in this paper is Nei's genetic distance 1972 [10].In addition, Power M e the pairwise linkage disequilibrium, where the output is displayed for each mark ey were inserted in the data file (Figure 7).Therefore, it is crucial that the marker re in the "right" order, which is important only if a genetic map with marker positio mosomes is available.vus rvus is primarily designed for the assignment of parents to their offspring using s.Nevertheless, it is sometimes used for genetic analysis.It is available for dow eldgenetics.com.The input data sheet is not as user-friendly as some of the other pr can be converted from GenePop, which in turn can be converted from GenAlEx.alculates the PIC value as per Botstein et al. [9] and He is unbiased.In crossed pop Table 9. Cervus output reporting the number of alleles per locus (k), number of individuals (N), observed (Hobs) and expected (Hexp) heterozygosity, PIC, combined non-exclusion probability for first parent (NE-1P), second parent (NE-2P), parent pair (NE-PP), identity (NE-I) and sib identity (NE-SI), the Hardy-Weinberg equilibrium significance (HW), and the F test (F).

Arlequin
Arlequin, available at http://cmpg.unibe.ch/software/arlequin3/,produces output displayed in a browser page, and thus is not ideal for conversion into a word document.On the other hand, the particular computation run by Arlequin is AMOVA (analyses of molecular variance) as described by Excoffier et al. [25].It considers haplotype, and with 90 genotypes, the total degree of freedom is 179 10).The AMOVA output is very similar to the GenAlEx one (Table 5).Where: T , F ST = 0.252 = F IT , F SC = 0.240 = F IS , and F CT = 0.015 = F ST .
He and Ho are reported for each locus within each population and produce the same average outcome as the GDA software.Linkage disequilibrium, where the deviation from random association between alleles at different loci [26], expressed as D = p ij − p i p j , is a potentially useful additional feature of Arlequin.However, although the instruction manual asserts the computation of the linkage disequilibrium coefficient (D) is possible, this seems not to be true.On the contrary, significance is reported as the P values of χ 2 with 1000 permutations.Moreover, the number of loci linked to each locus for each population analyzed is provided.Unfortunately, even when the locus name is inserted, it is not reflected in the output, where the loci are simply numbered starting at zero.Similarly, the populations are numbered as pop1#, pop2#, pop3#, etc. rather than using the given name.This could easily lead to mistakes and confusion.In addition, there are sometimes discrepancies between the data saved in the browser output file and that saved as an xls file.

Structure
Structure software [19] is available for downloaded at http://web.stanford.edu/group/pritchardlab/structure.html.Preparation of the data file in order to run Structure presents some problems.Conversion from GenAlEx is not straightforward since (i) an extra space is required at the end of the second row to allow the program to read the last number, and (ii) population names are not converted automatically.Moreover, particular care must be taken when dealing with missing data and their code, for doing so differs from other software packages (in Structure, "−9" is used as default, but it is possible to set it differently).However, with suitable modification, it is easy to convert files directly from Excel by saving it as a text file.
In Structure, the analysis should be set in agreement with the populations' information and the procedures used in the population sampling.Useful information to assist clustering includes three possible options: (i) Considering individuals with or without common ancestry, (ii) with or without use of sampling locations, and (iii) to set the allele frequencies as either independent or dependent in each population.
Fundamentally, Structure performs a K-mean cluster analyses [27].As with all K-mean cluster methods, in Structure, the analysis should be performed trying different values for the number of clusters (K).Clearly, some logical pre-cluster division can be argued in agreement with the data typologies, number of populations, regions, groups, etc.Nevertheless, several runs with different K values should be performed and compared.Moreover, it is sometimes useful to run single populations alone to test if they include different subpopulations.Evanno [28] has suggested the use of ∆K in order to aid determination of the correct number of clusters.This should help in most situations, but should not be used as an exclusive criterion.The STRUCTUREHARVESTER software, available on line [29], has been developed to determine Evanno computation.
The Structure output can be displayed as a "triangle plot" in which two clusters are plotted at two vertices and all the others at the third (Figure 8).When more than three clusters are obtained, it can require further classification.However, a more useful, frequently used output is the bar plot on which the clusters are shown using different colors that can then be divided to highlight populations, or can be sorted by the Q value (Figure 9a,b).The picture gives a clear idea of how the individuals are divided among clusters/populations, and hence, the population similarity and the collections structure.Structure can provide the histograms of Fst, alpha, and likelihood for each cluster, as well as a tree plot of the distance among clusters.It is also possible to plot the average proportion of the Q values directly on a geographic map [30].
In Structure, the analysis should be set in agreement with the populations' information and the procedures used in the population sampling.Useful information to assist clustering includes three possible options: (i) Considering individuals with or without common ancestry, (ii) with or without use of sampling locations, and (iii) to set the allele frequencies as either independent or dependent in each population.
Fundamentally, Structure performs a K-mean cluster analyses [27].As with all K-mean cluster methods, in Structure, the analysis should be performed trying different values for the number of clusters (K).Clearly, some logical pre-cluster division can be argued in agreement with the data typologies, number of populations, regions, groups, etc.Nevertheless, several runs with different K values should be performed and compared.Moreover, it is sometimes useful to run single populations alone to test if they include different subpopulations.Evanno [28] has suggested the use of ∆K in order to aid determination of the correct number of clusters.This should help in most situations, but should not be used as an exclusive criterion.The STRUCTUREHARVESTER software, available on line [29], has been developed to determine Evanno computation.
The Structure output can be displayed as a "triangle plot" in which two clusters are plotted at two vertices and all the others at the third (Figure 8).When more than three clusters are obtained, it can require further classification.However, a more useful, frequently used output is the bar plot on which the clusters are shown using different colors that can then be divided to highlight populations, or can be sorted by the Q value (Figure 9a,b).The picture gives a clear idea of how the individuals are divided among clusters/populations, and hence, the population similarity and the collections structure.Structure can provide the histograms of Fst, alpha, and likelihood for each cluster, as well as a tree plot of the distance among clusters.It is also possible to plot the average proportion of the Q values directly on a geographic map [30].

Conclusions
The different software packages available often use different methods and tools to describe populations.Table 11 provides an overview of the main function available for each program assessed in this paper.Overall, the author recommends using GenAlEx and/or Power Marker to insert data, subsequently exporting/importing and converting as required.In addition, either GDA and/or Power Marker can be used to perform most of the statistical analyses required for measuring genetic diversity, such as the percent of polymorphism, allele number, polymorphic allele number, and the expected and observed heterozygosity.In GDA, these parameters refer both to the loci and to the populations, while in Power Marker, several subsets of data should be run per population.Power Marker also computes PIC values, while GDA also computes private alleles.Both programs have different methodologies for computing population distance.Finally, GenAlEx and Arlequin are useful for determining analyses of molecular variance and Structure provides a clear illustration of population clustering.

Conclusions
The different software packages available often use different methods and tools to describe populations.Table 11 provides an overview of the main function available for each program assessed in this paper.Overall, the author recommends using GenAlEx and/or Power Marker to insert data, subsequently exporting/importing and converting as required.In addition, either GDA and/or Power Marker can be used to perform most of the statistical analyses required for measuring genetic diversity, such as the percent of polymorphism, allele number, polymorphic allele number, and the expected and observed heterozygosity.In GDA, these parameters refer both to the loci and to the populations, while in Power Marker, several subsets of data should be run per population.Power Marker also computes PIC values, while GDA also computes private alleles.Both programs have different methodologies for computing population distance.Finally, GenAlEx and Arlequin are

Figure 1 .
Figure 1.Diagram of the relationships between the gene diversity components.I = individual, S = subpopulation, T = total population.

Figure 1 .
Figure 1.Diagram of the relationships between the gene diversity components.I = individual, S = subpopulation, T = total population.

Figure 2 .
Figure 2. Structure of the data inserted by GenAlEx, the Excel macro for genetic analyses.(a) Template, (b) data in D sheet.

Figure 2 .
Figure 2. Structure of the data inserted by GenAlEx, the Excel macro for genetic analyses.(a) Template; (b) data in D sheet.

J
2018, 1, x FOR PEER REVIEW 14 of 19

Figure 6 .
Figure 6.Power Marker output for the genetic data.

Figure 6 .
Figure 6.Power Marker output for the genetic data.

Figure 6 .
Figure 6.Power Marker output for the genetic data.

Figure 7 .
Figure 7. Power Marker output for the pairwise linkage disequilibrium.

Figure 7 .
Figure 7. Power Marker output for the pairwise linkage disequilibrium.

Figure 8 .
Figure 8. Structure output of the triangle plot with the relationships among clusters.

Figure 8 .Figure 9 .
Figure 8. Structure output of the triangle plot with the relationships among clusters.

Figure 9 .
Figure 9. Structure output of the bar plot clusters either by populations (a) or sorted by the Q value (b) are reported with different colors.

Table 1 .
[8]elic situation and computation of the genetic parameters in three populations analyzed using two markers where each one has three possible alleles; adapted from Turpeinen et al.[8].

Table 2 .
GenAlEx output of the data in Figure2per locus.Sheet HFL.

Table 3 .
GenAlEx output of the data in Figure2per population.Sheet HFP.

Table 4 .
Computation of different parameters of distance between populations.Sheets NeiP, uNeiP, and F ST P. (A) Nei's genetic distance [10]; (B) Pairwise Population Matrix of Nei's Unbiased Genetic Distance; (C) Pairwise Population F ST Values.

Table 6 .
Descriptive statistics output of GDA per population (A) and per locus (B).Where n is the number of observations, P the polymorphism, A the alleles number, Ap the polymorphic alleles number, He the expected heterozygosity, and Ho the observed heterozygosity.

Table 8 .
Genetic distances computed by GDA.Above diagonal distance.

Table 11 .
Comparison of different characteristics of most frequently used software.