Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information

Divergent adaptation to highland and tropical environments in Bolivian Creole cattle

Bolivian Creole cattle populations evolved under low levels of breeding management and, during more than 500 years of natural selection, became adapted to various environments such as the contrasting highland and subtropical environments. Recently, highland Creole cattle were crossbred with Holstein to improve dairy production. The aim of this research was to evaluate the divergent adaptation through selection footprints of Bolivian Creole cattle from Andean highland and tropical lowlands, and to evaluate the effect of Holstein introgression in highland Creole. For this purpose, 130 Creole cattle (75 highland, 55 lowland) and 88 Holstein were genotyped using a microarray. The database was used to determine population structure and admixture and detect selection sweeps using FST, Rsb, XP-EHH, and ROH. Ancestry inference suggested that selection peaks were not due to Holstein introgression. The NCBI database was used to retrieve genes from the common regions and then perform gene ontology analysis. The most prominent selection peaks were on BTA20 and BTA23 and included the PRLR (slick phenotype) and Class I and IIa BoLA genes. Other windows contained candidate genes for hypoxia (ANXA2, NDUFA4L2), angiogenesis and haematological parameters (ANXA2, CPLANE1, NRP1, NRP2), immune response (IL7R, IL6ST, IL31RA, C6, C7, STAT6, NKG2A, IRAK4, KLR, CLEC), oxidative stress (GSTA, HSD17B6) and morphological traits (PLAG1, CHCHD7, CAP2, ARL15). GO analysis revealed enrichment terms and pathways related to immune response, glutathione and retinol metabolism and reported QTLs for coat characteristics, immune response and tick resistance. The results suggest the complex mechanism in the adaptation of Bolivian Creole cattle to the contrasting highland and subtropical environments.

Comparative genome analysis of international transboundary cattle breeds

Indian-origin Gir and Ongole cattle are international transboundary breeds that are reared in Brazil, The United States, Mexico, Malaysia, Panama, and other nations to provide meat and dairy products. These breeds have shown substantial genetic diversity in recent years and well suited to the ecological niche in Brazil. 90 cattle samples of Indian Gir (n = 15), Ongole (n = 17), Brazilian Gyr (n = 27), and Nellore (n = 31) breeds were genotyped using Illumina BovineHD BeadChip. Samples were analyzed to identify selection signatures using two complementing approaches: Integrated Haplotype Score (iHS) and Fixation Index (FST). Gir versus Gyr and Ongole versus Nellore revealed Pairwise FST differences of 2.85 % and 2.35 %, respectively. Using integrated haplotype score (iHS) method, 4004, 3322, 3437, and 3485 genes were found in Gir, Gyr, Ongole, and Nellore, respectively, underlying top 1 % of selected regions. Under top 1 % of selected regions, FST based method identified1897 genes for the Ongole-Nellore pair and 1966 genes for the Gir-Gyr pair. Runs of homozygosity (ROH) analysis revealed that both recent as well as ancient inbreeding in these breeds were in range of 2.6-4.5 % indicating populations to be less inbred. Numerous candidate genes, including IER5, MILR1 (immunity related traits) in Gir; and FGF12, SV2C, JMY (average daily gain, body size, reproduction related traits) in Ongole, were found under the top-selected regions. Nellore breed had carcass/growth traits (PARP2, and KCNJ11) and genes linked to mammary gland development, udder size, and carcass (MYO16, MYO1B) were found in Gyr. Present findings reveals that Brazilian cattle population (Gyr and Nellore) is more selected for carcass and growth traits along with milk production traits, whereas in Indian cattle population (Gir and Ongole) selection signature related to immunity and adaptation were more prominent. Further, sufficient genetic diversity exist within these cattle breeds for their genetic improvement.

Genomic analysis reveals population structure and selection signatures in plateau dairy cattle

BACKGROUND

To solve the problem of an insufficient supply of dairy products in Tibet, work has been carried out to improve native dairy cattle and introduce purebred dairy cattle from low-altitude areas. The harsh environment of the plateau not only severely limits the production performance of high-yielding dairy cattle, such as Holstein and Jersey cattle, but also challenges their survival. The population structure and plateau adaptation mechanism of plateau dairy cattle are rarely reported. In this study, key genes and pathways affecting plateau purebred and crossbred dairy cattle were explored using genetic chip information.

RESULTS

The results showed that the genetic diversity of the Tibet dairy cattle population was higher than that of the native cattle and plains dairy cattle. Purebred Holstein and Jersey cattle in Tibet were genetically closer to dairy cattle in the plains, and crossbred dairy cattle were admixed with more Tibet cattle and Apaijiza cattle. Based on the fixation index (FST), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches, 60 and 40 genes were identified in plateau Holstein and Jersey cattle, respectively. A total of 78 and 70 genes were identified in crossbred cattle compared to Holstein and Tibet cattle respectively. These genes are related to cardiac health and development, neuronal development and function, angiogenesis and hematopoietic, pigmentation, growth and development, and immune response.

CONCLUSIONS

Our results provide a glimpse into diverse selection signatures in plateau dairy cattle, which can be used to enhance our understanding of the genomic basis of plateau adaptation in dairy cattle. These results support further research on breeding strategies such as marker-assisted selection and gene editing in plateau dairy cattle populations.

Advancing the Indian cattle pangenome: characterizing non-reference sequences in Bos indicus

BACKGROUND

India harbors the world's largest cattle population, encompassing over 50 distinct Bos indicus breeds. This rich genetic diversity underscores the inadequacy of a single reference genome to fully capture the genomic landscape of Indian cattle. To comprehensively characterize the genomic variation within Bos indicus and, specifically, dairy breeds, we aim to identify non-reference sequences and construct a comprehensive pangenome.

RESULTS

