RNA-seq analysis identifies differentially expressed gene in different types of donkey skeletal muscles.

The main component of donkey meat is skeletal muscle, and different muscle fibers have been found to be associated with meat quality. However, RNA-seq technology has yet to be used to profile the transcriptomic changes of different muscles of the donkey. In this study, the characterizations of different muscles on the gene expression profiles of Dezhou donkey were obtained, the aim was to identify the important genes in donkey muscles, and aid in improving donkey meat quality via RNA-seq. In the donkey gluteus (DG) and donkey longissimus dorsi (DL) group, GO enrichment indicated that DEGs were mainly involved in the biological regulation and metabolic process, and KEGG analysis shows that a total of 427 DEGs were mapped to 216 KEGG pathways and 23 KEGG pathways were significantly enriched such as the ribosome, glycolysis/gluconeogenesis, glucagon signaling pathway and biosynthesis of amino acids pathways. Meanwhile, 504 DEGs were mapped to 223 KEGG pathways, in which 17 were significantly enriched including cardiac muscle contraction and oxytocin signaling pathway in donkey hamstring muscles (DH) and DL group. In addition, the tenderness in donkey meat might involve muscle fiber type and glucose metabolism, which might profit from the DEGs including MYH1, MYH7, TNNC1, TNNI3, TPM3, ALDOA, ENO3, and PGK1. The genes found in this study will provide some ideas for further understanding the molecular mechanism of donkey meat quality.

Macrophages Rapidly Seal off the Punctured Zebrafish Larval Brain through a Vital Honeycomb Network Structure.

There is accumulating evidence that macrophages play additional important roles in tissue damage besides their typical phagocytosis. Although the aggregation of macrophages on injured sites has long been observed, few researchers have focused on the role of the overall structure of macrophage aggregation. In this study, we developed a standardized traumatic brain injury (TBI) model in zebrafish larvae to mimic edema and brain tissue spillage symptoms after severe brain trauma. Using time-lapse imaging, we showed that macrophages/microglia in zebrafish larvae responded rapidly and dominated the surface of injured tissue, forming a meaningful honeycomb network structure through their compact aggregation and connection. Disrupting this structure led to fatal edema-like symptoms with severe loss of brain tissue. Using the RNA-Seq, together with the manipulation of in vitro cell lines, we found that collagen IV was indispensable to the formation of honeycomb network structures. Our study thus revealed a novel perspective regarding macrophages forming a protective compact structure with collagen IV. This honeycomb network structure acted as a physical barrier to prevent tissue loss and maintain brain homeostasis after TBI. This study may provide new evidence of macrophages' function for the rapid protection of brain tissue after brain injury.

Ikzf1 regulates embryonic T lymphopoiesis via Ccr9 and Irf4 in zebrafish.

Ikzf1 is a Krüppel-like zinc-finger transcription factor that plays indispensable roles in T and B cell development. Although the function of Ikzf1 has been studied extensively, the molecular mechanism underlying T lymphopoiesis remains incompletely defined during the embryonic stage. Here we report that the genetic ablation of ikzf1 in mutant zebrafish resulted in abrogated embryonic T lymphopoiesis. This was ascribed to impaired thymic migration, proliferation, and differentiation of hematopoietic stem/progenitor cells (HSPCs). Ccr9a and Irf4a, two indispensable factors in T lymphopoiesis, were the direct targets of Ikzf1 and were absent in the ikzf1 mutants. Genetic deletion of either ccr9a or irf4a in the corresponding mutant embryos led to obvious T cell development deficiency, which was mainly caused by disrupted thymic migration of HSPCs. Restoration of ccr9a in ikzf1 mutants obviously promoted HSPC thymus homing. However, the HSPCs then failed to differentiate into T cells. Additional replenishment of irf4a efficiently induced HSPC proliferation and T cell differentiation. Our findings further demonstrate that Ikzf1 regulates embryonic T lymphopoiesis via Ccr9 and Irf4 and provide new insight into the genetic network of T lymphocyte development.

Utilizing mobile digital radiography for detection of thoracolumbar vertebrae traits in live donkeys.

