Genomic diversity and affinities in population groups of North West India: An analysis of Alu insertion and a single nucleotide polymorphism

The North West region of India is extremely important to understand the peopling of India, as it acted as a corridor to the foreign invaders from Eurasia and Central Asia. A series of these invasions along with multiple migrations led to intermixture of variable populations, strongly contributing to genetic variations. The present investigation was designed to explore the genetic diversities and affinities among the five major ethnic groups from North West India; Brahmin, Jat Sikh, Bania, Rajput and Gujjar. A total of 327 individuals of the abovementioned ethnic groups were analyzed for 4 Alu insertion marker loci (ACE, PV92, APO and D1) and a Single Nucleotide Polymorphism (SNP) rs2234693 in the intronic region of the ESR1 gene. Statistical analysis was performed to interpret the genetic structure and diversity of the population groups. Genotypes for ACE, APO, ESR1 and PV92 loci were found to be in Hardy–Weinberg equilibrium in all the ethnic groups, while significant departures were observed at the D1 locus in every investigated population after Bonferroni's correction. The average heterozygosity for all the loci in these ethnic groups was fairly substantial ranging from 0.3927 ± 0.1877 to 0.4333 ± 0.1416. Inbreeding coefficient indicated an overall 10% decrease in heterozygosity in these North West Indian populations. The gene differentiation among the populations was observed to be of the order of 0.013. Genetic distance estimates revealed that Gujjars were close to Banias and Jat Sikhs were close to Rajputs. Overall the study favored the recent division of the populations of North West India into largely endogamous groups. It was observed that the populations of North West India represent a more or less homogenous genetic entity, owing to their common ancestral history as well as geographical proximity.     

Genetic Polymorphisms in VIR Genes among Indian Plasmodium vivax Populations

The vir genes are antigenic genes and are considered to be possible vaccine targets. Since India is highly endemic to Plasmodium vivax, we sequenced 5 different vir genes and investigated DNA sequence variations in 93 single-clonal P. vivax isolates. High variability was observed in all the 5 vir genes; the vir 1/9 gene was highly diverged across Indian populations. The patterns of genetic diversity do not follow geographical locations, as geographically distant populations were found to be genetically similar. The results in general present complex genetic diversity patterns in India, requiring further in-depth population genetic and functional studies.     

The origin of Indian Star tortoises (Geochelone elegans) based on nuclear and mitochondrial DNA analysis: A story of rescue and repatriation

The Indian Star tortoise (Geochelone elegans) belongs to the family Testunidae and is distributed in southwest India and Sri Lanka. In addition to facing loss of its natural habitat, the species is also illegally traded as food and as an exotic pet internationally. Here we report DNA-based analyses for identification and repatriation of these tortoises into their natural habitat. We have attempted to establish the geographical origin of these tortoises rescued from smugglers, by comparing the microsatellite and mitochondrial markers of rescued animals with animals of known provenance. Star tortoises exhibited strong genetic structure in India. The populations from western India were genetically distinct at microsatellite and mitochondrial loci from southern populations. The rescued individuals had similar multilocus genotypes and mitochondrial DNA haplotypes as the reference individuals from south India. However, the precise geographic origin of many of the rescued samples remains unresolved, because we could not assign them to southern populations and the Neighbor-Joining cluster analysis indicated that some of rescued tortoises formed distinct clusters. These data strongly suggest that the rescued group of tortoises is composed of a mix of individuals from differentiated source populations that are probably located in southern India and possibly Sri Lanka. Our study provides valuable information based on molecular markers for the assessment of genetic diversity in Indian Star tortoises.     

Low genetic differentiation in a sedentary bird: house sparrow population genetics in a contiguous landscape

The house sparrow Passer domesticus has been declining in abundance in many localities, including Finland. We studied the genetic diversity and differentiation of the house sparrow populations across Finland in the 1980s, at the onset of the species'' decline in abundance. We genotyped 472 adult males (the less dispersive sex) from 13 locations in Finland (covering a range of 400 × 800 km) and one in Sweden (Stockholm) for 13 polymorphic microsatellite markers. Our analysis of Finnish ringing records showed that natal dispersal distances are limited (90% <16 km), which confirmed earlier finding from other countries. The Finnish populations were panmictic, and genetically very homogeneous and the limited dispersal was sufficiently large to maintain their connectivity. However, all Finnish populations differed significantly from the Stockholm population, even though direct geographical distance to it was often smaller than among Finnish populations. Hence, the open sea between Finland and Sweden appears to form a dispersal barrier for this species, whereas dispersal is much less constrained across the Finnish mainland (which lacks geographical barriers). Our findings provide a benchmark for conservation biologists and emphasize the influence of landscape structure on gene flow.     

