European Cattle Research Articles

The Practical Model and Evidence of Organic Evolution Opposite to Darwin’s Theory of Natural Selection and Sexual Selection (Human Evolution)

Hybridization (breeding) is practical evidence and a model of Darwin’s theory. But it would be true, if hybridization between two plants or animal species is possible and produced a fertile, reproductively isolated offspring. However, hybridization between two plants or animal species is not possible due to structural, behavioural differences, and seasonal isolations. If imposed, the fertilization fails, if the fertilization is successful, the embryo may abort, or the young may die. If the hybrid is survived up to maturity, it must become sterile. However, a very rare case the hybrids become fertile but those produce so-called varieties /races only; those species that produce fertile hybrids (e.g. Indian cattle Bos indicus and European cattle Bos taurus) must merge into a species.to satisfy the modern definition of species. Moreover, the artificial selection is also a skilled sexual selection, as the breeders choose the fittest, most vigour, and most fertile/productive, beautiful, colourful ornamented organism. But breeders also failed to develop a reproductively isolated species/variety/race by Johnson’s pure line selection, cloning, genetic engineering and mutation breeding. Even, a new species is not evolved by the natural hybridization. Consequently, there is no evidence of evolution of a new species either artificially or naturally. So, recent research claims that sexual selection theory is fundamentally flawed and simply wrong. Hence, evolutionary biologists rejected sexual selection. Thus, sexual selection is opposite to the evolution of humans from the lower animal like a chimpanzee. It is assumed that macroevolution occurs through hybridization; so, such an assumption is not valid.

Whole Genome Resequencing Reveals Genetic Diversity and Selection Signatures of Ethiopian Indigenous Cattle Adapted to Local Environments

Cattle are among the most important domesticated bovid species in the world, of which Ethiopia possesses large populations adapted to different agro-ecologies and production systems. Though several molecular population genetic studies have been done on Ethiopian indigenous cattle, genomic diversity and selection signatures pertinent to adaptation to the different local environments have yet to be comprehensively characterized. Here, the whole genome sequences of 151 samples from 14 Ethiopian indigenous cattle populations were analyzed to assess genomic diversity and differentiation as well as signatures of positive selection (using Hp, iHS, FST, and XP-CLR) in comparison to Sudanese zebu, Asian zebu, Ankole, and African and European taurine cattle. High genomic differentiation was observed between Ethiopian and non-Ethiopian cattle populations, while low genomic differentiation and inbreeding were present between and within Ethiopian cattle populations. Sixteen genome regions overlapping with 40 candidate genes were commonly identified by at least three genome scan methods. High frequencies of missense SNPs in ITPR2, CHADL, GNAS, STING1, and KIT genes with high haplotype differentiations were observed in Ethiopian cattle compared to non-Ethiopian cattle. The candidate genes were significantly associated with several biological functions and molecular pathways responsible for nutrient metabolism, skeletal development, immune response, reproduction, water balance, coat color pigmentation, and circulatory homeostasis. Our results provide new insights into the adaptation of the Ethiopian indigenous cattle to the country’s diverse environments.

Genome-wide analysis emancipates genomic diversity and signature of selection in Altay white-headed cattle of Xinjiang, China.

Altay white-headed cattle have not received enough attention for several reasons. Due to irrational breeding and selection practices, the number of pure Altay white-headed cattle has decreased significantly and the breed is now on the eve of extinction. The genomic characterization will be a crucial step towards understanding the genetic basis of productivity and adaptability to survival under native Chinese agropastoral systems; nevertheless, no attempt has been made in Altay white-headed cattle. In the current study, we compared the genomes of 20 Altay white-headed cattle to the genomes of 144 individuals in representative breeds. Population genetic diversity revealed that the nucleotide diversity of Altay white-headed cattle was less than that of indicine breeds and comparable to that of Chinese taurus cattle. Using population structure analysis, we also found that Altay white-headed cattle carried the ancestry of the European and East Asian cattle lineage. In addition, we used three different methods (F ST, θπ ratio and XP-EHH) to investigate the adaptability and white-headed phenotype of Altay white-headed cattle and compared it with Bohai black cattle. We found EPB41L5, SCG5 and KIT genes on the list of the top one percent genes, these genes might have an association with environmental adaptability and the white-headed phenotype for this breed. Our research reveals the distinctive genomic features of Altay white-headed cattle at the genome-wide level.

