Skin transcriptomic analysis reveals candidate genes and pathways associated with thermotolerance in hair sheep

The skin plays an important role in thermoregulation. Identification of genes on the skin that contribute to increased heat tolerance can be used to select animals with the best performance in warm environments. Our objective was to identify candidate genes associated with the heat stress response in the skin of Santa Ines sheep. A group of 80 sheep assessed for thermotolerance was kept in a climatic chamber for 8 days at a stress level temperature of 36 °C (10 am to 04 pm) and a maintenance temperature of 28 °C (04 pm to 10 am). Two divergent groups, with seven animals each, were formed after ranking them by thermotolerance using rectal temperature. From skin biopsy samples, total RNA was extracted, quantified, and used for RNA-seq analysis. 15,989 genes were expressed in sheep skin samples, of which 4 genes were differentially expressed (DE; FDR < 0.05) and 11 DE (FDR 0.05-0.177) between the two divergent groups. These genes are involved in cellular protection against stress (HSPA1A and HSPA6), ribosome assembly (28S, 18S, and 5S ribosomal RNA), and immune response (IGHG4, GNLY, CXCL1, CAPN14, and SAA-4). The candidate genes and main pathways related to heat tolerance in Santa Ines sheep require further investigation to understand their response to heat stress in different climatic conditions and under solar radiation. It is essential to verify whether these genes and pathways are present in different breeds and to understand the relationship between heat stress and other genes identified in this study.

Exploring candidate genes for heat tolerance in ovine through liver gene expression

Thermotolerance has become an essential factor in the prevention of the adverse effects of heat stress, but it varies among animals. Identifying genes related to heat adaptability traits is important for improving thermotolerance and for selecting more productive animals in hot environments. The primary objective of this research was to find candidate genes in the liver that play a crucial role in the heat stress response of Santa Ines sheep, which exhibit varying levels of heat tolerance. To achieve this goal, 80 sheep were selected based on their thermotolerance and placed in a climate chamber for 10 days, during which the average temperature was maintained at 36 °C from 10 a.m. to 4 p.m. and 28 °C from 4 p.m. to 10 a.m. A subset of 14 extreme animals, with seven thermotolerant and seven non-thermotolerant animals based on heat loss (rectal temperature), were selected for liver sampling. RNA sequencing and differential gene expression analysis were performed. Thermotolerant sheep showed higher expression of genes GPx3, RGS6, GPAT3, VLDLR, LOC101108817, and EVC. These genes were mainly related to the Hedgehog signaling pathway, glutathione metabolism, glycerolipid metabolism, and thyroid hormone synthesis. These enhanced pathways in thermotolerant animals could potentially mitigate the negative effects of heat stress, conferring greater heat resistance.

Dynamic clade transitions and the influence of vaccine rollout on the spatiotemporal circulation of SARS-CoV-2 variants in São Paulo, Brazil

Since 2021, the emergence of variants of concern (VOC) has led Brazil to experience record numbers of in COVID-19 cases and deaths. The expanded spread of the SARS-CoV-2 combined with a low vaccination rate has contributed to the emergence of new mutations that may enhance viral fitness, leading to the persistence of the disease. Due to limitations in the real-time genomic monitoring of new variants in some Brazilian states, we aimed to investigate whether genomic surveillance, coupled with epidemiological data and SARS-CoV-2 variants spatiotemporal spread in a smaller region, can reflect the pandemic progression at a national level. Our findings revealed three SARS-CoV-2 variant replacements from 2021 to early 2022, corresponding to the introduction and increase in the frequency of Gamma, Delta, and Omicron variants, as indicated by peaks of the Effective Reproductive Number (Reff). These distinct clade replacements triggered two waves of COVID-19 cases, influenced by the increasing vaccine uptake over time. Our results indicated that the effectiveness of vaccination in preventing new cases during the Delta and Omicron circulations was six and eleven times higher, respectively, than during the period when Gamma was predominant, and it was highly efficient in reducing the number of deaths. Furthermore, we demonstrated that genomic monitoring at a local level can reflect the national trends in the spread and evolution of SARS-CoV-2.