Five representative genomes of prominent dairy breeds, including Gir, Kankrej, Tharparkar, Sahiwal, and Red Sindhi, were sequenced using 10X Genomics 'linked-read' technology. Assemblies generated from these linked-reads ranged from 2.70 Gb to 2.77 Gb, comparable to the Bos indicus Brahman reference genome. A pangenome of Bos indicus cattle was constructed by comparing the newly assembled genomes with the reference using alignment and graph-based methods, revealing 8 Mb and 17.7 Mb of novel sequence respectively. A confident set of 6,844 Non-reference Unique Insertions (NUIs) spanning 7.57 Mb was identified through both methods, representing the pangenome of Indian Bos indicus breeds. Comparative analysis with previously published pangenomes unveiled 2.8 Mb (37%) commonality with the Chinese indicine pangenome and only 1% commonality with the Bos taurus pangenome. Among these, 2,312 NUIs encompassing ~ 2 Mb, were commonly found in 98 samples of the 5 breeds and designated as Bos indicus Common Insertions (BICIs) in the population. Furthermore, 926 BICIs were identified within 682 protein-coding genes, 54 long non-coding RNAs (lncRNA), and 18 pseudogenes. These protein-coding genes were enriched for functions such as chemical synaptic transmission, cell junction organization, cell-cell adhesion, and cell morphogenesis. The protein-coding genes were found in various prominent quantitative trait locus (QTL) regions, suggesting potential roles of BICIs in traits related to milk production, reproduction, exterior, health, meat, and carcass. Notably, 63.21% of the bases within the BICIs call set contained interspersed repeats, predominantly Long Interspersed Nuclear Elements (LINEs). Additionally, 70.28% of BICIs are shared with other domesticated and wild species, highlighting their evolutionary significance.

CONCLUSIONS

This is the first report unveiling a robust set of NUIs defining the pangenome of Bos indicus breeds of India. The analyses contribute valuable insights into the genomic landscape of desi cattle breeds.

Selection signatures associated with adaptation in South African Drakensberger, Nguni, and Tuli beef breeds

In the present study 1,709 cattle, including 1,118 Drakensberger (DRB), 377 Nguni (NGI), and 214 Tuli (TUL), were genotyped using the GeneSeek® Genomic Profiler™ 150 K bovine SNP panel. A genomic data set of 122,632 quality-filtered single nucleotide polymorphisms (SNPs) were used to identify selection signatures within breeds based on conserved runs of homozygosity (ROH) and heterozygosity (ROHet) estimated with the detectRUNS R package. The mean number of ROH per animal varied across breeds ranging from 36.09 ± 12.82 (NGI) to 51.82 ± 21.01 (DRB), and the mean ROH length per breed ranged between 2.31 Mb (NGI) and 3.90 Mb (DRB). The smallest length categories i.e., ROH < 4 Mb were most frequent, indicating historic inbreeding effects for all breeds. The ROH based inbreeding coefficients (FROH) ranged between 0.033 ± 0.024 (NGI) and 0.081 ± 0.046 (DRB). Genes mapped to candidate regions were associated with immunity (ADAMTS12, LY96, WDPCP) and adaptation (FKBP4, CBFA2T3, TUBB3) in cattle and genes previously only reported for immunity in mice and human (EXOC3L1, MYO1G). The present study contributes to the understanding of the genetic mechanisms of adaptation, providing information for potential molecular application in genetic evaluation and selection programs.

Combined Use of Univariate and Multivariate Approaches to Detect Selection Signatures Associated with Milk or Meat Production in Cattle

OBJECTIVES

The aim of this study was to investigate the genomic structure of the cattle breeds selected for meat and milk production and to identify selection signatures between them.

METHODS

A total of 391 animals genotyped at 41,258 SNPs and belonging to nine breeds were considered: Angus (N = 62), Charolais (46), Hereford (31), Limousin (44), and Piedmontese (24), clustered in the Meat group, and Brown Swiss (42), Holstein (63), Jersey (49), and Montbéliarde (30), clustered in the Milk group. The population stratification was analyzed by principal component analysis (PCA), whereas selection signatures were identified by univariate (Wright fixation index, FST) and multivariate (canonical discriminant analysis, CDA) approaches. Markers with FST values larger than three standard deviations from the chromosomal mean were considered interesting. Attention was focused on markers selected by both techniques.

RESULTS

A total of 10 SNPs located on seven different chromosomes (7, 10, 14, 16, 17, 18, and 24) were identified. Close to these SNPs (±250 kb), 165 QTL and 51 genes were found. The QTL were grouped in 45 different terms, of which three were significant (Bonferroni correction < 0.05): milk fat content, tenderness score, and length of productive life. Moreover, genes mainly associated with milk production, immunity and environmental adaptation, and reproduction were mapped close to the common SNPs.

CONCLUSIONS

The results of the present study suggest that the combined use of univariate and multivariate approaches can help to better identify selection signatures due to directional selection.

Genome-wide scan for selection signatures in Mexican Sardo Negro Zebu cattle

The Sardo Negro cattle (SN) is the only zebu cattle breed developed in Mexico. Since its development, the selection could have led to an increase in the homozygosity level in some regions of the genome and made differentiation with other cattle populations. We aimed to identify and characterize selection signatures in SN using medium-density SNP data using four approaches: 1) Runs of homozygosity (ROH) 2) Nucleotide Diversity 3) Tajima's D and 4) the Wright's fixation index (FST). A sample of 555 SN animals genotyped for 65k SNPs was used to obtain ROH segments considered regions under selection. The FST values were estimated by comparing the sample of genotyped SN animals with samples of genotyped animals from the Gir, Brahman, and Ongole breeds. Only one region mapped to 35.78-42.51 Mb on BTA6 was considered a selection signature by the ROH method. This selection signature overlapped with the lowest diversity, negative values of Tajima's D and a diversification region between SN and the other Zebu breeds by FST. We found several candidate genes (LCORL, NCAPG, and SLIT2) related to growth and other economically important productive traits in this common region. Using the FST method, different regions, such as regions on BTA8 (8:93.4-93.9 Mb), BTA11 (11:99.2-99.7), and BTA14 (14: 26.1-26.8) related to growth and milk traits also were defined as candidate selection signatures. The selective signals identified in this study reflected the direction of the selection pressure that primarily involves the increase of live weight traits in the Sardo Negro cattle breeding program.

Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India