It has been well-established that the number of vertebrae is associated with body size and meat productivity. In current study we utilized a digital radiography (DR) technology to detect the number of thoracolumbar vertebrae in live donkeys. For this purpose, we introduced for the first time a groundbreaking device designed by our team for assessing thoracolumbar vertebrae number traits in equids, employing a sample of 1,000 donkeys sourced from five distinct donkey farms. This assessment incorporates a range of crucial body metrics, including body height, length, and various other measurements. Subsequently, our study determined the number of thoracolumbar vertebrae in 112 donkeys, utilizing the DR system. These findings were further validated through post-mortem evaluations conducted by slaughtering the donkeys. Our findings demonstrated a remarkable resemblance between the thoracolumbar vertebrae numbers visualized through the DR system in live donkeys and those obtained via slaughter verification. In conclusion, this research underscores the accuracy and effectiveness of the DR system for the detection of thoracolumbar vertebrae in live donkeys, which might be helpful for assessing the body size and meat productivity. We also recommended the utilization of DR system for counting thoracolumbar vertebrae in other animals in live state and could be a useful addition to livestock business industry for the prediction of body size and meat productivity efficiency.

Elucidating the Role of Transcriptomic Networks and DNA Methylation in Collagen Deposition of Dezhou Donkey Skin.

DNA methylation represents a predominant epigenetic modification with broad implications in various biological functions. Its role is particularly significant in the process of collagen deposition, a fundamental aspect of dermal development in donkeys. Despite its critical involvement, the mechanistic insights into how DNA methylation influences collagen deposition in donkey skin remain limited. In this study, we employed whole genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) to investigate the epigenetic landscape and gene expression profiles in the dorsal skin tissues of Dezhou donkeys across three developmental stages: embryonic (YD), juvenile (2-year-old, MD), and mature (8-year-old, OD). Our analysis identified numerous differentially methylated genes that play pivotal roles in skin collagen deposition and overall skin maturation, including but not limited to COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, GLUL, SFRP2, FOSL1, SERPINE1, MMP1, MMP2, MMP9, and MMP13. Notably, we observed an inverse relationship between gene expression and DNA methylation proximal to transcription start sites (TSSs), whereas a direct correlation was detected in regions close to transcription termination sites (TTSs). Detailed bisulfite sequencing analyses of the COL1A1 promoter region revealed a low methylation status during the embryonic stage, correlating with elevated transcriptional activity and gene expression levels. Collectively, our findings elucidate key genetic markers associated with collagen deposition in the skin of Dezhou donkeys, underscoring the significant regulatory role of DNA methylation. This research work contributes to the foundational knowledge necessary for the genetic improvement and selective breeding of Dezhou donkeys, aiming to enhance skin quality attributes.

Comprehensive transcriptomic analysis unveils the interplay of mRNA and LncRNA expression in shaping collagen organization and skin development in Dezhou donkeys.

The primary focus of donkey hide gelatin processing lies in the dermal layer of donkey hide due to its abundant collagen content. However, the molecular mechanism involved in collagen organization and skin development in donkey skin tissue across various developmental stages remains incomplete. The current study aims to investigate the transcriptomic screening of lncRNAs and mRNA associated with skin development and collagen organization across different ages in Dezhou donkeys' skin. In the pursuit of this objective, we used nine skin tissue samples obtained from Dezhou donkeys at various ages including 8-month fetal stage, followed by 2 and 8 years. RNA-seq analysis was performed for the transcriptomic profiling of differentially expressed genes (DEGs) and lncRNAs associated with skin development in different age groups. Our investigation revealed the presence of 6,582, 6,455, and 405 differentially expressed genes and 654, 789, and 29 differentially expressed LncRNAs within the skin tissues of Dezhou donkeys when comparing young donkeys (YD) vs. middle-aged donkeys (MD), YD vs. old donkeys (OD), and MD vs. OD, respectively. Furthermore, we identified Collagen Type I Alpha 1 Chain (COL1A1), Collagen Type III Alpha 1 Chain (COL3A1), and Collagen Type VI Alpha 5 Chain (COL6A5) as key genes involved in collagen synthesis, with COL1A1 being subject to cis-regulation by several differentially expressed LncRNAs, including ENSEAST00005041187, ENSEAST00005038497, and MSTRG.17248.1, among others. Interestingly, collagen organizational and skin development linked pathways including Protein digestion and absorption, metabolic pathways, Phosphatidylinositol 3-Kinase-Protein Kinase B signaling pathway (PI3K-Akt signaling pathway), Extracellular Matrix-Receptor Interaction (ECM-receptor interaction), and Relaxin signaling were also reported across different age groups in Dezhou donkey skin. These findings enhance our comprehension of the molecular mechanisms underlying Dezhou donkey skin development and collagen biosynthesis and organization, thus furnishing a solid theoretical foundation for future research endeavors in this domain.

A novel SNP in NKX1-2 gene is associated with carcass traits in Dezhou donkey.

