A small proportion of Zebu genetic background in crossbred calves may not be enough to improve resistance against natural bovine Babesia spp. infections.

The main objective of cattle breeders in tropical and subtropical regions is to acquire animals with taurine-productive traits adapted to the broad weather range of these regions. However, one of the main challenges on using taurine genetics in these areas is the high susceptibility of these animals to tick-borne diseases. Consequently, the present study evaluated from 10 November 2021-19 April 2022, the over 13 assessments, the Babesia bovis and Babesia bigemina DNA loads and the IgG anti-B. bovis and anti-B. bigemina levels in Angus (n = 17, 100% Taurine) and Ultrablack (n = 14, ∼82% taurine and 18% Zebu) calves. Data were analyzed using a multivariate mixed model with repeated measures of the same animal including the fixed effects of evaluation, genetic group, sex, Babesia spp., and their interactions. The repeatability values were estimated from the (co)variances matrix and expressed for each species. The correlations between the DNA loads (CNlog) and IgG titers (S/P) values for the two species were also estimated using the same model. Regarding the specific IgG antibody titers for both Babesia spp., no significant differences were observed between the two genetic groups. However, for B. bovis and B. bigemina DNA loads, Ultrablack calves presented significantly higher values than Angus calves. Under the conditions evaluated in this study, our findings suggest that the low percentage of Zebu genetic in the Ultrablack breed was insufficient to improve resistance against babesiosis. Further studies must demonstrate if the low percentages of Zebu genetics in Taurine breeds can modify the susceptibility to babesiosis infections.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep.

BACKGROUND: The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS: Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS: The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS: The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Detection and quantification of Babesia bovis and Babesia bigemina using different target genes.

Molecular assays have been widely used for the detection and quantification of bovine babesiosis due to their high sensitivity and specificity. However, variations in the sensitivity of pathogen detection may occur depending on the selected target gene. Thus, this study aimed to compare the detection sensitivity (DS) of Babesia bovis and B. bigemina infection levels in artificially and naturally infected cattle using quantitative PCR (qPCR) and six target genes. For B. bovis, the merozoite surface antigen genes 2b and 2c (msa-2b and msa-2c), and the mitochondrial cytochrome b gene (cybmt) were used. For B. bigemina, the genes encoding the proteins associated with rhoptry 1c (rap-1c), rap-1a, and cybmt were used. Six bovines, free of babesiosis, were artificially infected with 1 × 10-8 red blood cells infected (iRBC) with B. bovis (n = 3) or 1 × 10-6B. bigemina iRBC (n = 3). The animals were evaluated daily until parasitemia was confirmed (≥ 2.0%). The quantity of iRBC present in each animal was determined by examining blood smears. Blood samples were then subjected to DNA extraction, serial dilution, and qPCR analysis to determine the DS of each target gene. In addition, 30 calves naturally infected by Babesia spp. were also evaluated using the same six target genes. Regarding the artificial infection, B. bovis cybmt showed 25-fold higher sensitivity than the msa-2b and msa-2c genes, while the B. bigemina cybmt exhibited 5-fold and 25-fold higher sensitivity than the rap-1a and rap-1c genes, respectively. The rap-1a gene was found to be 5 times more sensitive than the rap-1c gene, while the B. bovis msa-2b and msa-2c genes exhibited similar DS. The positive frequencies of naturally infected calves for the target cybmt, msa-2b, and msa-2c genes (B. bovis) were: 100%, 33.3% and 50%, while cybmt, rap-1a, and rap-1c genes (B. bigemina) were 90%, 83.3%, and 63.3%, respectively. This study may contribute to the selection of suitable genes for molecular monitoring of bovine babesiosis. Mitochondrial genes could be considered as an alternative to improve the sensitivity of B. bovis and B. bigemina detection using qPCR.

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.

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.

Genome-Wide Association Study for Haemonchus contortus Resistance in Morada Nova Sheep.

