# Analysis of MTNR1A Genetic Polymorphisms and Their Association with the Reproductive Performance Parameters in Two Mediterranean Sheep Breeds

Asma Arjoune, Abrar B Alsaleh, Safia A Messaoudi, Hanen Chelbi, Refka Jelassi, Mourad Assidi, Taha Najar, Brahim Haddad, Marc-André Sirard

► To cite this version:

Asma Arjoune, Abrar B Alsaleh, Safia A Messaoudi, Hanen Chelbi, Refka Jelassi, et al.. Analysis of MTNR1A Genetic Polymorphisms and Their Association with the Reproductive Performance Parameters in Two Mediterranean Sheep Breeds. Animals, 2023, 13 (3), pp.448. 10.3390/ani13030448. pasteur-05011580

HAL Id: pasteur-05011580

https://riip.hal.science/pasteur-05011580v1

Submitted on 30 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Analysis of MTNR1A Genetic Polymorphisms and Their Association with the Reproductive Performance Parameters in Two Mediterranean Sheep Breeds

Asma Arjoune <fnref n="1" />,<fnref n="2" />, Abrar B. Alsaleh <fnref n="3" />, Safia A. Messaoudi ³ID, Hanen Chelbi ⁴ID, Refka Jelassi ⁴ID, Mourad Assidi ⁵ID, Taha Najar <fnref n="2" />, Brahim Haddad <fnref n="2" /> and Marc-André Sirard <fnref n="1" />,*ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Arjoune, A.; Alsaleh, A.B.; Messaoudi, S.A.; Chelbi, H.; Jelassi, R.; Assidi, M.; Najar, T.; Haddad, B.; Sirard, M.-A. Analysis of *MTNR1A* Genetic Polymorphisms and Their Association with the Reproductive Performance Parameters in Two Mediterranean Sheep Breeds.
*Animals* **2023**, *13*, 448. https://doi.org/10.3390/ani13030448

Academic Editors: Rosaura Pérez-Pe and Serafeim Papadopoulos

Received: 13 December 2022

Revised: 19 January 2023

Accepted: 24 January 2023

Published: 28 January 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** The aim of this research was to study the effect of *MTNR1A* gene polymorphisms on the reproductive performance of two Mediterranean sheep breeds. A total of 26 SNPs were found, and three SNPs caused amino acid changes. Two SNPs were totally linked and correlated with reproductive activity resumption. This study consists of a piece of new useful knowledge for the control of sheep reproductive seasonality and farm management that could improve genetic selection, especially in countries where agriculture relies mainly on sheep farming.

**Abstract:** Sheep farming plays an important economic role, and it contributes to the livelihoods of many rural poor in several regions worldwide and particularly in Tunisia. Therefore, the steady improvement of ewes' reproductive performance is a pressing need. The *MTNR1A* gene has been identified as an important candidate gene that plays a key role in sheep reproduction and its sexual inactivity. It is involved in the control of photoperiod-induced seasonality mediated by melatonin secretion. The aim of this study was to identify SNPs in the *MTNR1A* gene in two Tunisian breeds, Barbarine (B) and Queue Fine de l'Ouest (QFO). DNA extracted from the blood of 77 adult ewes was sequenced. Selected ewes were exposed to adult fertile rams. A total of 26 SNPs were detected; 15 SNPs in the promoter region and 11 SNPs in the exon II were observed in both (B) and (QFO) breeds. The SNP rs602330706 in exon II is a novel SNP detected for the first time only in the (B) breed. The SNPs rs430181568 and rs40738822721 (SNP18 and SNP20 in our study, respectively) were totally linked in this study and can be considered a single marker. DTL was associated with SNP18 and SNP20 in (B) ewes (*p* < 0.05); however, no significant difference was detected between the three genotypes (G/G, G/A, and A/A) at these two SNPs. Fertility rate and litter size parameters were not affected by SNP18 and SNP20. There was an association between these two polymorphisms and (B) lambs' birth weights (*p* < 0.05). Furthermore, the ewes with the A/A genotype gave birth to lambs with a higher weight compared to the other two genotypes for this breed (*p* < 0.05). There was not an association between SNP 18 and SNP20 and (QFO) ewes' reproductive parameters. These results might be considered in future sheep selection programs for reproductive genetic improvement.

**Keywords:** *MTNR1A*; SNPs; sheep; reproductive parameters

<footnote n="1">¹ Centre de Recherche en Reproduction, Développement et Santé Intergénérationnelle, Département des Sciences Animales, Faculté des Sciences de l'Agriculture et de l'Alimentation, Université Laval, Quebec, QC G1V 0A6, Canada</footnote>

<footnote n="2">² Department of Animal Production, National Agronomic Institute of Tunisia, University of Carthage, Tunis 1082, Tunisia</footnote>

<footnote n="3">³ Department of Forensic Sciences, Naïf Arab University for Security Sciences, Riyadh 14812, Saudi Arabia</footnote>

<footnote n="4">⁴ Laboratory of Medical Parasitology, Biotechnologies, and Biomolecules, Pasteur Institute of Tunis, Tunis 1002, Tunisia</footnote>

<footnote n="5">5 Center of Excellence in Genomic Medicine Research (CEGMR), Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21589, Saudi Arabia</footnote>

<footnote>* Correspondence: marc-andre.sirard@fsaa.ulaval.ca</footnote>
# 1. Introduction

The most expressed seasonal rhythm for small ruminants living in the Mediterranean areas is reproductive activity. It is a form of adaptation that several sheep breeds exert in the face of changes in environmental factors [1]. Ewes' reproduction seasonality is manifested by the existence of the seasonal anestrus period (the period in days from ram placement with ewes to lambing); it is controlled by the photoperiod [2–4]. Photoperiod annual variation is one of the main environmental factors that synchronizes reproductive activity and melatonin secretion [5]. Therefore, the short-days season is characterized by shorter photoperiods and promotes nocturnal synthesis of melatonin that stimulates their sexual activity [6]. Furthermore, the melatonin nocturnal production time in fall and winter is longer than in spring and winter [7]. It exerts a fundamental function in the modulation of sheep reproduction rhythm [8] through specific receptors located in the hypothalamic, pituitary, and gonadal axis [9]. Two high-affinity melatonin receptors have been identified in mammals: melatonin receptor 1A (MTNR1A, MT1) and melatonin receptor 1B (MTNR1B, MT2) [10,11]. MTNR1A is located on chromosome 26 in sheep. So far, several scientists have explored polymorphisms on the MTNR1A gene. In addition, these mutations were associated with the reproductive seasonality traits in different sheep breeds [12–18]. However, with MTNR1B, although it is less studied, it was reported that polymorphism in this receptor gene could affect ewes' reproductive seasonality and litter size [19]. Two single nucleotide polymorphisms (SNP) at positions 606 and 612 have been associated with different reproductive parameters [20–22]; however, in the Ile de France breed, there is no association between these SNPs and reproductive traits [23]. Discovering a polymorphic site in the MTNR1A gene has provided new useful knowledge for the control of sheep reproductive seasonality and farm management that could improve genetic selection, especially in countries where agriculture relies mainly on sheep farming. It will be beneficial to study the MTNR1A gene polymorphism in two main sheep breeds in Tunisia, Barbarine (B) and Queue Fine de l'ouest (QFO). Therefore, this research aims at studying MTNR1A polymorphisms in two Tunisian sheep breeds and investigating their association with the ewes' reproductive performance parameters.

# 2. Materials and Methods

## 2.1. Ethical Approval

