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.

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

BACKGROUND

Understanding the molecular underpinnings of phenotypic variations is critical for enhancing poultry breeding programs. The Brazilian broiler (TT) and laying hen (CC) lines exhibit striking differences in body weight, growth potential, and muscle mass. Our work aimed to compare the global transcriptome of wing and pectoral tissues during the early development (days 2.5 to 3.5) of these chicken lines, unveiling disparities in gene expression and regulation.

RESULTS

Different and bona-fide transcriptomic profiles were identified for the compared lines. A similar number of up- and downregulated differentially expressed genes (DEGs) were identified, considering the broiler line as a reference. Upregulated DEGs displayed an enrichment of protease-encoding genes, whereas downregulated DEGs exhibited a prevalence of receptors and ligands. Gene Ontology analysis revealed that upregulated DEGs were mainly associated with hormone response, mitotic cell cycle, and different metabolic and biosynthetic processes. In contrast, downregulated DEGs were primarily linked to communication, signal transduction, cell differentiation, and nervous system development. Regulatory networks were constructed for the mitotic cell cycle and cell differentiation biological processes, as their contrasting roles may impact the development of distinct postnatal traits. Within the mitotic cell cycle network, key upregulated DEGs included CCND1 and HSP90, with central regulators being NF-κB subunits (RELA and REL) and NFATC2. The cell differentiation network comprises numerous DEGs encoding transcription factors (e.g., HOX genes), receptors, ligands, and histones, while the main regulatory hubs are CREB, AR and epigenetic modifiers. Clustering analyses highlighted PIK3CD as a central player within the differentiation network.

CONCLUSIONS

Our study revealed distinct developmental transcriptomes between Brazilian broiler and layer lines. The gene expression profile of broiler embryos seems to favour increased cell proliferation and delayed differentiation, which may contribute to the subsequent enlargement of pectoral tissues during foetal and postnatal development. Our findings pave the way for future functional studies and improvement of targeted traits of economic interest in poultry.

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.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs

BACKGROUND

Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs.

RESULTS

The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes.

CONCLUSIONS

After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

Rumen and fecal microbiomes are related to diet and production traits in Bos indicus beef cattle

Background

Ruminants harbor a complex microbial community within their gastrointestinal tract, which plays major roles in their health and physiology. Brazil is one of the largest producers of beef in the world and more than 90% of the beef cattle herds are composed of pure and crossbred Nelore (Bos indicus). Despite its importance to the Brazilian economy and human feeding, few studies have characterized the Nelore microbiome. Therefore, using shotgun metagenomics, we investigated the impact of diet on the composition and functionality of the Nelore microbiome, and explored the associations between specific microbial taxa and their functionality with feed efficiency and methane emission.

Results

The ruminal microbiome exhibited significantly higher microbial diversity, distinctive taxonomic profile and variations in microbial functionality compared to the fecal microbiome, highlighting the distinct contributions of the microbiomes of these environments. Animals subjected to different dietary treatments exhibited significant differences in their microbiomes' archaeal diversity and in the abundance of 89 genera, as well as in the functions associated with the metabolism of components of each diet. Moreover, depending on the diet, feed-efficient animals and low methane emitters displayed higher microbial diversity in their fecal microbiome. Multiple genera were associated with an increase or decrease of the phenotypes. Upon analyzing the functions attributed to these taxa, we observed significant differences on the ruminal taxa associated with feed efficient and inefficient cattle. The ruminal taxa that characterized feed efficient cattle stood out for having significantly more functions related to carbohydrate metabolism, such as monosaccharides, di-/oligosaccharides and amino acids. The taxa associated with methane emission had functions associated with methanogenesis and the production of substrates that may influence methane production, such as hydrogen and formate.

Conclusion

Our findings highlight the significant role of diet in shaping Nelore microbiomes and how its composition and functionality may affect production traits such as feed efficiency and methane emission. These insights provide valuable support for the implementation of novel feeding and biotechnological strategies.

Ruminal and feces metabolites associated with feed efficiency, water intake and methane emission in Nelore bulls

