Adipogenesis is a complex and precisely orchestrated process mediated by a network of adipogenic regulatory factors. Several studies have highlighted the relevance of lncRNAs in adipocyte differentiation, but the precise molecular mechanism has largely remained elusive. In the present study, we performed Ribo-Zero RNA-Seq to investigate both the poly(A)+and poly(A)-lncRNAs of in vitro cultured bovine preadipocytes and differentiated adipocytes. A stringent set of 2882 lncRNAs was finally identified. A comparison of the lncRNAs expression profiles revealed that 16 lncRNAs are differentially expressed during adipocyte differentiation. We focused on the most downregulated lncRNA, which we named adipocyte differentiation-associated long noncoding RNA (ADNCR). Mechanistically, ADNCR inhibited adipocyte differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-204, thereby augmenting the expression of the miR-204 target gene, SIRT1, which is known to inhibit adipocyte differentiation and adipogenic gene expression by docking with NCoR and SMART to repress PPARγ activity. Our data not only provide a valuable genomic resource for the identification of lncRNAs with functional roles in adipocyte differentiation but also reveal new insights into understanding the mechanisms of adipogenic differentiation.

Muscle development is regulated under a series of complicate processes, and non-coding RNAs, such as microRNAs (miRNAs) and circular RNAs (circRNAs), have been reported to play important roles in regulating myoblast proliferation and differentiation. We found that miR-107 expression was high in skeletal muscle of Qinchuan cattle. Overexpression of miR-107 inhibited bovine myoblasts differentiation and protected cells from apoptosis. Wnt3a was identified as a target of miR-107 by luciferase activity, real-time qPCR, and western blotting assays. Knockdown of Wnt3a inhibited bovine myoblasts differentiation and apoptosis, and this effect was similar to miR-107 overexpression. We also found circFGFR4 to promote myoblasts differentiation and to induce cell apoptosis. Via luciferase screening and RNA pull-down assays, circFGFR4 was observed to sponge miR-107. Overexpression of circFGFR4 increased the expression of Wnt3a, whereas this effect was abolished by miR-107. These results demonstrated that circFGFR4 binding miR-107 promotes cell differentiation via targeting Wnt3a in bovine primary myoblasts.

The Mongolian cattle are one of the most widespread breeds with strictly Bos taurus morphological features in northern China. In our current study, we presented a diversity of mitochondrial DNA (mtDNA) D-loop region and Y chromosome SNP markers in 25 male and 8 female samples of Mongolian cattle from the Xinjiang Uygur autonomous region in Western China, and detected 21 B. taurus and four Bos indicus (zebu) mtDNA haplotypes. Among four B. indicus mtDNA haplotypes, two haplotypes belonged to I1 haplogroup and the remaining two haplotypes belonged to I2 haplogroup. In contrast, all 25 male Mongolian cattle samples revealed B. taurus Y chromosome haplotype and no B. indicus haplotypes were found. Historical and archeological records indicate that B. taurus was introduced to Xinjiang during the second millennium BC and B. indicus appeared in this region by the second century AD. The two types of cattle coexisted for many centuries in Xinjiang, as depicted in clay and wooden figurines unearthed in the Astana cemetery in Turfan (3rd-8th century AD). Multiple lines of evidence suggest that the earliest B. indicus introgression in the Mongolian cattle may have occurred during the 2nd-7th centuries AD through the Silk Road around the Xinjiang region. This conclusion differs from the previous hypothesis that zebu introgression to Mongolian cattle happened during the Mongol Empire era in the 13th century.