The care and the use of animals were under the control of the veterinarians of the Livestock and Pastures office (LPO) of Tunisia in accordance with the guidelines of the Animal Ethics Committee of the National School of Veterinary Medicine (Ethical approval no: CEEA-ENMV 57/2022). Blood samples were collected by a veterinarian during routine flock health control procedures.

## 2.2. Animals and Management

Ewes (n = 77) of meat-type breeding were used in this study. The ewes were of two Tunisian autochthonous breeds: Barbarine (n = 42) and Queue Fine de l'ouest (n = 35). These two breeds are the most typical purebred breeds used in sheep breeding. The animals were produced and located in the northeast of Tunisia (Zaghouan Governorate) at the LPO farms; the farming was based on a semi-extensive system located at a latitude of 36.3667 and a longitude: of 9.9 36°22′0″ North and 9°54′0″ East with a Mediterranean climate and a hot summer. All the ewes selected for the present study were multiparous aged from 3 to 5 years (3.825 ± 0.785) with body condition score (BCS) from 3.5 to 4 and had been maintained in a natural photoperiod condition from birth. The ewes were characterized by seasonal breeding; they show a period of anoestrus from March to September. They can lamb once a year from January to March. In addition, the ewes were clinically examined and their BCS on a 1 to 5 scale was assessed by a veterinarian.
## 2.3. Reproductive Data Collection

Rams of the two breeds were introduced in female flocks on the 1st of May of each year, so lambing can occur in October where there are more grazing areas, ensuring adequate milk production and lamb development. Rams were moved away from the ewe flocks after 90 days. For each breed, days to lambing (DTL), the number of lambs born for each ewe, and lambs' birth weights were recorded for 2 years/campaigns: 2017/2018 and 2018/2019. From the farm data records, the following reproductive variables for both breeds were collected: DTL, number of lambed ewes, litter size, and lambs' birth weights.

## 2.4. Sampling and Sequencing

Ten (10) mL of blood from the jugular vein of each ewe was collected using vacutainer collection EDTA tubes. Total DNA was extracted using the DNeasy blood and tissue kit (Qiagen, Mississauga, ON, Canada) according to the manufacturer's protocol. All the primers used for the amplifications of the promoter and exon II regions were those reported by Luridiana et al. [12]. The polymerase chain reaction (PCR) and the electrophoresis conditions were performed following Luridiana et al. [12]. All PCR products were sequenced in forward direction by a commercial service (Macrogen, Seoul, Korea). Information about primers and fragment amplifications are available in Supplementary Table S17.

## 2.5. Sequencing Data Analysis

DNA sequencing results were visualized using Series Data Collection Software 3130 (Thermo Fisher Scientific, Waltham, MA, USA) and analyzed using Sequencing Analysis software v5.3.1 (Thermo Fisher Scientific, USA). Genetic mutations were detected using Seqscape® software v2.6 (Thermo Fisher Scientific, USA) [24]. The sequences alignment was conducted via MEGA X: Molecular Evolutionary Genetics Analysis with the reference sequence of the sheep genome AY524665.1 and NM_001009725 from GenBank.

## 2.6. Statistical Analysis

Statistical analysis of the results was performed using the software SPSS Statistics version 28 (IBM, New York, NY, USA). Frequencies of genetic changes were compared using the Chi-squares ($X_2$) test and in case of non-validity the Fisher's exact tests. Accordingly, comparisons of variables between 2 groups were carried out using ANOVA test for normally distributed variables. As for the continuous variables, they were expressed in these mean ± standard deviation (SD). The value of $p \le 0.05$ was considered statistically significant.

# 3. Results

