Isolation and characterization of microsatellite loci in the highland papaya Vasconcellea × heilbornii V. Badillo (Caricaceae)

Nine polymorphic microsatellite loci were identified in the tropical plant Vasconcellea ×heilbornii and used to estimate allelic diversity in two populations of southern Ecuador. Allelic richness ranged from two to five alleles, observed heterozygosity ranged from 0.150 to 0.947 and expected heterozygosity ranged from 0.186 to 0.701. Most of these markers also amplified microsatellite loci in two other Vasconcellea species (Vasconcellea stipulata and Vasconcellea cundinamarcensis). Hence, these markers will be useful for population genetic analysis and the evaluation of genetic diversity and gene flow in these species.     

Conservation genetics of Butte County meadowfoam (Limnanthes floccosa ssp. californica Arroyo), an endangered vernal pool endemic

The endangered annual plant Limnanthes floccosa ssp. californica Arroyo is restricted to vernal pools in Butte County, California. To identify populations with unique genetic resources, guide reintroduction efforts, and design seed collection scenarios for long-term ex situ seed storage we determined extant genetic diversity and structure by surveying 457 individuals from 21 distinct populations using nine polymorphic microsatellite markers. We found low within population genetic diversity: low allelic diversity (1.9 [0.06 SE] alleles/locus); low heterozygosity (H _obs = 0.10 ± 0.018, H _exp = 0.19 ± 0.015), and a high fixation index (0.556 ± 0.044). The number of polymorphic loci ranged between 11 and 89%. Bayesian ordination determined 20 distinct populations and we found high genetic structure among these (F _st = 0.65, P < 0.0001). We identified notable gene flow barriers across populations, confirming regional structuring between three previously defined population density centers and two outlying populations (F _st = 0.21, P < 0.0001). Population size estimates ranged between ~50 and >5,000 extant plants per site. Our study confirms previous isozyme-based results and suggests that the loss of any population would represent a significant loss in the species’ genetic diversity. Recovery requires active restoration of existing populations and permanent habitat protection. We recommend close comparison of microhabitats of declining populations with genetically similar populations, to determine the potential for human assisted gene flow via seed movement to recover declining populations.     

Polyploid origin,genetic diversity and population structure in the tetraploid sea lavender <Emphasis Type="Italic">Limonium narbonense</Emphasis> Miller (Plumbaginaceae) from eastern Spain

Limonium narbonense Miller is a fertile tetraploid species with a sporophytic self-incompatibility system. This sea lavender is found in coastal salt marshes which have been under intense human pressure during the past decades resulting in significant habitat fragmentation. Eleven microsatellite loci specifically designed for this species were amplified in 135 individuals from five populations. These markers were used to investigate the polyploid nature, the levels of genetic diversity and population structure in this species. L. narbonense showed high levels of genetic diversity (A = 7.82, P = 100% H _T = 0.446), consistent with its likely autotetraploid origin revealed in this study and obligate outcrossing breeding system. Inbreeding (F _IS) values were low in the three southern populations (mean F _IS = 0.062), and higher in the northern populations (mean F _IS = 0.184). Bayesian analysis of population structure revealed that populations could be grouped into two genetic clusters, one including three southern populations and the other the two northernmost ones. Individuals from the two northernmost populations showed higher admixture of the two genetic clusters than individuals from the three southern ones. A thorough analysis of microsatellite electrophoretic patterns suggests an autotetraploid origin for L. narbonense. The genetic structure revealed in this study is attributed to a recent migration from the southern area. This result suggests a net gene flow from the south to the north, likely facilitated by migratory movements of birds visiting the temporary flooded ponds occupied by L. narbonense.     

Genetic divergence in nuclear genomes between populations of <Emphasis Type="Italic">Fagus crenata</Emphasis> along the Japan Sea and Pacific sides of Japan

