Field-based artificial crossings indicate partial compatibility of reciprocal crosses between <Emphasis Type="Italic">Pinus sylvestris</Emphasis> and <Emphasis Type="Italic">Pinus mugo</Emphasis> and unexpected chloroplast DNA inheritance

Crossability relationships between Scots pine (Pinus sylvestris L.) and mountain dwarf pine (Pinus mugo Turra) was studied, using artificial pollination approach. Partial compatibility of the reciprocal crossings of these species was proved experimentally, validating the idea of a spontaneous formation of their hybrid swarms under natural conditions. The hybrids were validated using organellar DNA markers and nuclear DNA microsatellites. Based on the percentage of filled seeds, the interspecific crossings were less efficient than the intraspecific cross-pollinations of P. sylvestris and P. mugo individuals. Both species were found to intercross readily with individuals of their putative hybrid swarm, P. mugo exhibiting a higher hybridological affinity towards putatively hybrid individuals than P. sylvestris. Validation of the hybrids confirmed the paternal inheritance of chloroplast DNA (cpDNA) in the combination P. sylvestris × P. mugo only. Surprisingly, in the reciprocal crossing P. mugo × P. sylvestris, maternal inheritance of cpDNA was revealed. Obtained results offer a new insight into the direction and intensity of gene flow within the hybrid swarms of Scots pine and mountain dwarf pine.     

Genetic diversity of the African poplar (<Emphasis Type="Italic">Populus ilicifolia</Emphasis>) populations in Kenya

We evaluated the genetic diversity of the African poplar (Populus ilicifolia) populations found in Kenya compared with reference samples of five poplar species from North America and one species introduced in Kenya from India (KEFRI-Kenya). Amplified fragment length polymorphism (AFLP) was used with the objective of providing important information for breeding and in situ/ex situ conservation of this species. Samples collected from three locations along the species’ natural range (Athi, Ewaso Nyiro, and Tana rivers) were compared with four samples of locally planted Populus deltoides stand introduced from India and ten reference samples from North America. Six AFLP primer combinations produced 521 clear bands for analysis. The percentage polymorphic loci were lowest in Tana (20.4 %) and highest in Athi (40.6 %). The average heterozygosity across the studied populations was between 0.07 and 0.3. AMOVA revealed more genetic variation partitioning within population (87 %; P?<?0.01) than among populations (13 %; P?<?0.01) suggesting significant genetic variation between populations. Further, UPGMA delineation showed two clusters of the Tana, Athi, and Ewaso Nyiro populations clustered together compared to the North America and India/KEFRI reference samples. Moreover, the study showed that the Athi population is more diverse than those of Tana and Ewaso Nyiro and may be important for conservation, domestication, and improvement studies. The genetic differentiation (F _ST ?=?0.134) among Kenyan P. ilicifolia populations suggests limited possibility of gene flow between these populations.     

Investigating the genetic diversity and differentiation patterns in the <Emphasis Type="Italic">Penstemon scariosus</Emphasis> species complex under different sample sizes using AFLPs and SSRs

Habitat fragmentation due to anthropogenic activities is the major cause of biodiversity loss. Endemic and narrowly distributed species are the most susceptible to habitat degradation. Penstemon scariosus is one of many species whose natural habitat is vulnerable to industrialization. All varieties of P. scariosus (P. scariosus var. albifluvis, P. scariosus var. cyanomontanus, P. scariosus var. garrettii, P. scariosus var. scariosus) have small distribution ranges, but only P. scariosus var. albifluvis is being considered for listing under the Endangered Species Act. We used eight microsatellites or simple sequence repeats (SSRs) loci and two amplified fragment length polymorphism (AFLP) primer combinations to investigate the population genetic structure and diversity of P. scariosus varieties. Moreover, we compared the utility of the two marker systems in conservation genetics and estimated an appropriate sample size in population genetic studies. Genetic differentiation among populations based on F_st ranged from low to moderate (F_st?=?0.056–0.157) and from moderate to high when estimated with D_es (D_es?=?0.15–0.32). Also, AMOVA analysis shows that most of the genetic variation is within populations. Inbreeding coefficients (F_is) were high in all varieties (0.20–0.56). The Bayesian analysis, STRUCTURE, identified three clusters from SSR data and four clusters from AFLPs. Clusters were not consistent between marker systems and did not represent the current taxonomy. MEMGENE revealed that a high proportion of the genetic variation is due to geographic distance (R²?=?0.38, P?=?0.001). Comparing the genetic measurements from AFLPs and SSRs, we found that AFLP results were more accurate than SSR results across sample size when populations were larger than 25 individuals. As sample size decreases, the estimates become less stable in both AFLP and SSR datasets. Finally, this study provides insight into the population genetic structure of these varieties, which could be used in conservation efforts.     

