A quasi-exclusive European ancestry in the Senepol tropical cattle breed highlights the importance of the slick locus in tropical adaptation

Selective signatures in composite MONTANA TROPICAL beef cattle reveal potential genomic regions for tropical adaptation

Genomic regions related to tropical adaptability are of paramount importance for animal breeding nowadays, especially in the context of global climate change. Moreover, understanding the genomic architecture of these regions may be very relevant for aiding breeding programs in choosing the best selection scheme for tropical adaptation and/or implementing a crossbreeding scheme. The composite MONTANA TROPICAL® population was developed by crossing cattle of four different biological types to improve production in harsh environments. Pedigree and genotype data (51962 SNPs) from 3215 MONTANA TROPICAL® cattle were used to i) characterize the population structure; ii) identify signatures of selection with complementary approaches, i.e. Integrated Haplotype Score (iHS) and Runs of Homozygosity (ROH); and iii) understand genes and traits related to each selected region. The population structure based on principal components had a weak relationship with the genetic contribution of the different biological types. Clustering analyses (ADMIXTURE) showed different clusters according to the number of generations within the composite population. Considering results of both selection signatures approaches, we identified only one consensus region on chromosome 20 (35399405-40329703 bp). Genes in this region are related to immune function, regulation of epithelial cell differentiation, and cell response to ionizing radiation. This region harbors the slick locus which is related to slick hair and epidermis anatomy, both of which are related to heat stress adaptation. Also, QTLs in this region were related to feed intake, milk yield, mastitis, reproduction, and slick hair coat. The signatures of selection detected here arose in a few generations after crossbreeding between contrasting breeds. Therefore, it shows how important this genomic region may be for these animals to thrive in tropical conditions. Further investigations on sequencing this region can identify candidate genes for animal breeding and/or gene editing to tackle the challenges of climate change.

Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat-Wool Sheep Breed

Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat-wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat-wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat-wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1-2 Mb were much higher than ROH2-4 Mb and ROH8-16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat-wool sheep breed.

Physiological responses of Holstein calves and heifers carrying the SLICK1 allele to heat stress in California and Florida dairy farms

Inheritance of the SLICK1 allele of the prolactin receptor gene improves thermotolerance of lactating Holstein cows under humid heat stress conditions. The aim of this study was to investigate whether pre- and postweaning Holstein heifers carrying the SLICK1 allele would show physiological responses indicative of higher tolerance to heat stress in high- and low-humidity climates. A total of 101 heifer calves of two age groups heterozygous for the SLICK1 allele and 103 wild-type half-siblings were evaluated during July 2020 in 3 dairy farms in central California and 2 in south Florida. Dry bulb temperature and relative humidity data were recorded during evaluation and used to calculate the temperature-humidity index (THI). Physiological measurements were obtained between 1600 and 1900 h in California, and 1200 and 1400 h in Florida and included rectal temperature, respiration rate, skin temperature, and sweating rate. Data were analyzed via Generalized Linear Mixed Models including the main effects of genotype, state, group, sire, farm within state, and interactions, with THI included as a covariate. The correlations between THI and dependent variables were analyzed via linear regression. The average 24-h THI was higher in Florida compared with California (90 vs. 72, respectively); the main driver of the higher THI in Florida was the high relative humidity (average 85.6% in Florida vs. 36.7% in California). In Florida, the rectal temperature of slick calves was 0.4°C lower than non-slick calves (39.5 ± 0.1 vs 39.9 ± 0.1°C); no differences were detected between slick and non-slick calves in California. Regardless of genotype, heifer calves in Florida had higher respiration rate, higher rectal and skin temperatures, and lower sweating rate than in California. This study is the first to evaluate physiological responses of calves carrying the SLICK1 allele under heat stress conditions in different climates. Our findings demonstrate that the presence of this allele is associated with lower rectal temperatures in pre- and post-weaning Holstein females. According to the physiological parameters evaluated, calves raised in Florida appeared to be under more severe heat stress; in those conditions, the SLICK1 allele was advantageous to confer thermotolerance as evidenced by lower rectal temperature in slick animals.

Coat score. A possible explanation for the zebuine selective sweep located on bovine chromosome 5: 47,670,001-48,100,000 bp

Over 65% of the world's cattle population resides in warm areas where heat stress conditions limit the breed of European taurine cattle. Composite breeds were developed to retain the main traits of both parental breeds. The skin plays a central role in animal response to heat stress. Research on the genetic architecture of skin traits has identified genes and regions related to warm resistance skin features. The aim of this study was to determine whether the indicine proportion accounted for coat type or whether there were genes of large effect segregating in Brangus. Bulls (n = 108) were genotyped using microarrays and their coat score and hair length were evaluated. Indicine-taurine genome-wide composition was estimated and GWAS was performed. Although significant correlations between indicine proportion and traits were not observed, four windows of SNPs on BTA4 and BTA5 explained more than 2% of the trait variance. The GWAS for coat score in summer showed the main peak on BTA5:46,941,446-48,030,219 bp, accounting for 4.65% of the variance. Our results suggest that the variation in coat score and undercoat hair length in Argentinian Brangus bulls is associated with the presence of some particular gene variants, rather than with the whole indicine genetic content.