The objectives of this study were twofold: (1) to identify potential differences in the ruminal and fecal metabolite profiles of Nelore bulls under different nutritional interventions; and (2) to identify metabolites associated with cattle sustainability related-traits. We used different nutritional interventions in the feedlot: conventional (Conv; n = 26), and by-product (ByPr, n = 26). Thirty-eight ruminal fluid and 27 fecal metabolites were significantly different (P < 0.05) between the ByPr and Conv groups. Individual dry matter intake (DMI), residual feed intake (RFI), observed water intake (OWI), predicted water intake (WI), and residual water intake (RWI) phenotypes were lower (P < 0.05) in the Conv group, while the ByPr group exhibited lower methane emission (ME) (P < 0.05). Ruminal fluid dimethylamine was significantly associated (P < 0.05) with DMI, RFI, FE (feed efficiency), OWI and WI. Aspartate was associated (P < 0.05) with DMI, RFI, FE and WI. Fecal C22:1n9 was significantly associated with OWI and RWI (P < 0.05). Fatty acid C14:0 and hypoxanthine were significantly associated with DMI and RFI (P < 0.05). The results demonstrated that different nutritional interventions alter ruminal and fecal metabolites and provided new insights into the relationship of these metabolites with feed efficiency and water intake traits in Nelore bulls.

Genetic parameters for production, quality, and colors from eggs in Brazilian lineages of chickens

Genetic parameters were estimated for egg production, egg quality, and eggshell colors in eight lineages of Brazilian laying hens. Age at first egg (AFE), total egg production up to the 45th week (PROD), egg weight (EW), albumen height (AH), yolk color (YC), the Haugh units (HU), eggshell strength (ESS), eggshell thickness (EST), yolk weight (YW), eggshell weight (ESW), and eggshell color (L*, a*, and b*) were measured in 2030 eggs obtained from 645 laying hens. Variance components were estimated from a mixed animal model, which included the fixed effects of contemporary groups, cage location, and hen line, and the additive genetic, permanent environmental, and residual as random effects. In general, heritabilities were low to moderate (h² = 0.11 to 0.48). Genetic correlations among eggshell quality traits were moderate to high (0.36 and 0.69). High genetic correlations were obtained between the eggshell color traits [r_g = -0.90 (L* and a*); r_g = -0.64 (L* and b*); and r_g = 0.65 (a* and b*)]. Results suggest that EW is strongly correlated with ESW, but the genetic correlations between EW and ESS and between EW and EST were low. Genetic correlations between L* and eggshell quality traits were low to moderate, suggesting that L* has little or no relation with external egg quality. However, genetic correlations between a* and b* values and eggshell quality traits were high. The genetic correlations between eggshell color and eggshell quality traits were low, suggesting that the eggshell color does not influence external egg quality. Genetic correlations between PROD and egg quality traits were negative and varied between -0.42 and -0.05. This antagonistic relationship emphasizes the importance of adopting breeding schemes that allow the simultaneous genetic progress of these traits by considering their genetic correlation and economic relevance, such as the selection index.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil

BACKGROUND

The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue.

RESULTS

Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified.

CONCLUSION

The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

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.

Repressive epigenetic mechanisms, such as the H3K27me3 histone modification, were predicted to affect muscle gene expression and its mineral content in Nelore cattle

Epigenetic repression has been linked to the regulation of different cell states. In this study, we focus on the influence of this repression, mainly by H3K27me3, over gene expression in muscle cells, which may affect mineral content, a phenotype that is relevant to muscle function and beef quality. Based on the inverse relationship between H3K27me3 and gene expression (i.e., epigenetic repression) and on contrasting sample groups, we computationally predicted regulatory genes that affect muscle mineral content. To this end, we applied the TRIAGE predictive method followed by a rank product analysis. This methodology can predict regulatory genes that might be affected by repressive epigenetic regulation related to mineral concentration. Annotation of orthologous genes, between human and bovine, enabled our investigation of gene expression in the Longissimus thoracis muscle of Bos indicus cattle. The animals under study had a contrasting mineral content in their muscle cells. We identified candidate regulatory genes influenced by repressive epigenetic mechanisms, linking histone modification to mineral content in beef samples. The discovered candidate genes take part in multiple biological pathways, i.e., impulse transmission, cell signalling, immunological, and developmental pathways. Some of these genes were previously associated with mineral content or regulatory mechanisms. Our findings indicate that epigenetic repression can partially explain the gene expression profiles observed in muscle samples with contrasting mineral content through the candidate regulators here identified.