Among the gastrointestinal nematodes affecting sheep, Haemonchus contortus is the most prevalent and virulent, resulting in health problems and production losses. Therefore, selecting sheep resistant to H. contortus is a suitable and sustainable strategy for controlling endoparasites in flocks. Here, 287 lambs of the native Brazilian Morada Nova hair sheep breed were subjected to two consecutive artificial infections with H. contortus and assessed for fecal egg count (FEC), packed cell volume (PCV), and live weight (LW). Forty-four animals ranked as having extreme resistance phenotypes were genotyped using the Illumina OvineSNP50v3 chip. A case−control genome-wide association study (GWAS) detected 37 significant (p < 0.001) markers in 12 ovine chromosomes in regions harboring quantitative trait loci (QTL) for FEC, Trichostrongylus spp. adults and larvae, weight, and fat; and candidate genes for immune responses, mucins, hematological parameters, homeostasis, and growth. Four single-nucleotide polymorphisms (SNP; OAR1_rs427671974, OAR2_rs419988472, OAR5_rs424070217, and OAR17_rs401006318) genotyped by qPCR followed by high-resolution melting (HRM) were associated with FEC and LW. Therefore, molecular markers detected by GWAS for H. contortus resistance in Morada Nova sheep may support animal selection programs aimed at controlling gastrointestinal nematode infections in flocks. Furthermore, genotyping of candidate genes using HRM qPCR may provide a rapid and efficient tool for animal identification.

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.

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.

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.

Differentially expressed genes in the femur cartilage transcriptome clarify the understanding of femoral head separation in chickens.

Locomotor problems are among one of the main concerns in the current poultry industry, causing major economic losses and affecting animal welfare. The most common bone anomalies in the femur are dyschondroplasia, femoral head separation (FHS), and bacterial chondronecrosis with osteomyelitis (BCO), also known as femoral head necrosis (FHN). The present study aimed to identify differentially expressed (DE) genes in the articular cartilage (AC) of normal and FHS-affected broilers by RNA-Seq analysis. In the transcriptome analysis, 12,169 genes were expressed in the femur AC. Of those, 107 genes were DE (FDR < 0.05) between normal and affected chickens, of which 9 were downregulated and 98 were upregulated in the affected broilers. In the gene-set enrichment analysis using the DE genes, 79 biological processes (BP) were identified and were grouped into 12 superclusters. The main BP found were involved in the response to biotic stimulus, gas transport, cellular activation, carbohydrate-derived catabolism, multi-organism regulation, immune system, muscle contraction, multi-organism process, cytolysis, leukocytes and cell adhesion. In this study, the first transcriptome analysis of the broilers femur articular cartilage was performed, and a set of candidate genes (AvBD1, AvBD2, ANK1, EPX, ADA, RHAG) that could trigger changes in the broiler´s femoral growth plate was identified. Moreover, these results could be helpful to better understand FHN in chickens and possibly in humans.

Differential Expression of Myogenic and Calcium Signaling-Related Genes in Broilers Affected With White Striping.

White Striping (WS) has been one of the main issues in poultry production in the last years since it affects meat quality. Studies have been conducted to understand WS and other myopathies in chickens, and some biological pathways have been associated to the prevalence of these conditions, such as extracellular calcium level, oxidative stress, localized hypoxia, possible fiber-type switching, and cellular repairing. Therefore, to understand the genetic mechanisms involved in WS, 15 functional candidate genes were chosen to be analyzed by quantitative PCR (qPCR) in breast muscle of normal and WS-affected chickens. To this, the pectoral major muscle (PMM) of 16 normal and 16 WS-affected broilers were collected at 42 days of age and submitted to qRT-PCR analysis. Out of the 15 genes studied, six were differentially expressed between groups. The CA2, CSRP3, and PLIN1 were upregulated, while CALM2, DNASE1L3, and MYLK2 genes were downregulated in the WS-affected when compared to the normal broilers. These findings highlight that the disruption on muscle and calcium signaling pathways can possibly be triggering WS in chickens. Improving our understanding on the genetic basis involved with this myopathy might contribute for reducing WS in poultry production.

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.

Exploring the genetic architecture of feed efficiency traits in chickens.