Identification and genetic diversity of two invasive <Emphasis Type="Italic">Pissodes</Emphasis> spp. Germar (Coleoptera: Curculionidae) in their introduced range in the southern hemisphere

During the first half of the twentieth century, two accidental cases of introduction of Pissodes weevils were recorded from the southern hemisphere. The weevils in South Africa were identified as the deodar weevil (Pissodes nemorensis) and those in South America as the small banded pine weevil (Pissodes castaneus). Wide distribution of the two species in their invasive range, general difficulty in identifying some Pissodes spp., and the varying feeding and breeding behaviours of the species in South Africa has necessitated better evidence of species identity and genetic diversity of both species and population structure of the species in South Africa. Barcoding and the Jerry-to-Pat region of the COI gene were investigated. Morphometric data of the South African species was analysed. Our results confirmed the introduction of only one Pissodes species of North American origin to South Africa. However, this species is not P. nemorensis, but an unrecognized species of the P. strobi complex or a hybrid between P. strobi and P. nemorensis. Only P. castaneus, of European origin, was identified from South America. We identified ten mitochondrial DNA haplotypes from South Africa with evidence of moderate genetic structure among geographic populations. Terminal leader and bole-feeding weevils did not differ at the COI locus. A single haplotype was identified from populations of P. castaneus in South America. Results of the present study will have implications on quarantine, research and management of these insect species.     

Genetic diversity and structure of the globally invasive tree, <Emphasis Type="Italic">Paraserianthes lophantha</Emphasis> subspecies <Emphasis Type="Italic">lophantha</Emphasis>, suggest an introduction history characterised by varying propagule pressure

An emerging insight in invasion biology is that intra-specific genetic variation, human usage, and introduction histories interact to shape genetic diversity and its distribution in populations of invasive species. We explore these aspects for the tree species Paraserianthes lophantha subsp. lophantha, a close relative of Australian wattles (genus Acacia). This species is native to Western Australia and is invasive in a number of regions globally. Using microsatellite genotype and DNA sequencing data, we show that native Western Australian populations of P. lophantha subsp. lophantha are geographically structured and are more diverse than introduced populations in Australia (New South Wales, South Australia, and Victoria), the Hawaiian Islands, Portugal, and South Africa. Introduced populations varied greatly in the amount of genetic diversity contained within them, from being low (e.g. Portugal) to high (e.g. Maui, Hawaiian Islands). Irrespective of provenance (native or introduced), all populations appeared to be highly inbred (F _IS ranging from 0.55 to 0.8), probably due to selfing. Although introduced populations generally had lower genetic diversity than native populations, Bayesian clustering of microsatellites and phylogenetic diversity indicated that introduced populations comprise a diverse array of genotypes, most of which were also identified in Western Australia. The dissimilarity in the distribution and number of genotypes in introduced regions suggests that non-native populations originated from different native sources and that introduction events differed in propagule pressure.     

The Role of Geographical and Ecological Factors on Population Divergence of the Neotropical otter <Emphasis Type="Italic">Lontra longicaudis</Emphasis> (Carnivora,Mustelidae)