The sequences of *MTNR1A* promoter and exon II regions from the total sample were compared with the reference sequence AY524665.1 and NM_001009725 from Gen-Bank. All genotype positions were inserted according to the sheep genome assembly version Oar_rambouillet_v1.0:26:17354220:17378573 (Oar_Rambouillet_v1.0, accession Ref-Seq GCF_002742125.1). As shown in Table 1, an overall total of 26 SNPs were detected; 15 SNPs in the promoter region, and 11 SNPs in the exon II were observed in both (B) and (QFO) breed. The SNPs number, breeds identifiers, location, locus, and amino acid change in the MTNR1A gene are summarized in Table 1. The allele and genotype frequency and the Odds Ratio (OR) are included in Table 2. All SNPs were in Hardy-Weinberg equilibrium except SNP6 (rs415456480) in both breeds ($p = 0.004$ and $p \le 0.0001$ for (B) and (QFO), respectively) and SNP7 (rs406184829) within the promoter region in (QFO) ewes ($p = 0.00027$) due to their high heterozygosity rate. Within exon II, the most studied region, five SNPs were missense, and the others were synonymous. Three SNPs caused amino acid changes, (rs407388227 G > A and rs404378206 G > A causing Val/Ile and rs416266900 A > C causing Ala/Asp substitution) in the amino acid sequence. The others twenty-three SNPs detected in our study were silent.
<table><caption>**Table 1.** SNPs number, dbSNPs identifiers, location (using the genome version Oar_rambouillet_v1.0 as a reference), *p*-values of the Hardy–Weinberg equilibrium, locus, and amino acid change in the *MTNR1A* gene. B: Barbarine, QFO: queue fine de l'ouest.</caption><thead><tr><th rowspan="2">SNP Number</th><th rowspan="2">dbSNPs</th><th rowspan="2">Position</th><th colspan="3">HWE (p-Value)</th><th rowspan="2">MTNR1A Region</th><th rowspan="2">A.A' Change</th></tr><tr><th>B</th><th>QFO</th><th>All Subject</th></tr></thead><tbody><tr><td>SNP1</td><td>rs428941001</td><td>Chr 26: 17379252</td><td>0.52</td><td>0.23</td><td>0.27</td><td>Promoter</td><td>–</td></tr><tr><td>SNP2</td><td>rs419424336</td><td>Chr 26: 17379237</td><td>1.00</td><td>0.31</td><td>0.15</td><td>Promoter</td><td>–</td></tr><tr><td>SNP3</td><td>rs429917252</td><td>Chr26: 17379178</td><td>1.00</td><td>0.31</td><td>0.15</td><td>Promoter</td><td>–</td></tr><tr><td>SNP4</td><td>rs405080439</td><td>Chr 26: 17379097</td><td>0.52</td><td>0.76</td><td>0.51</td><td>Promoter</td><td>–</td></tr><tr><td>SNP5</td><td>rs428880789</td><td>Chr 26: 17379083</td><td>1.00</td><td>0.76</td><td>1.00</td><td>Promoter</td><td>–</td></tr><tr><td>SNP6</td><td>rs415456480</td><td>Chr 26: 17379037</td><td>0.004</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>Promoter</td><td>–</td></tr><tr><td>SNP7</td><td>rs406184829</td><td>Chr 26: 17378992</td><td>0.19</td><td>0.0003</td><td>&lt;0.0001</td><td>Promoter</td><td>–</td></tr><tr><td>SNP8</td><td>rs412826644</td><td>Chr 26: 17378833</td><td>1.00</td><td>1.00</td><td>1.00</td><td>Promoter</td><td>–</td></tr><tr><td>SNP9</td><td>rs411931887</td><td>Chr 26: 17378874</td><td>1.00</td><td>0.35</td><td>0.67</td><td>Promoter</td><td>–</td></tr><tr><td>SNP10</td><td>rs402949406</td><td>Chr 26: 17378871</td><td>1.00</td><td>0.54</td><td>1.00</td><td>Promoter</td><td>–</td></tr><tr><td>SNP11</td><td>rs426266687</td><td>Chr 26: 17378842</td><td>1.00</td><td>0.54</td><td>0.83</td><td>Promoter</td><td>–</td></tr><tr><td>SNP12</td><td>rs399461430</td><td>Chr 26: 17378769</td><td>1.00</td><td>0.54</td><td>1.00</td><td>Promoter</td><td>–</td></tr><tr><td>SNP13</td><td>rs400561563</td><td>Chr 26: 17378728</td><td>1.00</td><td>0.57</td><td>0.35</td><td>Promoter</td><td>–</td></tr><tr><td>SNP14</td><td>rs419743392</td><td>Chr 26: 17378706</td><td>1.00</td><td>0.36</td><td>0.82</td><td>Promoter</td><td>–</td></tr><tr><td>SNP15</td><td>rs406334919</td><td>Chr 26: 17378624</td><td>0.55</td><td>0.51</td><td>0.15</td><td>Promoter</td><td>–</td></tr><tr><td>SNP16</td><td>rs419680097</td><td>Chr 26:17355611</td><td>0.59</td><td>0.13</td><td>0.13</td><td>Exon2</td><td>–</td></tr><tr><td>SNP17</td><td>rs406779174</td><td>Chr 26:17355458</td><td>0.20</td><td>0.31</td><td>0.018</td><td>Exon2</td><td>–</td></tr><tr><td>SNP18</td><td>rs430181568</td><td>Chr 26:17355452</td><td>0.59</td><td>0.59</td><td>0.27</td><td>Exon2</td><td>–</td></tr><tr><td>SNP19</td><td>rs420819884</td><td>Chr 26:17355389</td><td>1.00</td><td>0.17</td><td>0.086</td><td>Exon2</td><td>–</td></tr><tr><td>SNP20</td><td>rs407388227</td><td>Chr 26:17355358</td><td>0.95</td><td>0.13</td><td>0.13</td><td>Exon2</td><td>Val/lle</td></tr><tr><td>SNP21</td><td>rs427019119</td><td>Chr 26:17355281</td><td>0.70</td><td>0.57</td><td>1.00</td><td>Exon2</td><td>–</td></tr><tr><td>SNP22</td><td>rs417800445</td><td>Chr 26:17355263</td><td>1.00</td><td>0.57</td><td>0.51</td><td>Exon2</td><td>–</td></tr><tr><td>SNP23</td><td>rs404378206</td><td>Chr 26:17355190</td><td>1.00</td><td>1.00</td><td>1.00</td><td>Exon2</td><td>Val/lle</td></tr><tr><td>SNP24</td><td>rs429718221</td><td>Chr 26:17355173</td><td>1.00</td><td>0.74</td><td>0.65</td><td>Exon2</td><td>–</td></tr><tr><td>SNP25</td><td>rs416266900</td><td>Chr 26:17355171</td><td>1.00</td><td>0.58</td><td>0.72</td><td>Exon2</td><td>Ala/Asp</td></tr><tr><td>SNP26</td><td>rs602330706</td><td>Chr 26:17355149</td><td>1.00</td><td>1.00</td><td>1.00</td><td>Exon2</td><td>–</td></tr></tbody></table>
<table><caption>**Table 2.** Allele and genotype frequency of the *MTNR1A* gene SNPs in the two sheep breeds. All genotype positions were inserted according to the sheep genome assembly version Oar_rambouillet_v1.0:26:17354220:17378573 (Oar_Rambouillet_v1.0, accession RefSeq GCF_002742125.1).</caption><thead><tr><th rowspan="2">SNP Number</th><th rowspan="2">dbSNPs</th><th rowspan="2">Genotype</th><th colspan="2">Genotype Frequency (%)</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">p-Value</th><th rowspan="2">Allele</th><th colspan="2">Allele Frequency</th></tr><tr><th>B</th><th>Q</th><th>B</th><th>Q</th></tr></thead><tbody><tr><td rowspan="3">SNP1</td><td rowspan="3">rs428941001</td><td>GG</td><td>33 (13)</td><td>14 (6)</td><td>1.00</td><td rowspan="3">0.075</td><td>G</td><td>0.60</td><td>0.44</td></tr><tr><td>GT</td><td>54 (21)</td><td>60 (25)</td><td>2.58 (0.83-7.97)</td><td>T</td><td>0.4</td><td>0.56</td></tr><tr><td>TT</td><td>13 (5)</td><td>26 (11)</td><td>4.77<br>(1.14-19.98)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP2</td><td rowspan="3">rs419424336</td><td>CC</td><td>72 (28)</td><td>45 (19)</td><td>1.00</td><td rowspan="3">0.018</td><td>C</td><td>0.85</td><td>0.64</td></tr><tr><td>CT</td><td>26 (10)</td><td>38 (16)</td><td>2.36 (0.88-6.29)</td><td>T</td><td>0.15</td><td>0.36</td></tr><tr><td>TT</td><td>3 (1)</td><td>17 (7)</td><td>10.32<br>(1.17-90.78)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP3</td><td rowspan="3">rs429917252</td><td>AA</td><td>3 (1)</td><td>17 (7)</td><td>10.32<br>(1.17-90.78)</td><td rowspan="3">0.018</td><td>G</td><td>0.85<br>(66)</td><td>0.64 (54)</td></tr><tr><td>GA</td><td>26 (10)</td><td>38 (16)</td><td>2.36 (0.88-6.29)</td><td>A</td><td>0.15</td><td>0.36</td></tr><tr><td>GG</td><td>72 (28)</td><td>45 (19)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP4</td><td rowspan="3">rs405080439</td><td>CC</td><td>13 (33)</td><td>8 (19)</td><td>1.00</td><td rowspan="3">0.170</td><td>C</td><td>0.6</td><td>0.46</td></tr><tr><td>CT</td><td>54 (21)</td><td>55 (23)</td><td>1.78 (0.62-5.14)</td><td>T</td><td>0.4</td><td>0.54</td></tr><tr><td>TT</td><td>13 (5)</td><td>26 (11)</td><td>3.57<br>(0.90-14.15)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP5</td><td rowspan="3">rs428880789</td><td>CC</td><td>41 (16)</td><td>19 (8)</td><td>1.00</td><td rowspan="3">0.063</td><td>C</td><td>0.64</td><td>0.46</td></tr><tr><td>CT</td><td>46 (18)</td><td>55 (23)</td><td>2.56 (0.90-7.30)</td><td>T</td><td>0.36</td><td>0.54</td></tr><tr><td>TT</td><td>13 (5)</td><td>26 (11)</td><td>4.40<br>(1.13-17.07)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP6</td><td rowspan="3">rs415456480</td><td>CC</td><td>77 (30)</td><td>67 (28)</td><td>1.00</td><td rowspan="3">0.150</td><td>C</td><td>0.83</td><td>0.7</td></tr><tr><td>CT</td><td>13 (5)</td><td>7 (3)</td><td>0.64 (0.14-2.94)</td><td>T</td><td>0.17</td><td>0.3</td></tr><tr><td>TT</td><td>10 (4)</td><td>26 (11)</td><td>2.95<br>(0.84-10.34)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP7</td><td rowspan="3">rs406184829</td><td>AA</td><td>87 (34)</td><td>69 (29)</td><td>1.00</td><td rowspan="3">0.056</td><td>A</td><td>0.92</td><td>0.76</td></tr><tr><td>AG</td><td>10 (4)</td><td>14 (6)</td><td>1.76 (0.45-6.84)</td><td>G</td><td>0.08</td><td>0.24</td></tr><tr><td>GG</td><td>3 (1)</td><td>17 (7)</td><td>8.21<br>(0.95-70.67)</td><td></td><td></td><td></td></tr><tr><td rowspan="2">SNP8</td><td rowspan="2">rs412826644</td><td>CC</td><td>100 (42)</td><td>88 (38)</td><td>1.00</td><td rowspan="2">0.008</td><td>C</td><td>1</td><td>0.94</td></tr><tr><td>CT</td><td>0.00</td><td>12 (5)</td><td>N/A</td><td>T</td><td>0.00</td><td>0.06</td></tr><tr><td rowspan="3">SNP9</td><td rowspan="3">rs411931887</td><td>AA</td><td>14 (6)</td><td>28 (12)</td><td>5.33<br>(1.37-20.71)</td><td rowspan="3">0.026</td><td>G</td><td>0.62</td><td>0.43</td></tr><tr><td>GA</td><td>48 (20)</td><td>58 (25)</td><td>3.33<br>(1.10-10.09)</td><td>A</td><td>0.38</td><td>0.57</td></tr><tr><td>GG</td><td>38 (16)</td><td>14 (6)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP10</td><td rowspan="3">rs402949406</td><td>GG</td><td>14 (6)</td><td>28 (12)</td><td>4.86<br>(1.30-18.13)</td><td rowspan="3">0.032</td><td>T</td><td>0.63</td><td>0.44</td></tr><tr><td>TG</td><td>45 (19)</td><td>56 (24)</td><td>3.07 (1.06-8.91)</td><td>G</td><td>0.37</td><td>0.56</td></tr><tr><td>TT</td><td>40 (17)</td><td>16 (7)</td><td>1.00</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th rowspan="2">SNP Number</th><th rowspan="2">dbSNPs</th><th rowspan="2">Genotype</th><th colspan="2">Genotype Frequency (%)</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">p-Value</th><th rowspan="2">Allele</th><th colspan="2">Allele Frequency</th></tr><tr><th>B</th><th>Q</th><th>B</th><th>Q</th></tr></thead><tbody><tr><td rowspan="3">SNP11</td><td rowspan="3">rs426266687</td><td>AA</td><td>14 (6)</td><td>28 (12)</td><td>4.57<br>(1.22-17.16)</td><td rowspan="3">0.05</td><td>G</td><td>0.62</td><td>0.44</td></tr><tr><td>GA</td><td>48 (20)</td><td>56 (24)</td><td>2.74 (0.94-7.98)</td><td>A</td><td>0.38</td><td>0.56</td></tr><tr><td>GG</td><td>38 (16)</td><td>16 (7)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP12</td><td rowspan="3">rs399461430</td><td>AA</td><td>40 (17)</td><td>16 (7)</td><td>1.00</td><td rowspan="3">0.032</td><td>A</td><td>0.63</td><td>38</td></tr><tr><td>AG</td><td>45 (19)</td><td>56 (24)</td><td>3.07 (1.06-8.91)</td><td>G</td><td>0.37</td><td>0.56</td></tr><tr><td>GG</td><td>14 (6)</td><td>28 (12)</td><td>4.86<br>(1.30-18.13)</td><td></td><td></td><td></td></tr><tr><td rowspan="2">SNP13</td><td rowspan="2">rs400561563</td><td>GA</td><td>22 (9)</td><td>33 (14)</td><td>1.72 (0.65-4.59)</td><td rowspan="2">0.27</td><td>A</td><td>0.89</td><td>0.83</td></tr><tr><td>GG</td><td>78 (31)</td><td>67 (28)</td><td>1.00</td><td>G</td><td>0.11</td><td>0.17</td></tr><tr><td rowspan="3">SNP14</td><td rowspan="3">rs419743392</td><td>CC</td><td>38 (15)</td><td>14 (6)</td><td>1.00</td><td rowspan="3">0.038</td><td>C</td><td>0.61</td><td>0.43</td></tr><tr><td>CT</td><td>48 (19)</td><td>57 (24)</td><td>3.16 (1.03-9.70)</td><td>T</td><td>0.39</td><td>0.57</td></tr><tr><td>TT</td><td>15 (6)</td><td>29 (12)</td><td>5.00<br>(1.28-19.53)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP15</td><td rowspan="3">rs406334919</td><td>AA</td><td>75 (30)</td><td>43 (18)</td><td>1.00</td><td rowspan="3">0.0051</td><td>A</td><td>0.86</td><td>0.63</td></tr><tr><td>AG</td><td>22 (9)</td><td>40 (17)</td><td>3.15 (1.16-8.53)</td><td>G</td><td>0.14</td><td>0.37</td></tr><tr><td>GG</td><td>2 (1)</td><td>17 (7)</td><td>11.67<br>(1.33-102.72)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP16</td><td rowspan="3">rs419680097</td><td>GG</td><td>70 (30)</td><td>70 (30)</td><td>1.00</td><td rowspan="3">0.88</td><td>G</td><td>0.83</td><td>0.81</td></tr><tr><td>GT</td><td>26 (11)</td><td>23 (10)</td><td>0.91 (0.34-2.46)</td><td>T</td><td>0.17</td><td>0.19</td></tr><tr><td>TT</td><td>5 (2)</td><td>7 (3)</td><td>1.50 (0.23-9.63)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP17</td><td rowspan="3">rs406779174</td><td>TT</td><td>19 (8)</td><td>49 (21)</td><td>1.00</td><td rowspan="3">0.0014</td><td>T</td><td>0.37</td><td>0.67</td></tr><tr><td>TC</td><td>37 (16)</td><td>37 (16)</td><td>0.38 (0.13-1.11)</td><td>C</td><td>0.63</td><td>0.33</td></tr><tr><td>CC</td><td>44 (19)</td><td>14 (16)</td><td>0.12 (0.04-0.41)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP18</td><td rowspan="3">rs430181568</td><td>GG</td><td>70 (30)</td><td>70 (30)</td><td>1.00</td><td rowspan="3">n/a</td><td>G</td><td>0.83</td><td>0.83</td></tr><tr><td>GA</td><td>26 (11)</td><td>26 (11)</td><td>1.00 (0.38-2.66)</td><td>A</td><td>0.7</td><td>0.7</td></tr><tr><td>AA</td><td>5 (2)</td><td>5 (2)</td><td>1.00 (0.13-7.57)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP19</td><td rowspan="3">rs420819884</td><td>GG</td><td>95 (41)</td><td>77 (33)</td><td>1.00</td><td rowspan="3">0.024</td><td>G</td><td>0.98</td><td>0.86</td></tr><tr><td>GA</td><td>5 (2)</td><td>19 (8)</td><td>4.97<br>(0.99-25.01)</td><td>A</td><td>0.02</td><td>0.14</td></tr><tr><td>AA</td><td>0.00</td><td>5 (2)</td><td>n/a</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP20</td><td rowspan="3">rs407388227</td><td>GG</td><td>70 (30)</td><td>70 (30)</td><td>1.00</td><td rowspan="3">0.88</td><td>G</td><td>0.83</td><td>0.81</td></tr><tr><td>GA</td><td>26 (11)</td><td>23 (10)</td><td>0.91 (0.34-2.46)</td><td>A</td><td>0.17</td><td>0.19</td></tr><tr><td>AA</td><td>4 (2)</td><td>7 (3)</td><td>1.50 (0.23-9.63)</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP21</td><td rowspan="3">rs427019119</td><td>AA</td><td>9 (4)</td><td>0.00</td><td>n/a</td><td rowspan="3">0.041</td><td>A</td><td>0.28</td><td>0.16</td></tr><tr><td>AG</td><td>37 (16)</td><td>33 (14)</td><td>0.69 (0.28-1.71)</td><td>G</td><td>0.72</td><td>0.84</td></tr><tr><td>GG</td><td>53 (23)</td><td>67 (29)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP22</td><td rowspan="3">rs417800445</td><td>AA</td><td>4 (2)</td><td>0.00</td><td>n/a</td><td rowspan="3">0.14</td><td>A</td><td>0.24</td><td>0.15</td></tr><tr><td>AG</td><td>40 (17)</td><td>30 (13)</td><td>0.61 (0.25-1.50)</td><td>G</td><td>0.76</td><td>0.85</td></tr><tr><td>GG</td><td>56 (24)</td><td>70 (30)</td><td>1.00</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th rowspan="2">SNP Number</th><th rowspan="2">dbSNPs</th><th rowspan="2">Genotype</th><th colspan="2">Genotype Frequency (%)</th><th rowspan="2">OR (95% CI)</th><th rowspan="2">p-Value</th><th rowspan="2">Allele</th><th colspan="2">Allele Frequency</th></tr><tr><th>B</th><th>Q</th><th>B</th><th>Q</th></tr></thead><tbody><tr><td rowspan="3">SNP23</td><td rowspan="3">rs404378206</td><td>GG</td><td>98 (42)</td><td>79 (34)</td><td>1.00</td><td rowspan="3">0.0042</td><td>G</td><td>0.99</td><td>0.9</td></tr><tr><td>GA</td><td>2 (1)</td><td>21 (9)</td><td>11.12<br>(1.34–92.14)</td><td>A</td><td>0.01</td><td>0.1</td></tr><tr><td>–</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP24</td><td rowspan="3">rs429718221</td><td>TT</td><td>16 (7)</td><td>9 (4)</td><td>0.44 (0.11–1.83)</td><td rowspan="3">0.51</td><td>T</td><td>0.42</td><td>0.34</td></tr><tr><td>TC</td><td>51 (22)</td><td>49 (21)</td><td>0.74 (0.30–1.86)</td><td>C</td><td>0.58</td><td>0.66</td></tr><tr><td>CC</td><td>33 (14)</td><td>42 (18)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP25</td><td rowspan="3">rs416266900</td><td>AA</td><td>5 (2)</td><td>0.00</td><td>n/a</td><td rowspan="3">0.14</td><td>A</td><td>0.23</td><td>0.14</td></tr><tr><td>AC</td><td>37 (16)</td><td>28 (12)</td><td>0.60 (0.24–1.51)</td><td>C</td><td>0.77</td><td>0.86</td></tr><tr><td>CC</td><td>58 (25)</td><td>72 (31)</td><td>1.00</td><td></td><td></td><td></td></tr><tr><td rowspan="3">SNP26</td><td rowspan="3">rs602330706</td><td>CC</td><td>95 (41)</td><td>100 (43)</td><td>1.00</td><td rowspan="3">0.093</td><td>C</td><td>0.98</td><td>1</td></tr><tr><td>CT</td><td>5 (2)</td><td>0.00</td><td>n/a</td><td>T</td><td>0.02</td><td>0.00</td></tr><tr><td>–</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The most frequently evaluated SNPs in many ewes' breeds in previous study (rs406779174, rs430181568, rs407388227) for the association with seasonal reproduction were detected in our study. The SNP rs430181568 and rs407388227 (SNP18 and SNP20 in our study, respectively) were linked (D' = 0.99 and r² = 0.94), so that they can be considered as a single marker. Figure 1 shows the linkage disequilibrium plot among all the SNPs in MTNR1A gene detected in our study.

