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.

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.

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.

Pathogenic variability among Pasteurella multocida type A isolates from Brazilian pig farms.

BACKGROUND: Pasteurella multocida type A (PmA) is considered a secondary agent of pneumonia in pigs. The role of PmA as a primary pathogen was investigated by challenging pigs with eight field strains isolated from pneumonia and serositis in six Brazilian states. Eight groups of eight pigs each were intranasally inoculated with different strains of PmA (1.5 mL/nostril of 10e7 CFU/mL). The control group (n = 12) received sterile PBS. The pigs were euthanized by electrocution and necropsied by 5 dpi. Macroscopic lesions were recorded, and swabs and fragments of thoracic and abdominal organs were analyzed by bacteriological and pathological assays. The PmA strains were analyzed for four virulence genes (toxA: toxin; pfhA: adhesion; tbpA and hgbB: iron acquisition) by PCR and sequencing and submitted to multilocus sequence typing (MLST). RESULTS: The eight PmA strains were classified as follows: five as highly pathogenic (HP) for causing necrotic bronchopneumonia and diffuse fibrinous pleuritis and pericarditis; one as low pathogenic for causing only focal bronchopneumonia; and two as nonpathogenic because they did not cause injury to any pig. PCR for the gene pfhA was positive for all five HP isolates. Sequencing demonstrated that the pfhA region of the HP strains comprised four genes: tpsB1, pfhA1, tpsB2 and pfhA2. The low and nonpathogenic strains did not contain the genes tpsB2 and pfhA2. A deletion of four bases was observed in the pfhA gene in the low pathogenic strain, and an insertion of 37 kb of phage DNA was observed in the nonpathogenic strains. MLST clustered the HP isolates in one group and the low and nonpathogenic isolates in another. Only the nonpathogenic isolates matched sequence type 10; the other isolates did not match any type available in the MLST database. CONCLUSIONS: The hypothesis that some PmA strains are primary pathogens and cause disease in pigs without any co-factor was confirmed. The pfhA region, comprising the genes tpsB1, tpsB2, pfhA1 and pfhA2, is related to the pathogenicity of PmA. The HP strains can cause necrotic bronchopneumonia, fibrinous pleuritis and pericarditis in pigs and can be identified by PCR amplification of the gene pfhA2.

Genetic characterization of porcine circovirus type 2 in captive wild boars in southern Brazil.

Porcine circovirus type 2 (PCV2) has been identified in pig population in Brazil since 2000, but scarce studies involving wild boars with PCV2 infection are reported in the country. This study aimed to perform the genetic characterization of PCV2 detected in clinically healthy captive wild boars from farms located in Southern Brazil. Bronchial and mesenteric lymph nodes from 129 clinically healthy captive wild boars were tested by nested PCR for PCV2 detection. Six out of 38 positive samples (29.5%) were submitted to a quantitative real time PCR (qPCR) and genetic sequencing. Viral load up to 1.19 × 109 viral DNA copies/uL was detected in lymph nodes samples by qPCR. According to the ORF2 gene sequence analysis, all PCV2 samples were classified into PCV2b genotype. Comparisons based on a 702 nt region of the ORF2 of all six isolates revealed a high degree of similarity between these isolates. The ORF2 sequences characterized here share 97.1-100% of nucleotide identity and 95.7-100% of amino acid identity with other PCV2b isolated in Brazil from wild boars and feral pigs. This study reports the first detection and genetic characterization of PCV2b in captive wild boars in Brazil and provides important information on PCV2 infection in this domesticated species.

Archaeal Community Changes Associated with Cultivation of Amazon Forest Soil with Oil Palm.

This study compared soil archaeal communities of the Amazon forest with that of an adjacent area under oil palm cultivation by 16S ribosomal RNA gene pyrosequencing. Species richness and diversity were greater in native forest soil than in the oil palm-cultivated area, and 130 OTUs (13.7%) were shared between these areas. Among the classified sequences, Thaumarchaeota were predominant in the native forest, whereas Euryarchaeota were predominant in the oil palm-cultivated area. Archaeal species diversity was 1.7 times higher in the native forest soil, according to the Simpson diversity index, and the Chao1 index showed that richness was five times higher in the native forest soil. A phylogenetic tree of unclassified Thaumarchaeota sequences showed that most of the OTUs belong to Miscellaneous Crenarchaeotic Group. Several archaeal genera involved in nutrient cycling (e.g., methanogens and ammonia oxidizers) were identified in both areas, but significant differences were found in the relative abundances of Candidatus Nitrososphaera and unclassified Soil Crenarchaeotic Group (prevalent in the native forest) and Candidatus Nitrosotalea and unclassified Terrestrial Group (prevalent in the oil palm-cultivated area). More studies are needed to culture some of these Archaea in the laboratory so that their metabolism and physiology can be studied.