The geographic distribution of the populations of a species are influenced by the spatial structure of the ecosystems, the environmental factors and the presence of geographic barriers. The Neotropical otter, Lontra longicaudis, is widely distributed throughout the Americas, where a wide range of environmental conditions and geographical features could promote genetic and morphological variation on the three currently recognized subspecies. In this study, we combined phylogeographic, morphometric and environmental niche modelling analyses to examine whether: (1) genetic variation is associated with the presence of barriers to gene flow and/or hydrography; (2) genetic and morphologic variation are associated with environmental variation; and (3) the observed variation in L. longicaudis populations corresponds to the previously defined subspecies. We found strong phylogeographic structure between the northern (L. l. annectens) and the two-southern subspecies (L. l. longicaudis and L. l. enudris), and although shallower, we also detected genetic differentiation between the two South American subspecies. Such genetic differentiation corresponds to the hydrography and to the geographical barriers characteristic of the distributional area of the species. We found a correlation between the shape of the skull and mandible with the environmental variation through the distribution of the species, and we rejected the hypothesis of niche equivalency and similarity between the three identified genetic lineages, suggesting adaptations to different environmental conditions. Our results support that the variation in environmental conditions, in concert with geographical barriers to gene flow and hydrography, have led to population divergence of L. longicaudis along the Neotropics. These results have important taxonomic implications for the species and its conservation.     

Genomic clues to the parental origin of the wild flowering cherry <Emphasis Type="Italic">Prunus yedoensis</Emphasis> var. <Emphasis Type="Italic">nudiflora</Emphasis> (Rosaceae)

Prunus yedoensis Matsumura is one of the popular ornamental flowering cherry trees native to northeastern Asia, and its wild populations have only been found on Jeju Island, Korea. Previous studies suggested that wild P. yedoensis (P. yedoensis var. nudiflora) is a hybrid species; however, there is no solid evidence on its exact parental origin and genomic organization. In this study, we developed a total of 38 nuclear gene-based DNA markers that can be universally amplifiable in the Prunus species using 586 Prunus Conserved Orthologous Gene Set (Prunus COS). Using the Prunus COS markers, we investigated the genetic structure of wild P. yedoensis populations and evaluated the putative parental species of wild P. yedoensis. Population structure and phylogenetic analysis of 73 wild P. yedoensis accessions and 54 accessions of other Prunus species revealed that the wild P. yedoensis on Jeju Island is a natural homoploid hybrid. Sequence-level comparison of Prunus COS markers between species suggested that wild P. yedoensis might originate from a cross between maternal P. pendula f. ascendens and paternal P. jamasakura. Moreover, approximately 81% of the wild P. yedoensis accessions examined were likely F1 hybrids, whereas the remaining 19% were backcross hybrids resulting from additional asymmetric introgression of parental genotypes. These findings suggest that complex hybridization of the Prunus species on Jeju Island can produce a range of variable hybrid offspring. Overall, this study makes a significant contribution to address issues of the origin, nomenclature, and genetic relationship of ornamental P. yedoensis.     

The role of ecological factors in determining phylogeographic and population genetic structure of two sympatric island skinks (<Emphasis Type="Italic">Plestiodon kishinouyei</Emphasis> and <Emphasis Type="Italic">P. stimpsonii</Emphasis>)

We conducted comparative phylogeographic and population genetic analyses of Plestiodon kishinouyei and P. stimpsonii, two sympatric skinks endemic to islands in the southern Ryukyus, to explore different factors that have influenced population structure. Previous phylogenetic studies using partial mitochondrial DNA indicate similar divergence times from their respective closest relatives, suggesting that differences in population structure are driven by intrinsic attributes of either species rather than the common set of extrinsic factors that both presumably have been exposed to throughout their history. In this study, analysis of mtDNA sequences and microsatellite polymorphism demonstrate contrasting patterns of phylogeography and population structure: P. kishinouyei exhibits a lower genetic variability and lower genetic differentiation among islands than P. stimpsonii, consistent with recent population expansion. However, historical demographic analyses indicate that the relatively high genetic uniformity in P. kishinouyei is not attributable to recent expansion. We detected significant isolation-by-distance patterns among P. kishinouyei populations on the land bridge islands, but not among P. stimpsonii populations occurring on those same islands. Our results suggest that P. kishinouyei populations have maintained gene flows across islands until recently, probably via ephemeral Quaternary land bridges. The lower genetic variability in P. kishinouyei may also indicate smaller effective population sizes on average than that of P. stimpsonii. We interpret these differences as a consequence of ecological divergence between the two species, primarily in trophic level and habitat preference.     

