Molecular evolution and signatures of selective pressures on Bos, focusing on the Nelore breed (Bos indicus).

Evolutionary history leads to genome changes over time, especially for species that have experienced intense selective pressures over a short period. Here, we investigated the genomic evolution of Bos species by searching for potential selection signatures, focusing on Nelore, an economically relevant cattle breed in Brazil. We assessed the genomic processes determining the molecular evolution across Nelore and thirteen other related taxa by evaluating (i) amino acid sequence conservation, (ii) the dN/dS ratio, and (iii) gene families' turnover rate (λ). Low conserved regions potentially associated with fatty acid metabolism seem to reflect differences in meat fat content in taxa with different evolutionary histories. All Bos species presented genes under positive selection, especially B. indicus and Nelore, which include transport protein cobalamin, glycolipid metabolism, and hormone signaling. These findings could be explained by constant selective pressures to obtain higher immune resistance and efficient metabolism. The gene contraction rate across the Nelore + B. indicus branch was almost nine times higher than that in other lineages (λ = 0.01043 vs. 0.00121), indicating gene losses during the domestication process. Amino acid biosynthesis, reproductive and innate immune system-related pathways were associated with genes recognized within the most frequent rapidly evolving gene families and in genes under positive selection, supporting the substantial relevance of such traits from a domestication perspective. Our data provide new insights into how the genome may respond to intense artificial selection in distinct taxa, and reinforces the presence of selective pressures on traits potentially relevant for future animal breeding investments.

The Impact of Selective Dry Cow Therapy Adopted in a Brazilian Farm on Bacterial Diversity and the Abundance of Quarter Milk.

We aimed to evaluate the impact of selective dry cow therapy (SDCT) (protocol 1: antimicrobial combined with internal teat sealant (ITS); vs. protocol 2: ITS alone) on bacterial diversity and the abundance of quarter milk. Eighty high production cows (parity ≤ 3 and an average milk yield of 36.5 kg/cow/day) from the largest Brazilian dairy herd available were randomly selected; milk quarter samples were collected for microbiological culture (MC) on the day of drying-off (n = 313) and on day 7 post-calving (n = 313). Based on the results of the MC before and after calving, 240 quarters out of 313 were considered healthy, 38 were cured, 29 showed new infections and 6 had persistent infections. Mammary quarters were randomly selected based on intramammary information status and SDCT protocols for bacterial diversity analyses. The bacterial diversity was similar when comparing both healthy and cured quarters submitted to both drying-off protocols. Despite healthy cows that were treated at dry-off using only teat sealant showing no alteration in the alpha and beta bacterial diversity, they did show a higher abundance of bacterial groups that may be beneficial to or commensals of the mammary gland, which implies that antibiotic therapy should be reserved for mammary quarters with a history of mastitis.

Stool and Ruminal Microbiome Components Associated With Methane Emission and Feed Efficiency in Nelore Beef Cattle.

Background: The impact of extreme changes in weather patterns on the economy and human welfare is one of the biggest challenges our civilization faces. From anthropogenic contributions to climate change, reducing the impact of farming activities is a priority since it is responsible for up to 18% of global greenhouse gas emissions. To this end, we tested whether ruminal and stool microbiome components could be used as biomarkers for methane emission and feed efficiency in bovine by studying 52 Brazilian Nelore bulls belonging to two feed intervention treatment groups, that is, conventional and by-product-based diets. Results: We identified a total of 5,693 amplicon sequence variants (ASVs) in the Nelore bulls' microbiomes. A Differential abundance analysis with the ANCOM approach identified 30 bacterial and 15 archaeal ASVs as differentially abundant (DA) among treatment groups. An association analysis using Maaslin2 software and a linear mixed model indicated that bacterial ASVs are linked to the host's residual methane emission (RCH4) and residual feed intake (RFI) phenotype variation, suggesting their potential as targets for interventions or biomarkers. Conclusion: The feed composition induced significant differences in both abundance and richness of ruminal and stool microbial populations in ruminants of the Nelore breed. The industrial by-product-based dietary treatment applied to our experimental groups influenced the microbiome diversity of bacteria and archaea but not of protozoa. ASVs were associated with RCH4 emission and RFI in ruminal and stool microbiomes. While ruminal ASVs were expected to influence CH4 emission and RFI, the relationship of stool taxa, such as Alistipes and Rikenellaceae (gut group RC9), with these traits was not reported before and might be associated with host health due to their link to anti-inflammatory compounds. Overall, the ASVs associated here have the potential to be used as biomarkers for these complex phenotypes.