Genetic diversity and structure in Fagus crenata were studied by analyzing 14 nuclear microsatellite loci in 23 populations distributed throughout the species’ range. Although population differentiation was very low (F _ST = 0.027; R _ST = 0.041), both neighbor-joining tree and Bayesian clustering analyses provided clear evidence of genetic divergence between populations along the Japan Sea (Japan Sea lineage) and Pacific (Pacific lineage) sides of Japan, indicating that physical barriers to migration and gene flow, notably the mountain ranges separating the populations along the Japan Sea and Pacific sides, have promoted genetic divergence between these populations. The two lineages of the nuclear genome are generally consistent with those of the chloroplast genome detected in a previous study, with several discrepancies between the two genomes. Within-population genetic diversity was generally very high (average H _E = 0.839), but decreased in a clinal fashion from southwest to northeast, largely among populations of the Japan Sea lineage. This geographical gradient may have resulted from the late-glacial and postglacial recolonization to the northeast, which led to a loss of within-population genetic diversity due to cumulative founder effects.     

Genetic structure of the critically endangered Red‐headed Wood Pigeon Columba janthina nitens and its implications for the management of threatened island populations

The Red‐headed Wood Pigeon Columba janthina nitens is endemic to the Ogasawara Islands, an oceanic island chain located 1000 km south of the main islands of Japan. The subspecies is at high risk of extinction because of its small population size and restricted habitat range. We undertook genetic analyses of this pigeon using sequences of a portion of the mitochondrial control region and five microsatellite markers to estimate the genetic characteristics of two wild populations from the Bonin and Volcano Islands, as well as one captive breeding population. The genetic diversity of the wild individuals was exceptionally low in both the mitochondria (nucleotide diversity = 0.00105) and at the microsatellite (3.2 alleles per locus and H_E = 0.12) loci. Higher numbers of microsatellite genotypes were observed in the Volcano Islands population than in the Bonin Islands population, which may be because of the relatively low impact of human disturbance. The most common mitochondrial haplotypes and microsatellite alleles observed in the two wild populations were completely fixed in the captive population. Our results suggest that the genetic diversity of the captive population needs to be increased. However, introduction of a wild individual into a captive population can lead to a decreased genetic diversity in the wild population and therefore should be done with caution. The genetic differentiation between the Bonin and the Volcano island groups was low, and the populations of the two island groups should be regarded as a single evolutionarily significant unit. However, special consideration is required for habitat conservation in the Volcano Islands, which may be functioning as a sanctuary for the Red‐headed Wood Pigeon. For the long‐term conservation of threatened bird species that live on remote oceanic islands, determination of management units considering gene flow caused by their flying capacity and maintenance of genetically suitable wild and captive populations are essential.     

Microsatellite analysis of the natterjack toad (Bufo calamita) in Denmark: populations are islands in a fragmented landscape

The European natterjack toad (Bufo calamita) has declined rapidly in recent years, primarily due to loss of habitat, and in Denmark it is estimated that 50% of the isolated populations are lost each decade. To efficiently manage and conserve this species and its genetic diversity, knowledge of the genetic structure is crucial. Based on nine polymorphic microsatellite loci, the genetic diversity, genetic structure and gene flow were investigated at 12 sites representing 5–10% of the natterjack toad localities presently known in Denmark. The expected heterozygosity (H _E) within each locality was generally low (range: 0.18–0.43). Further analyses failed to significantly correlate genetic diversity with population size, degree of isolation and increasing northern latitude, indicating a more complex combination of factors in determining the present genetic profile. Genetic differentiation was high (overall θ = 0.29) and analyses based on a Bayesian clustering method revealed that the dataset constituted 11 genetic clusters, defining nearly all sampling sites as distinct populations. Contemporary gene flow among populations was undetectable in nearly all cases, and the failure to detect a pattern of isolation by distance within major regions supported this apparent lack of a gene flow continuum. Indications of a genetic bottleneck were found in three populations. The analyses suggest that the remaining Bufo calamita populations in Denmark are genetically isolated, and represent independent units in a highly fragmented gene pool. Future conservation management of this species is discussed in light of these results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Analysis of genetic diversity in the endangered Hucho taimen from China