Deciphering diversity in populations of various linguistic and ethnic affiliations of different geographical regions of India: Analysis based on 15 microsatellite markers

The extent of genetic polymorphism at fifteen autosomal microsatellite markers in 54 ethnically, linguistically and geographically diverse human populations of India was studied to decipher intrapopulation diversity. The parameters used to quantify intrapopulation diversity were average allele diversity, average heterozygosity, allele range (base pairs), and number of alleles. Multilocus genotype frequencies calculated for selected populations were utilized for testing conformity with the assumption of Hardy-Weinberg equilibrium. The exact test values, after Bonferroni correction, showed significant deviation amongst Gowda (vWA, Penta E); Dhangar, Satnami and Gounder (D8S1179); Hmar (FGA); Kuki and Balti (vWA) groups. Relatively low number of alleles and allelic diversity (base-pairs size) had been observed in populations of central India as compared with southern and northern regions of the country. The communities of Indo-Caucasoid ethnic origin and Indo-European linguistic family (Kshatriya of Uttar Pradesh) showed highest allelic diversity, as well as rare alleles, not reported in any other Indian populations. Analysis based on average heterozygosity was also found to be lowest among the populations of central India (0.729) and highest among the populations from north (0.777) and west (0.784) regions of the country, having Indo-Caucasoid ethnic origin and Austro-Asiatic linguistic affiliation. The maximum power of discrimination (85%-89%) had been observed at loci FGA, Penta E, D18S51 and D21S11, suggested high intrapopulation diversity in India. Genetic diversity revealed by STR markers was consistent with the known demographic histories of populations. Thus, the present study clearly demonstrated that the intrapopulation diversity is not only present at the national level, but also within smaller geographical regions of the country. This is the first attempt to understand the extent of diversity within populations of India at such a large scale at genomic level.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Population genetic structure of Penaeus monodon, in relation to monsoon current patterns in Southwest, East and Andaman coastal waters of India

The black tiger shrimp (Penaeus monodon), a commercially important penaeid species, is widely distributed across the Indo-Pacific region. Genetic diversity in P. monodon collected from eight geographical regions in Southwest, East and Andaman coastal waters of India (N = 418) was investigated using 10 polymorphic microsatellite loci. Average observed heterozygosity at sampled loci were high, ranging from 0.643 (Coromandel Coast) to 0.753 (South Andaman). Pairwise F_ST (ranged from 0.005 to 0.078) and R_ST (ranged from 0.005 to 0.171) estimates revealed surprisingly strong and statistically significant genetic structure among tiger shrimp populations. A synthetic map generated by multidimensional scaling shows an apparent cline in allele frequencies paralleling the roughly circular flow of surface currents in the Bay of Bengal. Significant heterozygote deficiencies were noted in most population samples at most loci. Andaman Island sites showed the highest diversity. Recognition of high genetic diversity and distinct population structuring of P. monodon in Indian seas has important implications for future domestication of this species in India, for two reasons: identification of the best wild founding stocks for aquaculture and, subsequently, the potential impacts of release of domesticates to the wild, either accidentally or deliberately (i.e. for stock enhancement).     

Haplotype analysis of genomic polymorphisms in and around the myotonic dystrophy locus in diverse populations of India

The frequencies of haplotypes based upon the (CTG)n repeat and three other biallelic markers in and around the myotonic dystrophy (DM) locus were estimated in 13 ethnically, linguistically and geographically diverse sub-populations of India. The range of CTG repeats in caste populations was 5-31, while in tribal populations the range was shorter (5-23). Extensive variation in frequencies of large (CTG)n alleles (> or =18 repeats) was found in Indian populations. The implications of this finding on DM epidemiology are discussed. Haplotype diversity was found to be very high in most populations. The majority of the Indian DM patients carried a haplotype that is commonly found among DM patients globally; this is the most common haplotype in the class of large (> or =18 repeats) CTG alleles. However, one haplotype was found to be present in particularly high frequency in Indian populations; this haplotype was also found among Indian DM patients. This haplotype may be a characteristic of Indian and possibly of other East Asian populations.     

Genetic variability and population structure in Sapindus emarginatus Vahl from India