A traveling SARS-CoV-2 laboratory as part of a pandemic response among vulnerable Brazilian populations

BACKGROUND

Brazil has been dramatically hit by the SARS-CoV-2 pandemic and is a world leader in COVID-19 morbidity and mortality. Additionally, the largest country of Latin America has been a continuous source of SARS-CoV-2 variants and shows extraordinary variability of the pandemic strains probably related to the country´s outstanding position as a Latin American economical and transportation hub. Not all regions of the country show sufficient infrastructure for SARS-CoV-2 diagnosis and genotyping which can negatively impact the pandemic response.

METHODS

Due to this reason and to disburden the diagnostic system of the inner São Paulo State, the Butantan Institute established the Mobile Laboratory (in Portuguese: LabMovel) for SARS-CoV-2 testing which started a trip of the most important "hotspots" of the most populous Brazilian region. The LabMovel initiated in two important cities of the State: Aparecida do Norte (an important religious center) and the Baixada Santista region which incorporates the port of Santos, the busiest in Latin America. The LabMovel was fully equipped with an automatized system for SARS-CoV-2 diagnosis and sequencing/genotyping. It also integrated the laboratory systems for patient records and results divulgation including in the Federal Brazilian Healthcare System.

RESULTS

Currently,16,678 samples were tested, among them 1,217 from Aparecida and 4,564 from Baixada Santista. We tracked the delta introductio in the tested regions with its high diversification. The established mobile SARS-CoV-2 laboratory had a major impact on the Public Health System of the included cities including timely delivery of the results to the healthcare agents and the Federal Healthcare system, evaluation of the vaccination status of the positive individuals in the background of exponential vaccination process in Brazil and scientific and technological divulgation of the fieldwork to the most vulnerable populations.

CONCLUSIONS

The SARS-CoV-2 pandemic has demonstrated worldwide the importance of science to fight against this viral agent and the LabMovel shows that it is possible to integrate researchers, clinicians, healthcare workers and patients to take rapid actions that can in fact mitigate this and other epidemiological situations.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

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.

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.

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 &lt; 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.

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.

Plant genetic diversity by DNA barcoding to investigate propolis origin

Identify the botanical origins of a certain type of propolis may be challenging and time demanding, since it involves bee's behavior observation, plant resins collection and chemical analysis. Thus, this study aimed to determine the plant genetic materials in propolis from southern Brazil using the DNA barcoding to investigate their botanical origins, as well as to compare it with the phytochemical composition determined by ultra-high-performance liquid chromatography coupled with high-resolution mass spectrometry (UHPLC-HRMS) and with the pollinic profile. As principal results, non-native Populus carolinensis Moench (Salicaceae) was almost the only DNA source in some propolis samples, which coincided with the presence of flavonoids typical from poplar exudates. Conversely, other propolis samples had DNA material coming mainly from native plant species, most of them characterized to the species level, although no specific chemical markers from those plants could be identified by UHPLC-HRMS. However, pollen from several plants identified by the DNA barcoding were extracted from some propolis samples. Despite the identification of typical diterpenes, DNA material from Araucaria angustifolia (Bertol.) Kuntze (Araucariaceae), which have been indicated as a major resin source for propolis from preservation areas in southern Brazil, was found in very small abundancies, likely because bees do not drag tissue material containing DNA when collecting resin from this native species. In conclusion, DNA barcoding analysis successfully provided information about the provenance of propolis, although, depending on the plant resin sources, this information is likely to come from pollen.

Simulation of group testing scenarios can boost COVID-19 screening power

The COVID-19 has severely affected economies and health systems around the world. Mass testing could work as a powerful alternative to restrain disease dissemination, but the shortage of reagents is a limiting factor. A solution to optimize test usage relies on ‘grouping’ or ‘pooling’ strategies, which combine a set of individuals in a single reaction. To compare different group testing configurations, we developed the poolingr package, which performs an innovative hybrid in silico/in vitro approach to search for optimal testing configurations. We used 6759 viral load values, observed in 2389 positive individuals, to simulate a wide range of scenarios. We found that larger groups (>100) framed into multi-stage setups (up to six stages) could largely boost the power to detect spreaders. Although the boost was dependent on the disease prevalence, our method could point to cheaper grouping schemes to better mitigate COVID-19 dissemination through identification and quarantine recommendation for positive individuals.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore 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.

