Integrative analysis of long non-coding RNA and mRNA in broilers with valgus-varus deformity

Genome-wide methylation and transcriptome differential analysis of skeletal muscle in broilers with valgus-varus deformity

1. Valgus-varus deformity (VVD) is a disease that severely affects leg function in broilers and for which there is no effective control method current available. Although DNA methylation has an important impact on most physiological and pathological processes, its involvement in skeletal muscle growth and development in VVD broilers is unknown. In this study, genome-wide DNA methylation was analysed in VVD-affected and normal broilers using whole genome resulphite sequencing.2. The results showed that in the cytosine-phosphoric acid-guanine (CG) sequence environment there was a methylation rate of about 55% and 4,265 differentially methylated regions (DMRs) were found in the CG. Of these, 550 were located in the promoter, 547 in the exon region, and 1,718 in the intron region.3. All differentially methylated genes (DMGs) were analysed for enrichment of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The GO was enriched in pathways related to protein degradation such as proteasome complex, endopeptidase complex and extracellular region. The KEGG pathways were enriched in signalling pathways related to protein degradation and catabolism such as proteasome, nitrogen metabolism, adherens junction and alanine.4. Protein interactions analysis revealed that FOS, MYL9, and FRAS1 had a high degree of interactions, in which the DNA methylation level of the MYL9 promoter region was negatively correlated with mRNA expression level. Further studies showed that 5-azacytidine (5-AzaC) inhibited DNMT1 and DNMT3A gene expression and promoted MYL9 expression.5. This study systematically investigated overall DNA methylation patterns in the leg muscle of VVD and normal broilers. It screened common differential genes in conjunction with transcriptomic data to further identify genes associated with muscle growth and development. This study provides new insights to better understand the pathogenesis of VVD from an epigenetic perspective.

Gut-X Axis and Its Role in Poultry Bone Health: A Review

The normal development and growth of bones are critical for poultry health. With the rapid increase in poultry growth rates achieved over the last few decades, juvenile meat-type poultry exhibit a high incidence of leg weakness and lameness. These issues are significant contributors to poor animal welfare and substantial economic losses. Understanding the potential etiology of bone problems in poultry will aid in developing treatments for bone diseases. The gut microbiota represents the largest micro-ecosystem in animals and is closely related to many metabolic disorders, including bone disease. It achieves this by secreting secondary metabolites and coordinating with various tissues and organs through the circulatory system, which leads to the concept of the gut-X axis. Given its importance, modulating gut microbiota to influence the gut-X axis presents new opportunities for understanding and developing innovative therapeutic approaches for poultry bone diseases. In light of the extensive literature on this topic, this review focuses on the effects of gut microbiota on bone density and strength in poultry, both directly and indirectly, through the regulation of the gut-X axis. Our aim is to provide scientific insights into the bone health problems faced by poultry.

Functional analysis of a novel pathogenic variant in CREBBP associated with bone development

BACKGROUND

CREBBP has been extensively studied in syndromic diseases associated with skeletal dysplasia. However, there is limited research on the molecular mechanisms through which CREBBP may impact bone development. We identified a novel pathogenic CREBBP variant (c.C3862T/p.R1288W, which is orthologous to mouse c.3789 C > T/p.R1289W) in a patient with non-syndromic polydactyly.

METHODS

We created a homozygous Crebbp p.R1289W mouse model and compared their skeletal phenotypes to wild-type (WT) animals. Bone marrow stem cells (BMSCs) were isolated and assessed for their proliferative capacity, proportion of apoptotic cells in culture, and differentiation to chondrocytes and osteocytes.

RESULTS

We observed a significant decrease in body length in 8-week-old homozygous Crebbp p.R1289W mice. The relative length of cartilage of the digits of Crebbp p.R1289W mice was significantly increased compared to WT mice. BMSCs derived from Crebbp p.R1289W mice had significantly decreased cell proliferation and an elevated rate of apoptosis. Consistently, cell proliferative capacity was decreased and the proportion of apoptotic cells was increased in the distal femoral growth plate of Crebbp p.R1289W compared to WT mice. Chemical induction of BMSCs indicated that Crebbp p.R1289W may promote chondrocyte differentiation.

CONCLUSION

The Crebbp p.R1289W variant plays a pathogenic role in skeletal development in mice.

IMPACT

