Comparative Transcriptomic Analysis of the Pituitary Gland between Cattle Breeds Differing in Growth: Yunling Cattle and Leiqiong Cattle

Insights into Adaption and Growth Evolution: Genome-Wide Copy Number Variation Analysis in Chinese Hainan Yellow Cattle Using Whole-Genome Re-Sequencing Data

Copy number variation (CNV) serves as a crucial source of genomic variation and significantly aids in the mining of genomic information in cattle. This study aims to analyze re-sequencing data from Chinese Hainan yellow cattle, to uncover breed CNV information, and to elucidate the resources of population genetic variation. We conducted whole-genome sequencing on 30 Chinese Hainan yellow cattle, thus generating 814.50 Gb of raw data. CNVs were called using CNVnator software, and subsequent filtering with Plink and HandyCNV yielded 197,434 high-quality CNVs and 5852 CNV regions (CNVRs). Notably, the proportion of deleted sequences (81.98%) exceeded that of duplicated sequences (18.02%), with the lengths of CNVs predominantly ranging between 20 and 500 Kb This distribution demonstrated a decrease in CNVR count with increasing fragment length. Furthermore, an analysis of the population genetic structure using CNVR databases from Chinese, Indian, and European commercial cattle breeds revealed differences between Chinese Bos indicus and Indian Bos indicus. Significant differences were also observed between Hainan yellow cattle and European commercial breeds. We conducted gene annotation for both Hainan yellow cattle and European commercial cattle, as well as for Chinese Bos indicus and Indian Bos indicus, identifying 206 genes that are expressed in both Chinese and Indian Bos indicus. These findings may provide valuable references for future research on Bos indicus. Additionally, selection signatures analysis based on Hainan yellow cattle and three European commercial cattle breeds identified putative pathways related to heat tolerance, disease resistance, fat metabolism, environmental adaptation, candidate genes associated with reproduction and the development of sperm and oocytes (CABS1, DLD, FSHR, HSD17B2, KDM2A), environmental adaptation (CNGB3, FAM161A, DIAPH3, EYA4, AAK1, ERBB4, ERC2), oxidative stress anti-inflammatory response (COMMD1, OXR1), disease resistance (CNTN5, HRH4, NAALADL2), and meat quality (EHHADH, RHOD, GFPT1, SULT1B1). This study provides a comprehensive exploration of CNVs at the molecular level in Chinese Hainan yellow cattle, offering theoretical support for future breeding and selection programs aimed at enhancing qualities of this breed.

Genome-wide association and functional genomic analyses for body conformation traits in North American Holstein cattle

Body conformation traits are directly associated with longevity, fertility, health, and workability in dairy cows and have been under direct genetic selection for many decades in various countries worldwide. The main objectives of this study were to perform genome-wide association studies and functional enrichment analyses for fourteen body conformation traits using imputed high-density single nucleotide polymorphism (SNP) genotypes. The traits analyzed include body condition score (BCS), body depth (BD), bone quality (BQ), chest width (CW), dairy capacity (DC), foot angle (FAN), front legs view (FLV), heel depth (HDe), height at front end (HFE), locomotion (LOC), rear legs rear view (RLRV), rear legs side view (RLSV), stature (ST), and a composite feet and legs score index (FL) of Holstein cows scored in Canada. De-regressed estimated breeding values from a dataset of 39,135 North American Holstein animals were used as pseudo-phenotypes in the genome-wide association analyses. A mixed linear model was used to estimate the SNP effects, which ranged from 239,533 to 242,747 markers depending on the trait analyzed. Genes and quantitative trait loci (QTL) located up to 100 Kb upstream or downstream of the significant SNPs previously cited in the Animal QTLdb were detected, and functional enrichment analyses were performed for the candidate genes identified for each trait. A total of 20, 60, 13, 17, 27, 8, 7, 19, 4, 10, 13, 15, 7, and 13 genome-wide statistically significant SNPs for Bonferroni correction based on independent chromosomal segments were identified for BCS, BD, BQ, CW, DC, FAN, FLV, HDe, HFE, LOC, RLRV, RLSV, ST, and FL, respectively. The significant SNPs were located across the whole genome, except on chromosomes BTA24, BTA27, and BTA29. Four markers (for BCS, BD, HDe, and RLRV) were statistically significant when considering a much stricter threshold for the Bonferroni correction for multiple tests. Moreover, the genomic regions identified overlap with various QTL previously reported for the trait groups of exterior, health, meat and carcass, milk, production, and reproduction. The functional enrichment analyses revealed 27 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously reported as linked to bone development, metabolism, as well as infectious and immunological diseases.

Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in southwestern China

