# Novel Candidate Genes Involved in an Initial Stage of White Striping Development in Broiler Chickens

Suelen Fernandes Padilha, Adriana Mércia Guaratini Ibelli, Jane Oliveira Peixoto, Maurício Egídio Cantão, Gabriel Costa Monteiro Moreira, Lana Teixeira Fernandes, Fernando Castro Tavernari, Marcos Antônio Zanella Morés, Ana Paula Almeida Bastos, Laila Talarico Dias, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Suelen Fernandes Padilha, Adriana Mércia Guaratini Ibelli, Jane Oliveira Peixoto, Maurício Egídio Cantão, Gabriel Costa Monteiro Moreira, et al., Novel Candidate Genes Involved in an Initial Stage of White Striping Development in Broiler Chickens. Animals, 2024, 14, 10.3390/ani14162379. hal-04683411

HAL Id: hal-04683411

https://hal.inrae.fr/hal-04683411v1

Submitted on 2 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Article

# Novel Candidate Genes Involved in an Initial Stage of White Striping Development in Broiler Chickens

Suelen Fernandes Padilha<fnref n="1" />, Adriana Mércia Guaratini Ibelli <fnref n="2" />,<fnref n="3" />,*,†, Jane Oliveira Peixoto <fnref n="2" />,<fnref n="3" />, Maurício Egídio Cantão <fnref n="2" />, Gabriel Costa Monteiro Moreira <fnref n="4" />, Lana Teixeira Fernandes <fnref n="2" />, Fernando Castro Tavernari <fnref n="2" />,<fnref n="5" />, Marcos Antônio Zanella Morés <fnref n="2" />, Ana Paula Almeida Bastos <fnref n="2" />,<fnref n="3" />, Laila Talarico Dias <fnref n="1" />, Rodrigo Almeida Teixeira <fnref n="1" /> and Mônica Corrêa Ledur <fnref n="2" />,<fnref n="5" />,*

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Padilha, S.F.; Ibelli, A.M.G.; Peixoto, J.O.; Cantão, M.E.; Moreira, G.C.M.; Fernandes, L.T.; Tavernari, F.C.; Morés, M.A.Z.; Bastos, A.P.A.; Dias, L.T.; et al. Novel Candidate Genes Involved in an Initial Stage of White Striping Development in Broiler Chickens. *Animals* **2024**, *14*, 2379. https://doi.org/10.3390/ani14162379

Academic Editor: Sylwester Świątkiewicz

Received: 1 June 2024
Revised: 13 August 2024
Accepted: 14 August 2024
Published: 16 August 2024

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** White striping (WS) is a condition frequently observed in poultry farming, which reduces consumer acceptability and causes economic losses for the industry. Although common, the development of WS during the animal's growth remains unclear. Therefore, the objective of this study was to identify genes involved in the occurrence of WS in broiler chickens at 35 days of age and through which processes they act. For this, breast samples were collected from normal and WS-affected chickens, and complete RNA sequencing was performed. Thirty genes showing an altered expression were identified, of which 23 were more expressed and 7 were less expressed in affected chickens compared to normal ones. Fourteen of these genes were related to WS for the first time and are associated with muscle development (CEPBD, DUSP8, METTL21EP, NELL2, UBE3D), lipid metabolism (PDK4, DDIT4, FKBP5, DGAT2, LIPG, TDH, and RGCC), and collagen (COL4A5 and COL4A6). Genes related to muscle, as well as cellular processes, are possibly involved with the initial phase of WS development. Genes linked to lipid metabolism and collagen may be more related to the progression of this condition.

**Abstract:** White striping (WS) is a myopathy characterized by the appearance of white stripes parallel to the muscle fibers in the breast of broiler chickens, composed of adipose and connective tissues. This condition causes economic losses and, although common, its etiology remains poorly understood. Hence, the objective was to identify genes and biological mechanisms involved in the early stages of WS using a paternal broiler line that grows slightly slower than commercial ones, at 35 days of age, through the RNA sequencing of the pectoralis major muscle. Thirty genes were differentially expressed between normal and WS-affected chickens, with 23 upregulated and 7 downregulated in the affected broilers. Of these, 14 genes are novel candidates for WS and are implicated in biological processes related to muscle development (CEPBD, DUSP8, METTL21EP, NELL2, and UBE3D), lipid metabolism (PDK4, DDIT4, FKBP5, DGAT2, LIPG, TDH, and RGCC), and collagen (COL4A5 and COL4A6). Genes related to changes in muscle fiber type and the processes of apoptosis, autophagy, proliferation, and differentiation are possibly involved with the initial stage of WS development. In contrast, the genes linked to lipid metabolism and collagen may have their expression altered due to the progression of the myopathy.

<footnote n="1">¹ Departamento de Zootecnia, Programa de Pós-Graduação em Zootecnia, Universidade Federal do Paraná, Curitiba 80035-050, PR, Brazil; suelenfpadilha7@gmail.com (S.F.P.); lailatalarico@gmail.com (L.T.D.); rteixeiraufpr@gmail.com (R.A.T.)</footnote>

<footnote n="2">² Embrapa Suínos e Aves, Concórdia 89715-899, SC, Brazil; jane.peixoto@embrapa.br (J.O.P.); mauricio.cantao@embrapa.br (M.E.C.); lanatf@yahoo.fr (L.T.F.); fernando.tavernari@embrapa.br (F.C.T.); marcos.mores@embrapa.br (M.A.Z.M.); ana.bastos@embrapa.br (A.P.A.B.)</footnote>

<footnote n="3">³ Programa de Pós-Graduação em Ciências Veterinárias, Universidade Estadual do Centro Oeste, Guarapuava 85040-167, PR, Brazil</footnote>

<footnote n="4">⁴ UVSQ, INRAE, BREED, Université Paris-Saclay, 78350 Jouy-en-Josas, France; gabriel.costa-monteiro-moreira@inrae.fr</footnote>

<footnote n="5">⁵ Programa de Pós-Graduação em Zootecnia, Universidade do Estado de Santa Catarina, UDESC-Oeste, Chapecó 89815-630, SC, Brazil</footnote>

<footnote>* Correspondence: adriana.ibelli@embrapa.br (A.M.G.I.); monica.ledur@embrapa.br (M.C.L.)</footnote>

<footnote>† Current address: Embrapa Pecuária Sudeste, São Carlos 13560-970, SP, Brazil.</footnote>
**Keywords:** gene expression; pectoral myopathy; RNA-Seq; transcriptome

# 1. Introduction

White striping (WS) myopathy is characterized by the appearance of white stripes, composed of adipose and connective tissue, that run parallel to the muscle fibers, predominantly in the pectoralis major muscle and less frequently in the thigh muscles [1,2]. This condition predominantly affects birds with high growth rates, greater breast yields, and heavier weights at slaughter [3,4]. This myopathy is undesirable as it impacts the visual appearance, nutritional value, and processing yield of the meat, as well as reduces consumer acceptability and leads to economic losses [1]. WS is considered one of the most common myopathies in poultry farming, with an incidence rate of 10 to 12% on commercial broiler farms [5,6].

Although it is a common condition, the etiology of WS remains poorly understood. It has been determined that the expression of this pathology is partially genetic, with heritability estimates ranging from low to high (0.19 to 0.65) [7-9]. These studies have also reported genetic correlation between WS and other traits, with estimates ranging from 0.08 to 0.33 for body weight and from 0.03 to 0.68 for breast yield [7,9]. Genome-wide association studies (GWASs) [10-12] have identified several candidate genes for WS and indicated that it is a polygenic trait with no specific genomic region contributing a large portion to the overall genetic variance [13]. However, candidate genes in chromosomes 1, 5, 11, 17, and 20 have been identified [10,14].

Additionally, it has been suggested that the occurrence of WS is related to processes such as hypoxia, oxidative stress, immune system activation, angiogenesis, disruptions in calcium pathways, endoplasmic reticulum (ER) stress, apoptosis, collagen metabolism, autophagy, insulin signaling, circulatory system development, cell response to stimulus, and cytokine production, among others [12,13,15-18]. However, conclusive evidence regarding the etiology of WS and its temporal progression remains elusive, indicating the need for studies to better understand the development of this myopathy over time [18].

Most research has concentrated on hybrid animals from commercial lines with high breast yields. Moreover, few studies have been conducted with animals younger than 42 days, which could further elucidate the development of this myopathy in broiler chickens of different lines, particularly in pure lines undergoing intense selection. Therefore, the objective of this study was to identify genes and biological mechanisms involved in the occurrence of WS in a Brazilian paternal broiler line at 35 days of age.

# 2. Materials and Methods

## 2.1. *Animals and Sample Collection*

Male broilers from a paternal broiler line called TT, developed by Embrapa Swine and Poultry National Research Center, were raised until 35 days of age on the experimental farm of the same institution, located in Concórdia, Santa Catarina State, Brazil. For this line, selection emphasis was on improving body weight, feed conversion, carcass and cuts yield, viability, fertility, and hatchability [19]. The TT line has been used as a reference population for genomic studies in chickens, and for more details on the selection and structure of this broiler line, see Marchesi et al. [20].

At the hatchery, the chicks were vaccinated against fowlpox and Marek's disease. The management conditions adhered to line-specific guidelines, providing water and feed (formulated according to Rostagno et al. [21]) ad libitum. At 35 days of age, the chickens were euthanized by cervical dislocation, following the ethical standards set by the Ethics Committee on Animal Use at Embrapa Swine and Poultry (approval no. 008/2019) and in line with international animal welfare guidelines.

After slaughter, the animals were necropsied, and the breast muscles were visually evaluated for the presence of WS according to the classification by Kuttappan et al. [22].
All affected animals exhibited a mild degree of myopathy. Subsequently, approximately 1 g samples were collected from the cranial region of the pectoralis major muscle (PMM) of eight animals: four samples from animals displaying mild WS (affected) and four from normal broilers (controls), with no evidence of WS or other myopathic changes. These samples were immediately frozen in liquid nitrogen and stored at -80 °C for gene expression analysis. Samples for histological analysis were collected from the opposite side of the breast and preserved in 4% paraformaldehyde.

## 2.2. *Histopathological Analysis*

The samples were sectioned into 5 m sections, dehydrated in alcohol, cleared in xylol, and embedded in paraffin. The tissues were then sectioned into 3 µm sections, mounted on slides, and stained with hematoxylin and eosin for the morphological evaluation and identification of myopathic lesions.

## 2.3. RNA Extraction, Preparation, and Sequencing of RNA-seq Libraries

Approximately 100 mg of PMM tissue was homogenized in a mortar containing liquid nitrogen and combined with 1 mL of Trizol reagent (Invitrogen, Waltham, MA, USA) for total RNA extraction, followed by purification using the RNeasy mini kit (Qiagen, Hilden, Germany). RNA purity and concentration were measured using a Biodrop spectrophotometer (Biochrom, St. Albans, UK), and integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Only samples with RIN > 8.0 were used to prepare libraries using the Illumina Stranded mRNA Prep kit (Illumina, San Diego, CA, USA) with 2 µg of total RNA following the manufacturer's protocol. Library size was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), and quantified by qPCR. Paired-end sequencing (2 × 100 bp) was performed on an Illumina NextSeq 2000 (Illumina, San Diego, CA, USA), at the NGS Soluções Genômicas facility (Piracicaba, SP, Brazil). The sequence files obtained were in fastq format.

## 2.4. RNA-seq Data Analysis