<figure><img src="image_5.png" /><figcaption>Figure 1. The linkage disequilibrium plot (p value) among all the SNPs in MTNR1A gene sequence found in our study.</figcaption></figure>

In the current study we found that the G allele was the most frequent in rs430181568 (SNP18) and rs407388227 (SNP20) loci in the two studied breeds (0.83 and 0.81 for (B) and (QFO), respectively), and consequently, the genotype G/G (70%) was the most frequent. Instead, for rs406779174, C (0.63) and T (0.67) alleles were the most frequent in (B) and (QFO) breeds, respectively (Table 2). Consequently, C/C (44%) and T/T were the most frequent genotypes in (B) and (QFO) breeds, respectively.

Of the possible 2²⁶ haplotypes, 15 haplotypes were found to be common and were included in the subsequent analysis. Haplotype names and frequencies are shown in Only eight haplotypes had a frequency greater than 5% and only five of them were significant (*p* < 0.05) (Supplementary Table S1).
Each SNP was analyzed to determine its relationship with DTL, fertility rate, litter size, and lambs' birth weights in both ewes' breeds. SNP17 (rs406779174), SNP18 (rs430181568), and SNP20 (rs407388227) were associated with DTL in the second year (2018/2019) (*p* < 0.05) but not in the first year (2017/208) of the study (Supplementary Tables S3 and S5). Furthermore, the ewes with G/G and G/A genotypes at SNP18 and SNP20 had a shorter DTL (*p* < 0.05) (Supplementary Table S5). Ewes with TT and CC genotypes at SNP17 had a shorter DTL (*p* < 0.05) (Supplementary Table S3).