The peculiarity of Indian cattle lies in milk quality, resistance to diseases and stressors as well as adaptability. The investigation addressed selection signatures in Gir and Tharparkar cattle, belonging to arid ecotypes of India. Double digest restriction-site associated DNA sequencing (ddRAD-seq) yielded nearly 26 million high-quality reads from unrelated seven Gir and seven Tharparkar cows. In all, 19,127 high-quality SNPs were processed for selection signature analysis. An approach involving within-population composite likelihood ratio (CLR) statistics and between-population FST statistics was used to capture selection signatures within and between the breeds, respectively. A total of 191 selection signatures were addressed using CLR and FST approaches. Selection signatures overlapping 86 and 73 genes were detected as Gir- and Tharparkar-specific, respectively. Notably, genes related to production (CACNA1D, GHRHR), reproduction (ESR1, RBMS3), immunity (NOSTRIN, IL12B) and adaptation (ADAM22, ASL) were annotated to selection signatures. Gene pathway analysis revealed genes in insulin/IGF pathway for milk production, gonadotropin releasing hormone pathway for reproduction, Wnt signalling pathway and chemokine and cytokine signalling pathway for adaptation. This is the first study where selection signatures are identified using ddRAD-seq in indicine cattle breeds. The study shall help in conservation and leveraging genetic improvements in Gir and Tharparkar cattle.

Assessment of genetic diversity, inbreeding and collection completeness of Jersey bulls in the US National Animal Germplasm Program

Genomic selection and extensive use of a few elite bulls through artificial insemination are leading to reduced genetic diversity in Jersey cattle. Conservation of genetic diversity through gene banks can protect a breed's genetic diversity and genetic gain, ensuring continued genetic advancement in the future. The availability of genomic information in the US National Animal Germplasm Program (NAGP) facilitates characterization of Jersey bulls in the germplasm collection. Therefore, in this study, we compared the genetic diversity and inbreeding between Jersey bulls in the NAGP and the national cooperator database (NCD). The NCD is maintained and curated by the Council on Dairy Cattle Breeding (CDCB). We found the genetic diversity to be marginally higher in NAGP (Ho = 0.34 ± 0.17) relative to the NCD population (Ho = 0.33 ± 0.16). The average pedigree and genomic inbreeding (FPED, FGRM, FROH > 2Mb) were similar between the groups, with estimates of 7.6% with FPED, 11.07% with FGRM and 20.13% with FROH > 2Mb. An increasing trend in inbreeding was detected, and a significantly higher level of inbreeding was estimated among the older bulls in the NAGP collection, suggesting an overrepresentation of the genetics from elite bulls. Results from principal component analyses (PCA) provided evidence that the NAGP collection is representative of the genetic variation found in the NCD population and a broad majority of the loci segregating (98.2%) in the NCD population were also segregating in the NAGP. Ward's clustering was used to assess collection completeness of Jerseys in the NAGP by comparison with top 1000 sires of bulls, top 1000 sires of cow, and bulls with high Lifetime Net Merit (NM$). All the clusters were represented in the NAGP suggesting that most of the genetic diversity in the US Jersey population is represented in the NAGP and confirmed the PCA results. The decade of birth was the major driver grouping bulls into clusters, suggesting the importance of selection over time. Selection signature analysis between the historic bulls in the NAGP with the newer bulls, born in the decade after implementation of genomic selection, identified selection for milk production, fat and protein yield, fertility, health, and reproductive traits. Cluster analysis revealed that the NAGP has captured allele frequency changes over time associated with selection, validating the strategy of repeated sampling and suggests that the continuation of a repeated sampling policy is essential for the germplasm collection to maintain its future utility. While NAGP should continue to collect bulls that have large influence on the population due to selection, care should be taken to include the entire breadth of bulls, including low merit bulls.

Historical museum samples reveal signals of selection and drift in response to changing insecticide use in an agricultural pest moth

In response to environmental and human-imposed selective pressures, agroecosystem pests frequently undergo rapid evolution, with some species having a remarkable capacity to rapidly develop pesticide resistance. Temporal sampling of genomic data can comprehensively capture such adaptive changes over time, for example, by elucidating allele frequency shifts in pesticide resistance loci in response to different pesticides. Here, we leveraged museum specimens spanning over a century of collections to generate temporal contrasts between pre- and post-insecticide populations of an agricultural pest moth, Helicoverpa armigera. We used targeted exon sequencing of 254 samples collected across Australia from the pre-1950s (prior to insecticide introduction) to the 1990s, encompassing decades of changing insecticide use. Our sequencing approach focused on genes that are known to be involved in insecticide resistance, environmental sensation, and stress tolerance. We found an overall lack of spatial and temporal population structure change across Australia. In some decades (e.g., 1960s and 1970s), we found a moderate reduction of genetic diversity, implying stochasticity in evolutionary trajectories due to genetic drift. Temporal genome scans showed extensive evidence of selection following insecticide use, although the majority of selected variants were low impact. Finally, alternating trajectories of allele frequency change were suggestive of potential antagonistic pleiotropy. Our results provide new insights into recent evolutionary responses in an agricultural pest and show how temporal contrasts using museum specimens can improve mechanistic understanding of rapid evolution.

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations

OBJECTIVE

The objective of this study was to identify genes associated with 305-day milk yield (MY) and fat yield (FY) that also influence the adaptability of the Thai multibreed dairy cattle population to tropical conditions.

METHODS

A total of 75,776 imputed and actual single nucleotide polymorphisms (SNPs) from 2,661 animals were used to identify genomic regions associated with MY and FY using the single-step genomic best linear unbiased predictions. Fixed effects included herd-yearseason, breed regression, heterosis regression and calving age regression effects. Random effects were animal additive genetic and residual. Individual SNPs with a p-value smaller than 0.05 were selected for gene mapping, function analysis, and quantitative trait loci (QTL) annotation analysis.