Metabolomics Changes in Meat and Subcutaneous Fat of Male Cattle Submitted to Fetal Programming

This study investigated changes in meat and subcutaneous fat metabolomes and possible metabolic pathways related to prenatal nutrition in beef cattle. For this purpose, 18 Nellore bulls were used for meat sampling and 15 for fat sampling. The nutritional treatments during the gestation were: NP-not programmed or control, without protein-energy supplementation; PP-partially programmed, with protein-energy supplementation (0.3% of body weight (BW)) only in the final third of pregnancy; and FP-full programming, with protein-energy supplementation (0.3% of BW) during the entire pregnancy. The meat and fat samples were collected individually 24 h after slaughter, and the metabolites were extracted using a combination of chemical reagents and mechanical processes and subsequently quantified using liquid chromatography or flow injection coupled to mass spectrometry. The data obtained were submitted to principal component analysis (PCA), analysis of variance (ANOVA), and functional enrichment analysis, with a significance level of 5%. The PCA showed an overlap between the treatments for both meat and fat. In meat, 25 metabolites were statistically different between treatments (p ≤ 0.05), belonging to four classes (glycerophospholipids, amino acids, sphingolipids, and biogenic amine). In fat, 10 significant metabolites (p ≤ 0.05) were obtained in two classes (phosphatidylcholine and lysophosphatidylcholine). The functional enrichment analysis showed alterations in the aminoacyl-tRNA pathway in meat (p = 0.030); however, there was no pathway enriched for fat. Fetal programming influenced the meat and fat metabolomes and the aminoacyl-tRNA metabolic pathway, which is an important candidate for the biological process linked to meat quality and related to fetal programming in beef cattle.

Genome-Wide Association (GWAS) Applied to Carcass and Meat Traits of Nellore Cattle

The meat market has enormous importance for the world economy, and the quality of the product offered to the consumer is fundamental for the success of the sector. In this study, we analyzed a database which contained information on 2470 animals from a commercial farm in the state of São Paulo, Brazil. Of this total, 2181 animals were genotyped, using 777,962 single-nucleotide polymorphisms (SNPs). After quality control analysis, 468,321 SNPs provided information on the number of genotyped animals. Genome-wide association analyses (GWAS) were performed for the characteristics of the rib eye area (REA), subcutaneous fat thickness (SFT), shear force at 7 days' ageing (SF7), and intramuscular fat (IMF), with the aid of the single-step genomic best linear unbiased prediction (ssGBLUP) method, with the purpose of identifying possible genomic windows (~1 Mb) responsible for explaining at least 0.5% of the genetic variance of the traits under analysis (≥0.5%). These genomic regions were used in a gene search and enrichment analyses using MeSH terms. The distributed heritability coefficients were 0.14, 0.20, 0.18, and 0.21 for REA, SFT, SF7, and IMF, respectively. The GWAS results indicated significant genomic windows for the traits of interest in a total of 17 chromosomes. Enrichment analyses showed the following significant terms (FDR ≤ 0.05) associated with the characteristics under study: for the REA, heat stress disorders and life cycle stages; for SFT, insulin and nonesterified fatty acids; for SF7, apoptosis and heat shock proteins (HSP27); and for IMF, metalloproteinase 2. In addition, KEGG (Kyoto encyclopedia of genes and genomes) enrichment analysis allowed us to highlight important metabolic pathways related to the studied phenotypes, such as the growth hormone synthesis, insulin-signaling, fatty acid metabolism, and ABC transporter pathways. The results obtained provide a better understanding of the molecular processes involved in the expression of the studied characteristics and may contribute to the design of selection strategies and future studies aimed at improving the productivity of Nellore cattle.

Retrospective Insights of the COVID-19 Epidemic in the Major Latin American City, São Paulo, Southeastern Brazil