Forkhead box A2 (Foxa2) has been recognized as one of the most potent transcriptional activators that is implicated in the control of feeding behavior and energy homeostasis. However, similar researches about the effects of genetic variations of Foxa2 gene on growth traits are lacking. Therefore, this study detected Foxa2 gene polymorphisms by DNA pool sequencing, PCR-RFLP and PCR-ACRS methods in 822 individuals from three Chinese cattle breeds. The results showed that four sequence variants (SVs) were screened, including two mutations (SV1, g. 7005 C>T and SV2, g. 7044 C>G) in intron 4, one mutation (SV3, g. 8449 A>G) in exon 5 and one mutation (SV4, g. 8537 T>C) in the 3'UTR. Notably, association analysis of the single mutations with growth traits in total individuals (at 24months) revealed that significant statistical difference was found in four SVs, and SV4 locus was highly significantly associated with growth traits throughout all three breeds (P<0.05 or P<0.01). Meanwhile, haplotype combination CCCCAGTC also indicated remarkably associated to better chest girth and body weight in Jiaxian Red cattle (P<0.05). We herein described a comprehensive study on the variability of bovine Foxa2 gene that was predictive of molecular markers in cattle breeding for the first time.

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Circular RNAs (circRNAs) have been identified from various tissues and species, but their regulatory functions during developmental processes are not well understood. We examined circRNA expression profiles of two developmental stages of bovine skeletal muscle (embryonic and adult musculus longissimus) to provide first insights into their potential involvement in bovine myogenesis. We identified 12 981 circRNAs and annotated them to the Bos taurus reference genome, including 530 circular intronic RNAs (ciRNAs). One parental gene could generate multiple circRNA isoforms, with only one or two isoforms being expressed at higher expression levels. Also, several host genes produced different isoforms when comparing development stages. Most circRNA candidates contained two to seven exons, and genomic distances to back-splicing sites were usually less than 50 kb. The length of upstream or downstream flanking introns was usually less than 105 nt (mean≈11 000 nt). Several circRNAs differed in abundance between developmental stages, and real-time quantitative PCR (qPCR) analysis largely confirmed differential expression of the 17 circRNAs included in this analysis. The second part of our study characterized the role of circLMO7-one of the most down-regulated circRNAs when comparing adult to embryonic muscle tissue-in bovine muscle development. Overexpression of circLMO7 inhibited the differentiation of primary bovine myoblasts, and it appears to function as a competing endogenous RNA for miR-378a-3p, whose involvement in bovine muscle development has been characterized beforehand. Congruent with our interpretation, circLMO7 increased the number of myoblasts in the S-phase of the cell cycle and decreased the proportion of cells in the G0/G1 phase. Moreover, it promoted the proliferation of myoblasts and protected them from apoptosis. Our study provides novel insights into the regulatory mechanisms underlying skeletal muscle development and identifies a number of circRNAs whose regulatory potential will need to be explored in the future.

Paired box 7 (PAX7) gene regulates the conversion of muscle satellite cells into myogenic cells and participates in multi-step processes in myogenesis. Expression levels of PAX7 are decisive for its regulatory function. Previous reports revealed that PAX7 were responsible for the developmental traits of muscle. The relationship of the PAX7 promoter variants and livestock phenotypic traits has not been fully elucidated. We detected a novel 10-bp insertion/deletion (indel) polymorphism in the bovine PAX7 promoter and revealed that the indel altered the binding of the transcriptional factor ZNF219. Luciferase reporter assay showed that deletion-deletion (Del-Del) genotype of the PAX7 gene showed 2.79-fold higher promoter activity than the insertion-insertion (Ins-Ins) genotype (P < 0.05), and ZNF219 overexpression significantly diminished the luciferase activity in Ins-Ins groups. Moreover, the expression of PAX7 and its down-stream genes were detected in fetal skeletal muscle of cattle with different PAX7 genotypes, where the Del-Del genotype also displayed high expression levels. Statistical association analysis demonstrated that this indel had significant effects on early growth traits in cattle. These findings provide a complete overview of the function of the PAX7 10-bp variant, which may have potential as a genetic marker for marker-assisted selection in improving economically significant traits of cattle.