Quality control was conducted using Trimmomatic software v0.39 [23], which removed short sequences (<70 bp), low-quality sequences (QPhred ≤20), and adapter sequences. After quality control, the sequences were mapped against the chicken reference genome (GRCg7b, Ensembl release 111) using STAR software v. 2.7.11a [24], and reads were counted with HTseq-count [25]. Differential expression analysis between normal and WS-affected animals was performed using the edgeR package [26] from R v 4.3.2 [27], and gene hierarchies were visualized in heatmap form. Genes were considered differentially expressed (DE) if they had a false discovery rate (FDR) < 0.05 after correction for multiple testing using the Benjamani-Hochberg (BH) method [28]. DE genes with negative and positive log2 fold change (Log2FC) were considered downregulated and upregulated, respectively, in WS-affected chickens compared to normal chickens.

## 2.5. Functional Annotation and Gene Interaction Networks

Gene ontology analysis was conducted using the annotation packages org.Gg.eg.db [29] and AnnotationDbi [30] from Bioconductor 3.18 [31], utilizing information from *Gallus gallus* (GRCg7b). The functional annotation of the DE genes and the enriched biological processes (BPs) were identified through the clusterProfiler package [32] from Bioconductor, using the MSigDB database [33] with information pertaining to *Gallus gallus* (GRCg7b). Gene interaction analysis was performed with gene expression data from the DE genes using the STRING database [34] with annotations for *Homo sapiens* to predict genetic interactions based on co-expression and co-localization. Finally, genes newly identified as related to WS were termed "novel candidate genes for WS." A search based on the initial and final positions of each of these genes was conducted to determine whether they were located in regions of quantitative trait loci (QTLs) previously identified in the chicken genome using the chickenQTLdb [35].
## 2.6. qPCR Confirmation Analysis

To confirm the RNA-Seq results, 20 samples (10 normal and 10 affected) were collected and underwent total RNA extraction with the Trizol standard protocol (Invitrogen, Waltham, MA, USA). The total RNA was submitted to first-strand cDNA synthesis using SuperScript III First-Strand Synthesis SuperMix (Invitrogen, Waltham, MA, USA), following the manufacturer's instructions.

Five candidate genes widely used as reference were chosen to be evaluated according to their stability to be used as normalizers in the gene expression analysis from the *pectoral major* muscle tissue: *HMBS* (Hydroxymethylbilane Synthase) [36], *MRPS30* (Mitochondrial Ribosomal Protein S30) [37], *RPL30* (Ribosomal Protein L30) [38], *RPL4* (Ribosomal Protein L4) [38], and *RPL5* (Ribosomal Protein L5) [37] (Supplementary File S1, Table S1). The qPCR reactions were performed with QuantStudio 6 Real-Time PCR equipment (Life Technologies, Carlsbad, CA, USA) using the following protocol: 1X GoTaq qPCR Master Mix (Promega, WI, USA), 0.133 uM forward primer, 0.133 uM reverse primer, cDNA in 1:10 dilution, and ultra-pure water (Promega, WI, USA) up to 15 µL of final volume. All reactions were performed in duplicate, using negative controls to detect possible contamination. The cycling conditions were 95 °C for 2 min, 40 cycles of 15 s at 95 °C and 30 s at 60 °C, and a melting curve analysis step was added following the manufacturer protocol. The cycle threshold (Ct), standard deviation (SD), and coefficient of variation (CV) were obtained for each sample. Primer's efficiencies were obtained with Web-based LinRegPCR [39–42]. Ct means were used to assess the gene stability performed with the automated pipeline endoGenes [43]. This pipeline uses the three most known programs to evaluate gene stability, geNorm, NormFinder, and BestKeeper, generating a final ranking using the RankAggreg package with the most stable genes. *HMBS* and *RPL30* were selected as the most stable genes found for pectoral major muscle tissue (Supplementary File S1: Figure S1).

Nine candidate genes were chosen for qPCR confirmation analysis. The primers for the genes: Pyruvate dehydrogenase kinase 4 (*PDK4*), FK506 binding protein 5 (*FKBP5*), CCAAT enhancer binding protein delta (*CEBPD*), Diacylglycerol O-acyltransferase 2 (*DGAT2*), regulator of cell cycle (*RGCC*), cholinergic receptor nicotinic gamma subunit (*CHRNG*), carbonic anhydrase 3A (*CA3A*), and uncoupling protein 3 (*UCP3*) were designed using the PRIMER BLAST program [44] with sequences downloaded from the NCBI [45] and Ensembl 111 [46] databases (Supplementary File S1: Table S2). Primer sequence for four and a half LIM domains 1 (*FHL1*) gene was obtained from a previous study [12] (Supplementary File S1: Table S2). For the qPCR reactions, 1X of GoTaq qPCR Master Mix (Madison, WI, USA), 0.133 uM of forward primer and 0.133 uM reverse primer, cDNA in 1:10 dilution, and ultra-pure water (Promega, WI, USA) were used. qPCR reactions were performed in duplicate in 15 µL of final volume. Subsequently, qPCR reactions were submitted to the QuantStudio 6 Real-Time PCR equipment (Applied Biosystems, Foster City, CA, USA). The reaction conditions to assess the specificities of the amplification were as follows: 95 °C for 2 min, 40 cycles for 15 s at 95 °C, and 30 s at 60 °C. The primer efficiencies were obtained with Web-based LinRegPCR [39–42].

The expression level of target genes was determined using Ct values, and changes in gene expression were calculated with the Pfaffl method, considering primer's efficiency [47], using the geometric mean of the endogenous genes *HMBS* and *RPL30*. A pairwise correlation using the ggpubr package from the R v. 4.3.2 environment [27] was performed to compare and confirm the expression levels between the RNA-Seq (log2 fold-change, log2FC) and qPCR (fold change, FC) data, according to Zhang et al. [48].

## 3. Results

The histopathological analysis and the macroscopic evaluation revealed two types of samples: (1) normal muscles, showing organized muscle fibers with homogeneous size and coloration (Figure 1A,C), and (2) samples affected by WS, characterized by fibers of varying
sizes, mild presence of hypereosinophilic, rounded fibers with hyaline degeneration, and infiltration of adipose tissue within the connective tissue (lipidosis) (Figure 1B,D).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Histopathological analysis of 35-day-old chicken breasts showing microscopic features of normal (A,C) and white striping (B,D) muscles. Presence of hypereosinophilic fibers (arrowhead), grade of muscular degeneration (*) and adipose infiltration (arrow)—hematoxylin and eosin stain.</figcaption></figure>

The sequencing of PMM samples generated approximately 18.2 million reads per sample (minimum: 15.9 million; maximum: 20.6 million). After quality control, about 17.2 million reads/sample remained, of which 88.1% were mapped to the chicken reference genome (GRCg7b, Ensembl release 111), and 75.5% of these reads were mapped to gene regions (Supplementary File S1: Table S3).

The hierarchical cluster analysis of DE genes showed that although distinct expression patterns between normal and affected chickens were observed for some genes, they were not homogeneous across all affected samples (Figure 2A). This profile was primarily influenced by the mild degree of WS present in the affected animals, resulting in less stark contrasts compared to normal broilers.

A total of 7434 genes were expressed in the PMM transcriptome (Supplementary File S1, Table S4), among which 30 were DE between normal and WS-affected chickens (Table 1). Of these DE genes, 23 were upregulated and 7 were downregulated in the affected group compared to the normal counterpart (Figure 2B; Table 1). Fourteen novel candidate genes for WS were identified: *PDK4*, *DDIT4*, *FKBP5*, *DGAT2*, *LIPG*, *TDH*, *RGCC*, *CEPBD*, *DUSP8*, *METTL21EP*, *NELL2*, *UBE3D*, *COL4A5*, and *COL4A6* (Table 1).
<table><caption>**Table 1.** Differentially expressed genes up- and downregulated in WS-affected broilers, according to FDR (novel genes associated with white striping are shown in bold).</caption><thead><tr><th colspan="7">Upregulated</th></tr><tr><th>Chromosome</th><th>Associated Gene Name</th><th>Gene Description</th><th>Ensembl Gene Id</th><th>logFC</th><th>FDR</th></tr></thead><tbody><tr><td>2</td><td><b>PDK4</b></td><td>pyruvate dehydrogenase kinase 4</td><td>ENSGALG00010002741</td><td>3.41421</td><td>1.571 × 10<sup>-24</sup></td></tr><tr><td>1</td><td><b>MYBPC1</b></td><td>myosin binding protein C1</td><td>ENSGALG00010011558</td><td>2.75445</td><td>8.349 × 10<sup>-17</sup></td></tr><tr><td>26</td><td><b>FKBP5</b></td><td>FK506 binding protein 5</td><td>ENSGALG00010026893</td><td>2.40214</td><td>1.345 × 10<sup>-12</sup></td></tr><tr><td>2</td><td><b>CEBPD</b></td><td>CCAAT enhancer binding protein delta</td><td>ENSGALG00010005405</td><td>2.01735</td><td>7.776 × 10<sup>-9</sup></td></tr><tr><td>2</td><td><b>CA3A</b></td><td>carbonic anhydrase 3A</td><td>ENSGALG00010013288</td><td>2.00895</td><td>7.776 × 10<sup>-9</sup></td></tr><tr><td>26</td><td><b>TNNT2</b></td><td>troponin T2, cardiac type</td><td>ENSGALG00010026743</td><td>1.57929</td><td>3.182 × 10<sup>-5</sup></td></tr><tr><td>2</td><td><b>ASNS</b></td><td>asparagine synthetase (glutamine-hydrolyzing)</td><td>ENSGALG00010003382</td><td>1.52297</td><td>9.686 × 10<sup>-5</sup></td></tr><tr><td>1</td><td><b>PFKFB3</b></td><td>6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3</td><td>ENSGALG00010006245</td><td>1.45674</td><td>0.0002375</td></tr><tr><td>1</td><td><b>NELL2</b></td><td>neural EGFL like 2</td><td>ENSGALG00010013749</td><td>1.43497</td><td>0.0006832</td></tr><tr><td>2</td><td><b>FBXO32</b></td><td>F-box protein 32</td><td>ENSGALG00010008245</td><td>1.37948</td><td>0.0006832</td></tr><tr><td>3</td><td><b>SESN1</b></td><td>sestrin 1</td><td>ENSGALG00010014244</td><td>1.34848</td><td>0.0009661</td></tr><tr><td>3</td><td><b>TDH</b></td><td>L-treonina desidrogenase</td><td>ENSGALG00010011278</td><td>1.23466</td><td>0.0064443</td></tr><tr><td>4</td><td><b>MYOZ2</b></td><td>myzozenin 2</td><td>ENSGALG00010005510</td><td>1.21752</td><td>0.0089447</td></tr><tr><td>4</td><td><b>COL4A6</b></td><td>collagen type IV alpha 6 chain</td><td>ENSGALG00010015733</td><td>1.20722</td><td>0.0107838</td></tr><tr><td>4</td><td><b>FHL1</b></td><td>four and a half LIM domains 1</td><td>ENSGALG00010014929</td><td>1.20013</td><td>0.0088713</td></tr><tr><td>9</td><td><b>CHRNG</b></td><td>cholinergic receptor nicotinic gamma subunit</td><td>ENSGALG00010017878</td><td>1.18978</td><td>0.0143605</td></tr><tr><td>5</td><td><b>ASB2</b></td><td>ankyrin repeat and SOCS box containing 2</td><td>ENSGALG00010016987</td><td>1.17026</td><td>0.0122158</td></tr><tr><td>4</td><td><b>COL4A5</b></td><td>collagen type IV alpha 5 chain</td><td>ENSGALG00010015632</td><td>1.14888</td><td>0.0228397</td></tr><tr><td>5</td><td><b>MYBPC3</b></td><td>myosin binding protein C3</td><td>ENSGALG00010022468</td><td>1.14728</td><td>0.0190541</td></tr><tr><td>28</td><td><b>ARRDC2</b></td><td>arrestin domain containing 2</td><td>ENSGALG00010027232</td><td>1.14249</td><td>0.0216413</td></tr><tr><td>1</td><td><b>UCP3</b></td><td>uncoupling protein 3</td><td>ENSGALG00010002543</td><td>1.08505</td><td>0.0337988</td></tr><tr><td>6</td><td><b>DDIT4</b></td><td>DNA damage inducible transcript 4</td><td>ENSGALG00010021509</td><td>1.0787</td><td>0.0370386</td></tr><tr><td>1</td><td><b>RGCC</b></td><td>regulator of cell cycle</td><td>ENSGALG00010003192</td><td>1.07199</td><td>0.0372903</td></tr></tbody></table><table><thead><tr><th colspan="7">Downregulated</th></tr><tr><th>Chromosome</th><th>Associated gene name</th><th>Gene description</th><th>Ensembl gene id</th><th>logFC</th><th>FDR</th></tr></thead><tbody><tr><td>1</td><td><b>METTL21EP</b></td><td>methyltransferase like 21E, pseudogene</td><td>ENSGALG00010006037</td><td>-2.2487</td><td>1.40 × 10<sup>-10</sup></td></tr><tr><td>3</td><td><b>UBE3D</b></td><td>ubiquitin protein ligase E3D</td><td>ENSGALG00010003115</td><td>-2.0996</td><td>4.26 × 10<sup>-9</sup></td></tr><tr><td>6</td><td><b>SCD</b></td><td>stearoyl-CoA desaturase</td><td>ENSGALG00010021083</td><td>-1.9290</td><td>3.38 × 10<sup>-8</sup></td></tr><tr><td>1</td><td><b>DGAT2</b></td><td>diacylglycerol O-acyltransferase 2</td><td>ENSGALG00010003035</td><td>-1.5278</td><td>9.69 × 10<sup>-5</sup></td></tr><tr><td>5</td><td><b>DUSP8</b></td><td>dual specificity phosphatase 8</td><td>ENSGALG00010024637</td><td>-1.3427</td><td>0.0017</td></tr><tr><td>Z</td><td><b>LIPG</b></td><td>lipase G, endothelial type</td><td>ENSGALG00010009013</td><td>-1.0878</td><td>0.0466</td></tr><tr><td>3</td><td><b>ATF3</b></td><td>activating transcription factor 3</td><td>ENSGALG00010019307</td><td>-1.0621</td><td>0.0431</td></tr></tbody></table>