Transcriptome analysis reveals dysregulated gene expression networks in Sertoli cells of cattle-yak hybrids

Cattle-yak, the hybrid offspring of yak and taurine cattle, exhibits male sterility with normal female fertility. Spermatogenesis is arrested in adult cattle-yak, and apoptosis is elevated in spermatogenic cells. Currently, the mechanisms underlying these defects remain elusive. Sertoli cells are the only somatic cells that directly interact with spermatogenic cells in the seminiferous tubules and play essential roles in spermatogenesis. The present study was designed to investigate gene expression signatures and potential roles of Sertoli cells in hybrid sterility in cattle-yak. Immunohistochemical analysis showed that the 5 mC and 5hmC signals in Sertoli cells of cattle-yaks were significantly different from those of age-matched yaks (P < 0.05). Transcriptome profiling of isolated Sertoli cells identified 402 differentially expressed genes (DEGs) between cattle-yaks and yaks. Notably, niche factor glial cell derived neurotrophic factor (GDNF) was upregulated, and genes involved in retinoic acid (RA) biogenesis were changed in Sertoli cells of cattle-yak, suggesting possible impairments of spermatogonial fate decisions. Further studies showed that the numbers of proliferative gonocytes and undifferentiated spermatogonia in cattle-yak were significantly higher than those in yak (P < 0.01). Exogenous GDNF significantly promoted the proliferation of UCHL1-positive spermatogonia in yaks. Therefore, we concluded that altered GDNF expression and RA signaling impacted the fate decisions of undifferentiated spermatogonia in cattle-yak. Together, these findings highlight the role of Sertoli cells and their derived factors in hybrid sterility.

Analysis of genetic diversity and selection signals in Chaling cattle of southern China using whole-genome scan.

China has diversified resources of indigenous cattle, which are classified into Northern, Central, and Southern groups according to their geographical distribution. Chaling cattle belong to Southern group. This breed is famous for the production of good quality meat with elite meat grades. To analyze the genetic diversity of Chaling cattle, 20 samples were sequenced using whole-genome resequencing technology, along with 138 published whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle as control. It was found that Chaling cattle originated from Chinese indicine cattle. The genetic diversity of Chaling cattle is higher than that of Indian indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle, but lower than that of Chinese indicine cattle and Xiangxi cattle. Annotating the selection signals obtained by composite likelihood ratio, θπ, FST , π-ratio, and XP-EHH methods, several genes associated with immunity, heat tolerance, reproduction, growth, and meat quality showed strong selection signals. In general, this study provides a theoretical basis for analyzing the genetic mechanism of Chaling cattle with excellent adaptability, rough feeding tolerance, good immune performance, and good meat quality. This work lays a foundation for genetic breeding of Chaling cattle in future.

Genomic analysis reveals a KIT-related chromosomal translocation associated with the white coat phenotype inyak.

White coat pigmentation is a striking phenotype of many domesticated species and has various genetic controls. The Tianzhu White yak, an indigenous breed with a complete white coat, has fascinated Tibetans for centuries. However, the genetic basis of this trait remains unknown. Here, we conducted population genomics analysis and genome-wide association study based on the whole-genome sequencing data of 38 white and 59 non-white-coated yak. The results revealed the presence of KIT-linked Cs alleles characterized by the translocations between chromosomes 6 and 29 in all-white yak. Furthermore, structural variations showed additional duplications of the Cs alleles in white yak compared with colour-sidedness cattle. Interestingly, the Cs alleles associated with the white coat phenotype in yak were found to have introgressed from taurine cattle. Our findings unveil the shared genetic control of the white coat phenotype and its evolution in closely related bovine species.

Bazı Anadolu Yerli Sığır Irklarının MtDNA D-loop Dizi Analizi: Genetik Çeşitlilik ve Popülasyon Geçmişinin Değerlendirilmesi

