Genetic position of Hungarian Grey among European cattle and identification of breed-specific markers

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Parentage testing and looking for single nucleotide markers associated with antler quality in deer (Cervus elaphus)

To provide a cost-efficient parentage testing kit for red deer (Cervus elaphus), a 63 SNP set has been developed from a high-density Illumina BovineHD BeadChip containing 777 962 SNPs after filtering of genotypes of 50 stags. The successful genotyping rate was 38.6 % on the chip. The ratio of polymorphic loci among effectively genotyped loci was 6.5 %. The selected 63 SNPs have been applied to 960 animals to perform parentage control. Thirty SNPs out of the 63 had worked on the OpenArray platform. Their combined value of the probability of identity and exclusion probability was 4.9×10-11 and 0.99803, respectively.

A search for loci linked with antler quality was also performed on the genotypes of the above-mentioned stags. Association studies revealed 14 SNPs associated with antler quality, where low-quality antlers with short and thin main beam antlers had values from 1 to 2, while high-quality antlers with long and strong main beams had values between 4 and 5. The chance for a stag to be correctly identified as having high-value antlers is expected to be over 88 %.

Genetic adaptation to urban living: molecular DNA analyses of wild boar populations in Budapest and surrounding area

Studies of wild boar, Sus scrofa Linnaeus 1758, in urban and suburban areas of Budapest, Hungary, have indicated that these populations do not have continuous contact. Based on the assumption that the city has a discrete population, we hypothesized that the urban wild boar would differ genetically from those in suburban areas. Analysis of single-nucleotide polymorphism (SNP) data using the GeneSeek Genomic Profiler (GGP) Porcine 50 K system (Neogen, Scotland, UK) differentiated three populations: Buda (B) from the Western bank of the Danube; Buda Surrounding (BS); and Valkó (V) from the Eastern bank of the Danube. The coefficient of genetic differentiation (FST) for the B and BS populations was low. The inbreeding coefficients of the populations BS and V were close to zero, while population B had a high positive value reflecting the influence of founders and the inbreeding of the continuous urban population. The genome regions that were most differentiated between the B and BS populations were analyzed based on the FST values of the SNP markers using a mixed linear multi-locus model and BayeScan software. The most differentiated marker, WU_10.2_18_56278226, was found on chromosome 18. The surrounding region contained several candidate genes that could play important roles in adaptations related to human-induced stress. Two of these, encoding the adenylate cyclase 1 (ADCY1) and inhibin beta A chain precursor (INHBA) genes, were sequenced. While IHBA gene did not display variation, the allele distribution of the ADCY1 gene in the B population was significantly different from that of the BS population supporting the parapatric differentiation of wild boar.

On the origin and diversification of Podolian cattle breeds: testing scenarios of European colonization using genome-wide SNP data

Background

During the Neolithic expansion, cattle accompanied humans and spread from their domestication centres to colonize the ancient world. In addition, European cattle occasionally intermingled with both indicine cattle and local aurochs resulting in an exclusive pattern of genetic diversity. Among the most ancient European cattle are breeds that belong to the so-called Podolian trunk, the history of which is still not well established. Here, we used genome-wide single nucleotide polymorphism (SNP) data on 806 individuals belonging to 36 breeds to reconstruct the origin and diversification of Podolian cattle and to provide a reliable scenario of the European colonization, through an approximate Bayesian computation random forest (ABC-RF) approach.

Results

Our results indicate that European Podolian cattle display higher values of genetic diversity indices than both African taurine and Asian indicine breeds. Clustering analyses show that Podolian breeds share close genomic relationships, which suggests a likely common genetic ancestry. Among the simulated and tested scenarios of the colonization of Europe from taurine cattle, the greatest support was obtained for the model assuming at least two waves of diffusion. Time estimates are in line with an early migration from the domestication centre of non-Podolian taurine breeds followed by a secondary migration of Podolian breeds. The best fitting model also suggests that the Italian Podolian breeds are the result of admixture between different genomic pools.

Conclusions

This comprehensive dataset that includes most of the autochthonous cattle breeds belonging to the so-called Podolian trunk allowed us not only to shed light onto the origin and diversification of this group of cattle, but also to gain new insights into the diffusion of European cattle. The most well-supported scenario of colonization points to two main waves of migrations: with one that occurred alongside with the Neolithic human expansion and gave rise to the non-Podolian taurine breeds, and a more recent one that favoured the diffusion of European Podolian. In this process, we highlight the importance of both the Mediterranean and Danube routes in promoting European cattle colonization. Moreover, we identified admixture as a driver of diversification in Italy, which could represent a melting pot for Podolian cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00639-w.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.

Genetic status of lowland-type Racka sheep colour variants

Lowland-type Racka is an indigenous sheep breed that beside Hungarian Grey cattle and Mangalitza pig is one of the national symbols of Hungary. However, the genetic description of Racka sheep has not yet been conducted based on whole-genome screening. By using the Geneseek Ovine SNP50 BeadChip, we have sampled the genome of 126 Black and 128 White Racka sheep. For comparative purposes, we used 134 Hungarian Merinos and further 3345 animals from 81 different breeds have been included from an available database. Performance of a multidimensional scaling plot showed that White and Black Rackas represent well-separated groups among other sheep breeds and clustered separately from each other. However, the number and total length of Runs of Homozygosity was similar to other sheep breeds, except Soay. The inbreeding coefficients (method-of-moments relatedness F coefficient) of Black and White Racka were 0.147 and 0.133, respectively. Based on multidimensional scaling and admixture analyses and on comparisons of genetic distances of the investigated 84 populations, we suggest considering the colour variants of Racka as genetically differentiated breeds. The most differentiated markers between Black and White Racka highlight several candidate genes including 5-Hydroxytryptamine Receptor 5A, Insulin Induced Gene 1, Cyclin Dependent Kinase 5 and Melanocortin 1 Receptor. The results of this study help the recognition of Racka as a unique genetic resource among sheep and pave the way of application of genome screens to guide the resolution of questions arising among breeders.

Insights into Genetic Diversity, Runs of Homozygosity and Heterozygosity-Rich Regions in Maremmana Semi-Feral Cattle Using Pedigree and Genomic Data

Semi-feral local livestock populations, like Maremmana cattle, are the object of renewed interest for the conservation of biological diversity and the preservation and exploitation of unique and potentially relevant genetic material. The aim of this study was to estimate genetic diversity parameters in semi-feral Maremmana cattle using both pedigree- and genomic-based approaches (FIS and FROH), and to detect regions of homozygosity (ROH) and heterozygosity (ROHet) in the genome. The average heterozygosity estimates were in the range reported for other cattle breeds (HE=0.261, HO=0.274). Pedigree-based average inbreeding (F) was estimated at 4.9%. The correlation was low between F and genomic-based approaches (r=0.03 with FIS, r=0.21 with FROH), while it was higher between FIS and FROH (r=0.78). The low correlation between F and FROH coefficients may be the result of the limited pedigree depth available for the animals involved in this study. The ROH islands identified in Maremmana cattle included candidate genes associated with climate adaptation, carcass traits or the regulation of body weight, fat and energy metabolism. The ROHet islands contained candidate genes associated with nematode resistance and reproduction traits in livestock. The results of this study confirm that genome-based measures like FROH may be useful estimators of individual autozygosity, and may provide insights on pedigree-based inbreeding estimates in cases when animals' pedigree data are unavailable, thus providing a more detailed picture of the genetic diversity.