Additionally, SNP18 (rs430181568) and SNP20 (rs407388227) were associated with lambs' birth weights in 2019 (*p* < 0.05) in all ewes from both breeds; ewes with G/G and A/A genotypes at SNP18 and G/A and A/A genotypes at SNP20 gave birth to lambs with a better birth weight (*p* < 0.05) (Supplementary Table S6). SNP18 and SNP20 were not associated with fertility rate and litter size in all ewes from both breeds in both years (Supplementary Table S4 and S5). However, SNP17 (rs406779174) was associated with litter size in the second year of the study (2018/2019) (*p* < 0.05) (Supplementary Table S3); it was also associated with lambs' birth weights in both 2018 and 2019 (*p* < 0.05) (Supplementary Table S7).

In the (B) ewes breed, the results from the statistical analysis indicated that SNP18 (rs430181568) and SNP20 (rs407388227) were associated with DTL in the first year of the study (2017/2018) (*p* < 0.05) but not in the second year (Supplementary Tables S8 and S9). Ewes with G/G and G/A at SNP18 (rs430181568) and SNP20 (rs407388227) positions had a shorter DTL (*p* < 0.05) (Supplementary Table S8). There was also an association between these two polymorphisms and lambs' birth weights in 2019 (*p* < 0.05). Furthermore, the ewes with the G/A and A/A genotypes gave birth to lambs with a better weight compared to the ewes with G/G genotypes for this breed (*p* < 0.05) (Supplementary Table S16). However, SNP17 was not associated with any reproductive traits in the (B) ewes (Supplementary Tables S10 and S11).

In the ewes of the (QFO) breed, the values of DTL, fertility rate, and litter size parameters in both years were not affected by the SNP18 and SNP20 (Supplementary Tables S12 and S13). In addition, there was no association between SNP18, SNP20, and lambs' birth weights (Supplementary Table S16). SNP17 was only associated with litter size in the second year of the study (2018/2019) (Supplementary Table S15).

**4. Discussion**