From the nine genes selected for qPCR confirmation, it was possible to verify a high correlation between the expression levels of the RNA-Seq and the qPCR analysis (Figure 3) (r = 0.71), evoking consistent results from the RNA-Seq obtained in this study; therefore, the qPCR analysis validated the RNA-Seq results (Supplementary File S1: Figure S2).

In the search for positional candidate genes, 13 of the novel candidate genes for WS were found in regions associated with QTLs for the following traits: body weight (including weight gain, average daily gain, carcass weight, conformation score, feed conversion,
and feed efficiency), pectoral muscle and PMM (including breast weight, percentage, and width), muscle fiber (including number, diameter, and density), and fat content/percentage (intramuscular and carcass) (Supplementary File S1: Table S5). The *PDK4*, *NELL2*, *DUSP8*, and *COL4A5* genes are notably located in QTL regions associated with many of the aforementioned traits.

<figure><img src="image_7.png" /><figcaption>**Figure 2.** (A) A heatmap showing a hierarchical clustering of the genes and samples of the DE genes between WS-affected and normal broilers. Genes are presented in the rows, and samples, in the columns. Downregulated genes in the affected group are shown in blue, and upregulated genes, in red. (B) A volcano plot showing the main differentially expressed genes in WS-affected broilers (-Log₁₀ *p* < 3.699 is equal to FDR < 0.05).</figcaption></figure>

Gene ontology analysis revealed that DE genes were enriched in 518 BPs, with 23 reaching statistical significance (*p* < 0.05; Supplementary File S1: Table S6). The most representative processes are illustrated in Figure 4.

In Figure 5, it is possible to identify which genes are associated with the most relevant BPs. Some genes are shown to participate in more than one BP, and these can be categorized into three major groups: (1) muscle-related genes (*CHRNG*, *FBXO32*, *MYBPC1*, *ATF3*, *MYBPC3*, *ASB2*, *MYOZ2*, *FHL1*, and *TNNT2*); (2) genes involved in lipid metabolism
and cellular processes (LIPG, UCP3, SESN1, ASNS, DGAT2, SCD, PDK4, and ATF3); and
(3) genes associated with collagen (COL4A5 and COL4A6).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Expression levels of RNA-Seq (log2FC) and qPCR (FC).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Main biological process enriched from combined differentially expressed genes (up- and downregulated) between WS-affected and normal broilers.</figcaption></figure>

The gene network (Figure 5), constructed with the expression data from DE genes, demonstrates interactions among most of the candidate genes, indicating functional associations that may contribute to the manifestation of WS. Additionally, there is a notable connection between this network and significant BPs, as the interacting genes (Figure 5) also participate in the same BPs (Figure 4), which reinforces the importance and functional relationships of these genes.

The interaction among genes from the three large aforementioned groups (related to muscle development, lipid metabolism and cellular processes, and collagen) are shown in
<figure><img src="image_10.png" /><figcaption>Figure 5. Gene enrichment analysis of differentially expressed genes and biological processes related to WS in broiler chickens. Gray circles represent DE genes in WS-affected compared to normal broilers. Flesh-colored circles represent biological processes; their size indicates the number of related genes. Colored lines denote the category of the biological process.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Gene network of differentially expressed genes between WS-affected and normal broilers. Upregulated genes in the affected group are shown in red, and downregulated genes, in blue. Lines indicate gene interactions.</figcaption></figure>
## 4. Discussion

In this study, we investigated the global gene expression profile of the pectoralis major muscle (PMM) of normal and WS-affected broilers from a Brazilian paternal line at 35 days of age. This line has a slightly lower growth rate compared to other commercial lines, which may contribute to a reduced incidence of WS and milder degrees of the myopathy. Furthermore, the samples were obtained from 35-day-old broilers, while most studies published so far have typically involved older animals. It is interesting to note that there are no gene expression studies on WS with broilers that are at ages less than 42 days old, differing from Wooden breast studies. Therefore, it was possible to investigate the gene expression profile of broilers at an earlier stage of the myopathy development, with the identification of 30 DE genes between WS-affected chickens and their normal counterpart. These genes are implicated in biological processes related to muscle development, lipid metabolism and cellular processes, and collagen. Among the DE genes, fourteen candidate genes for WS were identified for the first time: *CEPBD*, *METTL21EP*, *DUSP8*, *NELL2*, *UBE3D*, *PDK4*, *DDIT4*, *FKBP5*, *DGAT2*, *LIPG*, *TDH*, *RGCC*, *COL4A5*, and *COL4A6*. These genes likely emerged, possibly due to the early age and the specific broiler line used, which allowed the evaluation of initial stages of the myopathy. Moreover, these genes are located in regions of QTL associated with traits such as body weight, pectoral muscle and PMM, muscle fiber, and intramuscular and carcass fat content/percentage. These findings highlight the importance of the biological pathways involved (muscle development and lipid metabolism) and reinforce the suggested role of these genes in modifying the muscle fiber type, a condition described in the literature as linked to breast myopathies [49], and in the accumulation of fat between muscle fibers, as observed in the histopathological analysis (Figure 1B).

### 4.1. Muscle-Related Genes

The association of genes related to muscle development processes with WS myopathy in broiler chickens has been previously reported [14,18,50]. In our study, muscle-related BPs were significant (Figure 3 and Supplementary File S1: Table S6), involving genes such as *MYBPC1*, *MYBPC3*, *MYOZ2*, *FHL1*, and *TNNT2*. These genes interact with each other and with *CEBPD*, *ATF3*, and *FBXO32* (Figure 5), and, along with *CHRNG* and *ASB2*, participate in these BPs (Supplementary File S1: Table S6). Except for *ATF3*, which was downregulated, all other genes in this group were upregulated in broilers affected with WS.

The overexpression of genes such as *MYBPC1*, *MYBPC3*, *MYOZ2*, and *TNNT2* has already been linked to muscle growth and the development of pectoral myopathies, including WS [12,49,51–53]. The myosin-binding protein C (MyBP-C) family is a group of sarcomeric proteins that are crucial for the structure and function of striated muscle. This family is composed of three isoforms, cardiac, slow skeletal, and fast skeletal, which are encoded by *MYBPC3*, *MYBPC1*, and *MYBPC2*, respectively [54]. The upregulation of these genes may shift muscle fibers from fast (glycolytic) to slow (oxidative), altering muscle structure [49,51–53] and affecting the contractile, metabolic, and biochemical properties of the muscle [55]. It has been shown that higher levels of the *MYOZ2* gene can decrease the expression of calcineurin, which is involved in the development and differentiation of skeletal muscle as well as in controlling muscle fiber type [56,57]. The *TNNT2* gene also regulates calcium for the contraction of striated muscles [58]. The overexpression of this gene reduces the calcium load in muscles [59], and disruptions in calcium pathways can trigger WS in chickens [18].

Another gene that could influence muscle fiber type in response to metabolic imbalances is *FHL1*. The upregulation of this gene enhances the expression of the slow oxidative fiber type [60] and has also been associated with WS [12,61] and Wooden Breast (WB) in commercial broiler chickens [62]. The overexpression of the *CA3A* gene (Figure 2 and Table 1) supports the hypothesis of a shift in muscle fiber type in WS-affected chickens, as *CA3A* is a marker for slow-twitch fiber types [63]. *CA3A* can be regulated by hypoxia-inducible factor 1 (HIF-1) and has been reported to be upregulated in chickens affected
with WB in a pure line at 47 days [49] and with WS in the Cobb line at 42 and 56 days of age [12,61]. Hence, the upregulation of *CA3A* in WS-affected chickens could indicate a hypoxic condition that may have led to tissue damage, resulting in transitioning muscle fibers to a slow (oxidative) type, thereby increasing *CA3A* expression.

*METTL21EP* was identified in the present study as a novel candidate gene for WS, possibly due to the early age and the broiler line used, and its downregulation might also be a consequence of the shift in muscle fiber type. Originating as a pseudogene from *METTL21E* [64], which is predominant in type IIb myofibers, its reduced expression can lead to decreased myofiber size [65]. Thus, the reduced expression of *METTL21EP* in WS-affected chickens likely reflects a shift to predominantly slow (type I, oxidative) muscle fibers. Moreover, this gene was recently reported as DE in chickens with WB at 47 days [66]

The expression of the *DUSP8* gene can be altered in response to oxidative stress [67], and its inhibition shifts skeletal muscle fibers to an oxidative type [68]. *DUSP8* is also located in a region of QTLs associated with traits such as weight gain, carcass weight, conformation, and breast weight and yield, as well as muscle fiber diameter and density (Supplementary File S1: Table S5), highlighting its potential role in the onset of WS. The inhibition of *DUSP8* promotes autophagy, while its increased expression may inhibit cell proliferation and induce apoptosis [69,70]. Therefore, the observed downregulation of *DUSP8* in WS-affected chickens, described for the first time in this study, could be a response to oxidative stress, promoting the autophagy of damaged muscle fibers and proliferative processes typical of oxidative (slow) fibers. This response could enhance lipid oxidation, leading to fat accumulation in the muscles, a hallmark of the WS (Figure 1B). Additionally, in this study, the biological processes related to lipid metabolism were significant (Figure 3).

