Next Article in Journal
Thermal and Circulatory Changes in Diverse Body Regions in Dogs and Cats Evaluated by Infrared Thermography
Next Article in Special Issue
G0S2 Gene Polymorphism and Its Relationship with Carcass Traits in Chicken
Previous Article in Journal
The Utility of Serological Analysis for Neospora caninum Infection in Dairy Cattle Farms Management: Serological Investigation and Evaluation of the Effects on Reproductive and Productive Performances in Two Study Herds in Northern Italy
Previous Article in Special Issue
Associations of Transcription Factor 21 Gene Polymorphisms with the Growth and Body Composition Traits in Broilers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 6
10.3390/ani12060788
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens’ Femoral Head Separation

by
Iara Goldoni
1,
Adriana Mércia Guaratini Ibelli
1,2,*,
Lana Teixeira Fernandes
2,
Jane de Oliveira Peixoto
1,2,
Ludmila Mudri Hul
1,
Maurício Egídio Cantão
2,
João José de Simoni Gouveia
3 and
Mônica Corrêa Ledur
2,4,*
1
Programa de Pós-Graduação em Ciências Veterinárias, Universidade Estadual do Centro-Oeste, R. Salvatore Renna, 875, Guarapuava 85015-430, PR, Brazil
2
Embrapa Suínos e Aves, Rodovia BR 153, km 110, Concórdia 89715-899, SC, Brazil
3
Programa de Pós-Graduação em Ciências Veterinárias no Semiárido, Universidade Federal do Vale do São Francisco, UNIVASF, Rodovia BR 407, 12 Lote 543, Petrolina 56300-000, PE, Brazil
4
Programa de Pós-Graduação em Zootecnia, UDESC-Oeste, Rua Beloni Trombeta Zanin, 680E, Chapecó 89815-630, SC, Brazil
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(6), 788; https://doi.org/10.3390/ani12060788
Submission received: 28 February 2022 / Revised: 15 March 2022 / Accepted: 17 March 2022 / Published: 20 March 2022
(This article belongs to the Special Issue Genetic Selection and Associated Defects in Poultry Production)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Femoral head necrosis (FHN) and other locomotor problems cause severe impacts on the poultry industry due to huge economic losses and reduced animal welfare. Femoral head separation (FHS), the initial phase of FHN, is usually a subclinical condition characterized by the detachment of articular cartilage from the bone. In this study, we aimed to identify genes and biological processes involved with FHS in broilers. A better understanding of the FHS molecular mechanisms can help to develop strategies to reduce this condition in chickens. Here, we described several genes that have their expression altered in the articular cartilage and femur when normal and FHS-affected animals were compared. Furthermore, genetic variants were found differing between the studied groups. Therefore, performing an integrated analysis of these datasets, we were able to detect genes and variants related to FHS in chickens. Some of them, such as SLC4A1, RHAG, ANK1, MKNK2, SPTB, ADA, C7 and EPB420 genes were highlighted and should be further explored to validate them as candidates to FHS and FHN in chickens and possibly in humans.

Abstract

Femoral head separation (FHS) is usually a subclinical condition characterized by the detachment of articular cartilage from the bone. In this study, a comprehensive analysis identifying shared and exclusive expression profiles, biological processes (BP) and variants related to FHS in the femoral articular cartilage and growth plate in chickens was performed through RNA sequencing analysis. Thirty-six differentially expressed (DE) genes were shared between femoral articular cartilage (AC) and growth plate (GP) tissues. Out of those, 23 genes were enriched in BP related to ion transport, translation factors and immune response. Seventy genes were DE exclusively in the AC and 288 in the GP. Among the BP of AC, the response against bacteria can be highlighted, and for the GP tissue, the processes related to chondrocyte differentiation and cartilage development stand out. When the chicken DE genes were compared to other datasets, eight genes (SLC4A1, RHAG, ANK1, MKNK2, SPTB, ADA, C7 and EPB420) were shared between chickens and humans. Furthermore, 89 variants, including missense in the SPATS2L, PRKAB1 and TRIM25 genes, were identified between groups. Therefore, those genes should be more explored to validate them as candidates to FHS/FHN in chickens and humans.

1. Introduction

The poultry industry has experienced huge growth in recent decades, due to the advances in new research and technologies, which has brought benefits to the industry and producers. Chicken meat is one of the most consumed proteins in the world, due to the low cost, for being considered cost-effective and healthy, and because it is exempted from cultural impediments [1]. The constant investment in genetic improvement, nutrition and management has generated a better performance of the broilers, with slaughter age reduction, better feed conversion and less water consumption, always ensuring the animal welfare [2].
The genetic breeding programs have brought numerous benefits for poultry production, such as the increased body weight, greater thoracic musculature, and reduction of abdominal fat. In the last 50 years, chicken growth rates have risen more than 300%, from 25 g per day to approximately 100 g per day [3,4]. However, some negative consequences, such as locomotor problems and immunological and metabolic disorders have appeared due to the selection for fast-growing [3,5].
Locomotor problems cause severe impacts on the poultry industry, affecting about 6% of the animals in commercial flocks, which results in huge economic losses [6]. Furthermore, animals with bone disorders have their welfare affected because they cannot eat and drink water properly, being susceptible to opportunistic agents [7]. The bacterial chondronecrosis with osteomyelitis (BCO), also known as femoral head necrosis (FHN), is the most prevalent leg disorder in broilers [8]. The real incidence of FHN in broilers is difficult to estimate since most of the lesions often remain subclinical, and its etiological basis is not fully understood [9]. However, some authors have shown 20% of FHN incidence in chickens [10]. FHN occurs in the proximal part of the chickens’ femur head, starting with the femoral head separation (FHS), which consists in the separation of the growth plate (GP) from the articular cartilage (AC) [9,11]. This occurs due to the high growth rates and the mineralization deficiency of the chondrocytes that cause damage to the bone, which could allow for the colonization by opportunistic bacteria featuring the BCO [12].
Since the etiology of FHN remains unclear, studies are being conducted to clarify its pathogenesis and the molecular mechanisms involved, and to discover genes with a potential role in triggering FHS and subsequent BCO in the head of the femur and tibia [10,13,14,15,16,17]. Although most of the previous studies describe biological processes (BP) involved with this condition in bone tissue, there are no studies evaluating gene expression profiles between bone and articular cartilage of normal and FHS-affected chickens. The femur GP or physis is a point of huge stress during the broilers growing phase. The GP tissue presents a cartilaginous matrix and sequential layers of long columns of chondrocytes in different maturation stages [12]. Above the GP is the AC that consists of chondrocytes and extracellular matrix (ECM), which is mainly composed by proteoglycan, glycosaminoglycan, and collagen fibers [18]. The understanding on how these tissues are interacting might be one way to improve the knowledge of molecular mechanisms involved with FHS. Therefore, in this study, the femoral growth plate and articular cartilage transcriptomes from normal and FHS-affected broilers were compared to identify common and exclusive molecular mechanisms and differentially expressed (DE) genes between groups, and polymorphisms highlighting possible new candidate genes involved with FHS in chickens.

2. Material and Methods

Transcriptomes used in the present study were obtained from a project that evaluated femoral head separation (FHS) in the growth plate and articular cartilage, published respectively by Peixoto et al. [15] and Hul et al. [17], which are available in the SRA database BioProject # PRJNA350521 for articular cartilage (AC) and PRJNA352962 for growth plate (GP). Briefly, 16 AC and GP samples from 35-day-old male broilers (Cobb500), previously characterized as normal (4 AC and 4 GP) or FHS-affected (4 AC and 4 GP), were used in this study. At the time of collection, the normal samples had a good adhesion between the AC and GP, and the FHS-affected samples had the AC easily detached from the GP. For the normal group, the AC was separated from the GP using a scalpel. These tissues were submitted to RNA extraction using Trizol and a cleanup using RNeasy mini kit (Qiagen, Hilden, Germany). RNA quality was assessed in 1% agarose gel and in Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Samples with an RNA integrity number (RIN) higher than 8 were used to construct the RNA-Seq libraries. Eight samples from each tissue (4 normal and 4 FHS-affected), were prepared for RNA sequencing. For GP, the libraries were prepared from 2 µg of RNA using the “TruSeq RNA Sample Prep Kit v2” (Illumina, San Diego, CA, USA), while for AC the “TruSeq Total Stranded Sample Preparation” (Illumina, San Diego, CA, USA) kit was used. Libraries were sent to the Functional Genomics Center, ESALQ, University of São Paulo, Piracicaba, São Paulo State, Brazil, for sequencing in Illumina HiSeq2500 equipment (2 × 100 bp, Illumina, San Diego, CA, USA), with samples of each tissue in a different lane.

2.1. Bioinformatics, Differential Expression Analyses and Functional Annotation