The Barbarine and Queue Fine de l'ouest are the dominant breeds in Tunisia and are raised specially for meat production. They are characterized by a metabolic and digestive adaptation to the contrasting environmental conditions. However, a slow genetic improvement has been shown in sheep flocks of these two breeds. Melatonin, which is normally generated in the brain but can also be secreted by granulosa cells, is a known modulator of follicle development, oocyte maturation, and embryo development [25,26]. In sheep, ewes treated by melatonin displayed better follicle development in culture than those from untreated ewes. Furthermore, melatonin modulates gene expression related to steroidogenesis, differentiation, and granulosa cell luteinization [27]. This hormone protects the integrity of oocytes and granulosa cells via its specific receptors *MTRN1A* and *MTRN1B2* by scavenging reactive oxygen species (ROS) and regulating apoptosis-related genes to prevent apoptosis [28,29]. Our recent study showed that silencing of the *MTNR1A* receptor may be involved in the melatonin antiproliferative effect [30]. Another study reported that *MTNR1A* receptor silencing stimulates follicular atresia by increasing the expression of apoptosis-related genes [31] and reducing progesterone production.

The sequencing analysis of the melatonin receptor 1 gene (*MTNR1A*) in (B) and QFO ewes' breeds revealed 26 SNPs located in the promoter region and exon II. All the SNPs had been already reported in different breeds [12,13,16,20], except SNP rs602330706 in exon II (SNP 26 in our study), which is, to the best of our knowledge, a novel SNP detected for the first time in our study only in the (B) breed. When compared to the Sarda breed and Awassi ewes, only three polymorphic sites causing amino acid changes were identified among our (B) and QFO ewes [12]. The exon II of the *MTNR1A* gene is the
most studied because, within its nucleotide sequence, two specific mutations were found at positions 606 and 612 (SNP17 and SNP18 in our study, respectively) that have been associated with seasonality traits. However, these SNPs are synonymous mutations and therefore not causative in the Rasa Aragonesa breed [16], so it was difficult to explain how they influence the reproductive seasonality [13]. The analysis of the (B) and QFO exon II sequence showed that the SNP18 (old 612/MnII; rs430181568) was associated with the SNP20 (rs407388227), which was coherent with that found in other Mediterranean and European breeds [12,32,33]. This finding suggests that one or both SNPs may be implicated in the reproduction seasonality [13,32]. This amino acid change (Val > Ile) occurred in the fifth transmembrane domain (TM5), in position 220 of the protein chain [34]. TM5 is an essential site for *MTNR1A* function; amino acid change in this domain could be the origin of receptor signaling modification. Furthermore, other amino acid changes in the *MTNR1A* protein sequence at 195, 208, and 211 were considered as responsible for changing the binding capacity of the melatonin receptor. In addition, differences in cAMP inhibition were detected between Val 220 and Ile 220 at the *MTRN1A* receptor; this finding led us to suggest a potential modification in the molecular signal transduction of the melatonin signal in several sheep breeds with different genotypes at SNP20 [35]. Variations in melatonin signal perception may participate in the significant differences in the reproductive traits in our study.

In our study, SNP17, the most significative polymorphism, was not associated with DTL in both breeds, which is contradictory with findings in the Istrian Pramenka and Rasa Aragonesa breeds [13,15]. However, SNP18 (rs430181568) was associated with DTL for the first year studied only in Barbarine ewes. These results are in disagreement with those of Martinez-Royo et al. [36], where SNP17 was associated with DTL but not SNP18 (rs430181568) in the Rasa Aragonesa sheep breed. Our results also confirmed that SNP18 and SNP20 (rs407388227) are fully linked, and these data are consistent with the findings of several scientists [12,13,20,32]. In the present study, some reproductive traits in the Barbarine breed were affected by SNP18 (rs430181568) and SNP20 (rs407388227). Interestingly, there were no effects on the reproductive functions for the three genotypes in the (QFO) ewes. In the (QFO) breed, no significative effect has been identified for the reproductive traits such as DTL, fertility rate, and litter size at SNP18 (rs430181568) and SNP20 (rs407388227) in the studied period. However, SNP17 was associated with litter size in QFO ewes in the second year of the study. Furthermore, lambs' birth weights were affected by SNP18 (rs430181568) and SNP20 (rs407388227) in the (B) breed only; ewes with the G/A and A/A genotypes gave birth to lambs with better weights (p < 0.05).

In fact, Barbarine ewes carrying genotype A/A delayed the onset of the reproductive activity compared to the G/G and G/A genotypes. In both breeds, the different DTL recorded in the two years in the three genotypes may be associated with a different sensitivity to the photoperiod. Considering the mating trend in both the years, Barbarine ewes carrying the G/G and G/A phenotypes mated in the first 30 days after ram introduction; however, ewes with the A/A genotype have the most mating within 40 days from the ram introduction. In the (QFO) ewes, the highest number of the ewes carrying the G/G and G/A genotypes mated in the first 40 days after ram introduction. Therefore, the *MTRN1A* receptor seems to influence reproductive activity in sheep, so the different polymorphisms in SNP18 and SNP20 may modulate the melatonin signal transduction, and this could change the ewes' reproductive response.

The different *MTNR1A* genotypes seem to affect the melatonin signal transmission and seasonality, so we can suggest that G/G in (B) ewes was less sensitive to the photoperiod. Their reproductive activity resumption was advanced compared to A/A ewes in both the years. In fact, sheep carrying the G/G and G/A genotypes responded earlier to the ram effect as if they were in a less deep anestrous compared to the A/A genotype, so A/A ewes were surely in a deeper anestrous state. However, different results were reported in QFO ewes for both years. For the first year, the A/A genotypes were advanced in reproductive recovery compared to the G/G and G/A genotypes, but for the second year, the G/A
ewes were less sensitive to the photoperiod compared to the G/A and A/A ewes. This means that the highest number of G/G and G/A ewes suggesting that their anestrous state was not as deep or as shallow than in A/A ewes. Instead, ewes with the A/A genotype showed the mating peak later than the other genotypes so they need a longer stimulus to onset the reproductive activity. In fact, the depth of anestrous is the most important factor that affects the response of ewes to the ram introduction [37]; DTL length may depend also on the follicle development. In the (B) and (QFO) breeds, litter size and fertility rate were not associated with SNP18 (rs430181568) and SNP20 (rs407388227) in the *MTNR1A* gene sequence. SNP18 and SNP20 were associated with the body weight of lambs at birth for Barbarine ewes in the last year of the study but not with the (QFO) breed. In fact, Barbarine ewes with the A/A genotype recorded the highest body weight of lambs at birth (*p* < 0.05). This finding was also compatible with the results of Mahdi et al., who reported that mutations located in the *MTNR1A* gene sequence could affect lambs' weights at birth for the Awassi breed [38]. The effect of SNP18 and SNP20 in the (QFO) breed is considered marginal, which is in agreement with what is described by Cosso et al. and Hernandez et al. in the Awassi and Ile de France breeds, respectively [18,23].

In the present study, Barbarine ewes with G/G and G/A showed readiness to resume the reproductive activity earlier than the QFO ewes at the ram introduction. Our different results, from the first and the second year of the study, let us suggest that sheep breed, ewes' age and environmental factors may be the most important triggers that influence the ewes' response to the photoperiod; this fact could explain the lack of effect of these two mutations (rs430181568 and rs407388227) on QFO ewes' reproductive performance parameters.

**5. Conclusions**