RESULTS

A substantial number of QTLs associated with MY (9,334) and FY (8,977) were identified by integrating SNP genotypes and QTL annotations. Notably, we discovered 17 annotated QTLs within the health and exterior QTL classes, corresponding to nine unique genes. Among these genes, Rho GTPase activating protein 15 (ARHGAP15) and catenin alpha 2 (CTNNA2) have previously been linked to physiological traits associated with tropical adaptation in various cattle breeds. Interestingly, these two genes also showed signs of positive selection, indicating their potential role in conferring tolerance to trypanosomiasis, a prevalent tropical disease.

CONCLUSION

Our findings provide valuable insights into the genetic basis of MY and FY in the Thai multibreed dairy cattle population, shedding light on the underlying mechanisms of tropical adaptation. The identified genes represent promising targets for future breeding strategies aimed at improving milk and fat production while ensuring resilience to tropical challenges. This study significantly contributes to our understanding of the genetic factors influencing milk production and adaptability in dairy cattle, facilitating the development of sustainable genetic selection strategies and breeding programs in tropical environments.

Genetic diversity and positive signatures of selection in indigenous cattle breeds of Iran

Animal domestication, climate changes over time, and artificial selection have played significant roles in shaping the genome structure of various animal species, including cattle. These processes have led to the emergence of several indigenous cattle breeds with distinct genetic characteristics. This study focused on unraveling the genetic diversity and identifying candidate genomic regions in eight indigenous cattle breeds of Iran. The data consisted of ∼777 962 single nucleotide polymorphisms (SNPs) of 89 animals from Iranian indigenous cattle scattered throughout the country. We employed various methods, including integrated haplotype score, FST, and cross-population composite likelihood ratio, to conduct a genome scan for detecting selection signals within and between cattle populations. Average observed heterozygosity across the populations was 0.36, with a range of 0.32-0.40. In addition, negative and low rates of inbreeding (FIS) in the populations were observed. The genome-wide analysis revealed several genomic regions that harbored candidate genes associated with production traits (e.g., MFSD1, TYW5, ADRB2, BLK, and CRTC3), adaptation to local environmental constraints (CACNA2D1, CXCL3, and GRO1), and coat color (DYM). Finally, the study of the reported quantitative trait loci (QTL) regions in the cattle genome demonstrated that the identified regions were associated with QTL related to important traits such as milk composition, body weight, daily gain, feed conversion, and residual feed intake. Overall, this study contributes to a better understanding of the genetic diversity and potential candidate genes underlying important traits in Iranian indigenous cattle breeds, which can inform future breeding and conservation efforts.

Selective Sweeps in Cattle Genomes in Response to the Influence of Urbanization and Environmental Contamination

