Signatures of selection in Nelore cattle revealed by whole-genome sequencing data

Genome-wide analysis of genetic diversity and selection signatures in Fuzhou cattle

The Fuzhou cattle breed, native to northeast China, is widely recognized for its adaptability, disease resistance, and docility. Despite being known for these qualities, its population has declined recently, and there is a significant lack of genomic studies on this species. We sequenced 21 samples from a primary breeding farm to determine the genetic structure, diversity, and selection signature to address this. Additionally, we combined 100 published genomic datasets from diverse geographical regions to characterize the genomic variation of Fuzhou cattle. There were 53 752 978 bi-allelic SNPs retained for downstream analysis. In population structure analysis, Fuzhou cattle show a predominantly East Asian taurine ancestry, with strong genetic affinities to Hanwoo and Yanbian cattle. Despite high nucleotide diversity within the Bos taurine lineage, genetic diversity analysis also revealed significant levels of inbreeding in Fuzhou cattle populations, indicating the need for conservation. Utilizing various methods such as θπ, iHS, FST, π-ratio, and XP-EHH, we identified genes associated with traits like growth, meat quality, energy metabolism, and immunity. Several genes related to cold adaptation were identified, including PLIN5, PLB1, and CPT2. These findings provide a basis for conservation strategies to safeguard the genetic resources of Fuzhou cattle.

Genetic parameters, genomic prediction, and identification of regulatory regions located on chromosome 14 for weight traits in Nellore cattle

This study aimed to investigate functional variants in chromosome 14 (BTA14) and its impact in genomic selection for birth weight (BW), weaning weight (WW), and yearling weight (YW) in Nellore cattle. Genetic parameter estimation and the weighted single-step genomic best linear unbiased prediction (WssGBLUP) analyses were performed. Direct additive heritability estimates were high for WW and YW, and moderate for BW. Trait-associated variants distributed across multiple regions on BTA14 were observed in the weighted single-step genome-wide association studies (WssGWAS) results, implying a polygenic genetic architecture for weight in different ages. Several genes have been found in association with the weight traits, including the CUB And Sushi multiple domains 3 (CSMD3), thyroglobulin (TG), and diacylglycerol O-acyltransferase 1 (DGAT1) genes. The variance explained per SNP was higher in six functional classes of gene regulatory regions (5UTR, CpG islands, downstream, upstream, long non-coding RNA, and transcription factor binding sites (TFBS)), highlighting their importance for weight traits in Nellore cattle. A marginal increase in accuracy was observed when the selected functional variants (SV) information was considered in the WssGBLUP method, probably because of the small number of SV available on BTA14. The identified genes, pathways, and functions contribute to a better understanding of the genetic and physiological mechanisms regulating weight traits in the Nellore breed.

Whole genome resequencing reveals the correlation between selection signatures and adaptability of Micropterus salmoides to artificial fed

Largemouth bass (Micropterus salmoides, LMB) is an important aquaculture species due to its excellent flesh quality and environmental adaptability. It has been continuously introduced to many countries and cultured for decades. Here, an LMB population was used for selective breeding to improve growth rate and feed adaptability. After five generations of breeding, the growth rate improved by 38%, and feed adaptability improved by 22% compared to the non-breeding population. To study the underlying genetic mechanism, 100 LMB from the breeding population and 100 from the non-breeding population were sampled for whole-genome resequencing. The population genetics analysis shows that the breeding population has a higher inbreeding coefficient and linkage disequilibrium (LD) level, a lower nucleic acid diversity and effective population size (Ne). Using [Formula: see text] (fixation index), we found that the average [Formula: see text] value between the two populations was 0.07, with the highest [Formula: see text] value reaching 0.38, which overlaps with the trypsin gene. Additionally, other genes exhibiting high [Formula: see text] values are associated with functions such as neural development, glucose metabolism, and growth. Using [Formula: see text] and nucleic acid diversity as criteria, we identified 698 genes that are positively selected in the breeding population, and gene functional enrichment analysis shows that 36 genes are related to the olfactory receptor pathway. Overall, our study found that multiple genes were selected in the LMB breeding population. These genes may be associated with adaptation and digestion of artificial feed in fish.

