Phenotypic divergence between broiler and layer chicken lines is regulated at the molecular level during development.

BACKGROUND: Understanding the molecular underpinnings of phenotypic variations is critical for enhancing poultry breeding programs. The Brazilian broiler (TT) and laying hen (CC) lines exhibit striking differences in body weight, growth potential, and muscle mass. Our work aimed to compare the global transcriptome of wing and pectoral tissues during the early development (days 2.5 to 3.5) of these chicken lines, unveiling disparities in gene expression and regulation. RESULTS: Different and bona-fide transcriptomic profiles were identified for the compared lines. A similar number of up- and downregulated differentially expressed genes (DEGs) were identified, considering the broiler line as a reference. Upregulated DEGs displayed an enrichment of protease-encoding genes, whereas downregulated DEGs exhibited a prevalence of receptors and ligands. Gene Ontology analysis revealed that upregulated DEGs were mainly associated with hormone response, mitotic cell cycle, and different metabolic and biosynthetic processes. In contrast, downregulated DEGs were primarily linked to communication, signal transduction, cell differentiation, and nervous system development. Regulatory networks were constructed for the mitotic cell cycle and cell differentiation biological processes, as their contrasting roles may impact the development of distinct postnatal traits. Within the mitotic cell cycle network, key upregulated DEGs included CCND1 and HSP90, with central regulators being NF-κB subunits (RELA and REL) and NFATC2. The cell differentiation network comprises numerous DEGs encoding transcription factors (e.g., HOX genes), receptors, ligands, and histones, while the main regulatory hubs are CREB, AR and epigenetic modifiers. Clustering analyses highlighted PIK3CD as a central player within the differentiation network. CONCLUSIONS: Our study revealed distinct developmental transcriptomes between Brazilian broiler and layer lines. The gene expression profile of broiler embryos seems to favour increased cell proliferation and delayed differentiation, which may contribute to the subsequent enlargement of pectoral tissues during foetal and postnatal development. Our findings pave the way for future functional studies and improvement of targeted traits of economic interest in poultry.

Differential expression of miRNAs associated with pectoral myopathies in young broilers: insights from a comparative transcriptome analysis.

INTRODUCTION: White Striping (WS) and Wooden Breast (WB) pectoral myopathies are relevant disorders for contemporary broiler production worldwide. Several studies aimed to elucidate the genetic components associated with the occurrence of these myopathies. However, epigenetic factors that trigger or differentiate these two conditions are still unclear. The aim of this study was to identify miRNAs differentially expressed (DE) between normal and WS and WB-affected broilers, and to verify the possible role of these miRNAs in metabolic pathways related to the manifestation of these pectoral myopathies in 28-day-old broilers. RESULTS: Five miRNAs were DE in the WS vs control (gga-miR-375, gga-miR-200b-3p, gga-miR-429-3p, gga-miR-1769-5p, gga-miR-200a-3p), 82 between WB vs control and 62 between WB vs WS. Several known miRNAs were associated with WB, such as gga-miR-155, gga-miR-146b, gga-miR-222, gga-miR-146-5p, gga-miR- 29, gga-miR-21-5p, gga-miR-133a-3p and gga-miR-133b. Most of them had not previously been associated with the development of this myopathy in broilers. We also have predicted 17 new miRNAs expressed in the broilers pectoral muscle. DE miRNA target gene ontology analysis enriched 6 common pathways for WS and WB compared to control: autophagy, insulin signaling, FoxO signaling, endocytosis, and metabolic pathways. The WS vs control contrast had two unique pathways, ERBB signaling and the mTOR signaling, while WB vs control had 14 unique pathways, with ubiquitin-mediated proteolysis and endoplasmic reticulum protein processing being the most significant. CONCLUSIONS: We found miRNAs DE between normal broilers and those affected with breast myopathies at 28 days of age. Our results also provide novel evidence of the miRNAs role on the regulation of WS and in the differentiation of both WS and WB myopathies. Overall, our study provides insights into miRNA-mediated and pathways involved in the occurrence of WS and WB helping to better understand these chicken growth disorders in an early age. These findings can help developing new approaches to reduce these complex issues in poultry production possibly by adjustments in nutrition and management conditions. Moreover, the miRNAs and target genes associated with the initial stages of WS and WB development could be potential biomarkers to be used in selection to reduce the occurrence of these myopathies in broiler production.