Physicochemical composition, fatty acid profile and sensory attributes of meat (longissimus lumborum muscle) from Nellore and Nellore-cross bulls

This study aimed to compare the physicochemical characteristics, fatty acid composition, and sensory attributes of the meat from three genetic groups: Nellore (Nell), ½ Nellore × ½ Angus (NeAn), and ¼ Nellore × ¼ Angus × ½ Senepol (NASe). Longissimus lumborum muscle from 30 slaughtered bulls with a body weight of 549 ± 32.5 kg was used. The water holding capacity was greater for the Nell and NeAn groups than for the NASe group. Meat samples from the NASe group exhibited a higher L* index than those from the NeAn group and lower a* and b* color indexes than those from the Nell group. The meat fatty acid profiles showed that the Nell group had higher concentrations of 12:0, 14:0, 18:1 t11, 14:1 c9, 16:1 c9, 18:1 c9, 18:1 c11, 18:2 c9, t11 (conjugated linoleic acid (CLA)), and 20:3 n-6 polyunsaturated fatty acid (PUFA) than the NeAn and NASe groups. The total saturated (ΣSFA), unsaturated (ΣUFA), and monounsaturated (ΣMUFA) fatty acid concentrations were higher and the ΣPUFA:ΣSFA ratio was lower in the Nell group than in the NeAn group. The Δ⁹-desaturase C16 activity was significantly higher in the Nell and NASe groups than in the NeAn group. The atherogenicity index (AI) tended to be lower in the crossbreeds than in the Nell breed. The NASe group presented meat with better tenderness, juiciness, and overall acceptance than the Nell and NeAn groups and was therefore the best genetic group for beef production of the tested groups.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species

BACKGROUND

Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available.

RESULTS

The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep).

CONCLUSIONS

We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.

The Identification of Runs of Homozygosity Gives a Focus on the Genetic Diversity and Adaptation of the "Charolais de Cuba" Cattle

Simple Summary

The Charolais de Cuba cattle is a tropical adapted breed founded in Cuba around 120 years ago from Charolais French specimens. Nowadays, it is still a closed breed and remains as a small population. In this work, we analyzed the inbreeding and diversity patterns, as well as the population size, of this recent adapted breed via a run of homozygosity (ROH) analysis. We found that the genomic inbreeding levels are higher in the Charolais de Cuba breed compared to French and British Charolais populations. Nevertheless, we detected that the effective population size experienced a very similar decline during the last century in the three Charolais populations studied. Finally, a number of regions with exceptional patterns of long homozygosity were identified in this breed, and these could be related to processes of adaptation to tropical conditions.

Abstract

Inbreeding and effective population size (Ne) are fundamental indicators for the management and conservation of genetic diversity in populations. Genomic inbreeding gives accurate estimates of inbreeding, and the Ne determines the rate of the loss of genetic variation. The objective of this work was to study the distribution of runs of homozygosity (ROHs) in order to estimate genomic inbreeding (F_ROH) and an effective population size using 38,789 Single Nucleotide Polymorphisms (SNPs) from the Illumina Bovine 50K BeadChip in 86 samples from populations of Charolais de Cuba (n = 40) cattle and to compare this information with French (n = 20) and British Charolais (n = 26) populations. In the Cuban, French, and British Charolais populations, the average estimated genomic inbreeding values using the F_ROH statistics were 5.7%, 3.4%, and 4%, respectively. The dispersion measured by variation coefficient was high at 43.9%, 37.0%, and 54.2%, respectively. The effective population size experienced a very similar decline during the last century in Charolais de Cuba (from 139 to 23 individuals), in French Charolais (from 142 to 12), and in British Charolais (from 145 to 14) for the ~20 last generations. However, the high variability found in the ROH indicators and F_ROH reveals an opportunity for maintaining the genetic diversity of this breed with an adequate mating strategy, which can be favored with the use of molecular markers. Moreover, the detected ROH were compared to previous results obtained on the detection of signatures of selection in the same breed. Some of the observed signatures were confirmed by the ROHs, emphasizing the process of adaptation to tropical climate experienced by the Charolais de Cuba population.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Prospects for gene introgression or gene editing as a strategy for reduction of the impact of heat stress on production and reproduction in cattle

In cattle, genetic variation exists in regulation of body temperature and stabilization of cellular function during heat stress. There are opportunities to reduce the impact of heat stress on cattle production by identifying the causative mutations responsible for genetic variation in thermotolerance and transferring specific alleles that confer thermotolerance to breeds not adapted to hot climates. An example of a mutation conferring superior ability to regulate body temperature is the group of frame-sift mutations in the prolactin receptor gene (PRLR) that lead to a truncated receptor and development of cattle with a short, sleek hair coat. Slick mutations in PRLR have been found in several extant breeds derived from criollo cattle. The slick mutation in Senepol cattle has been introgressed into dairy cattle in Puerto Rico, Florida and New Zealand. An example of a mutation that confers cellular protection against elevated body temperature is a deletion mutation in the promoter region of a heat shock protein 70 gene called HSPA1L. Inheritance of the mutation results in amplification of the transcriptional response of HSPA1L to heat shock and increased cell survival. The case of PRLR provides a promising example of the efficacy of the genetic approach outlined in this paper. Identification of other mutations conferring thermotolerance at the whole-animal or cellular level will lead to additional opportunities for using genetic solutions to reduce the impact of heat stress.

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short-term evolutionary response to climate pressure induced by their post-domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome-wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.