A genomic study was conducted to identify the effects of urbanization and environmental contaminants with heavy metals on selection footprints in dairy cattle populations reared in the megacity of Bengaluru, South India. Dairy cattle reared along the rural-urban interface of Bengaluru with/without access to roughage from public lakeshores were selected. The genotyped animals were subjected to the cross-population-extended haplotype homozygosity (XP-EHH) methodology to infer selection sweeps caused by urbanization (rural, mixed, and urban) and environmental contamination with cadmium and lead. We postulated that social-ecological challenges contribute to mechanisms of natural selection. A number of selection sweeps were identified when comparing the genomes of cattle located in rural, mixed, or urban regions. The largest effects were identified on BTA21, displaying pronounced peaks for selection sweeps for all three urbanization levels (urban_vs_rural, urban_vs_mixed and rural_vs_mixed). Selection sweeps are located in chromosomal segments in close proximity to the genes lrand rab interactor 3 (RIN3), solute carrier family 24 member 4 (SLC24A4), tetraspanin 3 (TSPAN3), and proline-serine-threonine phosphatase interacting protein 1 (PSTPIP1). Functional enrichment analyses of the selection sweeps for all three comparisons revealed a number of gene ontology (GO) and KEGG terms, which were associated with reproduction, metabolism, and cell signaling-related functional mechanisms. Likewise, a number of the chromosomal segments under selection were observed when creating cattle groups according to cadmium and lead contaminations. Stronger and more intense positive selection sweeps were observed for the cadmium contaminated group, i.e., signals of selection on BTA 16 and BTA19 in close proximity to genes regulating the somatotropic axis (growth factor receptor bound protein 2 (GRB2) and cell ion exchange (chloride voltage-gated channel 6 (CLCN6)). A few novel, so far uncharacterized genes, mostly with effects on immune physiology, were identified. The lead contaminated group revealed sweeps which were annotated with genes involved in carcass traits (TNNC2, SLC12A5, and GABRA4), milk yield (HTR1D, SLCO3A1, TEK, and OPCML), reproduction (GABRA4), hypoxia/stress response (OPRD1 and KDR), cell adhesion (PCDHGC3), inflammatory response (ADORA2A), and immune defense mechanism (ALCAM). Thus, the findings from this study provide a deeper insight into the genomic regions under selection under the effects of urbanization and environmental contamination.

Population genomic structures and signatures of selection define the genetic uniqueness of several fancy and meat rabbit breeds

Following the recent domestication process of the European rabbit (Oryctolagus cuniculus), many different breeds and lines, distinguished primarily by exterior traits such as coat colour, fur structure and body size and shape, have been constituted. In this study, we genotyped, with a high-density single-nucleotide polymorphism panel, a total of 645 rabbits from 10 fancy breeds (Belgian Hare, Champagne d'Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex and Rhinelander) and three meat breeds (Italian White, Italian Spotted and Italian Silver). ADMIXTURE analysis indicated that breeds with similar phenotypic traits (e.g. coat colour and body size) shared common ancestries. Signatures of selection using two haplotype-based approaches (iHS and XP-EHH), combined with the results obtained with other methods previously reported that we applied to the same breeds, we identified a total of 5079 independent genomic regions with some signatures of selection, covering about 1777 Mb of the rabbit genome. These regions consistently encompassed many genes involved in pigmentation processes (ASIP, EDNRA, EDNRB, KIT, KITLG, MITF, OCA2, TYR and TYRP1), coat structure (LIPH) and body size, including two major genes (LCORL and HMGA2) among many others. This study revealed novel genomic regions under signatures of selection and further demonstrated that population structures and signatures of selection, left into the genome of these rabbit breeds, may contribute to understanding the genetic events that led to their constitution and the complex genetic mechanisms determining the broad phenotypic variability present in these untapped rabbit genetic resources.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.

Single-step genome-wide association studies and post-GWAS analyses for the number of oocytes and embryos in Gir cattle

Genome-Wide Association Studies (GWAS) are used for identification of quantitate trait loci (QTL) and genes associated with several traits. We aimed to identify genomic regions, genes, and biological processes associated with number of total and viable oocytes, and number of embryos in Gir dairy cattle. A dataset with 17,526 follicular aspirations, including the following traits: number of viable oocytes (VO), number of total oocytes (TO), and number of embryos (EMBR) from 1641 Gir donors was provided by five different stock farms. A genotype file with 2093 animals and 395,524 SNP markers was used to perform a single-step GWAS analysis for each trait. The top 10 windows with the highest percentage of additive genetic variance explained by 100 adjacent SNPs were selected. The genomic regions identified in our work were overlapped with QTLs from QTL database on chromosomes 1, 2, 5, 6, 7, 8, 9, 13, 17, 18, 20, 21, 22, 24, and 29. These QTLs were classified as External, Health, Meat and carcass, Production or Reproduction traits, and about 38% were related to Reproduction. In total, 117 genes were identified, of which 111 were protein-coding genes. Exclusively associations were observed for 42 genes with EMBR, and 1 with TO. Also, 42 genes were in common between VO and TO, 28 between VO and EMBR, and four genes were in common among all traits. In conclusion, great part of the identified genes plays a functional role in initial embryo development or general cell functions. The protein-coding genes ARNT, EGR1, HIF1A, AHR, and PAX2 are good markers for the production of oocytes and embryos in Gir cattle.

Dual-Purpose Cattle Raised in Tropical Conditions: What Are Their Shortcomings in Sound Productive and Reproductive Function?

Dual-purpose husbandry might well be the most commonly employed cattle management system in tropical regions worldwide. The advantages of producing both meat and milk, although in reduced quantities, gives an edge to the farmer in coping with the volatile economic conditions that prevail in the region. Herein, we discuss the different methods of cattle management under tropical conditions based on the financial and social structure of this system. An account of the sanitary and nutritional conditions available to the farmers and how these factors affect the profitability of the enterprise will also be given. Finally, we will discuss how these systems can take advantage of several biotechnological procedures, and how these tools (such as controlled natural mating, artificial insemination, and embryo transfer) affect reproductive outcomes. The present review will mainly concentrate on production systems located less than 1000 m above sea level, as the problems and shortcomings of cattle raised above this arbitrary landmark are quite different.

Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle

Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (F_ROH) and pedigree-based inbreeding (F_ped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH_1-2Mb: 58.17%), representing ancient inbreeding, ROH_2-4Mb (22.74%), ROH_4-8Mb (11.34%), ROH_8-16Mb (5.51%), and ROH_>16Mb (2.24%). Combining our results, we conclude that the increase in general F_ROH and F_ped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with F_ROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.

Genetic Architecture and Signatures of Selection in the Caqueteño Creole (Colombian Native Cattle)

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.

Genomic breed composition of pure registered Brazilian Gir

INTRODUCTION

Gir is a Bos indicus breed originally from India, first imported to Brazil in 1850. High-performance Dairy Gir has been systematically selected in Brazil from its arrival. Since the major phenotypic difference between Gir in India and Brazil is a higher milk production in the former, it is speculated that Brazilian Gir was strategically crossed with Holstein or another Bos taurus to improve milk yield. This study evaluated the purity of Brazilian Gir breed stocks from BASA Farms in Brazil, trying to identify possible admixture events with other cattle breeds based on DNA analysis.

MATERIALS AND METHODS

The population included 1061 pure registered individuals genotyped using two commercial platforms with 37 k and 25 k SNPs. Admixture analysis was performed individually to estimate levels of genomic composition derived from six different reference populations, three indicine and three taurine breeds.

RESULTS

A Gir ancestry of 99% or higher was found for 94.2% of the population, while the remaining showed levels of non-Gir ancestry up to 6.8%. Only five individuals were identified with possible taurine ancestry, all of them exhibiting levels lower than 2%. The remaining non-Gir ancestry identified was derived from indicine breeds. The levels of admixture observed in the population were from low to non-detectable. No consistent patterns of admixture were observed indicating sustained introgression of taurine lines as means of genetic improvement.

CONCLUSION

According to these results, genetic improvement achieved by Brazilian Gir breeders is the result of within-breed selection methods applied intensively over the past five decades, rather than the result from sustained introgression.

Genome-wide scan for signatures of selection in the Brangus cattle genome

Brangus is a composite cattle breed developed with the objective of combining the advantages of Angus and Zebuine breeds (Brahman, mainly) in tropical climates. The aim of this work was to estimate breed composition both genome-wide and locally, at the chromosome level, and to uncover genomic regions evidencing positive selection in the Argentinean Brangus population/nucleus. To do so, we analysed marker data from 478 animals, including Brangus, Angus and Brahman. Average breed composition was 35.0% ± 9.6% of Brahman, lower than expected according to the theoretical fractions deduced by the usual cross-breeding practice in this breed. Local ancestry analysis evidenced that breed composition varies between chromosomes, ranging from 19.6% for BTA26 to 56.1% for BTA5. Using approaches based on allelic frequencies and linkage disequilibrium, genomic regions with putative selection signatures were identified in several chromosomes (BTA1, BTA5, BTA6 and BTA14). These regions harbour genes involved in horn development, growth, lipid metabolism, reproduction and immune response. We argue that the overlapping of a chromosome segment originated in one of the parental breeds and over-represented in the sample with the location of a signature of selection constitutes evidence of a selection process that has occurred in the breed since its take off in the 1950s. In this regard, our results could contribute to the understanding of the genetic mechanisms involved in cross-bred cattle adaptation and productivity in tropical environments.

Selection Signatures in South African Nguni and Bonsmara Cattle Populations Reveal Genes Relating to Environmental Adaptation

Climate change is a major influencing factor in beef production. The greenhouse gases produced from livestock production systems contribute to the overall greenhouse gas emissions. The aim of this study was to identify selection signatures within and between Nguni and Bonsmara cattle in relation to production and adaptation. For this purpose, genomic 150 K single nucleotide polymorphism data from Nguni (n = 231) and Bonsmara (n = 252) cattle in South Africa were used. Extended haplotype homozygosity (EHH) based analysis was executed within each population using integrated haplotype score (iHS). The R package rehh was used for detecting selection signatures across the two populations with cross population EHH (XP-EHH). Total of 121 regions of selection signatures were detected (p < 0.0001) in the Bonsmara and Nguni populations. Several genes relating to DNA methylation, heat stress, feed efficiency and nitrogen metabolism were detected within and between each population. These regions also included QTLs associated with residual feed intake, residual gain, carcass weight, stature and body weight in the Bonsmara, while QTLs associated with conception rate, shear force, tenderness score, juiciness, temperament, heat tolerance, feed efficiency and age at puberty were identified in Nguni. Based on the results of the study it is recommended that the Nguni and Bonsmara be utilized in crossbreeding programs as they have beneficial traits that may allow them to perform better in the presence of climate change. Results of this study coincide with Nguni and Bonsmara breed characteristics and performance, and furthermore support informative crossbreeding programs to enhance livestock productivity in South Africa.

Comparative genomics and signatures of selection in North Atlantic eels

Comparative genomic approaches can identify putative private and shared signatures of selection. We performed a comparative genomic study of North Atlantic eels, European eel (Anguilla Anguilla) and American eel (A. rostrata). The two sister species are nearly undistinguishable at the phenotypic level and despite a wide non-overlapping continental distribution, they spawn in partial sympatry in the Sargasso Sea. Taking advantage of the newly assembled and annotated genome, we used genome wide RAD sequencing data of 359 individuals retrieved from Sequence Nucleotide Archive and state-of-the-art statistic tests to identify putative genomic signatures of selection in North Atlantic eels. First, using the F_ST and XP-EHH methods, we detected apparent islands of divergence on a total of 7 chromosomes, particularly on chromosomes 6 and 10. Gene ontology analyses suggested candidate genes mainly related to energy production, development and regulation, which could reflect strong selection on traits related to eel migration and larval duration time. Gene effect prediction using SNPeff showed a high number of SNPs in noncoding regions, pointing to a possible regulatory role. Second, using the iHS method we detected shared regions under selection on a total of 11 chromosomes. Several hypotheses might account for the detection of shared islands of selection in North Atlantic eels, including parallel evolution due to adaptation to similar environments and introgression. Future comparative genomic studies will be needed to further clarify the causes and consequences of introgression, including the directionality of these introgression events.

Unraveling Admixture, Inbreeding, and Recent Selection Signatures in West African Indigenous Cattle Populations in Benin

The Dwarf Lagune and the Savannah Somba cattle in Benin are typical representatives of the endangered West African indigenous Shorthorn taurine. The Lagune was previously exported to African and European countries and bred as Dahomey cattle, whereas the Somba contributed to the formation of two indigenous hybrids known as Borgou and Pabli cattle. These breeds are affected by demographic, economic, and environmental pressures in local production systems. Considering current and historical genomic data, we applied a formal test of admixture, estimated admixture proportions, and computed genomic inbreeding coefficients to characterize the five breeds. Subsequently, we unraveled the most recent selection signatures using the cross-population extended haplotype homozygosity approach, based on the current and historical genotypes. Results from principal component analyses and high proportion of Lagune ancestry confirm the Lagune origin of the European Dahomey cattle. Moreover, the Dahomey cattle displayed neither indicine nor European taurine (EUT) background, but they shared on average 40% of autozygosity from common ancestors, dated approximately eight generations ago. The Lagune cattle presented inbreeding coefficients larger than 0.13; however, the Somba and the hybrids (Borgou and Pabli) were less inbred (≤0.08). We detected evidence of admixture in the Somba and Lagune cattle, but they exhibited a similar African taurine (AFT) ancestral proportion (≥96%) to historical populations, respectively. A moderate and stable AFT ancestral proportion (62%) was also inferred for less admixed hybrid cattle including the Pabli. In contrast, the current Borgou samples displayed a lower AFT ancestral proportion (47%) than historical samples (63%). Irrespective of the admixture proportions, the hybrid populations displayed more selection signatures related to economic traits (reproduction, growth, and milk) than the taurine. In contrast, the taurine, especially the Somba, presented several regions known to be associated with adaptive traits (immunity and feed efficiency). The identified subregion of bovine leukocyte antigen (BoLA) class IIb (including DSB and BOLA-DYA) in Somba cattle is interestingly uncommon in other African breeds, suggesting further investigations to understand its association with specific adaptation to endemic diseases in Benin. Overall, our study provides deeper insights into recent evolutionary processes in the Beninese indigenous cattle and their aptitude for conservation and genetic improvement.

Assessing the Genetic Background and Selection Signatures of Huaxi Cattle Using High-Density SNP Array

Huaxi cattle, a specialized beef cattle breed in China, has the characteristics of fast growth, high slaughter rate, and net meat rate, good reproductive performance, strong stress resistance, and wide adaptability. In this study, we evaluated the genetic diversity, population structure, and genetic relationships of Huaxi cattle and its ancestor populations at the genome-wide level, as well as detecting the selection signatures of Huaxi cattle. Principal component analysis (PCA) and phylogenetic analysis revealed that Huaxi cattle were obviously separated from other cattle populations. The admixture analysis showed that Huaxi cattle has distinct genetic structures among all populations at K = 4. It can be concluded that Huaxi cattle has formed its own unique genetic features. Using integrated haplotype score (iHS) and composite likelihood ratio (CLR) methods, we identified 143 and 199 potentially selected genes in Huaxi cattle, respectively, among which nine selected genes (KCNK1, PDLIM5, CPXM2, CAPN14, MIR2285D, MYOF, PKDCC, FOXN3, and EHD3) related to ion binding, muscle growth and differentiation, and immunity were detected by both methods. Our study sheds light on the unique genetic feature and phylogenetic relationship of Huaxi cattle, provides a basis for the genetic mechanism analysis of important economic traits, and guides further intensive breeding improvement of Huaxi cattle.

Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains

Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.

Selection Response Due to Different Combination of Antagonistic Milk, Beef, and Morphological Traits in the Alpine Grey Cattle Breed

Selection in local dual-purpose breeds requires great carefulness because of the need to preserve peculiar traits and also guarantee the positive genetic progress for milk and beef production to maintain economic competitiveness. A specific breeding plan accounting for milk, beef, and functional traits is required by breeders of the Alpine Grey cattle (AG), a local dual-purpose breed of the Italian Alps. Hereditability and genetic correlations among all traits have been analyzed for this purpose. After that, different selection indexes were proposed to identify the most suitable for this breed. Firstly, a genetic parameters analysis was carried out with different datasets. The milk dataset contained 406,918 test day records of milk, protein, and fat yields and somatic cells (expressed as SCS). The beef dataset included performance test data conducted on 749 young bulls. Average daily gain, in vivo estimated carcass yields, and carcass conformation (SEUROP) were the phenotypes obtained from the performance tests. The morphological dataset included 21 linear type evaluations of 11,320 first party cows. Linear type traits were aggregated through factor analysis and three factors were retained, while head typicality (HT) and rear muscularity (RM) were analyzed as single traits. Heritability estimates (h2) for milk traits ranged from 0.125 to 0.219. Analysis of beef traits showed h2 greater than milk traits, ranging from 0.282 to 0.501. Type traits showed a medium value of h2 ranging from 0.238 to 0.374. Regarding genetic correlation, SCS and milk traits were strongly positively correlated. Milk traits had a negative genetic correlation with the factor accounting for udder conformations (-0.40) and with all performance test traits and RM. These latter traits showed also a negative genetic correlation with udder volume (-0.28). The HT and the factor accounting for rear legs traits were not correlated with milk traits, but negatively correlated with beef traits (-0.32 with RM). We argue that the consequence of these results is that the use of the current selection index, which is mainly focused on milk attitude, will lead to a deterioration of all other traits. In this study, we propose more appropriate selection indexes that account for genetic relationships among traits, including functional traits.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

Background

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne’s disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (−); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip.

Results

The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/−) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10^− 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids.

Conclusion

In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07487-4.

Genomic scans for selection signatures revealed candidate genes for adaptation and production traits in a variety of cattle breeds

Domestication and selection are the major driving forces responsible for the determinative genetic variability in livestock. These selection patterns create unique genetic signatures within the genome. BovineSNP50 chip data from 236 animals (seven indicine and five taurine cattle breeds) were analyzed in the present study. We implemented three complementary approaches viz. iHS (Integrated haplotype score), ROH (Runs of homozygosity), and F_ST, to detect selection signatures. A total of 179, 56, and 231 regions revealed 518, 277, and 267 candidate genes identified by iHS, ROH, and F_ST methods, respectively. We found several candidate genes (e.g., NCR3, ARID5A, HIST1H2BN, DEFB4, DEFB7, HSPA1L, HSPA1B, and DNAJB4) related to production traits and the adaptation of indigenous breeds to local environmental constraints such as heat stress and disease susceptibility. However, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping.

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.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genetic characterization and population structure of different coat colour variants of Badri cattle

The present study aimed to genetically characterize the Badri cattle and its three colour variants and assess their population structure using 24 microsatellite markers. Out of 96 animals analyzed, 32 each were collected from grey (GVBC), brown (BrVBC) and black (BVBC) colour variants of Badri cattle (BC). The genetic diversity parameters including allele frequencies, observed and effective number of alleles, observed and expected heterozygosity, PIC, Shannon's indices and F-statistics were estimated using POPGENE software. Bottleneck analysis was performed using both qualitative and quantitative approaches. A total of 274 alleles (50 private and 224 shared) were scored for BC, GVBC, BrVBC and BVBC with mean number of 11.417, 9.083, 9.125 and 9.083 alleles, respectively. All populations exhibited average heterozygosity estimate > 0.5 indicating existence of substantial genetic variability, concurrent with revelations from Shannon's indices. Observed mean PIC estimates (> 0.74) were indicative of optimum informativeness of used microsatellite markers. The mean inbreeding estimates (F) in GVBC, BrVBC and BVBC were 0.041, - 0.024 and 0.016, respectively. The pair wise genetic (> 0.91) pointed towards similarity between different colour variant populations. STRUCTURE analysis also revealed clear admixture for the three Badri colour variants indicating absence of genetic differentiation. The present study revealed first-hand information that populations of Badri cattle with different phenotypes with respect to coat colour are genetically related and can be considered as a single breed. The comprehensive knowledge generated for Badri cattle will help in designing breeding plan for its genetic improvement and deciding the conservation priorities.

Genome-Wide Population Genetic Analysis of Commercial, Indigenous, Game, and Wild Chickens Using 600K SNP Microarray Data

Following chicken domestication, diversified chicken breeds were developed by both natural and artificial selection, which led to the accumulation of abundant genetic and phenotypic variations, making chickens an ideal genetic research model. To better understand the genetic structure of chicken breeds under different selection pressures, we genotyped various chicken populations with specific selection targets, including indigenous, commercial, gamecock, and wild ancestral chickens, using the 600K SNP array. We analyzed the population structure, genetic relationships, run of homozygosity (ROH), effective population number (Ne), and other genetic parameters. The wild ancestral population, red junglefowl (RJF), possessed the highest diversity, in comparison with all other domesticated populations, which was supported by linkage disequilibrium decay (LD), effective population number, and ROH analyses. The gamecock breeds, which were subjected to stronger male-biased selection for fighting-related traits, also presented higher variation than the commercial and indigenous breeds. Admixture analysis also indicated that game breed is a relatively independent branch of Chinese local breeds. Following intense selection for reproductive and productive traits, the commercial lines showed the least diversity. We also observed that the European local chickens had lower genetic variation than the Chinese local breeds, which could be attributed to the shorter history of the European breed. ROH were present in a breed specific manner and 191 ROH island were detected on four groups (commercial, local, game and wild chickens). These ROH islands were involved in egg production, growth and silky feathers and other traits. Moreover, we estimated the effective sex ratio of these breeds to demonstrate the change in the ratio of the two sexes. We found that commercial chickens had a greater sex imbalance between females and males. The commercial lines showed the highest female-to-male ratios. Interestingly, RJF comprised a greater proportion of males than females. Our results show the population genetics of chickens under selection pressures, and can aid in the development of better conservation strategies for different chicken breeds.

Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color

Coat color is among the most distinctive phenotypes in cattle. Worldwide, several breeds share peculiar coat color features such as the presence of a fawn pigmentation of the calf at birth, turning over time to grey, and sexual dichromatism. The aim of this study was to search for polymorphisms under differential selection by contrasting grey cattle breeds displaying the above phenotype with non-grey cattle breeds, and to identify the underlying genes. Using medium-density SNP array genotype data, a multi-cohort F_ST-outlier approach was adopted for a total of 60 pair-wise comparisons of the 15 grey with 4 non-grey cattle breeds (Angus, Limousin, Charolais, and Holstein), with the latter selected as representative of solid and piebald phenotypes, respectively. Overall, more than 50 candidate genes were detected; almost all were either directly or indirectly involved in pigmentation, and some of them were already known for their role in phenotypes related with hair graying in mammals. Notably, 17 relevant genes, including SDR16C5, MOS, SDCBP, and NSMAF, were located in a signal on BTA14 convergently observed in all the four considered scenarios. Overall, the key stages of pigmentation (melanocyte development, melanogenesis, and pigment trafficking/transfer) were all represented among the pleiotropic functions of the candidate genes, suggesting the complex nature of the grey phenotype in cattle.

Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus)

Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one F_ST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and F_ST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.