The downregulation of genes encoding muscle proteins have a potential role in the development of scrotal hernia in pigs.

BACKGROUND: Testicular descent is a physiological process regulated by many factors. Eventually, disturbances in the embryological/fetal development path facilitate the occurrence of scrotal hernia, a congenital malformation characterized by the presence of intestinal portions within the scrotal sac due to the abnormal expansion of the inguinal ring. In pigs, some genes have been related to this anomaly, but the genetic mechanisms involved remain unclear. This study aimed to investigate the expression profile of a set of genes potentially involved with the manifestation of scrotal hernia in the inguinal ring tissue. METHODS AND RESULTS: Tissue samples from the inguinal ring/canal of normal and scrotal hernia-affected male pigs with approximately 30 days of age were used. Relative expression analysis was performed using qPCR to confirm the expression profile of 17 candidate genes previously identified in an RNA-Seq study. Among them, the Myosin heavy chain 1 (MYH1), Desmin (DES), and Troponin 1 (TNNI1) genes were differentially expressed between groups and had reduced levels of expression in the affected animals. These genes encode proteins involved in the formation of muscle tissue, which seems to be important for increasing the resistance of the inguinal ring to the abdominal pressure, which is essential to avoid the occurrence of scrotal hernia. CONCLUSIONS: The downregulation of muscular candidate genes in the inguinal tissue clarifies the genetic mechanisms involved with this anomaly in its primary site, providing useful information for developing strategies to control this malformation in pigs and other mammals.

A joint analysis using exome and transcriptome data identifiescandidate polymorphisms and genes involved with umbilical hernia in pigs.

BACKGROUND: Umbilical Hernia (UH) is characterized by the passage of part of the intestine through the umbilical canal forming the herniary sac. There are several potential causes that can lead to the umbilical hernia such as bacterial infections, management conditions and genetic factors. Since the genetic components involved with UH are poorly understood, this study aimed to identify polymorphisms and genes associated with the manifestation of umbilical hernia in pigs using exome and transcriptome sequencing in a case and control design. RESULTS: In the exome sequencing, 119 variants located in 58 genes were identified differing between normal and UH-affected pigs, and in the umbilical ring transcriptome, 46 variants were identified, located in 27 genes. Comparing the two methodologies, we obtained 34 concordant variants between the exome and transcriptome analyses, which were located in 17 genes, distributed in 64 biological processes (BP). Among the BP involved with UH it is possible to highlight cell adhesion, cell junction regulation, embryonic morphogenesis, ion transport, muscle contraction, within others. CONCLUSIONS: We have generated the first exome sequencing related to normal and umbilical hernia-affected pigs, which allowed us to identify several variants possibly involved with this disorder. Many of those variants present in the DNA were confirmed with the RNA-Seq results. The combination of both exome and transcriptome sequencing approaches allowed us to better understand the complex molecular mechanisms underlying UH in pigs and possibly in other mammals, including humans. Some variants found in genes and other regulatory regions are highlighted as strong candidates to the development of UH in pigs and should be further investigated.

Genome-wide detection of CNVs and their association with performance traits in broilers.

BACKGROUND: Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age. RESULTS: We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels. CONCLUSIONS: Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers.

BACKGROUND: Bacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens. RESULTS: Our data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB. CONCLUSIONS: Our data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

Stable reference genes for expression studies in breast muscle of normal and white striping-affected chickens.