Genomic Signatures of Divergent Ecological Strategies in a Recent Radiation of Neotropical Wild Cats.

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here, we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological, and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, for example a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

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

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

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation.

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Immune status, well-being and gut microbiota in military supplemented with synbiotic ice cream and submitted to field training: a randomised clinical trial.

Strenuous physical activity, sleep deprivation and psychological stress are common features of military field training. The present study aimed to evaluate the effects of supplementation with a synbiotic ice cream on salivary IgA, gastrointestinal symptoms, well-being indicators and gut microbiota in young military participants undergoing field training. Sixty-five military completed the study: one group was supplemented for 30 d with synbiotic ice cream containing: 2·1 × 108 CFU/g for Lactobacillus acidophilus LA-5 and 2·7 × 109 CFU/g for Bifidobacterium animalis BB-12 and 2·3 g of inulin in the 60 g of ice cream at manufacture, and the other with a placebo ice cream. Volunteers were evaluated at pre-supplementation (baseline), post-supplementation and after a 5-d military training. Bifidobacterium and Lactobacillus genera were measured in stool samples and both showed a higher differential abundance post-supplementation and training. Salivary IgA and gastrointestinal symptoms decreased at post-training in both groups (P < 0·05; main effect of time); however, supplementation with synbiotic did not mitigate this effect. Tenseness and sleepiness were decreased in the synbiotic-treated group, but not in the placebo group at post-military training (P = 0·01 and 0·009, respectively; group × time effect). The other well-being indicators were not affected by the synbiotic supplementation. In conclusion, 30 d of synbiotic ice cream supplementation containing inulin, L. acidophilus LA-5 and B. animalis BB-12 favourably modulated gut microbiota and improved tenseness and sleepiness in healthy young military undergoing a 5-d field training. These improvements may be relevant to this population as they may influence the decision-making process in an environment of high physical and psychological stress.

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

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

Proteome alterations associated with the oleic acid and cis-9, trans-11 conjugated linoleic acid content in bovine skeletal muscle.

Oleic acid (OA) and cis-9, trans-11 conjugated linoleic acid (c9t11-CLA) are fatty acids found in beef with beneficial effects in human health. This study investigated differentially abundant proteins (DAPs) in skeletal muscle of bovines with extreme values of OA, and c9t11-CLA. For each one of the fatty acids, twenty muscle samples were divided into two groups (N = 10_High; N = 10_Low) and analyzed by high definition mass spectrometry. We identified 103 and 133 DAPs between the groups for each fatty acid. We found 64 and 45 up-regulated and 39 and 68 down-regulated proteins for OA and c9t11-CLA, respectively. Comparative analysis between proteomic and transcriptomic data revealed eight and ten genes with a consistent between mRNA expression levels and protein abundance for OA and c9t11-CLA, respectively. Unconventional myosin-Id (MYO1D), mineralocorticoid receptor (NR3C2), geranylgeranyl transferase type-2 subunit-alpha (RABGGTA), and uveal autoantigen with coiled-coil domains and ankyrin repeats (UACA) were found as putative candidate proteins for OA content. Fatty acid synthase (FASN), tubulin alpha-4A chain (TUBA4A), vinculin (VCL), NADH dehydrogenase 1 alpha subcomplex 5 (NDUFA5), and prefoldin subunit 6 (PFDN6) for c9t11-CLA. Our findings contribute to a deeper understanding of the molecular mechanisms behind the regulation of the OA and c9t11-CLA content in cattle skeletal muscle. SIGNIFICANCE: Questions about the association between meat intake and disease incidence in humans has driven animal scientist to pursue a better understanding of the biological processes associated with differences in the intramuscular fat composition. The beneficial effects of oleic acid and conjugated linoleic acid in human health have been demonstrated by improving the immune system and preventing atherosclerosis, different types of cancers, hypertension, and diabetes. Previous genome-wide association and gene expression studies identified genomic regions and differentially expressed genes associated with the fatty acid profile in skeletal muscle. In this work, differences were evaluated at the protein level. The use of a label-free quantitative proteomic approach, compared with muscle transcriptome results obtained by RNA-sequencing, allowed us to earn new insights into the variability in fatty acid deposition in skeletal muscle of farm animals. This study opens new avenues to explore the effect of the fatty acids in the skeletal muscle of livestock animals, which is associated with nutritional values of the meat, and perhaps to understand the mechanisms correlated with metabolic diseases in other species.