Mitochondrial DNA is the molecular indicator of maternal inheritance and is one of the frequently preferred markers in population genomics, phylogenetic and phylogeographic studies due to its features such as evolving faster than genomic DNA, lack of recombination and showing differences according to the geographical distribution of species. In this research, mitochondrial DNA D-Loop region sequence analyzes were performed in order to reveal the genetic diversity and phylogenetic relationships of some domestic cattle breeds bred in Anatolia. First of all, the sampling study of the cattle breeds that were the subject of the research, was completed. Mitochondrial DNA D-Loop region was amplified by Polymerase Chain Reaction from the DNA isolated samples using the Standard Phenol/Chloroform Method. After DNA sequencing, data of domestic cattle breeds were analyzed together with reference mitochondrial DNA sequences from GenBANK, and their haplotype and nucleotide diversity, and genetic and phylogenetic relationships within and between populations were evaluated. With the data obtained, it was determined that the nucleotide and haplotype diversity, haplotype numbers, intra and inter population variation were quite high in Anatolian native cattle breeds. It was also concluded that the native cattle breeds bred in Anatolia are phylogenetically between Asian and European cattle breeds.

Giardia duodenalis Styles, 1902 Prevalence in Cattle (Bos taurus Linnaeus, 1758) in Europe: A Systematic Review.

Giardia duodenalis Styles, 1902 is an infectious agent which can cause enteritic disease in cattle (Bos taurus Linnaeus, 1758) worldwide. As a zoonotic protozoan, it is important to acknowledge Giardia prevalence and assemblages found in cattle and risk factors associated with the Giardia infection in herds. This systematic review aims to estimate the prevalence of G. duodenalis and its assemblages in cattle and to identify the risk factors associated with Giardia infection in cattle in Europe. A systematic review was performed according to Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines to collect information from published studies in Europe. A total of 1414 studies were identified and 17 relevant studies were included in this review. Mean Giardia prevalence in cattle in Europe was 35.1%, with the highest prevalence found in neonatal animals (39.6%), but mean heard prevalence was 67.0%. Mixed infections of Giardia assemblages A and E were found most frequently (55.6%), while assemblages A and B were found more frequently in animals under 24 months old. Risk factors, such as deep litter with run-out, prolonged calf contact with the dam, and seasonality, such as winter and spring, were found to be potential risk factors for Giardia presence in the herds.

Comparison between indicine and taurine cattle DNA methylation reveals epigenetic variation associated to differences in morphological adaptive traits

ABSTRACT Indicine and taurine subspecies present distinct morphological traits as a consequence of environmental adaptation and artificial selection. Although the two subspecies have been characterized and compared at genome-wide level and at specific loci, their epigenetic diversity has not yet been explored. In this work, Reduced Representation Bisulphite Sequencing (RRBS) profiling of the taurine Angus (A) and indicine Nellore (N) cattle breeds was applied to identify methylation differences between the two subspecies. Genotyping by sequencing (GBS) of the same animals was performed to detect single nucleotide polymorphisms (SNPs) at cytosines in CpG dinucleotides and remove them from the differential methylation analysis. A total of 660,845 methylated cytosines were identified within the CpG context (CpGs) across the 10 animals sequenced (5 N and 5 A). A total of 25,765 of these were differentially methylated (DMCs). Most DMCs clustered in CpG stretches nearby genes involved in cellular and anatomical structure morphogenesis. Also, sequences flanking DMC were enriched in SNPs compared to all other CpGs, either methylated or unmethylated in the two subspecies. Our data suggest a contribution of epigenetics to the regulation and divergence of anatomical morphogenesis in the two subspecies relevant for cattle evolution and sub-species differentiation and adaptation.

Genome-Wide DNA Methylation Differences between Bos indicus and Bos taurus.

Disease risk is a persistent problem in domestic cattle farming, while economic traits are the main concern. This study aimed to reveal the epigenetic basis for differences between zebu (Bos indicus) and taurine cattle (Bos taurus) in disease, disease resistance, and economic traits, and provide a theoretical basis for the genetic improvement of domestic cattle. In this study, whole genome bisulfite sequencing (WGBS) was used to analyze the whole-genome methylation of spleen and liver samples from Yunnan zebu and Holstein cattle. In the genome-wide methylation pattern analysis, it was found that the methylation pattern of all samples was dominated by the CG type, which accounted for >94.9%. The DNA methylation levels of different functional regions and transcriptional elements in the CG background varied widely. However, the methylation levels of different samples in the same functional regions or transcriptional elements did not differ significantly. In addition, we identified a large number of differentially methylation region (DMR) in both the spleen and liver groups, of which 4713 and 4663 were annotated to functional elements, and most of them were annotated to the intronic and exonic regions of genes. GO and KEGG functional analysis of the same differentially methylation region (DMG) in the spleen and liver groups revealed that significantly enriched pathways were involved in neurological, disease, and growth functions. As a result of the results of DMR localization, we screened six genes (DNM3, INPP4B, PLD, PCYT1B, KCNN2, and SLIT3) that were tissue-specific candidates for economic traits, disease, and disease resistance in Yunnan zebu. In this study, DNA methylation was used to construct links between genotypes and phenotypes in domestic cattle, providing useful information for further screening of epigenetic molecular markers in zebu and taurine cattle.