The normalization with proper reference genes is a crucial step to obtain accurate mRNA expression levels in quantitative PCR (qPCR) studies. Therefore, in this study, 10 reference candidate genes were evaluated to determine their stability in normal pectoralis major muscle of broilers and those counterparts affected with White Striping (WS) myopathy at 42 days age. Four different tools were used for ranking the most stable genes: GeNorm, NormFinder, BestKeeper and Comparative Ct (ΔCt), and a general ranking was performed using the RankAggreg tool to select the best reference genes among all tools. From the 10 genes evaluated in the breast muscle of broilers, 8 were amplified. Most of the algorithms/tools indicated the same two genes, RPL30 and RPL5, as the most stable in the broilers breast muscle. In addition, there was agreement among the tools for the least stable genes: MRPS27, GAPDH and RPLP1 in the broilers breast muscle. Therefore, it is interesting to note that even with different tools for evaluating gene expression, there was consensus on the most and least stable genes. These results indicate that the Ribosomal protein L30 (RPL30) and Ribosomal protein L5 (RPL5) can be recommended for accurate normalization in qPCR studies with chicken pectoralis major muscle affected with White Striping and other myopathies.

Genome-wide association scan for QTL and their positional candidate genes associated with internal organ traits in chickens.

BACKGROUND: Poultry breeding programs have been focused on improvement of growth and carcass traits, however, this has resulted in correlated changes in internal organ weights and increased incidence of metabolic disorders. These disorders can affect feed efficiency or even cause death. We used a high density SNP array (600 K, Affymetrix) to estimate genomic heritability, perform genome-wide association analysis, and identify genomic regions and positional candidate genes (PCGs) associated with internal organ traits in an F2 chicken population. We integrated knowledge of haplotype blocks, selection signature regions and sequencing data to refine the list of PCGs. RESULTS: Estimated genomic heritability for internal organ traits in chickens ranged from low (LUNGWT, 0.06) to high (GIZZWT, 0.45). A total of 20 unique 1 Mb windows identified on GGA1, 2, 4, 7, 12, 15, 18, 19, 21, 27 and 28 were significantly associated with intestine length, and weights or percentages of liver, gizzard or lungs. Within these windows, 14 PCGs were identified based on their biological functions: TNFSF11, GTF2F2, SPERT, KCTD4, HTR2A, RB1, PCDH7, LCORL, LDB2, NR4A2, GPD2, PTPN11, ITGB4 and SLC6A4. From those genes, two were located within haplotype blocks and three overlapped with selection signature regions. A total of 13,748 annotated sequence SNPs were in the 14 PCGs, including 156 SNPs in coding regions (124 synonymous, 26 non-synonymous, and 6 splice variants). Seven deleterious SNPs were identified in TNFSF11, NR4A2 or ITGB4 genes. CONCLUSIONS: The results from this study provide novel insights to understand the genetic architecture of internal organ traits in chickens. The QTL detection performed using a high density SNP array covered the whole genome allowing the discovery of novel QTL associated with organ traits. We identified PCGs within the QTL involved in biological processes that may regulate internal organ growth and development. Potential functional genetic variations were identified generating crucial information that, after validation, might be used in poultry breeding programs to reduce the occurrence of metabolic disorders.

Identification of selection signatures involved in performance traits in a paternal broiler line.