About the inconvenience of handling mixed-breed herds; aspects of social behavior as a potential source of stress and economic losses

The objective of this study was to assess aspects of the social behavior of a mixed-breed herd of beef cows as a potential source for stress and economic losses. Angus (AN; N = 10), Brahman (BR; N = 10), and Senepol (SE; N = 10) cows were assigned to two groups (N = 15 each containing equal breed numbers) on separate pastures. Agonistic interactions (win/loss) during feeding were recorded daily for 45 days. Dominance values were estimated as the proportion of individuals dominated to total herdmates. From this, individuals were placed into social categories based upon linear ranking as follows: dominants (D), intermediate (I), and subordinates (S). Breed influenced (P < 0.01) social category, with SE cows being dominants (P < 0.05) over AN and BR cows. Interactions between AN and BR cows were less (P < 0.0005) than interactions between AN and SE (53 vs 140, respectively). Within breeds, BR (152) and SE (182) cows had more (P < 0.0005) agonistic interactions than AN (107) cows. Although apparently influenced by breed, agonistic interactions occurred more frequently (P < 0.005) between social categories than within social categories (814 vs 310, respectively). Dominant cows were involved in more agonistic interactions with cows from different social categories than were intermediate and subordinate cows (P < 0.0005). However, intermediate (100) and subordinate (157) cows generated more (P < 0.0005) agonistic interactions within their own social category than dominant cows (53). It was concluded that, in mixed-breed herds, breed influences both social organization and agonistic interactions which could be considered as potential sources of stress and economic losses.

Polledness in Argentinean Creole cattle, five centuries surviving

Polledness has been shown to have autosomal Mendelian inheritance, with the polled locus being dominant to the horned locus. This trait was mapped to the BTA1 centromeric end in several breeds. One of the distinctive attributes of Creole cattle, such as the Argentinean Creole, is the presence of long, lyre-shaped horns. However, polled native animals were reported before the introduction of modern selected European breeds. Here, we studied the origin of the polled mutation, either independent or introgressed, in a Creole line from the Creole cattle founder group at the IIACS-INTA Leales Experimental Station (northwest Argentina). The study sample (65 animals: 26 horned and 39 polled) was genotyped using high-density SNP microarrays and three previously reported genetic markers (P_{202 ID} , P_{80kb ID} and P_G ). A genome-wide association study, selection signatures, linkage disequilibrium analysis and copy number variations were used to detect the responsible region and the segregating haplotypes/alleles. The interval mapped in the Leales herd (1.23-2.13 Mb) overlapped with the region previously reported in several European cattle breeds, suggesting that the same locus could be segregating in this population. The previously reported variants P_F and P_G were not detected, thus dismissing the Holstein-Friesian and Nellore origins of the polled phenotype in this native breed. Conversely, the presence of the Celtic variant P_C suggests an almost complete co-segregation. The cluster analysis rejected the hypothesis of recent introgression, which is compatible with the historical record of polled Creole cattle in northwest Argentina.

Origin of prehistoric cattle excavated from four archaeological sites in central and northeastern Thailand

Cattle have been domesticated in Southeast Asia, including Thailand, for thousands of years, but the history of cattle domestication in the region remains unclear. To date, only genetic studies of modern Thai cattle DNA have been reported. To gain some insight into cattle domestication in the country, a total of 56 cattle remains excavated from four archaeological sites (dated to between 3550 and 1700 years before present (YBP)) in northeastern and central Thailand were analysed in this study. Of 56, the 157-bp D-loop fragment was successfully generated from 26 samples, all of which belonged to Bos taurus in haplogroup T/T3. One haplotype contained 19 members from all four archaeological sites and clustered with the ancient B. taurus from Iran, Turkey and China. Other haplotypes have not shared haplotype with B. taurus from other countries but they showed close relationship to those from China. This represents the first genetic evidence that B. taurus was domesticated in Thailand between 3550 and 1700 YBP. In addition, the close relationship among ancient Thai, Iranian and Chinese taurines suggests that cattle from the Near East were introduced into North China, and were subsequently brought into Thailand thousands of years ago.

Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics

The introduction of Iberian cattle in the Americas after Columbus' arrival imposed high selection pressures on a limited number of animals over a brief period of time. Knowledge of the genomic regions selected during this process may help in enhancing climatic resilience and sustainable animal production. We first determined taurine and indicine contributions to the genomic structure of modern Creole cattle. Second, we inferred their demographic history using approximate Bayesian computation (ABC), linkage disequilibrium (LD) and N e Slope (NeS) analysis. Third, we performed whole genome scans for selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) and population differentiation (F ST) to disentangle the genetic mechanisms involved in adaptation and phenotypic change by a rapid and major environmental transition. To tackle these questions, we combined SNP array data (~54,000 SNPs) in Creole breeds with their modern putative Iberian ancestors. Reconstruction of the population history of Creoles from the end of the 15th century indicated a major demographic expansion until the introduction of zebu and commercial breeds into the Americas ~180 years ago, coinciding with a drastic N e contraction. NeS analysis provided insights into short-term complexity in population change and depicted a decrease/expansion episode at the end of the ABC-inferred expansion, as well as several additional fluctuations in N e with the attainment of the current small N e only towards the end of the 20th century. Selection signatures for tropical adaptation pinpointed the thermoregulatory slick hair coat region, identifying a new candidate gene (GDNF), as well as novel candidate regions involved in immune function, behavioural processes, iron metabolism and adaptation to new feeding conditions. The outcomes from this study will help in future-proofing farm animal genetic resources (FAnGR) by providing molecular tools that allow selection for improved cattle performance, resilience and welfare under climate change.

Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle

Background

This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection.

Results

Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (F_ST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4.

Conclusions

The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0381-2) contains supplementary material, which is available to authorized users.

Convergent Evolution of Slick Coat in Cattle through Truncation Mutations in the Prolactin Receptor

Evolutionary adaptations are occasionally convergent solutions to the same problem. A mutation contributing to a heat tolerance adaptation in Senepol cattle, a New World breed of mostly European descent, results in the distinct phenotype known as slick, where an animal has shorter hair and lower follicle density across its coat than wild type animals. The causal variant, located in the 11^th exon of prolactin receptor, produces a frameshift that results in a truncated protein. However, this mutation does not explain all cases of slick coats found in criollo breeds. Here, we obtained genome sequences from slick cattle of a geographically distinct criollo breed, namely Limonero, whose ancestors were originally brought to the Americas by the Spanish. These data were used to identify new causal alleles in the 11^th exon of the prolactin receptor, two of which also encode shortened proteins that remove a highly conserved tyrosine residue. These new mutations explained almost 90% of investigated cases of animals that had slick coats, but which also did not carry the Senepol slick allele. These results demonstrate convergent evolution at the molecular level in a trait important to the adaptation of an animal to its environment.

Climate Adaptation of Tropical Cattle

There is sustained growth in the number of tropical cattle, which represent more than half of all cattle worldwide. By and large, most research in tropical areas is still focused on breeds of cattle, their particular advantages or disadvantages in tropical areas, and the tropical forages or feeds that could be usefully fed to them. A consistent issue for adaptation to climate is the heat of tropical environments. Changing the external characteristics of the animal, such as color and coat characteristics, is one way to adapt, and there are several major genes for these traits. However, further improvement in heat tolerance and other adaptation traits will need to use the entire genome and all physical and physiological systems. Apart from the response to heat, climate forcing through methane emission identifies dry season weight loss as an important if somewhat neglected trait in climate adaptation of cattle. The use of genome-estimated breeding values in tropical areas is in its infancy and will be difficult to implement, but will be essential for rapid, coordinated genetic improvement. The difficulty of implementation cannot be exaggerated and may require major improvements in methodology.

BREEDING AND GENETICS SYMPOSIUM:Breeding heat tolerant dairy cattle: the case for introgression of the "slick" prolactin receptor variant into dairy breeds

Increasing environmental temperatures are a threat to the sustainability of livestock production and, because of the high metabolic demands of lactation, to dairy production in particular. Summer heat waves in temperate climates reduce feed intake, milk production, and cow comfort. In extreme heat events, there is an increase in cow mortality. In tropical climates, dairy cattle are mostly (zebu) type or zebu crossbred with temperate dairy breeds. Crossbreeding is undertaken to combine the heat tolerance and tick resistance of zebu with the productivity of temperate dairy breeds. In the absence of improved heat tolerance, milk production and fertility of temperate cattle is severely impaired. We have recently identified a key role for the prolactin pathway in regulating heat tolerance. A de novo mutation in prolactin that impairs prolactin activity was discovered in hairy and heat intolerant, New Zealand dairy cattle. The phenotypes produced were remarkably similar to those seen in fescue toxicosis, a syndrome seen in grazing cattle in the U.S. where ingestion of ergovaline, a fungal toxin from infected pasture, inhibits prolactin secretion. Recognition of the role of prolactin in hairy cattle led us to identify a deletion in exon 10 of the long-form of the prolactin receptor in Senepol cattle that causes truncation of the protein and determines the slick coat and heat tolerance traits found in this , beef breed. The short form of the prolactin receptor is predicted to be unaffected by the deletion. Knowledge of this dominant mutation has provided the impetus to begin a crossbreeding program to investigate performance and heat tolerance of temperate dairy cattle carrying the slick, prolactin receptor variant. The perceived opportunity is to introgress this variant into temperate dairy cattle to enable performance and welfare improvement in hot climates. Heat tolerance of cattle with slick coats appears to be mostly associated with coat type although sweating ability may also be enhanced. Further investigation is required of performance traits in cows homozygous for the slick variant because the published data are almost exclusively from heterozygous animals. Combination of the slick mutation with other favorable genes for heat tolerance, especially those for coat color, will be particularly enabled by gene editing technologies, offering opportunities for further improvement in bovine thermotolerance.

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F st ) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.