Further discussing muscle regulation and functioning, additional genes—*FBXO32*, *CEBPD*, *ASB2*, and *CHRNG*—have also been linked to apoptosis and differentiation, processes previously associated with the development of WS [11,50]. The *FBXO32* gene is a key mediator of apoptosis following ischemia injury [71], and apoptosis can lead to tissue damage that requires the activation of repair mechanisms [12]. However, the upregulation of this gene suppresses differentiation and inhibits myotube formation. As this process is essential for repair after injury, it is suggested that the overexpression of *FBXO32* adversely impacts the regeneration of skeletal muscle [72]. In addition, Marchesi et al. [12] observed the downregulation of the *FBXO32* gene in the breast muscles of Cobb chickens with WS at 42 days. In contrast, in the current study involving the TT line chickens with 35 days of age, the *FBXO32* gene was found to be upregulated. The opposite expression levels of this gene observed in both studies may be attributed to the line used, the progression of the myopathy and the degree of fiber regeneration. Moreover, Li et al. [73] also reported a difference in the expression of this gene between pure breeds and a commercial line of broiler chickens.

In our study, the *CEBPD* gene was identified as a novel candidate gene for WS. Although it is not directly linked to skeletal muscle BP (Supplementary File S1: Table S6), it is located in a QTL region associated with breast muscle weight (Supplementary File S1: Table S6). This gene encodes a protein involved in various BPs, including cell differentiation, proliferation, apoptosis, and immune response [74,75]. *CEBPD* is also highly expressed under hypoxic conditions [76] and plays a crucial role in mediating inflammatory responses [77]. Thus, it is hypothesized that WS-affected broilers experienced muscular hypoxia, leading to the overexpression of the *CEBPD* gene. Nonetheless, the specific mechanisms through which this gene contributes to myopathy development need further elucidation.

On the other hand, although the *ASB2* and *CHRNG* genes do not interact with other muscle-related genes, they are involved in BPs related to muscle processes (Figure 4; Supplementary File S1: Table S5). The *ASB2* gene regulates myogenic differentiation through the ubiquitination and degradation of filamin B (FLNb) [78], while the *CHRNG* gene influences myotube differentiation [79]. Therefore, we suggest that in WS-affected animals, both genes (*ASB2* and *CHRNG*), which were upregulated, enhance myogenic differentiation and favor the regeneration of tissue damaged by the possible hypoxic conditions.
Some genes, such as *NELL2* and *ARRDC2* (upregulated), as well as *UBE3D* (downregulated), are not involved in muscle-related BPs (Supplementary File S1: Table S6). However, this study is the first to identify *NELL2* and *UBE3D* as candidate genes for WS. *NELL2* is found in the QTL region associated with body weight, carcass weight, breast weight, intramuscular fat content, and muscle fiber (Supplementary File S1: Table S5). It is known that *NELL2* overexpression may promote cell proliferation and inhibit apoptosis [80]. *ARRDC2* expression increases following catabolic stimuli and might be involved in autophagy activation [81]. Hence, it is suggested that the upregulation of *ARRDC2* in WS-affected broilers was a response to muscle injuries from apoptotic processes and potential oxidative stress. The upregulation of *NELL2* is thought to promote the proliferation and differentiation of muscle cells, aiding in the regeneration of damaged tissue. Finally, although the *UBE3D* gene has been minimally studied, it is located in QTL regions associated with body weight (Supplementary File S1: Table S5). Some reports indicate that its low expression in zebrafish embryos is linked to oxidative damage and inflammatory responses, and that inhibiting this gene can lead to apoptosis and retinal degeneration [82]. Therefore, it is assumed that the downregulation of this gene in the PMM of chickens could have caused oxidative damage and inflammatory reactions, inducing apoptosis and muscle fiber degeneration, which are hallmark histopathological features of WS.

Changes in the expression profiles of genes associated with muscular processes may indicate alterations that are potentially involved in the initial stage of WS development. The PMM is predominantly composed of fast-twitch-type IIb glycolytic muscle fibers [83]. Nonetheless, several factors such as activity rate, mechanical stress, hormones, aging, responses to muscle damage, and/or changes in muscle oxygenation can lead to fiber-type switching [84]. In this way, the upregulation of the *MYBPC1*, *MYBPC3*, *MYOZ2*, *TNNT2*, *FHL1*, and *CA3A* genes, along with the downregulation of *METTL21EP* and *DUSP8*, which are associated with changes in muscle fiber type, are important findings for understanding the alterations in the muscle tissue affected by WS. This study is the first to identify genes associated with muscle fiber type changes in tissues exclusively affected by WS.

## 4.2. Genes Related to Lipid Metabolism and Cellular Processes

The genetic selection applied to improve broiler chickens has led to reduced muscle oxygenation (hypoxia), which, combined with oxidative stress, contributes to disorders in the oxidation of fatty acids and the tricarboxylic acid cycle, resulting in muscle problems [15]. Significant BPs related to lipid metabolism and cellular processes were identified in this study (Figure 3 and Supplementary File S1: Table S6). These processes involve *PDK4* and *DDIT4*, novel candidate genes for WS, which interact with several other genes (*UCP3*, *PFKFB3*, *FKBP5*, *SCD*, *DGAT2*, *LIPG*, *SESN1*, *ASNS*, and *ATF3*) implicated in these BPs (Figures 4 and 5). With the exception of *FKBP5*, all have been linked to breast myopathies [12,52,85]. These findings support the hypotheses of Lake and Abasht [86] and Pizzol et al. [17], who suggested that WS may be caused by dysregulation in glucose and lipid metabolism. Additionally, cellular processes such as autophagy, apoptosis, and differentiation have also been linked to WS [11,17,18,50] and are likely activated in response to the onset of the myopathy.

Several genes identified in the current study are regulated by hypoxic conditions, e.g., *PDK4*, *DDIT4*, *PFKFB3*, *SESN1*, and *CEBPD*. Changes in muscle oxygenation can induce fiber-type switching [84] and activate various cellular signaling pathways that control metabolism. These pathways stimulate the expression of the *PDK4* gene via the nuclear receptor ERRγ [87,88]. In our study, *PDK4* was found to be more expressed in WS-affected chickens, enriched in the BP for fatty acid response (Figure 4), and located in QTL regions for intramuscular fat percentage (Supplementary File S1: Table S5). The increased expression of this gene suggests the enhanced oxidation of fatty acids [89]. Given that WS-affected chickens exhibit high lipid content in the pectoral muscle [90], it is posited that the excessive oxidation of fatty acids could have led to mitochondrial overload and incomplete oxidation, resulting in the intramuscular accumulation of metabolites derived
from fatty acids [91]. Moreover, as oxygen is necessary for lipid oxidation, the likely hypoxic condition could have led to incomplete oxidation and intensified the accumulation of fatty acids in the muscle, as evidenced by histopathological analysis (Figure 1B).

The *UCP3* gene interacts with *PDK4* (Figure 5), and, in some way, the overexpression of *PDK4* accounts for the observed increase in *UCP3* expression. This upregulation decreases mitochondrial membrane potential, enhances the ATP/ADP ratio, and favors fatty acid oxidation over glucose utilization [92,93]. *UCP3* primarily facilitates the transport of fatty acid anions from muscle mitochondria, protecting against the detrimental effects of high fatty acid levels [94]. This gene was found to be downregulated in Ross chickens exhibiting severe levels of WS [11]. However, unlike in our study, the comparison was between chickens with severe and mild WS.

Additionally, the *SCD* and *DGAT2* genes also interact with *PDK4* (Figure 5), with *DGAT2* implicated in the BP related to fatty acid response (Figure 4). Nevertheless, these genes were downregulated in WS-affected chickens (Figure 2 and Table 1). *SCD* inhibition enhances β-oxidation in muscles [95], while *DGAT2* inhibition in skeletal myotubes reduces glucose uptake and the partitioning of fatty acids into triacylglycerols (TAGs) [96,97]. Thus, *DGAT2* downregulation may reduce TAG formation and elevate harmful fatty acid levels in muscle tissue (Figure 1B), leading to lipotoxicity and the development of pectoral myopathy [85]. Similarly, *LIPG* downregulation disrupts lipid metabolism, potentially triggering pectoral myopathy alongside other lipid-related genes [52]. Both *DGAT2* and *LIPG* have been previously associated with WB pectoral myopathy [52,85], but are linked to WS for the first time in the present study.

The *TDH* and *RGCC* genes were upregulated in the PMM of WS-affected chickens (Figure 2; Table 1) and, although not involved in the BPs identified in this study (Figure 4), were related to WS for the first time, establishing them as key functional candidates for myopathy development. The *TDH* gene regulates fatty acid degradation, biosynthesis, and ketobody formation [98] and may influence the development of chicken breast muscle fibers [99]. Positioned in QTL for carcass fat content (Supplementary File S1: Table S5), the *RGCC* gene regulates cell differentiation [100], and when upregulated, *RGCC* promotes fatty acid oxidation [101,102], contributing to metabolic imbalance and lipid accumulation in muscle (Figure 1B), which is indicative of WS.

In addition to activating various metabolic control pathways, hypoxic conditions induce the expression of *DDIT4*, whose overexpression can promote apoptosis [103] and inhibit mTOR signaling [104], limiting the autophagy [104,105]. Therefore, the upregulation of *DDIT4* in WS-affected chickens, now identified as a novel candidate gene, could suppress the mTOR pathway. This suppression may enhance apoptosis and autophagy processes to mitigate damage from hypoxia and oxidative stress, serving as an important mechanism for degrading damaged cells in pectoral muscles affected by WS [106].

The *PFKFB3* and *FKBP5* genes interact with both *PDK4* and *DDIT4* (Figure 5). *FKBP5* has emerged as a novel candidate gene for WS. During hypoxic conditions, the hypoxia-inducible factor 1 (HIF-1) complex activates *PFKFB3*, promoting its expression [107]. The overexpression of *PFKFB3* increases cell proliferation [108] and modulates autophagy [109]. Although *FKBP5* does not participate in lipid metabolism or cellular processes (Figure 4 and Supplementary File S1: Table S6), it can enhance autophagy [110] and trigger p53-dependent apoptosis [111]. Furthermore, *FKBP5* upregulation increases muscle cell differentiation and tissue regeneration [112]. Thus, the likely hypoxic environment in PMM affected by WS could stimulate increased expression of both *PFKFB3* and *FKBP5*, which may lead to the apoptosis and autophagy of injured cells while enhancing cell proliferation and differentiation, aiding in tissue regeneration.

In addition to the *PFKFB3* and *FKBP5* genes, *DDIT4* formed a cluster including the *SESN1*, *ASNS*, and *ATF3* genes (Figure 5). Sestrin family genes, such as *SESN1*, are typically induced under conditions of DNA damage, hypoxia, and oxidative stress [113,114]. The overexpression of *SESN1* in chicken skeletal muscle cells stimulates cell proliferation and myoblast differentiation [115]. Similarly, the overexpression of *ASNS* promotes proliferation
and resistance to oxidative stress in human stem cells [116]. Therefore, we hypothesize that the overexpression of *SESN1* and *ASNS* could promote muscle proliferation and differentiation in order to control oxidative stress and repair the damaged tissue.

