Whole transcriptome analysis of the pectoralis major muscle reveals molecular mechanisms involved with white striping in broiler chickens.

White striping (WS) is one of the most common myopathies identified in broiler chickens leading to substantial production losses, where the incidence reaches 12% in commercial chickens. It occurs primarily in heavier chickens being a modification of the breast muscle characterized by the presence of pale parallel streaks in the same orientation of the muscle fibers. Since the WS etiology remains unclear, we aimed to identify the biological and genetic mechanisms involved in its occurrence through the whole transcriptome analysis of WS in affected and unaffected chicken breast muscles. A total of 11,177 genes were expressed in the pectoralis major muscle. Out of those, 1,441 genes were differentially expressed (FDR ≤ 0.01) between the two analyzed groups, being, respectively, 772 genes upregulated and 669 downregulated in the WS affected group. A total of 36 significantly overrepresented GO terms related to WS myopathy were enriched, and the most relevant biological processes were activation of immune system, angiogenesis, hypoxia, cell death, and striated muscle contraction. The unbalance of those biological processes may trigger the occurrence of the WS phenotype in broilers. The possible lack of capillary blood supply homogeneously in the muscle triggers the hypoxia, following the activation of glycolysis, calcium signaling and apoptosis related genes facilitating the tissue damage and WS incidence.

Relationship of runs of homozygosity with adaptive and production traits in a paternal broiler line.

Genomic regions under high selective pressure present specific runs of homozygosity (ROH), which provide valuable information on the genetic mechanisms underlying the adaptation to environment imposed challenges. In broiler chickens, the adaptation to conventional production systems in tropical environments lead the animals with favorable genotypes to be naturally selected, increasing the frequency of these alleles in the next generations. In this study, ~1400 chickens from a paternal broiler line were genotyped with the 600 K Affymetrix® Axiom® high-density (HD) genotyping array for estimation of linkage disequilibrium (LD), effective population size (N e ), inbreeding and ROH. The average LD between adjacent single nucleotide polymorphisms (SNPs) in all autosomes was 0.37, and the LD decay was higher in microchromosomes followed by intermediate and macrochromosomes. The N e of the ancestral population was high and declined over time maintaining a sufficient number of animals to keep the inbreeding coefficient of this population at low levels. The ROH analysis revealed genomic regions that harbor genes associated with homeostasis maintenance and immune system mechanisms, which may have been selected in response to heat stress. Our results give a comprehensive insight into the relationship between shared ROH regions and putative regions related to survival and production traits in a paternal broiler line selected for over 20 years. These findings contribute to the understanding of the effects of environmental and artificial selection in shaping the distribution of functional variants in the chicken genome.

Preservation and reactivation of Candidatus Jettenia asiatica and Anammoxoglobus propionicus using different preservative agents.

Anaerobic ammonium oxidation (anammox) bacteria have peculiar characteristics that make them difficult to cultivate. The conservation of these microorganisms in culture collections or laboratories requires successful preservation and reactivation techniques. Furthermore, studies have shown that successful reactivation may be preservative dependent. Considering this, the present study aimed to evaluate the preservation and reactivation of anammox consortia enriched from swine manure treatment lagoons, by using different preservative agents at different temperatures: KNO3 (at 4 °C), glycerol (-20 °C, -80 °C), and skimmed cow milk (-20 °C, -80 °C, -200 °C). After 4 months, the biomass was thawed (except for KNO3), and the reestablishment of anammox activity was evaluated by stoichiometric coefficients. Microbial community transformation during the reactivation process was also studied by 16S rDNA sequence analysis. The results showed that the anammox biomass preserved with glycerol or skimmed cow milk at -80 °C recovered activity, while the biomass preserved with other methodologies did not reestablish activity during the studied time (90 days). The bacterial community from the biomass with anammox activity was characterized and showed the presence of Candidatus Brocadia anammoxidans, Candidatus Jettenia asiatica, and Candidatus Anammoxoglobus propionicus. Preservation with skimmed cow milk at -80 °C favored the selection of Candidatus Anammoxoglobus propionicus, while preservation with glycerol at -80 °C was successful for Candidatus Jettenia asiatica. The present study was effective on anammox sludge preservation and reactivation using low-cost processes for anammox cultures preservation, which is important for biomass transport and deammonification reactor start up.

Association of RUNX2 and TNFSF11 genes with production traits in a paternal broiler line.

