Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle.

BACKGROUND: Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19). RESULTS: A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima's D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle. CONCLUSIONS: We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle.

BACKGROUND: The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data. RESULT: The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1). CONCLUSION: We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Using expression data to fine map QTL associated with fertility in dairy cattle.

BACKGROUND: Female fertility is an important trait in dairy cattle. Identifying putative causal variants associated with fertility may help to improve the accuracy of genomic prediction of fertility. Combining expression data (eQTL) of genes, exons, gene splicing and allele specific expression is a promising approach to fine map QTL to get closer to the causal mutations. Another approach is to identify genomic differences between cows selected for high and low fertility and a selection experiment in New Zealand has created exactly this resource. Our objective was to combine multiple types of expression data, fertility traits and allele frequency in high- (POS) and low-fertility (NEG) cows with a genome-wide association study (GWAS) on calving interval in Australian cows to fine-map QTL associated with fertility in both Australia and New Zealand dairy cattle populations. RESULTS: Variants that were significantly associated with calving interval (CI) were strongly enriched for variants associated with gene, exon, gene splicing and allele-specific expression, indicating that there is substantial overlap between QTL associated with CI and eQTL. We identified 671 genes with significant differential expression between POS and NEG cows, with the largest fold change detected for the CCDC196 gene on chromosome 10. Our results provide numerous candidate genes associated with female fertility in dairy cattle, including GYS2 and TIGAR on chromosome 5 and SYT3 and HSD17B14 on chromosome 18. Multiple QTL regions were located in regions with large numbers of copy number variants (CNV). To identify the causal mutations for these variants, long read sequencing may be useful. CONCLUSIONS: Variants that were significantly associated with CI were highly enriched for eQTL. We detected 671 genes that were differentially expressed between POS and NEG cows. Several QTL detected for CI overlapped with eQTL, providing candidate genes for fertility in dairy cattle.

Definition of metafounders based on population structure analysis.

BACKGROUND: Limitations of the concept of identity by descent in the presence of stratification within a breeding population may lead to an incomplete formulation of the conventional numerator relationship matrix ( A ). Combining A with the genomic relationship matrix ( G ) in a single-step approach for genetic evaluation may cause inconsistencies that can be a source of bias in the resulting predictions. The objective of this study was to identify stratification using genomic data and to transfer this information to matrix A , to improve the compatibility of A and G . METHODS: Using software to detect population stratification (ADMIXTURE), we developed an iterative approach. First, we identified 2 to 40 strata ( k ) with ADMIXTURE, which we then introduced in a stepwise manner into matrix A , to generate matrix A Γ using the metafounder methodology. Improvements in consistency between matrix G and A Γ were evaluated by regression analysis and through the comparison of the overall mean and mean diagonal values of both matrices. The approach was tested on genotype and pedigree information of European and North American Brown Swiss animals (85,249). Analyses with ADMIXTURE were initially performed on the full set of genotypes (S1). In addition, we used an alternative dataset where we avoided sampling of closely related animals (S2). RESULTS: Results of the regression analyses of standard A on G were - 0.489, 0.780 and 0.647 for intercept, slope and fit of the regression. When analysing S1 data results of the regression for A Γ on G corresponding values were - 0.028, 1.087 and 0.807 for k =7, while there was no clear optimum k . Analyses of S2 gave a clear optimal k =24, with - 0.020, 0.998 and 0.817 as results of the regression. For this k differences in mean and mean diagonal values between both matrices were negligible. CONCLUSIONS: The derivation of hidden stratification information based on genotyped animals and its integration into A improved compatibility of the resulting A Γ and G considerably compared to the initial situation. In dairy breeding populations with large half-sib families as sub-structures it is necessary to balance the data when applying population structure analysis to obtain meaningful results.

Whole genome sequences of 70 indigenous Ethiopian cattle.