The *ATF3* gene was found to be downregulated in an initial stage of WS development (Figure 2; Table 1). Studies have shown that increased *ATF3* expression can reduce β-oxidation rates, while decreased expression may lessen apoptosis and suppress both the production of reactive oxygen species (ROS) and inflammatory responses [117–119]. Therefore, the downregulation of *ATF3* in the pectoral muscles of WS-affected chickens might elevate β-oxidation rates. Alongside the other genes discussed, this could deregulate lipid metabolism, favoring lipid accumulation in muscle fibers (observed in Figure 1B) and contributing to the development of myopathy. Furthermore, the downregulation of *ATF3* might reflect an adaptive response to stabilize oxidative stress, apoptosis, and inflammation during the development of the myopathy, thereby aiding in the repair of muscle tissue damage. However, this gene was previously reported to be upregulated in a study from our group comparing the transcriptome of normal and WS-affected commercial broilers at 42 days of age [12]. This discrepancy indicates that the expression of this gene is modulated by the lipid accumulation and tissue regeneration.

### 4.3. Collagen-Related Genes

Research has documented an association between collagen genes and both WS and WB in broiler chickens [14,16,120]. Additionally, WS-affected chicken breast filets exhibit higher collagen content, including type IV collagen [121–123]. In our study, collagen-related BPs were significant (Figures 3 and 4 and Supplementary File S1: Table S6), and the *COL4A5* and *COL4A6* genes, found to be upregulated in WS-affected chickens (Figure 2 and Table 1), were identified as novel candidate genes for WS.

Type IV collagen is the primary structural component of the basement membrane, as it forms a supportive network beneath epithelial and endothelial cells and serves as a barrier between tissue compartments [124]. The upregulation of the *COL4A5* and *COL4A6* genes in affected animals could likely result from myofibril necrosis [125], which may lead to collagen accumulation. The accumulation of type IV collagen can induce ER stress, causing alterations in the structure of the extracellular matrix and apoptosis in the pectoral muscles affected by myopathies [120].

Therefore, it is plausible that an imbalance of hypoxia and oxidative stress mechanisms in the pectoralis major muscle may have facilitated the emergence of tissue injuries. For tissue regeneration, apoptosis and autophagy could be induced, along with proliferation and myogenic differentiation. This could result in a transition of muscle fibers from fast (glycolytic) to slow (oxidative) types. The predominance of slow-type fibers could have disrupted metabolism, leading to enhanced fatty acid degradation and/or incomplete oxidation. Excessive or incomplete lipid oxidation may cause lipid accumulation in muscle, a phenomenon partly characteristic of WS myopathy. On the other hand, tissue injuries might also trigger excessive collagen production in an attempt to rebuild the muscle's basement membrane. This collagen accumulation could induce endoplasmic reticulum stress, further contributing to the development of WS. Further studies with this broiler line at an earlier age could help elucidate the mechanisms activating muscle genes and determine whether there was indeed a state of hypoxia and oxidative stress.

## 5. Conclusions

In the initial stages of WS, observed in the evaluated broiler line at 35 days of age, a set of 30 differentially expressed genes was identified. Among these, *CEBPD*, *METTL21EP*, *DUSP8*, *NELL2*, *UBE3D*, *PDK4*, *DDIT4*, *FKBP5*, *DGAT2*, *LIPG*, *TDH*, *RGCC*, *COL4A5*, and *COL4A6* were associated with this myopathy for the first time. Thirteen of these novel candidates for white striping were located in QTL regions linked to the traits of body weight, carcass and breast weight, muscle fibers, and intramuscular fat content. Moreover, 15 genes were associated with muscle biological processes, including 8 linked to changes
in muscle fiber type. These, along with other genes involved in apoptosis, autophagy, proliferation, and differentiation, are highlighted as likely involved in the early stages of WS development. This finding demonstrates the potential shift from glycolytic to oxidative muscle fiber types in WS-affected chickens and suggests mechanisms for lipid accumulation in muscles exhibiting WS. In this sense, eight genes involved in lipid metabolism and two related to collagen were also identified, whose altered expression may be due to the development of the myopathy.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ani14162379/s1: Supplementary File S1: Table S1: Primers for the reference candidate genes evaluated in pectoral major muscle tissue; Supplementary File S1: Figure S1: Overall ranking of reference candidate genes based on the Genorm, Normfinder, and Bestkeeper software ranked using RankAggreg tool; Supplementary File S1: Table S2: Primers for the candidate genes evaluated in pectoral major muscle tissue; Supplementary File S1: Table S3: Number of raw reads, after quality control and mapped reads; Supplementary File S1: Table S4: Genes expressed in the pectoralis major muscle of normal and WS-affect broilers; Supplementary File S1: Figure S2: Pearson correlation (r) obtained between RNA-Seq and qPCR of the nine selected genes. The straight line is the linear regression, the genes are shown by the black dots, and the gray area represents the 95% confidence interval; Supplementary File S1: Table S5: List of QTLs located in the region of novel candidate genes for white striping in broilers; Supplementary File S1: Table S6: Biological processes enriched with differentially expressed genes between white striping-affected and normal broilers.

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

**Funding:** This study was financed by Project Number 10.19.03.034.00.00 from the Brazilian Agricultural Research Corporation (EMBRAPA), National Council for Scientific and Technological Development (CNPq) through Project Number 423159/2021-6, and The Santa Catarina State Research Foundation through Project Number 2021TR001463.

**Institutional Review Board Statement:** The animal study protocol was approved by the Ethics Committee on the Use of Animals (CEUA) from the Embrapa Swine and Poultry National Research Center (case no. 008/2019, approved on 1 August 2019) for studies involving animals.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data are available in the article.

**Acknowledgments:** SFP was sponsored by a scholarship from the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil. LTF received a visiting specialist fellowship from the National Council of Scientific and Technological Development (CNPq)/Embrapa. MCL, FCT, and AMGI are recipients of a productivity fellowship from CNPq. We thank CAPES (Finance Code 001) for the free access to the journals used in the literature review.

**Conflicts of Interest:** The authors declare no conflicts of interest.
## References

1. Praud, C.; Jimenez, J.; Pampouille, E.; Couroussé, N.; Godet, E.; Le Bihan-Duval, E.; Berri, C. Molecular Phenotyping of White Striping and Wooden Breast Myopathies in Chicken. *Front. Physiol.* **2020**, *11*, 549704. [CrossRef] [PubMed]

2. Velleman, S.G. Recent Developments in Breast Muscle Myopathies Associated with Growth in Poultry. *Annu. Rev. Anim. Biosci.* **2019**, *7*, 289–308. [CrossRef] [PubMed]

3. Kuttappan, V.A.; Brewer, V.B.; Apple, J.K.; Waldroup, P.W.; Owens, C.M. Influence of Growth Rate on the Occurrence of White Striping in Broiler Breast Fillets. *Poult. Sci.* **2012**, *91*, 2677–2685. [CrossRef] [PubMed]

4. Lorenzi, M.; Mudalal, S.; Cavani, C.; Petracci, M. Incidence of White Striping under Commercial Conditions in Medium and Heavy Broiler Chickens in Italy. *J. Appl. Poult. Res.* **2014**, *23*, 754–758. [CrossRef]

5. Ferreira, T.Z.; Casagrande, R.A.; Vieira, S.L.; Driemeier, D.; Kindlein, L. An Investigation of a Reported Case of White Striping in Broilers. *J. Appl. Poult. Res.* **2014**, *23*, 748–753. [CrossRef]

6. Petracci, M.; Mudalal, S.; Bonfiglio, A.; Cavani, C. Occurrence of White Striping under Commercial Conditions and Its Impact on Breast Meat Quality in Broiler Chickens. *Poult. Sci.* **2013**, *92*, 1670–1675. [CrossRef] [PubMed]

7. Alnahhas, N.; Berri, C.; Chabault, M.; Chartrin, P.; Boulay, M.; Bourin, M.C.; Bihan-Duval, E. Le Genetic Parameters of White Striping in Relation to Body Weight, Carcass Composition, and Meat Quality Traits in Two Broiler Lines Divergently Selected for the Ultimate PH of the Pectoralis Major Muscle. *BMC Genet.* **2016**, *17*, 61. [CrossRef]

8. Bailey, R.A.; Souza, E.; Avendano, S. Characterising the Influence of Genetics on Breast Muscle Myopathies in Broiler Chickens. *Front. Physiol.* **2020**, *11*, 547817. [CrossRef]

9. Bailey, R.A.; Watson, K.A.; Bilgili, S.F.; Avendano, S. The Genetic Basis of Pectoralis Major Myopathies in Modern Broiler Chicken Lines. *Poult. Sci.* **2015**, *94*, 2870–2879. [CrossRef]

10. Lake, J.A.; Dekkers, J.C.M.; Abasht, B. Genetic Basis and Identification of Candidate Genes for Wooden Breast and White Striping in Commercial Broiler Chickens. *Sci. Rep.* **2021**, *11*, 6785. [CrossRef]

11. Malila, Y.; Uengwetwanit, T.; Arayamethakorn, S.; Srimarut, Y.; Thanatsang, K.V.; Soglia, F.; Strasburg, G.M.; Rungrassamee, W.; Visessanguan, W. Transcriptional Profiles of Skeletal Muscle Associated with Increasing Severity of White Striping in Commercial Broilers. *Front. Physiol.* **2020**, *11*, 521734. [CrossRef] [PubMed]

12. Marchesi, J.A.P.; Ibelli, A.M.G.; Peixoto, J.O.; Cantão, M.E.; Pandolfi, J.R.C.; Marciano, C.M.M.; Zanella, R.; Settles, M.L.; Coutinho, L.L.; Ledur, M.C. Whole Transcriptome Analysis of the Pectoralis Major Muscle Reveals Molecular Mechanisms Involved with White Striping in Broiler Chickens. *Poult. Sci.* **2019**, *98*, 590–601. [CrossRef] [PubMed]

13. Vanderhout, R.J.; Abdalla, E.A.; Leishman, E.M.; Barbut, S.; Wood, B.J.; Baes, C.F. Genetic Architecture of White Striping in Turkeys (*Meleagris gallopavo*). *Sci. Rep.* **2024**, *14*, 9007. [CrossRef] [PubMed]

14. Pampouille, E.; Berri, C.; Boitard, S.; Hennequet-Antier, C.; Beaucercq, S.A.; Godet, E.; Praud, C.; Jégo, Y.; Le Bihan-Duval, E. Mapping QTL for White Striping in Relation to Breast Muscle Yield and Meat Quality Traits in Broiler Chickens. *BMC Genom.* **2018**, *19*, 202. [CrossRef]

15. Boerboom, G.; Van Kempen, T.; Navarro-Villa, A.; Pérez-Bonilla, A. Unraveling the Cause of White Striping in Broilers Using Metabolomics. *Poult. Sci.* **2018**, *97*, 3977–3986. [CrossRef] [PubMed]

16. Bordini, M.; Soglia, F.; Davoli, R.; Zappaterra, M.; Petracci, M.; Meluzzi, A. Molecular Pathways and Key Genes Associated with Breast Width and Protein Content in White Striping and Wooden Breast Chicken Pectoral Muscle. *Front. Physiol.* **2022**, *13*, 936768. [CrossRef]

17. Pizzol, M.S.D.; Ibelli, A.M.G.; Cantão, M.E.; Campos, F.G.; de Oliveira, H.C.; de Oliveira Peixoto, J.; Fernandes, L.T.; de Castro Tavernari, F.; Morés, M.A.Z.; Bastos, A.P.A.; et al. Differential Expression of MiRNAs Associated with Pectoral Myopathies in Young Broilers: Insights from a Comparative Transcriptome Analysis. *BMC Genom.* **2024**, *25*, 104. [CrossRef]