Enhanced progesterone support during stimulated cycles of transvaginal follicular aspiration improves bovine in vitro embryo production

The in vitro production (IVP) of cattle embryos requires that germinal-vesicle stage oocytes undergo a period of maturation in vitro prior to fertilization and culture to the blastocyst stage. Success of IVP in taurine cattle is enhanced following ovarian stimulation prior to oocyte retrieval (OPU), particularly if preceded by a short period of FSH withdrawal (‘coasting’). However, evidence regarding the importance of progesterone (P4) support during OPU-IVP is equivocal. The current study, therefore, determined the effects of increased peripheral P4 concentrations during FSH-stimulated (‘coasted’) cycles of OPU. Progesterone support was provided by either an active corpus luteum (CL) and/or one of two intravaginal P4 releasing devices (i.e., CIDR® [1.38 g P4] or PRID® Delta [1.55 g P4]). Expt. 1 established an initial estrus prior to OPU, allowing CL formation (single luteal phase) spanning the first two of five cycles of OPU; the remaining three cycles were supported by either a CIDR® or PRID® Delta. Expt. 2 commenced with two cycles of dominant follicle removal (including prostaglandin F2α) undertaken seven days apart prior to six cycles of OPU. The absence of a CL meant that these cycles were supported only by a CIDR® or PRID® Delta. As each experiment involved several sequential cycles of OPU, the cumulative effects of device use on vaginal discharges were also assessed. Each experiment involved 10 sexually mature Holstein heifers. In the absence of a CL, peak plasma P4 concentrations were greater (P = 0.002) for the PRID® Delta (4.3 ± 0.22) than for the CIDR® (2.9 ± 0.22). In Expt. 1 there was an interaction (P < 0.05) between CL presence at OPU and P4 device on Day 8 blastocyst yields, indicating an effect of P4 device only when the CL was absent. The percentage hatching/hatched blastocysts of matured oocytes for the CIDR® and PRID® Delta was 44.3 ± 5.04 and 41.0 ± 5.40 in the presence, and 17.1 ± 3.48 and 42.2 ± 3.76 in the absence, of a CL (P = 0.018). Combined analyses of data from Expt. 1 and 2, when no CL was present, confirmed that Day 8 blastocyst yields were greater (P = 0.022) for the PRID® Delta than the CIDR®. Vaginal discharge scores were higher (P < 0.001) for the PRID® Delta than the CIDR® in Expt. 1 but not in Expt 2; however scores were low, did not increase with repeated use, and thus were deemed of no clinical or welfare concern. In conclusion, enhanced P4 support during FSH-stimulated cycles of OPU-IVP can improve in vitro embryo development.

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1–13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

Understanding the uptake of diagnostics for sustainable gastrointestinal nematode control by European dairy cattle farmers: a multi-country cross-sectional study.

To mitigate emerging anthelmintic resistance (AR) in cattle, sustainable gastrointestinal nematode control strategies should be adopted. A multi-centre study was set up to understand the factors affecting European dairy cattle farmers' adoption of diagnostics and to gauge for differences between regions. The data were collected through a multi-lingual survey by participating countries of the European Co-operation in Science and Technology (COST) action COMbatting Anthelmintic Resistance in ruminants (COMBAR). Four countries provided sufficient data to be included in the data analysis: Norway, Italy, Germany and Austria. Three models were estimated and validated through structural equation modelling. Norway, along with Germany and Austria (pooled dataset) showed similar trends that align with previous studies. AR risk perception had no influence on the adoption intention of diagnostics, a positive influence was found for attitude towards diagnostics and subjective norms (i.e., perceived opinion of others), and a negative influence of attitudes towards anthelminthics. Additionally, routine (i.e., perception of the current treatment) had an indirect effect on adoption intention through attitudes. Italy's data deviated from these findings, presenting a positive effect of the perceived severity of AR, and perceived behavioural control (i.e., perceived ability to perform a specific behaviour) on adoption intention of diagnostics. Finally, Norway's data set allowed for inclusion of a measurement of current behaviour in the model, identifying a direct positive effect of the perceived actual behaviour of other farmers on their own behaviour.