The use of an alternative feed additive, containing benzoic acid, thymol, eugenol, and piperine, improved growth performance, nutrient and energy digestibility, and gut health in weaned piglets.

This research evaluated a feed additive (benzoic acid, eugenol, thymol, and piperine), associated or not with colistin, in weaned piglets feeding. The parameters evaluated were growth performance, apparent total tract digestibility (ATTD) of nutrients, diarrhea incidence, intestinal morphology, relative weights of digestive organs, microbial diversity, and the percentages of operational taxonomic units of microorganisms in the cecum content of pigs. One-hundred and eight crossbred piglets (5.3 ± 0.5 kg) were used in a three-phase feeding program (21 to 35, 36 to 50, 51 to 65 d of age) and fed a control diet with no inclusion of growth promoter feed additive, a diet with 40 ppm of colistin, a diet with 0.3% of alternative additive, and a diet with 0.3% of alternative additive and 40 ppm of colistin. The diets were based on corn, soybean meal, dairy products, and spray-dried blood plasma and formulated to provide 3.40, 3.38, and 3.20 Mcal of ME/kg and 14.5, 13.3, and 10.9 g/kg of digestible lysine, in phases 1, 2, and 3, respectively. The piglets were housed three per pen, with nine replicates per diet, in a complete randomized block design based on initial BW. The data were submitted to ANOVA and means were separated by Tukey test (5%), using SAS. Pigs fed diets with the alternative feed additive had greater (P < 0.05) ADG (114.3 vs. 91.8 g) and ADFI (190.1 vs. 163.3 g) in phase 1 than pigs fed diets without the product. The alternative additive improved (P < 0.05) ATTD of crude protein (CP) in phase 1 (71.0% vs. 68.6%), gross energy in phases 1 (77.4% vs. 75.2%) and 3 (79.0% vs. 77.1%), and dry matter in phase 3 (79.1% vs. 77.1%). The antibiotic inclusion in the diets increased (P < 0.05) ATTD of CP in phase 1 (71.5% vs. 68.2%). The alternative feed additive tended (P = 0.06) to increase (46%) normal feces frequency, decreased (P < 0.05) goblet cells count (104.3 vs. 118.1) in the jejunum, and decreased (P < 0.05) small intestine (4.60% vs. 4.93%) and colon (1.41% vs. 1.65%) relative weights, compared with pigs not fed with the alternative additive. There was a tendency (P = 0.09) for a lower concentration of Escherichia-Shigella (1.46% vs. 3.5%) and lower (P < 0.05) percentage of Campylobacter (0.52% vs. 10.21%) in the cecum content of piglets fed diets containing essential oils and benzoic acid compared with pigs fed diets without the alternative feed additive. The alternative feed additive was effective in improving growth performance, diets digestibility, and gut health in piglets soon after weaning.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses.

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

The structure of microbial populations in Nelore GIT reveals inter-dependency of methanogens in feces and rumen.