18. Marciano, C.M.M.; Ibelli, A.M.G.; Marchesi, J.A.P.; de Oliveira Peixoto, J.; Fernandes, L.T.; Savoldi, I.R.; do Carmo, K.B.; Ledur, M.C. Differential Expression of Myogenic and Calcium Signaling-Related Genes in Broilers Affected with White Striping. *Front. Physiol.* **2021**, *12*, 712464. [CrossRef]

19. Fornari, M.B.; Zanella, R.; Ibelli, A.M.G.; Fernandes, L.T.; Cantão, M.E.; Thomaz-Soccol, V.; Ledur, M.C.; Peixoto, J.O. Unraveling the Associations of Osteoprotegerin Gene with Production Traits in a Paternal Broiler Line. *Springerplus* **2014**, *3*, 682. [CrossRef]

20. Marchesi, J.A.P.; Buzanskas, M.E.; Cantão, M.E.; Ibelli, A.M.G.; Peixoto, J.O.; Joaquim, L.B.; Moreira, G.C.M.; Godoy, T.F.; Sbardella, A.P.; Figueiredo, E.A.P.; et al. Relationship of Runs of Homozygosity with Adaptive and Production Traits in a Paternal Broiler Line. *Animal* **2018**, *12*, 1126–1134. [CrossRef]

21. Rostagno, H.S.; Albino, L.T.; Donzele, J.L.; Gomes, P.C.; Oliveira, R.D.; Lopes, D.C.; Ferreira, A.S.; Barreto, S.D.T.; Euclides, R.F. *Tabelas Brasileiras Para Aves e Suínos: Composição de Alimentos e Exigências Nutricionais*; Rostagno, H.S., Ed.; Universidade Federal de Viçosa: Viçosa, Brasil, 2017.

22. Kuttappan, V.A.; Lee, Y.S.; Erf, G.F.; Meullenet, J.F.C.; Mckee, S.R.; Owens, C.M. Consumer Acceptance of Visual Appearance of Broiler Breast Meat with Varying Degrees of White Striping. *Poult. Sci.* **2012**, *91*, 1240–1247. [CrossRef]

23. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. *Bioinformatics* **2014**, *30*, 2114–2120. [CrossRef] [PubMed]

24. 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. [CrossRef] [PubMed]
25. Anders, S.; pyl, P.T.; Huber, W. HTSeq—A Python Framework to Work with High-Throughput Sequencing Data. *Bioinformatics* **2015**, *31*, 166–169. [CrossRef] [PubMed]

26. Chen, Y.; Chen, L.; Lun, A.T.L.; Baldoni, P.L.; Smyth, G.K. EdgeR 4.0: Powerful Differential Analysis of Sequencing Data with Expanded Functionality and Improved Support for Small Counts and Larger Datasets. *bioRxiv* **2024**. [CrossRef]

27. R Core Team. R: A Language and Environment for Statistical Computing. Available online: https://www.r-project.org (accessed on 30 January 2024).

28. 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. [CrossRef]

29. Carlson, M. Org.Gg.Eg.Db: Genome Wide Annotation for Chicken. R Package Version 3.8.2. 2019. Available online: https://doi.org/10.18129/B9.bioc.org.Gg.eg.db (accessed on 3 August 2024).

30. Pagès, H.; Carlson, M.; Falcon, S.; Li, N. AnnotationDbi: Manipulation of SQLite-Based Annotations in Bioconductor. 2023. Available online: https://doi.org/10.18129/B9.bioc AnnotationDbi (accessed on 3 August 2024).

31. Gentleman, R.C.; Carey, V.J.; Bates, D.M.; Bolstad, B.; Dettling, M.; Dudoit, S.; Ellis, B.; Gautier, L.; Ge, Y.; Gentry, J.; et al. Bioconductor: Open Software Development for Computational Biology and Bioinformatics. *Genome Biol.* **2004**, *5*, R80. [CrossRef]

32. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. ClusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. *Innovation* **2021**, *2*, 100141. [CrossRef]

33. Liberzon, A.; Subramanian, A.; Pinchback, R.; Thorvaldsdóttir, H.; Tamayo, P.; Mesirov, J.P. Molecular Signatures Database (MSigDB) 3.0. *Bioinformatics* **2011**, *27*, 1739–1740. [CrossRef]

34. 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. [CrossRef]

35. Hu, Z.L.; Park, C.A.; Reecy, J.M. Bringing the Animal QTLdb and CorrDB into the Future: Meeting New Challenges and Providing Updated Services. *Nucleic Acids Res.* **2022**, *50*, D956–D961. [CrossRef] [PubMed]

36. 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. [CrossRef] [PubMed]

37. Nascimento, C.S.; Barbosa, L.T.; Brito, C.; Fernandes, R.P.M.; Mann, R.S.; Pinto, A.P.G.; Oliveira, H.C.; Dodson, M.V.; Guimarães, S.E.F.; Duarte, M.S. Identification of Suitable Reference Genes for Real Time Quantitative Polymerase Chain Reaction Assays on Pectoralis Major Muscle in Chicken (*Gallus gallus*). *PLoS ONE* **2015**, *10*, e0127935. [CrossRef]

38. 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. [CrossRef]

39. Rausch, T.; Hsi-Yang Fritz, M.; Korbel, J.O.; Benes, V. Alfred: Interactive Multi-Sample BAM Alignment Statistics, Feature Counting and Feature Annotation for Long- and Short-Read Sequencing. *Bioinformatics* **2019**, *35*, 2489–2491. [CrossRef]

40. Rausch, T.; Fritz, M.H.Y.; Untergasser, A.; Benes, V. Tracy: Basecalling, Alignment, Assembly and Deconvolution of Sanger Chromatogram Trace Files. *BMC Genom.* **2020**, *21*, 230. [CrossRef]

41. Untergasser, A.; Ruijter, J.M.; Benes, V.; van den Hoff, M.J.B. Web-Based LinRegPCR: Application for the Visualization and Analysis of (RT)-QPCR Amplification and Melting Data. *BMC Bioinform.* **2021**, *22*, 398. [CrossRef]

42. GEAR: Genome Analysis Server. Available online: https://www.gear-genomics.com/ (accessed on 3 August 2024).

43. EndoGenes—Análise de Genes de Referência. Available online: https://hanielcedraz.shinyapps.io/endoGenes/ (accessed on 3 August 2024).

44. 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. [CrossRef]

45. Início Gene. NCBI. Available online: https://www.ncbi.nlm.nih.gov/gene (accessed on 3 August 2024).

46. Navegador do Genoma Ensembl 112. Available online: https://www.ensembl.org/index.html (accessed on 3 August 2024).

47. Pfaffl, M.W. A New Mathematical Model for Relative Quantification in Real-Time RT-PCR. *Nucleic Acids Res.* **2001**, *29*, e45. [CrossRef]

48. Zhang, J.; Kaiser, M.G.; Deist, M.S.; Gallardo, R.A.; Bunn, D.A.; Kelly, T.R.; Dekkers, J.C.M.; Zhou, H.; Lamont, S.J. Transcriptome Analysis in Spleen Reveals Differential Regulation of Response to Newcastle Disease Virus in Two Chicken Lines. *Sci. Rep.* **2018**, *8*, 1278. [CrossRef] [PubMed]

49. Mutryn, M.F.; Brannick, E.M.; Fu, W.; Lee, W.R.; Abasht, B. Characterization of a Novel Chicken Muscle Disorder through Differential Gene Expression and Pathway Analysis Using RNA-Sequencing. *BMC Genom.* **2015**, *16*, 399. [CrossRef]

50. Zambonelli, P.; Zappaterra, M.; Soglia, F.; Petracci, M.; Sirri, F.; Cavani, C.; Davoli, R. Detection of Differentially Expressed Genes in Broiler Pectoralis Major Muscle Affected by White Striping—Wooden Breast Myopathies. *Poult. Sci.* **2016**, *95*, 2771–2785. [CrossRef] [PubMed]

51. Kong, B.W.; Hudson, N.; Seo, D.; Lee, S.; Khatri, B.; Lassiter, K.; Cook, D.; Piekarski, A.; Dridi, S.; Anthony, N.; et al. RNA Sequencing for Global Gene Expression Associated with Muscle Growth in a Single Male Modern Broiler Line Compared to a Foundational Barred Plymouth Rock Chicken Line. *BMC Genom.* **2017**, *18*, 82. [CrossRef]
52. Papah, M.B.; Abasht, B. Dysregulation of Lipid Metabolism and Appearance of Slow Myofiber-Specific Isoforms Accompany the Development of Wooden Breast Myopathy in Modern Broiler Chickens. *Sci. Rep.* **2019**, *9*, 17170. [CrossRef]

53. Zhang, Z.; Du, H.; Yang, C.; Li, Q.; Qiu, M.; Song, X.; Yu, C.; Jiang, X.; Liu, L.; Hu, C.; et al. Comparative Transcriptome Analysis Reveals Regulators Mediating Breast Muscle Growth and Development in Three Chicken Breeds. *Anim. Biotechnol.* **2019**, *30*, 233–241. [CrossRef] [PubMed]

54. McNamara, J.W.; Sadayappan, S. Skeletal Myosin Binding Protein-C: An Increasingly Important Regulator of Striated Muscle Physiology. *Arch. Biochem. Biophys.* **2018**, *660*, 121–128. [CrossRef]

55. Qaisar, R.; Bhaskaran, S.; Van Remmen, H. Muscle Fiber Type Diversification during Exercise and Regeneration. *Free Radic. Biol. Med.* **2016**, *98*, 56–67. [CrossRef]

56. Frey, N.; Barrientos, T.; Shelton, J.M.; Frank, D.; Rütten, H.; Gehring, D.; Kuhn, C.; Lutz, M.; Rothermel, B.; Bassel-Duby, R.; et al. Mice Lacking Calsarcin-1 Are Sensitized to Calcineurin Signaling and Show Accelerated Cardiomyopathy in Response to Pathological Biomechanical Stress. *Nat. Med.* **2004**, *10*, 1336–1343. [CrossRef]

57. Schulz, R.A.; Yutzey, K.E. Calcineurin Signaling and NFAT Activation in Cardiovascular and Skeletal Muscle Development. *Dev. Biol.* **2004**, *266*, 1–16. [CrossRef] [PubMed]

58. Wei, B.; Jin, J.P. TNNT1, TNNT2, and TNNT3: Isoform Genes, Regulation, and Structure–Function Relationships. *Gene* **2016**, *582*, 1–13. [CrossRef]

59. Yang, J.; Gong, Y.; Cai, J.; Liu, Q.; Zhang, Z. Lnc-3215 Suppression Leads to Calcium Overload in Selenium Deficiency-Induced Chicken Heart Lesion via the Lnc-3215-MiR-1594-TNN2 Pathway. *Mol. Ther. Nucleic Acids.* **2019**, *18*, 1–15. [CrossRef]

60. Cowling, B.S.; McGrath, M.J.; Nguyen, M.A.; Cottle, D.L.; Kee, A.J.; Brown, S.; Schessl, J.; Zou, Y.; Joya, J.; Bönnemann, C.G.; et al. Identification of FHL1 as a Regulator of Skeletal Muscle Mass: Implications for Human Myopathy. *J. Cell Biol.* **2008**, *183*, 1033–1048. [CrossRef] [PubMed]

61. Kong, B.; Owens, C.; Bottje, W.; Shakeri, M.; Choi, J.; Zhuang, H.; Bowker, B. Proteomic Analyses on Chicken Breast Meat with White Striping Myopathy. *Poult. Sci.* **2024**, *103*, 103682. [CrossRef]