SARS-CoV-2 epidemic in Brazil: how the displacement of variants has driven distinct epidemic waves

Brazil ranks as third in terms of total number of reported SARS-CoV-2 cases globally. The COVID-19 epidemic in Brazil was characterised by the co-circulation of multiple variants as a consequence of multiple independent introduction events occurring through time. Here, we describe the SARS-CoV-2 variants that are currently circulating and co-circulating in the country, with the aim to highlight which variants have driven the different epidemic waves. For this purpose, we retrieved metadata information of Coronavirus sequences collected in Brazil and available at the GISAID database. SARS-CoV-2 lineages have been identified along with eleven variants, labelled as VOCs (Alpha, Gamma, Beta, Delta and Omicron) VOIs (Lambda and Mu) VUMs (B.1.1.318) and FMVs (Zeta, Eta and B.1.1.519). Here we show that, in the Brazilian context, after 24 months of sustained transmission and evolution of SARS-CoV-2, local variants (among them the B.1.1.28 and B.1.1.33) were displaced by recently introduced VOCs firstly with the Gamma, followed by Delta and more recently Omicron. The rapid spread of some of those VOCs (such as Gamma and Omicron) was also mirror by a large increase in the number of cases and deaths in the country. This in turn reinforces that, due to the emergence of variants that appear to induce a substantial evasion against neutralizing antibody response, it is important to strengthen genomic effort within the country and how vaccination still remains a critical process to protect the vulnerable population, still at risk of infection and death.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs

Simple Summary

Findings from the analysis of the pig transcriptome may help to better understand the biological mechanisms that can be modulated by the diet. Thus, the aim of this study was to identify the differentially expressed genes from the skeletal muscle and liver samples of pigs fed diets with two different levels of soybean oil (1.5 or 3%). The FA profile in the tissues was modified by the diet mainly related to monounsaturated (MUFA) and polyunsaturated (PUFA). This nutrigenomics study verified the effect of different levels of soybean oil in the pig diet on the transcriptome profile of skeletal muscle and liver, where the higher level of soybean oil added to the diet led to a higher expression of genes targeting biological processes related to lipid oxidation and consequently to metabolic diseases and inflammation.

Abstract

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Land use and roles of soil bacterial community in the dissipation of atrazine

Atrazine (ATZ) is one of the most widely used herbicides in the world even though it is classified as a carcinogenic endocrine disruptor. This study focused on how land use (grazing versus cultivation in parallel soils, the latter under no-till with a seven-year history of ATZ application) and bacterial community diversity affected ATZ dissipation. Proteobacteria, Bacteroidetes, Firmicutes, Actinobacteria, Acidobacteria, Verrucomicrobia, Planctomycetes, and Gemmatimonadetes were the dominant phyla in both soils. The mineralization of ATZ was much higher in soils under cultivation up to the onset of moderate diversity depletion (dilution =10^-3), corresponding to 44-52% of the amount applied (< 5% in the grazed soil). This was attributed to the higher diversity and complexity of the soils´ bacterial communities which consist of microbial groups that were more adapted as a result of previous exposure to ATZ. In these cases, ATZ dissipation was attributed mainly to mineralization (DT₅₀ = 4-11 d). However, formation of non-extractable ATZ residues was exceptionally important in the other cases (DT₅₀ = 17-44 d). The cultivated soils also presented a higher number of bacterial genera correlated with ATZ dissipation, in which Acidothermus, Aquicela, Arenimonas, Candidatus_Koribacter, Hirschia, MND1, Nitrospira, Occallatibacter, OM27_clade, and Ralstonia are suggested as potential ATZ-degraders. Finally, ATZ dissipation was mostly associated with an abundance of microbial functions related to energy supply and N-metabolism, suggesting co-metabolism is its first biodegradation step.

Effect of dietary soybean oil inclusion on liver-related transcription factors in a pig model for metabolic diseases

