Genome- wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi- arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunis ia n sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of Europe an gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58–49.61 Mbp). These results were confirmed using FST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds

Genetic structure of Tunisian sheep breeds as inferred from genome-wide SNP markers

Assessing the status of genetic variability of native sheep breeds could provide important clues for research and policy makers to devise better strategies for the conservation and management of genetic resources. In this study, a genetic investigation of Tunisian sheep breeds using a genome-wide scan of approximately 50,000 SNPs was performed. To reconstruct genetic structure and relationships among four sheep breeds, 40 samples belonging to fat-tailed Barbarine, Queue Fine de l'Ouest, Noire de Thibar and D'Man breeds were genotyped using Illumina Ovine SNP50 BeadChip. Tunisian breeds averaged 96 % polymorphic loci with an expected heterozygosity (He = 0.36). Genetic analysis of relationship between breeds using Bayesian clustering, MDS and Neighbor-Network analysis, and estimation of FST genetic structure, highlighted the genetic differentiation of Noire de Thibar breed from the other local breeds, reflecting the effect of past events of introgression of European gene pool. The Queue Fine de l'Ouest breed showed a genetic heterogeneity and was close to Barbarine and D'Man breeds, as evidenced by MDS and the lowest level of differentiation with Barbarine breed (FST = 1.8 %). The D'Man breed shared a considerable gene flow with the thin-tailed Queue Fine de l'Ouest breed. Possible factors explaining the genetic patterns observed, such as considerable gene flow probably due to anthropogenic activities in the light of population management and conservation programs

Genetic parameters of somatic cell scores using random regression test-day models with Legendre polynomials in Tunisian dairy cattle

9 Pág.The aim of this study was to analyse the genetic parameters of test-day somatic cell scores (TDSCS) in Tunisian dairy cattle using random regression models (RRM) with homogeneous residual variances. Data included 43,647 test-day somatic cell count (TDSCC) records collected on 4825 Holstein-Friesian cows with parities up to the fifth, and calving dates between 2000 and 2014. Records of the same animal belonging to different lactations were treated as repeated records. Additive genetic (AG) and permanent environmental (PE) random effects were fitted using Legendre polynomials functions (LP) up to the fifth-order. The models included fixed regressions on days-in-milk (DIM), modelling the average phenotypic curve using Ali-Schaeffer's lactation function. Overall, three models (LP3, LP4 and LP5) where compared using the -2log marginal density for Bayes factor (-2logp) and the deviance information criterion (DIC). Models were also compared by similarities between the estimated breeding values and genetic trends. For all models, the total variances tended to decrease during the first month of the lactation, then they maintained a nearly constant trend. All models showed the largest AG and PE variances at the beginning and at the end of lactation with an oscillating pattern in the middle. The heritability estimates of from the three models ranged from 0.01 to 0.08 with the lowest heritability estimates being in the mid-lactation period. Genetic correlations between adjacent DIM were high and then decreased as the intervals widened. All the best statistical criteria were obtained with the LP5 model, however, no differences were observed between the estimated 305 average breeding values (EBV305) from the three models. This assumes that the LP3 model can be selected as a parsimonious model since it contains the smallest number of parameters. This RRM study provided a good insight into the changes in genetic merit of the somatic cell score trait in Tunisian dairy cattle and that an EBV for the SCS should be included in the selection programs. It also showed the benefits given of using RRMs to perform national genetic evaluations of dairy cattle for milk production traits in Tunisia.Peer reviewe

Random regression models to estimate genetic parameters for milk yield, fat, and protein contents in Tunisian Holsteins