Yanbian cattle is inhabitant of North of China, exhibiting many phenotypic features, such as long, dense body hair, and abundant intramuscular fat, designed to combat the extreme cold climate adaption. In the current study, we studied the cold tolerance of nine Yanbian cattle by whole genome resequencing and compared with African tropical cattle, N'Dama, as a control group. Yanbian cattle was aligned to the Bos taurus reference genome (ARS-UCD1.2) yielding an average of 10.8 fold coverage. The positive selective sweep analysis for the cold adaption in Yanbian cattle were analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 292 overlapped genes. The strongest selective signal was found on BTA16 with potential mutation in CORT gene, a regulatory gene of primary hormone in the hypothalamic-pituitary-adrenal (HPA) axis, is reported to be associated with the cold stress, representedfour missense mutations (c.269C > T, p.Lys90Ile; c.251A > G, p.Glu84Gly; c.112C > T, p.Pro38Ser; c.86G > A, p.Pro29His). Meanwhile another gene on BTA6, showed significantly higher selective sweep signals for a cold adapted trait for hair follicle and length development, FGF5 (fibroblast growth factor 5) with a missense mutation (c.191C > T, p.Ser64Phe). Moreover, cold adapted Yanbian cattle was statistically compared with the hot adapted N'Dama cattle, a taurine cattle reported to show superior heat tolerance than zebu cattle, making them better adapted to the hot regions of Africa. XP-CLR, Fst, and θπ ratio were used to compare both breeds, yielding 487, 924, and 346 genes respectively. Among the 12 overlapped genes, (CD36) (c.638A > G, p.Lys 213Arg) involved in fat digestion and absorption plays an important role in membrane transport of long-chain fatty acid and its expression could increase in cold exposure. Henceforth, our study provides a novel genetic insights into the cold climate adaptation of Yanbian cattle and identified three candidate genes (CORT, FGF5, and CD36), which can add to an understanding of the cold climate adaptation of Yanbian cattle.

A new direction for the treatment of skeletal myopathies, which are mainly caused by abnormal mitochondrial metabolism, is the application of drugs and active substances to relieve oxidative stress in mitochondria. Phloretin, a dihydrochalcone active substance widely present in succulent fruits, has attracted attention for its strong antioxidant activity. This study aimed to investigate the potential antioxidant effects of phloretin and its potential mechanism of action in C2C12 mouse myoblasts. Under oxidative stress caused by 500 μmol/L H2O2, the addition of 10 μmol/L phloretin ameliorated the high level of reactive oxygen species, increased CuZn/Mn-dependent superoxide dismutase activities, and restored the loss of mitochondrial membrane potential. Additionally, apoptosis, necrocytosis, and the inhibition of cell proliferation caused by H2O2 stimulation were alleviated by phloretin. Moreover, phloretin significantly increased the expression of cyclin D1 and alleviated the stagnation trend of the G1 phase of cell proliferation caused by H2O2. Furthermore, the addition of phloretin simultaneously significantly increased the protein and mRNA expression of heme oxygenase-1 (HO-1) and alleviated the inhibitory phosphorylation of p-nuclear factor erythroid 2-related factor 2 (Nrf2), p-AMP-activated protein kinase (AMPK), and p-liver kinase B1 (LKB1) induced by H2O2. Moreover, the expression of nuclear Nrf2 was higher with phloretin treatment than without phloretin treatment. Overall, phloretin alleviated the proliferation inhibition and apoptosis induced by H2O2 and exerted antioxidant effects via the LKB1/AMPK/Nrf2/HO-1 pathway in C2C12 cells. These results provide insight for the application of phloretin to alleviate oxidative damage to muscle.