BACKGROUND: Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits. Selection signatures may be identified using different methodologies, of which some are based on detecting contiguous sequences of homozygous identical-by-descent haplotypes, called runs of homozygosity (ROH), or estimating fixation index (FST) of genomic windows that indicates genetic differentiation. This study aimed to identify selection signatures in a paternal broiler TT line at generations 7th and 16th of selection and to investigate the genes annotated in these regions as well as the biological pathways involved. For such purpose, ROH and FST-based analysis were performed using whole genome sequence of twenty-eight chickens from two different generations. RESULTS: ROH analysis identified homozygous regions of short and moderate size. Analysis of ROH patterns revealed regions commonly shared among animals and changes in ROH abundance and size between the two generations. Results also suggest that whole genome sequencing (WGS) outperforms SNPchip data avoiding overestimation of ROH size and underestimation of ROH number; however, sequencing costs can limited the number of animals analyzed. FST-based analysis revealed genetic differentiation in several genomic windows. Annotation of the consensus regions of ROH and FST windows revealed new and previously identified genes associated with traits of economic interest, such as APOB, IGF1, IGFBP2, POMC, PPARG, and ZNF423. Over-representation analysis of the genes resulted in biological terms of skeletal muscle, matrilin proteins, adipose tissue, hyperglycemia, diabetes, Salmonella infections and tyrosine. CONCLUSIONS: Identification of ROH and FST-based analyses revealed selection signatures in TT line and genes that have important role in traits of economic interest. Changes in the genome of the chickens were observed between the 7th and 16th generations showing that ancient and recent selection in TT line may have acted over genomic regions affecting diseases and performance traits.

Proximal femoral head transcriptome reveals novel candidate genes related to epiphysiolysis in broiler chickens.

BACKGROUND: The proximal femoral head separation (FHS) or epiphysiolysis is a prevalent disorder affecting the chicken femur epiphysis, being considered a risk factor to infection which can cause bacterial chondronecrosis with osteomyelitis in broilers. To identify the genetic mechanisms involved in epiphysiolysis, differentially expressed (DE) genes in the femur of normal and FHS-affected broilers were identified using RNA-Seq technology. Femoral growth plate (GP) samples from 35-day-old commercial male broilers were collected from 4 healthy and 4 FHS-affected broilers. Sequencing was performed using an Illumina paired-end protocol. Differentially expressed genes were obtained using the edgeR package based on the False Discovery Rate (FDR < 0.05). RESULTS: Approximately 16 million reads/sample were generated with 2 × 100 bp paired-end reads. After data quality control, approximately 12 million reads/sample were mapped to the reference chicken genome (Galgal5). A total of 12,645 genes were expressed in the femur GP. Out of those, 314 were DE between groups, being 154 upregulated and 160 downregulated in FHS-affected broilers. In the functional analyses, several biological processes (BP) were overrepresented. Among them, those related to cell adhesion, extracellular matrix (ECM), bone development, blood circulation and lipid metabolism, which are more related to chicken growth, are possibly involved with the onset of FHS. On the other hand, BP associated to apoptosis or cell death and immune response, which were also found in our study, could be related to the consequence of the FHS. CONCLUSIONS: Genes with potential role in the epiphysiolysis were identified through the femur head transcriptome analysis, providing a better understanding of the mechanisms that regulate bone development in fast-growing chickens. In this study, we highlighted the importance of cell adhesion and extracellular matrix related genes in triggering FHS. Furthermore, we have shown new insights on the involvement of lipidemia and immune response/inflammation with FHS in broilers. Understanding the changes in the GP transcriptome might support breeding strategies to address poultry robustness and to obtain more resilient broilers.

Unraveling genomic associations with feed efficiency and body weight traits in chickens through an integrative approach.

BACKGROUND: Feed efficiency and growth rate have been targets for selection to improve chicken production. The incorporation of genomic tools may help to accelerate selection. We genotyped 529 individuals using a high-density SNP chip (600 K, Affymetrix®) to estimate genomic heritability of performance traits and to identify genomic regions and their positional candidate genes associated with performance traits in a Brazilian F2 Chicken Resource population. Regions exhibiting selection signatures and a SNP dataset from resequencing were integrated with the genomic regions identified using the chip to refine the list of positional candidate genes and identify potential causative mutations. RESULTS: Feed intake (FI), feed conversion ratio (FC), feed efficiency (FE) and weight gain (WG) exhibited low genomic heritability values (i.e. from 0.0002 to 0.13), while body weight at hatch (BW1), 35 days-of-age (BW35), and 41 days-of-age (BW41) exhibited high genomic heritability values (i.e. from 0.60 to 0.73) in this F2 population. Twenty unique 1-Mb genomic windows were associated with BW1, BW35 or BW41, located on GGA1-4, 6-7, 10, 14, 24, 27 and 28. Thirty-eight positional candidate genes were identified within these windows, and three of them overlapped with selection signature regions. Thirteen predicted deleterious and three high impact sequence SNPs in these QTL regions were annotated in 11 positional candidate genes related to osteogenesis, skeletal muscle development, growth, energy metabolism and lipid metabolism, which may be associated with body weight in chickens. CONCLUSIONS: The use of a high-density SNP array to identify QTL which were integrated with whole genome sequence signatures of selection allowed the identification of candidate genes and candidate causal variants. One novel QTL was detected providing additional information to understand the genetic architecture of body weight traits. We identified QTL for body weight traits, which were also associated with fatness in the same population. Our findings form a basis for further functional studies to elucidate the role of specific genes in regulating body weight and fat deposition in chickens, generating useful information for poultry breeding programs.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines.