Graduate Student Literature Review: The DRB3 gene of the bovine major histocompatibility complex-Discovery, diversity, and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine major histocompatibility complex (MHC), also known as the bovine leukocyte antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. The BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and Creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

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 Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

Whole-genome resequencing deciphers patterns of genetic diversity, phylogeny, and evolutionary dynamics in Kashmir cattle

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

Genome-wide association study of Nelore and Angus heifers with low and high ovarian follicle counts

The number of antral follicles is considered an important fertility trait because animals with a high follicle count (HFC) produce more oocytes and embryos per cycle. Identification of these animals by genetic markers such as single nucleotide polymorphisms (SNPs) can accelerate selection of future generations. The aim of this study was to perform a genome wide association study (GWAS) on Nelore and Angus heifers with HFC and low (LFC) antral follicle counts. The groups HFC and LFC for genotyping were formed based on the average of total follicles (≥ 3 mm) counted in each breed consistently ± standard deviation. A total of 72 Nelore heifers (32 HFC and 40 LFC) and 48 Angus heifers (21 HFC and 27 LFC) were selected and the DNA was extracted from blood and hair bulb. Genotyping was done using the Illumina Bovine HD 770K BeadChip. The GWAS analysis showed 181 and 201 SNPs with genotype/phenotype association (P ≤ 0.01) in Nelore and Angus heifers, respectively. Functional enrichment analysis was performed on candidate genes that were associated with SNPs. A total of 97 genes were associated to the 181 SNPs in the Nelore heifers and the functional analysis identified genes (ROBO1 and SLIT3) in the ROBO-SLIT pathway that can be involved in the control of germ cell migration in the ovary as it is involved in lutheal cell migration and fetal ovary development. In the Angus heifers, 57 genes were associated with the 201 SNPs, highlighting Fribilin 1 (FBN1) gene, involved in regulation of growth factors directly involved in follicle activation and development. In summary, GWAS for Nelore and Angus heifers showed SNPs associated with higher follicle count phenotype. Furthermore, these findings offer valuable insights for the further investigation of potential mechanism involved in follicle formation and development, important for breeding programs for both breeds.

Evidence for selective sweeps in the MHC gene repertoire of various cattle breeds

Major Histocompatibility Complex (MHC) genes are among the immune genes that have been extensively studied in vertebrates and are necessary for adaptive immunity. In the immunological response to infectious diseases, they play several significant roles. This research paper provides the selection signatures in the MHC region of the bovine genome as well as how certain genes related to innate immunity are undergoing a positive selective sweep. Here, we investigated signatures of historical selection on MHC genes in 15 different cattle populations and a total of 427 individuals. To identify the selection signatures, we have used three separate summary statistics. The findings show potential selection signatures in cattle from whom we isolated genes involved in the MHC. The most significant regions related to the bovine MHC are BOLA, non-classical MHC class I antigen (BOLA-NC1), Microneme protein 1 (MIC1) , Cluster of Differentiation 244 (CD244), Gap Junction Alpha-5 Protein (GJA5). It will be possible to gain new insight into immune system evolution by understanding the distinctive characteristics of MHC in cattle.

Sequence level genome-wide associations for bull production and fertility traits in tropically adapted bulls

BACKGROUND

The genetics of male fertility is complex and not fully understood. Male subfertility can adversely affect the economics of livestock production. For example, inadvertently mating bulls with poor fertility can result in reduced annual liveweight production and suboptimal husbandry management. Fertility traits, such as scrotal circumference and semen quality are commonly used to select bulls before mating and can be targeted in genomic studies. In this study, we conducted genome-wide association analyses using sequence-level data targeting seven bull production and fertility traits measured in a multi-breed population of 6,422 tropically adapted bulls. The beef bull production and fertility traits included body weight (Weight), body condition score (CS), scrotal circumference (SC), sheath score (Sheath), percentage of normal spermatozoa (PNS), percentage of spermatozoa with mid-piece abnormalities (MP) and percentage of spermatozoa with proximal droplets (PD).