Adipogenesis is a complex cellular process, which needs a series of molecular events, including long non-coding RNA (lncRNA). In the present study, a novel lncRNA named BADLNCR1 was identified as a regulator during bovine adipocyte differentiation, which plays an inhibitory role in lipid droplet formation and adipogenic marker gene expression. CHIPR-seq data demonstrated a potential competitive binding motif between BADLNCR1 and sterol regulatory element-binding proteins 1 and 2 (SREBP1/2). Dual-luciferase reporter assay indicated target relationship between KLF2 and BADLNCR1. Moreover, after the induction of KLF2, the expression of adipogenic gene reduced, while the expression of BADLNCR1 increased. Real-time quantitative PCR (qPCR) showed that BADLNCR1 negatively regulated mRNA expression of GLRX5 gene, a stimulator of genes that promoted formation of lipid droplets and expression of adipogenic genes. GLRX5 could partially reverse the effect of BADLNCR1 in bovine adipocyte differentiation. Dual-luciferase reporter assay stated that BADLNCR1 significantly reduced the enhancement of C/EBPα on promoter activity of GLRX5 gene. Furthermore, CHIP-PCR and CHIRP-PCR confirmed the suppressing effect of BADLNCR1 on binding of C/EBPα to GLRX5 promoter. Collectively, this study revealed the molecular mechanisms underlying the negative regulation of BADLNCR1 in bovine adipogenic differentiation.

Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.

Testes produce sperm, and investigations into gene expression in the testes will enhance the understanding of the roles of testicular genes in male reproduction. Cryptorchidism, the failure of one or both testes to descend into the scrotal sac, is a common congenital malformation in horses. The major clinical consequence of this abnormality is impaired fertility. The aim of this study was to analyze the expression patterns of testicular genes and to identify the differentially expressed genes (DEGs) in testes between cryptorchid and normal horses. In this study, the gene expression patterns in equine testes and the DEGs between mature descended testes (DTs) and undescended testes (UDTs) were identified by RNA-seq and validated by real-time qPCR. Our results provide comprehensive transcriptomic data on equine testes. The transcriptomic analysis revealed 11 affected genes that were downregulated in UDTs, possibly as a result of the higher temperature in the abdomen than in the scrotal sac. These 11 genes have previously been associated with male reproduction, and their downregulation might explain the impaired fertility of cryptorchid horses. Two homozygous missense mutations detected in horses with cryptorchidism were absent in normal horses and were listed as potential pathogenic mutations; these mutations should be verified in the future.

Copy number variation is a part of genomic structural variation and has caused widespread concern. According to the results of high-throughput screening of the MLLT10 gene, we found that the copy number variation region of the MLLT10 gene was correlated with bovine growth traits. We aimed to detect the MLLT10 gene copy number variation and provide materials for the Chinese yellow cattle breed. In this study, the SPSS software was used to analyze the correlation among the copy number type of six different cattle breeds (i.e., Qinchuan, Xianan, Jiaxian, Yanbian, Sinan, Yunling) and the corresponding growth traits. The results showed the following: In Qinchuan cattle, the copy number duplication type was greater than the deletion and normal types; in Xianan cattle, the copy number duplication and normal types were less as compared with the deletion type; and in Yunling cattle, the frequency of the duplication type was dominant among the three types of copy number variants. The correlation analysis result showed that there is a significant correlation between the copy number variation (CNV) of the MLLT10 gene and the growth traits of three cattle breeds. Furthermore, correlation analysis showed that MLLT10 CNV had positive effects on growth traits such as hip width, rump length, hucklebone width, and cannon bone circumference (p < 0.05). This study provides a basis for the molecular-assisted marker breeding of cattle and contributes to the breeding of cattle.