Modern State of Populations of Endemic <Emphasis Type="Italic">Oxytropis</Emphasis> Species from Baikal Siberia and Their Phylogenetic Relationships Based on Chloroplast DNA Markers

The genetic diversity and population structure of the endemic species of Baikal Siberia Oxytropis triphylla, O. bargusinensis, and O. interposita were studied for the first time on the basis of the nucleotide polymorphism of intergenic spacers psbA–trnH, trnL–trnF, and trnS–trnG of chloroplast DNA. All populations of these species were characterized by a high haplotype (0.762–0.924) and relatively low nucleotide (0.0011–0.0022) diversity. Analysis of the distribution of variability in O. triphylla and O. bargusinensis showed that there was no significant genetic differentiation between populations of each species; the gene flow was 4.43 and 8.91, respectively. The high level of genetic diversity in the studied populations indicates a relatively stable state of these populations. A study of the phylogenetic relationships of closely related species confirms the concept of the origin of O. bargusinensis and O. tompudae as a result of intersectional hybridization of the species of the sections Orobia and Verticillares.     

Molecular genetic characteristic of dinucleotide microsatellite loci in parthenogenetic lizards <Emphasis Type="Italic">Darevskia unisexualis</Emphasis>

In the present study, the first molecular genetic investigation of dinucleotide (GT) _n microsatellite loci in parthenogenetic lizards Darevskia unisexualis was performed. New polymorphic locus, Du214, (GenBank Ac. No. EU252542) was identified and characterized in detail. It was demonstrated that allele of this locus differed in the size and structure of microsatellite locus, as well as in point mutations, the combinations of which enabled the isolation of stabile fixed double nucleotide substitutions A-A (alleles 2 and 4) and G-T (alleles 1, 3, 5, and 6). Double nucleotide substitutions described were also identified in the orthlogous loci of the parental species genomes, D. raddei (G-T) and D. valentine (A-A). Based on the analysis of allele distribution pattern at this locus in all populations of parthenospecies D. unisexualis, mathematic model was elaborated and realized. Using this model, frequencies of allelic variants for all populations of the species of interest were calculated and population genetic structure of D. unisexualis was characterized. Genetic contribution of each population to the species gene pool was determined. The data obtained demonstrated that microsatellite variation was one of the factors of clonal and genetic diversity of a parthenospecies.     

Phylogeographical structure inferred from cpDNA sequence variation of <Emphasis Type="Italic">Zygophyllum xanthoxylon</Emphasis> across north-west China

Zygophyllum xanthoxylon, a desert species, displaying a broad east–west continuous distribution pattern in arid Northwestern China, can be considered as a model species to investigate the biogeographical history of this region. We sequenced two chloroplast DNA spacers (psbK-psbI and rpl32-trnL) in 226 individuals from 31 populations to explore the phylogeographical structure. Median-joining network was constructed and analysis of AMOVA, SMOVA, neutrality tests and distribution analysis were used to examine genetic structure and potential range expansion. Using species distribution modeling, the geographical distribution of Z. xanthoxylon was modeled during the present and at the Last Glacial Maximum (LGM). Among 26 haplotypes, one was widely distributed, but most was restricted to either the eastern or western region. The populations with the highest levels of haplotype diversity were found in the Tianshan Mountains and its surroundings in the west, and the Helan Mountains and Alxa Plateau in the east. AMOVA and SAMOVA showed that over all populations, the species lacks phylogeographical structure, which is speculated to be the result of its specific biology. Neutrality tests and mismatch distribution analysis support past range expansions of the species. Comparing the current distribution to those cold and dry conditions in LGM, Z. xanthoxylon had a shrunken and more fragmented range during LGM. Based on the evidences from phylogeographical patterns, distribution of genetic variability, and paleodistribution modeling, Z. xanthoxylon is speculated most likely to have originated from the east and migrated westward via the Hexi Corridor.     

Demographic history and population genetic structure of <Emphasis Type="Italic">Hagenia abyssinica</Emphasis> (Rosaceae), a tropical tree endemic to the Ethiopian highlands and eastern African mountains