Intense selection for production traits has improved the genetic gain of important economic traits. However, selection for performance and carcass traits has led to the onset of locomotors problems and decreasing bone strength in broilers. Thus, genes associated with bone integrity traits have become candidates for genetic studies in order to reduce the impact of bone disorders in broilers. This study investigated the association of the RUNX2 and TNFSF11 genes with 79 traits related to performance, carcass composition, organs, and bone integrity in a paternal broiler line. Analyses of genetic association between single-nucleotide polymorphisms (SNPs) and traits were carried out using the maximum likelihood procedures for mixed models. Genetic associations (P < 0.05) were found between SNP g.124,883A>G in the RUNX2 gene and chilled femur weight (additive plus dominance deviation effects within sex) and with performance traits (additive within sex and additive effects). The SNP g.14,862T>C in the TNFSF11 gene presented genetic associations (P < 0.05) with additive plus dominance deviation effects within sex for performance traits. Suggestive genetic associations (P < 0.10) were found with abdominal fat and its yield. Selection based on SNPs g.14,862T>C in TNFSF11 and g.124,883A>G in RUNX2 could be used to improve performance and carcass quality traits in the population studied, although SNP g.14,862T>C was not in Hardy-Weinberg equilibrium because it was not undergoing a selection process. Furthermore, it is important to validate these markers in an unrelated population for use in the selection process.

The involvement of RUNX2 and SPARC genes in the bacterial chondronecrosis with osteomyelitis in broilers.

Economic losses due to an increase of leg disorders in broilers have become a major concern of the poultry industry. Despite the efforts to reduce skeletal abnormalities in chickens, insufficient progress has been made. Bacterial chondronecrosis with osteomyelitis (BCO) is one of the main disorders that affect bone integrity in broilers. However, the genetic pathways and genes involved in most bone problems, including BCO, remains unclear. In this study, femoral samples from male broilers with 45 days of age affected or not with BCO were used to compare the relative expression with a reverse transcription real time PCR approach of 13 candidate genes: SPP1 (osteopontin), TNFRSF11B (osteoprotegerin), SPARC (osteonectin), CALB1 (calbidin 1), CALM (Calmodulin 2), IBSP (sialoprotein), COL1A2 (collagen, type I, α 2), BMP2 (bone morphogenetic protein 2), BMP3 (bone morphogenetic protein 3), RANKL (κ-B nuclear factor ligand), SMAD1 (SMAD family member 1), LEPR (leptin receptor) and RUNX2 (related transcription factor Runt 2). Differential expression test between affected and non-affected groups was performed using the REST software. The RUNX2 and SPARC genes were downregulated (P<0.05) in the affected group, with reduced expression of fourfold when compared with the non-affected group. This result indicates that the downregulation of RUNX2 and SPARC can contribute to an increased incidence of BCO in broilers.

Allele- and parent-of-origin-specific effects on expression of the KCNJ11 gene: A candidate for meat tenderness in cattle.

In contrast to the Mendelian inheritance model, parental alleles can contribute unequally to gene expression, which may result in phenotypic variance among individuals and bias in the predicted additive effect of molecular markers associated with production traits. Given the need to understand the effects of allelic variation and parent-of-origin effects on the expression of genes with a commercial interest in cattle, we analyzed the expression of KCNJ11 (potassium inwardly rectifying channel, subfamily J, member 11), which was previously described as a functional candidate gene for meat tenderness. Allele-specific and parent-of-origin-dependent expression of this gene were assessed in bovine muscle using the rs379610823 single nucleotide polymorphism as a reference. Biallelic expression was observed; however, the T allele was expressed at significantly higher levels than the C allele. Furthermore, increased expression of KCNJ11 was found in animals harboring the maternal T allele. This study is the first to describe the differential allelic expression of bovine KCNJ11. Our findings are important for understanding the mechanisms that underlie the pattern of KCNJ11 expression and its potential impact on the phenotypic variation of meat tenderness in Nelore beef cattle. This reinforces the need for further investigation of allelic- and parent-of-origin expression deviation in genetic markers eligible for the selection of target traits.

Genome-wide association study of periweaning failure-to-thrive syndrome (PFTS) in pigs.

Porcine periweaning-failure-to-thrive syndrome (PFTS) is a condition that affects newly weaned piglets. It is characterised by a progressive debilitation leading to death, in the absence of infectious, nutritional, management or environmental factors. In this study, we present the first report of PFTS in South America and the results of a genome-wide association study to identify the genetic markers associated with the appearance of this condition in a crossbred swine population. Four chromosomal regions were associated with PFTS predisposition, one located on SSCX, one on SSC8, and the two other regions on SSC14. Regions on SSC8 and SSC14 harbour important functional candidate genes involved in human depression and might have an important role in PFTS. Our findings contribute to the increasing knowledge about this syndrome, which has been investigated since 2007, and to the identification of the aetiology of this disease.

