SNP rs403212791 in exon 2 of the MTNR1A gene is associated with reproductive seasonality in the Rasa Aragonesa sheep breed

The aim of this study was to characterize and identify causative SNPs in the MTNR1A gene responsible for the reproductive seasonality traits in the Rasa aragonesa sheep breed. A total of 290 ewes (155, 84 and 51 mature, young and ewe lambs, respectively) from one flock were controlled from January to August. The following three reproductive seasonality traits were considered: the total days of anoestrus (TDA) and the progesterone cycling months (P4CM); both ovarian function seasonality traits based on blood progesterone levels; and the oestrus cycling months (OCM) based on oestrous detection, which indicate behavioural signs of oestrous. We have sequenced the total coding region plus 733 and 251 bp from the promoter and 3'-UTR regions, respectively, from the gene in 268 ewes. We found 9 and 4 SNPs associated with seasonality traits in the promoter (for TDA and P4CM) and exon 2 (for the three traits), respectively. The SNPs located in the gene promoter modify the putative binding sites for various trans-acting factors. In exon 2, two synonymous SNPs affect RFLP sites, rs406779174/RsaI (for the three traits) and rs430181568/MnlI (for OCM), and they have been related with seasonal reproductive activity in previous association studies with other breeds. SNP rs400830807, which is located in the 3'-UTR, was associated with the three traits, but this did not modify the putative target sites for ovine miRNAs according to in silico predictions. Finally, the SNP rs403212791 (NW_014639035.1: g.15099004G > A), which is also associated with the three seasonality phenotypes, was the most significant SNP detected in this study and was a non-synonymous polymorphism, leading a change from an Arginine to a Cysteine (R336C). Haplotype analyses confirmed the association results and showed that the effects found for the seasonality traits were caused by the SNPs located in exon 2. We have demonstrated that the T allele in the SNP rs403212791 in the MNTR1A gene is associated with a lower TDA and higher P4CM and OCM values in the Rasa Aragonesa breed.