Inference of Genetic Structure and the Process of Population Formation in Nepalese Native Goats Using Uniparental and Genome-Wide Markers.

Abstract

Nepal is a small, landlocked country with marked elevational variation from the Terai plains to the Himalayas. Here, four indigenous goat populations (Chyangra, Sinhal, Khari, and Terai) are raised at different elevations. This study aimed to clarify the genetic structure of these populations and how they are formed and propagated across the Himalayan region. We analyzed 136 Nepalese goats using mitochondrial (mt) DNA D-loop and sex-determining region Y (SRY) 3'-untranslated region (UTR) sequences, as well as 50 K SNP array data. The mtDNA haplogroups D (0.162) and G (0.03) were detected only in Chyangra, whereas haplogroup B was predominant in Sinhal (0.42), followed by Khari (0.260). Regarding SRY haplotypes, Y2B was detected in all populations, whereas Y1AB (0.42) was found only in Chyangra. Genome-wide SNP analysis showed that Chyangra was genetically related to Tibetan and Central Asian goats, while Terai resembled South Asian goats. Interestingly, Sinhal formed a distinct cluster, whereas Khari exhibited an admixed genetic structure. These findings suggest that Nepalese goats originate from at least three ancestral lineages and that an additional migration route may have existed through the southern Himalayas.