BACKGROUND: At present, donkey meat in the market shows an imbalance between supply and demand, and there is an urgent need to cultivate a meat-type Dezhou donkey breed. On the one hand, it can improve the imbalance in the market, and on the other hand, it can promote the rapid development of the donkey industry. This study aimed to reveal significant genetic variation in the NK1 homeobox 2 gene (NKX1-2) of Dezhou donkeys and investigate the association between genotype and body size in Dezhou donkeys. RESULTS: In this study, a SNP (g.54704925 A > G) was identified at the exon4 by high-depth resequencing of the Dezhou donkey NKX1-2 gene. The AA genotype is the dominant genotype. The g.54704925 A > G site was significantly associated with body length, thoracic girth, and hide weight (P < 0.05), while it was highly significantly associated with body height and carcass weight (P < 0.01) in Dezhou donkeys. CONCLUSION: Overall, the results of this study showed that the NKX1-2 gene could be a candidate gene for breeding meat-type Dezhou donkeys, and the g.54704925 A > G locus could be used as a marker locus for selection and breeding.

Polymorphism detection of PRKG2 gene and its association with the number of thoracolumbar vertebrae and carcass traits in Dezhou donkey.

BACKGROUND: Previous studies have shown that the protein kinase cGMP-dependent 2 (PRKG2) gene is associated with dwarfism in humans, dogo Argentines, and Angus cattle, as well as with height and osteoblastogenesis in humans. Therefore, the PRKG2 gene was used as the target gene to explore whether this gene is associated with several thoracolumbar vertebrae and carcass traits in Dezhou donkeys. RESULTS: In this study, fifteen SNPs were identified by targeted sequencing, all of which were located in introns of the PRKG2 gene. Association analysis illustrated that the g.162153251 G > A, g.162156524 C > T, g.162158453 C > T and, g.162163775 T > G were significantly different from carcass weight. g.162166224 G > A, g.162166654 T > A, g.162167165 C > A, g.162167314 A > C and, g.162172653 G > C were significantly associated with the number of thoracic vertebrae. g.162140112 A > G was significantly associated with the number and the length of lumbar vertebrae, and g.162163775 T > G was significantly associated with the total number of thoracolumbar vertebrae. CONCLUSION: Overall, the results of this study suggest that PRKG2 gene polymorphism can be used as a molecular marker to breed high-quality Dezhou donkeys.

Metabolic Alternations During Gestation in Dezhou Donkeys and the Link to the Gut Microbiota.

The maternal intestinal microbial community changes dramatically during pregnancy and plays an important role in animal growth, metabolism, immunity and reproduction. However, our understanding of microbiota compositional dynamics during the whole pregnancy period in donkey is incomplete. This study was carried out to evaluate gut microbiota alterations as well as the correlation with serum biochemical indices, comparing pregnant donkeys to non-pregnant donkeys. A total of 18 pregnant (including EP, early-stage pregnancy; MP, middle-stage pregnancy and LP, late-stage pregnancy) and six non-pregnant (C as a control) donkey blood samples and rectum contents were collected. The results showed that pregnant donkeys had higher microbial richness than non-pregnant donkeys and that the lowest microbial diversity occurred at the EP period. Moreover, the relative abundances of the families Clostridiaceae and Streptococcaceae were significantly higher in the EP group (p < 0.05) than that in the C and MP groups, while the relative abundances of the families Lachnospiraceae and Rikenellaceae were significantly lower in the EP group (p < 0.05) than that in the C group. The predicted microbial gene functions related to the inflammatory response and apoptosis, such as Staphylococcus aureus infection, the RIG-1-like receptor signaling pathway and apoptosis, were mainly enriched in EP. Furthermore, pregnant donkeys had higher glucose levels than non-pregnant donkeys, especially at EP period. EP donkeys had lower triglyceride, total protein and albumin levels but higher malondialdehyde, interleukin 1ß, interleukin 6 and tumor necrosis factor-α levels than those in the C and MP groups. Additionally, there were strong correlations between inflammatory cytokine levels and the relative abundances of genera belonging to the Clostridiaceae and Streptococcaceae families. This is the first comparative study performed in donkeys that indicates that pregnancy status (especially in the early pregnancy period) alters the gut microbiota composition, which was correlated with serum biochemical parameters. These results could provide useful information for improving the reproductive management in Dezhou donkeys.

Dynamic alterations in the donkey fecal bacteria community and metabolome characteristics during gestation.