Inbreeding and pedigree analysis of the European red dairy cattle

BackgroundRed dairy cattle breeds have an important role in the European dairy sector because of their functional characteristics and good health. Extensive pedigree information is available for these breeds and provides a unique opportunity to examine their population structure, such as effective population size, depth of the pedigree, and effective number of founders and ancestors, and inbreeding levels. Animals with the highest genetic contributions were identified. Pedigree data included 9,073,403 animals that were born between 1900 and 2019 from Denmark, Finland, Germany, Latvia, Lithuania, the Netherlands, Norway, Poland, and Sweden, and covered 32 breeds. The numerically largest breeds were Red Dairy Cattle and Meuse-Rhine-Yssel.ResultsThe deepest average complete generation equivalent (9.39) was found for Red Dairy Cattle in 2017. Mean pedigree completeness ranged from 0.6 for Finncattle to 7.51 for Red Dairy Cattle. An effective population size of 166 animals was estimated for the total pedigree and ranged from 35 (Rotes Höhenvieh) to 226 (Red Dairy Cattle). Average generation intervals were between 5 and 7 years. The mean inbreeding coefficient for animals born between 1960 and 2018 was 1.5%, with the highest inbreeding coefficients observed for Traditional Angler (4.2%) and Rotes Höhenvieh (4.1%). The most influential animal was a Dutch Meuse-Rhine-Yssel bull born in 1960. The mean inbreeding level for animals born between 2016 and 2018 was 2% and highest for the Meuse-Rhine-Yssel (4.64%) and Rotes Hohenvieh breeds (3.80%).ConclusionsWe provide the first detailed analysis of the genetic diversity and inbreeding levels of the European red dairy cattle breeds. Rotes Höhenvieh and Traditional Angler have high inbreeding levels and are either close to or below the minimal recommended effective population size, thus it is necessary to implement tools to monitor the selection process in order to control inbreeding in these breeds. Red Dairy Cattle, Vorderwälder, Swedish Polled and Hinterwälder hold more genetic diversity. Regarding the Meuse-Rhine-Yssel breed, given its decreased population size, increased inbreeding and low effective population size, we recommend implementation of a breeding program to prevent further loss in its genetic diversity.

Meta-analysis and meta-regression of core body temperature in taurine and zebuine cattle under different environmental conditions

Heat stress affects livestock productivity and is often determined by an increase in the core body temperature (CBT). To evaluate the latter, both vaginal and rectal temperatures can be used. The aims of this study were (1) to perform a meta-analysis in order to evaluate the impact of heat stress in comparison to thermoneutral conditions on CBT in Bos taurus and Bos taurus indicus cattle, and (2) to evaluate the effects of genetic group, sex, productive purpose, climate zone, and temperature-humidity index (THI class) heat stress on the vaginal and rectal temperatures of cattle. Eighty-nine studies were included in this meta-analysis. Systematic review and meta-analysis were performed following the PRISMA guidelines. The standardized mean difference (SMD) in CBT between the heat stress (HS) and thermoneutral (TN) environments was used as a measure of effect size, considering a 95% confidence interval (CI). Cochran's Q test and I² statistic indicated high variability among studies; thus, a random effects model was used. Subgroup and meta-regression analysis were performed to investigate the source of heterogeneity. Genetic group, sex, productive purpose, climate zone, and THI class were included as covariates. The SMD in CBT between the HS and TN environments was 7.45 °C (95%CI: 5.89;9.02 °C, P = 0.0001), indicating a physiological response to the heat stress condition. The SMD in CBT was 3.92 °C higher in Bos taurus than Bos taurus indicus animals (P < 0.05). Animals in climate zone D (warm temperate zone) showed a lower SMD in CBT than those in climate zone B (arid zone) (SMD = 4.04 ºC vs 12.38 ºC). The THI class affected the SMD in CBT (P < 0.05). The genetic group, climate zone, THI class and interaction between genetic group and climate zone explained 29.89% of the heterogeneity observed among studies. In conclusion, heat stress exposure affects the vaginal and rectal temperatures of cattle. The CBT measured by rectal and vaginal temperatures can be used to detect heat stress-induced changes in CBT. Bos taurus and Bos taurus indicus animals seem to respond differently to heat stress in terms of CBT.

Review on Bovine Respiratory Syncytial Virus Characteristic, Pathogenesis and Control Methods Applied for the Disease