The distribution of genetic diversity among natural populations is significantly shaped by geographical and environmental heterogeneity. The key objectives of this study were to outline the population genetic structure and to investigate the effects of historical and current factors in shaping the population structure of an endemic tropical tree, Hagenia abyssinica. We used 11 polymorphic microsatellites to estimate genetic variability and evaluate gene flow among natural populations of H. abyssinica. Further, we employed ecological niche modeling approaches, to analyze the demographic history and map potential distributions of H. abyssinica during the Last Glacial Maximum and the present. Significant levels of genetic diversity (H _O = 0.477, H _E = 0.439) were observed among the sampled locations. High coefficient of genetic differentiation (F _ST = 0.32) and considerable genetic variation within the sampled locations (68.01%) were detected. Our results indicated the existence of three genetic groups with limited gene exchange and revealed positive correlations (r = 0.425, P < 0.05) between genetic diversity and geo-graphic distance. The ecological niche modeling (ENM) results support the existence of three distribution zones during the Last Glacial Maximum (LGM), with high probability of occurrence (0.8–1.0), and indicated slight distribution disturbances during and after the LGM. The fundamental patterns of genetic diversity and population structuring of H. abyssinica result from a combination of both environmental and geographical factors, including long-term isolation by distance and characteristic life history of this species. Our ENM results identified three zones that could have served as glacial refugia for this species and lay a foundation for further studies, outlining demographic histories and population structures of Afromontane species.     

Genetic diversity and population structure of the roughskin sculpin (<Emphasis Type="Italic">Trachidermus fasciatus</Emphasis> Heckel) inferred from microsatellite analyses: implications for its conservation and management

Identification of population units is crucial for management and monitoring programs, especially for endangered wild species. The roughskin sculpin (Trachidermus fasciatus Heckel) is a small catadromous fish and has been listed as a second class state protected aquatic animal since 1988 in China. To achieve sustainable conservation of this species, it is necessary to clarify the existing genetic structure both between and within populations. Here, population genetic structure among eight populations of T. fasciatus were analyzed by using 16 highly polymorphic microsatellites. High levels of genetic variation were observed in all populations. All pairwise F _ST estimates were significant after false discovery rate correction (overall average F _ST = 0.054). Furthermore, both STRUCTURE and discriminant analysis of principal components (DAPC) analysis showed that the eight populations were grouped into six clusters. BAYESASS analysis showed generally low recent and asymmetric migration among populations. All these results suggested significant genetic structure across populations. However, there was no isolation by distance relationship among populations, likely resulting from barriers to gene flow created by habitat fragmentation. Our results highlight the need for in situ conservation efforts for T. fasciatus across its entire distribution range, through maximizing habitat size and quality to preserve overall genetic diversity and evolutionary potential.     

Genetic Diversity and Phylogenetic Relationships of Two Closely Related Northeast China <Emphasis Type="Italic">Vicia</Emphasis> Species Revealed with RAPD and ISSR Markers

RAPD and ISSR analyses revealed genetic diversity and relationships among 11 populations of two closely related northeast China Vicia species, Vicia ramuliflora and V. unijuga. Both methods yielded similar and complementary results, showing high genetic diversity. Vicia ramuliflora had 100% polymorphic loci in both RAPD and ISSR, and V. unijuga had 100% polymorphic loci for RAPD and 98.96% for ISSR. Genetic differentiation was moderate among populations of each species. Genetic variation was distributed mainly within populations for the two species. The high level of gene flow was important for the allocation of genetic variation. The UPGMA dendrogram and principal coordinates analysis at the level of individuals and populations showed that V. ramuliflora and V. unijuga were more closely related than either of them was to the outgroup species, V. cracca. The small molecular variance of V. ramuliflora and V. unijuga supports the conclusion that these two species had a common ancestor.     

Geography,geology and ecology influence population genetic diversity and structure in the endangered endemic Azorean <Emphasis Type="Italic">Ammi</Emphasis> (Apiaceae)