RESULTS

After quality control, 13,398,171 polymorphisms were tested for their associations with each trait in a mixed-model approach, fitting a multi-breed genomic relationship matrix. A Bonferroni genome-wide significance threshold of 5 × 10^- 8 was imposed. This effort led to identifying genetic variants and candidate genes underpinning bull fertility and production traits. Genetic variants in Bos taurus autosome (BTA) 5 were associated with SC, Sheath, PNS, PD and MP. Whereas chromosome X was significant for SC, PNS, and PD. The traits we studied are highly polygenic and had significant results across the genome (BTA 1, 2, 4, 6, 7, 8, 11, 12, 14, 16, 18, 19, 23, 28, and 29). We also highlighted potential high-impact variants and candidate genes associated with Scrotal Circumference (SC) and Sheath Score (Sheath), which warrants further investigation in future studies.

CONCLUSION

The work presented here is a step closer to identifying molecular mechanisms that underpin bull fertility and production. Our work also emphasises the importance of including the X chromosome in genomic analyses. Future research aims to investigate potential causative variants and genes in downstream analyses.

Evaluation of the Genetic Diversity, Population Structure and Selection Signatures of Three Native Chinese Pig Populations

Indigenous pig populations in Hainan Province live in tropical climate conditions and a relatively closed geographical environment, which has contributed to the formation of some excellent characteristics, such as heat tolerance, strong disease resistance and excellent meat quality. Over the past few decades, the number of these pig populations has decreased sharply, largely due to a decrease in growth rate and poor lean meat percentage. For effective conservation of these genetic resources (such as heat tolerance, meat quality and disease resistance), the whole-genome sequencing data of 78 individuals from 3 native Chinese pig populations, including Wuzhishan (WZS), Tunchang (TC) and Dingan (DA), were obtained using a 150 bp paired-end platform, and 25 individuals from two foreign breeds, including Landrace (LR) and Large White (LW), were downloaded from a public database. A total of 28,384,282 SNPs were identified, of which 27,134,233 SNPs were identified in native Chinese pig populations. Both genetic diversity statistics and linkage disequilibrium (LD) analysis indicated that indigenous pig populations displayed high genetic diversity. The result of population structure implied the uniqueness of each native Chinese pig population. The selection signatures were detected between indigenous pig populations and foreign breeds by using the population differentiation index (F_ST) method. A total of 359 candidate genes were identified, and some genes may affect characteristics such as immunity (IL-2, IL-21 and ZFYVE16), adaptability (APBA1), reproduction (FGF2, RNF17, ADAD1 and HIPK4), meat quality (ABCA1, ADIG, TLE4 and IRX5), and heat tolerance (VPS13A, HSPA4). Overall, the findings of this study will provide some valuable insights for the future breeding, conservation and utilization of these three Chinese indigenous pig populations.

Genome-wide scans for selection signatures in indigenous pigs revealed candidate genes relating to heat tolerance

Heat stress is a major problem that constrains pig productivity. Understanding and identifying adaptation to heat stress has been the focus of recent studies, and the identification of genome-wide selection signatures can provide insights into the mechanisms of environmental adaptation. Here, we generated whole-genome re-sequencing data from six Chinese indigenous pig populations to identify genomic regions with selection signatures related to heat tolerance using multiple methods: three methods for intra-population analyses (Integrated Haplotype Score, Runs of Homozygosity and Nucleotide diversity Analysis) and three methods for inter-population analyses (Fixation index (F_ST), Cross-population Composite Likelihood Ratio and Cross-population Extended Haplotype Homozygosity). In total, 1 966 796 single nucleotide polymorphisms were identified in this study. Genetic structure analyses and F_ST indicated differentiation among these breeds. Based on information on the location environment, the six breeds were divided into heat and cold groups. By combining two or more approaches for selection signatures, outlier signals in overlapping regions were identified as candidate selection regions. A total of 163 candidate genes were identified, of which, 29 were associated with heat stress injury and anti-inflammatory effects. These candidate genes were further associated with 78 Gene Ontology functional terms and 30 Kyoto Encyclopedia of Genes and Genomes pathways in enrichment analysis (P < 0.05). Some of these have clear relevance to heat resistance, such as the AMPK signalling pathway and the mTOR signalling pathway. The results improve our understanding of the selection mechanisms responsible for heat resistance in pigs and provide new insights of introgression in heat adaptation.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