Dietary fatty acids (FA) are components of the lipids, which contribute to membrane structure, energy input, and biological functions related to cellular signaling and transcriptome regulation. However, the consumers still associate dietary FA with fat deposition and increased occurrence of metabolic diseases such as obesity and atherosclerosis. Previous studies already demonstrated that some fatty acids are linked with inflammatory response, preventing metabolic diseases. To better understand the role of dietary FA on metabolic diseases, for the first time, a study to identify key transcription factors (TF) involved in lipid metabolism and inflammatory response by transcriptome analysis from liver samples of animal models was performed. The key TF were identified by functional enrichment analysis from the list of differentially expressed genes identified in liver samples between 35 pigs fed with 1.5% or 3.0% soybean oil. The functional enrichment analysis detected TF linked to lipid homeostasis and inflammatory response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key TF related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive food technology and science.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus

BACKGROUND

Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus.

RESULTS

Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity.

CONCLUSIONS

We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

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.

Replacement of the Gamma by the Delta variant in Brazil: Impact of lineage displacement on the ongoing pandemic

The coronavirus disease 2019 (COVID-19) epidemic in Brazil was driven mainly by the spread of Gamma (P.1), a locally emerged variant of concern (VOC) that was first detected in early January 2021. This variant was estimated to be responsible for more than 96 per cent of cases reported between January and June 2021, being associated with increased transmissibility and disease severity, a reduction in neutralization antibodies and effectiveness of treatments or vaccines, and diagnostic detection failure. Here we show that, following several importations predominantly from the USA, the Delta variant rapidly replaced Gamma after July 2021. However, in contrast to what was seen in other countries, the rapid spread of Delta did not lead to a large increase in the number of cases and deaths reported in Brazil. We suggest that this was likely due to the relatively successful early vaccination campaign coupled with natural immunity acquired following prior infection with Gamma. Our data reinforce reports of the increased transmissibility of the Delta variant and, considering the increasing concern due to the recently identified Omicron variant, argues for the necessity to strengthen genomic monitoring on a national level to quickly detect the emergence and spread of other VOCs that might threaten global health.

Transcriptomic characterization of the molecular mechanisms induced by RGMa during skeletal muscle nuclei accretion and hypertrophy

Background

The repulsive guidance molecule a (RGMa) is a GPI-anchor axon guidance molecule first found to play important roles during neuronal development. RGMa expression patterns and signaling pathways via Neogenin and/or as BMP coreceptors indicated that this axon guidance molecule could also be working in other processes and diseases, including during myogenesis. Previous works from our research group have consistently shown that RGMa is expressed in skeletal muscle cells and that its overexpression induces both nuclei accretion and hypertrophy in muscle cell lineages. However, the cellular components and molecular mechanisms induced by RGMa during the differentiation of skeletal muscle cells are poorly understood. In this work, the global transcription expression profile of RGMa-treated C2C12 myoblasts during the differentiation stage, obtained by RNA-seq, were reported.

Results

RGMa treatment could modulate the expression pattern of 2,195 transcripts in C2C12 skeletal muscle, with 943 upregulated and 1,252 downregulated. Among them, RGMa interfered with the expression of several RNA types, including categories related to the regulation of RNA splicing and degradation. The data also suggested that nuclei accretion induced by RGMa could be due to their capacity to induce the expression of transcripts related to ‘adherens junsctions’ and ‘extracellular-cell adhesion’, while RGMa effects on muscle hypertrophy might be due to (i) the activation of the mTOR-Akt independent axis and (ii) the regulation of the expression of transcripts related to atrophy. Finally, RGMa induced the expression of transcripts that encode skeletal muscle structural proteins, especially from sarcolemma and also those associated with striated muscle cell differentiation.

Conclusions

These results provide comprehensive knowledge of skeletal muscle transcript changes and pathways in response to RGMa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08396-w.

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.

GBS-MeDIP: A protocol for parallel identification of genetic and epigenetic variation in the same reduced fraction of genomes across individuals

The GBS-MeDIP protocol combines two previously described techniques, Genotype-by-Sequencing (GBS) and Methylated-DNA-Immunoprecipitation (MeDIP). Our method allows for parallel and cost-efficient interrogation of genetic and methylomic variants in the DNA of many reduced genomes, taking advantage of the barcoding of DNA samples performed in the GBS and the subsequent creation of DNA pools, then used as an input for the MeDIP. The GBS-MeDIP is particularly suitable to identify genetic and methylomic biomarkers when resources for whole genome interrogation are lacking.