The baqcom pipeline (https://github.com/hanielcedraz/BAQCOM, accessed on 31 December 2021) was used for quality control (QC) and mapping. QC for each sequencing file was performed with the Trimmomatic v. 0.38 [19] to remove short reads (<70 bp), low-quality reads (QPhred < 20) and adapter sequences. AC and GP sequences were mapped against the chicken reference genome (GRCg6a, Ensembl 101) using the STAR software v 2.7 [20], and reads were counted with HTSeq-counts [21]. The EdgeR package [22] was used to identify the DE genes between normal and affected groups in each tissue separately, considering differentially expressed those with false discovery rate (FDR) ≤ 0.05 after correcting for the Benjamini–Hochberg (BH) multiple tests [23]. Positive and negative log2 fold-change corresponded to DE genes upregulated or downregulated, respectively, in FHS-affected compared to the normal group. A heatmap was generated to visualize patterns of expression across samples from both tissues using R. The expressed genes were annotated using the Biomart database (https://www.ensembl.org/biomart, accessed on 31 December 2021). Visualization and Integrated Discovery (DAVID 6.8) [24] and Panther (http://pantherdb.org/, accessed on 31 December 2021) [25] databases were used for assessing functional profiles of DE genes based on the biological processes (BP), cellular components (CC), molecular functions (MF) and metabolic pathways categories of gene ontology (GO). The enrichment analyses were performed using the chicken genome information available in those databases.

2.2. qPCR Validation

A relative quantification using qPCR was performed to confirm the expression profile of eight chosen DE genes in the GP using eight samples from the four normal and four from the FHS-affected group. The genes selected were solute carrier family 4 member 1 (SLC4A1), Rh- associated glycoprotein (RHAG), ankyrin 1 (ANK1), MAP kinase interacting serine/threonine kinase 2 (MKNK2), spectrin beta, erythrocytic (SPTB), adenosine deaminase (ADA), complement C7 (C7) and erythrocyte membrane protein band 4.2 (EPB42). Gene sequences were downloaded from Gallus gallus on Genebank (http://www.ncbi.nlm.nih.gov/gene/, accessed on 31 January 2022) and Ensembl (www.ensembl.org, accessed on 31 January 2022). Primers were designed in exon–exon junctions using the Primer-Blast online tool [26], and their quality was evaluated in the NetPrimer online software (http://www.premierbiosoft.com/netprimer/, accessed on 31 January 2022) (Table 1). For AC, the analysis was also performed in eight samples for five genes (SPTB, MKNK2, SLC4A1, C7 and EPB42), since the RHAG, ANK1 and ADA genes have already been confirmed in the analyzed population by Hul et al. [17]. Ribosomal protein 4 (RPL4) and Ribosomal protein 30 (RPL30) genes were used as reference for GP, and Ribosomal protein 5 (RPL5) and Ribosomal Protein Lateral Stalk Subunit P1 (RPLP1) for AC, as described by Peixoto et al. [15] and Hul et al. [17], respectively. The cDNA was synthetized using SuperScript III First-Strand Synthesis SuperMix (Invitrogen, Waltham (MA), USA), and qPCR reactions were carried out in QuantStudio 6 (Applied Biosystems, Waltham (MA), USA) equipment. The reactions had a final volume of 15 μL with 1×GoTaq qPCR Master Mix with BRYT Green (Promega, Madison (WI), USA), 0.13 µM of each primer and 2 μL of cDNA. Cycle threshold (Ct) mean for each replicate sample was obtained and normalized using the geometric mean of the reference genes, calculating 2−ΔCt [27]. Statistical analyses for group comparisons were performed using Mann–Whitney–Wilcoxon test in the R environment, considering DE as the genes with p ≤ 0.05.

2.3. Integrated Analysis

To identify the main molecular mechanisms involved in FHS, the common and exclusive DE genes in the AC and GP datasets were obtained, and the David [24] database was used to find BP involved with FHS based on the gene ontology database. The REVIGO tool [28] was used to reduce and highlight the most abundant BP in our data. Furthermore, gene networks with the common DE genes between the two tissues and those obtained with DE genes only detected in AC or GP were constructed using the STRING database [29].
Furthermore, we used the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds, accessed on 9 November 2020) to search datasets with the terms “femoral head necrosis” and “epiphysiolysis” to better understand the role of the main DE genes involved with FHS across different species. This search returned two main datasets: identification of potential biomarkers for improving the precision of early detection of steroid-induced osteonecrosis of the femoral head (GSE123568) and gene expression profile of hip cartilage with necrosis of femoral head (GSE74089). The first one evaluated the gene expression on blood from 40 patients (10 normal and 30 affected with femoral head necrosis) and the second one, the hip cartilage in 24 patients (12 normal and 12 affected with femoral head necrosis). DE genes were obtained using the GEOR tool from the GEO dataset (https://www.ncbi.nlm.nih.gov/geo/info/geo2r.html, accessed on 9 November 2020). The list of DE genes from the two datasets obtained in the GEO database and the DE genes found in the GP and AC from our study were submitted to the InteractiVenn (http://www.interactivenn.net/, accessed on 10 November 2020) [30] to find common DE genes among the experiments.

2.4. Polymorphism Identification Using the RNA-Seq Data

The Genome Analysis Tool kit 3.6 (GATK) [31] was used for polymorphism identification. We followed the best practices guidelines standard parameters for transcriptome variant analysis available on the GATK website (https://gatk.broadinstitute.org/, accessed on 1 February 2021). The Picard tools 2.5 (https://broadinstitute.github.io/picard/index.html, accessed on 17 May 2017) was used to generate the genome index, to assign read groups and mark duplicates. In the GATK, the CIGAR strings determination (SplitNCigarReads), qualities reassigning mapping, base recalibration, variants calling, and filtering were performed. To minimize the false positive variants, the following filters were used in GATK to filter and select variants: FS > 30.0, MQRankSum < −12.5, SNPcluster considering 3 variants in a 35 bp window, QD < 5.0, MQ < 50.0, GQ < 5.0, QUAL ≥ 30.0, ReadPosRankSum < −8.0 and DP ≥ 100.0. Furthermore, variants with less than 10 reads per sample were removed from the downstream analysis. Once the polymorphisms were identified in the GP and AC datasets, a filter was performed to obtain the variants that differed between normal and FHS-affected groups. The Variant Effect Predictor (VEP) tool [32] was used for variant annotation, effect and consequence predictions using the Gallus gallus genome (GRCg6a) with the Ensembl 103 annotation database version. Furthermore, the Stringdb [29] and EnrichR [33,34] databases were used to verify the interactions between DE genes and variants identified in genes.

3. Results

3.1. RNA-Sequencing, Mapping and Characterization

The RNA sequencing generated about 345 million paired end reads for all samples (N = 16), and after the quality control, approximately 306 million reads (88.78%) remained for further analyses (Supplementary File S1: Table S1), with an average of 16,355,490 reads for GP and 26,819,697 for AC. A clear separation between the samples from the two tissues, as well as between the normal and affected groups can be verified in the heatmap (Supplementary File S2: Figure S1). As expected, most of the genes were characterized as protein coding (95%), and the remaining 5% were classified in IncRNA, pseudogenes, miRNAs and other noncoding RNAs (Table 2).

3.2. Differentially Expressed Genes in AC and GP Transcriptomes

Using the Ensembl annotation 101, in the AC transcriptome, 12,182 genes were identified and 106 genes were DE; as well, 99 (93.44%) were upregulated and seven (6.56%) were downregulated in the AC of the FHS-affected compared to the normal group (Supplementary File S1: Table S2). Regarding the GP tissue, 12,632 genes were expressed, where 324 were DE between the analyzed groups. From those, 174 were upregulated (53.7%) and 150 (46.3%) downregulated in the affected compared to the normal group (Supplementary File S1: Table S3).
Comparing the DE genes between the AC and the GP, 36 DE genes were common to both transcriptomes (Table 3). From these 36, 33 were upregulated and three were downregulated in the FHS-affected group in both transcriptomes. From the 70 DE genes remaining in the AC data, four were downregulated and 66 were upregulated in the FHS-affected broilers (Table S2). Finally, 141 genes were upregulated and 147 were downregulated in the affected group when considering the DE genes exclusive to the GP tissue (Table S3).
The enrichment analysis performed based on the common DE genes between the two tissues showed that the main biological processes found were related to ion transport, inflammatory and defense response, homeostasis, protein biogenesis and immune response processes (Figure 1).
The most relevant BP for the AC tissue was the defense response to bacterium, cell chemotaxis and regulation of cell motility (Table S4). Conversely, the regulation of lipid metabolic process, cytokine-mediated signaling pathway, T cell activation, inflammatory response, interleukin-1 beta production, small GTPase-mediated signal transduction, and negative regulation of chondrocyte differentiation and cartilage development were the BP most prominent on GP tissues (Table S5).
A gene network was constructed based on the 36 common DE genes. From those, 23 were recognized by Stringdb tool, where it was possible to observe a main network of grouped genes related to ion transport and translation factors (Figure 2). The transferrin (TF) gene is involved in antimicrobial activity, while some genes such as Hemoglobin subunit gamma-2 (HBG2) and Rh-associated glycoprotein (RHAG) are involved in oxygen and ion transport. Another group connected the PTPRC, CCL10, ISG12-2 and SAMD9L genes that are most related to the immune response BP (Figure 2).
Another two gene networks were constructed based on the DE genes exclusive for AC or for GP. For AC, 61 genes were recognized by the Stringdb tool, which highlighted a group of genes related to the response against bacteria, with all genes involved in this BP being more expressed in the FHS-affected group (Figure 3). The second network was constructed with 266 GP DE genes, and the processes related to chondrocyte differentiation and cartilage development stood out, grouping genes that had the lowest expression in the affected group (Figure 4).

3.3. qPCR Validation

In the qPCR analysis, from the eight candidate genes evaluated in GP, four (C7, EPB42, MKNK2 and SPTB) were DE. The MKNK2 and SPTB genes were also DE in the AC. Furthermore, all evaluated genes were upregulated in the FHS-affected group (Figure 5A,B), having the same expression profile and therefore confirming the RNA-Seq results.

3.4. Comparison of the DE Genes from AC and GP Tissues with Different Datasets

After the differentially expressed analysis, we compared the DE genes from the AC and GP tissues with other datasets obtained in the GEO database, as previously described. According to the GEOR analysis, the comparison performed between normal and FHN-affected individuals (GSE123568) evinced 5250 DE genes in the blood of those individuals, while 5978 DE genes were found in the hip cartilage of normal and FHN-affected individuals (GSE74089). Considering these four datasets, four genes were DE in all studies: SLC4A1 (Solute Carrier Family 4 Member 1), EPB42 (Erythrocyte Membrane Protein Band 4.2), SPTB (Spectrin Beta, Erythrocytic) and ANK1 (Ankyrin 1) (Figure 6A), all of them related to erythrocyte and ion transport functions.
We also performed a comparison using DE genes from cartilage (chicken and human) and bone datasets to verify the presence of shared genes in the tissues affected by FHN/FHS. In this comparison, eight shared genes were found: the four previously described (SLC4A1, EPB42, ANK1 and SPTB) and four additional genes: CCL26 (Chemokine (C-C motif) ligand 26), ADA (Adenosine Deaminase), SLC25A37 (Mitoferrin-1) and MKNK2 (MAPK Interacting Serine/Threonine Kinase 2), all involved with immune response (Figure 6B).

3.5. Variant Identification and Annotation

A total of 89,823 variants (SNPs and InDels) were found in the analyzed samples with the GATK tool. Out of those, 89 differed between the normal and FHS-affected groups (Supplementary File S1: Table S6). These 89 variants were processed in the VEP tool and approximately 3% from all of them (three) were first described in the current study, while 86 (~97%) have already been identified in chickens. One of the three novel SNPs was located in the 5′ UTR region of the NET1 (Neuroepithelial Cell Transforming 1) gene, another SNP overlapped the downstream region of two genes, CCM2 Scaffold Protein (CCM2) and NAC Alpha Domain Containing (NACAD, ENSGALG00000050931—predicted), and the other was a missense variant in the Staphylococcal Nuclease and Tudor Domain Containing 1 (SND1, ENSGALG00000054297—predicted). Variants were mainly classified as downstream gene variant, synonymous variant, upstream gene variant, intron variant, 3′ prime UTR variant and missense variants (Figure 7).
Out of the described coding variants, 90% was classified in synonymous variants, 10% was missense (Supplementary File S1: Table S7). Considering the four missense variants, three were previously described in the genes PRKAB1 (rs733851944), SPATS2L (rs737569658) and TRIM25 (rs314175171), with a Sorting Intolerant from Tolerant (SIFT) score of 0.36 and 0.16. The novel variant (chr1:1049560) in the SND1 had a SIFT = 0, being predicted as deleterious.
After identifying the variants that differ between the normal and affected groups in the RNA sequences, a search was performed to verify if some of the variants were present in the DE genes in both AC and GP tissues previously found between normal and FHS-affected samples. None of the variants were found in the common DE genes; however, two polymorphisms were found in the DE genes in the GP: a SNP (rs737937191) in the downstream region of the LOXL2 gene and the other in the amino acid position 52 of the SND1 (ENSGALG00000054297) gene.
Some SNPs that differed between normal and affected broilers were found in genes that have already been described in genomic regions associated with related phenotypes in the chicken QTLdb (https://www.animalgenome.org/cgi-bin/QTLdb/GG/index, accessed on 31 December 2021), such as a SNP in the RPS24 gene, located in a QTL region associated with tibia volume (135,882, Animal QTLdb) (Supplementary File S1: Table S8). Furthermore, some of the genes that have variants differing between the normal and affected groups, such as Solute Carrier Family 25 Member 3 (SLC25A3), Yip1 Domain Family Member 3 (YIPF3), DEAD-Box Helicase 55 (DDX55), Exportin 5 (XPO5), RNA Polymerase I And III Subunit C (POLR1C), R3H Domain And Coiled-Coil Containing 1 (R3HCC1), Transmembrane P24 Trafficking Protein 2 (TMED2), Ras Responsive Element Binding Protein 1 (RREB1), charged Multivesicular Body Protein 7 (CHMP7), DNA Polymerase Eta (POLH) and Lysyl oxidase homolog 2 (LOXL2) have already been associated with waist-to-hip ratio adjusted for body mass index (BMI) in humans (Supplementary File S1: Table S9).
The interactions between the DE genes and genes with variants were verified through a gene network, where it was possible to observe that those genes have physical or functional interactions and were grouped in five main subnetworks (Figure 8). For example, the ADA gene was DE in AC and GP and grouped with several genes, which have SNPs differing between normal and affected groups, such as the TRIM25, LUM, MYH9, RPS24 and NACAD.

4. Discussion

The mechanical stress caused by the overload of the skeleton that modern chickens sustain due to the high growth rates has been recognized as one of the main causes of locomotor problems [11,12,35]. Among these locomotor problems, the FHN/BCO is a disorder that has a large impact on the poultry production due to the huge economic losses and its negative effect on animal welfare [8,9,15]. One of the biggest issues to study FHN is its early detection, since the animals do not show clinical signs often, being visible only after slaughter [16].
There are few studies related to FHN and other bone integrity problems in chickens, especially those approaching the molecular mechanisms involved with these conditions. Previous studies from our group have shown some biological processes (BP) and genes associated with FHN/BCO in different ages and lines of chickens [14,15,16,17]. Paludo et al. [13] pointed out that the downregulated expression of RUNX2 and SPARC genes may be associated with reduced vascularization and poor bone mineralization, increasing the risk of skeletal problems in chickens. The impaired collagen formation, connective tissue cell adhesion, and bone extracellular matrix (ECM) were also identified as predisposing factors for BCO and other leg disorders [14,15,16,17]. The translocation of bacteria from blood to bone has also been suggested as a secondary condition [16], whereas the bacterial profile of the microbiome in the blood of animals affected with BCO was different from healthy ones [36]. However, most of these studies focus only on the bone tissue, and there are no studies evaluating the articular cartilage and the femoral head jointly. Therefore, our study aimed to find common DE genes and BP between those tissues as well as genes exclusively DE in each tissue (AC and GP), and some of these genes were validated by qPCR (Figure 5). We also identified variants in the AC and GP transcriptomes between normal and FHS-affected broilers through an integrated analysis of the RNA sequencing datasets previously published by Peixoto et al. [15] and Hul et al. [17].
In the global transcriptome characterization of the GP and AC tissues, about 12,000 genes were expressed in both tissues. From those, 11,141 were shared and 891 and 441 genes were exclusively to the GP and AC, respectively. As expected, due to the RNA sequencing methodology used, approximately 95% of the genes were coding, although lncRNAs, miRNAS, mitochondrial and miscRNAs were also identified in the analyzed datasets (Table 2).
A total of 36 DE genes were shared by both tissues, and 91.6% of the common genes were upregulated in the FHS-affected group (Table 3). Some of those genes have already been highlighted in other studies [14,16,17]. However, when the common BPs were evaluated, most of them were related to ion transport, defense response and those related with tissue homeostasis (Figure 1).
Two main branches were observed in the common DE genes network (Figure 2): one was composed of the TF, HBG2, ANK1, SLC25A37, EPB42, RHAG, ADD2, KEL, MKNK2 genes and another by the PTPRC, CCLI10, ISG12-2 AND SAMD9L genes. In general, the first group of genes are related to blood type and to cartilage regeneration, and the second is related to inflammatory response [15,37]. Moreover, a set of variants was identified differing between normal and affected groups, where none of them were located in the common DE genes. It was observed that some of these polymorphisms were in genes that are co-expressed or co-regulated with DE genes (Figure 7). The possible role of some of these genes in the FHS are discussed below.
The transferrin (TF) was in the main branch of the gene network and its primary function is to transport iron to tissues that require this mineral [38], being involved in defense against systemic infection [38]. There are two forms of transferrin in birds: ovotransferrin, found in oviducts, and serum transferrin, secreted by the liver. Serum transferrin may stimulate cell proliferation and is regulated by iron levels, while ovotransferrin has a bacteriostatic role independent of iron levels [39,40]. Furthermore, the TF has a function in endochondral ossification, being considered the angiogenic molecule released by the hypertrophic cartilage [41]. TF presented strong interactions with the HBG2 (Hemoglobin Subunit Gamma 2) gene, which is involved in the transport of oxygen from the lungs to the peripheral tissues, where the beta chain is a component of adult hemoglobin A and D [8]. HBG2, as the TF, was upregulated in the FHS-affected broilers (Table 3). The gamma globin genes (HBG1 and HBG2) are expressed in several tissues, such as fetal liver, spleen and bone marrow [42]. Chicken has multiple types of globins, and erythropoiesis occurs in two waves: one primitive that acts in blood cells during early embryonic development, and the other producing definitive erythrocytes in late embryonic and post-hatching development [43]. In the bone and cartilage, globin function is not clear, and no information is available on its involvement with FHS in chickens.
The SLC25A37 (Solute Carrier Family 25 Member 37) works as a solute transporter located in the internal mitochondrial membrane, importing iron, which is essential for the synthesis of mitochondrial heme and iron–sulfur clusters [44]. In general, solute carrier proteins are essential to transport a wide set of molecules, such as ions, amino acids and vitamins to the tissues. Recently, it has been demonstrated that they are sensitive to extracellular levels of phosphate [45] and also regulate the pH of the cells [46]. In humans, several bone morphogenic proteins such as BMP13, BMP14 and BMP15 were responsible for regulating genes of the SLC family [47].
The ANK1 gene encodes a binding protein of the cytoskeleton, helping to bind other membrane proteins to the actin–spectrin cytoskeleton [48]. It also acts on contact activation, maintenance and proliferation of specialized membrane domains [48]. Multiple ankyrin isoforms with different affinities for various target proteins are expressed in a tissue-specific manner, regulated by development. The ANK1 is usually found in erythrocytes, but it has already been found in the brain and muscles [49]. In our study, the high expression of ANK1 can possibly affect the actin structure of the cytoskeleton, changing the structural integrity of the femoral articular cartilage and contributing to the occurrence of FHS. In addition, its expression is related to cell damage, which can be a consequence of FHS, because when FHS starts, positive regulation can act as a sign of attempt to combat the progression of this condition. This gene is also involved in the upregulation of inflammatory cytokines in osteoarthritic lesions [50]. The ANK1 gene is regulated by the EPB42, another gene that was upregulated in the FHS-affected group. The ANK1, EPB42 and RHAG (Figure 2) are closely related to some genes that have SNPs differing between the two analyzed groups, such as TMED2, DDX55 (Figure 7) and POLRC1.
The MKNK2 (Serine/threonine-protein kinase 2) gene interacts with MAP kinase and may respond to environmental stress and cytokines. Among its related pathways are the interleukin 1 signaling pathway. MAP/ERK signaling acts as a mediator of the suppressive effects of IFN-gamma in hematopoiesis [51]. The MKNK2 is not normally expressed in rats, but appears in stressful situations, showing its adaptive function and can also be a signal for cellular apoptosis [52]. Bringing it to our study, this gene may be involved in the initial phase of FHN with the function of initial signaling.
Another DE gene, the MYH15, has also been associated with the adrenergic signaling process in cardiomyocytes, in cases of imbalance in the O2 supply and removal of CO2 from tissues [53]. This gene could be associated with FHS as a consequence of local inflammation.
The EPX gene is released during the immune response by eosinophils, having a cytotoxic effect on cells [54,55] and bactericidal activity [56]. The eosinophil peroxidase, an extremely cytotoxic molecule that has anti-inflammatory and pro-inflammatory properties, regulates inflammation by combat invading microorganisms [57,58]. In chickens, the EPX gene has been studied as a biomarker for inflammatory events in the gastrointestinal tract [59] and has been related to leg disorders [15,17]. The EPX upregulation in broilers affected with FHS may indicate the immune system response to inflammation and, in severe cases, may be related to local necrosis.
The extracellular matrix (ECM) is a stable structural component that is located under the epithelium and close to connective tissue [60], and is responsible in providing support to the tissues and organs throughout the body. It acts in biochemical processes of the body assisting in the differentiation, morphogenesis and homeostasis of tissues [61]. The upregulated genes ADA and RHAG also participate in the process of metabolic glycosaminoglycans (GAGs) and aminoglycans involved in the metabolism of ECM [62], indicating that the body tries to repair the damage caused by FHS through tissue remodeling. Furthermore, the ADA enzyme acts as an endogenous regulator of the adaptive immune system, especially on proliferation and differentiation of T lymphocytes, regulating cell metabolism and triggering several physiological effects on cell proliferation [63,64].
Inflammation is an essential component of the immune system; however, an excess of inflammation can cause tissue damage [65]. The ADA gene also acts as a sensor providing information to the immune system regarding tissue damage, protecting cells from excessive tissue damage associated with inflammation [66]. The upregulation of ADA has a regulatory role in immune responses, acting on the activation and regulation of lymphocyte and neutrophil levels [64,67]. In our study, the ADA was upregulated in AC and GP in the FHS-affected group (Table 3), possibly due to the large tissue inflammation that occurs in the affected animals. The presence of specific variants in this gene between the studied groups also highlights that those genetic mechanisms could be involved with FHS in chickens. The ADA gene is in a branch of the gene network with several genes that have SNPs with different genotypes between the normal and affected groups, such as TRIM25, RPS24, NACAD, and others (Figure 7). The TRIM25 and ADA are important genes to the RNA machinery [66,68], and the involvement of the ADA in bone metabolism has been observed in humans [69]. Deficiency in ADA expression has led to a reduction of bone volume and downregulation of RANKL expression, while ADA upregulation can lead to skeletal abnormalities such as scapular spurring [70,71].
The ADA and IFI6 genes also play an important role in the regulation of apoptosis, which is an essential physiological mechanism in the development and tissue homeostasis [14,72]. The IFI6 gene, also known as ISG12, regulates cellular metabolism during the differentiation of osteoblasts and apoptosis [14,73]. The upregulation of the IFI6 gene may be related to a causative factor, stimulating apoptosis in the articular cartilage, leading the animal to be more susceptible to FHN.
Furthermore, once the exclusively DE genes in the GP were evaluated, the main biological processes were those related to chondrocyte differentiation and cartilage development, grouping three main genes that had the lowest expression in the affected group: CHADL, GDF5 and PTHLH (Figure 6, Table S5). Chondroadherin-like (CHADL) is a member of a family of collagen-associated small leucine-rich proteins (SLRPs) responsible for signaling of differentially regulated collagen fibrils during development, homeostasis, or pathogenesis. CHADL is expressed in cartilaginous tissues, influences collagen fibrillogenesis and negatively modulates chondrocyte differentiation [74].
In humans, increased GDF5 gene expression is related to joint tissue remodeling. After joint damage, its high expression in chondrocytes is observed, both in the new cartilage recovery tissue and in the adjacent damaged cartilage [75]. The PTHLH gene promotes the proliferation of chondrocytes and inhibits their hypertrophy and terminal differentiation. In patients with osteoarthritis, this gene is highly expressed in chondrocytes and synovial fluid [76].
Among the BP of the AC, the response against bacteria can be highlighted, with all genes involved in this BP being more expressed in the FHS-affected group (Figure 1). Among them, cathelicidin-2 (CAMP) and cathelicidin-3 (CATH3) stand out, which have a broad spectrum of antimicrobial activity and the ability to limit inflammation by inhibiting the activation of Toll-like receptor 2 (TLR2) and TLR4 [77].
These results suggest that the growth plate is not fully efficient in the process of recovery and remodeling of damaged cartilage. Conversely, articular cartilage can partially control bacterial infection and the intensity of inflammation that can intensify tissue damage. Other mechanisms should be evoked.
In addition to the identification of the DE genes, biological processes and SNPs found in our study, another interesting result presented here was the comparison of the chicken transcriptomes with other datasets. When we compared the DE genes obtained in the AC and GP tissues with other gene expression datasets obtained in the GEO, we found four genes in common, SLC4A1 (Solute Carrier Family 4 Member 1), EPB42 (Erythrocyte Membrane Protein Band 4.2), SPTB (Spectrin Beta. Erythrocytic) and ANK1 (Ankyrin 1) (Figure 6A), all of them related to erythrocyte and ion transport functions. In the second comparison, using the DE genes of bone and cartilage datasets, besides the four previously described genes (SLC4A1, EPB42, ANK1 and SPTB), four additional genes were found (CCL26, ADA, SLC25A37 and MKNK2), all involved with immune response. Furthermore, two genes, the CCL26 and SLC4A1, have not been enriched in the gene network. CCL26 (chemokine (C-C motif) ligand 26) is from the cytokine family of secreted proteins involved in immunoregulatory and inflammatory processes, displaying chemotactic activity for normal peripheral blood eosinophils and basophils, and it may contribute to the eosinophil accumulation in atopic diseases. Finally, the SLC4A1 gene encodes a protein that is part of the anion exchanger family and is expressed in the plasma membrane of erythrocytes, involved in the transport of carbon dioxide from the tissues to the lungs. Furthermore, this gene is located in QTL regions previously described as associated with femur mineral content and femur weight, which reinforces its importance as a good candidate gene to FHS.
Based on the functions and the high expression of the genes that were DE in both tissues, we had two hypotheses: the first is that the high expression of these genes may be a consequence of the FHS, due to cell damage and from the attempt to combat the progression of this condition. The second hypothesis is that the high expression and possible mutations of these genes cause FHN, as they can alter the structure of the actin cytoskeleton, affecting the structural integrity of the femoral articular cartilage, contributing to the occurrence of FHS. Therefore, using the two datasets in this study, we could highlight common genes that were regulated in the two tissues, as well as those that were exclusively DE in each one. With this information, new genes involved with FHS in broilers were described, where most of them have never been cited as candidates to this condition in chickens. We have also shown conserved mechanisms among huma and chicken osteonecrosis. Furthermore, genetic variants were identified; some of them were described for the first time in this study, and some of them differed between the normal and FHS-affected groups and could be further investigated as molecular markers for this condition. Finally, our results can also contribute to understand the onset of this condition in other species since the molecular mechanisms seem to be at least in part conserved among them.

5. Conclusions

With the integrated analysis of gene expression of the two tissues and variant identification, we were able to detect strong candidate genes and SNPs related to FHS/FHN in chickens. Comparing the AC and GP DE genes with the other datasets related to FHN allowed us to identify genes possibly involved with FHN, evincing the shared mechanisms across different species. Finally, the role of SLC4A1, RHAG, ANK1, MKNK2, SPTB, ADA, C7 and EPB420 genes should be more explored in order to validate them as candidates for FHS/FHN.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12060788/s1, Supplementary File S1: Table S1. Input paired read numbers, reads number after quality control and mapping information of each sample. A table containing the number of raw reads and after quality control and mapping information. Table S2. Expressed and DE genes between normal and FHS-affected broilers with 35 days of age in the AC. An excel sheet with the expressed and DE genes between normal and FHS-affected broilers with 35 days of age. Table S3. Expressed and DE genes between normal and FHS-affected broilers with 35 days of age in the GP. An excel sheet with the expressed and DE genes between normal and FHS-affected broilers with 35 days of age. Table S4. Main biological processes involved with FHS considering the DE genes exclusive. Table S5. Main biological processes involved with FHS considering the DE genes exclusive in the growth plate using the DAVID database. Table S6. Variants identified between normal and FHS-affected broilers with 35 days of age in the AC and GP transcriptomes according to the location in the Gallus gallus reference genome (GRCg6). Table S7. Variants annotated by the VEP tool (Ensembl 103). Table S8. Variants annotated in QTL regions using the VEP tool (Ensembl 103). Table S9. Candidate genes in human QTLs according to EnrichR. The genes identified in the common differentially expressed and variant analysis were used to verify their presence in human GWAS. Supplementary File S2: Figure S1. Heatmap gene cluster classification for AC and GP normal and FHS-affected samples. In the heatmap, the expression for each gene is presented in the rows, and samples are visualized in the columns, showing a hierarchical clustering of genes and samples. In red, DE genes were downregulated, and in green, upregulated in the affected samples.

Author Contributions

I.G., formal analysis, investigation, methodology, visualization, writing—original draft, writing—review and editing; A.M.G.I., conceptualization, data curation, formal analysis, investigation, methodology, project administration, supervision, writing—original draft, writing—review and editing; L.T.F., methodology, writing—original draft, writing—review and editing; J.d.O.P., conceptualization, data curation, formal analysis, investigation, methodology, project administration, supervision, funding acquisition, writing—original draft, writing—review and editing; J.J.d.S.G., methodology, formal analysis, writing—review and editing; M.E.C., conceptualization, formal analysis, methodology, writing—review and editing; L.M.H., formal analysis, investigation, visualization, writing—review and editing; M.C.L., conceptualization, data curation, investigation, methodology, validation, visualization, project administration, writing—original draft, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Brazilian Agricultural Research Corporation (Embrapa), project number 01.11.07.002.04.03, from the Brazilian Government.

Institutional Review Board Statement

All animal procedures were performed in accordance with the Ethics Committee on Animal Utilization of the Embrapa Swine and Poultry National Research Center under protocol number 012/2012.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used or analyzed during the current study are available from the corresponding author on reasonable request. The transcriptome sequences are available in the SRA database with BioProject number PRJNA352962 and PRJNA350521.

Acknowledgments

We thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil—Finance Code 001 for the free access to the journals used in the literature review. L.T.F. is supported by a visiting specialist fellowship from the National Council of Scientific and Technological Development (CNPq) grant #380282/2021-6. M.C.L. and A.M.G.I. are recipients of a productivity fellowship from CNPq.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Voila, M.; Triches, D. A cadeia de carne de frango: Uma análise dos mercados brasileiro e mundial de 2002 a 2012. Rev. Teor. E Evidência Econômica 2015, 21, 126–148. [Google Scholar] [CrossRef] [Green Version]
  2. Pimentel, D.; Berger, B.; Filiberto, D.; Newton, M.; Wolfe, B.; Karabinakis, E.; Clark, S.; Poon, E.; Abbett, E.; Nandagopal, S. Water resources: Agricultural and environmental issues. Bioscience 2004, 54, 909–918. [Google Scholar] [CrossRef] [Green Version]
  3. Barbut, S.; Sosnicki, A.A.; Lonergan, S.M.; Knapp, T.; Ciobanu, D.C.; Gatcliffe, L.J.; Huff-Lonergan, E.; Wilson, E.W. Progress in reducing the pale, soft and exudative (PSE) problem in pork and poultry meat. Meat Sci. 2008, 79, 46–63. [Google Scholar] [CrossRef]
  4. Knowles, T.G.; Kestin, S.C.; Haslam, S.M.; Brown, S.N.; Green, L.E.; Butterworth, A.; Pope, S.J.; Pfeiffer, D.; Nicol, C.J. Leg Disorders in Broiler Chickens: Prevalence, Risk Factors and Prevention. PLoS ONE 2008, 3, e1545. [Google Scholar] [CrossRef] [Green Version]
  5. Tallentire, C.W.; Leinonen, I.; Kyriazakis, I. Breeding for efficiency in the broiler chicken: A review. Agron. Sustain. Dev. 2016, 36, 66. [Google Scholar] [CrossRef] [Green Version]
  6. Paz, I.C.D.L.A. Problemas locomotores em frangos de corte—Revisão./locomotor problems in broiler chickens—A review. Rev. Bras. Eng. Biossistemas 2008, 2, 263. [Google Scholar] [CrossRef] [Green Version]
  7. Bradshaw, R.H.; Kirkden, R.D.; Broom, D.M. A Review of the Aetiology and Pathology of Leg Weakness in Broilers in Relation to Welfare. Avian Poult. Biol. Rev. 2002, 13, 45–103. [Google Scholar] [CrossRef] [Green Version]
  8. Wideman, R.F. Bacterial chondronecrosis with osteomyelitis and lameness in broilers: A review. Poult. Sci. 2016, 95, 325–344. [Google Scholar] [CrossRef]
  9. Packialakshmi, B.; Rath, N.C.; Huff, W.E.; Huff, G.R. Poultry Femoral Head Separation and Necrosis: A Review. Avian Dis. 2015, 59, 349–354. [Google Scholar] [CrossRef]
  10. Li, P.F.; Zhou, Z.L.; Shi, C.Y.; Hou, J.F. Downregulation of basic fibroblast growth factor is associated with femoral head necrosis in broilers. Poult. Sci. 2015, 94, 1052–1059. [Google Scholar] [CrossRef]
  11. Durairaj, V.; Okimoto, R.; Rasaputra, K.; Clark, F.D.; Rath, N.C. Histopathology and Serum Clinical Chemistry Evaluation of Broilers With Femoral Head Separation Disorder. Avian Dis. 2009, 53, 21–25. [Google Scholar] [CrossRef] [PubMed]
  12. Wideman Bone circulatory disturbances in the development of spontaneous bacterial chondronecrosis with osteomyelitis: A translational model for the pathogenesis of femoral head necrosis. Front. Endocrinol. (Lausanne) 2013, 3, 183. [CrossRef] [Green Version]
  13. Paludo, E.; Ibelli, A.M.G.; Peixoto, J.O.; Tavernari, F.C.; Lima-Rosa, C.A.V.; Pandolfi, J.R.C.; Ledur, M.C. The involvement of RUNX2 and SPARC genes in the bacterial chondronecrosis with osteomyelitis in broilers. Animal 2017, 11, 1063–1070. [Google Scholar] [CrossRef] [PubMed]
  14. Petry, B.; Savoldi, I.R.; Ibelli, A.M.G.; Paludo, E.; de Oliveira Peixoto, J.; Jaenisch, F.R.F.; de Córdova Cucco, D.; Ledur, M.C. New genes involved in the Bacterial Chondronecrosis with Osteomyelitis in commercial broilers. Livest. Sci. 2018, 208, 33–39. [Google Scholar] [CrossRef]
  15. Peixoto, J.D.O.; Savoldi, I.R.; Ibelli, A.M.G.; Cantão, M.E.; Jaenisch, F.R.F.; Giachetto, P.F.; Settles, M.L.; Zanella, R.; Marchesi, J.A.P.; Pandolfi, J.R.; et al. Proximal femoral head transcriptome reveals novel candidate genes related to epiphysiolysis in broiler chickens. BMC Genom. 2019, 20. [Google Scholar] [CrossRef] [Green Version]
  16. de Peixoto, J.O. Improving skeletal structure in meat-type chickens given fast weight gain and large breast yield through genetic selection. CAB Rev. Perspect. Agric. Vet. Sci. Nutr. Nat. Resour. 2020, 15. [Google Scholar] [CrossRef]
  17. Hul, L.M.; Ibelli, A.M.G.; Savoldi, I.R.; Marcelino, D.E.P.; Fernandes, L.T.; Peixoto, J.O.; Cantão, M.E.; Higa, R.H.; Giachetto, P.F.; Coutinho, L.L.; et al. Differentially expressed genes in the femur cartilage transcriptome clarify the understanding of femoral head separation in chickens. Sci. Rep. 2021, 11, 17965. [Google Scholar] [CrossRef]
  18. Gao, Y.; Liu, S.; Huang, J.; Guo, W.; Chen, J.; Zhang, L.; Zhao, B.; Peng, J.; Wang, A.; Wang, Y.; et al. The ECM-cell interaction of cartilage extracellular matrix on chondrocytes. BioMed Res. Int. 2014, 2014, 648459. [Google Scholar] [CrossRef] [Green Version]
  19. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  20. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  21. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  22. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  24. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  25. Mi, H.; Ebert, D.; Muruganujan, A.; Mills, C.; Albou, L.P.; Mushayamaha, T.; Thomas, P.D. PANTHER version 16: A revised family classification, tree-based classification tool, enhancer regions and extensive API. Nucleic Acids Res. 2021, 49, D394–D403. [Google Scholar] [CrossRef]
  26. Ye, J.; Coulouris, G.; Zaretskaya, I.; Cutcutache, I.; Rozen, S.; Madden, T.L. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012, 13, 134. [Google Scholar] [CrossRef] [Green Version]
  27. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef] [PubMed]
  28. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. Revigo summarizes and visualizes long lists of gene ontology terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Heberle, H.; Meirelles, V.G.; da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef]
  31. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’ayan, A. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [Green Version]
  34. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [Green Version]
  35. Kestin, S.C.; Gordon, S.; Su, G.; Sørensen, P. Relationships in broiler chickens between lameness, liveweight, growth rate and age. Vet. Rec. 2001, 148, 195–197. [Google Scholar] [CrossRef]
  36. Jiang, T.; Mandal, R.K.; Wideman, R.F.; Khatiwara, A.; Pevzner, I.; Min Kwon, Y. Molecular Survey of Bacterial Communities Associated with Bacterial Chondronecrosis with Osteomyelitis (BCO) in Broilers. PLoS ONE 2015, 10, e0124403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Fan, Q.; Liu, Z.; Shen, C.; Li, H.; Ding, J.; Jin, F.; Sha, L.; Zhang, Z. Microarray study of gene expression profile to identify new candidate genes involved in the molecular mechanism of leptin-induced knee joint osteoarthritis in rat. Hereditas 2018, 155, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Alcantara, J.; Schryvers, A.B. Transferrin binding protein two interacts with both the N-lobe and C-lobe of ovotransferrin. Microb. Pathog. 1996, 20, 73–85. [Google Scholar] [CrossRef] [PubMed]
  39. Bullen, J.J.; Rogers, H.J.; Griffiths, E. Role of iron in bacterial infection. Curr. Top. Microbiol. Immunol. 1978, 80, 1–35. [Google Scholar] [PubMed]
  40. Giansanti, F.; Leboffe, L.; Angelucci, F.; Antonini, G. The nutraceutical properties of ovotransferrin and its potential utilization as a functional food. Nutrients 2015, 7, 9105–9115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Carlevaro, M.F.; Albini, A.; Ribatti, D.; Gentili, C.; Benelli, R.; Cermelli, S.; Cancedda, R.; Cancedda, F.D. Transferrin promotes endothelial cell migration and invasion: Implication in cartilage neovascularization. J. Cell Biol. 1997, 136, 1375–1384. [Google Scholar] [CrossRef]
  42. HBG2 Gene—GeneCards | HBG2 Protein | HBG2 Antibody.
  43. Alev, C.; Shinmyozu, K.; McIntyre, B.A.S.; Sheng, G. Genomic organization of zebra finch alpha and beta globin genes and their expression in primitive and definitive blood in comparison with globins in chicken. Dev. Genes Evol. 2009, 219, 353–360. [Google Scholar] [CrossRef] [PubMed]
  44. SLC25A37 Gene—GeneCards | MFRN1 Protein | MFRN1 Antibody.
  45. Pizzagalli, M.D.; Bensimon, A.; Superti-Furga, G. A guide to plasma membrane solute carrier proteins. FEBS J. 2021, 288, 2784–2835. [Google Scholar] [CrossRef] [PubMed]
  46. Chigaev, A. Does aberrant membrane transport contribute to poor outcome in adult acute myeloid leukemia? Front. Pharmacol. 2015, 6, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Zhang, L.; Luo, Q.; Shu, Y.; Zeng, Z.; Huang, B.; Feng, Y.; Zhang, B.; Wang, X.; Lei, Y.; Ye, Z.; et al. Transcriptomic landscape regulated by the 14 types of bone morphogenetic proteins (BMPs) in lineage commitment and differentiation of mesenchymal stem cells (MSCs). Genes Dis. 2019, 6, 258–275. [Google Scholar] [CrossRef] [PubMed]
  48. Michaely, P.; Tomchick, D.R.; Machius, M.; Anderson, R.G.W. Crystal structure of a 12 ANK repeat stack from human ankyrinR. EMBO J. 2002, 21, 6387–6396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Bagnato, P.; Barone, V.; Giacomello, E.; Rossi, D.; Sorrentino, V. Binding of an ankyrin-1 isoform to obscurin suggests a molecular link between the sarcoplasmic reticulum and myofibrils in striated muscles. J. Cell Biol. 2003, 160, 245–253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Nummenmaa, E.; Hämäläinen, M.; Pemmari, A.; Moilanen, L.J.; Tuure, L.; Nieminen, R.M.; Moilanen, T.; Vuolteenaho, K.; Moilanen, E. Transient Receptor Potential Ankyrin 1 (TRPA1) Is Involved in Upregulating Interleukin-6 Expression in Osteoarthritic Chondrocyte Models. Int. J. Mol. Sci. 2020, 22, 87. [Google Scholar] [CrossRef] [PubMed]
  51. Scheper, G.C.; Parra, J.L.; Wilson, M.; van Kollenburg, B.; Vertegaal, A.C.O.; Han, Z.-G.; Proud, C.G. The N and C Termini of the Splice Variants of the Human Mitogen-Activated Protein Kinase-Interacting Kinase Mnk2 Determine Activity and Localization. Mol. Cell. Biol. 2003, 23, 5692–5705. [Google Scholar] [CrossRef] [Green Version]
  52. Chrestensen, C.A.; Eschenroeder, A.; Ross, W.G.; Ueda, T.; Watanabe-Fukunaga, R.; Fukunaga, R.; Sturgill, T.W. Loss of MNK function sensitizes fibroblasts to serum-withdrawal induced apoptosis. Genes Cells 2007, 12, 1133–1140. [Google Scholar] [CrossRef] [PubMed]
  53. Daneshyar, M.; Kermanshahi, H.; Golian, A. The effects of turmeric supplementation on antioxidant status, blood gas indices and mortality in broiler chickens with T3-induced ascites. Br. Poult. Sci. 2012, 53, 379–385. [Google Scholar] [CrossRef]
  54. Acharya, K.R.; Ackerman, S.J. Eosinophil Granule Proteins: Form and Function. J. Biol. Chem. 2014, 289, 17406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Davoine, F.; Lacy, P. Eosinophil cytokines, chemokines, and growth factors: Emerging roles in immunity. Front. Immunol. 2014, 5, 570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Jong, E.C.; Henderson, W.R.; Klebanoff, S.J. Bactericidal activity of eosinophil peroxidase. J. Immunol. 1980, 124, 1378–1382. [Google Scholar] [PubMed]
  57. Johnston, R.B.; Baehner, R.L. Chronic granulomatous disease: Correlation between pathogenesis and clinical findings. Pediatrics 1971, 48, 730–739. [Google Scholar] [CrossRef] [PubMed]
  58. Yamaguchi, Y.; Zhang, D.E.; Sun, Z.; Albee, E.A.; Nagata, S.; Tenen, D.G.; Ackerman, S.J. Functional characterization of the promoter for the gene encoding human eosinophil peroxidase. J. Biol. Chem. 1994, 269, 19410–19419. [Google Scholar] [CrossRef]
  59. Celi, P.; Verlhac, V.; Pérez Calvo, E.; Schmeisser, J.; Kluenter, A.M. Biomarkers of gastrointestinal functionality in animal nutrition and health. Anim. Feed Sci. Technol. 2019, 250, 9–31. [Google Scholar] [CrossRef]
  60. Hay, E.D. Collagen and Other Matrix Glycoproteins in Embryogenesis. In Cell Biology of Extracellular Matrix; Springer US: New York, NY, USA, 1991; pp. 419–462. [Google Scholar]
  61. Frantz, C.; Stewart, K.M.; Weaver, V.M. The extracellular matrix at a glance. J. Cell Sci. 2010, 123, 4195–4200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Velleman, S.G. The Role of the Extracellular Matrix in Skeletal Muscle Development. Poult. Sci. 2000, 79, 985–989. [Google Scholar] [CrossRef] [PubMed]
  63. Aran, J.M.; Colomer, D.; Matutes, E.; Vives-Corrons, J.L.; Franco, R. Presence of adenosine deaminase on the surface of mononuclear blood cells: Immunochemical localization using light and electron microscopy. J. Histochem. Cytochem. 1991, 39, 1001–1008. [Google Scholar] [CrossRef] [PubMed]
  64. Desrosiers, M.D.; Cembrola, K.M.; Fakir, M.J.; Stephens, L.A.; Jama, F.M.; Shameli, A.; Mehal, W.Z.; Santamaria, P.; Shi, Y. Adenosine Deamination Sustains Dendritic Cell Activation in Inflammation. J. Immunol. 2007, 179, 1884–1892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Sitkovsky, M.V.; Lukashev, D.; Apasov, S.; Kojima, H.; Koshiba, M.; Caldwell, C.; Ohta, A.; Thiel, M. Physiological Control of Immune Response and Inflammatory Tissue Damage by Hypoxia Inducible Factors and Adenosine A2A Receptors. Annu. Rev. Immunol. 2004, 22, 657–682. [Google Scholar] [CrossRef] [PubMed]
  66. Kumar, V.; Sharma, A. Adenosine: An endogenous modulator of innate immune system with therapeutic potential. Eur. J. Pharmacol. 2009, 616, 7–15. [Google Scholar] [CrossRef] [PubMed]
  67. Mills, J.H.; Alabanza, L.M.; Mahamed, D.A.; Bynoe, M.S. Extracellular adenosine signaling induces CX3CL1 expression in the brain to promote experimental autoimmune encephalomyelitis. J. Neuroinflammation 2012, 9, 683. [Google Scholar] [CrossRef] [Green Version]
  68. Wang, Y.; Wang, Y.; Luo, W.; Song, X.; Huang, L.; Xiao, J.; Jin, F.; Ren, Z.; Wang, Y. Roles of long non-coding RNAs and emerging RNA-binding proteins in innate antiviral responses. Theranostics 2020, 10, 9407–9424. [Google Scholar] [CrossRef]
  69. Strazzulla, L.C.; Cronstein, B.N. Regulation of bone and cartilage by adenosine signaling. Purinergic Signal. 2016, 12, 583–593. [Google Scholar] [CrossRef] [Green Version]
  70. Manson, D.; Diamond, L.; Oudjhane, K.; Hussain, F.B.; Roifman, C.; Grunebaum, E. Characteristic scapular and rib changes on chest radiographs of children with ADA-deficiency SCIDS in the first year of life. Pediatr. Radiol. 2013, 43, 589–592. [Google Scholar] [CrossRef]
  71. Sauer, A.V.; Mrak, E.; Hernandez, R.J.; Zacchi, E.; Cavani, F.; Casiraghi, M.; Grunebaum, E.; Roifman, C.M.; Cervi, M.C.; Ambrosi, A.; et al. ADA-deficient SCID is associated with a specific microenvironment and bone phenotype characterized by RANKL/OPG imbalance and osteoblast insufficiency. Blood 2009, 114, 3216–3226. [Google Scholar] [CrossRef] [PubMed]
  72. ADA Gene—GeneCards | ADA Protein | ADA Antibody.
  73. Ananko, E.A.; Kondrakhin, Y.V.; Merkulova, T.I.; Kolchanov, N.A. Recognition of interferon-inducible sites, promoters, and enhancers. BMC Bioinform. 2007, 8, 56. [Google Scholar] [CrossRef] [Green Version]
  74. Tillgren, V.; Ho, J.C.S.; Önnerfjord, P.; Kalamajski, S. The Novel Small Leucine-rich Protein Chondroadherin-like (CHADL) Is Expressed in Cartilage and Modulates Chondrocyte Differentiation. J. Biol. Chem. 2015, 290, 918–925. [Google Scholar] [CrossRef] [Green Version]
  75. Kania, K.; Colella, F.; Riemen, A.H.K.; Wang, H.; Howard, K.A.; Aigner, T.; Dell’Accio, F.; Capellini, T.D.; Roelofs, A.J.; De Bari, C. Regulation of Gdf5 expression in joint remodelling, repair and osteoarthritis. Sci. Rep. 2020, 10, 157. [Google Scholar] [CrossRef] [PubMed]
  76. Duan, Z.X.; Li, Y.S.; Tu, C.; Xie, P.; Li, Y.H.; Qi, L.; Li, Z.H. Identification of a potential gene target for osteoarthritis based on bioinformatics analyses. J. Orthop. Surg. Res. 2020, 15, 1–10. [Google Scholar] [CrossRef] [PubMed]
  77. Coorens, M.; Banaschewski, B.J.H.; Baer, B.J.; Yamashita, C.; van Dijk, A.; Haagsman, H.P.; Veldhuizen, R.A.W.; Veldhuizen, E.J.A. Killing of Pseudomonas aeruginosa by Chicken Cathelicidin-2 Is Immunogenically Silent, Preventing Lung Inflammation In Vivo. Infect. Immun. 2017, 85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Animals 12 00788 g001 550
Figure 1. Main biological process involved with FHS considering the DE genes shared between articular cartilage and femoral growth plate using the REVIGO tool.
Figure 1. Main biological process involved with FHS considering the DE genes shared between articular cartilage and femoral growth plate using the REVIGO tool.
Animals 12 00788 g001
Animals 12 00788 g002 550
Figure 2. Gene network of the common DE genes between the articular cartilage and femoral growth plate tissues from normal and FHS-affected groups using STRING. Colored circles represent genes, and lines represent the predicted interactions between genes.
Figure 2. Gene network of the common DE genes between the articular cartilage and femoral growth plate tissues from normal and FHS-affected groups using STRING. Colored circles represent genes, and lines represent the predicted interactions between genes.
Animals 12 00788 g002
Animals 12 00788 g003 550
Figure 3. Gene network of unique DE genes in the articular cartilage from normal and FHS-affected broilers using STRING. Colored circles represent genes with similar functions or acting in similar pathways, and lines represent the predicted interactions between genes.
Figure 3. Gene network of unique DE genes in the articular cartilage from normal and FHS-affected broilers using STRING. Colored circles represent genes with similar functions or acting in similar pathways, and lines represent the predicted interactions between genes.
Animals 12 00788 g003
Animals 12 00788 g004 550
Figure 4. Gene network of unique DE genes in the femoral growth plate of normal and FHS-affected broilers using STRING. Colored circles represent genes with similar functions or acting in similar pathways, and lines represent the predicted interactions between genes.
Figure 4. Gene network of unique DE genes in the femoral growth plate of normal and FHS-affected broilers using STRING. Colored circles represent genes with similar functions or acting in similar pathways, and lines represent the predicted interactions between genes.
Animals 12 00788 g004
Animals 12 00788 g005 550
Figure 5. Expression profile (2−ΔCt) and statistical significance between normal and FHS-affected broilers obtained with qPCR for GP (A) and AC (B) for the candidate genes evaluated.
Figure 5. Expression profile (2−ΔCt) and statistical significance between normal and FHS-affected broilers obtained with qPCR for GP (A) and AC (B) for the candidate genes evaluated.
Animals 12 00788 g005
Animals 12 00788 g006 550
Figure 6. (A) Venn diagram constructed with DE genes from normal and femoral head necrosis and/or FHS-affected samples in four evaluated datasets. AC: chicken femoral articular cartilage; GP: chicken femoral growth plate; Hip_AC: human articular cartilage; Blood: human peripheral blood. (B) Venn diagram constructed with DE genes from normal and femoral head necrosis and/or FHS affected samples in three evaluated datasets. AC: chicken femoral articular cartilage; GP: femoral growth plate; Hip_AC: human articular cartilage.
Figure 6. (A) Venn diagram constructed with DE genes from normal and femoral head necrosis and/or FHS-affected samples in four evaluated datasets. AC: chicken femoral articular cartilage; GP: chicken femoral growth plate; Hip_AC: human articular cartilage; Blood: human peripheral blood. (B) Venn diagram constructed with DE genes from normal and femoral head necrosis and/or FHS affected samples in three evaluated datasets. AC: chicken femoral articular cartilage; GP: femoral growth plate; Hip_AC: human articular cartilage.
Animals 12 00788 g006
Animals 12 00788 g007 550
Figure 7. Consequences of all variants identified the articular cartilage and growth plate transcriptomes between normal and affected broilers using the VEP tool.
Figure 7. Consequences of all variants identified the articular cartilage and growth plate transcriptomes between normal and affected broilers using the VEP tool.
Animals 12 00788 g007
Animals 12 00788 g008 550
Figure 8. Gene network constructed using the Stringdb with DE genes and genes with SNPs differing between control and affected groups. The line thickness indicates the strength of data support. The circles represent genes, and different colors represents group of genes with similar functions or acting in similar pathways.
Figure 8. Gene network constructed using the Stringdb with DE genes and genes with SNPs differing between control and affected groups. The line thickness indicates the strength of data support. The circles represent genes, and different colors represents group of genes with similar functions or acting in similar pathways.
Animals 12 00788 g008
Table
Table 1. Primers used in the qPCR analysis of articular cartilage and growth plate of normal and FHS-affected broilers.
Table 1. Primers used in the qPCR analysis of articular cartilage and growth plate of normal and FHS-affected broilers.
GeneEnsembl IDPrimer (5′ to 3′)
C7ENSGALG00000014835F: GTGGTCCTTCCCTGGACATC
Gallus gallus complement C7 R: GTTGGATCGCGCTTCACATT
EPB42ENSGALG00000021230F: GCTACAGACCTGGTACTTGGA
Gallus gallus erythrocyte membrane protein band 4.2 R: GACTTGCAGGGTTCTAACTTCG
MKNK2ENSGALG00000003845F: CATCCTGCAGAGGAACAGCA
Gallus gallus MAPK interacting serine/threonine kinase 2 R: CCTGGAGGCTGCTTTGATGA
SLC4A1ENSGALG00000039978F:TCGAGCTCAAATGCCTGGAC
Gallus gallus solute carrier family 4 member 1 (Diego blood group) R: GGCATCTGCACCTCGTTGTA
SPTBENSGALG00000036805F: CTGGATGCTGATGAGGCCAA
Gallus gallus spectrin beta, erythrocytic R: CTTTCCTCGTCCTTCGGCTT
ANK1 *ENSGALG00000029534F: CCACCATCCCACCATTCAGT
Ankyrin 1 R: ACGGTCACAAACTCCAGCAT
ADA *ENSGALG00000004170F:TTCGGCAAGAAAAGAGGGGT
adenosine deaminase R: GTGTTTGGTAGCTGACGTGC
RHAG *ENSGALG00000016684F:TCTGGAGATCACGGCCTTTG
Rh associated glycoprotein R:GCTCCAATATCTGTGGCCTGA
* Hul et al. [17].
Table
Table 2. Characterization of the transcripts identified in the articular cartilage (AC) and femoral growth plate (GP) tissues.
Table 2. Characterization of the transcripts identified in the articular cartilage (AC) and femoral growth plate (GP) tissues.
ACGP
Annotated TranscriptsTranscript Number%Transcript Number%
IG_V_gene30.02%200.16%
LncRNA4133.40%4543.59%
MiRNA200.16%220.17%
misc_RNA10.01%20.02%
Mt_Rrna20.02%20.02%
protein_coding11,57795.03%11,94294.54%
Pseudogene1551.27%1621.28%
rRNA10.01%00
scaRNA10.01%00
snoRNA80.07%260.21%
snRNA10.01%10.01%
sRNA 10.01%
Total annotated transcripts12,182 12,632
Table
Table 3. Common differentially expressed genes between articular cartilage (AC) and femoral growth plate (GP) with information of the Ensembl ID, gene name, log2 fold change (logFC) and false discovery rate (FDR).
Table 3. Common differentially expressed genes between articular cartilage (AC) and femoral growth plate (GP) with information of the Ensembl ID, gene name, log2 fold change (logFC) and false discovery rate (FDR).
ACGP
Ensembl IDGene NamelogFCFDRlogFCFDR
ENSGALG00000043254EPX4.254.11 × 10−92.590.0009
ENSGALG00000003212TSPO23.667.20 × 10−51.540.005
ENSGALG00000040279RHCE2.020.00081.290.002
ENSGALG00000026518RUNDC3A2.440.0031.100.02
ENSGALG00000039978SLC4A12.700.0031.200.01
ENSGALG00000007447TUBB11.990.0031.280.01
ENSGALG00000011190PLACL22.540.0051.100.0002
ENSGALG00000042105PLCB21.810.0050.670.03
ENSGALG00000014585CCL262.360.0081.120.04
ENSGALG00000021230EPB422.690.0081.330.02
ENSGALG00000004170ADA1.710.0081.010.001
ENSGALG00000015358MYH15−1.910.01−1.380.03
ENSGALG00000016684RHAG2.430.011.340.006
ENSGALG00000029857GIMAP61.790.010.760.01
ENSGALG00000013575IFI62.870.012.220.01
ENSGALG00000009479SAMD9L3.100.011.780.01
ENSGALG00000000378SLC25A371.760.011.320.001
ENSGALG00000041693 2.810.011.330.04
ENSGALG00000003845MKNK21.400.010.630.009
ENSGALG00000026948ADD22.120.011.271.85
ENSGALG00000036805SPTB2.190.011.090.005
ENSGALG00000051068 2.380.021.410.0005
ENSGALG00000003594ANK12.110.020.990.01
ENSGALG00000014736KEL1.990.021.210.005
ENSGALG00000043671 −1.240.02−1.030.02
ENSGALG00000028273HBE13.740.022.240.04
ENSGALG00000039269RNF2132.060.0271.060.03
ENSGALG00000014835C72.840.0271.550.001
ENSGALG00000028357LOC4284212.500.0361.490.006
ENSGALG00000002192PTPRC1.220.0370.610.03
ENSGALG00000045776CPN21.770.0370.900.02
ENSGALG00000014686FBN2−1.080.037−1.720.0005
ENSGALG00000030247TMOD42.050.041.060.01
ENSGALG00000009552LOC4232771.780.041.140.04
ENSGALG00000044326LOC4268201.540.040.800.01
ENSGALG00000006453TF3.040.040.670.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Goldoni, I.; Ibelli, A.M.G.; Fernandes, L.T.; Peixoto, J.d.O.; Hul, L.M.; Cantão, M.E.; Gouveia, J.J.d.S.; Ledur, M.C. Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens’ Femoral Head Separation. Animals 2022, 12, 788. https://doi.org/10.3390/ani12060788

AMA Style

Goldoni I, Ibelli AMG, Fernandes LT, Peixoto JdO, Hul LM, Cantão ME, Gouveia JJdS, Ledur MC. Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens’ Femoral Head Separation. Animals. 2022; 12(6):788. https://doi.org/10.3390/ani12060788

Chicago/Turabian Style

Goldoni, Iara, Adriana Mércia Guaratini Ibelli, Lana Teixeira Fernandes, Jane de Oliveira Peixoto, Ludmila Mudri Hul, Maurício Egídio Cantão, João José de Simoni Gouveia, and Mônica Corrêa Ledur. 2022. "Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens’ Femoral Head Separation" Animals 12, no. 6: 788. https://doi.org/10.3390/ani12060788

APA Style

Goldoni, I., Ibelli, A. M. G., Fernandes, L. T., Peixoto, J. d. O., Hul, L. M., Cantão, M. E., Gouveia, J. J. d. S., & Ledur, M. C. (2022). Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens’ Femoral Head Separation. Animals, 12(6), 788. https://doi.org/10.3390/ani12060788

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop