Bovine Genome Research Articles

Detection of putative loci affecting milk yield in Turkish Awassi sheep using microsatellite markers.

On the basis of comparisons between bovine and ovine genome mapping information, the aim of the study was to analyze the genetic diversity of selected DNA microsatellites from the bovine genome and to investigate their correlation with the average daily milk yield in Awassi sheep. 18 informative microsatellite markers were selected from the significant QTL regions affecting milk yield identified in the bovine genome in previous studies. The selected microsatellite markers were then amplified by PCR as reciprocal amplifications on the genomic DNA of Awassi sheep, with standard daily milk yield records. Thus, in this study, 18 microsatellite markers associated with milk yield in the bovine genome were examined for both determination of genetic polymorphism within the flock and the effects of marker loci on average daily milk yield in Awassi sheep. Allele frequencies of markers were determined based on the results of fragment analysis. The analysis of variance showed that the 123bp allele at the marker locus BMS1341 on BTA2 significantly influenced the average daily milk yield of Ivesi sheep (P < 0.01). On the other hand, the BMS381 locus with a 115bp allele on BTA2, the MCM140 locus with a 185bp allele on BTA6, the BMS2721 locus with a 155bp allele, the BM1237 locus with 174 and 180bp alleles on BTA7, and finally, the BMS1967 locus with a 117bp allele, the BM4208 locus with 176 and 182bp alleles, and the INRA locus with a185 bp allele on BTA8 showed moderately significant effects on the average daily milk yield of Ivesi ewes (P < 0.05).

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application