Bovine Respiratory Syncytial Virus (BRSV) is an envelope, non-segmented, negative-stranded RNA virus belonging to the Pneumovirus genus within the subfamily Pneumovirinae, family Paramyxoviridae. The virus virion consists of a lipid envelope, derived from the host plasma membrane, containing three virally encoded trans-membrane surface glycoproteins, which are organized separately into spikes on the surface of the virion. The virus causes disease in cattle, sheep, goats and camelids where cattle are supposed to be a natural host and all other animal species listed may play as an epidemiological role in certain circumstances. Human Respiratory Syncytial Virus (HRSV) is also an important respiratory pathogen in infants and young children. The distribution of BRSV is worldwide and the virus has been isolated from cattle in Europe, America and Asia. The virus causes regular winter outbreaks of respiratory disease in cattle and a sero prevalence of 30–70% has been detected in this species. BRSV is one of the main causes of severe pneumonia, interstitial edema, and emphysema in cattle which spreads by infected animal aerosols, direct contact and transmission through objects (fomites). It can vary in its ability to cause disease, ranging from fatal to no clinical signs shown. BRSV can also make the animal more vulnerable to secondary infections, where the virus weakens the immune system so that bacteria that are usually harmless can cause disease. BRSV replicates primarily in the superficial layer of the respiratory ciliated epithelium and replication can also be detected in type II pneumocytes. Various types of proteins of this virus can contribute for its pathogenicity. Since the peak incidence of severe BRSV disease is between 2 and 6 months, an effective BRSV vaccine must be capable of stimulating an effective immune response within the first months of life. However, the presence of maternallyderived antibodies poses a major obstacle and deletion of non-essential genes represents an attractive option for production of a live, attenuated virus vaccine.

Favored single nucleotide variants identified using whole genome Re-sequencing of Austrian and Chinese cattle breeds

Cattle have been essential for the development of human civilization since their first domestication few thousand years ago. Since then, they have spread across vast geographic areas following human activities. Throughout generations, the cattle genome has been shaped with detectable signals induced by various evolutionary processes, such as natural and human selection processes and demographic events. Identifying such signals, called selection signatures, is one of the primary goals of population genetics. Previous studies used various selection signature methods and normalized the outputs score using specific windows, in kbp or based on the number of SNPs, to identify the candidate regions. The recent method of iSAFE claimed for high accuracy in pinpointing the candidate SNPs. In this study, we analyzed whole-genome resequencing (WGS) data of ten individuals from Austrian Fleckvieh (Bos taurus) and fifty individuals from 14 Chinese indigenous breeds (Bos taurus, Bos taurus indicus, and admixed). Individual WGS reads were aligned to the cattle reference genome of ARS. UCD1.2 and subsequently undergone single nucleotide variants (SNVs) calling pipeline using GATK. Using these SNVs, we examined the population structure using principal component and admixture analysis. Then we refined selection signature candidates using the iSAFE program and compared it with the classical iHS approach. Additionally, we run Fst population differentiation from these two cattle groups. We found gradual changes of taurine in north China to admixed and indicine to the south. Based on the population structure and the number of individuals, we grouped samples to Fleckvieh, three Chinese taurines (Kazakh, Mongolian, Yanbian), admixed individuals (CHBI_Med), indicine individuals (CHBI_Low), and a combination of admixed and indicine (CHBI) for performing iSAFE and iHS tests. There were more significant SNVs identified using iSAFE than the iHS for the candidate of positive selection and more detectable signals in taurine than in indicine individuals. However, combining admixed and indicine individuals decreased the iSAFE signals. From both within-population tests, significant SNVs are linked to the olfactory receptors, production, reproduction, and temperament traits in taurine cattle, while heat and parasites tolerance in the admixed individuals. Fst test suggests similar patterns of population differentiation between Fleckvieh and three Chinese taurine breeds against CHBI. Nevertheless, there are genes shared only among the Chinese taurine, such as PAX5, affecting coat color, which might drive the differences between these yellowish coated breeds, and those in the greater Far East region.

Seven Shades of Grey: A Follow-Up Study on the Molecular Basis of Coat Colour in Indicine Grey Cattle Using Genome-Wide SNP Data.