São Paulo is the financial center of Brazil, with a population of over 12 million, that receives travelers from all over the world for business and tourism. It was the first city in Brazil to report a case of COVID-19 that rapidly spread across the city despite the implementation of the restriction measures. Despite many reports, much is still unknown regarding the genomic diversity and transmission dynamics of this virus in the city of São Paulo. Thus, in this study, we provide a retrospective overview of the COVID-19 epidemic in São Paulo City, Southeastern, Brazil, by generating a total of 9995 near-complete genome sequences from all the city's different macro-regions (North, West, Central, East, South, and Southeast). Our analysis revealed that multiple independent introduction events of different variants (mainly Gamma, Delta, and Omicron) occurred throughout time. Additionally, our estimates of viral movement within the different macro-regions further suggested that the East and the Southeast regions were the largest contributors to the Gamma and Delta viral exchanges to other regions. Meanwhile, the North region had a higher contribution to the dispersion of the Omicron variant. Together, our results reinforce the importance of increasing SARS-CoV-2 genomic monitoring within the city and the country to track the real-time evolution of the virus and to detect earlier any eventual emergency of new variants of concern that could undermine the fight against COVID-19 in Brazil and worldwide.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (−0.47), ribeye area (REA) and protein 4 (prot4) (−0.33), REA and protein 2 (prot2) (−0.3), carcass and prot4 (−0.31), carcass and prot2 (−0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (−0.25). Positive correlations were observed among the four protein factors (0.45–0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Correction: Lesbon et al. Nucleocapsid (N) Gene Mutations of SARS-CoV-2 Can Affect Real-Time RT-PCR Diagnostic and Impact False-Negative Results. Viruses 2021, 13, 2474

The authors hereby request the inclusion of two authors (Olivia Teixeira and Maria Cristina Nonato) in the recently published article in Viruses entitled “Nucleocapsid (N) gene mutations of SARS-CoV-2 can affect real-time RT-PCR diagnostic and impact false-negative results” [...]

Genomic epidemiology of the SARS-CoV-2 epidemic in Brazil

The high numbers of COVID-19 cases and deaths in Brazil have made Latin America an epicentre of the pandemic. SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, but important gaps remain in our understanding of virus transmission dynamics at a national scale. We use 17,135 near-complete genomes sampled from 27 Brazilian states and bordering country Paraguay. From March to November 2020, we detected co-circulation of multiple viral lineages that were linked to multiple importations (predominantly from Europe). After November 2020, we detected large, local transmission clusters within the country. In the absence of effective restriction measures, the epidemic progressed, and in January 2021 there was emergence and onward spread, both within and abroad, of variants of concern and variants under monitoring, including Gamma (P.1) and Zeta (P.2). We also characterized a genomic overview of the epidemic in Paraguay and detected evidence of importation of SARS-CoV-2 ancestor lineages and variants of concern from Brazil. Our findings show that genomic surveillance in Brazil enabled assessment of the real-time spread of emerging SARS-CoV-2 variants.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR &lt; 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

The Divergent Pattern of SARS-CoV-2 Variant Predominance and Transmission Dynamics in the Brazilian Island of Ilhabela

Our effort in SARS-CoV-2 genomic surveillance in Brazil has detected the Alpha Variant of Concern with a predominance higher than 75% in the population of Ilhabela island (São Paulo State) at a time when the Gamma VOC was already predominating the mainland raised concerns for closer surveillance on this island. Therefore, we intensified the surveillance for 24 weeks by generating data from 34% of local positive cases. Our data show that the patterns of VOC predominance dynamics and infection rates were in general distinct from the mainland. We report here the first known case of Alpha predominance in a Brazilian population, a delay greater than 3 months for the Gamma to dominate the previous variants compared to the mainland, and a faster dispersion rate of Gamma and Delta VOCs compared to the mainland. Phylogenetic analysis revealed the SARS-CoV-2 transmission dynamics in Ilhabela were characterized by multiple independent introduction events of Gamma and Delta, with a few events of Alpha introduction, two of them followed by community transmission. This study evidenced the peculiar behavior of SARS-CoV-2 variants in an isolated population and brought to light the importance of specific programs for SARS-CoV-2 genomic surveillance in isolated populations.