BackgroundStructural variants (SVs) such as deletions, duplications, and insertions are known to contribute to phenotypic variation but remain challenging to identify and genotype. A more complete, accessible, and assessable collection of SVs will assist efforts to study SV function in cattle and to incorporate SV genotyping into animal evaluation.ResultsIn this work we produced a large and deeply characterized collection of SVs in Holstein cattle using two popular SV callers (Manta and Smoove) and publicly available Illumina whole-genome sequence (WGS) read sets from 310 samples (290 male, 20 female, mean 20X coverage). Manta and Smoove identified 31 K and 68 K SVs, respectively. In total the SVs cover 5% (Manta) and 6% (Smoove) of the reference genome, in contrast to the 1% impacted by SNPs and indels. SV genotypes from each caller were confirmed to accurately recapitulate animal relationships estimated using WGS SNP genotypes from the same dataset, with Manta genotypes outperforming Smoove, and deletions outperforming duplications. To support efforts to link the SVs to phenotypic variation, overlapping and tag SNPs were identified for each SV, using genotype sets extracted from the WGS results corresponding to two bovine SNP chips (BovineSNP50 and BovineHD). 9% (Manta) and 11% (Smoove) of the SVs were found to have overlapping BovineHD panel SNPs, while 21% (Manta) and 9% (Smoove) have BovineHD panel tag SNPs. A custom interactive database (https://svdb-dc.pslab.ca) containing the identified sequence variants with extensive annotations, gene feature information, and BAM file content for all SVs was created to enable the evaluation and prioritization of SVs for further study. Illustrative examples involving the genes POPDC3, ORM1, G2E3, FANCI, TFB1M, FOXC2, N4BP2, GSTA3, and COPA show how this resource can be used to find well-supported genic SVs, determine SV breakpoints, design genotyping approaches, and identify processed pseudogenes masquerading as deletions.ConclusionsThe resources developed through this study can be used to explore sequence variation in Holstein cattle and to develop strategies for studying SVs of interest. The lack of overlapping and tag SNPs from commonly used SNP chips for most of the SVs suggests that other genotyping approaches will be needed (for example direct genotyping) to understand their potential contributions to phenotype. The included SV genotype assessments point to challenges in characterizing SVs, especially duplications, using short-read data and support ongoing efforts to better characterize cattle genomes through long-read sequencing. Lastly, the identification of previously known functional SVs and additional CDS-overlapping SVs supports the phenotypic relevance of this dataset.

Telomere-to-telomere assemblies of cattle and sheep Y-chromosomes uncover divergent structure and gene content.

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. Here, we assemble complete and gapless telomere to telomere (T2T) Y chromosomes for these species. We find that the pseudo-autosomal regions are similar in length, but the total chromosome size is substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity is accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 19MYA. The centromeres also differ dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosomes have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

DeepGenomeScan of 15 Worldwide Bovine Populations Detects Spatially Varying Positive Selection Signals.

Identifying genomic regions under selection is essential for understanding the genetic mechanisms driving species evolution and adaptation. Traditional methods often fall short in detecting complex, spatially varying selection signals. Recent advances in deep learning, however, present promising new approaches for uncovering subtle selection signals that traditional methods might miss. In this study, we utilized the deep learning framework DeepGenomeScan to detect spatially varying selection signatures across 15 bovine populations worldwide. Our analysis uncovered novel insights into selective sweep hotspots within the bovine genome, revealing key genes associated with physiological and adaptive traits that were previously undetected. We identified significant quantitative trait loci linked to milk protein and fat percentages. By comparing the selection signatures identified in this study with those reported in the Bovine Genome Variation Database, we discovered 38 novel genes under selection that were not identified through traditional methods. These genes are primarily associated with milk and meat yield and quality. Our findings enhance our understanding of spatially varying selection's impact on bovine genomic diversity, laying a foundation for future research in genetic improvement and conservation. This is the first deep learning-based study of selection signatures in cattle, offering new insights for evolutionary and livestock genomics research.

Whole genome re-sequencing reveals high altitude adaptation signatures and admixture in Ladakhi cattle

Ladakhi cattle, native to the high-altitude region of Ladakh in northern India (ranging from 3,000 to 5,000 m above sea level), have evolved unique genetic adaptations to thrive in harsh environmental conditions, such as hypoxia, extreme cold, and low humidity. This study explored the genome of Ladakhi cattle to investigate genetic structure, selection signatures, and adaptive mechanisms. Whole genome sequencing reads, generated on Illumina NovaSeq 6000 platform, were aligned to the Bos taurus reference genome with BWA-MEM. SNPs were identified and filtered using GATK and bcftools, and functionally annotated with SnpEff. For population genomic analysis, PCA and admixture modeling assessed genetic structure, while Neighbor-Joining trees, LD decay, nucleotide diversity (π), and FST evaluated phylogenetic relationships and genetic variation. Selective sweeps were detected using RAiSD, and gene-set enrichment and protein–protein interaction analyses were conducted to explore functional pathways related to adaptation. The study revealed 3,759,279 unique SNPs and demonstrated that Ladakhi cattle form a distinct genetic cluster with an estimated admixture of 68 % Bos indicus and 32 % Bos taurus ancestry. Key findings include rapid linkage disequilibrium decay, low inbreeding level, and the identification of selection signatures and genes associated with hypoxia response, energy metabolism, and cold adaptation. Mean nucleotide diversity (π, 0.0037) and FST values indicated moderate genetic differentiation from other breeds. The analysis highlighted selection signatures for genes like HIF1A, ENO4, ANGPT1, EPO, NOS3, MAPK3, HMOX1, BCL2,CAMK2D, MTOR, AKT2,PIK3CB, and MAP2K1, among others, including various keratin and heat shock proteins. The interaction between genes associated with hypoxia signaling (HIF-1) and other enriched pathways such as PI3K, mTOR, NFκB, ERK, and ER stress, reveals a complex mechanism for managing hypoxic stress in Ladakhi cattle. These findings offer valuable insights for breeding programs aimed at enhancing livestock resilience in extreme environments and enhance understanding of mammalian adaptation to high-altitude conditions.

Development and validation of a 5K low-density SNP chip for Hainan cattle

BackgroundThis study aimed to design and develop a 5K low-density liquid chip for Hainan cattle utilizing targeted capture sequencing technology. The chip incorporates a substantial number of functional single nucleotide polymorphism (SNP) loci derived from public literature, including SNP loci significantly associated with immunity, heat stress, meat quality, reproduction, and other traits. Additionally, SNPs located in the coding regions of immune-related genes from the Bovine Genome Variation Database (BGVD) and Hainan cattle-specific SNP loci were included.ResultsA total of 5,293 SNPs were selected, resulting in 9,837 DNA probes with a coverage rate of 85.69%, thereby creating a Hainan cattle-specific 5K Genotyping by Target Sequencing (GBTS) liquid chip. Evaluation with 152 cattle samples demonstrated excellent clustering performance and a detection rate ranging from 96.60 to 99.07%, with 94.5% of SNP sites exhibiting polymorphism. The chip achieved 100% gender coverage and displayed a heterozygosity rate between 14.20% and 29.65%, with a repeatability rate of 99.65–99.85%. Analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed the potential regulatory roles of exonic SNPs in immune response pathways.ConclusionThe development and validation of the 5K GBTS liquid chip for Hainan cattle represent a valuable tool for genome analysis and genetic diversity assessment. Furthermore, it facilitates breed identification, gender determination, and kinship analysis, providing a foundation for the efficient utilization and development of local cattle genetic resources.

PSVII-10 Genome-wide discovery and characterization of short interspersed nuclear elements (SINEs) in the bovine genome

Abstract Short interspersed nuclear elements (SINEs) constitute a class of nonautonomous retrotransposons endowed with the capability to multiply and spread within a genome. SINEs serve as reliable indicators for assessing genetic variability in dairy populations, aiding in precise determination of breed composition. However, SINEs within the bovine genome remain insufficiently characterized, hindering the understanding of their impact on dairy production and health. This study aimed to systematically identify, classify and unveil the characteristics of SINEs within the bovine reference genome (ARS-UCD1.3). Multiple de novo identification tools, including RepeatModeler, SINE_scan, LongRepMarker, and SINE-finder, were employed for SINE identification and characterization. A total of 1.6 million copies of SINEs were identified, constituting approximately 10.8% of the genome. The SINEs were stratified into seven distinct families based on their structural attributes, denoted as BosSINEL, BosSINE1, BosSINE2, BosSINE3, BosSINE4, BosSINE5 and BosSINE6. BosSINE1 represents a novel family, whereas the remaining six families were previously documented in Repbase. Among these families, BosSINE3 emerges as the most prolific, boasting 254,586 complete copies (over 80% of the length of the consensus sequence), followed by BosSINEL (164,538 copies), BosSINE2 (162,207 copies) and BosSINE1 (150,799 copies). Conversely, the remaining families exhibit lower abundance, each harboring no more than 100,000 copies within the genome. Structural comparisons among SINE families elucidate distinctive features. Notably, BosSINEL lacks a tRNA-related region, including boxA and boxB, which is the hallmark element of the classic SINE structure according to SINEBase. Conversely, all other families exhibit tRNA-related regions. BosSINE1, the lengthiest at 328 bp (excluding Poly-A), possesses two tRNA-related regions, distinguishing it from other families. BosSINEL, BosSINE1, BosSINE2 and BosSINE3 exhibit long interspersed nuclear elements (LINE)-related regions, suggesting their retrotransposition by engaging with the enzymatic machinery encoded by LINEs. In contrast, BosSINE6, the shortest at 109 bp, features a GC-rich region resulting in a high GC content of 67%. Additionally, BosSINE4 contains an AT-rich region. Evolutionary analysis disclosed that BosSINEL is the youngest SINE family in the bovine genome, potentially still harboring active roles and contributing to the genetic diversity related to different dairy production and health performance. Conversely, BosSINE1 to BosSINE6 are comparatively aged in terms of evolution, likely fixed in the bovine genome with diminished or lost activities. In summary, this investigation unveils the comprehensive profile of SINE in the bovine genome, delineating both structural and evolutionary characteristics. Particularly noteworthy is the emergence of BosSINEL as the youngest SINE family, suggesting active roles in the bovine genome. Since lineage-specific TEs (such as Bov-A2 elements) have been reported to regulate interferon-inducible gene expression in bovine, our data warrants further investigations to reveal the potential of SINEs as markers for improving dairy production and health.

PSVI-24 Gene expression profiles of granulosa and peripheral white blood cells from fertile and sub-fertile beef heifers

Abstract Beef heifer fertility issues contribute to a major economic loss in the cow-calf production industry. Therefore, identifying beef heifers with superior genetic potential for improved fertility would increase profitability. This study aimed to identify differences in the transcriptome profiles from granulosa and peripheral white blood cells (pWBCs) of beef heifers with varying reproductive potential. For this, Angus-Simmental crossbred heifers were subjected to an estrus synchronization and fixed-time artificial insemination (AI) protocol (7-D CO-Synch + CIDR) followed by exposure for 60-d to a fertile bull. Depending on the presence or absence of conceptus 120 d post-AI, heifers were classified as fertile (pregnant by AI) or sub-fertile (non-pregnant by AI or bull-breeding). Pregnancies were terminated, and all animals in both groups (fertile, n = 8; and sub-fertile, n = 5) were cycling when the blood and ovaries were collected from each heifer. Total RNA was extracted from the pWBCs and granulosa cells and subjected to library preparation and sequencing on the Nova-Seq platform. The read counts were obtained after data quality control using FastQC v0.11.9 and MultiQC v1.12 and alignment to the Ensemble’s ARS UCD1.2 Bos taurus genome reference using STAR aligner v2.7.5. The filtered data were subjected to differential expression analysis using DESeq2. We identified 1,061 and 72 significantly differentially expressed genes (DEGs) with P-values ≤ 0.05 and absolute (log2 fold change) ≥ 0.5 from pWBC and granulosa, respectively. Notably, 12 targets including 9 protein coding genes (PLCL1, DNER, GNAS, CDH3, PER1, ITGA2B, CXCL12, ENSBTAG00000048613, ENSBTAG00000051519), bta-mir-2887-1, 5-8S-rRNA and U5 were found as DEGs in both the tissues. Based on a differential co-expression analysis using PCIT, we identified GNAS and DNER as hub genes in pWBC and the granulosa cells of the sub-fertile heifer group. The 12 shared genes were over-represented for pathways such as NF kappa B and chemokine signaling, regulation of actin cytoskeleton, and platelet activation. Some of the identified genes have been previously associated with fertility, while others are novel. A detailed understanding of the underlying biological mechanisms of the top genes and a follow-up study with a larger sample size at different time points could validate the candidates identified in this study for their role as potential therapeutic targets.

403 Impact of lysine, methionine, and histidine deficiencies on gene expression and upstream transcriptional regulators in bovine mammary epithelial cells: An RNA-sequencing analysis

Abstract Essential (E) amino acids (AA) may alter bovine mammary epithelial cell functions by impacting mRNA expression of genes through their upstream transcription factors. The objective of this research was to study the impact of EAA deficiencies on gene expression and to identify upstream transcriptional factors mediating their effects in bovine mammary epithelial cells. Epithelial cells were collected from the mammary glands of lactating dairy cows and cultured in combined media of DMEM/F12 and Medium 170, differentiated for 5 d, and subsequently treated in triplicate for 60 h with a control medium containing all EAA at normal plasma concentrations for cows (CTL), or a medium with 20 μM Lys (LY), 5 μM Met (LM), or 10 μM His (LH), which are each 25% of normal concentration. RNA was isolated from the cells and sequenced on an Illumina NovaSeq 6000 producing 150 base pair paired-end reads with a minimum sequencing depth of 150 million reads. Reads were cleaned with ‘fastp’ (Chen et al., 2018), aligned to the bovine reference genome (ARS-UCD1.2) using ‘STAR’ aligner (Dobin et al., 2013), and genes counted using ‘featureCounts’ (Liao et al., 2014). Differential expression of genes between treatment pairs was determined using the ‘DESeq2’ package (Love et al., 2014) in R. Genes with an absolute log2 fold change ≥ 1 and a false discovery rate of &amp;lt; 0.05 were considered differentially expressed. The number of differentially expressed genes relative to CTL was 780 for LY, 383 for LH, and 536 for LM. QIAGEN Ingenuity Pathway Analysis (IPA) identified four pathways consistently affected across LY, LM, and LH: kinetochore metaphase signaling, mitotic roles of polo-like kinase, cell cycle: G2/M DNA damage checkpoint regulation, and interferon signaling (Figure 1-3). Notably, the kinetochore metaphase signaling pathway was identified as the most significant across all treatments, suggesting a strong link to cell proliferation and apoptosis during EAA deficiency. The interferon signaling might be related to the Janus kinase 2/signal transducer and activator of transcription 5 (Jak2/STAT5) pathway which can promote the expression of genes encoding milk proteins such as β-casein and β-lactoglobulin. Furthermore, the IPA identified upstream transcription factors related to cell proliferation, apoptosis, autophagy, protein translation, endoplasmic reticulum biogenesis, and amino acid metabolism including NUPR1, FOXM1, MYC, CCND1, TP53, DDIT3, and ATF4 across all treatments. The findings suggest that AA deprivation regulates bovine mammary epithelial cell function by instigating amino acid response pathways and altering secretory cell number.

204 Changes in DNA methylation 5 days after exposure to acute heat stress in beef cattle skeletal muscle

Abstract The objective of this study was to investigate the effects of acute heat stress on the skeletal muscle epigenome using reduced representation bisulfite sequencing (RRBS). March-born crossbred (5/8 Red Angus, 3/8 Simmental) steers (n = 14; average body weight = 308.7 kg) were exposed to heat stress (HS) conditions of 35°C + 35% relative humidity (Temperature Humidity Index = 81) for 12h/d for 2 d (48 h). Longissimus dorsi samples were collected via biopsy and flash-frozen in liquid nitrogen immediately following the HS period at 48 h and from the contralateral longissimus dorsi 5 d later. DNA was isolated and underwent RRBS using 150 bp paired-end reads at an achieved average depth of 15.3 million reads/sample. Reads were trimmed and mapped to a bisulfite-converted bovine genome (ARSUCD1.2). Methylation calls were extracted and used to identify differentially methylated regions (DMRs) with methylkit in R comparing samples collected 5 d post-HS to those taken immediately post-48 h of HS. Differentially methylated regions (1,000 bp) had q &amp;lt; 0.05 and at least a 15% difference in methylation. There were 1,015 DMRs in genes. Hypomethylated DMRs 5 d post-HS included IGF1R, IGF2, IGF2BP1, IGFBP4, IGF2R, IKBKB, and IL34. Hypermethylated DMRs 5 d post-HS included FGF5, FGFR1, and FGFR3. There were 81 DMRs in putative promoter regions (2000 bp 5’ of a gene). MYLK2, SIGIRR, and 25 other genes had hypomethylated promoter regions. Genes that contained DMRs were used as input for DAVID to explore KEGG pathways. Enriched pathways (P &amp;lt; 0.1) included multiple hormone pathways, MAPK, PI3K-Akt, HIF-1, and VEGF. Hypomethylation of IGF-related genes and hypermethylation of FGF-related genes suggests impaired growth capacity post-HS, as does the enrichment of MAPK and PI3K-Akt pathways. Hypomethylation of SIGIRR’s promoter region, IL34 and IKBKB, indicates enhanced anti-inflammatory activity given their known regulation of inflammatory cytokine production. This, coupled with enrichment of the HIF-1 and VEGF pathways, suggests regulatory dampening of HS-induced oxidative stress. These findings reveal key regulatory changes in muscle growth, inflammatory response, and oxidative stress activity for cattle recovering from acute HS. Previous cattle skeletal muscle RNA studies predicted that HS altered metabolism, inflammation, and oxidative stress pathways, similar to the methylation results of this study. Overall, evidence suggests growth dysregulation, which may help explain observed HS-induced poor growth due to HS, as well as subsequent anti-inflammatory activity during recovery from HS.

Temporal stability of the rumen microbiome and its longitudinal associations with performance traits in beef cattle

The rumen microbiome is the focus of a growing body of research, mostly based on investigation of rumen fluid samples collected once from each animal. Exploring the temporal stability of rumen microbiome profiles is imperative, as it enables evaluating the reliability of findings obtained through single-timepoint sampling. We explored the temporal stability of rumen microbiomes considering taxonomic and functional aspects across the 7-month growing-finishing phase spanning 6 timepoints. We identified a temporally stable core microbiome, encompassing 515 microbial genera (e.g., Methanobacterium) and 417 microbial KEGG genes (e.g., K00856—adenosine kinase). The temporally stable core microbiome profiles collected from all timepoints were strongly associated with production traits with substantial economic and environmental impact (e.g., average daily gain, daily feed intake, and methane emissions); 515 microbial genera explained 45–83%, and 417 microbial genes explained 44–83% of their phenotypic variation. Microbiome profiles influenced by the bovine genome explained 54–87% of the genetic variation of bovine traits. Overall, our results provide evidence that the temporally stable core microbiome identified can accurately predict host performance traits at phenotypic and genetic level based on a single timepoint sample taken as early as 7 months prior to slaughter.

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing.

Structural variants play an important role in evolutionary processes. Besides, they constitute a large source of inter individual genetic variation that might represent a major factor in the aetiology of complex, multifactorial traits. Their importance in adaptation is becoming increasingly evident in literature. Yet, the characterization of the genomic landscape of structural variants in local breeds remains scarce to date. Herein, we investigate patterns and gene annotation of structural variants in the Creole cattle from Guadeloupe breed using whole genome sequences from 23 bulls representative of the population. In total, we detected 32821 ascertained SV defining 15258 regions, representing ~ 17% of the Creole cattle genome. Among these, 6639 regions have not been previously reported in the Database of Genomic Variants archive. Average number of structural variants detected per individual in the studied population is in the same order of magnitude of that observed in indicine populations and higher than that reported in taurine breeds. We observe an important within-individual variability where approximately half of the detected structural variants have low frequency (MAF < 0.25). Most of the detected structural variants (55%) occurred in intergenic regions. Genic structural variants overlapped with 7793 genes and the predicted effect of most of them is ranked as "modifier". Among the structural variants that were predicted to have a high functional impact on the protein, a 5.5 Kb in length, highly frequent deletion on chromosome 2, affects ALPI, a gene associated with the interaction between gut microbiota and host immune system. The 6639 newly identified structural variants regions include three deletions and three duplications shared by more than 80% of individuals that are significantly enriched for genes related to tRNA threonylcarbamoyladenosine metabolic process, important for temperature adaptation in thermophilic organisms, therefore suggesting a potential role in the thermotolerance of Creole cattle from Guadeloupe cattle to tropical climate. Overall, highly frequent structural variants that are specific to the Creole cattle population encompass olfactory receptor and immunity genes as well as genes involved in muscle tone, muscle development and contraction. Beyond mapping and characterizing structural variants in the Creole cattle from Guadeloupe breed, this study provides valuable information for a better understanding of the potential role of chromosomal rearrangements in adaptive traits in cattle.

Analysis of Runs of Homozygosity in Aberdeen Angus Cattle.

A large number of cattle breeds have marked phenotypic differences. They are valuable models for studying genome evolution. ROH analysis can facilitate the discovery of genomic regions that may explain phenotypic differences between breeds affecting traits of economic importance. This paper investigates genome-wide ROH of 189 Aberdeen Angus bulls using the Illumina Bovine GGP HD Beadchip150K to structurally and functionally annotate genes located within or in close ROH of the Aberdeen Angus cattle genome. The method of sequential SNP detection was used to determine the ROH. Based on this parameter, two ROH classes were allocated. The total length of all ROH islands was 11,493 Mb. As a result of studying the genomic architecture of the experimental population of Aberdeen Angus bulls, nine ROH islands and 255 SNPs were identified. Thirteen of these overlapped with regions bearing 'selection imprints' previously identified in other breeds of cattle, and five of these regions were identified in other Aberdeen Angus populations. The total length of the ROH islands was 11,493 Mb. The size of individual islands ranged from 0.038 to 1.812 Mb. Structural annotation showed the presence of 87 genes within the identified ROH islets.

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle.

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China, which has the advantages of fast growth, excellent meat quality, improved tolerance ability, and important landscape value. Copy number variation (CNV) is a significant source of gene structural variation and plays a crucial role in evolution and phenotypic diversity. Based on the latest reference genome ARS-UCD2.0, this study analyzed the genome-wide distribution of CNVs in Yunling cattle using short-read whole-genome sequencing data (n = 129) and single-molecule long-read sequencing data (n = 1), and a total of 16,507 CNVs were detected. After merging CNVs with overlapping genomic positions, 3,728 CNV regions (CNVRs) were obtained, accounting for 0.61% of the reference genome. The functional analysis indicated significant enrichment of CNVRs in 96 GO terms and 57 KEGG pathways, primarily related to cell adhesion, signal transduction, neuromodulation, and nutritional metabolism. Additionally, 111 CNVRs overlapped with 76 quantitative trait loci (QTLs), including Subcutaneous fat thickness QTL, Longissimus muscle area QTL, and Marbling score QTL. Several CNVR-overlapping genes, including BZW1, AOX1, and LOC100138449, overlap with regions associated with meat color and quality QTLs. Furthermore, Vst analysis showed that PSMB4, ERICH1, SMC2, and PPP4R3A were highly divergent between Yunling and Brahman cattle. In summary, we have constructed the genomic CNV map of Yunling cattle for the first time using whole-genome resequencing. This provides valuable genetic variation resources for the study of the Yunling cattle genome and contributes to the study of economic traits in Yunling cattle.

Beef Cattle Genome Project: Advances in Genome Sequencing, Assembly, and Functional Genes Discovery.

Beef is a major global source of protein, playing an essential role in the human diet. The worldwide production and consumption of beef continue to rise, reflecting a significant trend. However, despite the critical importance of beef cattle resources in agriculture, the diversity of cattle breeds faces severe challenges, with many breeds at risk of extinction. The initiation of the Beef Cattle Genome Project is crucial. By constructing a high-precision functional annotation map of their genome, it becomes possible to analyze the genetic mechanisms underlying important traits in beef cattle, laying a solid foundation for breeding more efficient and productive cattle breeds. This review details advances in genome sequencing and assembly technologies, iterative upgrades of the beef cattle reference genome, and its application in pan-genome research. Additionally, it summarizes relevant studies on the discovery of functional genes associated with key traits in beef cattle, such as growth, meat quality, reproduction, polled traits, disease resistance, and environmental adaptability. Finally, the review explores the potential of telomere-to-telomere (T2T) genome assembly, structural variations (SVs), and multi-omics techniques in future beef cattle genetic breeding. These advancements collectively offer promising avenues for enhancing beef cattle breeding and improving genetic traits.

Gene Editing Cattle for Enhancing Heat Tolerance: A Welfare Review of the “PRLR-SLICK Cattle” Case

In March 2022 the US Food and Drug Administration (FDA) published a risk assessment of a recent animal gene editing proposal submitted by Acceligen™. The proposal concerned the possibility of changing the cattle genome to obtain a slicker, shorter hair coat. Using CRISPR-Cas9 it was possible to introduce an intentional genomic alteration (IGA) to the prolactin receptor gene (PRLR), thereby producing PRLR-SLICK cattle. The goal was to diminish heat stress in the cattle by enhancing their heat-tolerance. With regard to unintended alterations (i.e., off-target effects), the FDA stated that the IGA posed a low, but still present, risk to animal safety. The aim of this article is to present some initial insights into the welfare issues raised by PRLR-SLICK cattle by addressing the question: Do SLICK cattle have better welfare than non-SLICK cattle when exposed to heat stress? Two potential welfare concerns are examined. The first is pleiotropy, an issue that arises when one gene affects multiple traits. Given the pleiotropic nature of prolactin, it has been suggested that the IGA for SLICK cattle may also affect their hepatic and other functions. The second concern relates not primarily to direct effects on cattle health, but rather to the indirect risk that this more heat-tolerant animal would just be used in the livestock sector under farming conditions that are such that the net welfare improvement would be non-existent.

Single-cell transcriptomics across 2,534 microbial species reveals functional heterogeneity in the rumen microbiome.

Deciphering the activity of individual microbes within complex communities and environments remains a challenge. Here we describe the development of microbiome single-cell transcriptomics using droplet-based single-cell RNA sequencing and pangenome-based computational analysis to characterize the functional heterogeneity of the rumen microbiome. We generated a microbial genome database (the Bovine Gastro Microbial Genome Map) as a functional reference map for the construction of a single-cell transcriptomic atlas of the rumen microbiome. The atlas includes 174,531 microbial cells and 2,534 species, of which 172 are core active species grouped into 12 functional clusters. We detected single-cell-level functional roles, including a key role for Basfia succiniciproducens in the carbohydrate metabolic niche of the rumen microbiome. Furthermore, we explored functional heterogeneity and reveal metabolic niche trajectories driven by biofilm formation pathway genes within B. succiniciproducens. Our results provide a resource for studying the rumen microbiome and illustrate the diverse functions of individual microbial cells that drive their ecological niche stability or adaptation within the ecosystem.

Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle.

Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.

Nuclear sequences of mitochondrial origin in domestic yak

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation.

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.