•

Our protocol combines two well-known techniques, GBS and MeDIP.

•

The GBS-MeDIP allows for parallel interrogation of genetic/methylomic variants.

•

Allows for cost-efficient analysis of reduced genomes from multiple individuals.

•

Ideal for the identification of genetic and/or methylomic markers at reduced cost.

CNV detection and their association with growth, efficiency and carcass traits in Santa Inês sheep

Copy number variations (CNV) are an important source of genetic variation. CNV has been increasingly studied and frequently associated with diseases and productive traits in livestock animals. However, CNV-based genome-wide association studies (GWAS) in Santa Inês sheep, one of the principal sheep breeds in Brazil, have not yet been reported. Thus, the aim of this study was to investigate the association between CNV and growth, efficiency and carcass traits in sheep. The Illumina OvineSNP50 BeadChip array was used to detect CNV in 491 Santa Inês individuals. Then, CNV-based GWAS was performed with a linear mixed model approach considering a genomic relationship matrix, for ten traits: (1) growth: body weight at three (W3) and six (W6) months of age; (2) efficiency: residual feed intake (RFI) and feed efficiency (FE) and (3) carcass: external carcass length (ECL), leg length (LL), carcass yield (CY), commercial cuts weight (CCW), loin eye area (LEA) and subcutaneous fat thickness (SFT). We identified 1,167 autosomal CNV in 438 sheep, with 294 non-redundant CNV, ranging from 21.8 to 861.9 kb, merged into 216 distinct copy number variation regions (CNVRs). One significant CNV segment (p_FDR -value<0.05) in OAR3 was associated with CY, while another significant CNV in OAR6 was associated with RFI. Additionally, another 5 CNV segments were considered relevant for investigation in the future studies. The significant segments overlapped 4 QTLs and spanned 8 genes, including the SPAST, TGFA and ADGRL3 genes, involved in cell differentiation and energy metabolism. Therefore, the results of the present study increase knowledge about CNV in sheep, their possible impacts on productive traits, and provide information for future investigations, being especially useful for those interested in structural variations in the sheep genome.

Genetic effects of heat stress on milk fatty acids in a Brazilian Holstein cattle

The present study aimed to estimate covariance components of milk fatty acids (FA) and to compare the genomic estimated breeding values under general and heat-stress effects. Data consisted of 38,762 test-day records from 6,344 Holstein cows obtained from May 2012 through January 2018 on 4 dairy herds from Brazil. Single-trait repeatability test-day models with random regressions as a function of temperature-humidity index values were used for genetic analyses. The models included contemporary groups, parity order (1-6), and days in milk classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, differences in heritability estimates between environments (general and heat stress) increased (0.03 to 0.06) for unsaturated FA traits, such as unsaturated, monounsaturated, and polyunsaturated, at higher heat-stress levels. In contrast, heritability estimated between environments for saturated FA traits, including saturated FA, palmitic acid (C16:0), and stearic acid (C18:0) did not observe significant differences between environments. In addition, our study revealed negative genetic correlations between general and heat-stress additive genetic effects (antagonistic effect) for the saturated FA, C16:0, C18:0, and C18:1, which ranged from -0.007 to -0.32. Spearman's ranking correlation between genomic estimated breeding values ranged from -0.27 to 0.99. Results indicated a moderate to strong interaction of genotype by the environment for most FA traits comparing a heat-stress environment with thermoneutral conditions. Our findings point out novel opportunities to explore the use of FA milk profile and heat-stress models.

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.

Liquid Biopsy as a Diagnostic and Prognostic Tool for Women and Female Dogs with Breast Cancer

INTRODUCTION

Breast cancer (BC) is the malignant neoplasm with the highest mortality rate in women and female dogs are good models to study BC.

OBJECTIVE

We investigated the efficacy of liquid biopsy to detect gene mutations in the diagnosis and follow-up of women and female dogs with BC.

MATERIALS AND METHODS

In this study, 57 and 37 BC samples were collected from women and female dogs, respectively. After core biopsy and plasma samples were collected, the DNA and ctDNA of the tumor fragments and plasma were processed for next generation sequencing (NGS) assay. After preprocessing of the data, they were submitted to the Genome Analysis ToolKit (GATK).