Hucho taimen are listed as endangered in China. The population size has declined recently, prompting an increase in the level of listing from grade three in 2002 to grade five in 2006. We analyzed the genetic diversity of wild populations using 17 microsatellite markers to establish a scientific basis for conservation of this species. We collected tissue samples from four populations in the Heilongjiang River basin: Huma River (HM), Hutou (HT), Haiqing (HQ), and Zhuaji (ZJ). A total of 21 loci were amplified, 18 of which were polymorphic. The number of alleles per locus ranged from 2 to 9 (mean: 4.1905). There were 13 highly polymorphic loci and 5 moderately polymorphic loci. Analysis of five genetic diversity parameters (Na, Ne, Ho, He, and PIC) suggested moderate levels of diversity within the populations. The populations were ranked HT > HQ > ZJ > HM, but the differences in diversity were not statistically significant (P > 0.05). A comparison of variation among all four populations suggested Hardy–Weinberg disequilibrium at 20% of the loci. Genetic differentiation (Fst) was 0.0644 and the gene flow among populations was estimated at 3.36 individuals per generation. The majority of diversity (93.88%) occurred among individuals within a population. In contrast, relatively little (6.12%) of the genetic diversity was distributed between the populations. An analysis of genetic differentiation and genetic distance between pairs of populations revealed that both parameters were higher in comparisons of the HM population to the HT, HQ, and ZJ populations than among the three latter populations. This suggests that the HM population has a distinct genetic structure. We hypothesize that habitat degradation and excessive fishing, not low genetic diversity, has caused the decline in H. taimen populations. However, this species should be protected from further declines in genetic diversity.     

Conservation genetics of the alligator snapping turtle: cytonuclear evidence of range-wide bottleneck effects and unusually pronounced geographic structure

A previous mtDNA study indicated that female-mediated gene flow was extremely rare among alligator snapping turtle populations in different drainages of the Gulf of Mexico. In this study, we used variation at seven microsatellite DNA loci to assess the possibility of male-mediated gene flow, we augmented the mtDNA survey with additional sampling of the large Mississippi River System, and we evaluated the hypothesis that the consistently low within-population mtDNA diversity reflects past population bottlenecks. The results show that dispersal between drainages of the Gulf of Mexico is rare (F _STmsat  = 0.43, Φ_STmtDNA = 0.98). Past range-wide bottlenecks are indicated by several genetic signals, including low diversity for microsatellites (1.1–3.9 alleles/locus; H _e = 0.06–0.53) and mtDNA (h = 0.00 for most drainages; π = 0.000–0.001). Microsatellite data reinforce the conclusion from mtDNA that the Suwannee River population might eventually be recognized as a distinct taxonomic unit. It was the only population showing fixation or near fixation for otherwise rare microsatellite alleles. Six evolutionarily significant units are recommended on the basis of reciprocal mtDNA monophyly and high levels of microsatellite DNA divergence.     

Genetic diversity in fragmented populations of the critically endangered spider orchid Caladenia huegelii: implications for conservation