Sapindus emarginatus is an economically important tropical tree species sparsely distributed in different geographical provinces like Gangetic Plains, Western Ghats, and Deccan Plateau in India. In the present paper estimation of genetic variability within and among 41 accessions representing five populations was carried out using 3 single primer amplification reaction (SPAR) methods viz. RAPD, DAMD and ISSR. The cumulative data analysis was carried out for all three SPAR methods, and showed 82.32% polymorphism across all the accessions of S. emarginatus. Jaccard's similarity values among 41 accessions ranged from 0.15 to 0.49 with an average value of 0.37. The intra-population genetic diversity revealed highest values of Nei's genetic diversity (0.19,) Shannon information index (0.29) and polymorphic loci (55.18%), among the accessions of Gujarat (GJ) population, while the corresponding lowest values were (0.10), (0.15) and (26.40%) respectively among the accessions of Rajasthan (RJ) population. The maximum inter-population average genetic distance (0.20) was between Karnataka (KA) and RJ, while the corresponding least genetic distance (0.06) was between Allahabad (AL) and Varanasi (VS) populations. The analysis of molecular variance (AMOVA) revealed maximum percentage of variation among individuals of populations (72%) followed by 16% among regions and 12% among populations. Principal coordinate analysis (PCA) of cumulative data also supported the clustering pattern in the UPGMA dendrogram. These results suggest that genetic diversity is corroborating with the geographical diversity. Mantel's test was performed which revealed a highly significant correlation between cumulative vs RAPD, and showed the maximum (0.93) correlation coefficient, followed by cumulative vs ISSR (0.78) and cumulative vs DAMD (0.91) respectively, and this clearly indicates that the SPAR methods (RAPD, DAMD and ISSR) are sufficiently informative and are suitable to analyze the genetic variability within and among the populations of S. emarginatus.     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.     

Biased morph ratios and skewed mating success contribute to loss of genetic diversity in the distylous Pulmonaria officinalis

Background and Aims

In heterostylous plant species, skewed morph ratios are not uncommon and may arise from a range of factors. Despite the recognized importance of skewed morph ratios on overall reproductive success within populations, little is known about the impact of skewed morph ratios on population genetic diversity and differentiation in heterostylous species. This study specifically aimed to clarify the effect of population size and morph bias on population genetic diversity and differentiation in the temperate forest herb Pulmonaria officinalis. This species is characterized by a distylous breeding system and shows morph-specific differences in reproductive success.

Methods

Genetic diversity was determined for 27 P. officinalis populations in northern Belgium by using eight recently developed microsatellite markers. Multiple regressions were used to assess the relationship between genetic diversity, morph bias and population size, and F_ST-values were calculated for short- and long-styled morphs separately to study genetic differentiation as a function of morph type.

Key Results

For all genetic measures used, morph bias was more important in explaining patterns of genetic diversity than population size, and in all cases patterns of population genetic diversity followed a quadratic function, which showed a symmetrical decrease in genetic diversity with increasing morph bias. However, probably due to the reproductive advantage of L-morphs relative to S-morphs, maximum genetic diversity was found in populations showing an excess of L-morphs (60·7 % L-morph). On the other hand, no significant difference in pairwise genetic distances between populations was observed between L- (0·107) and S-morphs (0·106).

Conclusions

Our results indicate that significant deviations from equal morph ratios not only affect plant reproductive success but also population genetic diversity of heterostylous plant species. Hence, when defining conservation measures for populations of heterostylous plant species, morph ratios should be considered as an important trait affecting their long-term population viability.     

Biodiversity patterns in Indian Ocean corals, and effects of taxonomic error in data

Corals from 26 sites in the Indian Ocean, from numerous taxonomic sources, are analysed for distribution patterns after applying a consistent synonymy. The data set contains double the amount of distribution information used previously. Cluster analysis shows several contiguous regional groupings within this ocean, and no geographically dispersed groupings exist. Between-site species similarity correlates negatively with distance between sites. Coral species diversity and genus diversity plotted with latitude and longitude confirm that a band of high diversity stretches across the entire Indian Ocean, but that in this Ocean, unlike the Pacific and greater Caribbean area, there is no clear gradient with latitude at least up to the latitudinal limits of coral growth. About half the species are widespread throughout the Indian Ocean. Species frequency/distribution curves and cumulative frequency curves show that 150 species occur at only 2–4 sites. Multi-dimensional scaling (MDS) analysis using only these species confirms that they are important in forming sub-regional groupings which are superimposed on a general Indian Ocean homogeneity.The perennial problem of error in taxonomic data sets is examined. Increasing taxonomic error is introduced into the data using random methods. It is found that about 25% more error can be added to this data set before the relationship between similarity coefficient and geographical distance between pairs of sites is lost. Measures extracted from the clustering procedure using the original data and the data sets with added taxonomic error, show a sharp loss in cluster formation after addition of about 10 or 20% more error.     

An updated checklist of the ants of India with their specific distributions in Indian states (Hymenoptera,Formicidae)