RESULTS

In women, 1788 variants were identified in tumor fragments and 221 variants in plasma; 66 variants were simultaneously detected in tumors and plasma. Conversely, in female dogs, 1430 variants were found in plasma and 695 variants in tumor fragments; 59 variants were simultaneously identified in tumors and plasma. The most frequently mutated genes in the tumor fragments of women were USH2A, ATM, and IGF2R; in female dogs, they were USH2A, BRCA2, and RRM2. Plasma of women showed the most frequent genetic variations in the MAP3K1, BRCA1, and GRB7 genes, whereas plasma from female dogs had variations in the NF1, ERBB2, and KRT17 genes. Mutations in the AKT1, PIK3CA, and BRIP genes were associated with tumor recurrence, with a highly pathogenic variant in PIK3CA being particularly prominent. We also detected a gain-of-function mutation in the GRB7, MAP3K1, and MLH1 genes.

CONCLUSION

Liquid biopsy is useful to identify specific genetic variations at the beginning of BC manifestation and may be accompanied over the entire follow-up period, thereby supporting the clinicians in refining interventions.

SAP30 Gene Is a Probable Regulator of Muscle Hypertrophy in Chickens

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

Does algae β-glucan affect the fecal bacteriome in dairy calves?

β-glucans has been reported to be associated with many health-promoting and improvements in animal performance, however, information about their effects on the bacterial community remains unknown. This study aimed to investigate how the addition of β-glucans can affect the fecal bacterial community with possible consequences on animal growth and health. For this, newborn Holstein calves (n = 14) were individually housed in tropical shelters and blocked according to sex, date, and weight at birth and randomly assigned to one of the following treatments: (1) Control: milk replacer (14% solids, 24% CP, 18.5% fat); (2) β-glucans: milk replacer supplemented with β-glucans (2 g/d). All calves were bucket fed 6 L/d of milk replacer and received water and starter concentrate ad libitum starting on d 2. To evaluate the bacteriome, fecal samples were collected at weeks 1, 2, 4, and 8. The bacterial community was assessed through sequencing of the V3-V4 region of the 16S rRNA gene on the Illumina MiSeq platform and analyzed using the DADA2 pipeline. No differences for Shannon and Chao1 indexes were observed for treatments, but both indexes increased with age (P < 0.001). There were dissimilarities in the structure of the bacterial community during the pre-weaning period (P = 0.01). In a deeper taxonomic level, Collinsella (Actinobacteriota), Prevotella (Bacteroidota), and Lactobacillus (Firmicutes) were the most abundant genera (9.84, 9.54, and 8.82% of the sequences, respectively). β-glucans promoted a higher abundance of Alloprevotella and Holdemanella, which may indicate a beneficial effect of supplementation on dairy calves. The bacterial community was highly correlated with the fecal score at weeks 1 and 2 and with starter concentrate intake at week 8. In conclusion, algae β-glucan supplementation could be beneficial to fecal bacteriome and consequently to the health and performance of dairy calves.

Genome assembly and annotation of the tambaqui (Colossoma macropomum): an emblematic fish of the Amazon River Basin

Colossoma macropomum, known as "tambaqui", is the largest Characiformes fish in the Amazon River Basin and a leading species in Brazilian aquaculture and fisheries. Good quality meat and excellent adaptability to culture systems are some of its remarkable farming features. To support studies into the genetics and genomics of the tambaqui, we have produced the first high-quality genome for the species. We combined Illumina and PacBio sequencing technologies to generate a reference genome, assembled with 39× coverage of long reads and polished to a consensus quality value (QV) of 36 with 130× coverage of short reads. The genome was assembled into 1269 scaffolds (a total of 1,221,847,006 bases), with a scaffold N50 size of 40 Mb, where 93% of all assembled bases were placed in the largest 54 scaffolds corresponding to the diploid karyotype of the tambaqui. Furthermore, the NCBI Annotation Pipeline annotated genes, pseudogenes, and non-coding transcripts using the RefSeq database as evidence, guaranteeing a high-quality annotation. A Genome Data Viewer for the tambaqui was produced, which will benefit groups interested in exploring the unique genomic features of the species. The availability of a highly accurate genome assembly for tambaqui provides the foundation for the discovery of novel ecological and evolutionary insights, and is a helpful resource for aquaculture.