The Orchidaceae is characterised by a diverse range of life histories, reproductive strategies and geographic distribution, reflected in a variety of patterns in the population genetic structure of different species. In this study, the genetic diversity and structure was assessed within and among remnant populations of the critically endangered sexually deceptive orchid, Caladenia huegelii. This species has experienced severe recent habitat loss in a landscape marked by ancient patterns of population fragmentation within the Southwest Australian Floristic Region, a global biodiversity hotspot. Using seven polymorphic microsatellite loci, high levels of within-population diversity (mean alleles/locus = 6.73; mean H _E = 0.690), weak genetic structuring among 13 remnant populations (F _ST = 0.047) and a consistent deficit of heterozygotes from Hardy–Weinberg expectation were found across all populations (mean F _IS = 0.22). Positive inbreeding coefficients are most likely due to Wahlund effects and/or inbreeding effects from highly correlated paternity and typically low fruit set. Indirect estimates of gene flow (Nm = 5.09 using F _ST; Nm = 3.12 using the private alleles method) among populations reflects a historical capacity for gene flow through long distance pollen dispersal by sexually deceived wasp pollinators and/or long range dispersal of dust-like orchid seed. However, current levels of gene flow may be impacted by habitat destruction, fragmentation and reduced population size. A genetically divergent population was identified, which should be a high priority for conservation managers. Very weak genetic differentiation indicates that the movement and mixing of seeds from different populations for reintroduction programs should result in minimal negative genetic effects.     

Development of Microsatellite Markers and Analysis of Genetic Diversity and Population Structure of Colletotrichum gloeosporioides from Ethiopia

Twenty three polymorphic microsatellite markers were developed for citrus plant pathogenic fungus, Colletotrichum gloeosporioides, and were used to analyze genetic diversity and population structure of 163 isolates from four different geographical regions of Ethiopia. These loci produced a total of 118 alleles with an average of 5.13 alleles per microsatellite marker. The polymorphic information content values ranged from 0.104 to 0.597 with an average of 0.371. The average observed heterozygosity across all loci varied from 0.046 to 0.058. The gene diversity among the loci ranged from 0.106 to 0.664. Unweighted Neighbor-joining and population structure analysis grouped these 163 isolates into three major groups. The clusters were not according to the geographic origin of the isolates. Analysis of molecular variance showed 85% of the total variation within populations and only 5% among populations. There was low genetic differentiation in the total populations (F_ST = 0.049) as evidenced by high level of gene flow estimate (N_m = 4.8 per generation) among populations. The results show that Ethiopian C. gloeosporioides populations are generally characterized by a low level of genetic diversity. The newly developed microsatellite markers were useful in analyzing the genetic diversity and population structure of the C. gloeosporioides populations. Information obtained from this study could be useful as a base to design strategies for better management of leaf and fruit spot disease of citrus in Ethiopia.     

Genetic impacts of habitat loss on the rare ironstone endemic Tetratheca paynterae subsp. paynterae

Expansion of mining in the banded ironstone ranges of southern Western Australia has focussed attention on the genetic impacts of habitat loss on rare endemic taxa. One example is Tetratheca paynterae subsp. paynterae (Elaeocarpaceae), an insect-pollinated, perennial shrub confined to 4 ha of banded ironstone outcrops in the Windarling Range. Mining removed 1,900 of the 7,700 recorded plants in 2004. Further reductions could occur if it can be demonstrated that the viability of the remaining population is not threatened. To investigate the potential impact of reductions in population size due to mining we first used Bayesian clustering and principal coordinate analysis to define population boundaries based on differentiated gene pools. The level of genetic diversity and spatial genetic structuring was then compared among populations that ranged in size from 46 to 4400 individuals. Analysis with 11 microsatellite loci revealed lower genetic diversity in small populations (A _R = 4.5–4.8) than a large population (A _R = 6.3) together with significant pair-wise differences among populations separated by distances of 80 m or more. Spatial autocorrelation analysis showed the extent of spatial genetic structure differed among populations of different size, consistent with near-neighbour mating and limited dispersal. Fine scale spatial structure was consistent with historically low gene flow. Analysis of the impact of possible expansions in mining revealed small, isolated populations of T. paynterae subsp. paynterae were of high conservation value. While their removal would reduce plant numbers and genetic diversity by less than 5%, unique genotypes will be lost resulting in a 30% decline in genetic differentiation.     

Genetic diversity and population differentiation of the Chinese soft-shelled turtle (Pelodiscus sinensis) in three geographical populations