Chicken feed efficiency (FE) traits are the most important economic traits in broiler production. Several studies evaluating genetic factors affecting food consumption in chickens are available. However, most of these studies identified genomic regions containing putative quantitative trait loci for each trait separately. It is still a challenge to find common gene networks related to these traits. Therefore, here, a genome-wide association study (GWAS) was conducted to explore candidate genomic regions responsible for Feed Intake (FI), Body Weight Gain (BWG) and Feed Conversion Ratio (FCR) traits and their gene networks. A total of 1430 broilers from an experimental population was genotyped with the high density Affymetrix 600K SNP array. A total of 119 associated SNPs located in 20 chromosomes were identified, where some of them were common in more than one FE trait. In addition, novel genomic regions were prospected considering the SNPs dominance effects and sex interaction, identifying putative candidate genes only when these effects were fit in the model. Relevant candidate genes such as ATRNL1, PIK3C2A, PTPRN2, SORCS3 and gga-mir-1759 were highlighted in this study helping to elucidate the genomic architecture of feed efficiency traits. These results provide new insights on the mechanisms underlying the consumption and utilization of food in chickens.

Genes and SNPs Involved with Scrotal and Umbilical Hernia in Pigs.

Hernia is one of the most common defects in pigs. The most prevalent are the scrotal (SH), inguinal (IH) and umbilical (UH) hernias. We compared the inguinal ring transcriptome of normal and SH-affected pigs with the umbilical ring transcriptome of normal and UH-affected pigs to discover genes and pathways involved with the development of both types of hernia. A total of 13,307 transcripts was expressed in the inguinal and 13,302 in the umbilical ring tissues with 94.91% of them present in both tissues. From those, 35 genes were differentially expressed in both groups, participating in 108 biological processes. A total of 67 polymorphisms was identified in the inguinal ring and 76 in the umbilical ring tissue, of which 11 and 14 were novel, respectively. A single nucleotide polymorphism (SNP) with deleterious function was identified in the integrin α M (ITGAM) gene. The microtubule associated protein 1 light chain 3 γ (MAP1LC3C), vitrin (VIT), aggrecan (ACAN), alkaline ceramidase 2 (ACER2), potassium calcium-activated channel subfamily M α 1 (KCNMA1) and synaptopodin 2 (SYNPO2) genes are highlighted as candidates to trigger both types of hernia. We generated the first comparative study of the pig umbilical and inguinal ring transcriptomes, contributing to the understanding of the genetic mechanism involved with these two types of hernia in pigs and probably in other mammals.

Putative Epigenetic Biomarkers of Stress in Red Blood Cells of Chickens Reared Across Different Biomes.

Production animals are constantly subjected to early adverse environmental conditions that influence the adult phenotype and produce epigenetic effects. CpG dinucleotide methylation in red blood cells (RBC) could be a useful epigenetic biomarker to identify animals subjected to chronic stress in the production environment. Here we compared a reduced fraction of the RBC methylome of chickens exposed to social isolation to non-exposed. These experiments were performed in two different locations: Brazil and Sweden. The aim was to identify stress-associated DNA methylation profiles in RBC across these populations, in spite of the variable conditions to which birds are exposed in each facility and their different lineages. Birds were increasingly exposed to a social isolation treatment, combined with food and water deprivation, at random periods of the day from weeks 1-4 after hatching. We then collected the RBC DNA from individuals and compared a reduced fraction of their methylome between the experimental groups using two bioinformatic approaches to identify differentially methylated regions (DMRs): one using fixed-size windows and another that preselected differential peaks with MACS2. Three levels of significance were used (P ≤ 0.05, P ≤ 0.005, and P ≤ 0.0005) to identify DMRs between experimental groups, which were then used for different analyses. With both of the approaches more DMRs reached the defined significance thresholds in BR individuals compared to SW. However, more DMRs had higher fold change values in SW compared to BR individuals. Interestingly, ChrZ was enriched above expectancy for the presence of DMRs. Additionally, when analyzing the locations of these DMRs in relation to the transcription starting site (TSS), we found three peaks with high DMR presence: 10 kb upstream, the TSS itself, and 20-40 kb downstream. Interestingly, these peaks had DMRs with a high presence (>50%) of specific transcription factor binding sites. Three overlapping DMRs were found between the BR and SW population using the most relaxed p-value (P ≤ 0.05). With the most stringent p-value (P ≤ 0.0005), we found 7 and 4 DMRs between treatments in the BR and SW populations, respectively. This study is the first approximation to identify epigenetic biomarkers of long-term exposure to stress in different lineages of production animals.