BACKGROUND: The success of different species of ruminants in the colonization of a diverse range of environments is due to their ability to digest and absorb nutrients from cellulose, a complex polysaccharide found in leaves and grass. Ruminants rely on a complex and diverse microbial community, or microbiota, in a unique compartment known as the rumen to break down this polysaccharide. Changes in microbial populations of the rumen can affect the host's development, health, and productivity. However, accessing the rumen is stressful for the animal. Therefore, the development and use of alternative sampling methods are needed if this technique is to be routinely used in cattle breeding. To this end, we tested if the fecal microbiome could be used as a proxy for the rumen microbiome due to its accessibility. We investigated the taxonomic composition, diversity and inter-relations of two different GIT compartments, rumen and feces, of 26 Nelore (Bos indicus) bulls, using Next Generation Sequencing (NGS) metabarcoding of bacteria, archaea and ciliate protozoa. RESULTS: We identified 4265 Amplicon Sequence Variants (ASVs) from bacteria, 571 from archaea, and 107 from protozoa, of which 143 (96 bacteria and 47 archaea) were found common between both microbiomes. The most prominent bacterial phyla identified were Bacteroidetes (41.48%) and Firmicutes (56.86%) in the ruminal and fecal microbiomes, respectively, with Prevotella and Ruminococcaceae UCG-005 the most relatively abundant genera identified in each microbiome. The most abundant archaeal phylum identified was Euryarchaeota, of which Methanobrevibacter gottschalkii, a methanogen, was the prevalent archaeal species identified in both microbiomes. Protozoa were found exclusively identified in the rumen with Bozasella/Triplumaria being the most frequent genus identified. Co-occurrence among ruminal and fecal ASVs reinforces the relationship of microorganisms within a biological niche. Furthermore, the co-occurrence of shared archaeal ASVs between microbiomes indicates a dependency of the predominant fecal methanogen population on the rumen population. CONCLUSIONS: Co-occurring microorganisms were identified within the rumen and fecal microbiomes, which revealed a strong association and inter-dependency between bacterial, archaeal and protozoan populations of the same microbiome. The archaeal ASVs identified as co-occurring between GIT compartments corresponded to the methanogenic genera Methanobrevibacter and Methanosphaera and represented 26.34% of the overall archaeal sequencesdiversity in the rumen and 42.73% in feces. Considering that these archaeal ASVs corresponded to a significant part of the overall diversity of both microbiomes, which is much higher if one includes the interactions of these co-occurring with other rumen archaea ASVs, we suggest that fecal methanogens could be used as a proxy of ruminal methanogens.

Gene expression profile suggesting immunological dysregulation in two Brazilian Bloom's syndrome cases.

BACKGROUND: Bloom syndrome (BS) is a rare autosomal recessive chromosome instability disorder. The main clinical manifestations are growth deficiency, telangiectasic facial erythema, immunodeficiency, and increased risk to develop neoplasias at early age. Cytogenetic test for sister chromatid exchanges (SCEs) is used as a diagnostic marker for BS. In addition, most patients also present mutations in the BLM gene, related to defects in the DNA repair mechanism. However, the molecular mechanism behind the pathogenicity of BS is still not completely understood. METHODS: We describe two patients confirmed with BS by SCE and molecular analysis. Also, we performed the gene expression profile by the RNA-seq methodology in mRNA transcripts for differential gene expression analysis using as a biological condition for comparison BS versus health controls. RESULTS: We detected 216 differentially expressed genes related to immunological pathways such as positive regulation and activation of B cells, immune effector process and absence of difference of DNA repair genes expression. In addition; we also observed differentially expressed genes associated with apoptosis control, such as BCL2L1, CASP7, CDKN1A, E2F2, ITPR, CD274, TNFAIP6, TNFRSF25, TNFRSF13C, and TNFRSF17. CONCLUSION: Our results suggest that the combination of altered expression of genes involved in signaling pathways of immune response and apoptosis control may contribute directly to the main characteristics observed in BS, such as recurrent infections, growth failure, and high risk of cancer. Transcriptome studies of other instability syndromes could allow a more accurate analysis of the relevant gene interactions associated with the destabilization of the genome. This is a first description of the profile of differential gene expression related to immunological aspects detected in patients with BS by RNA-seq.

Role of neuroendocrine modulation and biochemistry in the sepsis in Piaractus mesopotamicus.