Signatures of Selection for Environmental Adaptation and Zebu × Taurine Hybrid Fitness in East African Shorthorn Zebu

The East African Shorthorn Zebu (EASZ) cattle are ancient hybrid between Asian zebu × African taurine cattle preferred by local farmers due to their adaptability to the African environment. The genetic controls of these adaptabilities are not clearly understood yet. Here, we genotyped 92 EASZ samples from Kenya (KEASZ) with more than 770,000 SNPs and sequenced the genome of a pool of 10 KEASZ. We observe an even admixed autosomal zebu × taurine genomic structure in the population. A total of 101 and 165 candidate regions of positive selection, based on genome-wide SNP analyses (meta-SS, Rsb, iHS, and ΔAF) and pooled heterozygosity (Hp) full genome sequence analysis, are identified, in which 35 regions are shared between them. A total of 142 functional variants, one novel, have been detected within these regions, in which 30 and 26 were classified as of zebu and African taurine origins, respectively. High density genome-wide SNP analysis of zebu × taurine admixed cattle populations from Uganda and Nigeria show that 25 of these regions are shared between KEASZ and Uganda cattle, and seven regions are shared across the KEASZ, Uganda, and Nigeria cattle. The identification of common candidate regions allows us to fine map 18 regions. These regions intersect with genes and QTL associated with reproduction and environmental stress (e.g., immunity and heat stress) suggesting that the genome of the zebu × taurine admixed cattle has been uniquely selected to maximize hybrid fitness both in terms of reproduction and survivability.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

WIDDE: a Web-Interfaced next generation database for genetic diversity exploration, with a first application in cattle

BACKGROUND

The advent and democratization of next generation sequencing and genotyping technologies lead to a huge amount of data for the characterization of population genetic diversity in model and non model-species. However, efficient storage, management, cross-analyzing and exploration of such dense genotyping datasets remain challenging. This is particularly true for the bovine species where many SNP datasets have been generated in various cattle populations with different genotyping tools.

DESCRIPTION

