Discovery and computational exploration of SNPs in GNRHR gene and their influence on protein structure and function in Indian goat breeds

This study focuses on identifying functional single nucleotide polymorphisms (SNPs) within the gonadotropin-releasing hormone receptor (GNRHR) gene and conducting subsequent in-silico analysis of their effects on protein structure and function in two distinct South Indian goat breeds, namely Attapady Black (n = 120) and Malabari goats (n = 180), known for their divergent prolificacy traits. Utilizing a DNA pool sequencing assay, ten SNPs were uncovered in the study population: c.-1129T>G, c.-1069A>G, c.-978A>C, c.-605A>G, c.-33A>G, c.-29T>G, c.48G>A, c.75G>A, c.209T>G, and c.*212A>G. Notably, two polymorphisms, c.-1129T>G and c.-33A>G, were novel. Additionally, two polymorphisms, c.-33A>G and c.-978A>C, were exclusive to Malabari goats. Analysis of upstream variants revealed modifications to transcription factor and micro-RNA (miRNA) binding sites, suggesting potential alterations in GNRHR expression. Of particular significance was the non-synonymous exonic variant at c.209T>G locus, resulting in methionine to arginine substitution at the 70^th position within the first intracellular loop of the receptor protein. This amino acid change may have implications for the functional dynamics of the receptor as GnRHR intracellular loops are involved in G protein coupling thereby facilitation of downstream signalling pathways. The identified SNPs and their in-silico impact analysis contribute to our understanding of the molecular mechanisms underlying reproductive traits in these goat populations, with implications for future breeding strategies and genomic selection programs.