62. Zhang, X.; Antonelo, D.; Hendrix, J.; To, V.; Campbell, Y.; Von Staden, M.; Li, S.; Suman, S.P.; Zhai, W.; Chen, J.; et al. Proteomic Characterization of Normal and Woody Breast Meat from Broilers of Five Genetic Strains. *Meat Muscle Biol.* **2020**, *4*, 1–17. [CrossRef]

63. Dowling, P.; Gargan, S.; Zweyer, M.; Sabir, H.; Swandulla, D.; Ohlendieck, K. Proteomic Profiling of Carbonic Anhydrase CA3 in Skeletal Muscle. *Expert Rev. Proteom.* **2021**, *18*, 1073–1086. [CrossRef]

64. Qian, S.H.; Chen, L.; Xiong, Y.L.; Chen, Z.X. Evolution and Function of Developmentally Dynamic Pseudogenes in Mammals. *Genome Biol.* **2022**, *23*, 235. [CrossRef] [PubMed]

65. Wang, C.; Zhang, B.; Ratliff, A.C.; Arlington, J.; Chen, J.; Xiong, Y.; Yue, F.; Nie, Y.; Hu, K.; Jin, W.; et al. Methyltransferase-like 21e Inhibits 26S Proteasome Activity to Facilitate Hypertrophy of Type IIb Myofibers. *FASEB J.* **2019**, *33*, 9672. [CrossRef]

66. Bordini, M.; Wang, Z.; Soglia, F.; Petracci, M.; Schmidt, C.J.; Abasht, B. RNA-Sequencing Revisited Data Shed New Light on Wooden Breast Myopathy. *Poult. Sci.* **2024**, *103*, 103902. [CrossRef]

67. Oehrl, W.; Cotsiki, M.; Panayotou, G. Differential Regulation of M3/6 (DUSP8) Signaling Complexes in Response to Arsenite-Induced Oxidative Stress. *Cell Signal.* **2013**, *25*, 429–438. [CrossRef]

68. Boyer, J.G.; Prasad, V.; Song, T.; Lee, D.; Fu, X.; Grimes, K.M.; Sargent, M.A.; Sadayappan, S.; Molkentin, J.D. ERK1/2 Signaling Induces Skeletal Muscle Slow Fiber-Type Switching and Reduces Muscular Dystrophy Disease Severity. *JCI Insight* **2019**, *4*, e127356. [CrossRef]

69. Almasi, S.; Crawford Parks, T.E.; Ravel-Chapuis, A.; MacKenzie, A.; Côté, J.; Cowan, K.N.; Jasmin, B.J. Differential Regulation of Autophagy by STAU1 in Alveolar Rhabdomyosarcoma and Non-transformed Skeletal Muscle Cells. *Cell. Oncol.* **2021**, *44*, 851–870. [CrossRef]

70. Zhang, X.; Nan, H.; Li, S. Dual Specificity Phosphatases 8 Directly Interacts with Mitogen-Activated Protein Kinase 1 to Regulate the Proliferation and Apoptosis of Rheumatoid Arthritis Fibroblast-like Synovial Cells. *Chin. J. Rheumatol.* **2022**, *433–438*, C7-1. [CrossRef]

71. Xie, P.; Guo, S.; Fan, Y.; Zhang, H.; Gu, D.; Li, H. Atrogin-1/MAFbx Enhances Simulated Ischemia/Reperfusion-Induced Apoptosis in Cardiomyocytes through Degradation of MAPK Phosphatase-1 and Sustained JNK Activation. *J. Biol. Chem.* **2009**, *284*, 5488–5496. [CrossRef] [PubMed]

72. Tintignac, L.A.; Lagirand, J.; Batonnet, S.; Sirri, V.; Leibovitch, M.P.; Leibovitch, S.A. Degradation of MyoD Mediated by the SCF (MAFbx) Ubiquitin Ligase. *J. Biol. Chem.* **2005**, *280*, 2847–2856. [CrossRef]

73. Li, J.; Hu, Y.; Lan, H.; Li, L.; Hu, X.; Li, N. P3014 The Study on the Genetic Mechanism of Varied Atrogin-1 Expression in Different Chicken Lines. *J. Anim. Sci.* **2016**, *94*, 58–59. [CrossRef]

74. Balamurugan, K.; Sterneck, E. The Many Faces of C/EBPδ and Their Relevance for Inflammation and Cancer. *Int. J. Biol. Sci.* **2013**, *9*, 917–933. [CrossRef]

75. Yu, Y.J.; Xu, Y.Y.; Lan, X.O.; Liu, X.Y.; Zhang, X.L.; Gao, X.H.; Geng, L. Shikonin Induces Apoptosis and Suppresses Growth in Keratinocytes via CEBP-δ Upregulation. *Int. Immunopharmacol.* **2019**, *72*, 511–521. [CrossRef]

76. Mao, X.G.; Xue, X.Y.; Lv, R.; Ji, A.; Shi, T.Y.; Chen, X.Y.; Jiang, X.F.; Zhang, X. CEBPD Is a Master Transcriptional Factor for Hypoxia Regulated Proteins in Glioblastoma and Augments Hypoxia Induced Invasion through Extracellular Matrix-Integrin Mediated EGFR/PI3K Pathway. *Cell Death Dis.* **2023**, *14*, 269. [CrossRef]
77. Ko, C.Y.; Chang, W.C.; Wang, J.M. Biological Roles of CCAAT/Enhancer-Binding Protein Delta during Inflammation. *J. Biomed. Sci.* **2015**, *22*, 6. [CrossRef] [PubMed]

78. Bello, N.F.; Lamsoul, I.; Heuzé, M.L.; Métais, A.; Moreaux, G.; Calderwood, D.A.; Duprez, D.; Moog-Lutz, C.; Lutz, P.G. The E3 Ubiquitin Ligase Specificity Subunit ASB2β Is a Novel Regulator of Muscle Differentiation That Targets Filamin B to Proteasomal Degradation. *Cell Death Differ.* **2009**, *16*, 921. [CrossRef]

79. Lu, P.Y.; Taylor, M.; Jia, H.T.; Ni, J.H. Muscle LIM Protein Promotes Expression of the Acetylcholine Receptor γ-Subunit Gene Cooperatively with the Myogenin-E12 Complex. *Cell. Mol. Life Sci.* **2004**, *61*, 2386-2392. [CrossRef] [PubMed]

80. Liu, J.; Liu, D.; Zhang, X.; Li, Y.; Fu, X.; He, W.; Li, M.; Chen, P.; Zeng, G.; DiSanto, M.E.; et al. NELL2 Modulates Cell Proliferation and Apoptosis via ERK Pathway in the Development of Benign Prostatic Hyperplasia. *Clin. Sci.* **2021**, *135*, 1591-1608. [CrossRef] [PubMed]

81. Gordon, B.S.; Rossetti, M.L.; Kimball, S.R. Alpha Arrestin Domain Containing 2 and Alpha Arrestin Domain Containing 3 Expression Is Reduced in Skeletal Muscle by Anabolic Stimuli and Increased by Catabolic Stimuli. *FASEB J.* **2019**, *33*, 700.10. [CrossRef]

82. Xia, H.; Zhang, Q.; Shen, Y.; Bai, Y.; Ma, X.; Zhang, B.; Qi, Y.; Zhang, J.; Hu, Q.; Du, W.; et al. Ube3d, a New Gene Associated with Age-Related Macular Degeneration, Induces Functional Changes in Both In Vivo and In Vitro Studies. *Mol. Ther. Nucleic Acids* **2020**, *20*, 217-230. [CrossRef] [PubMed]

83. Lilburn, M.S.; Griffin, J.R.; Wick, M. From Muscle to Food: Oxidative Challenges and Developmental Anomalies in Poultry Breast Muscle. *Poult. Sci.* **2019**, *98*, 4255-4260. [CrossRef] [PubMed]

84. Pette, D.; Staron, R.S. Transitions of Muscle Fiber Phenotypic Profiles. *Histochem. Cell Biol.* **2001**, *115*, 359-372. [CrossRef] [PubMed]

85. Abasht, B.; Papah, M.B.; Qiu, J. Evidence of Vascular Endothelial Dysfunction in Wooden Breast Disorder in Chickens: Insights through Gene Expression Analysis, Ultra-Structural Evaluation and Supervised Machine Learning Methods. *PLoS ONE* **2021**, *16*, e0243983. [CrossRef] [PubMed]

86. Lake, J.A.; Abasht, B. Glucolipotoxicity: A Proposed Etiology for Wooden Breast and Related Myopathies in Commercial Broiler Chickens. *Front. Physiol.* **2020**, *11*, 528245. [CrossRef] [PubMed]

87. Lee, J.H.; Kim, E.J.; Kim, D.K.; Lee, J.M.; Park, S.B.; Lee, I.K.; Harris, R.A.; Lee, M.O.; Choi, H.S. Hypoxia Induces PDK4 Gene Expression through Induction of the Orphan Nuclear Receptor ERRγ. *PLoS ONE* **2012**, *7*, e46324. [CrossRef]

88. Wang, X.; Shen, X.; Yan, Y.; Li, H. Pyruvate Dehydrogenase Kinases (PDKs): An Overview toward Clinical Applications. *Biosci. Rep.* **2021**, *41*, 20204402. [CrossRef]

89. Pettersen, I.K.N.; Tusubira, D.; Ashrafi, H.; Dyrstad, S.E.; Hansen, L.; Liu, X.Z.; Nilsson, L.I.H.; Løvsletten, N.G.; Berge, K.; Wergedahl, H.; et al. Upregulated PDK4 Expression Is a Sensitive Marker of Increased Fatty Acid Oxidation. *Mitochondrion* **2019**, *49*, 97-110. [CrossRef]

90. Kuttappan, V.A.; Brewer, V.B.; Mauromoustakos, A.; McKee, S.R.; Emmert, J.L.; Meullenet, J.F.; Owens, C.M. Estimation of Factors Associated with the Occurrence of White Striping in Broiler Breast Fillets. *Poult. Sci.* **2013**, *92*, 811-819. [CrossRef]

91. Koves, T.R.; Ussher, J.R.; Noland, R.C.; Slentz, D.; Mosedale, M.; Ilkayeva, O.; Bain, J.; Stevens, R.; Dyck, J.R.B.; Newgard, C.B.; et al. Mitochondrial Overload and Incomplete Fatty Acid Oxidation Contribute to Skeletal Muscle Insulin Resistance. *Cell Metab.* **2008**, *7*, 45-56. [CrossRef]

92. Bezaire, V.; Spriet, L.L.; Campbell, S.; Sabet, N.; Gerrits, M.; Bonen, A.; Harper, M.-E.; Harper, M.-E. Constitutive UCP3 Overexpression at Physiological Levels Increases Mouse Skeletal Muscle Capacity for Fatty Acid Transport and Oxidation. *FASEB J.* **2005**, *19*, 977-979. [CrossRef]

93. García-Martínez, C.; Sibille, B.; Solanes, G.; Darimont, C.; Macé, K.; Villarroya, F.; Gómez-Foix, A.M. Overexpression of UCP3 in Cultured Human Muscle Lowers Mitochondrial Membrane Potential, Raises ATP/ADP Ratio, and Favors Fatty Acid versus Glucose Oxidation. *FASEB J.* **2001**, *15*, 2033-2035. [CrossRef]

94. Schrauwen, P.; Hinderling, V.; Hesselink, M.K.C.; Schaart, G.; Kornips, E.; Saris, W.H.M.; Westerterp-Plantenga, M.; Langhans, W. Etomoxir-Induced Increase in UCP3 Supports a Role of Uncoupling Protein 3 as a Mitochondrial Fatty Acid Anion Exporter. *FASEB J.* **2002**, *16*, 1688-1690. [CrossRef] [PubMed]