We developed WIDDE, a Web-Interfaced Next Generation Database that stands as a generic tool applicable to a wide range of species and marker types ( http://widde.toulouse.inra.fr). As a first illustration, we hereby describe its first version dedicated to cattle biodiversity, which includes a large and evolving cattle genotyping dataset for over 750,000 SNPs available on 129 (89 public) different cattle populations representative of the world-wide bovine genetic diversity and on 7 outgroup bovid species. This version proposes an optional marker and individual filtering step, an export of genotyping data in different popular formats, and an exploration of genetic diversity through a principal component analysis. Users can also explore their own genotyping data together with data from WIDDE, assign their samples to WIDDE populations based on distance assignment method and supervised clustering, and estimate their ancestry composition relative to the populations represented in the database.

CONCLUSION

The cattle version of WIDDE represents to our knowledge the first database dedicated to cattle biodiversity and SNP genotyping data that will be very useful for researchers interested in this field. As a generic tool applicable to a wide range of marker types, WIDDE is overall intended to the genetic diversity exploration of any species and will be extended to other species shortly. The structure makes it easy to include additional output formats and new tools dedicated to genetic diversity exploration.

The genomics of ecological vicariance in threespine stickleback fish

Populations occurring in similar habitats and displaying similar phenotypes are increasingly used to explore parallel evolution at the molecular level. This generally ignores the possibility that parallel evolution can be mimicked by the fragmentation of an ancestral population followed by genetic exchange with ecologically different populations. Here we demonstrate such an ecological vicariance scenario in multiple stream populations of threespine stickleback fish divergent from a single adjacent lake population. On the basis of demographic and population genomic analyses, we infer the initial spread of a stream-adapted ancestor followed by the emergence of a lake-adapted population, that selective sweeps have occurred mainly in the lake population, that adaptive lake–stream divergence is maintained in the face of gene flow from the lake into the streams, and that this divergence involves major inversion polymorphisms also important to marine-freshwater stickleback divergence. Overall, our study highlights the need for a robust understanding of the demographic and selective history in evolutionary investigations.

Threespine stickleback fish are adapted to lake and stream habitats in Central Europe. Here, the authors show colonization of a lake basin by a stream-adapted ancestor, followed by the emergence of a lake-adapted population in the face of gene flow across lake–stream boundaries.

Signatures of positive selection in East African Shorthorn Zebu: A genome-wide single nucleotide polymorphism analysis

The small East African Shorthorn Zebu (EASZ) is the main indigenous cattle across East Africa. A recent genome wide SNP analysis revealed an ancient stable African taurine x Asian zebu admixture. Here, we assess the presence of candidate signatures of positive selection in their genome, with the aim to provide qualitative insights about the corresponding selective pressures. Four hundred and twenty-five EASZ and four reference populations (Holstein-Friesian, Jersey, N’Dama and Nellore) were analysed using 46,171 SNPs covering all autosomes and the X chromosome. Following F_ST and two extended haplotype homozygosity-based (iHS and Rsb) analyses 24 candidate genome regions within 14 autosomes and the X chromosome were revealed, in which 18 and 4 were previously identified in tropical-adapted and commercial breeds, respectively. These regions overlap with 340 bovine QTL. They include 409 annotated genes, in which 37 were considered as candidates. These genes are involved in various biological pathways (e.g. immunity, reproduction, development and heat tolerance). Our results support that different selection pressures (e.g. environmental constraints, human selection, genome admixture constrains) have shaped the genome of EASZ. We argue that these candidate regions represent genome landmarks to be maintained in breeding programs aiming to improve sustainable livestock productivity in the tropics.

Revisiting demographic processes in cattle with genome-wide population genetic analysis

The domestication of the aurochs took place approximately 10,000 years ago giving rise to the two main types of domestic cattle known today, taurine (Bos taurus) domesticated somewhere on or near the Fertile Crescent, and indicine (Bos indicus) domesticated in the Indus Valley. However, although cattle have historically played a prominent role in human society the exact origin of many extant breeds is not well known. Here we used a combination of medium and high-density Illumina Bovine SNP arrays (i.e., ~54,000 and ~770,000 SNPs, respectively), genotyped for over 1300 animals representing 56 cattle breeds, to describe the relationships among major European cattle breeds and detect patterns of admixture among them. Our results suggest modern cross-breeding and ancient hybridisation events have both played an important role, including with animals of indicine origin. We use these data to identify signatures of selection reflecting both domestication (hypothesized to produce a common signature across breeds) and local adaptation (predicted to exhibit a signature of selection unique to a single breed or group of related breeds with a common history) to uncover additional demographic complexity of modern European cattle.

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.

Genomic data as the "hitchhiker's guide" to cattle adaptation: tracking the milestones of past selection in the bovine genome

The bovine species have witnessed and played a major role in the drastic socio-economical changes that shaped our culture over the last 10,000 years. During this journey, cattle "hitchhiked" on human development and colonized the world, facing strong selective pressures such as dramatic environmental changes and disease challenge. Consequently, hundreds of specialized cattle breeds emerged and spread around the globe, making up a rich spectrum of genomic resources. Their DNA still carry the scars left from adapting to this wide range of conditions, and we are now empowered with data and analytical tools to track the milestones of past selection in their genomes. In this review paper, we provide a summary of the reconstructed demographic events that shaped cattle diversity, offer a critical synthesis of popular methodologies applied to the search for signatures of selection (SS) in genomic data, and give examples of recent SS studies in cattle. Then, we outline the potential and challenges of the application of SS analysis in cattle, and discuss the future directions in this field.

Functionally reciprocal mutations of the prolactin signalling pathway define hairy and slick cattle

Lactation, hair development and homeothermy are characteristic evolutionary features that define mammals from other vertebrate species. Here we describe the discovery of two autosomal dominant mutations with antagonistic, pleiotropic effects on all three of these biological processes, mediated through the prolactin signalling pathway. Most conspicuously, mutations in prolactin (PRL) and its receptor (PRLR) have an impact on thermoregulation and hair morphology phenotypes, giving prominence to this pathway outside of its classical roles in lactation.

The hormone prolactin is a known modulator of mammalian lactation and hair growth. Here, the authors describe two dominant mutations in bovine prolactin and its receptor, demonstrating antagonistic effects on these traits and highlighting a role for this pathway in sweat gland function and thermoregulation.

Genomic signatures reveal new evidences for selection of important traits in domestic cattle

We investigated diverse genomic selections using high-density single nucleotide polymorphism data of five distinct cattle breeds. Based on allele frequency differences, we detected hundreds of candidate regions under positive selection across Holstein, Angus, Charolais, Brahman, and N'Dama. In addition to well-known genes such as KIT, MC1R, ASIP, GHR, LCORL, NCAPG, WIF1, and ABCA12, we found evidence for a variety of novel and less-known genes under selection in cattle, such as LAP3, SAR1B, LRIG3, FGF5, and NUDCD3. Selective sweeps near LAP3 were then validated by next-generation sequencing. Genome-wide association analysis involving 26,362 Holsteins confirmed that LAP3 and SAR1B were related to milk production traits, suggesting that our candidate regions were likely functional. In addition, haplotype network analyses further revealed distinct selective pressures and evolution patterns across these five cattle breeds. Our results provided a glimpse into diverse genomic selection during cattle domestication, breed formation, and recent genetic improvement. These findings will facilitate genome-assisted breeding to improve animal production and health.

On the History of Cattle Genetic Resources

Cattle are our most important livestock species because of their production and role in human culture. Many breeds that differ in appearance, performance and environmental adaptation are kept on all inhabited continents, but the historic origin of the diverse phenotypes is not always clear. We give an account of the history of cattle by integrating archaeological record and pictorial or written sources, scarce until 300 years ago, with the recent contributions of DNA analysis. We describe the domestication of their wild ancestor, migrations to eventually all inhabited continents, the developments during prehistory, the antiquity and the Middle Ages, the relatively recent breed formation, the industrial cattle husbandry in the Old and New World and the current efforts to preserve the cattle genetic resources. Surveying the available information, we propose three main and overlapping phases during the development of the present genetic diversity: (i) domestication and subsequent wild introgression; (ii) natural adaptation to a diverse agricultural habitat; and (iii) breed development.

The SLICK hair locus derived from Senepol cattle confers thermotolerance to intensively managed lactating Holstein cows

The SLICK haplotype (http://omia.angis.org.au/OMIA001372/9913/) in cattle confers animals with a short and sleek hair coat. Originally identified in Senepol cattle, the gene has been introduced into Holsteins. The objectives of the current study were to determine (1) whether lactating Holsteins with the slick hair phenotype have superior ability for thermoregulation compared with wild-type cows or relatives not inheriting the SLICK haplotype, and (2) whether seasonal depression in milk yield would be reduced in SLICK cows. In experiment 1, diurnal variation in vaginal temperature in the summer was monitored for cows housed in a freestall barn with fans and sprinklers. Vaginal temperatures were lower in slick-haired cows than in relatives and wild-type cows. In experiment 2, acute responses to heat stress were monitored after cows were moved to a dry lot in which the only heat abatement was shade cloth. The increases in rectal temperature and respiration rate caused by heat stress during the day were lower for slick cows than for relatives or wild-type cows. Moreover, sweating rate was higher for slick cows than for cows of the other 2 types. In experiment 3, effects of season of calving (summer vs. winter) on milk yield and composition were determined. Compared with milk yield of cows calving in winter, milk yield during the first 90 d in milk was lower for cows calving in the summer. However, this reduction was less pronounced for slick cows than for wild-type cows. In conclusion, Holsteins with slick hair have superior thermoregulatory ability compared with non-slick animals and experience a less drastic depression in milk yield during the summer.

Adaptive admixture in the West African bovine hybrid zone: insight from the Borgou population

Understanding the adaptive response to environmental fluctuations represents a central issue in evolutionary biology. Population admixture between divergent ancestries has often been considered as an efficient short-term adaptation strategy. Cattle populations from the West African Bos taurus × Bos indicus hybrid zone represent a valuable resource to characterize the effect of such adaptive admixture at the genome level. We here provide a detailed assessment of the global and local genome ancestries of the Borgou breed, one of the most representative cattle of this hybrid zone. We analysed a large data set consisting of 38,100 SNPs genotyped on 203 Borgou and 591 individuals representative of all the different cattle ancestries. At the global genomic level, we show that Borgou is a stabilized admixed breed whose origin (c. 130 years ago) traces back to the great African rinderpest pandemic, several centuries after the last admixture events, the West African zebus originate from (c. 500 years ago). To identify footprints of adaptive admixture, we combined the identification of signatures of selection and the functional annotation of the underlying genes using systems biology tools. The detection of the SILV coat coloration gene likely under artificial selection may be viewed as a validation of our approach. Overall, our results suggest that the long-term presence of pathogens and the intermediate environmental conditions are the main acting selective pressures. Our analytical framework can be extended to other model or nonmodel species to understand the process that shapes the patterns of genetic variability in hybrid zones.

Genome-wide association study and ancestral origins of the slick-hair coat in tropically adapted cattle

The slick hair coat (SLICK) is a dominantly inherited trait typically associated with tropically adapted cattle that are from Criollo descent through Spanish colonization of cattle into the New World. The trait is of interest relative to climate change, due to its association with improved thermo-tolerance and subsequent increased productivity. Previous studies localized the SLICK locus to a 4 cM region on chromosome (BTA) 20 and identified signatures of selection in this region derived from Senepol cattle. The current study compares three slick-haired Criollo-derived breeds including Senepol, Carora, and Romosinuano and three additional slick-haired cross-bred lineages to non-slick ancestral breeds. Genome-wide association (GWA), haplotype analysis, signatures of selection, runs of homozygosity (ROH), and identity by state (IBS) calculations were used to identify a 0.8 Mb (37.7-38.5 Mb) consensus region for the SLICK locus on BTA20 in which contains SKP2 and SPEF2 as possible candidate genes. Three specific haplotype patterns are identified in slick individuals, all with zero frequency in non-slick individuals. Admixture analysis identified common genetic patterns between the three slick breeds at the SLICK locus. Principal component analysis (PCA) and admixture results show Senepol and Romosinuano sharing a higher degree of genetic similarity to one another with a much lesser degree of similarity to Carora. Variation in GWA, haplotype analysis, and IBS calculations with accompanying population structure information supports potentially two mutations, one common to Senepol and Romosinuano and another in Carora, effecting genes contained within our refined location for the SLICK locus.

Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep

Background

Discerning the traits evolving under neutral conditions from those traits evolving rapidly because of various selection pressures is a great challenge. We propose a new method, composite selection signals (CSS), which unifies the multiple pieces of selection evidence from the rank distribution of its diverse constituent tests. The extreme CSS scores capture highly differentiated loci and underlying common variants hauling excess haplotype homozygosity in the samples of a target population.

Results

The data on high-density genotypes were analyzed for evidence of an association with either polledness or double muscling in various cohorts of cattle and sheep. In cattle, extreme CSS scores were found in the candidate regions on autosome BTA-1 and BTA-2, flanking the POLL locus and MSTN gene, for polledness and double muscling, respectively. In sheep, the regions with extreme scores were localized on autosome OAR-2 harbouring the MSTN gene for double muscling and on OAR-10 harbouring the RXFP2 gene for polledness. In comparison to the constituent tests, there was a partial agreement between the signals at the four candidate loci; however, they consistently identified additional genomic regions harbouring no known genes. Persuasively, our list of all the additional significant CSS regions contains genes that have been successfully implicated to secondary phenotypic diversity among several subpopulations in our data. For example, the method identified a strong selection signature for stature in cattle capturing selective sweeps harbouring UQCC-GDF5 and PLAG1-CHCHD7 gene regions on BTA-13 and BTA-14, respectively. Both gene pairs have been previously associated with height in humans, while PLAG1-CHCHD7 has also been reported for stature in cattle. In the additional analysis, CSS identified significant regions harbouring multiple genes for various traits under selection in European cattle including polledness, adaptation, metabolism, growth rate, stature, immunity, reproduction traits and some other candidate genes for dairy and beef production.

Conclusions

CSS successfully localized the candidate regions in validation datasets as well as identified previously known and novel regions for various traits experiencing selection pressure. Together, the results demonstrate the utility of CSS by its improved power, reduced false positives and high-resolution of selection signals as compared to individual constituent tests.

Information compression exploits patterns of genome composition to discriminate populations and highlight regions of evolutionary interest

Background

Genomic information allows population relatedness to be inferred and selected genes to be identified. Single nucleotide polymorphism microarray (SNP-chip) data, a proxy for genome composition, contains patterns in allele order and proportion. These patterns can be quantified by compression efficiency (CE). In principle, the composition of an entire genome can be represented by a CE number quantifying allele representation and order.

Results

We applied a compression algorithm (DEFLATE) to genome-wide high-density SNP data from 4,155 human, 1,800 cattle, 1,222 sheep, 81 dogs and 49 mice samples. All human ethnic groups can be clustered by CE and the clusters recover phylogeography based on traditional fixation index (F_ST) analyses. CE analysis of other mammals results in segregation by breed or species, and is sensitive to admixture and past effective population size. This clustering is a consequence of individual patterns such as runs of homozygosity. Intriguingly, a related approach can also be used to identify genomic loci that show population-specific CE segregation. A high resolution CE ‘sliding window’ scan across the human genome, organised at the population level, revealed genes known to be under evolutionary pressure. These include SLC24A5 (European and Gujarati Indian skin pigmentation), HERC2 (European eye color), LCT (European and Maasai milk digestion) and EDAR (Asian hair thickness). We also identified a set of previously unidentified loci with high population-specific CE scores including the chromatin remodeler SCMH1 in Africans and EDA2R in Asians. Closer inspection reveals that these prioritised genomic regions do not correspond to simple runs of homozygosity but rather compositionally complex regions that are shared by many individuals of a given population. Unlike F_ST, CE analyses do not require ab initio population comparisons and are amenable to the hemizygous X chromosome.

Conclusions

We conclude with a discussion of the implications of CE for a complex systems science view of genome evolution. CE allows one to clearly visualise the evolution of individual genomes and populations through a formal, mathematically-rigorous information space. Overall, CE makes a set of biological predictions, some of which are unique and await functional validation.

Genome-wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26-0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise FST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome-wide smoothed FST , 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty-four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.

Comparison of tick resistance of crossbred Senepol × Limousin to purebred Limousin cattle

The comparison of resistance to natural tick infestation by Rhipicephalus microplus (Canestrini, 1887) of crossbred Senepol × Limousin and purebred Limousin cattle was investigated. The Senepol breed, originated from St Croix Island in the Caribbean is considered as a Bos taurus breed adapted to tropical conditions. Despite its B. taurus genetic background, it is believed to have a good tick resistance, but this resistance has never been assessed previously. Tick counts under natural infestation were carried out to investigate the difference of susceptibility between crossbred Senepol × Limousin and purebred Limousin cattle. Mixed-effect models were used to assess the effect of the breed on the number of ticks. Results show that Senepol × Limousin are five times less infested by ticks than purebred Limousin. These results underline the opportunity to use Senepol cattle for crossing with susceptible B. taurus breeds in tick infested areas, to combine tick resistance with beef production abilities.