Genomic epidemiology reveals the impact of national and international restrictions measures on the SARS-CoV-2 epidemic in Brazil

Brazil has experienced some of the highest numbers of COVID-19 cases and deaths globally and from May 2021 made Latin America a pandemic epicenter. Although SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, important gaps remain in our understanding of virus transmission dynamics at the national scale. Here, we describe the genomic epidemiology of SARS-CoV-2 using near-full genomes sampled from 27 Brazilian states and a bordering country - Paraguay. We show that the early stage of the pandemic in Brazil was characterised by the co-circulation of multiple viral lineages, linked to multiple importations predominantly from Europe, and subsequently characterized by large local transmission clusters. As the epidemic progressed under an absence of effective restriction measures, there was a local emergence and onward international spread of Variants of Concern (VOC) and Variants Under Monitoring (VUM), including Gamma (P.1) and Zeta (P.2). In addition, we provide a preliminary genomic overview of the epidemic in Paraguay, showing evidence of importation from Brazil. These data reinforce the usefulness and need for the implementation of widespread genomic surveillance in South America as a toolkit for pandemic monitoring that provides a means to follow the real-time spread of emerging SARS-CoV-2 variants with possible implications for public health and immunization strategies.

SARS-COV-2 GENOMIC MONITORING IN THE STATE OF SÃO PAULO UNVEILS TWO EMERGING AY.43 SUBLINEAGES

Delta VOC is highly diverse with more than 120 sublineages already described as of November 30, 2021. In this study, through active monitoring of circulating severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) variants in the state of São Paulo, southeast Brazil, we identified two emerging sublineages from the ancestral AY.43 strain which were classified as AY.43.1 and AY.43.2. These sublineages were defined by the following characteristic nonsynonymous mutations ORF1ab:A4133V and ORF3a:T14I for the AY.43.1 and ORF1ab:G1155C for the AY.43.2 and our analysis reveals that they might have a likely‐Brazilian origin. Much is still unknown regarding their dissemination in the state of São Paulo and Brazil as well as their potential impact on the ongoing vaccination process. However, the results obtained in this study reinforce the importance of genomic surveillance activity for timely identification of emerging SARS‐CoV‐2 variants which can impact the ongoing SARS‐CoV‐2 vaccination and public health policies.

Genomic monitoring of the SARS-CoV-2 B1.1.7 (WHO VOC Alpha) in the Sao Paulo state, Brazil

The SARS-CoV-2 alpha VOC (also known as lineage B.1.1.7) initially described in the autumn, 2020 in UK, rapidly became the dominant lineage across much of Europe. Despite multiple studies reporting molecular evidence suggestive of its circulation in Brazil, much is still unknown about its genomic diversity in the state of São Paulo, the main Brazilian economic and transportation hub. To get more insight regarding its transmission dynamics into the State we performed phylogenetic analysis on all alpha VOC strains obtained between February and August 2021 from the Sao Paulo state Network for Pandemic Alert of Emerging SARS-CoV-2 variants. The performed phylogenetic analysis showed that most of the alpha VOC genomes were interspersed with viral strains sampled from different Brazilian states and other countries suggesting that multiple independent Alpha VOC introductions from Brazil and overseas have occurred in the São Paulo State over time. Nevertheless, large monophyletic clusters were also observed especially from the Central-West part of the São Paulo State (the city of Bauru) and the metropolitan region of the São Paulo city. Our results highlight the Alpha VOC molecular epidemiology in the São Paulo state and reinforce the need for continued genomic surveillance strategies for the real-time monitoring of potential emerging SARS-CoV-2 variants during the ever-growing vaccination process.

Allele and genotype frequency for milk beta-casein in dairy cattle in the northern region of Tocantins State, Brazil