Impact of long-term cropping of glyphosate-resistant transgenic soybean [Glycine max (L.) Merr.] on soil microbiome.

The transgenic soybean [Glycine max (L.) Merrill] occupies about 80 % of the global area cropped with this legume, the majority comprising the glyphosate-resistant trait (Roundup Ready(®), GR or RR). However, concerns about possible impacts of transgenic crops on soil microbial communities are often raised. We investigated soil chemical, physical and microbiological properties, and grain yields in long-term field trials involving conventional and nearly isogenic RR transgenic genotypes. The trials were performed at two locations in Brazil, with different edaphoclimatic conditions. Large differences in physical, chemical and classic microbiological parameters (microbial biomass of C and N, basal respiration), as well as in grain production were observed between the sites. Some phyla (Proteobacteria, Actinobacteria, Acidobacteria), classes (Alphaproteobacteria, Actinomycetales, Solibacteres) and orders (Rhizobiales, Burkholderiales, Myxococcales, Pseudomonadales), as well as some functional subsystems (clustering-based subsystems, carbohydrates, amino acids and protein metabolism) were, in general, abundant in all treatments. However, bioindicators related to superior soil fertility and physical properties at Londrina were identified, among them a higher ratio of Proteobacteria:Acidobacteria. Regarding the transgene, the metagenomics showed differences in microbial taxonomic and functional abundances, but lower in magnitude than differences observed between the sites. Besides the site-specific differences, Proteobacteria, Firmicutes and Chlorophyta were higher in the transgenic treatment, as well as sequences related to protein metabolism, cell division and cycle. Although confirming effects of the transgenic trait on soil microbiome, no differences were recorded in grain yields, probably due to the buffering capacity associated with the high taxonomic and functional microbial diversity observed in all treatments.

Influenza A Viruses of Human Origin in Swine, Brazil.

The evolutionary origins of the influenza A(H1N1)pdm09 virus that caused the first outbreak of the 2009 pandemic in Mexico remain unclear, highlighting the lack of swine surveillance in Latin American countries. Although Brazil has one of the largest swine populations in the world, influenza was not thought to be endemic in Brazil's swine until the major outbreaks of influenza A(H1N1)pdm09 in 2009. Through phylogenetic analysis of whole-genome sequences of influenza viruses of the H1N1, H1N2, and H3N2 subtypes collected in swine in Brazil during 2009-2012, we identified multiple previously uncharacterized influenza viruses of human seasonal H1N2 and H3N2 virus origin that have circulated undetected in swine for more than a decade. Viral diversity has further increased in Brazil through reassortment between co-circulating viruses, including A(H1N1)pdm09. The circulation of multiple divergent hemagglutinin lineages challenges the design of effective cross-protective vaccines and highlights the need for additional surveillance.

Assessment of bacterial and archaeal community structure in Swine wastewater treatment processes.

Microbial communities from two field-scale swine wastewater treatment plants (WWTPs) were assessed by pyrosequencing analyses of bacterial and archaeal 16S ribosomal DNA (rDNA) fragments. Effluent samples from secondary (anaerobic covered lagoons and upflow anaerobic sludge blanket [UASB]) and tertiary treatment systems (open-pond natural attenuation lagoon and air-sparged nitrification-denitrification tank followed by alkaline phosphorus precipitation process) were analyzed. A total of 56,807 and 48,859 high-quality reads were obtained from bacterial and archaeal libraries, respectively. Dominant bacterial communities were associated with the phylum Firmicutes, Bacteroidetes, Proteobacteria, or Actinobacteria. Bacteria and archaea diversity were highest in UASB effluent sample. Escherichia, Lactobacillus, Bacteroides, and/or Prevotella were used as indicators of putative pathogen reduction throughout the WWTPs. Satisfactory pathogen reduction was observed after the open-pond natural attenuation lagoon but not after the air-sparged nitrification/denitrification followed by alkaline phosphorus precipitation treatment processes. Among the archaeal communities, 80% of the reads was related to hydrogeno-trophic methanogens Methanospirillum. Enrichment of hydrogenotrophic methanogens detected in effluent samples from the anaerobic covered lagoons and UASB suggested that CO2 reduction with H2 was the dominant methanogenic pathway in these systems. Overall, the results served to improve our current understanding of major microbial communities' changes downgradient from the pen and throughout swine WWTP as a result of different treatment processes.

A human-like H1N2 influenza virus detected during an outbreak of acute respiratory disease in swine in Brazil.