BACKGROUND: Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection. RESULTS: A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development. CONCLUSIONS: In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

A genome-wide association study reveals novel genomic regions and positional candidate genes for fat deposition in broiler chickens.

BACKGROUND: Excess fat content in chickens has a negative impact on poultry production. The discovery of QTL associated with fat deposition in the carcass allows the identification of positional candidate genes (PCGs) that might regulate fat deposition and be useful for selection against excess fat content in chicken's carcass. This study aimed to estimate genomic heritability coefficients and to identify QTLs and PCGs for abdominal fat (ABF) and skin (SKIN) traits in a broiler chicken population, originated from the White Plymouth Rock and White Cornish breeds. RESULTS: ABF and SKIN are moderately heritable traits in our broiler population with estimates ranging from 0.23 to 0.33. Using a high density SNP panel (355,027 informative SNPs), we detected nine unique QTLs that were associated with these fat traits. Among these, four QTL were novel, while five have been previously reported in the literature. Thirteen PCGs were identified that might regulate fat deposition in these QTL regions: JDP2, PLCG1, HNF4A, FITM2, ADIPOR1, PTPN11, MVK, APOA1, APOA4, APOA5, ENSGALG00000000477, ENSGALG00000000483, and ENSGALG00000005043. We used sequence information from founder animals to detect 4843 SNPs in the 13 PCGs. Among those, two were classified as potentially deleterious and two as high impact SNPs. CONCLUSIONS: This study generated novel results that can contribute to a better understanding of fat deposition in chickens. The use of high density array of SNPs increases genome coverage and improves QTL resolution than would have been achieved with low density. The identified PCGs were involved in many biological processes that regulate lipid storage. The SNPs identified in the PCGs, especially those predicted as potentially deleterious and high impact, may affect fat deposition. Validation should be undertaken before using these SNPs for selection against carcass fat accumulation and to improve feed efficiency in broiler chicken production.

CYP24A1 and TRPC3 Gene Expression in Kidneys and Their Involvement in Calcium and Phosphate Metabolism in Laying Hens.

Ca and P homeostasis across the egg-laying cycle is a complex process involving absorption in the small intestine, reabsorption/excretion in the kidneys, and eggshell gland secretion. Diets with inadequate calcium and phosphorus can interfere with their absorption and digestibility, resulting in eggshell quality losses and reduced productive life, affecting egg production and welfare. A better understanding of gene expression profiles in the kidneys of laying hens during the late egg-laying period could clarify the renal role in mineral metabolism at this late stage. Therefore, the performance, egg quality and bone integrity-related traits, and expression profiles of kidney candidate genes were evaluated in 73-week-old laying hens receiving different Ca and P ratios in their diet: a high Ca/P ratio (HR, 22.43), a low ratio (LR, 6.71), and a medium ratio (MR, 11.43). The laying hens receiving the HR diet had improved egg production and eggshell quality traits compared to the other two groups. Humerus length was shorter in the HR than in the other groups. The CYP24A1 and TRPC3 genes were differentially expressed (p.adj ≤ 0.05) among the groups. Therefore, their expression profiles could be involved in calcium and phosphate transcellular transport in 73-week-old laying hens as a way to keep mineral absorption at adequate levels.

Whole-Exome Sequencing Indicated New Candidate Genes Associated with Unilateral Cryptorchidism in Pigs.

INTRODUCTION: Cryptorchidism is a hereditary anomaly characterized by the incomplete descent of one or both testicles to the scrotum. One of the challenges of this anomaly is that the retained testicle maintains its endocrine function. As a consequence, cryptorchid animals produce hormone-tainted meat in comparison to castrated animals and are likely to be more aggressive. Cryptorchidism can lead to reduced animal welfare outcomes and cause economic losses. Identifying genetic markers for cryptorchidism is an essential step toward mitigating these negative outcomes and may facilitate genome manipulation to reduce the occurrence of cryptorchidism. Attempts to identify such markers have used genome-wide association studies. Using whole-exome sequencing, we aimed to identify single nucleotide polymorphisms (SNPs) in the coding regions of cryptorchid pigs and to characterize functional pathways concerning these SNPs. METHODS: DNA was extracted and sequenced from 5 healthy and 5 cryptorchid animals from the Landrace breed, using the Illumina HiSeq 2500 platform. Data were pre-processed using the SeqyClean tool and further mapped against the swine reference genome (Sus scrofa 11.1) using BWA software. GATK was used to identify polymorphisms (SNPs and InDels), which were annotated using the VEP tool. Network prediction and gene ontology enrichment analysis were conducted using the Cytoscape platform, and STRING software was used for visualization. RESULTS: A total of 63 SNPs were identified across the genes PIGB, CCPG1, COMMD9, LDLRAD3, TRIM44, MYLPF, SEPTIN, ZNF48, TIA1, FAIM2, KRT18, FBP1, FBP2, CTSL, DAPK1, DHX8, GPR179, DEPDC1B, ENSSSCG00000049573, ENSSSCG00000016384, ENSSSCG00000022657, ENSSSCG00000038825, and ENSSSCG00000001229. Using pathway enrichment analyses and network prospection, we have identified the following significant adjusted p value threshold of 0.001 involved with the biological function pathways of estrogen signaling, cytoskeleton organization, and the pentose phosphate pathway. CONCLUSION: Our data suggest the involvement of new SNPs and genes in developing cryptorchidism in pigs. However, further studies are needed to validate our results in a larger cohort population. Variations in the GPR179 gene, with implications at the protein level, may be associated with the appearance of this anomaly in the swine. Finally, we are showing that the estrogen signaling pathway may be involved in the pathophysiological mechanisms of this congenital anomaly as previously reported in GWAS.

Upregulated genes in articular cartilage may help to counteract femoral head separation in broilers with 21 days of age.

Femoral head necrosis (FHN) is one of the most common conditions in fast growing broilers, being characterized by separation of articular cartilage from epiphysis and classified as femoral head separation (FHS) or FHS with laceration (FHSL) depending on severity. Although molecular mechanisms involved with this disorder have been observed, its etiology is still unclear. Therefore, the expression of 15 candidate genes, chosen based on previous transcriptomic studies, was evaluated in the articular cartilage (AC) of normal and FHS-affected broilers at 21 days of age. Samples were collected based on the absence or presence of FHS for physical-chemical and qPCR analysis. The AvBD2, RHAG, COL28A1, ADA and ANGPTL7 were upregulated in FHS-affected broilers compared to the healthy group. These genes are involved in immune response, defense against pathogens, inflammation, cellular migration and adhesion, indicating different molecular mechanisms to control FHS progression at early age. Our results can contribute to improve the knowledge on FHN etiology in chickens and other species, such as horse and pigs that are severely affected by bone disorders.

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

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.

Downregulation of growth plate genes involved with the onset of femoral head separation in young broilers.