Reference genes for proximal femoral epiphysiolysis expression studies in broilers cartilage.

The use of reference genes is required for relative quantification in gene expression analysis and the stability of these genes can be variable depending on the experimental design. Therefore, it is indispensable to test the reliability of endogenous genes previously to their use. This study evaluated nine candidate reference genes to select the most stable genes to be used as reference in gene expression studies with the femoral cartilage of normal and epiphysiolysis-affected broilers. The femur articular cartilage of 29 male broilers with 35 days of age was collected, frozen and further submitted to RNA extraction and quantitative PCR (qPCR) analysis. The candidate reference genes evaluated were GAPDH, HMBS, HPRT1, MRPS27, MRPS30, RPL30, RPL4, RPL5, and RPLP1. For the gene stability evaluation, three software were used: GeNorm, BestKeeper and NormFinder, and a global ranking was generated using the function RankAggreg. In this study, the RPLP1 and RPL5 were the most reliable endogenous genes being recommended for expression studies with femur cartilage in broilers with epiphysiolysis and possible other femur anomalies.

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.

Transcriptome analysis identifies genes involved with the development of umbilical hernias in pigs.

Umbilical hernia (UH) is one of the most frequent defects affecting pig production, however, it also affects humans and other mammals. UH is characterized as an abnormal protrusion of the abdominal contents to the umbilical region, causing pain, discomfort and reduced performance in pigs. Some genomic regions associated to UH have already been identified, however, no study involving RNA sequencing was performed when umbilical tissue is considered. Therefore, here, we have sequenced the umbilical ring transcriptome of five normal and five UH-affected pigs to uncover genes and pathways involved with UH development. A total of 13,216 transcripts were expressed in the umbilical ring tissue. From those, 230 genes were differentially expressed (DE) between normal and UH-affected pigs (FDR <0.05), being 145 downregulated and 85 upregulated in the affected compared to the normal pigs. A total of 68 significant biological processes were identified and the most relevant were extracellular matrix, immune system, anatomical development, cell adhesion, membrane components, receptor activation, calcium binding and immune synapse. The results pointed out ACAN, MMPs, COLs, EPYC, VIT, CCBE1 and LGALS3 as strong candidates to trigger umbilical hernias in pigs since they act in the extracellular matrix remodeling and in the production, integrity and resistance of the collagen. We have generated the first transcriptome of the pig umbilical ring tissue, which allowed the identification of genes that had not yet been related to umbilical hernias in pigs. Nevertheless, further studies are needed to identify the causal mutations, SNPs and CNVs in these genes to improve our understanding of the mechanisms of gene regulation.

Inguinal Ring RNA Sequencing Reveals Downregulation of Muscular Genes Related to Scrotal Hernia in Pigs.

Scrotal hernias (SH) are common congenital defects in commercial pigs, characterized by the presence of abdominal contents in the scrotal sac, leading to considerable production and animal welfare losses. Since the etiology of SH remains obscure, we aimed to identify the biological and genetic mechanisms involved in its occurrence through the whole transcriptome analysis of SH affected and unaffected pigs' inguinal rings. From the 22,452 genes annotated in the pig reference genome, 13,498 were expressed in the inguinal canal tissue. Of those, 703 genes were differentially expressed (DE, FDR < 0.05) between the two groups analyzed being, respectively, 209 genes upregulated and 494 downregulated in the SH-affected group. Thirty-seven significantly overrepresented GO terms related to SH were enriched, and the most relevant biological processes were muscular system, cell differentiation, sarcome reorganization, and myofibril assembly. The calcium signaling, hypertrophic cardiomyopathy, dilated cardiomyopathy, and cardiac muscle contraction were the major pathways possibly involved in the occurrence of the scrotal hernias. The expression profile of the DE genes was associated with the reduction of smooth muscle differentiation, followed by low calcium content in the cell, which could lead to a decreased apoptosis ratio and diminished muscle contraction of the inguinal canal region. We have demonstrated that genes involved with musculature are closely linked to the physiological imbalance predisposing to scrotal hernia. According to our study, the genes MYBPC1, BOK, SLC25A4, SLC8A3, DES, TPM2, MAP1CL3C, and FGF1 were considered strong candidates for future evaluation.