Indigenous animal genetic resources play a crucial role in preserving global genetic diversity and supporting the livelihoods of millions of people. In Ethiopia, the majority of the cattle population consists of indigenous breeds. Understanding the genetic architecture of these cattle breeds is essential for effective management and conservation efforts. In this study, we sequenced DNA samples from 70 animals from seven indigenous cattle breeds, generating about two terabytes of pair-end reads with an average coverage of 14X. The sequencing data were pre-processed and mapped to the cattle reference genome (ARS-UCD1.2) with an alignment rate of 99.2%. Finally, the variant calling process produced approximately 35 million high-quality SNPs. These data provide a deeper understanding of the genetic landscape, facilitate the identification of causal mutations, and enable the exploration of evolutionary patterns to assist cattle improvement and sustainable utilization, particularly in the face of unpredictable climate changes.

Conservation and upgradation of indigenous Cholistani cattle breed of Pakistan: a pathway to sustainable livestock.

Pakistan is endowed with many established indigenous zebu Bos indicus type (humped) cattle breeds including Sahiwal, Red Sindhi, Bhagnari and Cholistani. Amongst these indigenous cattle breeds, Sahiwal and Red Sindhi have extensively been navigated and hence these two are acclaimed as internationally recognized breeds. However, research work on Cholistani cattle breed actually initiated in 2010 and has attained a steady pace. This breed was a new entrant in Livestock Census of Pakistan since 2006. Cholistani is a hardy, tick-resistant, adaptable cattle breed being reared under pastoral nomadism of the Cholistan desert, Pakistan. The present narrative review is the first of its kind intended to sum-up all the research work conducted about this indigenous cattle breed, and to put forth research gaps for this formerly neglected cattle breed. The review discusses the research work conducted on Cholistani cattle breed under five major research subjects/domains i.e. production attributes, theriogenology-related attributes, hematochemical attributes, disease, epidemiologic and therapeutic attributes, and genetic attributes. Future horizon for research avenues has also been given. It is the dire need of time that specific breed-oriented conservation and propagation programs may be initiated in the country so that sustained livestock and enhance socioeconomic profiling of rural communities may be attained.

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip.

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1-5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Equivalence of variance components between standard and recursive genetic models using LDL' transformations.

BACKGROUND: Recursive models are a category of structural equation models that propose a causal relationship between traits. These models are more parameterized than multiple trait models, and they require imposing restrictions on the parameter space to ensure statistical identification. Nevertheless, in certain situations, the likelihood of recursive models and multiple trait models are equivalent. Consequently, the estimates of variance components derived from the multiple trait mixed model can be converted into estimates under several recursive models through LDL' or block-LDL' transformations. RESULTS: The procedure was employed on a dataset comprising five traits (birth weight-BW, weight at 90 days-W90, weight at 210 days-W210, cold carcass weight-CCW and conformation-CON) from the Pirenaica beef cattle breed. These phenotypic records were unequally distributed among 149,029 individuals and had a high percentage of missing data. The pedigree used consisted of 343,753 individuals. A Bayesian approach involving a multiple-trait mixed model was applied using a Gibbs sampler. The variance components obtained at each iteration of the Gibbs sampler were subsequently used to estimate the variance components within three distinct recursive models. CONCLUSIONS: The LDL' or block-LDL' transformations applied to the variance component estimates achieved from a multiple trait mixed model enabled inference across multiple sets of recursive models, with the sole prerequisite of being likelihood equivalent. Furthermore, the aforementioned transformations simplify the handling of missing data when conducting inference within the realm of recursive models.

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics.

BACKGROUND: Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs. RESULTS: Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample. CONCLUSIONS: The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

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.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle.

BACKGROUND: Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle. RESULTS: We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits. CONCLUSIONS: This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Lifetime milk production of Holstein cattle in the humid tropics compared to Holstein-Gyr and Holstein-Brahman crosses.