Previous, studies have shown that the dynein transporter compound has a role in diseases such as intellectual disability and cerebral malformations. However, the study of CNV in DYNC1I2 gene has not been reported. Q-PCR and data association analysis were used for DYNC1I2 gene copy in this study. In this study, blood samples were collected from five breeds of Chinese cattle (Qingchuan cattle, Xianan cattle, Yunling cattle, Pinan cattle and Guyuan cattle) for DYNC1I2 gene CNV type detection. SPSS 20.0 software and method of ANOVA were used to analyzed the association between types of CNV and growth traits. Results reveal that the distribution of different copy number types in different cattle breeds is different. Association analysis indicate that CNV of DYNC1I2 gene showed a positive effect in cattle growth: in XN cattle, individuals with deletion types showed better performance on height at hip cross (P < 0.05); individuals with duplication types have better performance on body length (P < 0.05) in PN cattle; individuals with deletion types was significantly correlated with chest width and Hucklebone width (P < 0.05) in QC cattle; individuals with duplication types in Yunling cattle were better than the normal types, and there was a significant correlation between copy number variant and chest depth (P < 0.05). The results showed that CNV markers closely related to cattle production traits were detected at DNA level, which could be used as an important candidate molecular marker for marker-assisted selection of growth traits in Chinese cattle, and provided a new research basis for genetics and breeding of Chinese beef cattle.

As a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the function of circRNAs in bovine muscle development. Based on existing sequencing data, we identified circINSR. The localization of circINSR in bovine myoblasts was investigated by fluorescence in situ hybridization. Molecular and biochemical assays were used to confirm the role of circINSR in myoblast proliferation and the cell cycle. Mitochondrial membrane potential and annexin V-PE/7-AAD staining assays were performed to assess cell apoptosis. Additionally, interactions between circINSR, miR-34a, and target mRNAs were examined using bioinformatics, a luciferase assay, and RNA immunoprecipitation. We found that circINSR was highly expressed in embryonic muscle tissue. Overexpression of circINSR significantly promoted proliferation and reduced apoptosis of embryonic myoblasts. Our data suggested that circINSR may act as a sponge of miR-34a and could function through de-repression of target genes in muscle cells. This study proposes that circINSR may function as a regulator of embryonic muscle development. circINSR regulates cells proliferation and apoptosis through miR-34a-modulated Bcl-2 and CyclinE2 expression.

The Longlin goat is one of the most valuable livestock species in Guangxi Autonomous Region of China, but its genomic diversity and selective signals are not clearly elucidated. Here we compared 20 genomes of Longlin goat to 66 genomes of other seven goat breeds worldwide to analyze patterns of Longlin goat genetic variation. We found the lowest linkage disequilibrium at the large distances between SNPs associated with the highest effective population size in the recent generations ago in Longlin goat. The eight goat breeds could be divided into Euro-African and East Asian goat population. Interestingly, like East Asian taurine, the same two migration phases might have occurred in the history of East Asian goat. More importantly, we identified selective signals implicated in immune resistance to disease, especially for skin disease, in Longlin goat. Our findings will not only help understand the evolutionary history and breed characteristic but can provide valuable resources for conservation of germplasm resources and implementation of crossbreeding programs.

Generally, copy number variation (CNV) is a large-scale structural variation between 50 bp and 1 kb of the genome. It can affect gene expression and is an important reason for genetic diversity and phenotypic trait diversity. Studies have shown that the eukaryotic translation initiation factor 4A2 (EIF4A2) gene plays an essential role in muscle development in both humans and pigs. However, the influence of bovine EIF4A2's copy number change on phenotypic traits has not been reported.

Domestic yak, cattle, and their hybrids are fundamental to herder survival at high altitudes on the Tibetan Plateau. However, little is known about their history. Bos remains are uncommon in this region, and ancient domestic yak have not been securely identified. To identify Bos taxa and investigate their initial management, we conducted zooarchaeological analyses of 193 Bos specimens and sequenced five nuclear genomes from recently excavated assemblages at Bangga. Morphological data indicated that more cattle than yak were present. Ancient mitochondrial DNA and nuclear genome sequences identified taurine cattle and provided evidence for domestic yak and yak-cattle hybridization ~2500 years ago. Reliance on diverse Bos species and their hybrid has increased cattle adaptation and herder resilience to plateau conditions. Ancient cattle and yak at Bangga were closely related to contemporary livestock, indicating early herder legacies and the continuity of cattle and yak husbandry on the Tibetan Plateau.