The Sichuan-Yunnan region is the main production area of yaks in southwestern China, with rich genetic resources of Yaks. Nevertheless, there have been limited study on the genetic characteristics of the entire yak populations in Tibet and southwestern China. In this study, we performed whole-genome resequencing to identify genetic variation information in a total of 198 individuals from six yak breeds (populations) in Sichuan (Muli yak, Jinchuan yak, Changtai yak, Maiwa yak), Yunnan (Zhongdian yak), and Tibet (Tibetan yak). The aim was to investigate the whole-genome genetic diversity, population genetic structure, and genome selection signatures. We observed that all six populations exhibit abundant genetic diversity. Except for Tibetan yaks, which showed low nucleotide diversity (0.00104), the remaining yak populations generally displayed high nucleotide diversity (0.00129-0.00153). Population genetic structure analysis revealed that, among the six yak populations, Muli yak exhibited greater differentiation from other yak populations and formed a distinct cluster independently. The Maiwa yak population displayed a complex genetic structure and exhibited gene exchange with Jinchuan and Changtai yaks. Positive selection signals were detected in candidate genes associated with growth (GNB4, HMGA2, TRPS1, and LTBP1), reproduction (PI4KB, DYNC1I1, and GRIP1), immunity (CD200 and IL1RAP), lactation (SNX13 and CPM), hypoxia adaptation (NDUFB6, PRKN, and MRPS9), hair (KRT24, KRT25, and KRT26), meat quality (SUCLG2), digestion and absorption (CLDN1), and pigment deposition (OCA2) using the integrated Pi and F ST methods. This study provides significant insights into understanding the whole-genome genetic characteristics of yak populations in Tibet and southwestern China.

Chromosome-level genome assembly and annotation of the Yunling cattle with PacBio and Hi-C sequencing data

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China. It is bred by crossing the Brahman, the Murray Grey and the Yunnan Yellow cattle. Yunling cattle can adapt to the tropical and subtropical climate environment, and has good reproductive ability and growth speed under high temperature and high humidity conditions, it also has strong resistance to internal and external parasites and with good beef performance. In this study, we generated a high-quality chromosome-level genome assembly of a male Yunling cattle using a combination of short reads sequencing, PacBio HiFi sequencing and Hi-C scaffolding technologies. The genome assembly(3.09 Gb) is anchored to 31 chromosomes(29 autosomes plus one X and Y), with a contig N50 of 35.97 Mb and a scaffold N50 of 112.01 Mb. It contains 1.62 Gb of repetitive sequences and 20,660 protein-coding genes. This first construction of the Yunling cattle genome provides a valuable genetic resource that will facilitate further study of the genetic diversity of bovine species and accelerate Yunling cattle breeding efforts.

Integrated Analysis of Transcriptome and Metabolome Profiles in the Longissimus Dorsi Muscle of Buffalo and Cattle

Buffalo meat is gaining popularity for its nutritional properties, such as its low fat and cholesterol content. However, it is often unsatisfactory to consumers due to its dark color and low tenderness. There is currently limited research on the regulatory mechanisms of buffalo meat quality. Xinglong buffalo are raised in the tropical Hainan region and are undergoing genetic improvement from draught to meat production. For the first time, we evaluated the meat quality traits of Xinglong buffalo using the longissimus dorsi muscle and compared them to Hainan cattle. Furthermore, we utilized a multi-omics approach combining transcriptomics and metabolomics to explore the underlying molecular mechanism regulating meat quality traits. We found that the Xinglong buffalo had significantly higher meat color redness but lower amino acid content and higher shear force compared to Hainan cattle. Differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were identified, with them being significantly enriched in nicotinic acid and nicotinamide metabolic and glycine, serine, and threonine metabolic pathways. The correlation analysis revealed that those genes and metabolites (such as: GAMT, GCSH, PNP, L-aspartic acid, NADP+, and glutathione) are significantly associated with meat color, tenderness, and amino acid content, indicating their potential as candidate genes and biological indicators associated with meat quality. This study contributes to the breed genetic improvement and enhancement of buffalo meat quality.