The study of selection signature and its applications on identification of candidate genes using whole genome sequencing data in chicken - a review

Chicken is a major source of protein for the increasing human population and is useful for research purposes. There are almost 1,600 distinct regional breeds of chicken across the globe, among which a large body of genetic and phenotypic variations has been accumulated due to extensive natural and artificial selection. Moreover, natural selection is a crucial force for animal domestication. Several approaches have been adopted to detect selection signatures in different breeds of chicken using whole genome sequencing (WGS) data including integrated haplotype score (iHS), cross-populated extend haplotype homozygosity test (XP-EHH), fixation index (F_ST), cross-population composite likelihood ratio (XP-CLR), nucleotide diversity (Pi), and others. In addition, gene enrichment analyses are utilized to determine KEGG pathways and gene ontology (GO) terms related to traits of interest in chicken. Herein, we review different studies that have adopted diverse approaches to detect selection signatures in different breeds of chicken. This review systematically summarizes different findings on selection signatures and related candidate genes in chickens. Future studies could combine different selection signatures approaches to strengthen the quality of the results thereby providing more affirmative inference. This would further aid in deciphering the importance of selection in chicken conservation for the increasing human population.

Identification of copy number variations in the genome of Dairy Gir cattle

Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.

Towards Sustainable Sources of Omega-3 Long-Chain Polyunsaturated Fatty Acids in Northern Australian Tropical Crossbred Beef Steers through Single Nucleotide Polymorphisms in Lipogenic Genes for Meat Eating Quality

This study aimed to identify single nucleotide polymorphisms (SNP) in lipogenic genes of northern Australian tropically adapted crossbred beef cattle and to evaluate associations with healthy lipid traits of the Longissimus dorsi (loin eye) muscle. The hypothesis tested was that there are significant associations between SNP loci encoding for the fatty acid binding protein 4 (FABP4), stearoyl-CoA desaturase (SCD) and fatty acid synthase (FASN) genes and human health beneficial omega-3 long-chain polyunsaturated fatty acids (ω3 LC-PUFA) within the loin eye muscle of northern Australian crossbred beef cattle. Brahman, Charbray, and Droughtmaster crossbred steers were fed on Rhodes grass hay augmented with desmanthus, lucerne, or both, for 140 days and the loin eye muscle sampled for intramuscular fat (IMF), fat melting point (FMP), and fatty acid composition. Polymorphisms in FABP4, SCD, and FASN genes with significant effects on lipid traits were identified with next-generation sequencing. The GG genotype at the FABP4 g.44677239C>G locus was associated with higher proportion of linoleic acid than the CC and CG genotypes (p < 0.05). Multiple comparisons of genotypes at the SCD g.21266629G>T locus indicated that the TT genotype had significantly higher eicosapentaenoic, docosapentaenoic, and docosahexaenoic acids than GG genotype (p < 0.05). Significant correlations (p < 0.05) between FASN SNP and IMF, saturated and monounsaturated fatty acids were observed. These results provide insights into the contribution of lipogenic genes to intramuscular fat deposition and SNP marker-assisted selection for improvement of meat-eating quality, with emphasis on alternate and sustainable sources of ω3 LC-PUFA, in northern Australian tropical crossbred beef cattle, hence an acceptance of the tested hypothesis.