The Chinese soft-shelled turtle (Pelodiscus sinensis) is one of the most important economical chelonians in the world. To understand the genetic variations of the Chinese soft-shelled turtle in China, 62 individuals were sampled from three localities and 18 polymorphic microsatellite loci tested were used to detect genetic diversity and population structure. Results showed that the genetic diversity of the wild P. sinensis was high. Except for the Wuhu populations, the majority of microsatellite loci are not deviation from Hardy–Weinberg equilibrium in the other two populations. AMOVA analysis indicated that genetic variations occurred mainly within populations (97.4%) rather than among populations (2.6%). The gene flow estimates (Nm) among three geographic populations demonstrated that strong gene flow existed (Nm > 1, mean 6). The present study supported that different habitats, breed turtles escaped, multiple paternity and long evolutionary history may be responsible for the current genetic diversity and differentiation in the wild Chinese soft-shelled turtle.     

Microsatellite multiplexes for high‐throughput genotyping of French grunts (Haemulon flavolineatum,Pisces: Haemulidae) and their utility in other grunt species

Grunts of the genus Haemulon are abundant in the western Atlantic and are an important component of regional fisheries. Management efforts for these and other Caribbean reef fish require an understanding of population structure and the extent of dispersal among populations. We characterized 11 polymorphic microsatellite loci for French grunts (Haemulon flavolineatum) and grouped them into three PCR multiplexes to facilitate high‐throughput genotyping. These loci will be useful for population genetic studies of French grunts. Cross‐species amplifications suggest these loci are also suitable for population genetic studies of H. aurolineatum, H. plumieri, and H. sciurus.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Genetic variation of wild litchi (<Emphasis Type="Italic">Litchi chinensis</Emphasis> Sonn. subsp. <Emphasis Type="Italic">chinensis</Emphasis>) revealed by microsatellites

Understanding the amount and distribution of genetic diversity in natural populations can inform the conservation strategy for the species in question. In this study, genetic variation at eight nuclear microsatellite loci was used to investigate genetic diversity and population structure of wild litchi (Litchi chinensis Sonn. subsp. chinensis). Totally 215 individuals were sampled, representing nine populations of wild litchi. All eight loci were polymorphic, with a total of 51 alleles. The expected heterozygosity in the nine populations ranged from 0.367 to 0.638 with an average value of 0.526. Inbreeding within wild litchi populations was indicated by a strong heterozygote defect. Significant bottleneck events were detected in the populations from Yunnan and Vietnam, which could be responsible for lower levels of genetic diversity in these populations. Measures of genetic differentiation (F _ST = 0.269) indicated strong differentiation among wild litchi populations. Significant correlation was found between genetic differentiation and geographical distance (r = 0.655, P = 0.002), indicating a strong isolation by distance in these populations. Bayesian clustering suggested genetic separation among three regional groups, namely, the western group, the central group and the eastern group. Some conservation strategies for wild litchi populations were also proposed based on our results.     

Demographic bottlenecks and low gene flow in remnant populations of the critically endangered <Emphasis Type="Italic">Berchemiella wilsonii</Emphasis> var<Emphasis Type="Italic">. pubipetiolata</Emphasis> (Rhamnaceae) inferred from microsatellite markers

Berchemiella wilsonii var. pubipetiolata (Rhamnaceae) is an endangered plant with only four remnant populations in eastern China. Population genetic information is essential for understanding population history and formulating conservation strategies for this species. Thirteen microsatellite loci were used to investigate genetic variation and population structure of the four remnant populations. Moderate levels of expected heterozygosity (H _E = 0.466–0.543) and low allelic diversity (A = 3.1–3.6 and A _R = 2.2–2.4, respectively) were observed within populations. Bottleneck tests found three out of four populations to deviate from mutation-drift equilibrium under the two-phase model (TPM), suggesting a recent population decline, which is congruent with known demographic history. The evolutionary history of the species seems dominated by genetic drift rather than gene flow. Low historical gene flow was inferred from several different approaches and N _m ranged from 0.582 by the private allele method to 0.783 by the coalescent method. Contemporary gene flow was also found to be even lower for only one first generation migrant was detected with individual-based assignment analysis. Restricted pollen and seed dispersal as well as a recent decline in population size associated with habitat fragmentation may have contributed to low levels of historical and contemporary gene flow, and resulted in a high genetic differentiation. Under this scenario, Berchemiella wilsonii var. pubipetiolata populations are expected to display more pronounced population genetic structure in the future as a result of increased inbreeding and genetic drift.     