Sepsis is a systemic process with multifactorial pathophysiology that affects most animal species. It is responsible for high rates of morbidity and mortality. This work aimed to study the biochemical and neuroendocrine changes of the sepsis process in Piaractus mesopotamicus after Aeromonas hydrophila inoculation analyzing changes in blood leukocyte and differences in neuroendocrine-biochemical modulation using RNA-seq. Fish showed hypercortisolemia, inhibition of glucose absorption, followed by hypocortisolemia and then hyperglycemia. Thyroid hormones (T3 and T4) showed immediate decrease in serum and T4 increased 6 h post-inoculation (HPI). Sepsis-induced hormonal alterations triggered changes in the metabolic pathways increasing protein and lipid catabolism, use of transient anaerobic glycolysis and liver injury. A reference transcriptome was constructed based on blood leukocytes from P. mesopotamicus. The assembly resulted in total 266,272 contigs with a N50 of 2786 bp. There was a reorganization of plasma membrane of leukocytes at the beginning of the septic process with increased expression of neuroendocrine receptors and with continuous flow of neurotransmitters, hormones and solutes with compensatory regulation at 6 HPI. Three and nine HPI seemed to be critical, the expression of a number of transcription factors was increased, including the modulatory DEGs related to glucocorticoid and thyroid hormones induced and suppressed (FDR < 0.05). Neuroendocrine modulation can regulate leukocytes and biochemical parameters of peripheral blood, being important sources for the study of the pathophysiology of sepsis. These finding highlights the importance of further studies focusing on biochemical-neuroendocrine changes in blood leukocytes and systemic sepsis.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin.

BACKGROUND: Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance. RESULTS: Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively. CONCLUSION: This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development.

Sexual dimorphism is a relevant factor in animal science, since it can affect the gene expression of economically important traits. Eventually, the interest in the prenatal phase in a transcriptome study may not comprise the period of development in which male and female conceptuses are phenotypically divergent. Therefore, it would be interesting if sex differentiation could be performed using transcriptome data, with no need for extra techniques. In this study, the sex of pig conceptuses (embryos at 25 days-old and fetuses at 35 days-old) was determined by reads counts per million (CPM) of Y chromosome-linked genes that were discrepant among samples. Thus, ten genes were used: DDX3Y, KDM5D, ZFY, EIF2S3Y, EIF1AY, LOC110255320, LOC110257894, LOC396706, LOC100625207, and LOC110255257. Conceptuses that presented reads CPM sum for these genes (ΣCPMchrY) greater than 400 were classified as males and those with ΣCPMchrY below 2 were classified as females. It was demonstrated that the sex identification can be performed at early stages of pig development from RNA-sequencing analysis of genes mapped on Y chromosome. Additionally, these results reinforce that sex determination is a mechanism conserved across mammals, highlighting the importance of using pigs as an animal model to study sex determination during human prenatal development.

Genome-wide association scan for QTL and their positional candidate genes associated with internal organ traits in chickens.

BACKGROUND: Poultry breeding programs have been focused on improvement of growth and carcass traits, however, this has resulted in correlated changes in internal organ weights and increased incidence of metabolic disorders. These disorders can affect feed efficiency or even cause death. We used a high density SNP array (600 K, Affymetrix) to estimate genomic heritability, perform genome-wide association analysis, and identify genomic regions and positional candidate genes (PCGs) associated with internal organ traits in an F2 chicken population. We integrated knowledge of haplotype blocks, selection signature regions and sequencing data to refine the list of PCGs. RESULTS: Estimated genomic heritability for internal organ traits in chickens ranged from low (LUNGWT, 0.06) to high (GIZZWT, 0.45). A total of 20 unique 1 Mb windows identified on GGA1, 2, 4, 7, 12, 15, 18, 19, 21, 27 and 28 were significantly associated with intestine length, and weights or percentages of liver, gizzard or lungs. Within these windows, 14 PCGs were identified based on their biological functions: TNFSF11, GTF2F2, SPERT, KCTD4, HTR2A, RB1, PCDH7, LCORL, LDB2, NR4A2, GPD2, PTPN11, ITGB4 and SLC6A4. From those genes, two were located within haplotype blocks and three overlapped with selection signature regions. A total of 13,748 annotated sequence SNPs were in the 14 PCGs, including 156 SNPs in coding regions (124 synonymous, 26 non-synonymous, and 6 splice variants). Seven deleterious SNPs were identified in TNFSF11, NR4A2 or ITGB4 genes. CONCLUSIONS: The results from this study provide novel insights to understand the genetic architecture of internal organ traits in chickens. The QTL detection performed using a high density SNP array covered the whole genome allowing the discovery of novel QTL associated with organ traits. We identified PCGs within the QTL involved in biological processes that may regulate internal organ growth and development. Potential functional genetic variations were identified generating crucial information that, after validation, might be used in poultry breeding programs to reduce the occurrence of metabolic disorders.

Co-Expression Networks Reveal Potential Regulatory Roles of miRNAs in Fatty Acid Composition of Nelore Cattle.