Babesia bovis and Babesia bigemina infection levels estimated by qPCR in Angus cattle from an endemic area of São Paulo state, Brazil.

The levels of infection by Babesia bovis and Babesia bigemina were estimated by absolute quantification through the quantitative PCR technique (qPCR). Fifty-one contemporaneous Angus cattle were evaluated on two occasions. The number of standard female Rhipicephalus microplus ticks present on the left side of the body was counted and blood samples were drawn from the tail vein into tubes containing the anticoagulant EDTA. The blood samples were submitted to DNA extraction and used to quantify the number of copies (NC) of DNA from B. bovis and B. bigemina by qPCR. The data on tick count and number of DNA copies were transformed for normalization and analyzed by a mixed model method. A multivariate model with repeated measures of the same animal, including the effects of collection, parasite species and their interaction, was used. The repeatability values were obtained from the matrix of (co)variances and were expressed for each species. The correlations between the counts of different species on the same animal, in the same collection or different collections, were also estimated. The results showed the qPCR could distinguish the two between infection by the two Babesia species. Infection levels by B. bovis and B. bigemina were detected in 100% and 98% of the animals, respectively. Significant differences were found (P<0.05) between the NC of the two Babesia species, B. bovis 1.49±0.07 vs. B. bigemina 0.82±0.06. Low repeatabilities were found for the counts of R. microplus and NC of B. bovis and B. bigemina: 0.05, 0.10 and 0.02, respectively. The correlations between R. microplus count and NC of B. bovis and B. bigemina were both very near zero. However, an association was observed between the NC of the two species, with a correlation coefficient of 0.30 for measures from the same collection. The absence of associations between the quantity of DNA from B. bovis and B. bigemina and the tick counts suggests that the variation of parasitemia by the hemoparasites did not depend on the tick infestation levels at the moment of each collection. The repeatability values estimated indicate that under the study conditions, the variations in the tick infestation levels and of parasitemia by B. bovis and B. bigemina depend more on factors related to each collection than on intrinsic factors of the animal.

Amblyomma americanum tick saliva serine protease inhibitor 6 is a cross-class inhibitor of serine proteases and papain-like cysteine proteases that delays plasma clotting and inhibits platelet aggregation.

We previously demonstrated that Amblyomma americanum tick serine protease inhibitor 6 (AamS6) was secreted into the host during tick feeding and that both its mRNA and protein were ubiquitously and highly expressed during the first 3 days of tick feeding. This study demonstrates that AamS6 is a cross-class inhibitor of both serine- and papain-like cysteine proteases that has apparent antihaemostatic functions. Consistent with the typical inhibitory serpin characteristics, enzyme kinetics analyses revealed that Pichia pastoris-expressed recombinant (r) AamS6 reduced initial velocities of substrate hydrolysis (V0) and/or maximum enzyme velocity (V(max)) of trypsin, chymotrypsin, elastase, chymase, and papain in a dose-response manner. We speculate that rAamS6 inhibited plasmin in a temporary fashion in that while rAamS6 reduced V0 of plasmin by up to ∼53%, it had no effect on V(max). Our data also suggest that rAmS6 has minimal or no apparent effect on V0 or V(max) of thrombin, factor Xa, and kallikrein. We speculate that AamS6 is apparently involved in facilitating blood meal feeding in that various amounts of rAamS6 reduced platelet aggregation by up to ∼47% and delayed plasma clotting time in the recalcification time assay by up to ∼210 s. AamS6 is most likely not involved with the tick's evasion of the host's complement defense mechanism, in that rAamS6 did not interfere with the complement activation pathway. Findings in this study are discussed in the context of expanding our understanding of tick proteins that control bloodmeal feeding and hence tick-borne disease transmission by ticks.

Resistance of cattle of various genetic groups to the tick Rhipicephalus microplus and the relationship with coat traits.