Crossbred cattle are commonly used for milk production in the tropics, combining the potential benefits of pure breeds with the heterosis effects of the offspring. However, no comprehensive assessment of lifetime productivity for crossbred versus purebred cattle in low-altitude tropical environments has been carried out. The present study compares the lifetime productivity of purebred Holstein (HO, n = 17,269), Gyr (GY4, n = 435), and Brahman (BR4, n = 622) with crossbreds Gyr × Holstein (GY × HO, n = 5521) and Brahman×Holstein (BR × HO, n = 5429) cows from dairy farms located in low and medium altitude tropical regions in Costa Rica. The production traits of interest were age at first calving (AFC), days open (DO), milk production per lactation (TMP), lactation length (LLEN), age at culling (ACUL), and number of lactations (NLAC). Estimates of heterosis were also calculated. The AFC for GY × HO crosses (33-34 months) was not significantly different (p > .05) from HO (33.8 months). For BR × HO crosses, a significant (p < .05) decrease in AFC (BR3HO1 35.6 months, BR2HO2 34.5 months, and BR1H03 33.3 months) was observed as the fraction of HO breed increased. Estimates of heterosis for AFC were favourable for both crosses, of a magnitude close to 3%. The DO for F1 crosses (GY2HO2 94 days; BR2HO2 96 days) was significantly (p < .05) lower than HO (123 days). Estimates of heterosis for DO were also favourable and above 15% for both crosses. The TMP and LLEN were higher for HO (TMP = 5003 kg; LLEN = 324 days) compared with GY × HO (TMP = 4428 to 4773 kg; LLEN = 298 to 312 days) and BR × HO (TMP = 3950 to 4761 kg; LLEN = 273 to 313 days) crosses. Heterosis for TMP was favourable but low for both crosses, with a magnitude below 3.0%. The NLAC for HO (4.6 lactations) was significantly (p < .05) lower than F1 (GY2HO2, 5.8 lactations; BR2HO2, 5.4 lactations). Heterosis for NLAC was above 6.0% for both crosses. Overall, estimates of lifetime income over feed costs per cow on average were USD 2637 (30.3%) and USD 734 (8.4%) higher in F1 GY × HO and BR × HO, respectively, compared to HO. In conclusion, crossbred animals, specifically those with Gyr and Brahman genetics, extend the productive lifespan, increasing economic returns.

Nutritional and productive performance of purebred Nellore heifers and crossed with Brangus and Braford finished in a feedlot system.

This study evaluated the nutritional and productive performance of Nellore purebred heifers and crossbred Brangus x Nellore (BGNE) and Braford x Nellore (BFNE) in a feedlot system. Thirty heifers (10 of each genetic group) with an average age of 18 months and an initial body weight of 261 kg were used. The experiment was structured and conducted according to a completely randomized design, with three treatments. Heifers received two diets (60 days each) during the experimental period. The experiment lasted 120 days with four experimental periods. Nellore heifers had a lower intake than crossbred heifers (P <0.05). There were no differences between BGNE and BFNE heifers, which had higher final body weight, average daily gain, feed efficiency, hot carcass weight and carcass length than NE heifers. Crossed heifers presented better fat cover than NE heifers. However, NE heifers had higher carcass dressing Despite presenting lower carcass yields than Nellore heifers, crossed heifers are more efficient and have higher performance and better fat cover on the carcass than purebred Nellore heifers. Crossbreeding synthetic breeds, such as Brangus and Braford breeds, with the Nellore breed is an effective way to increase the productivity and efficiency of feedlot heifers in tropical regions.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Deciphering genetic diversity in conserved cattle bulls to achieve sustainable development goals.

The primary objective of Sustainable Development Goal target 2.5 established by the United Nations is to ensure the preservation of genetic diversity in domesticated animals. The ICAR-National Bureau of Animal Genetic Resources in India has been actively engaged in the conservation of cattle and buffalo bull semen for long-term storage. This present study aimed to assess the genetic diversity present in the conserved cattle bull semen, which would aid in determining the most suitable strategy for future conservation management. A total of 192 bull semen belonging to 19 cattle breeds were selected to evaluate genetic diversity using 17 pairs of FAO recommended microsatellite primers. Total 267 alleles were detected across all the samples which indicates substantial amount of allelic variation is being maintained in conserved bulls. Further, all cattle bulls semen conserved showed higher observed heterozygosity than expected heterozygosity which indicates excess genetic diversity in all the populations. The FST, F IT and FIS value across the loci and population is 0.146 ± 0.009, 0.054 ± 0.038, and - 0.105 ± 0.035, respectively, which suggests lack of inbreeding in conserved cattle bull semen. This study has established genetic diversity in conserved cattle semen samples to achieve sustainable development goals. In addition, it provides compelling evidence that the current approach for conserving cattle bull semen is heading in the correct direction.