CREBBP has been extensively studied in syndromic diseases characterized by skeletal dysplasia. There is limited research regarding the molecular mechanism through which CREBBP may affect bone development. To our knowledge, we generated the first animal model of a novel Crebbp variant, which is predicted to be pathogenic for skeletal diseases. Certain pathogenic variants, such as Crebbp p.R1289W, can independently lead to variant-specific non-syndromic skeletal dysplasia.

Plasma exosome-derived miR-455-5p targets RPS6KB1 to regulate cartilage homeostasis in valgus-varus deformity (Gallus gallus)

Valgus-varus deformity (VVD) is a common long bone deformity in broilers. Imbalance in cartilage homeostasis is the main feature of leg disease. Exosomes act as an important intercellular communication vector that regulates chondrogenesis by encapsulating specific nucleic acids and proteins. However, the exact mechanism of how plasma exosomal miRNAs regulate cartilage homeostasis in VVD broilers remains unclear. This study first demonstrated the structural disorder, growth retardation, and reduced proliferative capacity of VVD cartilage in vitro and in vivo. Subsequently, VVD and Normal broiler plasma exosomes were collected for miRNA sequencing. Cartilage-specific miR-455-5p was extraordinarily emphasized by performing bioinformatics analysis on differential miRNA target genes and further validated by tissue expression profiling. PKH67 fluorescently labeled plasma exosomes were shown to be taken up by chondrocytes, deliver miR-455-5p, inhibit chondrocyte proliferation, and disrupt their homeostasis, and these effects could be inhibited by the miR-inhibitors. Mechanistically, MiR-455-5p targets Ribosomal Protein S6 Kinase B1 (RPS6KB1) to inhibit RPS6 phosphorylation and reduce the synthesis of key proteins for cartilage proliferation, which in turn inhibits cartilage proliferation and disrupts its homeostasis. In conclusion, the present study identified abnormalities in VVD cartilage tissue and clarified the specific mechanism by which plasma exosome-derived miR-455-5p regulates cartilage homeostasis.

Transcriptomics reveals the molecular regulation of Chinese medicine formula on improving bone quality in broiler

Skeletal disorder is of concern to the poultry industry as it affects animal welfare and production performance. Traditional Chinese medicine could improve bone quality and reduce the incidence of bone disease, but the molecular regulation of Chinese medicine formula (CMF) on improving bone quality in broilers is still unclear. This study was performed to research the effects of CMF on skeletal performance of Cobb broilers and reveal the molecular regulation. A total of 120 one-day-old Cobb broilers were randomly allocated into 4 equal groups of 30 chickens, with 5 replicates and 6 chickens in each replicate. The control (CON) group was fed a diet without CMF, while the CMF1, CMF2, and CMF3 groups were supplemented with different CMF at 6,000 mg/kg diet, respectively. The broilers were raised to 60 d of age, then bone tissues were collected for biomechanical properties, micro-CT detection and transcriptomic sequencing analysis. The results showed that CMF3 improved the biomechanical properties of broiler tibia, via increasing the elastic modulus (P < 0.05), yield strength (P > 0.05), maximum stress (P < 0.05) and fracture stress (P < 0.05) of the tibia. Micro-CT analysis indicated that CMF3 increased the bone mineral density (BMD), bone volume/total volume (BV/TV), bone surface density (BS/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), and decreased the trabecular separation (Tb.Sp) of femur cancellous bone (P < 0.05). RNA-seq analysis revealed 2,177 differentially expressed genes (DEGs) (|log2FoldChange| ≥ 1, FDR < 0.05) between the CMF3 group and CON group. Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) analysis showed 13 pathways mostly associated with bone growth and development and bone metabolism, and we identified 39 bone-related DEGs. This study suggests that CMF3 could improve bone strength and bone microstructure of broilers, and showed a positive effect on bone performance. Our research could provide a theoretical reference for the development of pollution-free feed additives to improve the skeletal performance of broilers, which could help promote healthy farming of chickens.

Carcass composition, meat quality, leg muscle status, and its mRNA expression profile in broilers affected by valgus-varus deformity

