Search

Background Availability of single nucleotide polymorphism (SNP) genotyping arrays and progress in statistical analyses have allowed the identification of genomic regions and genes under selection in chicken. In this study, SNP data from the 600 K Affymetrix chicken array were used to detect signatures of selection in 23 local Italian chicken populations. The populations were categorized into four groups for comparative analysis based on live weight (heavy vs light) and geographical area (Northern vs Southern Italy). Putative signatures of selection were investigated by combining three extended haplotype homozygosity (EHH) statistical approaches to quantify excess of haplotype homozygosity within (iHS) and between (Rsb and XP-EHH) groups. Presence of runs of homozygosity (ROH) islands was also analysed for each group. Results After editing, 541 animals and 313,508 SNPs were available for statistical analyses. In total, 15 candidate genomic regions that are potentially under selection were detected among the four groups: eight within a group by iHS and seven by combining the results of Rsb and XP-EHH, which revealed divergent selection between the groups. The largest overlap between genomic regions identified to be under selection by the three approaches was on chicken chromosome 8. Twenty-one genomic regions were identified with the ROH approach but none of these overlapped with regions identified with the three EHH-derived statistics. Some of the identified regions under selection contained candidate genes with biological functions related to environmental stress, immune responses, and disease resistance, which indicate local adaptation of these chicken populations. Conclusions Compared to commercial lines, local populations are predominantly reared as backyard chickens, and thus, may have developed stronger resistance to environmental challenges. Our results indicate that selection can play an important role in shaping signatures of selection in local chicken populations and can be a starting point to identify gene mutations that could have a useful role with respect to climate change.

About three decades of breeding and selection in the Valle del Belìce sheep are expected to have left several genomic footprints related to milk production traits. In this study, we have assembled a dataset with 451 individuals of the Valle del Belìce sheep breed: 184 animals that underwent directional selection for milk production and 267 unselected animals, genotyped for 40,660 single-nucleotide polymorphisms (SNPs). Three different statistical approaches, both within (iHS and ROH) and between (Rsb) groups, were used to identify genomic regions potentially under selection. Population structure analyses separated all individuals according to their belonging to the two groups. A total of four genomic regions on two chromosomes were jointly identified by at least two statistical approaches. Several candidate genes for milk production were identified, corroborating the polygenic nature of this trait and which may provide clues to potential new selection targets. We also found candidate genes for growth and reproductive traits. Overall, the identified genes may explain the effect of selection to improve the performances related to milk production traits in the breed. Further studies using high-density array data, would be particularly relevant to refine and validate these results.

Admixture is an evolutionary process that enables short-term adaptation. The Creole cattle from Guadeloupe is a tropically adapted breed. Its three-way admixture and long-term isolation offer a unique opportunity for understanding the genetic determinants of adaptive admixture in livestock. Here, we sequenced 23 Creole cattle from Guadeloupe (GUA) and combined our data with sequenced genomes of 99 cattle from 25 breeds representative of European, African and indicine groups to provide the most detailed exploration, to date, of patterns of genetic variation and to detect selection signatures in this population. We detect 17 228 983 single nucleotide polymorphisms (SNPs) and we confirmed the higher level of African and indicine ancestries, compared to the European ancestry, in the GUA genome. We show that, unlike Criollo cattle, GUA population originates directly from West Africa with indicine ancestry inherited via West African ancestors. We relied on consistency of signals across various methods based on excess of haplotype homozygosity, differences in allele frequencies and excess ⁄ deficiency of local ancestry to identify five strong candidate regions showing an excess of indicine ancestry. These encompass immune-, heat-tolerance- and physical exercise-related genes. Moreover, we found that a previously identified horn-related gene, RXFP2 is under strong selective pressure in GUA genome likely owing to human-driven (socio-cultural) pressure. Our study highlights the role played by population admixture for driving rapid adaptive response to local environmental constraints.

Background To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. Results The results indicate that a large part of the Merino’s genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. Conclusions To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as “headgear,” which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