As one of the 17 megadiverse countries of the world and with four biodiversity hotspots represented in its borders, India is home to an impressive diversity of life forms. However, much work remains to document and catalogue the species of India and their geographic distributions, especially for diverse invertebrate groups. In the present study, a comprehensive and critical list of Indian ant species is provided with up-to-date state-wise distribution. A total of 828 valid species and subspecies names belonging to 100 genera are listed from India. Potential erroneous data, misidentifications and dubious distributional records that may exist in the literature are also identified. The present exhaustive listing of Indian ants will provide a holistic view about diversity and distribution and will also help to identify major undersampled areas where future sampling and taxonomic efforts should be directed.     

Genetic relationship among the musk shrews,Suncus murinus insectivora,inhabiting islands and the continent based on mitochondrial DNA types

Cleavage patterns of mitochondrial DNA (mtDNA) by restriction endonucleases were examined in musk shrews collected at six trapping sites on two Japanese and two Indonesian islands, on Sri Lanka, situated close to the Indian subcontinent, and on the mainland of East Bengal in Bangladesh. No variation of mtDNAs was found among the Japanese and Indonesian shrews, despite their geographical isolation by the sea. In contrast, at least six mtDNA types were present in the Sri Lanka and the Bangladesh populations (three types for each), and these two populations seemed to be differentiated to the extent, which could be compared to the mice-intersubspecific differences. These populations were also differentiated from the Japanese-Indonesian type. Furthermore, a similar level of differentiation was also estimated between two mtDNA types within these respective populations. This feral species might be considered unique because of its high emigration rate caused by human movements and its high rate of subspecific hybridization.     

Variability in abundance and fluxes of dimethyl sulphide in the Indian Ocean

Dimethyl sulphide (DMS) is a biogenic gas of climatic significance on which limited information is available from the Indian Ocean. To fill this gap, we collected data on DMS and total dimethylsulphoniopropionate (DMSP_t) by participating in a dozen cruises. Here, we discuss the variability in DMS and DMSP_t in the north and central Indian Ocean in terms of their spatial and temporal variation. DMS and DMSP_t exhibited significant spatial and temporal variability. Apart from the concentration gradients in DMS within the Arabian Sea, Bay of Bengal and Central Indian Ocean basins, differences in average abundances were conspicuous between these basins. The Arabian Sea contained more DMS (mixed layer average was 7.8 nM) followed by the Bay of Bengal (2.8 nM) and the Central Indian Ocean (2.7 nM). The highest concentrations of DMS and DMSP_t (525 nM and 916 nM, respectively) were found in upwelling regimes along the west coast of India during the Southwest monsoon and fall intermonsoon seasons. Average surface DMS was the highest in the Arabian Sea. On the other hand observed sea-to-air fluxes of DMS were higher in the Bay of Bengal due to the prevalence of turbulent conditions. In the Arabian Sea wind speeds were low and hence the sea-to-air fluxes. The total diffusive flux of DMS from the study area to atmosphere is estimated to be about 1.02 × 10¹² g S y⁻¹, which contributes to 4.1–6.3% of the global DMS emission     

Genetic diversity in inter‐simple sequence repeat profiles across natural populations of Indian pomegranate (Punica granatum L.)

Pomegranate (Punica granatum L.), in the monogeneric family Punicaceae, is found in Iran, Afghanistan, India and Mediterranean countries. Iran is considered to be its primary centre of origin. In India, pomegranate occurs naturally only in the Western Himalayan regions of Jammu and Kashmir, Himachal Pradesh and Uttarakhand States. However, there is no information about genetic variation in wild pomegranate at population level. In this paper, we describe genetic diversity across natural populations of Indian pomegranate based on inter‐simple sequence repeat (ISSR) markers. Forty‐nine accessions representing eight populations from two regions were analysed using ISSR. Seventeen ISSR primers resulted in 268 polymorphic bands, with 87.01% polymorphism throughout the accessions. Pair‐wise population genetic distances ranged from 0.05 to 0.45, with a mean of 0.25 between populations. amova and Nei’s genetic diversity analyses revealed higher genetic variation within populations than among populations. A higher genetic differentiation (G_ST) was observed between the spatially distant populations, indicating a low level of genetic exchange (Nm) among these populations. However, clustering of populations was not in accordance with their geographical affiliations in the tree. The results indicate that the ISSR method is sufficiently informative and powerful to assess genetic variability in pomegranate, and that patterns of genetic variability observed among populations of wild pomegranate from the Western Himalaya differ. Estimation of genetic variation reported here provides a significant insight for in situ conservation and exploitation of genetic resources for this economically important species as potential breeding material.     