At present, there is a concern about the quality of milk and diseases related to its consumption, as it can generate discomfort and allergic reactions in some individuals due to its protein components. Thus, the present study was developed to identify the allele and genotype frequencies of genes for β casein, A1 and A2, in dairy herds in the region of Araguaína-TO, Brazil. Genetic material from 421 animals (crossbred dairy cattle in lactation) was used. All animals were numbered for identification, and DNA samples were extracted from hair bulbs. Samples for two markers from the polymorphic regions were characterized and confirmed by real time PCR using the ABI Prism® 7500 Sequence Detection System (Applied Biosystems). Allele and genotype frequencies were determined using the TaqMan™ detection system, where the primer and probe release different fluorescence signals for each allele of the polymorphism. The sampled herd showed frequencies of 28.27% for the A1 allele and 71.73% for the A2 allele. Genotype frequencies were 52.96% (223/421) for A2A2; 37.53% (158/421) for the A1A2 genotype; and 9.50% (40/421) for the A1A1 genotype. The frequency of the A1 allele for β-casein in dairy herds from the northern region of Tocantins was low and is per the results of previous studies. Although the A2A2 genotype of β-casein had a high relative frequency, the A1A2 genotype is still rather frequent, warranting greater selection pressure.

Nucleocapsid (N) Gene Mutations of SARS-CoV-2 Can Affect Real-Time RT-PCR Diagnostic and Impact False-Negative Results

The current COVID-19 pandemic demands massive testing by Real-time RT-PCR (Reverse Transcription Polymerase Chain Reaction), which is considered the gold standard diagnostic test for the detection of the SARS-CoV-2 virus. However, the virus continues to evolve with mutations that lead to phenotypic alterations as higher transmissibility, pathogenicity or vaccine evasion. Another big issue are mutations in the annealing sites of primers and probes of RT-PCR diagnostic kits leading to false-negative results. Therefore, here we identify mutations in the N (Nucleocapsid) gene that affects the use of the GeneFinder COVID-19 Plus RealAmp Kit. We sequenced SARS-CoV-2 genomes from 17 positive samples with no N gene detection but with RDRP (RNA-dependent RNA polymerase) and E (Envelope) genes detection, and observed a set of three different mutations affecting the N detection: a deletion of 18 nucleotides (Del28877-28894), a substitution of GGG to AAC (28881-28883) and a frameshift mutation caused by deletion (Del28877-28878). The last one cause a deletion of six AAs (amino acids) located in the central intrinsic disorder region at protein level. We also found this mutation in 99 of the 14,346 sequenced samples by the Sao Paulo state Network for Pandemic Alert of Emerging SARS-CoV-2 variants, demonstrating the circulation of the mutation in Sao Paulo, Brazil. Continuous monitoring and characterization of mutations affecting the annealing sites of primers and probes by genomic surveillance programs are necessary to maintain the effectiveness of the diagnosis of COVID-19.

Introduction of SARS-CoV-2 C.37 (WHO VOI lambda) in the Sao Paulo State, Southeast Brazil

The Lambda variants of interest (VOI) (C37/GR/452Q.V1/21G) was initially reported in Lima, Peru but has gained rapid dissemination through other Latin American countries. Nevertheless, the dissemination and molecular epidemiology of the Lambda VOI in Brazil is unknown apart from a single case report. In this respect, we characterized the circulation of the SARS‐CoV‐2 Lambda VOI (C37/GR/452Q.V1/21G) in Sao Paulo State, Brazil. From March to June 2021, we identified seven Lambda isolates in a set of approximately 8000 newly sequenced genomes of the Network for Pandemic Alert of Emerging SARS‐CoV‐2 variants from Sao Paulo State. Interestingly, in three of the positive patients, the Lambda VOI infection was probably related to a contact transmission. These individuals were fully vaccinated to COVID‐19 and presented mild symptoms. The remaining positive for Lambda VOI individuals showed different levels of COVID‐19 symptoms and one of them needed hospitalization (score 5, WHO). In our study, we present a low level of Lambda VOI circulation in the Sao Paulo State. This reinforces the essential role of molecular surveillance for the effective SARS‐CoV‐2 pandemic response, especially in regard to circulating variants.

Our study shows the low circulation of the Lambda variants of interest (VOI) in the São Paulo State.