International audience; Natural-driven selection is supposed to have left detectable signatures on the genome of North African cattle which are often characterized by the fixation of genetic variants associated with traits under selection pressure and/or an outstanding genetic differentiation with other populations at particular loci. Here, we investigate the population genetic structure and we provide a first outline of potential selection signatures in North African cattle using single nucleotide polymorphism genotyping data. After comparing our data to African, European and indicine cattle populations, we identified 36 genomic regions using three extended haplotype homozygosity statistics and 92 outlier markers based on Bayescan test. The 13 outlier windows detected by at least two approaches, harboured genes (e.g. GH1, ACE, ASIC3, HSPH1, MVD, BCL2, HIGD2A, CBFA2T3) that may be involved in physiological adaptations required to cope with environmental stressors that are typical of the North African area such as infectious diseases, extended drought periods, scarce food supply, oxygen scarcity in the mountainous areas and high-intensity solar radiation. Our data also point to candidate genes involved in transcriptional regulation suggesting that regulatory elements had also a prominent role in North African cattle response to environmental constraints. Our study yields novel insights into the unique adaptive capacity in these endangered populations emphasizing the need for the use of whole genome sequence data to gain a better understanding of the underlying molecular mechanisms.

The availability of genotyping assays has allowed the detailed evaluation of cattle genetic diversity worldwide. However, these comprehensive studies did not include some local European populations, including autochthonous Italian cattle. In this study, we assembled a large-scale, genome-wide dataset of single nucleotide polymorphisms scored in 3,283 individuals from 205 cattle populations worldwide to assess genome-wide autozygosity and understand better the genetic relationships among these populations. We prioritized European cattle, with a special focus on Italian breeds. Moderate differences in estimates of molecular inbreeding calculated from runs of homozygosity (FROH) were observed among domesticated bovid populations from different geographic areas, except for Bali cattle. Our findings indicated that some Italian breeds show the highest estimates of levels of molecular inbreeding among the cattle populations assessed in this study. Patterns of genetic differentiation, shared ancestry, and phylogenetic analysis all suggested the occurrence of gene flow, particularly among populations originating from the same geographical area. For European cattle, we observed a distribution along three main directions, reflecting the known history and formation of the analyzed breeds. The Italian breeds are split into two main groups, based on their historical origin and degree of conservation of ancestral genomic components. The results pinpointed that also Sicilian breeds, much alike Podolian derived-breeds, in the past experienced a similar non-European influence, with African and indicine introgression.

Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case–control genome-wide association study (GWAS), the Bayesian population differentiation analysis (FST) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and FST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (

International audience; Mutations in the FecL locus are associated with large variation in ovulation rate and litter size in the French Lacaune sheep breed. It has been shown that the B4GALNT2 gene within the FecL locus is most likely responsible for the high fecundity in the French breed. In this study, we have highlighted the segregation of the FecLL mutation within the B4GALNT2 gene in North African sheep breeds and notably in the highly prolific D'man breed. Genotyping of a sample of 183 Tunisian D'man individuals revealed a high frequency (0.65) of the prolific allele FecLLwhich was attributed to the adoption of a decades‐old breeding strategy based on the selection of ewe lambs born from large litter size. Homozygous LL ewes showed a significantly increased litter size compared to heterozygous and non‐carrier ewes (FecLL/FecLL = 2.47 ± 0.09 vs. FecLL/FecL+ = 2.23 ± 0.09, p < 0.05 and FecL+/FecL+ = 1.93 ± 0.18, p < 0.01). The presence of the FecLLpolymorphism in both D'man and Lacaune breeds argues for an ancestral origin of this mutation and brings an answer to the old question of the genetic determinism of the extreme prolificacy of the D'man ewes. The results of this study can help to establish planned genotype‐based mating allowing both higher profit for the breeders and an optimal management of the FecLL mutation in D'man sheep populations.

International audience; To carry out effective genome-wide association studies, information about linkage disequilibrium (LD) is essential. Here, we used medium-density SNP chips to provide estimates of LD in native Tunisian cattle. The two measures of LD that were used, mean r(2) and D', decreased from 0.26 to 0.05 and from 0.73 to 0.40, respectively, when the distance between markers increased from less than 20 Kb to 200 Kb. The decay in LD over physical distance occurred at a faster rate than that reported for European and other indigenous breeds, and reached background levels at less than 500 Kb distance. This is consistent with the absence of strong selective pressure within the Tunisian population and suggests that, in order to be effective, any potential genome-wide association mapping studies will need to use chips with higher marker density. An analysis of effective population size (Ne) based on LD data showed a decline in past Ne, with a sudden drop starting about eight generations ago. This finding, combined with the high levels of recent inbreeding revealed by runs of homozygosity (ROH) analysis, indicate that this population is endangered and may be in urgent need of a conservation plan that includes a well-designed genetic management program.

Mutations in the FecL locus are associated with large variation in ovulation rate and litter size in the French Lacaune sheep breed. It has been shown that the B4GALNT2 gene within the FecL locus is most likely responsible for the high fecundity in the French breed. In this study, we have highlighted the segregation of the FecL

International audience; In the present study, we used medium density SNP chips and comparisons with worldwide cattle to explore the phylogenetic history of Tunisian and Iberian cattle. We have pointed out an admixed origin of Tunisian populations and low African taurine proportions in the Iberian cattle genome. Our results suggest that iberian and Tunisian cattle derived from different origins and that a past admixture occurred between them

International audience; Understanding between and within-breeds genetic variability is essential in the choice of conservation management decisions for threatened populations. In this study we assessed the genetic diversity of the Algerian Guelmoise cattle (GUE) by analyzing data on 24 GUE individuals genotyped for the Illumina BovineSNP50 BeadChipv2. We also provided a detailed description of the population structure of GUE using comparisons with 23 worldwide cattle populations, selected as being representative of African, South European and indicine populations, in addition to four North African populations. We show that GUE is an admixed population which has strong genetic similarity to Tunisian cattle owing both to a common ancestral origin and a continuous gene flow. GUE was found to be less subject to foreign breed introgression than native Tunisian cattle. We also found support for a past admixture between GUE and the Moroccan Oulmes Zaer breed. The observed heterozygosity of GUE population was similar to that observed for the South European breeds and the other Maghreb populations (Ho similar to 0.30). By investigating the occurrence of Runs of Homozysosity (ROHs) in a subset of nine populations, we found that GUE had the lowest cumulative ROH length per animal with the largest proportion of their genome covered with medium and long ( > 5 Mb) ROHs. This likely indicates that GUE is derived from a large ancestral population which has subsequently undergone a recent increasing levels of inbreeding. ROH islands within GUE harbored genes enriched in immune defense and transcriptional regulation functions in addition to genes related to the Pregnancy-associated proteins family. Our study is the first of its kind providing insights into the genomic structure and the phylogenetic relationships of Guelmoise cattle population. It points out the possibility to establish a common genetic management program between Algeria and Tunisia aimed at safeguarding the endangered indigenous North African cattle.

Background Korean native chicken (KNC) is a well-known breed due to its superior meat taste. This breed, however, owing to a low growth rate, has a high market price. In order to overcome this disadvantage, the National Institute of Animal Science (NIAS) in Korea developed a commercial KNC breed, named Woorimatdag version 2 (WM2), an upgraded version of the Woorimatdag (WM1) breed and the WM2 was created by crossing the KNC with meat type breeds. This study aims to discriminate between WM2 and other chicken breeds using microsatellite (MS) markers. Methods A total of 302 individuals from eight Korean chicken populations were examined. The genetic diversity and population structure analysis were investigated using Cervus, API-CALC, STRUCTURE, PowerMarker programs. Results Based on heterozygosity and polymorphic information content (PIC) values, 30 MS markers were initially selected from 150 markers. The identified average number of alleles (Na), expected heterozygosity, and PIC values for the WM2 samples were 7.17, 0.741, and 0.682, respectively. Additionally, the paternity of individuals was assigned with a success rate of greater than 99% using 12 markers, the best minimum number of markers. The 12 selected markers contained heterozygosity and PIC values above 0.7 and probability of identity values around zero. Using these markers, the determined probability of identity (PI), PIhalf-sibs, and PIsibs values were 3.23E-33, 5.03E-22, and 8.61E-08, respectively. Conclusions WM2 is well differentiated with respect to other chicken breeds based on estimated genetic distances. The results presented here will contribute to the identification of commercial WM2 chicken in the market. Electronic supplementary material The online version of this article (doi:10.1186/s40781-015-0044-6) contains supplementary material, which is available to authorized users.

Background Tunisian local cattle populations are at risk of extinction as they were massively crossed with imported breeds. Preservation of indigenous livestock populations is important because each of them comprises a unique set of genes resulting from a local environment-driven selection that occurred over hundreds of years. The diversity and genetic structure of Tunisian local cattle populations are poorly understood. However, such information is crucial to the conservation and sustainable use of genetic resources. In addition, comparing the genomic structure of population sets from different parts of the world could help yield insight into their origin and history. In the present study, we provide a detailed assessment of the population structure of the three Tunisian local cattle populations using various methods, and we highlight their origin and history by investigating approximately ~38,000 SNPs in a broad panel of 878 individuals from 37 worldwide cattle breeds representative of African, European and indicine populations. Results Our study revealed a low level of divergence and high genetic diversity in Tunisian local cattle reflecting low levels of genetic drift. A Comparison with the worldwide cattle panel pinpointed the admixed origin of the genome of the three Tunisian populations with the two main European and African ancestries. Our results were in agreement with previous historical and archaeological reports about the past gene flow that existed between North African and South European breeds, in particular with Iberian cattle. We also detected a low-level indicine introgression in the three Tunisian populations and we inferred that indicine ancestry was inherited via African ancestors. Conclusions Our results represent the first study providing genetic evidence about the origin and history of Tunisian local cattle. The information provided by the fine-scale genetic characterization of our study will enhance the establishment of a national conservation strategy for these populations. These results may enable the identification of genetic variants involved in adaptation to harsh environmental conditions. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1638-6) contains supplementary material, which is available to authorized users.