Expression profile and bioinformatics analysis of circRNA and its associated ceRNA networks in longissimus dorsi from Lufeng cattle and Leiqiong cattle

This paper aims to explore the role of circRNA expression profiles and circRNA-associated ceRNA networks in the regulation of myogenesis in the longissimus dorsi of cattle breeds surviving under subtropical conditions in southern China by RNA sequencing and bioinformatics analysis. It also aims to provide comprehensive understanding of the differences in muscle fibers in subtropical cattle breeds and to expand the knowledge of the molecular networks that regulate myogenesis. With regard to meat quality indicators, results showed that the longissimus dorsi of LQC had lower pH (P < 0.0001), lower redness (P < 0.01), lower shear force (P < 0.05), and higher brightness (P < 0.05) than the longissimus dorsi of LFC. With regard to muscle fiber characteristics, the longissimus dorsi of LQC had a smaller diameter (P < 0.0001) and higher density of muscle fibers (P < 0.05). The analysis results show that the function of many circRNA-targeted mRNAs was related to myogenesis and metabolic regulation. Furthermore, in the analysis of the function of circRNA source genes, we hypothesized that btacirc_00497 and btacirc_034497 may regulate the function and type of myofibrils by affecting the expression of MYH6, MYH7, and NEB through competitive linear splicing.

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

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

The Association between Gut Microbiome Diversity and Composition and Heat Tolerance in Cattle

Cattle are raised around the world and are frequently exposed to heat stress, whether in tropical countries or in regions with temperate climates. It is universally acknowledged that compared to those in temperate areas, the cattle breeds developed in tropical and subtropical areas have better heat tolerance. However, the underlying mechanism of heat tolerance has not been fully studied, especially from the perspective of intestinal microbiomics. The present study collected fecal samples of cattle from four representative climatic regions of China, namely, the mesotemperate (HLJ), warm temperate (SD), subtropical (HK), and tropical (SS) regions. Then, the feces were analyzed using high-throughput 16S rRNA sequencing. The results showed that with increasing climatic temperature from HLJ to SS, the abundance of Firmicutes increased, accompanied by an increasing Firmicutes to Bacteroidota ratio. Proteobacteria showed a trend of reduction from HLJ to SS. Patescibacteria, Chloroflexi, and Actinobacteriota were particularly highest in SS for adapting to the tropical environment. The microbial phenotype in the tropics was characterized by an increase in Gram-positive bacteria and a decrease in Gram-negative bacteria, aerobic bacteria, and the forming of_biofilms. Consistently, the functional abundances of organismal systems and metabolism were decreased to reduce the material and energy demands in a hot environment. Genetic information processing and information storage and processing may be how gut flora deals with hot conditions. The present study revealed the differences in the structure and function of gut microbes of cattle from mesotemperate to tropical climates and provided an important reference for future research on the mechanism of heat tolerance regulated by the gut microbiota and a potential microbiota-based target to alleviate heat stress.

Genome-Wide Association Study and Selective Sweep Analysis Reveal the Genetic Architecture of Body Weights in a Chicken F2 Resource Population

Rapid growth is one of the most important economic traits in broiler breeding programs. Identifying markers and genes for growth traits may not only benefit marker-assisted selection (MAS)/genomic selection (GS) but also provide important information for understanding the genetic architecture of growth traits in broilers. In the present study, an F₂ resource population derived from a cross between the broiler and Baier yellow chicken (a Chinese local breed) was used and body weights from 1 to 12 weeks of age [body weight (BW) 1–BW12)] were measured. A total of 519 F₂ birds were genome re-sequenced, and a combination of genome-wide association study (GWAS) and selective sweep analysis was carried out to characterize the genetic architecture affecting chicken body weight comprehensively. As a result, 1,539 SNPs with significant effects on body weights at different weeks of age were identified using a genome-wide efficient mixed-model association (GEMMA) package. These SNPs were distributed on chromosomes 1 and 4. Besides, windows under selection identified for BW1–BW12 varied from 1,581 to 2,265. A total of 42 genes were also identified with significant effects on BW1–BW12 based on both GWAS and selective sweep analysis. Among these genes, diacylglycerol kinase eta (DGKH), deleted in lymphocytic leukemia (DLEU7), forkhead box O17 (FOXO1), karyopherin subunit alpha 3 (KPNA3), calcium binding protein 39 like (CAB39L), potassium voltage-gated channel interacting protein 4 (KCNIP4), and slit guidance ligand 2 (SLIT2) were considered as important genes for broiler growth based on their basic functions. The results of this study may supply important information for understanding the genetic architecture of growth traits in broilers.