Two distinct Lambda VOI introductions may have occurred in the Sao Paulo State.

Epidemiological data revealed Lambda VOI transmission in three cases of fully vaccinated individuals.

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.

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

Background

Feed efficiency and growth rate have been targets for selection to improve chicken production. The incorporation of genomic tools may help to accelerate selection. We genotyped 529 individuals using a high-density SNP chip (600 K, Affymetrix®) to estimate genomic heritability of performance traits and to identify genomic regions and their positional candidate genes associated with performance traits in a Brazilian F₂ Chicken Resource population. Regions exhibiting selection signatures and a SNP dataset from resequencing were integrated with the genomic regions identified using the chip to refine the list of positional candidate genes and identify potential causative mutations.

Results

Feed intake (FI), feed conversion ratio (FC), feed efficiency (FE) and weight gain (WG) exhibited low genomic heritability values (i.e. from 0.0002 to 0.13), while body weight at hatch (BW1), 35 days-of-age (BW35), and 41 days-of-age (BW41) exhibited high genomic heritability values (i.e. from 0.60 to 0.73) in this F₂ population. Twenty unique 1-Mb genomic windows were associated with BW1, BW35 or BW41, located on GGA1–4, 6–7, 10, 14, 24, 27 and 28. Thirty-eight positional candidate genes were identified within these windows, and three of them overlapped with selection signature regions. Thirteen predicted deleterious and three high impact sequence SNPs in these QTL regions were annotated in 11 positional candidate genes related to osteogenesis, skeletal muscle development, growth, energy metabolism and lipid metabolism, which may be associated with body weight in chickens.

Conclusions

The use of a high-density SNP array to identify QTL which were integrated with whole genome sequence signatures of selection allowed the identification of candidate genes and candidate causal variants. One novel QTL was detected providing additional information to understand the genetic architecture of body weight traits. We identified QTL for body weight traits, which were also associated with fatness in the same population. Our findings form a basis for further functional studies to elucidate the role of specific genes in regulating body weight and fat deposition in chickens, generating useful information for poultry breeding programs.

Genetic variants in glucocorticoid and mineralocorticoid receptors are associated with concentrations of plasma cortisol, muscle glycogen content, and meat quality traits in male Nellore cattle

The glucocorticoid receptor (GR) and mineralocorticoid receptor (MR) are key components in the regulation of the hypothalamic-pituitary-adrenal neuroendocrine axis and coordinate the physiological response to stress agents to reestablish homeostasis. Genetic variations of GR (NR3C1) and MR (NR3C2) genes could explain the alterations in animals to adapt to challenges, and therefore, their influence on production traits. The present study aimed to identify single-nucleotide polymorphisms (SNPs) in the bovine NR3C1 and NR3C2 genes and explore their associations to relevant traits of beef cattle production. Genotypes and phenotypes were collected from 241 male Nellore cattle (119 noncastrated and 122 castrated surgically) with an average of 24 ± 1.2 mo of age and live weight of 508 ± 39 kg. The traits evaluated were concentrations of plasma adrenocorticotropic hormone (ACTH) and cortisol, muscle glycogen and lactate content, and pH, color, cooking loss, and shear force of longissimus thoracis measured on the 1st, 7th, and 14th days postmortem. Five SNPs were identified, 2 in the NR3C1 gene and 3 in the NR3C2 gene. There was an associative relationship between the SNP NR3C1_1 g.3293A>G and postmortem plasma concentration of cortisol (P = 0.0008). The SNPs NR3C2_1 g.115T>C and NR3C2_2 g.570T>C were associated with muscle glycogen content (P = 0.0306 and P = 0.0158), postmortem plasma concentration of ACTH (P = 0.0118 and P = 0.0095), and cooking loss of the steak aged 1 d (P = 0.0398 and P = 0.0423). Haplotype analysis showed associations of GR haplotypes with postmortem plasma concentrations of cortisol and MR haplotypes with meat color, cooking losses, muscle glycogen content, and plasma concentrations of ACTH. The associations observed in the present study show that SNPs in GR and MR genes are related with changes of hypothalamic-pituitary-adrenal axis activity and metabolic profile in cattle, leading to individual variation in meat quality traits.