In donkeys, the gestation period is a dynamic and precisely coordinated process involving systemic and local alterations. Both the gut microbiota and its link with blood metabolites are thought to play significant roles in maintaining maternal health and supporting fetal development during the gestation period. This study was conducted to evaluate gut microbiota changes and the correlation with plasma metabolites in Dezhou donkeys during the gestation period. The donkeys were divided into the four following groups according to their pregnancy stages: the non-pregnant (NP), early stage of pregnancy (P1), middle stage of pregnancy (P2), and late stage of pregnancy (P3) groups. A total of 24 (n = 6 per group) samples of donkey feces and plasma were collected. The results showed that the diversity (Shannon index) of fecal bacteria significantly increased throughout the gestation period. The phyla Spirochaetota and Fibrobacterota varied significantly according to the stages of pregnancy (p < 0.05). At the genus level, the abundance of Treponema in pregnant donkeys was greater than that in non-pregnant donkeys (p < 0.05), and the genus Streptococcus reached its maximum abundance in the P2 period (p < 0.05). The abundance of Ruminococcaceae_NK4A214_group and norank_f_norank_o_WCHB1-41 linearly increased with the progression of pregnancy (p < 0.05). In addition, the host plasma metabolome was altered significantly during the gestation period. Testolic acid, estradiol-17beta 3-sulfate, equol 7'-o-glucuronide, equol 4'-o-glucuronide, estrone, estrone 3-glucuronide, and estradiol were the most significant differential enriched metabolites, and they increased gradually as gestation progressed. The altered metabolites were mainly enriched in pathways matched to bile secretion, ABC transporters, amino acid metabolism, protein digestion and absorption, mineral absorption, fatty acid degradation, glycerophospholipid metabolism, and steroid hormone biosynthesis. We also found a significant correlation between the shifts in donkey fecal bacteria and changes in the host metabolism. In summary, this study provided systematic data on the fecal bacterial changes and host plasma metabolism of donkeys throughout pregnancy. The results indicated that host-bacteria interactions during the gestation period influence the host metabolism. These interactions benefit the pregnant donkeys by providing a sufficient supply of nutrients and energy for fetal growth and maternal health.

Multi-Thoracolumbar Variations and NR6A1 Gene Polymorphisms Potentially Associated with Body Size and Carcass Traits of Dezhou Donkey.

The number of thoracolumbar vertebrae is a quantitative trait positively correlated with the economic traits of livestock. More thoracolumbar vertebrae individuals could genetically be used to improve the livestock population, as more thoracolumbar vertebrae means a longer carcass, which could bring more meat production. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate gene for effecting the number of vertebrae in livestock. The purposes of this study are as follows: (a) Analyzing the effect of TLN variation on body size and carcass traits of Dezhou donkey; (b) Studying the distribution of seven single nucleotide variants (SNVs) in NR6A1 gene of Dezhou donkey; (c) Exploring the relationship between latent SNVs and TLN, the body size and carcass traits. We examined the thoracic and lumbar vertebrae number and seven SNVs in NR6A1 gene of 455 Dezhou donkeys, and analyzed the relationships between them. Five types of thoracolumbar combinations (T17L5 (individual with 17 thoracic and five lumbar vertebrae) 2.4%, T18L5 75.8%, T19L5 1.1%, T17L6 11.9%, and T18L6 8.8%) of Dezhou donkeys were detected in this study. For one thoracolumbar vertebra added, the body length of Dezhou donkey increases by 3 cm and the carcass weight increases by 6 kg. Seven SNVs (g.18093100G > T, g.18094587G > T, g.18106043G > T, g.18108764G > T, g.18110615T > G, g.18112000C > T and g.18114954T > G) of the NR6A1 gene were found to have a significant association with the TLN, body size and carcass traits of Dezhou donkey (p < 0.05), respectively. For instance, g.18114954C > T is significantly associated with lumber vertebrae number, the total number of thoracolumbar, and carcass weight, and individuals with TT genotype had significantly larger value than CC genotype (p < 0.05). Using these 7SNVs, 16 different haplotypes were estimated. Compared to Hap3Hap3, individuals homozygous for Hap2Hap2 showed significantly longer length in one thoracic spine (STL), the total thoracic vertebrae and one thoracolumbar spine. Our study will not only extend the understanding of genetic variation in the NR6A1 gene of Dezhou donkey, but also provide useful information for marker assisted selection in donkey breeding program.

Comparing the Profiles of Raw and Cooked Donkey Meat by Metabonomics and Lipidomics Assessment.