Ruminal and Fecal Bacteriome of Dairy Calves Fed Different Levels and Sources of NDF

A starter concentrate containing different levels and sources of NDF can modify the gastrointestinal bacteriome. This study evaluated 18 Holstein calves housed in un-bedded suspended individual cages, fed one of three treatments: 22NDF: a conventional starter containing 22% NDF (n = 7); 31NDF: a starter with 31% NDF, replacing part of the corn by soybean hull (n = 6); and 22Hay: diet 22NDF plus coast-cross hay ad libitum (n = 5). All animals received 4 L of milk replacer daily, weaned at 8th week of age, and housed in wood shelters until week 10. To evaluate the bacteriome, the bacterial community of ruminal fluid and fecal samples was determined by sequencing V3 and V4 region amplicons of the 16S rRNA gene. Bacterial diversity in rumen was not affected by diet or age. The phyla Firmicutes and Bacteroidota, and Prevotella' genus were the most abundant in ruminal fluid and fecal samples. In feces, the α-diversity indices were higher for 22Hay. All indices were significantly affected by age. We believe that the ruminal bacteriome was affected by basal diet components, but not affected by NDF levels or sources. The supply of hay was effective in modifying the fecal bacteriome of dairy calves due to hind gut fermentation.

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.

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

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

Liquid Biopsy Can Detect Brca2 Gene Variants In Female Dogs With Mammary Neoplasia

Mammary tumors (MT) is one of the most prevalent malignancies in female dogs and women. Currently, molecular analyzes have shown that each tumor type presents its own genetic signature. In this context, liquid biopsy allows a comprehensive genetic characterization of the tumor, enabling early diagnosis and personalized treatment of patients. In women, deleterious mutations inherited in BRCA2 gene are associated with an increased risk of breast cancer, resistance to therapies, and worse prognosis. In female dogs, there are many divergent data on the involvement of BRCA2 gene with mammary carcinogenesis and what its pathogenic potential is. Therefore, the objective was to identify BRCA2 gene variants in 20 plasma DNA samples, from 10 newly diagnosed dogs with mammary cancer (RD), five control (CTR), and five mastectomized (MAST) patients. Eleven single nucleotide polymorphisms (SNPs) were detected, most of them in the exon 11, and two indels (deletion / insertion) in the BRCA2 gene. However, there was no statistically significant difference in the SNPs/indels detected between the groups. In addition, only one SNP (p.T1425P) and one deletion (p.L2307del) were considered deleterious using in silico computational models. Interestingly, most common variants were present in the plasma of all groups, except for the Ile2614Thr, Ile2614Val, Thr1425Pro, and p.L2307del variants. Thus, we observed that SNPs are common in the BRCA2 gene of female dogs with mammary tumors, with a similar condition identified in women with breast cancer. Liquid biopsy approach in dogs with MT is useful for genetic and therapeutic proposals. This article is protected by copyright. All rights reserved.

Comparative untargeted metabolome analysis of ruminal fluid and feces of Nelore steers (Bos indicus)

We conducted a study to identify the fecal metabolite profile and its proximity to the ruminal metabolism of Nelore steers based on an untargeted metabolomic approach. Twenty-six Nelore were feedlot with same diet during 105 d. Feces and rumen fluid were collected before and at slaughter, respectively. The metabolomics analysis indicated 49 common polar metabolites in the rumen and feces. Acetate, propionate, and butyrate were the most abundant polar metabolites in both bio-samples. The rumen presented significantly higher concentrations of the polar compounds when compared to feces (P < 0.05); even though, fecal metabolites presented an accentuated representability of the ruminal fluid metabolites. All fatty acids present in the ruminal fluid were also observed in the feces, except for C20:2n6 and C20:4n6. The identified metabolites offer information on the main metabolic pathways (higher impact factor and P < 0.05), as synthesis and degradation of ketone bodies; the alanine, aspartate and glutamate metabolisms, the glycine, serine; and threonine metabolism and the pyruvate metabolism. The findings reported herein on the close relationship between the ruminal fluid and feces metabolic profiles may offer new metabolic information, in addition to facilitating the sampling for metabolism investigation in animal production and health routines.