Fatty acid (FA) content affects the sensorial and nutritional value of meat and plays a significant role in biological processes such as adipogenesis and immune response. It is well known that, in beef, the main FAs associated with these biological processes are oleic acid (C18:1 cis9, OA) and conjugated linoleic acid (CLA-c9t11), which may have beneficial effects on metabolic diseases such as type 2 diabetes and obesity. Here, we performed differential expression and co-expression analyses, weighted gene co-expression network analysis (WGCNA) and partial correlation with information theory (PCIT), to uncover the complex interactions between miRNAs and mRNAs expressed in skeletal muscle associated with FA content. miRNA and mRNA expression data were obtained from skeletal muscle of Nelore cattle that had extreme genomic breeding values for OA and CLA. Insulin and MAPK signaling pathways were identified by WGCNA as central pathways associated with both of these fatty acids. Co-expression network analysis identified bta-miR-33a/b, bta-miR-100, bta-miR-204, bta-miR-365-5p, bta-miR-660, bta-miR-411a, bta-miR-136, bta-miR-30-5p, bta-miR-146b, bta-let-7a-5p, bta-let-7f, bta-let-7, bta-miR 339, bta-miR-10b, bta-miR 486, and the genes ACTA1 and ALDOA as potential regulators of fatty acid synthesis. This study provides evidence and insights into the molecular mechanisms and potential target genes involved in fatty acid content differences in Nelore beef cattle, revealing new candidate pathways of phenotype modulation that could positively benefit beef production and human consumption.

Mutation dynamics of CpG dinucleotides during a recent event of vertebrate diversification.

DNA methylation in CpGs dinucleotides is associated with high mutability and disappearance of CpG sites during evolution. Although the high mutability of CpGs is thought to be relevant for vertebrate evolution, very little is known on the role of CpG-related mutations in the genomic diversification of vertebrates. Our study analysed genetic differences in chickens, between Red Junglefowl (RJF; the living closest relative to the ancestor of domesticated chickens) and domesticated breeds, to identify genomic dynamics that have occurred during the process of their domestication, focusing particularly on CpG-related mutations. Single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) between RJF and these domesticated breeds were assessed in a reduced fraction of their genome. Additionally, DNA methylation in the same fraction of the genome was measured in the sperm of RJF individuals to identify possible correlations with the mutations found between RJF and the domesticated breeds. Our study shows that although the vast majority of CpG-related mutations found relate to CNVs, CpGs disproportionally associate to SNPs in comparison to CNVs, where they are indeed substantially under-represented. Moreover, CpGs seem to be hotspots of mutations related to speciation. We suggest that, on the one hand, CpG-related mutations in CNV regions would promote genomic 'flexibility' in evolution, i.e., the ability of the genome to expand its functional possibilities; on the other hand, CpG-related mutations in SNPs would relate to genomic 'specificity' in evolution, thus, representing mutations that would associate with phenotypic traits relevant for speciation.

Detection of Co-expressed Pathway Modules Associated With Mineral Concentration and Meat Quality in Nelore Cattle.

Meat quality is a complex trait that is influenced by genetic and environmental factors, which includes mineral concentration. However, the association between mineral concentration and meat quality, and the specific molecular pathways underlying this association, are not well explored. We therefore analyzed gene expression as measured with RNA-seq in Longissimus thoracis muscle of 194 Nelore steers for association with three meat quality traits (intramuscular fat, meat pH, and tenderness) and the concentration of 13 minerals (Ca, Cr, Co, Cu, Fe, K, Mg, Mn, Na, P, S, Se, and Zn). We identified seven sets of co-expressed genes (modules) associated with at least two traits, which indicates that common pathways influence these traits. From pathway analysis of module hub genes, we further found an over-representation for energy and protein metabolism (AMPK and mTOR signaling pathways) in addition to muscle growth, and protein turnover pathways. Among the identified hub genes FASN, ELOV5, and PDE3B are involved with lipid metabolism and were affected by previously identified eQTLs associated to fat deposition. The reported hub genes and over-represented pathways provide evidence of interplay among gene expression, mineral concentration, and meat quality traits. Future studies investigating the effect of different levels of mineral supplementation in the gene expression and meat quality traits could help us to elucidate the regulatory mechanism by which the genes/pathways are affected.