Valgus-varus deformity (VVD) is a common leg disease in commercial broilers, which seriously affects animal welfare and causes economic losses. Up to now, most of the studies on VVD have been on skeleton, whereas there are fewer studies on VVD muscle. In this study, carcass composition and meat quality of 35-day-old normal and VVD Cobb broilers assess the effect of VVD on broiler growth. Molecular biology, morphology, and RNA sequencing (RNA-seq) were used to study the difference between normal and VVD gastrocnemius muscle. In comparison with the normal broilers, the breast muscle and leg muscle of the VVD broilers had lower shear force, notably lower crude protein, lower water content, cooking loss, and deeper meat color (P < 0.05). The morphological results showed that the weight of skeletal muscle was significantly higher in the normal broilers than that in the VVD broilers (P < 0.01), the diameter and area of myofibrils in the affected VVD were smaller than in the normal broilers (P < 0.01). Quantitative real-time PCR (qPCR) of gastrocnemius muscle revealed that the expression of myasthenic marker genes, fast myofiber marker genes, and apoptosis-related factors were significantly higher in the VVD broilers than in the normal broilers (P < 0.01). In total, 736 differentially expressed genes (DEGs) were identified firstly in the normal and VVD leg muscle by RNA-seq. Gene ontology (GO) enrichment indicated that these DEGs were mainly involved in the multicellular organismal process and anatomical structure development. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that DEGs are significantly enriched in proteasome. Protein interaction analysis obtained that DEGs with high interaction were proteasome-related coding genes and ubiquitin-related genes, these DEGs were closely associated with muscle atrophy. These show that VVD has an adverse effect on growth characteristics, slaughter characteristics, and meat quality in broilers, which may cause leg muscle atrophy. This study provides some reference values and basis for studying the pathogenesis of VVD in broilers.

Valgus-varus deformity induced abnormal tissue metabolism, inflammatory damage and apoptosis in broilers

1. This study explored the tissue metabolic status and the relationship with inflammation in valgus-valgus deformity (VVD) broilers with increasing age.2. Tissue and blood from VVD and healthy broilers were collected at two, four and five weeks old. A fully automated biochemical analyser, real-time PCR, HE staining and enzyme-linked immunosorbent assay were used to detect tissue metabolic indexes, mRNA levels of inflammation and apoptosis cytokines in immune organs, histological changes and serum inflammation and immune-related protein contents in VVD broilers.3. The results showed that VVD increased the levels of total protein, albumin, alanine aminotransferase at five weeks of age, aspartate aminotransferase, urea and creatine kinase in blood at two weeks of age. It upregulated the gene expression of inflammatory factors IL-1β, IL-6, IL-8, TNF-α, NF-κB and TGF-β and apoptotic factors FAS, Bcl-2, caspase-3 and 9 in immune organs; increased levels of serum proteins TNF-α, IL-1β and IL-6 and decreased levels of serum immunoglobulins IgY and CD3⁺.4. In addition, with increasing age, IL-10 gene expression gradually increased in the BF and decreased in the spleen.5. In conclusion, VVD broilers have disorders of liver and kidney metabolism, inflammation and apoptosis of immune organs and increased levels of serum inflammatory factor proteins.

Identification and genetic analysis of major gene ST3GAL4 related to serum alkaline phosphatase in chicken

Alkaline phosphatase (ALP) is a marker of osteoblast maturation and an important indicator of bone metabolism. The activity of ALP can reflect the bone metabolism and growth traits of animals, so the polymorphism affecting ALP expression deserves further study. In this study, we identified an SNP site in ST3GAL4 found by genome-wide association studies (GWAS) in previous studies, 8 SNPs were also identified by DNA sequencing. Interestingly, there were 4 SNPs (rs475471G > A, rs475533C > T, rs475621A > G, rs475647C > A) completely linked by haplotype analysis. Therefore, we selected a tag SNP rs475471G > A to further analyze the ALP level of different genotypes in Hubbard leg disease population and an F₂ chicken resource population produced by Anka and Gushi chickens and carried out population genetic analysis in 18 chicken breeds. Association analysis showed that this QTL within ST3GAL4 was highly correlated with ALP level. The mutant individuals with genotype AA had the highest ALP level, followed by GA and GG carriers. The mutant individual carriers of AA and GA genotype had higher values for body weight (BW), chest width (CW), body slanting length (BSL), pelvis width (PW) at 4-week, the semi-evisceration weight (SEW), evisceration weight (EW) and Leg weight (LW) than GG genotypes. The amplification and typing of 4852 DNA samples from 18 different breeds showed GG genotype mainly existed in egg-type chickens and dual-type chickens, while the AA genotype was mainly distributed in commercial broilers and F₂ resource population. The individual carriers of the AA genotype had the highest ALP and showed better growth performance. Besides, tissue expression analysis used Cobb broiler showed significant differences between different genotypes in the spleen and duodenum. Taken together, this was the first time to determine 9 SNPs within ST3GAL4 related to ALP in chickens, 4 of them were complete linkage with each other, which provides useful information on the mutation of ST3GAL4 and could predict the serum ALP level of chicken early and as an effective potential molecular breeding marker for chickens.