Heat cooking of meat gives it a specific taste and flavor which are favored by many consumers. While the characteristic taste components of chicken, duck, pig, and seafood have been studied, there is a lack of information about the molecular components that give donkey meat its unique taste. Here, the characterization profiles of raw donkey meat (RDM) and cooked donkey meat (CDM) meat by metabonomics and lipidomics. The results showed that a total of 186 metabolites belonging to 8 subclasses were identified in CDM and RDM, including carbohydrates (27.42%), amino acids (17.20%), lipids (13.44%), and nucleotides (9.14%). In total, 37 differential metabolites were identified between CDM and RDM. Among these, maltotriose, L-glutamate, and L-proline might predominantly contribute to the unique umami and sweet taste of donkey meat. Comprehensive biomarker screening detected 9 potential metabolite markers for the discrimination among RDM and CDM, including L-glutamate, gamma-aminobutyric acid, and butane-1, 2, 3, 4-tetrol. Moreover, a total of 992 and 1,022 lipids belonging to 12 subclasses were identified in RDM and CDM, respectively, mainly including triglycerides (TGs) and glycerophospholipids (GPs). Of these lipids, 116 were significantly different between CDM and RDM. The abundances of 61 TGs rich in saturated and monounsaturated fatty acids were retained in CDM, whereas the abundances of 37 GPs rich in polyunsaturated fatty acids were reduced, suggesting that TGs and GPs might be the predominant lipids for binding and generating aroma compounds, respectively. A total of 13 lipids were determined as potential markers for the discrimination among RDM and CDM, including PC(O-16:2/2:0), LPE(22:5/0:0), and PC(P-16:0/2:0). In conclusion, this study provided useful information about the metabolic and lipid profiles of donkey meat which may explain its unique taste and flavor, which could serve as a basis for the development and quality control of donkey meat and its products.

Estimation of Liveweight from Body Measurements through Best Fitted Regression Model in Dezhou Donkey Breed.

Dezhou donkey is an excellent local hide, meat and milk breed in Shandong Province. To accelerate the speed of breeding, reproduction and conversation, correlation and multiple regression analysis between body weight (BW) and body size of Dezhou donkey rearing under intensive farms was made by SAS 9.4 software (Statistical Analysis for Windows). A total of 162 donkeys of both gender of age 2 ~ 10 years old were used to investigate the relationships between BW and body dimensional traits (cm) including height of withers, body length; thoracic depth, thoracic girth (TG), thoracic width (TW), circumference of cannon bone (CB), height of rump, rump length (RL) and rump width (RW). The results showed that BW and body measurements have positive and great correlations with R2 value ranged from 0.58 to 0.88 (P <0.01). The R2 values from the single-parameter equations showed that the TG was highly related to BW (0.72, P < 0.01). The stepwise regression equations were applied to obtain the best prediction equations, and the results indicated that the prediction accuracy for BW was improved with the inclusion of more body measurement variables. The "best fit" models were: BW (kg) = 1.88 × TG + 1.27 × BL + 2.55 × TW + 4.61 × CB + 2.18 × RW + 1.78 × RL - 422.8 (R2 = 0.906, P < .01). The predicted BW from the present equations showed the nearest value to the real BW (R2 > 0.94, P < .01). In addition, the equations derived to predict the BW of donkeys in Britain, Morocco were less satisfactory for use with the present Dezhou donkey breed because they overestimated or underestimated the BW due to the different donkey breed.

The caudal dorsal artery generates hematopoietic stem and progenitor cells via the endothelial-to-hematopoietic transition in zebrafish.

Zebrafish hematopoietic stem and progenitor cells (HSPCs) originate from the hemogenic endothelium of the ventral wall of the dorsal aorta (DA) through the endothelial-to-hematopoietic transition (EHT) from approximately 30 to 60 hours post fertilization (hpf). However, whether other artery sites can generate HSPCs de novo remains unclear. In this study, using live imaging and lineage tracing, we found that the caudal dorsal artery (CDA) in the caudal hematopoietic tissue directly gave rise to HSPCs through EHT. This process initiated from approximately 60 hpf and terminated at approximately 156 hpf. Compared with that in the DA, fewer EHT events were observed in the CDA. The EHT events in the DA and CDA were similarly regulated by Runx1 but differentially influenced by blood flow (i.e., the EHT frequency in CDA was affected to a lesser extent when circulation was compromised in the tnnt2a-/- mutant). Therefore, the whole artery, including both DA and CDA, was endowed with the ability to produce HSPCs during a much longer time period. Coincidently, the lineage tracing results indicated that adult hematopoietic cells originated from the embryonic endothelium, and those produced later preferentially colonized the adult thymus. Collectively, our study revealed that the CDA serves as an additional source of hematopoiesis, and it shows similar but not identical properties with the DA.