The effect of miR-23b-3p on regulating GH by targeting POU1F1 in Yanbian yellow cattle.

This study aimed to evaluate the effect of miR-23b-3p on growth hormone (GH) in pituitary cells of Yanbian yellow cattle. The mRNA and protein levels of GH and miR-23b-3p target genes were measured by real time fluorescence quantitative PCR (qPCR) and Western blot, respectively. The target relationship of miR-23b-3p was validated by double luciferase reporter gene system. The results showed that GH mRNA and protein levels in pituitary cells of Yanbian yellow cattle were significantly lower in the miR-23b-3p-mi group than in the NC group (P<0.01), while GH mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.05). The result of bioinformatics analysis and double luciferase reporter gene system validation proved that miR-23b-3p targeted 3'UTR of pituitary specific transcription factor 1 (POU1F1). POU1F1 mRNA and protein levels were lower miR-23b-3p-mi group than in the NC group (P<0.01), while POU1F1 mRNA and protein levels were higher in the miR-23b-3p-in group than in the iNC group (P<0.01). These results demonstrated that miR-23b-3p could regulate GH expression in pituitary cells by regulating POU1F1 gene.

Spatially resolved transcriptomic profiling of placental development in dairy cow.

The placenta plays a crucial role in successful mammalian reproduction. Ruminant animals possess a semi-invasive placenta characterized by a highly vascularized structure formed by maternal endometrial caruncles and fetal placental cotyledons, essential for full-term fetal development. The cow placenta harbors at least two trophoblast cell populations: uninucleate (UNC) and binucleate (BNC) cells. However, the limited capacity to elucidate the transcriptomic dynamics of the placental natural environment has resulted in a poor understanding of both the molecular and cellular interactions between trophoblast cells and niches, and the molecular mechanisms governing trophoblast differentiation and functionalization. To fill this knowledge gap, we employed Stereo-seq to map spatial gene expression patterns at near single-cell resolution in the cow placenta at 90 and 130 days of gestation, attaining high-resolution, spatially resolved gene expression profiles. Based on clustering and cell marker gene expression analyses, key transcription factors, including YBX1 and NPAS2, were shown to regulate the heterogeneity of trophoblast cell subpopulations. Cell communication and trajectory analysis provided a framework for understanding cell-cell interactions and the differentiation of trophoblasts into BNCs in the placental microenvironment. Differential analysis of cell trajectories identified a set of genes involved in regulation of trophoblast differentiation. Additionally, spatial modules and co-variant genes that help shape specific tissue structures were identified. Together, these findings provide foundational insights into important biological pathways critical to the placental development and function in cows.

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens).

BACKGROUND: There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity. RESULTS: A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks. CONCLUSIONS: Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak.

BACKGROUND: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY). RESULTS: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak. CONCLUSIONS: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Genetic parameters for milk yield and reproductive traits in Honduran Holstein cattle.

The aim of this study was to estimate the (co)variance components and genetic parameters for milk yield adjusted to 305d (MY305), calving-to-conception interval (CCI), number of services per conception (NSC) and calving interval (CI) of Honduran Holstein cows, by fitting a bivariate animal model using Maximum Restricted Likelihood procedures. Model included the fixed effects of calving number, the contemporary calving group (farm-season-year of calving and the cow age as covariate). The estimated means and standard deviations for MY, CCI, NSC and CI were, 5098.60 ± 1564.32 kg, 168.27 ± 104.71 days, 2.46 ± 1.69 services, and 448.73 ± 109.16 days, respectively; and their estimated heritabilities were 0.21 ± 0.05, 0.03 ± 0.028, 0.02 ± 0.024 and 0.06 ± 0.04, respectively. The genetic correlations between MY305 and CCI, NSC and CI were positive and antagonist, with values of 0.64 ± 0.52, 0.99 ± 0.56, and 0.32 ± 0.24 respectively. Even though moderate to low heritability was estimated for MY305, systematic selection for milk yield, with a reduction in reproductive efficiency, if considered as the only selection criterion is important to be considered. By including reproductive traits and considering permanent environment effects into the breeding program, might yield a slow, but constant and permanent improvement over time.