Genome-wide association study identifies SNPs for growth performance and serum indicators in Valgus-varus deformity broilers (Gallus gallus) using ddGBS sequencing

Background

Valgus-varus deformity (VVD) is a lateral or middle deviation of the tibiotarsus or tarsometatarsus, which is associated with compromised growth, worse bone quality and abnormal changes in serum indicators in broilers. To investigate the genetic basis of VVD, a genome wide association study (GWAS) was performed to identify candidate genes and pathways that are responsible for VVD leg disease, serum indicators and growth performance in broilers.

Results

In total, VVD phenotype, seven serum indicators and three growth traits were measured for 126 VVD broilers (case group) and 122 sound broilers (control group) based on a high throughput genome wide genotyping-by-sequencing (GBS) method. After quality control 233 samples (113 sound broilers and 120 VVD birds) and 256,599 single nucleotide polymorphisms (SNPs) markers were used for further analysis. As a result, a total of 5 SNPs were detected suggestively significantly associated with VVD and 70 candidate genes were identified that included or adjacent to these significant SNPs. In addition, 43 SNPs located on Chr24 (0.22 Mb - 1.79 Mb) were genome-wide significantly associated with serum alkaline phosphatase (ALP) and 38 candidate genes were identified. Functional enrichment analysis showed that these genes are involved in two Gene Ontology (GO) terms related to bone development (cartilage development and cartilage condensation) and two pathways related to skeletal development (Toll−like receptor signaling pathway and p53 signaling pathway). BARX2 (BARX homeobox 2) and Panx3 (Pannexin 3) related to skeleton diseases and bone quality were obtained according to functional analysis. According to the integration of GWAS with transcriptome analysis, HYLS1 (HYLS1 centriolar and ciliogenesis associated) was an important susceptibility gene.

Conclusions

The results provide some reference for understanding the relationship between metabolic mechanism of ALP and pathogenesis of VVD, which will provide a theoretical basis for disease-resistant breeding of chicken leg soundness.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08236-3.

Identification of key miRNAs affecting broilers with valgus-varus deformity by RNA sequencing and analysis of miRNA-mRNA interactions

Valgus-varus Deformity (VVD) leg disease often affects chickens raised in modern large-scale breeding operations. Losses due to VVD are costly to farmers, and the condition also causes significant suffering in affected birds. In this study, we profiled RNAs from the spleens of VVD (BS) and healthy (JS) broilers using high-throughput sequencing to identify miRNAs that might be involved in the development of the disease. Fifty differentially expressed miRNAs (DEMs) were found, of which 30 were up-regulated and 20 were down-regulated in VVD-affected birds (|log 2 Fold Change| ≥ 1 and q-value < 0.05). DEMs were matched with putative target genes and 864 target genes were found. Gene Ontology (GO) analyses of these target genes showed that they were significantly enriched in the "cytoplasm" term (q-value < 0.05), and most of the target genes were enriched in "cellular component". Kyoto encyclopedia of genes and genomes (KEGG) analysis showed that they were significantly enriched in 11 signaling pathways (P-value < 0.05), including metabolic pathways, 2-oxocarboxylic acid metabolism, regulation of actin cytoskeleton, purine metabolism, endocytosis and so on. And we found that they were enriched in immune-related pathways in which MAPK, Notch, JAK-Stat, Toll-like receptor, p53 and other single pathways were involved in the development of skeletal diseases. Differentially expressed mRNAs obtained from our previous study were used to construct an interaction network consisting of 16 DEMs and 21 differentially expressed mRNAs (|log 2 Fold Change| ≥ 1 and q-value < 0.05). We found that miR-12247-5p, miR-15c-5p, miR-15b-5p, FKBP5 and HSP90AB1 were at the center of network interaction. This study provides a foundation for further investigations of the pathogenesis and genetic mechanisms underlying VVD.