Expression of candidate genes for residual feed intake in tropically adapted Bos taurus and Bos indicus bulls under thermoneutral and heat stress environmental conditions

The objectives of this study were to measure the relative expression of the ATP1A1, NR3C1, POMC, NPY, and LEP genes in Caracu (Bos taurus) and Nelore (Bos indicus) bulls submitted to feed efficiency tests at high environmental temperatures, and to evaluate differences in adaptability to tropical conditions between breeds. Thirty-five Caracu and 30 Nelore bulls were submitted to a feed efficiency test using automated feeding stations. At the end of the test, the animals were subjected to thermoneutral (TN) and heat stress (HS) conditions. Blood samples were collected after the exposure to the TN and HS conditions and the relative expression of genes was measured by qPCR. The bulls exhibited lower expression of ATP1A1 in the HS condition than in the TN condition (1.98 ± 0.27 and 2.86 ± 0.26, P = 0.02), while the relative expression of NR3C1, POMC, and LEP did not differ (P > 0.05) between climatic conditions. The breed and feed intake influenced NPY and LEP expression levels (P < 0.05). Different climate conditions associated with residual feed intake can modify the gene expression patterns of ATP1A1 and NPY. The association observed among all genes studied shows that they are involved in appetite control. Bos taurus and Bos indicus bulls exhibited similar adaptability to tropical climate conditions.

First report of the prevalence and genetic characterization of Giardia duodenalis and Cryptosporidium spp. in Yunling cattle in Yunnan Province, southwestern China

Yunling cattle is an unique cattle breed distributed in Yunnan Province, southwestern China. It is yet to know whether Yunling cattle are infected with Giardia duodenalis and Cryptosporidium spp.. The objectives of the present study were to investigate the prevalence and characterize the assemblages of G. duodenalis and species of Cryptosporidium spp. in Yunling cattle in Yunnan province. The overall prevalence of G. duodenalis and Cryptosporidium spp. were 10.49% (41/391) and 0.77% (3/391), respectively. The age was considered as the risk factor for Yunling cattle infection with G. duodenalis (χ² = 8.082, OR = 2.56, P = 0.004). Two assemblages of G. duodenalis, assemblage A (n = 1) and assemblage E (n = 40), were identified by amplification of the β-giardin (bg) and glutamate dehydrogenase (gdh) gene loci using the nested PCR methods. Furthermore, Cryptosporidium andersoni (n = 1) and Cryptosporidium ryanae (n = 2) were detected by nested PCR targeting the small subunit (SSU) rRNA gene. This is the first report of G. duodenalis and Cryptosporidium spp. in Yunling cattle in China, which provided baseline date for further studies of the prevalence, genetic identity, and public health potential of these parasites in Yunling cattle.

Control of Protein and Energy Metabolism in the Pituitary Gland in Response to Three-Week Running Training in Adult Male Mice

It is assumed that crosstalk of central and peripheral tissues plays a role in the adaptive response to physical activity and exercise. Here, we wanted to study the effects of training and genetic predisposition in a marathon mouse model on mRNA expression in the pituitary gland. Therefore, we used a mouse model developed by phenotype selection for superior running performance (DUhTP) and non-inbred control mice (DUC). Both mouse lines underwent treadmill training for three weeks or were kept in a sedentary condition. In all groups, total RNA was isolated from the pituitary gland and sequenced. Molecular pathway analysis was performed by ingenuity pathway analysis (IPA). Training induced differential expression of 637 genes (DEGs) in DUC but only 50 DEGs in DUhTP mice. Genetic selection for enhanced running performance strongly affected gene expression in the pituitary gland and identified 1732 DEGs in sedentary DUC versus DUhTP mice. Training appeared to have an even stronger effect on gene expression in both lines and comparatively revealed 3828 DEGs in the pituitary gland. From the list of DEGs in all experimental groups, candidate genes were extracted by comparison with published genomic regions with significant effects on training responses in mice. Bioinformatic modeling revealed induction and coordinated expression of the pathways for ribosome synthesis and oxidative phosphorylation in DUC mice. By contrast, DUhTP mice were resistant to the positive effects of three-week training on protein and energy metabolism in the pituitary gland.