Three Azorean endemic Ammi species were initially described: Ammi trifoliatum (Wats.) Trel., Ammi seubertianum (Wats.) Trel. and Ammi huntii (Wats.) Trel. Many taxonomic changes have been conducted, and one to three species have been considered. Two species are currently accepted: A. trifoliatum, which occurs in almost every island, and A. seubertianum, with a narrower distribution. In this research, the population genetic diversity and structure of the Azorean Ammi species were assessed using eight specifically designed SSR markers. A wide sampling of A. seubertianum and A. trifoliatum was conducted in seven Azorean islands, and four A. huntii herbarium samples were also included to further contribute to the taxonomy of this genus in Azores. Flores populations showed the highest genetic diversity, while North of Topo, in São Jorge, showed the lowest. None of the populations analysed displayed signs of putative inbreeding. The population genetic structure analyses conducted partially provided support for the two currently accepted species, but other possible cryptic taxa may also be present. The complex clustering obtained seems to result from a combined action of geography, geology and ecology, and although island-specific genetic patterns were found, environmental conditions connected to different altitudes and the existence of micro-niches may also play an important role. A thorough morphological revision and ecophysiological studies should be conducted to clarify the number of endemic taxonomic entities present in the Azores archipelago.     

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

Mesoscale investigations based on microsatellite analysis of the freshwater sponge <Emphasis Type="Italic">Ephydatia fluviatilis</Emphasis> in the River-Sieg system (Germany) reveal a genetic divergence

Freshwater sponges play a major role in freshwater ecological system as important filter-feeding organisms and bioindicators. There are only few data about their ecological diversity and population genetic structure available, though a deeper knowledge is needed to propose proper conservation and effective management. The aim of this study was to assess data on distribution patterns of freshwater sponges to study the connectivity of genotypes of Ephydatia fluviatilis in a river system. We sampled specimens from River-Sieg system (River Agger and River Sieg, Germany). We hypothesized that strong anthropogenic influence would cause a uniform distribution of population structures. The genetic structure of E. fluviatilis populations was analysed with a set of eleven microsatellite loci from seven locations in River-Sieg system. Besides of E. fluviatilis, three other species co-occurred (Ephydatia mülleri, Spongilla lacustris, Eunapius fragilis). In contrast to our hypothesis, we observed an overall correlation between genetic and geographic distances among populations of this sessile species, which follows a clear isolation-by-distance pattern. A significant microsatellite polymorphism and high levels of genetic divergence between populations (F_ST) in upstream reaches were present. These results will provide important information for conservation management of populations with limited dispersal ability in connected river systems.     

Genetic diversity and population structure of <Emphasis Type="Italic">Melia azedarach</Emphasis> in North-Western Plains of India

Key message

Genetic structure among M. azedarach populations was detected and two subpopulations were present among them. A significant ‘isolation by distance’ was found in M. azedarach population in North-Western Plains of India.

Abstract

Melia azedarach is an important forest tree with pharmaceutical, insecticidal, pesticidal, and commercial significance. It is a good reforestation tree because of its fast growth and drought hardy nature. Genetic variation in a species allows itself to adapt, evolve and respond to environmental stress. It provides the basis for survival of a species and critically influences its evolutionary potential. Assessment of genetic diversity is necessary for improvement and conservation of a species. For this, microsatellite markers are of particular interest given the attributes like co-dominance, reproducibility, hyper variability and abundance throughout the genome. In the present study, we analyzed the genetic diversity and population structure of M. azedarach, an ecologically imperative species growing in the North-Western Plains of India. We developed 43 microsatellite markers, of which 20 were subsequently employed for analysis of diversity and population structure among 33 populations encompassing 318 genotypes representing North-Western Plains of India. A moderate level of diversity (Na = 5.1, Ho = 0.506, He = 0.712, I = 1.386) was assessed. The highest value of ΔK estimated using STRUCTURE indicated 2 subpopulations (K = 2). AMOVA exhibited 73 % variation within populations and 12 % variation was found among regions. Significant positive correlation between geographical and genetic distance was found (Rxy = 0.365, P = 0.010). The present study lays a foundation on a better understanding of genetic dynamics of the species and reveals its diversity and population structure in North-Western Plains of India.