Evidence for independent domestication of sheep mtDNA lineage A in India and introduction of lineage B through Arabian sea route

India ranks the second in the world in terms of its sheep population with approximately 74.26 million represented by 44 well-described breeds in addition to several non-descript populations. Genetic diversity and phylogeography of Indian sheep breeds remain poorly understood, particularly for south Indian breeds. To have a comprehensive view of the domestication history of Indian sheep, we sequenced the mitochondrial DNA (mtDNA) control region (D-loop) and cytochrome b gene (CYTB) of 16 Indian domestic sheep breeds, most of them (13) from the south India. We analysed these sequences along with published data of domestic and wild sheep from different countries, including India. The haplotype diversity was relatively high in Indian sheep, which were classified into the three known mtDNA lineages, namely A, B and C. Lineage A was predominant among Indian sheep whereas lineages B and C were observed at low frequencies but C was restricted to the breeds of north and east India. The median joining network showed five major expanding haplogroups of lineage A (A1–A5). Out of which, A2, A4 and A5 were more frequent in Indian sheep in contrast to breeds from other parts of the world. Among the 27 Indian sheep breeds analysed, Mandya and Sonadi breeds were significantly different from other Indian breeds in the MDS analyses. This was explained by a very high contribution of lineage B into these two breeds. The Approximate Bayesian Computation (ABC) provided evidence for the domestication of lineage A sheep in the Indian subcontinent. Contrary to the current knowledge, we also found strong support for the introduction of lineage B into Indian subcontinent through sea route rather than from the Mongolian Plateau. The neighbour-joining tree of domestic and wild sheep revealed the close genetic relationship of Indian domestic sheep with Pakistani wild sheep O. vignei blanfordi. Based on our analyses and archaeological evidences, we suggest the Indian subcontinent as one of the domestication centres of the lineage A sheep, while lineage B sheep might have arrived into India from elsewhere via Arabian sea route. To the best of our knowledge, this is the first comprehensive study on Indian sheep where we have analysed more than 740 animals belonging to 27 sheep breeds raised in various regions of India. Our study provides insight into the understanding of the origin and migratory history of Indian sheep.

Association analysis of novel SNPs in exon 10 of the growth hormone receptor gene with growth traits in Indian sheep

The growth hormone receptor (GHR) gene encodes the type I cytokine receptor that helps in joining the growth hormone to this receptor, thus promoting receptor-dimerization, leading to up-regulating growth. The ovine GHR gene located on chromosome 16, which consists of 10 exons and 9 introns, along with untranslated regions on either side, comprise a total size of 178.09 kb. However, earlier reports about polymorphism have mainly dealt with exon 10 which is also a larger fragment of this gene comprising 1102 bp. Hence, this study was carried out to detect polymorphism in exon 10 of the GHR gene and its association with growth traits. Genomic DNA was isolated from blood samples of Madras Red and Mecheri sheep breeds from India. Part of exon 10 (895 bp) of the GHR gene was amplified and sent for sequencing. The sequence analysis revealed transition of nucleotide G>A at loci G177624A and G177878A in both sheep breeds. Populations were screened by Tetra-primer ARMS-PCR. The genotype frequencies of GG, GA and AA were 0.276, 0.519 and 0.205 at 177624 G>A, and 0.307, 0.444 and 0.149 at 177878 G>A in Madras Red sheep; whereas in Mecheri they were 0.476, 0.372 and 0.152 at 177624 G>A, and 0.629, 0.314 and 0.057 at 177878 G>A, respectively. Likewise the estimated allele frequencies of G and A were 0.5355 and 0.4645 at 177624 G>A, and 0.5790 and 0.4210 at 177878 G>A in Madras Red sheep; whereas in Mecheri they were 0.6620 and 0.3380 at 177624 G>A, and 0.7860 and 0.2140 at 177878 G>A, respectively. The effect of sex was significant for birth, six and nine month weight; but non-significant for three and 12 month weight in Mecheri sheep. However, in the Madras Red breed the effect of sex was significant for all body weights except weaning weight. The effect of variations on growth traits, viz., birth weight, weight at weaning, and weight at six, nine and twelve months in both breeds were analysed for their association, and they were found non-significant. Since these SNPs are salient findings of GHR gene polymorphism in Indian sheep breeds, further investigation is required into the significant effects of these novel SNPs, which could be useful for genetic improvement based on marker assisted selection.

Analysis of selection signatures reveals important insights into the adaptability of high-altitude Indian sheep breed Changthangi

Changthangi is a high-altitude sheep breed of India that is adapted to cold and hypoxic climate of Himalayas. In the present study, we analysed population structure of Changthangi and contrasted it with selected Indian and European commercial sheep breeds to detect genomic regions under positive selection. The studied domesticated sheep breeds included Changthangi, Indian Garole, Deccani, Tibetan, Rambouillet and Australian Merino. While the PCA results revealed Changthangi clustered with Tibetan sheep; TREEMIX and ADMIXTURE results also detected the introgression of lowland Indian sheep inheritance in Changthangi. Cross-population comparisons of Changthangi using XP-EHH showed multiple functional regions present on OAR 7, 15 and 16, to be under selection in Changthangi sheep. These regions are related to adaptation to climatic and hypoxic stressors, nervous system and mitochondrial functioning. The genes present in these regions are suitable candidates for future studies on the genetic mechanisms underlying high-altitude adaptation.

Differential expression of miRNAs in skeletal muscles of Indian sheep with diverse carcass and muscle traits

The study presents the miRNA profiles of two Indian sheep populations with divergent carcass and muscle traits. The RNA sequencing of longissimus thoracis muscles from the two populations revealed a total of 400 known miRNAs. Myomirs or miRNAs specific to skeletal muscles identified in our data included oar-miR-1, oar-miR-133b, oar-miR-206 and oar-miR-486. Comparison of the two populations led to identification of 100 differentially expressed miRNAs (p < 0.05). A total of 45 miRNAs exhibited a log2 fold change of ≥ ( ±) 3.0. Gene Ontology analysis revealed cell proliferation, epithelial to mesenchymal transition, apoptosis, immune response and cell differentiation as the most significant functions of the differentially expressed miRNAs. The differential expression of some miRNAs was validated by qRT-PCR analysis. Enriched pathways included metabolism of proteins and lipids, PI3K-Akt, EGFR and cellular response to stress. The microRNA-gene interaction network revealed miR-21, miR-155, miR-143, miR-221 and miR-23a as the nodal miRNAs, with multiple targets. MicroRNA-21 formed the focal point of the network with 42 interactions. The hub miRNAs identified in our study form putative regulatory candidates for future research on meat quality traits in Indian sheep. Our results provide insight into the biological pathways and regulatory molecules implicated in muscling traits of sheep.