Our data confirm the presence of polymorphisms in the *MTNR1A* gene sequence in the Barbarine and Queue Fine de l'Ouest breeds. Twenty-six polymorphic sites were detected, three of which caused amino acid changes. The total linkage of SNP18 and SNP20, which are considered as a single marker, revealed that these polymorphisms can influence the reproductive activity resumption in Barbarine ewes with the G/G and G/A genotypes but not in QFO ewes. These results might be considered in sheep selection programs for reproductive genetic improvement. To explain the association more precisely between polymorphisms in the *MTNR1A* gene sequence and the reproductive traits, other studies with a higher number of ewes from both breeds need to be conducted. In addition, it would be useful to study the effect of the several genotypes at SNP18 and SNP20 on ovarian cells to demystify the intricate role of the *MTNR1A* gene at the ovary.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13030448/s1, Supplementary Table S1. Haplotypes and their frequency in all the ewes; Supplementary Table S2. Results of the first year in the studied ewes from both breeds: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S3. Results of the second year in the studied ewes from both breeds: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S4. Results of the first year in the studied ewes from both breeds: DTL, fertility rate, and litter size for each genotype at SNP18/SNP20; Supplementary Table S5. Results of the second year in the studied ewes from both breeds: DTL, fertility rate, and litter size for each genotype at SNP18/SNP20; Supplementary Table S6. Results of lambs' birth weights for the three years in the studied ewes from both breeds; Supplementary Table S7. Results of lambs' birth weights for the three years in the studied ewes from both breeds; Supplementary Table S8. Results of the first year in the studied Barbarine ewes: DTL, fertility rate, litter size, and lambs' birth weights for each genotype at SNP18/SNP20; Supplementary Table S9. Results of the second year in the studied Barbarine ewes: DTL, fertility rate, litter size, and lambs' birth weights for each genotype at SNP18/SNP20; Supplementary Table S10. Results of the first year in the studied ewes from Barbarine breed: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S11. Results of the second year in the studied ewes from *Barbarine breed*: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S12. Results of the second year in the studied Queue Fine de l'Ouest ewes: DTL, fertility rate, litter size, and lambs' birth weights
for each genotype at SNP18/SNP20; Supplementary Table S13. Results of the second year in the studied Queue Fine de l'Ouest ewes: DTL, fertility rate, litter size, and lambs' birth weights for each genotype at SNP18/SNP20; Supplementary Table S14. Results of the first year in the studied ewes from QFO breed: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S15. Results of the second year in the studied ewes from QFO breed: DTL, fertility rate, and litter size for each genotype at SNP17; Supplementary Table S16. Results of lambs' birth weights for the three years in the studied Barbarine and QFO ewes breed; Supplementary Table S17. Primer sequences information and PCRs.

**Author Contributions:** Conceptualization, M.A. and A.A.; methodology, H.C. and R.J.; software, A.B.A.; validation, M.A., S.A.M. and B.H.; formal analysis, A.A.; investigation, A.A.; resources, A.A.; data curation, A.A.; writing—original draft preparation, A.A.; writing—review and editing, A.A., S.A.M., M.A. and M.-A.S.; visualization, A.A.; supervision, M.A., M.-A.S. and B.H.; project administration, T.N.; funding acquisition, M.-A.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by NSERC # 2014-94532 Canada grant to MASirard.

**Institutional Review Board Statement:** The animal study protocol was approved by the Ethics Committee of the National School of Veterinary Medicine (CEEA-ENMV 57/2022) for studies involving animals.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author. The data are not publicly available to preserve privacy of the data.

**Acknowledgments:** We would like to acknowledge the support of the Livestock and Pastures office of Tunisia in addition to the financial support of NSERC Canada.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Anderson, G.M.; Hardy, S.L.; Valent, M.; Billings, H.J.; Connors, J.M.; Goodman, R.L. Evidence that Thyroid Hormones Act in the Ventromedial Preoptic Area and the Premammillary Region of the Brain to Allow the Termination of the Breeding Season in the Ewe. *Endocrinology* **2003**, *144*, 2892-2901. [CrossRef] [PubMed]

2. Malpaux, B.; Viguié, C.; Thiéry, J.C.; Chemineau, P. Contrôle photopériodique de la reproduction. *INRA Prod. Anim.* **1996**, *9*, 9-23. [CrossRef]

3. Belkadi, S.; Safsaf, B.; Heleili, N.; Tlidjane, M.; Belkacem, L.; Oucheriah, Y. Seasonal influence on sperm parameters, scrotal measurements, and serum testosterone in Ouled Djellal breed rams in Algeria. *Vet. World* **2017**, *10*, 1486. [CrossRef] [PubMed]

4. Karsch, F.J.; Bittman, E.L.; Foster, D.L.; Goodman, R.L.; Legan, S.J.; Robinson, J.E. Neuroendocrine Basis of Seasonal Reproduction A2-GREEP, ROY O. In Proceedings of the 1983 Laurentian Hormone Conference, Mont-Tremblant, QC, Canada, 1 February 1983; Academic Press: Boston, MA, USA, 1984; Volume 40, pp. 185-232.

5. Balasse, M.; Keller, M.; Balasescu, A.; Tresset, A.; Chemineau, P. *Cycle de Reproduction des Ovins et Rythme Saisonnier de L'élevage: Un Schéma Plurimillénaire Bouleversé par la Recherche Agronomique*; Editions Des Archives Contemporaines: Paris, France, 2017; pp. 25-50.

6. Ortavant, R.; Pelletier, J.; Ravault, J.P.; Thimonier, J.; Volland-Nail, P. Photoperiod: Main proximal and distal factor of the circannual cycle of reproduction in farm mammals. *Oxf. Rev. Reprod.* **1985**, *7*, 305-345.

7. Hastings, M.H.; Herbert, J.; Martensz, N.D.; Roberts, A.C. Melatonin and the brain in photoperiodic mammals. *Ciba Found Symp.* **1985**, *117*, 57-77. [CrossRef]

8. Carcangiu, V.; Mura, C.; Parmeggiani, A.; Piccione, G.; Bini, P.; Cosso, G.; Luridiana, S. Daily rhythm of blood melatonin concentrations in sheep of different ages. *Biol. Rhythm Res.* **2013**, *44*, 908-915. [CrossRef]

9. Sliwowska, J.H.; Billings, H.J.; Goodman, R.L.; Coolen, L.M.; Lehman, M.N. The Premammillary Hypothalamic Area of the Ewe: Anatomical Characterization of a Melatonin Target Area Mediating Seasonal Reproduction. *Biol. Reprod.* **2004**, *70*, 1768-1775. [CrossRef]

10. Weaver, D.R.; Liu, C.; Reppert, S.M. Nature's knockout: The Mel1b receptor is not necessary for reproductive and circadian responses to melatonin in Siberian hamsters. *Mol. Endocrinol.* **1996**, *10*, 1478-1487. [CrossRef]

11. Dubocovich, M.L.; Rivera-Bermudez, M.A.; Gerdin, M.J.; Masana, M.I. Molecular pharmacology, regulation and function of mammalian melatonin receptors. *Front. Biosci.-Landmark* **2003**, *8*, 1093-1108. [CrossRef]

12. Luridiana, S.; Cosso, G.; Pulinas, L.; Di Stefano, M.V.; Curone, G.; Carcangiu, V.; Mura, M.C. New polymorphisms at MTNR1A gene and their association with reproductive resumption in sarda breed sheep. *Theriogenology* **2020**, *158*, 438-444. [CrossRef]
13. Starič, J.; Farci, F.; Luridiana, S.; Mura, M.C.; Pulinas, L.; Cosso, G.; Carcangiu, V. Reproductive performance in three Slovenian sheep breeds with different alleles for the MTNR1A gene. *Anim. Reprod. Sci.* **2020**, *216*, 106352. [CrossRef] [PubMed]