Heat shock and structural proteins associated with meat tenderness in Nellore beef cattle, a Bos indicus breed

Nellore beef cattle, a Bos indicus (Zebu) breed, is well adapted to tropical conditions and has allowed Brazil to become one of the largest producers of red meat. Nevertheless, B. indicus breeds are reported to have less tender meat than Bos taurus. This study was designed to identify genes associated with meat tenderness and thus provides important information for breeding programs. A group of 138 animals was evaluated for longissimus thoracis muscle shear force (SF). Animals with the highest and lowest SF values (six animals each) were then selected for protein abundance studies. Samples were subjected to two-dimensional gel electrophoresis (2-DE) followed by peptide sequencing through mass spectrometry (MS) to identify differentially expressed proteins associated with SF values. Seventeen differentially expressed spots were observed (p<0.05) between the two groups. The 13 proteins identified included structural proteins (alpha actin-1, MLC1, MLC3, MLC2F and tropomyosin), related to cell organization (HSPB1 and HSP70), metabolism (beta-LG, ACBD6 and Complex III subunit I) and some uncharacterized proteins. Results confirm the existence of differentially expressed proteins associated with SF, which can lead to a better understanding of mechanisms involved in meat tenderness.

Relationship between expression of voltage-dependent anion channel (VDAC) isoforms and type of hexokinase binding sites on brain mitochondria

Voltage-dependent anion channels (VDAC) are pore-forming proteins found in the outer mitochondrial membrane of eukaryotes. VDACs are known to play an essential role in cellular metabolism and in early stages of apoptosis. In mammals, three VDAC isoforms have been identified. A proteomic approach was exploited to study the expression of VDAC isoforms in rat, bovine, and chicken brain mitochondria. Given the importance of mitochondrially bound hexokinase in regulation of aerobic glycolysis in brain, we studied the possibility that differences in the relative expression of VDAC isoforms may be a factor in determining the species-dependent ratio of type A/type B hexokinase binding sites on brain mitochondria. The spots were characterized, and the signal intensities among spots were compared. VDAC1 was the most abundantly expressed of the three isoforms. Moreover the expression of VDAC1 plus VDAC2 was significantly higher in bovine than in rat brain. Chicken brain mitochondria showed the highest VDAC1 expression and the lowest of VDAC2. Bovine brain mitochondria had the highest VDAC2 levels. We concluded that the nature of hexokinase binding site is not determined by the expression of a single VDAC isoform.

Effect of Cysticercus cellulosae fractions on the respiratory burst of pig neutrophils

Neutrophils, eosinophils and macrophages are cells that interact with invading parasites and naive hosts have been shown to have anti-parasitic activity. The initial reaction of these leukocytes is the generation of reactive oxygen species (ROS) to play in parasite expulsion. The present work was carried out to study the effect of total extract, scolex and membrane fractions from Cysticercus cellulosae on respiratory burst by pig neutrophils. Hydrogen peroxide (H2O2) production by neutrophils incubated with metacestode fractions from C. cellulosae showed an increase of: 190% (total extract), 120% (scolex) and 44% (membrane). High antioxidant catalatic activity (33%, 28%, 28% by total extract, scolex and membrane, respectively) was observed in neutrophils incubated with metacestode fractions, which could be an attempt at self-protection. Scolex and membrane fractions increased the phagocytic capacity of neutrophils (44% and 28%, respectively). On the other hand, total cysticerci did not alter the phagocytosis, possibly due to modifications in membrane function, caused by high ROS production from neutrophils in the presence of total cysticerci. Total fraction from C. cellulosae is toxic for neutrophils as shown by the decrease in phagocytic capacity, probably caused by high levels of ROS formation. The difference in toxicity of total extract, scolex and membrane fractions on neutrophils can be explained by the presence of an antigenic effect of the vesicular fluid in the total extract of C. cellulosae.