Hirschsprung's disease (HSCR) is a congenital disease in neonates characterized by the absence of the enteric ganglia in a variable length of the distal colon. This disease results from multiple genetic interactions that modulate the ability of enteric neural crest cells to populate developing gut. We previously reported that three rat strains with different backgrounds (susceptible AGH-Ednrbsl/sl, resistant F344-Ednrbsl/sl, and LEH-Ednrbsl/sl) but the same null mutation of Ednrb show varying severity degrees of aganglionosis. This finding suggests that strain-specific genetic factors affect the severity of HSCR. Consistent with this finding, a quantitative trait locus (QTL) for the severity of HSCR on chromosome (Chr) 2 was identified using an F2 intercross between AGH and F344 strains. In the present study, we performed QTL analysis using an F2 intercross between the susceptible AGH and resistant LEH strains to identify the modifier/resistant loci for HSCR in Ednrb-deficient rats. A significant locus affecting the severity of HSCR was also detected within the Chr 2 region. These findings strongly suggest that a modifier gene of aganglionosis exists on Chr 2. In addition, two potentially causative SNPs (or mutations) were detected upstream of a known HSCR susceptibility gene, Gdnf. These SNPs were possibly responsible for the varied length of gut affected by aganglionosis.

Abundant miRNAs have been identified in milk and mammary gland tissues of different species. Typically, RNA in milk can be extracted from different fractions including fat, whey and cells and the mRNA transcriptome of milk could serve as an indicator of the transcriptome of mammary gland tissue. However, it has not been adequately validated if the miRNA transcriptome of any milk fraction could be representative of that of mammary gland tissue. The objectives of this study were to (1) characterize the miRNA expression spectra from three milk fractions- fat, whey and cells; (2) compare miRNome profiles of milk fractions (fat, whey and cells) with mammary gland tissue miRNome, and (3) determine which milk fraction miRNome profile could be a better representative of the miRNome profile of mammary gland tissue. Milk from four healthy Canadian Holstein cows in mid lactation was collected and fractionated. Total RNA extracted from each fraction was used for library preparation followed by small RNA sequencing. In addition, miRNA transcripts of mammary gland tissues from twelve Holstein cows in our previous study were used to compare our data. We identified 210, 200 and 249 known miRNAs from milk fat, whey and cells, respectively, with 188 universally expressed in the three fractions. In addition, 33, 31 and 36 novel miRNAs from milk fat, whey and cells were identified, with 28 common in the three fractions. Among 20 most highly expressed miRNAs in each fraction, 14 were expressed in common and 11 were further shared with mammary gland tissue. The three milk fractions demonstrated a clear separation from each other using a hierarchical cluster analysis with milk fat and whey being most closely related. The miRNome correlation between milk fat and mammary gland tissue (rmean = 0.866) was significantly higher than the other two pairs (p < 0.01), whey/mammary gland tissue (rmean = 0.755) and milk cell/mammary gland tissue (rmean = 0.75), suggesting that milk fat could be an alternative non-invasive source of RNA in assessing miRNA activities in bovine mammary gland. Predicted target genes (1802) of 14 highly expressed miRNAs in milk fractions were enriched in fundamental cellular functions, infection, organ and tissue development. Furthermore, some miRNAs were highly enriched (FDR <0.05) in milk whey (3), cells (11) and mammary gland tissue (14) suggesting specific regulatory functions in the various fractions. In conclusion, we have obtained a comprehensive miRNA profile of the different milk fractions using high throughput sequencing. Our comparative analysis showed that miRNAs from milk fat accurately portrayed the miRNome of mammary gland tissue. Functional annotation of the top expressed miRNAs in milk confirmed their critical regulatory roles in mammary gland functions and potentially to milk recipients.