Molecular characterization of GPR50 gene and study of its comparative genetic variability in sheep breeds adapted to different thermo-contrasting climatic regimens.

Abstract

GPR50, formerly known as a melatonin-related receptor, is one of the three subtypes of melatonin receptor subfamily, together with MTNR1A and MTNR1B. GPR50, despite its high identity with the melatonin receptor family, does not bind melatonin and is considered to be an ortholog of MTNR1C in mammals. GPR50-expressing cells have been found in the dorsomedial nucleus of the hypothalamus, the periventricular nucleus, and the median eminence. Genetic and functional evidence have been recently investigated linking GPR50 to adaptive thermogenesis and torpor, but still, it is an orphan receptor and is yet to be studied conclusively. The aims of the study were to characterize the GPR50 gene of sheep and to study the sequence variability of the gene in Indian sheep breeds of two different thermo-varied agroclimatic conditions. Genomic DNA isolation was done and a 791-bp sequence was amplified using self-designed primers and SNP profiling done out of samples of all the breeds to study the relative frequency of SNPs in each of the breed. Five important non-synonymous mutations were observed in the various breeds studied. T698G, G1097A, G1270A, G1318A, and C1334G lead to the following substitution: valine by glycine, arginine by glutamine, threonine by alanine, isoleucine by valine, and serine by cytosine, respectively. Two synonymous mutations (T663G and C888T) were also observed in some of the studied breeds. G1270A and C888T were the most prevalent SNPs observed in nearly all of the breeds. C888T SNPs were observed in higher prevalence in Chokla, Marwari, and Magra in comparison to Gaddi and Bharat Merino. A PolyPhen-2 analysis, which is used to assess the potential damaging nature of an SNP, revealed that mutation T698G and G1270A were benign while G1097A, G1318A, and C1334G were damaging with a score of 0.987, 0.993, and 0.739, respectively. A 3-D homology model of the protein was prepared using c4zwjA (UniProt sequence ID) as a template using the online version of Phyre2 protein modeling software. The structure demonstrated closed similarity with other G-coupled receptor and it had a 45% α-helical content. G1270A and C888T may be taken up for SNP correlation in a larger population study for their association with heat stress protection.