14. Carcangiu, V.; Luridiana, S.; Vacca, G.M.; Daga, C.; Mura, M.C. A polymorphism at the melatonin receptor 1A (<i>MTNR1A</i>/*i*) gene in Sarda ewes affects fertility after AI in the spring. *Reprod. Fertil. Dev.* **2011**, *23*, 376–380. [CrossRef] [PubMed]

15. Martínez-Royo, A.; Lahoz, B.; Alabart, J.L.; Folch, J.; Calvo, J.H. Characterisation of the Melatonin Receptor 1A (MTNR1A) gene in the Rasa Aragonesa sheep breed: Association with reproductive seasonality. *Anim. Reprod. Sci.* **2012**, *133*, 169–175. [CrossRef] [PubMed]

16. Calvo, J.H.; Serrano, M.; Martinez-Royo, A.; Lahoz, B.; Sarto, P.; Ibañez-Deler, A.; Folch, J.; Alabart, J.L. SNP rs403212791 in exon 2 of the <i>em</i>*MTNR1A*/*em* gene is associated with reproductive seasonality in the Rasa aragonesa sheep breed. *Theriogenology* **2018**, *113*, 63–72. [CrossRef] [PubMed]

17. Mura, C.; Luridiana, S.; Bodano, S.; Daga, C.; Cosso, G.; Diaz, M.L.; Bini, P.P.; Carcangiu, V. Influence of melatonin receptor 1A gene polymorphisms on seasonal reproduction in Sarda ewes with different body condition scores and ages. *Anim. Reprod. Sci.* **2014**, *149*, 173–177. [CrossRef] [PubMed]

18. Cosso, G.; Nehme, M.; Luridiana, S.; Pulinas, L.; Curone, G.; Hosri, C.; Carcangiu, V.; Mura, M.C. Detection of Polymorphisms in the MTNR1A Gene and Their Association with Reproductive Performance in Awassi Ewes. *Animals* **2021**, *11*, 583. [CrossRef] [PubMed]

19. He, X.; Zhang, Z.; Chu, M. The effect of SNP rs400827589 in exon 2 of the MTNR1B gene on reproductive seasonality and litter size in sheep. *Vet. Med. Sci.* **2020**, *2020*, 804–812. [CrossRef] [PubMed]

20. Pulinas, L.; Starič, J.; Cosso, G.; Curone, G.; Mura, M.C.; Carcangiu, V.; Luridiana, S. MTNR1A gene polymorphisms and reproductive recovery after seasonal anoestrus in different Mediterranean sheep breeds. *Anim. Reprod. Sci.* **2022**, *236*, 106905. [CrossRef]

21. Meena, A.; Kumar, R.; Kumari, R.; Jyotsana, B.; Prince, L.; Kumar, S. Genetic polymorphism of melatonin receptor 1A (MTNR1A) gene in Indian sheep breeds. *Indian J. Anim. Sci.* **2013**, *83*, 1284–1287.

22. Carcangiu, V.; Mura, M.C.; Vacca, G.M.; Pazzola, M.; Dettori, M.L.; Luridiana, S.; Bini, P.P. Polymorphism of the melatonin receptor MT1 gene and its relationship with seasonal reproductive activity in the Sarda sheep breed. *Anim. Reprod. Sci.* **2009**, *116*, 65–72. [CrossRef]

23. Hernandez, X.; Bodin, L.; Chesneau, D.; Guillaume, D.; Chemineau, P.; Malpaux, B.; Migaud, M. Relationship between MT1 melatonin receptor gene polymorphism and seasonal physiological responses in Île-de-France ewes. *Reprod. Nutr. Dev.* **2005**, *45*, 151–162. [CrossRef] [PubMed]

24. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

25. Suofu, Y.; Li, W.; Jean-Alphonse, F.; Jia, J.; Khattar, N.; Li, J.; Baranov, S.; Leronni, D.; Mihalik, A.; He, Y.; et al. Dual role of mitochondria in producing melatonin and driving GPCR signaling to block cytochrome c release. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, E7997–E8006. [CrossRef]

26. Slominski, R.M.; Reiter, R.J.; Schlabritz-Loutsevitch, N.; Ostrom, R.S.; Slominski, A.T. Melatonin membrane receptors in peripheral tissues: Distribution and functions. *Mol. Cell. Endocrinol.* **2012**, *351*, 152–166. [CrossRef] [PubMed]

27. Arjoune, A.; Sirard, M.-A. Melatonin Signaling Pathways Implicated in Metabolic Processes in Human Granulosa Cells (KGN). *Int. J. Mol. Sci.* **2022**, *23*, 2988.

28. Cruz, M.H.C.; Leal, C.L.V.; Cruz, J.F.; Tan, D.X.; Reiter, R.J. Essential actions of melatonin in protecting the ovary from oxidative damage. *Theriogenology* **2014**, *82*, 925–932. [CrossRef] [PubMed]

29. Tanabe, M.; Tamura, H.; Taketani, T.; Okada, M.; Lee, L.; Tamura, I.; Maekawa, R.; Asada, H.; Yamagata, Y.; Sugino, N. Melatonin protects the integrity of granulosa cells by reducing oxidative stress in nuclei, mitochondria, and plasma membranes in mice. *J. Reprod. Dev.* **2015**, *61*, 35–41. [CrossRef] [PubMed]

30. Arjoune, A.; Sirard, M.-A. The genomic response of human granulosa cells (KGN) to melatonin and specific agonists/antagonists to the melatonin receptors. *Sci. Rep.* **2022**, *12*, 17539. [CrossRef]

31. Wang, S.-J.; Liu, W.-J.; Wang, L.-K.; Pang, X.-S.; Yang, L.-G. The role of Melatonin receptor MTNR1A in the action of Melatonin on bovine granulosa cells. *Mol. Reprod. Dev.* **2017**, *84*, 1140–1154. [CrossRef]

32. Mura, M.C.; Luridiana, S.; Pulinas, L.; Di Stefano, M.V.; Carcangiu, V. Reproductive response to male joining with ewes with different allelic variants of the MTNR1A gene. *Anim. Reprod. Sci.* **2019**, *200*, 67–74. [CrossRef] [PubMed]

33. Saxena, V.; Jha, B.; Meena, A.; Narula, H.; Kumar, D.; Naqvi, S. Assessment of Genetic Variability in the Coding Sequence of Melatonin Receptor Gene (MTNR1A) in Tropical Arid Sheep Breeds of India. *Reprod. Domest. Anim.* **2015**, *50*, 517–521. [CrossRef] [PubMed]

34. Barrett, P.; Schuster, C.; Mercer, J.; Morgan, P.J. Sensitization: A Mechanism for Melatonin Action in the Pars Tuberalis. *J. Neuroendocrinol.* **2003**, *15*, 415–421. [CrossRef] [PubMed]

35. Trecherel, E.; Batailler, M.; Chesneau, D.; Delagrange, P.; Malpaux, B.; Chemineau, P.; Migaud, M. Functional characterization of polymorphic variants for ovine MT1 melatonin receptors: Possible implication for seasonal reproduction in sheep. *Anim. Reprod. Sci.* **2010**, *122*, 328–334. [CrossRef]

36. Martinez-Royo, A.; Alabart, J.L.; Sarto, P.; Serrano, M.; Lahoz, B.; Folch, J.; Calvo, J.H. Genome-wide association studies for reproductive seasonality traits in Rasa Aragonesa sheep breed. *Theriogenology* **2012**, *99*, 21–29. [CrossRef]
37. Rosa, H.J.D.; Bryant, M.J. Seasonality of reproduction in sheep. *Small Rumin. Res.* **2003**, *48*, 155–171. [CrossRef]

38. Salih Mahdi, A.; Al-Khafaji, M.; Al-Wazeer, A.; Al-Khuzai, H. Combined effects of melatonin receptor gene polymorphism and vitamin E supplementation on litter size and lambs body weight in Awassi ewes. *Ann. Trop. Med. Public Health* **2020**, *23*, 64–69.
[CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.