Genetic variation within and between populations of hermaphroditic <Emphasis Type="Italic">Bulinus truncatus</Emphasis> tetraploid freshwater snails of the Albertine Rift,East Africa

Genetic variability within and among Bulinus truncatus of the Albertine Rift freshwater bodies were assessed to investigate the degree of inbreeding and gene flow in the snail populations. The effect of ploidy on the genetic structuring of B. truncatus is also described. We characterized the genetic structure of seven B. truncatus populations from Lake Albert, Lake Kivu, and Katosho swamp in Tanzania using five polymorphic microsatellite loci. Genetic differentiation was quantified using pairwise FST values and Nei’s standard genetic distances. Different alleles were observed across all loci and genetic diversity was low although it varied greatly across populations; observed heterozygosity was, however, higher than the expected heterozygosity in three of the populations studied. Significant heterozygote deficiencies were observed coupled with significant linkage disequilibria in five populations for all the five loci examined in this study. We found significant genetic differentiation among the seven freshwater bodies; private alleles were observed across all loci indicating restricted or absence of gene flow between populations. Limited snail dispersal and the reproductive biology of B. truncatus are the major forces shaping the genetic variation observed. Low genetic variation within B. truncatus populations exposes them to a high parasite infection risk as predicted in the Red Queen hypothesis.     

Genetic structure of the threatened Phaedranassa schizantha (Amaryllidaceae)

Phaedranassa schizantha (Amaryllidaceae) is an endangered species endemic to Ecuador and two varieties have been described: P. schizantha var. schizantha and P. schizantha var. ignea. We assessed population genetic structure and demographic patterns in 11 populations across the range of the species using 13 microsatellite loci. Our data show that genetic diversity was generally lower in the southern part of the range and was especially low in populations closest to cities. We found significant population differentiation (F_ST = 0.14, D_EST = 0.34) and evidence of a genetic bottleneck. Genetic variation did not show isolation by distance. Instead, results suggest genetic barriers around two main cities. Bayesian analysis identified two genetic groups, neither of which represents either of the two varieties previously recognized. Coalescent analysis indicates a relatively recent colonization pattern between the two genetic groups (< 3000 generations). Conservation efforts need to be taken to facilitate genetic exchange between the groups, especially between locations that seem to be genetically isolated.     

Genetic diversity and population structure of three Indian horse breeds

The genetic relationships of three Indian horse breeds—Marwari, Spiti, and Kathiawari were studied by genotyping 96 individuals with 20 polymorphic microsatellite markers. A total of 157 alleles were detected across 20 polymorphic loci. The Marwari population showed the highest allelic diversity (A = 5.7 and Ar = 5.14), followed by Spiti (A = 4.9 and Ar = 4.74) and Kathiawari (A = 4.1 and Ar = 3.82). The gene diversity was highest in the Spiti population (He = 0.67), followed by Marwari (He = 0.66) and Kathiawari (He = 0.59). Within population inbreeding estimates (f) in Marwari, Spiti and Kathiawari breeds were 0.18, 0.08, and 0.07, respectively, suggesting high level of inbreeding in these breeds. Analysis of bottleneck revealed evidence of recent bottleneck in Spiti and Kathiawari populations. Pair-wise Fst analysis, AMOVA and assignment tests demonstrated high genetic differentiation and low gene flow between populations. The information about genetic diversity and population structure will be useful for the future development of effective breeding management in order to preserve these Indian horse breeds.