Shades of grey and brown are a dominant component in mammal coat colours, representing a fundamental trait involved in a great number of processes including cryptism, sexual selection and signalling. The genetic mechanisms of the grey colouration in mammals are very complex and controlled by hundreds of genes whose effects and interactions are still largely unclear. In this study, we adopted a robust multi-cohort Fst outlier approach based on pairwise contrasts between seven grey indicine cattle breeds and both taurine and indicine non-grey cattle breeds in order to find genomic regions potentially related to the grey colouration. On the basis of three main drawn settings, built in order to control both the effect of the sample size and the genetic structure, we have identified some signals common to those obtained in a previous work employing only taurine cattle. In particular, using the top 1% Fst approach, we detected a candidate region (22.6–23.8 megabases) on chromosome 14 in which genes related to pigmentation have been already documented. In addition, when we constructed a phylogenetic tree using the significant markers identified in this study and including also the genotyping data at these loci of both the grey taurine and the extinct wild auroch, we found a topological repartition consistent with breed colour pattern rather than with the known bovine evolutionary history. Thus, on the basis of this evidence, together with the geographical distribution of the current taurine grey cattle, an ancestral indicine origin for the grey phenotype would seem to be a conceivable interpretation. In this context, a higher thermo-tolerance and less UV-induced damage of the grey phenotype might have favoured the retention of advantageous genes into the taurine genome during the post-Neolithic human-mediated cattle expansions.

Long read genome assemblies complemented by single cell RNA-sequencing reveal genetic and cellular mechanisms underlying the adaptive evolution of yak

Wild yak (Bos mutus) and domestic yak (Bos grunniens) are adapted to high altitude environment and have ecological, economic, and cultural significances on the Qinghai-Tibetan Plateau (QTP). Currently, the genetic and cellular bases underlying adaptations of yak to extreme conditions remains elusive. In the present study, we assembled two chromosome-level genomes, one each for wild yak and domestic yak, and screened structural variants (SVs) through the long-read data of yak and taurine cattle. The results revealed that 6733 genes contained high-FST SVs. 127 genes carrying special type of SVs were differentially expressed in lungs of the taurine cattle and yak. We then constructed the first single-cell gene expression atlas of yak and taurine cattle lung tissues and identified a yak-specific endothelial cell subtype. By integrating SVs and single-cell transcriptome data, we revealed that the endothelial cells expressed the highest proportion of marker genes carrying high-FST SVs in taurine cattle lungs. Furthermore, we identified pathways which were related to the medial thickness and formation of elastic fibers in yak lungs. These findings provide new insights into the high-altitude adaptation of yak and have important implications for understanding the physiological and pathological responses of large mammals and humans to hypoxia.

Assessment of genotyping array performance for genome-wide association studies and imputation in African cattle

BackgroundIn cattle, genome-wide association studies (GWAS) have largely focused on European or Asian breeds, using genotyping arrays that were primarily designed for European cattle. Because there is growing interest in performing GWAS in African breeds, we have assessed the performance of 23 commercial bovine genotyping arrays for capturing the diversity across African breeds and performing imputation. We used 409 whole-genome sequences (WGS) spanning global cattle breeds, and a real cohort of 2481 individuals (including African breeds) that were genotyped with the Illumina high-density (HD) array and the GeneSeek bovine 50 k array.ResultsWe found that commercially available arrays were not effective in capturing variants that segregate among African indicine animals. Only 6% of these variants in high linkage disequilibrium (LD) (r2 > 0.8) were on the best performing arrays, which contrasts with the 17% and 25% in African and European taurine cattle, respectively. However, imputation from available HD arrays can successfully capture most variants (accuracies up to 0.93), mainly when using a global, not continent-specific, reference panel, which partially reflects the unusually high levels of admixture on the continent. When considering functional variants, the GGPF250 array performed best for tagging WGS variants and imputation. Finally, we show that imputation from low-density arrays can perform almost as well as HD arrays, if a two-stage imputation approach is adopted, i.e. first imputing to HD and then to WGS, which can potentially reduce the costs of GWAS.ConclusionsOur results show that the choice of an array should be based on a balance between the objective of the study and the breed/population considered, with the HD and BOS1 arrays being the best choice for both taurine and indicine breeds when performing GWAS, and the GGPF250 being preferable for fine-mapping studies. Moreover, our results suggest that there is no advantage to using the indicus-specific arrays for indicus breeds, regardless of the objective. Finally, we show that using a reference panel that better represents global bovine diversity improves imputation accuracy, particularly for non-European taurine populations.