This study evaluated the resistance of cattle of different genetic groups to the tick Rhipicephalus microplus and the relationship with traits of the animals' hair and coat. Cows of the Senepol×Nelore (SN), Angus×Nelore (AN) and Nelore (NX) genetic groups were submitted to four consecutive artificial infestations, at 14-day intervals, each one with approximately 20,000 tick larvae placed on the animals' lumbar region. From the 19th to 23rd day of each infestation five counts of the number of ticks were performed on each animal's left body side. The tick count data (TTC) were transformed into log(10) (n+1), and also into percentage of return (PR), where n is the total number of ticks counted at each infestation. Hair samples were collected 24h after the last infestation with flat-nosed pliers. Measures of the average hair length (HL), coat thickness (CT), number of hairs per cm(2) (NHCM2) and weight of the samples (SW) were obtained. Pearson's correlation coefficients were calculated within genetic group to measure association between PR and the hair and coat data. There was a significant difference among genetic groups for the number of ticks, with the AN group having higher counts than the SN and NX groups. For the hair and coat traits, the NX and SN groups had lower values of HL and SW than did the AN group. The SN genetic group had lower NHCM2 counts than the NX and AN groups. There were positive correlations between TTC and CT (P<0.05) and SW (P<0.05) in the SN group. No significant correlation was found for the AN genetic group (P>0.05).

mRNA profile of Nellore calves after primary infection with Haemonchus placei.

Haemonchus parasites are responsible for many losses in animal production. However, few studies are available, especially of zebu cattle. In this study, we investigated mRNA differences of immune response genes in naïve Nellore calves infected with Haemonchus placei, relating these differences to patterns of cellular infiltrate. Calves were infected with 15,000 H. placei L3 larvae and after 7 days lymph node and abomasum tissues were collected. IL-2, IL-4, IL-8, IL-12, IL-13, IFN-γ, MCP-1, lysozyme, pepsinogen and TNF-α genes were evaluated by qPCR. Mast cells, eosinophils and globular leukocytes were counted by abomasum histology. In the infected group, IL-4, IL-13 and TNF-α were up-regulated in the abomasal lymph node. In the abomasum, IL-13 increased and TNF-α was down-regulated (p<0.05). No differences were detected for mast cells and eosinophil counts in abomasal tissue (p>0.05). We conclude that for this infection time, there was Th2 polarization but that cellular infiltrate in abomasal tissue takes longer to develop.

Prospecting candidate SNPs for backfat in Canchim beef cattle.

Canchim is a composite cattle breed developed in Brazil for beef production. One of the breeding objectives is to increase fat deposition. QTLs for fat thickness and/or marbling have been reported on BTA4 and BTA14. The IGFBP3 and DDEF1 genes, mapped to BTA4 and BTA14, respectively, affect adipogenesis. We looked for SNPs in the IGFBP3 and DDEF1 genes that could be associated with backfat thickness in Canchim beef cattle. For SNP identification, sires with the highest accuracy were ranked according to expected breeding value for fat thickness; the 12 extremes (six sires with the highest and six with the lowest expected breeding value for the trait) were chosen. Six regions of the IGFBP3 and 14 regions of the DDEF1 were sequenced using the Sanger method. Nine SNPs were identified in IGFBP3 and 76 in the DDEF1. After an initial analysis, two SNPs were selected to be genotyped for the whole population; these were DDEF1g.279401A>G and IGFBP3c.4394T>C(Trp>Arg). We found a significant effect (P ≤ 0.05) of allele substitution on backfat thickness; however, the IGFBP3 SNP did not significantly affect this trait.

On the search for markers of tick resistance in bovines.

Genetic differences in susceptibility to ticks (Rhipicephalus (Boophilus) microplus) are considerable in bovines. Here, mapping, association and gene expression approaches were employed to further advance our understanding of the molecular basis of tick resistance. A B. taurus x B. indicus F2 population was developed by Embrapa and 382 individuals were measured for parasitic load. Scanning of all chromosomes is in progress. Quantitative trait loci (QTL) for tick load were mapped to chromosomes 4, 5, 7, 10, 14, 18 and 23 out of the 20 chromosomes scanned and were dependent on the season in which the phenotype was scored. In the candidate gene approach, females from the genetic groups Nelore (NE--184), Canchim x Nelore (CN--153), Aberdeen Angus x Nelore (AN--123) and Simmental x Nelore (SN--120) were evaluated under natural infestation. Microsatellite markers close to the genes for interleukin 2 (IL2), interleukin 4 (IL4) and interferon gamma (IFNG) were analysed. Tick counts were associated with the marker for interleukin 4 (P < 0.05) in three genetic groups. Differences in cytokine mRNA levels of naive versus infested Nelore calves as well as between resistant versus susceptible cows from NE, CN and AN genetic groups were also investigated. Comparison of cytokines from infested and naïve animals showed downregulation of IL2. When resistant cows were compared to susceptible animals, IL8 was downregulated. These results reinforce the multiloci nature of tick resistance and the need to consider QTL and environment interactions.