95. Dobrzyn, A.; Dobrzyn, P.; Lee, S.H.; Miyazaki, M.; Cohen, P.; Asilmaz, E.; Hardie, D.G.; Friedman, J.M.; Ntambi, J.M. Stearoyl-CoA Desaturase-1 Deficiency Reduces Ceramide Synthesis by Downregulating Serine Palmitoyltransferase and Increasing β-Oxidation in Skeletal Muscle. *Am. J. Physiol. Endocrinol. Metab.* **2005**, *288*, 599-607. [CrossRef]

96. Bu, S.Y. Role of Acyl-CoA:Diacylglycerol Acyltransferase 2 (DGAT2) in Energy Metabolism in Skeletal Muscle Cells. *FASEB J.* **2015**, *29*, LB176. [CrossRef]

97. Bu, S.Y. Role of Dgat2 in Glucose Uptake and Fatty Acid Metabolism in C2C12 Skeletal Myotubes. *J. Microbiol. Biotechnol.* **2023**, *33*, 1563. [CrossRef]

98. Guerranti, R.; Pagani, R.; Neri, S.; Errico, S.V.; Leoncini, R.; Marinello, E. Inhibition and Regulation of Rat Liver L-Threonine Dehydrogenase by Different Fatty Acids and Their Derivatives. *Biochim. Biophys. Acta (BBA)-Gen. Subj.* **2001**, *1568*, 45-52. [CrossRef]

99. Li, Y.; Yuan, P.; Fan, S.; Zhai, B.; Jin, W.; Li, D.; Li, H.; Sun, G.; Han, R.; Liu, X.; et al. Weighted Gene Co-Expression Network Indicates That the DYNLL2 Is an Important Regulator of Chicken Breast Muscle Development and Is Regulated by MiR-148a-3p. *BMC Genom.* **2022**, *23*, 258. [CrossRef]
100. Li, F.; Luo, Z.; Huang, W.; Lu, Q.; Wilcox, C.S.; Jose, P.A.; Chen, S. Response Gene to Complement 32, a Novel Regulator for Transforming Growth Factor-β-Induced Smooth Muscle Differentiation of Neural Crest Cells. *J. Biol. Chem.* **2007**, *282*, 10133–10137.
[CrossRef]

101. Cheng, S.; Wan, X.; Yang, L.; Qin, Y.; Chen, S.; Liu, Y.; Sun, Y.; Qiu, Y.; Huang, L.; Qin, Q.; et al. RGCC-Mediated PLK1 Activity Drives Breast Cancer Lung Metastasis by Phosphorylating AMPKα2 to Activate Oxidative Phosphorylation and Fatty Acid Oxidation. *J. Exp. Clin. Cancer Res.* **2023**, *42*, 342. [CrossRef]

102. Tang, R.; Zhang, G.; Chen, S.Y. Response Gene to Complement 32 Protein Promotes Macrophage Phagocytosis via Activation of Protein Kinase C Pathway. *J. Biol. Chem.* **2014**, *289*, 22715–22722. [CrossRef]

103. Shoshani, T.; Faerman, A.; Mett, I.; Zelin, E.; Tenne, T.; Gorodin, S.; Moshel, Y.; Elbaz, S.; Budanov, A.; Chajut, A.; et al. Identification of a Novel Hypoxia-Inducible Factor 1-Responsive Gene, RTP801, Involved in Apoptosis. *Mol. Cell Biol.* **2002**, *22*, 2283. [CrossRef]

104. Brugarolas, J.; Lei, K.; Hurley, R.L.; Manning, B.D.; Reiling, J.H.; Hafen, E.; Witters, L.A.; Ellisen, L.W.; Kaelin, W.G. Regulation of MTOR Function in Response to Hypoxia by REDD1 and the TSC1/TSC2 Tumor Suppressor Complex. *Genes Dev.* **2004**, *18*, 2893. [CrossRef]

105. Laplante, M.; Sabatini, D.M. MTOR Signaling at a Glance. *J. Cell Sci.* **2009**, *122*, 3589–3594. [CrossRef]

106. Ayansola, H.; Liao, C.; Dong, Y.; Yu, X.; Zhang, B.; Wang, B. Prospect of Early Vascular Tone and Satellite Cell Modulations on White Striping Muscle Myopathy. *Poult. Sci.* **2021**, *100*, 100945. [CrossRef]

107. Obach, M.; Navarro-Sabaté, À.; Caro, J.; Kong, X.; Duran, J.; Gómez, M.; Perales, J.C.; Ventura, F.; Rosa, J.L.; Bartrons, R. 6-Phosphofructo-2-Kinase (Pfkfb3) Gene Promoter Contains Hypoxia-Inducible Factor-1 Binding Sites Necessary for Transactivation in Response to Hypoxia. *J. Biol. Chem.* **2004**, *279*, 53562–53570. [CrossRef]

108. Yalcin, A.; Clem, B.F.; Simmons, A.; Lane, A.; Nelson, K.; Clem, A.L.; Brock, E.; Siow, D.; Wattenberg, B.; Telang, S.; et al. Nuclear Targeting of 6-Phosphofructo-2-Kinase (PFKFB3) Increases Proliferation via Cyclin-Dependent Kinases. *J. Biol. Chem.* **2009**, *284*, 24223. [CrossRef]

109. Yalcin, A.; Clem, B.F.; Imbert-Fernandez, Y.; Ozcan, S.C.; Peker, S.; O'Neal, J.; Klarer, A.C.; Clem, A.L.; Telang, T.; Chesney, J. 6-Phosphofructo-2-Kinase (PFKFB3) Promotes Cell Cycle Progression and Suppresses Apoptosis via Cdk1-Mediated Phosphorylation of P27. *Cell Death Dis.* **2014**, *5*, e1337. [CrossRef]

110. Gassen, N.C.; Hartmann, J.; Schmidt, M.V.; Rein, T. FKBP5/FKBP51 Enhances Autophagy to Synergize with Antidepressant Action. *Autophagy* **2015**, *11*, 578–580. [CrossRef]

111. Chen, Y.; Liu, Z.; Wang, Y.; Zhuang, J.; Peng, Y.; Mo, X.; Chen, J.; Shi, Y.; Yu, M.; Cai, W.; et al. FKBP51 Induces P53-Dependent Apoptosis and Enhances Drug Sensitivity of Human Non-Small-Cell Lung Cancer Cells. *Exp. Ther. Med.* **2020**, *19*, 2236–2242. [CrossRef]

112. Gao, S.; Huang, S.; Zhang, Y.; Fang, G.; Liu, Y.; Zhang, C.; Li, Y.; Du, J. The Transcriptional Regulator KLF15 Is Necessary for Myoblast Differentiation and Muscle Regeneration by Activating FKBP5. *J. Biol. Chem.* **2023**, *299*, 105226. [CrossRef]

113. Lee, J.H.; Budanov, A.V.; Karin, M. Sestrins Orchestrate Cellular Metabolism to Attenuate Aging. *Cell Metab.* **2013**, *18*, 792–801. [CrossRef]

114. Peeters, H.; Debeer, P.; Bairoch, A.; Wilquet, V.; Huysmans, C.; Parthoens, E.; Fryns, J.P.; Gewillig, M.; Nakamura, Y.; Niikawa, N.; et al. PA26 Is a Candidate Gene for Heterotaxia in Humans: Identification of a Novel PA26-Related Gene Family in Human and Mouse. *Hum. Genet.* **2003**, *112*, 573–580. [CrossRef]

115. Cai, B.; Ma, M.; Chen, B.; Li, Z.; Abdalla, B.A.; Nie, Q.; Zhang, X. MiR-16-5p Targets SESN1 to Regulate the P53 Signaling Pathway, Affecting Myoblast Proliferation and Apoptosis, and Is Involved in Myoblast Differentiation. *Cell Death Dis.* **2018**, *9*, 367. [CrossRef]

116. Thomas, T.M.; Miyaguchi, K.; Edwards, L.A.; Wang, H.; Wollebo, H.; Aiguo, L.; Murali, R.; Wang, Y.; Braas, D.; Michael, J.S.; et al. Elevated Asparagine Biosynthesis Drives Brain Tumor Stem Cell Metabolic Plasticity and Resistance to Oxidative Stress. *Mol. Cancer Res.* **2021**, *19*, 1375–1388. [CrossRef]

117. Aung, H.H.; Lame, M.W.; Gohil, K.; An, C.I.; Wilson, D.W.; Rutledge, J.C. Induction of ATF3 Gene Network by Triglyceride-Rich Lipoprotein Lipolysis Products Increases Vascular Apoptosis and Inflammation. *Arterioscler. Thromb. Vasc. Biol.* **2013**, *33*, 2088–2096. [CrossRef] [PubMed]

118. Feng, J.; Li, Y.; Jin, X.; Gong, R.; Xia, Z. ATF3 Regulates Oxidative Stress and Extracellular Matrix Degradation via P38/Nrf2 Signaling Pathway in Pelvic Organ Prolapse. *Tissue Cell* **2021**, *73*, 101660. [CrossRef] [PubMed]

119. Kim, J.Y.; Park, K.J.; Hwang, J.Y.; Kim, G.H.; Lee, D.Y.; Lee, Y.J.; Song, E.H.; Yoo, M.G.; Kim, B.J.; Suh, Y.H.; et al. Activating Transcription Factor 3 Is a Target Molecule Linking Hepatic Steatosis to Impaired Glucose Homeostasis. *J. Hepatol.* **2017**, *67*, 349–359. [CrossRef]

120. Bordini, M.; Zappaterra, M.; Soglia, F.; Petracci, M.; Davoli, R. Weighted Gene Co-Expression Network Analysis Identifies Molecular Pathways and Hub Genes Involved in Broiler White Striping and Wooden Breast Myopathies. *Sci. Rep.* **2021**, *11*, 1776. [CrossRef]

121. Lee, B.; Kim, J.Y.; Ryu, Y.C.; Lee, K.; Choi, Y.M. Research Note: Expression Levels of Collagen-Related Genes in PSE Conditions and White Striping Features of Broiler Pectoralis Major Muscle. *Poult. Sci.* **2023**, *102*, 102471. [CrossRef]
122. Malila, Y.; U-Chupaj, J.; Srimarut, Y.; Chaiwiwattrakul, P.; Uengwetwanit, T.; Arayamethakorn, S.; Punyapornwithaya, V.; Sansamur, C.; Kirschke, C.P.; Huang, L.; et al. Monitoring of White Striping and Wooden Breast Cases and Impacts on Quality of Breast Meat Collected from Commercial Broilers (*Gallus gallus*). *Asian-Australas J. Anim. Sci.* **2018**, *31*, 1807. [CrossRef]

123. Petracci, M.; Mudalal, S.; Babini, E.; Cavani, C. Effect of White Striping on Chemical Composition and Nutritional Value of Chicken Breast Meat. *Ital. J. Anim. Sci.* **2014**, *13*, 179-183. [CrossRef]

124. Wang, D.; Mohammad, M.; Wang, Y.; Tan, R.; Murray, L.S.; Ricardo, S.; Dagher, H.; van Agtmael, T.; Savige, J. The Chemical Chaperone, PBA, Reduces ER Stress and Autophagy and Increases Collagen IV A5 Expression in Cultured Fibroblasts from Men With X-Linked Alport Syndrome and Missense Mutations. *Kidney Int. Rep.* **2017**, *2*, 739-748. [CrossRef]

125. Valentine, B.A. Skeletal Muscle. In *Pathologic Basis of Veterinary Disease*, 6th ed.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 908-953. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.