Femoral head separation (FHS) is characterized by the detachment of growth plate (GP) and articular cartilage, occurring in tibia and femur. However, the molecular mechanisms involved with this condition are not completely understood. Therefore, genes and biological processes (BP) involved with FHS were identified in 21-day-old broilers through RNA sequencing of the femoral GP. 13,487 genes were expressed in the chicken femoral head transcriptome of normal and FHS-affected broilers. From those, 34 were differentially expressed (DE; FDR ≤0.05) between groups, where all of them were downregulated in FHS-affected broilers. The main BP were enriched in receptor signaling pathways, ossification, bone mineralization and formation, skeletal morphogenesis, and vascularization. RNA-Seq datasets comparison of normal and FHS-affected broilers with 21, 35 and 42 days of age has shown three shared DE genes (FBN2, C1QTNF8, and XYLT1) in GP among ages. Twelve genes were exclusively DE at 21 days, where 10 have already been characterized (SHISA3, FNDC1, ANGPTL7, LEPR, ENSGALG00000049529, OXTR, ENSGALG00000045154, COL16A1, RASD2, BOC, GDF10, and THSD7B). Twelve SNPs were associated with FHS (p < 0.0001). Out of those, 5 were novel and 7 were existing variants located in 7 genes (RARS, TFPI2, TTI1, MAP4K3, LINK54, and AREL1). We have shown that genes related to chondrogenesis and bone differentiation were downregulated in the GP of FHS-affected young broilers. Therefore, these findings evince that candidate genes pointed out in our study are probably related to the onset of FHS in broilers.

SAP30 Gene Is a Probable Regulator of Muscle Hypertrophy in Chickens.

Animals with muscle hypertrophy phenotype are targeted by the broiler industry to increase the meat production and the quality of the final product. Studies characterizing the molecular machinery involved with these processes, such as quantitative trait loci studies, have been carried out identifying several candidate genes related to this trait; however, validation studies of these candidate genes in cell culture is scarce. The aim of this study was to evaluate SAP30 as a candidate gene for muscle development and to validate its function in cell culture in vitro. The SAP30 gene was downregulated in C2C12 muscle cell culture using siRNA technology to evaluate its impact on morphometric traits and gene expression by RNA-seq analysis. Modulation of SAP30 expression increased C2C12 myotube area, indicating a role in muscle hypertrophy. RNA-seq analysis identified several upregulated genes annotated in muscle development in treated cells (SAP30-knockdown), corroborating the role of SAP30 gene in muscle development regulation. Here, we provide experimental evidence of the involvement of SAP30 gene as a regulator of muscle cell hypertrophy.

A Missense Mutation in the MYBPH Gene Is Associated With Abdominal Fat Traits in Meat-Type Chickens.

Chicken is an important source of protein for human nutrition and a model system for growth and developmental biology. Although the genetic architecture of quantitative traits in meat-type chickens has been the subject of ongoing investigation, the identification of mutations associated with carcass traits of economic interest remains challenging. Therefore, our aim was to identify predicted deleterious mutation, which potentially affects protein function, and test if they were associated with carcass traits in chickens. For that, we performed a genome-wide association analysis (GWAS) for breast, thigh and drumstick traits in meat-type chickens and detected 19 unique quantitative trait loci (QTL). We then used: (1) the identified windows; (2) QTL for abdominal fat detected in a previous study with the same population and (3) previously obtained whole genome sequence data, to identify 18 predicted deleterious single nucleotide polymorphisms (SNPs) in those QTL for further association with breast, thigh, drumstick and abdominal fat traits. Using the additive model, a predicted deleterious SNP c.482C > T (SIFT score of 0.4) was associated (p-value < 0.05) with abdominal fat weight and percentage. This SNP is in the second exon of the MYBPH gene, and its allele frequency deviates from Hardy-Weinberg equilibrium. In conclusion, our study provides evidence that the c.482C > T SNP in the MYBPH gene is a putative causal mutation for fat deposition in meat-type chickens.