The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics.

BACKGROUND: The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí). RESULTS: A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions. CONCLUSIONS: Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

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.

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.

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.

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.

Transcriptome profile of skeletal muscle using different sources of dietary fatty acids in male pigs.

Pork is of great importance in world trade and represents the largest source of fatty acids in the human diet. Lipid sources such as soybean oil (SOY), canola (CO), and fish oil (FO) are used in pig diets and influence blood parameters and the ratio of deposited fatty acids. In this study, the main objective was to evaluate changes in gene expression in porcine skeletal muscle tissue resulting from the dietary oil sources and to identify metabolic pathways and biological process networks through RNA-Seq. The addition of FO in the diet of pigs led to intramuscular lipid with a higher FA profile composition of C20:5 n-3, C22:6 n-3, and SFA (C16:0 and C18:0). Blood parameters for the FO group showed lower cholesterol and HDL content compared with CO and SOY groups. Skeletal muscle transcriptome analyses revealed 65 differentially expressed genes (DEG, FDR 10%) between CO vs SOY, and 32 DEG for CO vs FO, and 531 DEG for SOY vs FO comparison. Several genes, including AZGP1, PDE3B, APOE, PLIN1, and LIPS, were found to be down-regulated in the diet of the SOY group compared to the FO group. The enrichment analysis revealed DEG involved in lipid metabolism, metabolic diseases, and inflammation between the oil groups, with specific gene functions in each group and altered blood parameters. The results provide mechanisms to help us understand the behavior of genes according to fatty acids.

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.

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.

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 (h2 = 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 [rg = -0.90 (L* and a*); rg = -0.64 (L* and b*); and rg = 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.

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.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers.

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

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.

Variants in myostatin and MyoD family genes are associated with meat quality traits in Santa Inês sheep.

Myostatin and MyoD family genes play vital roles in myogenesis and this study aimed to identify association of variants in MyoD1, MyoG, MyF5, MyF6, and MSTN genes with meat traits in Santa Inês sheep. A dataset with 44 variants and records of seven meat traits in 192 lambs (pH0, pH24, a*, b*, L*, tenderness assessed by shear force, and water-holding capacity) was used. Single-locus and haplotype association analyses were performed, and the significance threshold was established according to Bonferroni's method. Single-locus analysis revealed two associations at a Bonferroni level, where the variant c.935-185C > G in MyoD1 had an additive effect (-4.31 ± 1.08 N) on tenderness, while the variant c.464 + 185G > A in MyoG had an additive effect (-2.86 ± 0.64) on a*. Additionally, the haplotype replacement GT>AC in MSTN was associated with pH0 (1.26 ± 0.31), pH24 (1.07 ± 0.27), a* (-1.40 ± 0.51), and tenderness (3.83 ± 1.22 N), while the replacement GT > AG in MyoD1 was associated with pH0 (1.43 ± 0.26), pH24 (1.25 ± 0.22), b* (-1.06 ± 0.39), and tenderness (-4.13 ± 1.16 N). Our results have demonstrated that some variants in MyoG, MyF6, MyoD1, and MSTN can be associated with physicochemical meat traits in Santa Inês sheep.

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 (pFDR -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.

Genome-wide detection of CNVs and their association with performance traits in broilers.

BACKGROUND: Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age. RESULTS: We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels. CONCLUSIONS: Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Interplay among miR-29 family, mineral metabolism, and gene regulation in Bos indicus muscle.

An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.\\An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers.

BACKGROUND: Bacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens. RESULTS: Our data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB. CONCLUSIONS: Our data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

Identification of selection signatures involved in performance traits in a paternal broiler line.

BACKGROUND: Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits. Selection signatures may be identified using different methodologies, of which some are based on detecting contiguous sequences of homozygous identical-by-descent haplotypes, called runs of homozygosity (ROH), or estimating fixation index (FST) of genomic windows that indicates genetic differentiation. This study aimed to identify selection signatures in a paternal broiler TT line at generations 7th and 16th of selection and to investigate the genes annotated in these regions as well as the biological pathways involved. For such purpose, ROH and FST-based analysis were performed using whole genome sequence of twenty-eight chickens from two different generations. RESULTS: ROH analysis identified homozygous regions of short and moderate size. Analysis of ROH patterns revealed regions commonly shared among animals and changes in ROH abundance and size between the two generations. Results also suggest that whole genome sequencing (WGS) outperforms SNPchip data avoiding overestimation of ROH size and underestimation of ROH number; however, sequencing costs can limited the number of animals analyzed. FST-based analysis revealed genetic differentiation in several genomic windows. Annotation of the consensus regions of ROH and FST windows revealed new and previously identified genes associated with traits of economic interest, such as APOB, IGF1, IGFBP2, POMC, PPARG, and ZNF423. Over-representation analysis of the genes resulted in biological terms of skeletal muscle, matrilin proteins, adipose tissue, hyperglycemia, diabetes, Salmonella infections and tyrosine. CONCLUSIONS: Identification of ROH and FST-based analyses revealed selection signatures in TT line and genes that have important role in traits of economic interest. Changes in the genome of the chickens were observed between the 7th and 16th generations showing that ancient and recent selection in TT line may have acted over genomic regions affecting diseases and performance traits.

Proximal femoral head transcriptome reveals novel candidate genes related to epiphysiolysis in broiler chickens.

BACKGROUND: The proximal femoral head separation (FHS) or epiphysiolysis is a prevalent disorder affecting the chicken femur epiphysis, being considered a risk factor to infection which can cause bacterial chondronecrosis with osteomyelitis in broilers. To identify the genetic mechanisms involved in epiphysiolysis, differentially expressed (DE) genes in the femur of normal and FHS-affected broilers were identified using RNA-Seq technology. Femoral growth plate (GP) samples from 35-day-old commercial male broilers were collected from 4 healthy and 4 FHS-affected broilers. Sequencing was performed using an Illumina paired-end protocol. Differentially expressed genes were obtained using the edgeR package based on the False Discovery Rate (FDR < 0.05). RESULTS: Approximately 16 million reads/sample were generated with 2 × 100 bp paired-end reads. After data quality control, approximately 12 million reads/sample were mapped to the reference chicken genome (Galgal5). A total of 12,645 genes were expressed in the femur GP. Out of those, 314 were DE between groups, being 154 upregulated and 160 downregulated in FHS-affected broilers. In the functional analyses, several biological processes (BP) were overrepresented. Among them, those related to cell adhesion, extracellular matrix (ECM), bone development, blood circulation and lipid metabolism, which are more related to chicken growth, are possibly involved with the onset of FHS. On the other hand, BP associated to apoptosis or cell death and immune response, which were also found in our study, could be related to the consequence of the FHS. CONCLUSIONS: Genes with potential role in the epiphysiolysis were identified through the femur head transcriptome analysis, providing a better understanding of the mechanisms that regulate bone development in fast-growing chickens. In this study, we highlighted the importance of cell adhesion and extracellular matrix related genes in triggering FHS. Furthermore, we have shown new insights on the involvement of lipidemia and immune response/inflammation with FHS in broilers. Understanding the changes in the GP transcriptome might support breeding strategies to address poultry robustness and to obtain more resilient broilers.

Didelphis albiventris: an overview of unprecedented transcriptome sequencing of the white-eared opossum.

BACKGROUND: The white-eared opossum (Didelphis albiventris) is widely distributed throughout Brazil and South America. It has been used as an animal model for studying different scientific questions ranging from the restoration of degraded green areas to medical aspects of Chagas disease, leishmaniasis and resistance against snake venom. As a marsupial, D. albiventris can also contribute to the understanding of the molecular mechanisms that govern the different stages of organogenesis. Opossum joeys are born after only 13 days, and the final stages of organogenesis occur when the neonates are inside the pouch, depending on lactation. As neither the genome of this opossum species nor its transcriptome has been completely sequenced, the use of D. albiventris as an animal model is limited. In this work, we sequenced the D. albiventris transcriptome by RNA-seq to obtain the first catalogue of differentially expressed (DE) genes and gene ontology (GO) annotations during the neonatal stages of marsupial development. RESULTS: The D. albiventris transcriptome was obtained from whole neonates harvested at birth (P0), at 5 days of age (P5) and at 10 days of age (P10). The de novo assembly of these transcripts generated 85,338 transcripts. Approximately 30% of these transcripts could be mapped against the amino acid sequences of M. domestica, the evolutionarily closest relative of D. albiventris to be sequenced thus far. Among the expressed transcripts, 2077 were found to be DE between P0 and P5, 13,780 between P0 and P10, and 1453 between P5 and P10. The enriched GO terms were mainly related to the immune system, blood tissue development and differentiation, vision, hearing, digestion, the CNS and limb development. CONCLUSIONS: The elucidation of opossum transcriptomes provides an out-group for better understanding the distinct characteristics associated with the evolution of mammalian species. This study provides the first transcriptome sequences and catalogue of genes for a marsupial species at different neonatal stages, allowing the study of the mechanisms involved in organogenesis.