Genetic Diversity of Hesperis persica, Based On Inter-Retrotransposon Amplified Polymorphism
Pages 322-336
https://doi.org/10.61882/JCT.16.4.322
somayeh esfandani bozchaloyi
Abstract Introduction: Knowledge of spatial genetic structures provides a valuable tool for inferring the evolutionary forces such as selective pressures and drift. Low gene flow due to spatial isolation of populations may even increase the degree of local differentiation. Nevertheless, phenotypic plasticity rather than genetic differentiation may be an alternative way of matching genotypes to environment; indeed increasing environmental variation favors higher levels of plasticity.Genetic diversity is one aspect of biological diversity that is extremely important for conservation strategies, especially in rare and narrowly endemic species. The genus Hesperis L. (Brassicaceae) comprises biennial and perennial herbs and consists of 46 species worldwide, mainly occurring in different parts of Europe, Caucasus, Transcaucasia, and to a lesser extent in northern and central Asia, and mostly in Turkey with 28 species. The genus is represented by 11 or six species belonging to sections Hesperis Dvořák, Diaplictos Dvořák and Pachycarpos Fourn. in Iran.
Aim: Moreover, due to extensive morphological variability of this species in the country, there is possibility of having infra-specific taxonomic forms in this species. Therefore, we carried out population genetic analysis and morphometric study of 11 geographical populations for the first time in the country. For genetic study, we used the inter-retrotransposon amplified polymorphism (IRAP) method that displays insertional polymorphisms by amplifying the segments of DNA between two retrotransposons. It has been used in numerous studies of genetic diversity. The objectives of this research were to study genetic diversity Hesperis persica with a different geographical origin by inter-retrotransposon amplified polymorphism (IRAP) method.
Materials and Methods: A total of 73 individuals were sampled representing 11 natural populations of H. persica Boiss. subsp. persica and H. persica subsp. kurdica (F. Dvořák & Hadac) F. Dvořák, in Mazandaran, East Azerbaijan, Kohgilouye-Boirahmad, Chaharmahal Bakhtiari, Fars, Zanjan, Tehran, Kordestan, Provinces of Iran during July-Agust 2019-2024.
Fresh leaves were used randomly from 6-10 plants in each of the studied populations. These were dried by silica gel powder. CTAB activated charcoal protocol was used to extract genomic DNA. The quality of extracted DNA was examined by running on 0.8% agarose gel. A set of six outward-facing LTR primers were used for IRAP analysis. We also used 15 different combinations of outward-facing LTR pair primers. PCR reactions were carried in a 25μl volume containing 10 mM Tris-HCl buffer at pH 8; 50 mM KCl; 1.5 mM MgCl2; 0.2 mM of each dNTP (Bioron, Germany); 0.2 μM of a single primer; 20 ng genomic DNA and 3 U of Taq DNA polymerase (Bioron, Germany). The thermal program was carried out with an initial denaturation for 1 min at 94°C, followed by 40 cycles in three segments: 35 s at 95°C, 40s at 47°C and 55s at 72°C. Final extension was performed at 72°C for 5 min. The amplification products were observed by running on 1% agarose gel, followed by the ethidium bromide staining. The fragment size was estimated by using a 100 bp molecular size ladder (Fermentas, Germany).The IRAP profiles obtained for each samples were scored as binary characters. Parameter like Nei’s gene diversity (H), Shannon information index (I), number of effective alleles, and percentage of polymorphism were determined.Nei’s genetic distance among populations was used for Neighbor Joining (NJ) clustering and Neighbor-Net networking. Mantel test checked the correlation between geographical and genetic distance of the studied populations . These analyses were done by PAST ver. 2.17 , DARwin ver. 5 (2012) and SplitsTree4 V4.13.1 (2013) software. AMOVA (Analysis of molecular variance) test (with 1000 permutations) as implemented in GenAlex 6.4, and Nei,s Gst analysis as implemented in GenoDive ver.2 (2013) were used to show genetic difference of the populations. Moreover, populations, genetic differentiation was studied by G'ST est = standardized measure of genetic differentiation, and D_est = Jost measure of differentiation. The genetic structure of populations was studied by Bayesian based model STRUCTURE analysis , and maximum likelihood-based method of K-Means clustering of GenoDive ver. 2. (2013). For STRUCTURE analysis, data were scored as dominant markers. The Evanno test was performed on STRUCTURE result to determine proper number of K by using delta K value. In K-Means clustering, two summary statistics, pseudo-F, and Bayesian Information Criterion (BIC), provide the best fit for k.
Results: The highest value of percentage polymorphism (57.41%) was observed in Mazandaran, 30 km S. of Ramsar between Kash-e Chal mountain and Miankuh (population No.8, H. persica subsp. kurdica) which shows high value for gene diversity (0.34) and Shanon, information index (0.43). Population Chaharmahal Bakhtiari, Shahr-e Kurd, tang-e Sayyad protected area, Pir kuh mountain (No.3, H. persica subsp. persica) has the lowest value for percentage of polymorphism (28.11%) and the lowest value for Shanon, information index (0.088), and He (0.022).AMOVA (PhiPT = 0.98, P = 0.010), and Gst analysis (0.654, p = 0.001) revealed significant difference among the studied populations. It also revealed that, 40% of total genetic variability was due to within population diversity and 60% was due to among population genetic differentiation. Pairwise AMOVA produced significant difference among the studied populations. Moreover, we got high values for Hedrick standardized fixation index after 999 permutation (G’st = 0.654, P = 0.001) and Jost, differentiation index (D-est = 0.769, P = 0.001). These results indicate that the geographical populations of Hesperis persica are genetically differentiated from each other. The mean Nm = 0.455 was obtained for all IRAP loci, which indicates low amount of gene flow among the populations. Population assignment test also agreed with Nm result and could not identify significant gene flow among these populations. However, reticulogram obtained based on the least square method, revealed some amount of shared alleles among populations 2 and 3, and between 7 and 8, also between 1, and 4. This result is in agreement with grouping we obtained with PCoA plot, as these populations were placed close to each other. As evidenced by STRUCTURE plot based on admixture model, these shared alleles comprise very limited part of the genomes in these populations and all these results are in agreement in showing high degree of genetic stratification within Hesperis persica populations. In total 76 IRAP bands (loci) were obtained, out of which 14 bands were private. Populations 2 and 5-7 contained 2-5 private bands.
Conclusion: PCoA plot of populations was in agreement with WARD clustering of molecular data. These results indicated that geographical populations of Hesperis persica are well differentiated based on (IRAP) markers.

This work is licensed under a