Molecular characterization of Indian Sugarcane streak mosaic virus isolates reveals recombination and negative selection in the P1 gene

Sugarcane streak mosaic virus (SCSMV), a member of the genus Poacevirus is an important viral pathogen affecting sugarcane production in India. The P1 gene of ten Indian isolates was sequenced and compared with previously reported SCSMV isolates. Comparative sequence analysis revealed a high level of diversity in the P1 gene (83–98% nucleotide sequence identity; 87–100% amino acid sequence identity), and the Indian SCSMV isolates were found to be the most variable (up to 9% diversity at the amino acid level). Phylogenetic tree analysis showed clustering of 17 SCSMV isolates into two groups: group I included isolates from India (except SCSMV-TPT) and Pakistan, and group II consisted of isolates from Japan, Indonesia, Thailand and SCSMV-TPT. The results obtained from phylogenetic study were further supported by the different in silico analysis viz. SNPs (single nucleotide polymorphism), INDELs (insertion and deletion) and evolutionary distance analysis. A significant proportion of recombination sites were observed at the N terminal region of P1 gene. Analysis of selection pressure indicated that the P1 gene of the Indian SCSMV isolates is under strong negative or purifying selection. It is likely that recombination identified in Indian SCSMV isolates, along with strong purifying selection, enhances the speed of elimination of deleterious mutations in the P1 gene. The evolutionary processes (recombination and selection pressure) together contributed to the observed genetic diversity and population structure of Indian SCSMV isolates.     

Conservation genetics of endangered medicinal plant Commiphora wightii in Indian Thar Desert

To ascertain the conservation priorities and strategies for Commiphora wightii, an endangered medicinal plant of Indian Thar Desert, genetic diversity was estimated within and among different populations. The total of 155 amplification products were scored using ten each of RAPD and ISSR primers, exhibiting an overall 86.72% polymorphism across 45 individuals representing eight populations. The cumulative data of two markers were used to compute pair-wise distances. The Neighbor-Joining tree revealed high genetic differentiation among populations except Kiradu population. Nei's gene diversity (h) ranged between 0.082 and 0.193 with total diversity at species level is 0.294. Shannon's information index (I) ranged between 0.118 and 0.275 with an overall diversity of 0.439. Analysis of molecular variance showed more diversity among population level (56.65%) than at within population level (43.35%). The low gene flow value (Nm = 0.349) and high coefficient of genetic differentiation (G_ST = 0.589) and high fixation index (F_ST = 0.566) demonstrated elevated genetic differentiation among the population and can be predicted that these populations are not in Hardy–Weinberg proportions. Principal Co-ordinate Analysis confirms that Akal population has become phylogenetically more distinct and less diverse than the rest of the samples. Mantel's test revealed no correlation between genetic and geographical distances of populations (R² = 0.122). Overall highest diversity was observed in the population of Machiya Safari Park and Kiradu, while lowest in Akal population, later may constitute an evolutionary significant unit, having merit for special management.     

Extensive Variation and Sub-Structuring in Lineage A mtDNA in Indian Sheep: Genetic Evidence for Domestication of Sheep in India

Previous studies on mitochondrial DNA analysis of sheep from different regions of the world have revealed the presence of two major- A and B, and three minor- C, D and E maternal lineages. Lineage A is more frequent in Asia and lineage B is more abundant in regions other than Asia. We have analyzed mitochondrial DNA sequences of 330 sheep from 12 different breeds of India. Neighbor-joining analysis revealed lineage A, B and C in Indian sheep. Surprisingly, multidimensional scaling plot based on F_ST values of control region of mtDNA sequences showed significant breed differentiation in contrast to poor geographical structuring reported earlier in this species. The breed differentiation in Indian sheep was essentially due to variable contribution of two major lineages to different breeds, and sub- structuring of lineage A, possibly the latter resulting from genetic drift. Nucleotide diversity of this lineage was higher in Indian sheep (0.014 ± 0.007) as compared to that of sheep from other regions of the world (0.009 ± 0.005 to 0.01 ± 0.005). Reduced median network analysis of control region and cytochrome b gene sequences of Indian sheep when analyzed along with available published sequences of sheep from other regions of the world showed that several haplotypes of lineage A were exclusive to Indian sheep. Given the high nucleotide diversity in Indian sheep and the poor sharing of lineage A haplotypes between Indian and non-Indian sheep, we propose that lineage A sheep has also been domesticated in the east of Near East, possibly in Indian sub-continent. Finally, our data provide support that lineage B and additional lineage A haplotypes of sheep might have been introduced to Indian sub-continent from Near East, probably by ancient sea trade route.