Genetic Diversity and Signatures of Selection in a Native Italian Horse Breed Based on SNP Data

Simple Summary

The Bardigiano horse is a native Italian breed bred for living in rural areas, traditionally used in agriculture. The breed counts about 3000 horses, and it is nowadays mainly used for recreational purposes. The relatively small size and the closed status of the breed raise the issue of monitoring genetic diversity. We therefore characterized the breed’s genetic diversity based on molecular data. We showed a critical reduction of genetic variability mainly driven by past bottlenecks. We also highlighted homozygous genomic regions that might be the outcome of directional selection in recent years, in line with the conversion of Bardigiano horses from agricultural to riding purposes.

Abstract

Horses are nowadays mainly used for sport and leisure activities, and several local breeds, traditionally used in agriculture, have been exposed to a dramatic loss in population size and genetic diversity. The loss of genetic diversity negatively impacts individual fitness and reduces the potential long-term survivability of a breed. Recent advances in molecular biology and bioinformatics have allowed researchers to explore biodiversity one step further. This study aimed to evaluate the loss of genetic variability and identify genomic regions under selection pressure in the Bardigiano breed based on GGP Equine70k SNP data. The effective population size based on Linkage Disequilibrium (N_e) was equal to 39 horses, and it showed a decline over time. The average inbreeding based on runs of homozygosity (ROH) was equal to 0.17 (SD = 0.03). The majority of the ROH were relatively short (91% were ≤2 Mbp long), highlighting the occurrence of older inbreeding, rather than a more recent occurrence. A total of eight ROH islands, shared among more than 70% of the Bardigiano horses, were found. Four of them mapped to known quantitative trait loci related to morphological traits (e.g., body size and coat color) and disease susceptibility. This study provided the first genome-wide scan of genetic diversity and selection signatures in an Italian native horse breed.

Effect direction meta-analysis of GWAS identifies extreme, prevalent and shared pleiotropy in a large mammal

In genome-wide association studies (GWAS), variants showing consistent effect directions across populations are considered as true discoveries. We model this information in an Effect Direction MEta-analysis (EDME) to quantify pleiotropy using GWAS of 34 Cholesky-decorrelated traits in 44,000+ cattle with sequence variants. The effect-direction agreement between independent bull and cow datasets was used to quantify the false discovery rate by effect direction (FDRed) and the number of affected traits for prioritised variants. Variants with multi-trait p < 1e–6 affected 1∼22 traits with an average of 10 traits. EDME assigns pleiotropic variants to each trait which informs the biology behind complex traits. New pleiotropic loci are identified, including signals from the cattle FTO locus mirroring its bystander effects on human obesity. When validated in the 1000-Bull Genome database, the prioritized pleiotropic variants consistently predicted expected phenotypic differences between dairy and beef cattle. EDME provides robust approaches to control GWAS FDR and quantify pleiotropy.

Xiang et al. developed an Effect Direction Meta-analysis (EDME) approach to identify true pleiotropy. They used Cholesky-transformation to decorrelate the traits and identified many pleiotropic variants that consistently predicted phenotypic differences in cattle.