Passive monitoring for detection of influenza A viruses (IAVs) in pigs has been carried out in Brazil since 2009, detecting mostly the A(H1N1)pdm09 influenza virus. Since then, outbreaks of acute respiratory disease suggestive of influenza A virus infection have been observed frequently in Brazilian pig herds. During a 2010-2011 influenza monitoring, a novel H1N2 influenza virus was detected in nursery pigs showing respiratory signs. The pathologic changes were cranioventral acute necrotizing bronchiolitis to subacute proliferative and purulent bronchointerstitial pneumonia. Lung tissue samples were positive for both influenza A virus and A(H1N1)pdm09 influenza virus based on RT-qPCR of the matrix gene. Two IAVs were isolated in SPF chicken eggs. HI analysis of both swine H1N2 influenza viruses showed reactivity to the H1δ cluster. DNA sequencing was performed for all eight viral gene segments of two virus isolates. According to the phylogenetic analysis, the HA and NA genes clustered with influenza viruses of the human lineage (H1-δ cluster, N2), whereas the six internal gene segments clustered with the A(H1N1)pdm09 group. This is the first report of a reassortant human-like H1N2 influenza virus derived from pandemic H1N1 virus causing an outbreak of respiratory disease in pigs in Brazil. The emergence of a reassortant IAV demands the close monitoring of pigs through the full-genome sequencing of virus isolates in order to enhance genetic information about IAVs circulating in pigs.

Wild-type lytic bacteriophages against Salmonella Heidelberg: Further characterization and effect of prophylactic therapy in broiler chickens.

To characterize wild-type bacteriophages and their effect on Salmonella Heidelberg intestinal colonization in broilers, phages combined in a cocktail were continuously delivered via drinking water since the first day after hatching. After challenge with a field strain, broilers were evaluated at regular intervals for S. Heidelberg and bacteriophages in tissues and cecum, and gross and microscopic lesions in organs. Phages were highly virulent against S. Heidelberg by efficiency of plating. One-step growth curves exhibited eclipse period from 20 to 25 min, whereas the lowest latent period and higher burst size found were 45 min and 54 PFU/cell, respectively. Bacteriophage whole genomic sequencing analyses revealed a lack of genes related to lysogeny, antimicrobial resistance, and virulence factors. Relevant gross or microscopic lesions were absent in tissues analyzed from treated broilers. Although numerically stable bacteriophage concentrations were detected in the cecal contents of treated broilers, no significant difference was found for the S. Heidelberg cecal load in comparison to the untreated group and for the prevalence of positive tissues throughout the evaluated period. The phages produced turbid plaques against some S. Heidelberg re-isolated from treated broilers, suggesting the evolving of a resistant subpopulation. Overall, the results provide new evidence of the safety and in vitro replication of such phages in S. Heidelberg. Nevertheless, continuous administration of the phage suspension most likely induced the development of bacteriophage-resistant mutants, which might have affected the in vivo effect. Therefore, a putative administration protocol should be based on other strategies, such as short-term therapy at pre-harvest age.

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.

Analysis of composition-based metagenomic classification.

BACKGROUND: An essential step of a metagenomic study is the taxonomic classification, that is, the identification of the taxonomic lineage of the organisms in a given sample. The taxonomic classification process involves a series of decisions. Currently, in the context of metagenomics, such decisions are usually based on empirical studies that consider one specific type of classifier. In this study we propose a general framework for analyzing the impact that several decisions can have on the classification problem. Instead of focusing on any specific classifier, we define a generic score function that provides a measure of the difficulty of the classification task. Using this framework, we analyze the impact of the following parameters on the taxonomic classification problem: (i) the length of n-mers used to encode the metagenomic sequences, (ii) the similarity measure used to compare sequences, and (iii) the type of taxonomic classification, which can be conventional or hierarchical, depending on whether the classification process occurs in a single shot or in several steps according to the taxonomic tree. RESULTS: We defined a score function that measures the degree of separability of the taxonomic classes under a given configuration induced by the parameters above. We conducted an extensive computational experiment and found out that reasonable values for the parameters of interest could be (i) intermediate values of n, the length of the n-mers; (ii) any similarity measure, because all of them resulted in similar scores; and (iii) the hierarchical strategy, which performed better in all of the cases. CONCLUSIONS: As expected, short n-mers generate lower configuration scores because they give rise to frequency vectors that represent distinct sequences in a similar way. On the other hand, large values for n result in sparse frequency vectors that represent differently metagenomic fragments that are in fact similar, also leading to low configuration scores. Regarding the similarity measure, in contrast to our expectations, the variation of the measures did not change the configuration scores significantly. Finally, the hierarchical strategy was more effective than the conventional strategy, which suggests that, instead of using a single classifier, one should adopt multiple classifiers organized as a hierarchy.

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.

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.

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.

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.