22 Pág.This study aimed to find the parsimonious random regression model (RRM) to evaluate the genetic potential for milk yield (MY), fat content (FC), and protein content (PC) in Tunisian Holstein cows. For this purpose, 551,139; 331,654; and 302,396 test day records for MY, FC, and PC were analysed using various RRMs with different Legendre polynomials (LP) orders on additive genetic (AG) and permanent environmental (PE) effects, and different types of residual variances (RV). The statistical analysis was performed in a Bayesian framework with Gibbs sampling, and the model performances were assessed, mainly, on the predictive ability criteria. The study found that the optimal model for evaluating these traits was an RRM with a third LP order and nine classes of heterogeneous RV. In addition, the study found that heritability estimates for MY, FC, and PC ranged from 0.11 to 0.22, 0.11 to 0.17, and 0.12 to 0.18, respectively, indicating that genetic improvement should be accompanied by improvements in the production environment. The study also suggested that new selection rules could be used to modify lactation curves by exploiting the canonical transformation of the random coefficient covariance (RC) matrix or by using the combination of slopes of individual lactation curves and expected daily breeding values.This research has no funding source.Peer reviewe

Sequencing and Characterization of αs2-Casein Gene (CSN1S2) in the Old-World Camels Have Proven Genetic Variations Useful for the Understanding of Species Diversification

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3'-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels

Marker-assisted introgression of fecundity mutation into Barbarine sheep breed under low input production system

Fecundity gene introgression has become an economical method to improve sheep prolificacy in developing countries. The FecXBar and FecGH mutations are variants of the BMP15 and GDF9 sheep genes, respectively, identified in Tunisian Barbarine sheep prolific line “W-INRAT,” created through a prolificacy-based selection program. The first mutation caused increased prolificacy in heterozygous state and sterility in homozygous ewes. The aim of this work was to increase the number of effective carriers by the introgression of fecundity mutation into non-carrier conventional flocks based on a marker-assisted breeding program. The genotyping was carried out to follow up on the segregation of prolificacy mutation. The conventional ewes, raised in state farms of Tunisia and inseminated by “W-INRAT” rams, yielded 100% carriers females at heterozygous state. These females were selected to be inseminated with conventional rams and produced carriers’ lambs with 66.7% males and 33.3% females. The prolific males will be kept in breeding centers to disseminate the fecundity mutation in commercial flocks. This approach facilitates the dissemination of fecundity genes and contributes to livelihood improvement in communities raising Barbarine sheep

CA.RA.VA.N project: toward implementation of a modern camel selection system in Northern Africa?

A project entitled “CA.RA.VA.N., Toward a CAmel tRAnsnational VAlue chain” (http://www.arimnet2.net/index.php/researchprojects/projects-2nd-call-2/ca-ra-va-n), involving research institutions from Spain, Italy, France, Algeria, Morocco, and Tunisia, has been recently financed by EU through the Arimnet2 instrument. Additional partners from eligible (Tunisia) and non-eligible areas (Mauritania and Austria) have been associated to the project due to relevant interest in, and expected contribution to, the topic of the project. General goals and activities planned within the CA.RA.VA.N. project, together with an overview of the already launched initiatives, known constraints and emerging opportunities for the development of the dromedary sector in Northern African countries will be discussed

Social network analysis of the stakeholders involved in the dromedary sector in the Mediterranean region

Marginal arid zones in the south Mediterranean are faced with the dramatic departure of their labor forces through migration. Interest in the capacity of the dromedary species to enhance desert ecosystems and to be a potential lever of economic development of these marginal zones only started to grow in the last two decades. Based on an empirical survey of 179 stakeholders in four Mediterranean countries, we explored the links of the stakeholders in the dromedary sector in two dimensions: Horizontal links with peers and vertical links along the value chain and in resource management. Both descriptive statistics and social network analysis highlight the original organization of the dromedary sector around herders and their social and cultural organization at the territorial level. Therefore, even if milk production and processing start to constitute an opportunity for the young generation who do not necessarily have the financial capacity to invest in a large dromedary herd, this change towards milk valorization can only happen if it is linked with the traditional system based on mobility. Using a systemic approach and working toward multiple valorizations of dromedary products instead of only targeting milk productivity should be explored