Fertility quantitative trait loci (QTL) are of high interest in dairy cattle since insemination failure has dramatically increased in some breeds such as Holstein. High-throughput SNP analysis and SNP microarrays give the opportunity to genotype many animals for hundreds SNPs per chromosome. In this study, due to these techniques a dense SNP marker map was used to fine map a QTL underlying nonreturn rate measured 90 days after artificial insemination previously detected with a low-density microsatellite marker map. A granddaughter design with 17 Holstein half-sib families (926 offspring) was genotyped for a set of 437 SNPs mapping to BTA3. Linkage analysis was performed by both regression and variance components analysis. An additional analysis combining both linkage analysis and linkage-disequilibrium information was applied. This method first estimated identity-by-descent probabilities among base haplotypes. These probabilities were then used to group the base haplotypes in different clusters. A QTL explaining 14% of the genetic variance was found with high significance (P < 0.001) at position 19 cM with the linkage analysis and four sires were estimated to be heterozygous (P < 0.05). Addition of linkage-disequilibrium information refined the QTL position to a set of narrow peaks. The use of the haplotypes of heterozygous sires offered the possibility to give confidence in some peaks while others could be discarded. Two peaks with high likelihood-ratio test values in the region of which heterozygous sires shared a common haplotype appeared particularly interesting. Despite the fact that the analysis did not fine map the QTL in a unique narrow region, the method proved to be able to handle efficiently and automatically a large amount of information and to refine the QTL position to a small set of narrow intervals. In addition, the QTL identified was confirmed to have a large effect (explaining 13.8% of the genetic variance) on dairy cow fertility as estimated by nonreturn rate at 90 days.