In cattle, molecular control of oviduct receptivity to the embryo is poorly understood. Here, we used a bovine model for receptivity based on size of the pre-ovulatory follicle to compare oviductal ...global and candidate gene transcript abundance on day 4 of the estrous cycle. Growth of the pre-ovulatory follicle (POF) of Nelore (Bos indicus) cows was manipulated to produce two groups: large POF large corpus luteum (CL) group (LF-LCL; greater receptivity) and small POF-small CL group (SF-SCL). Oviductal samples were collected four days after GnRH-induced ovulation. Ampulla and isthmus transcriptome was obtained by RNA-seq, regional gene expression was assessed by qPCR, and PGR and ERa protein distribution was evaluated by immunohistochemistry. There was a greater abundance of PGR and ERa in the oviduct of LF-LCL animals thus indicating a greater availability of receptors and possibly sex steroids stimulated signaling in both regions. Transcriptomic profiles indicated a series of genes associated with functional characteristics of the oviduct that are regulated by the periovulatory sex steroid milieu and that potentially affect oviductal receptivity and early embryo development. They include tissue morphology changes (extra cellular matrix remodeling), cellular changes (proliferation), and secretion changes (growth factors, ions and metal transporters), and were enriched for the genes with increased expression in the LF-LCL group. In conclusion, differences in the periovulatory sex steroid milieu lead to different oviductal gene expression profiles that could modify the oviductal environment to affect embryo survival and development.
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.
Intramuscular fat (IMF) content is related to insulin resistance, which is an important prediction factor for disorders, such as cardiovascular disease, obesity and type 2 diabetes in human. At the ...same time, it is an economically important trait, which influences the sensorial and nutritional value of meat. The deposition of IMF is influenced by many factors such as sex, age, nutrition, and genetics. In this study Nellore steers (Bos taurus indicus subspecies) were used to better understand the molecular mechanisms involved in IMF content. This was accomplished by identifying differentially expressed genes (DEG), biological pathways and putative regulatory factors. Animals included in this study had extreme genomic estimated breeding value (GEBV) for IMF. RNA-seq analysis, gene set enrichment analysis (GSEA) and co-expression network methods, such as partial correlation coefficient with information theory (PCIT), regulatory impact factor (RIF) and phenotypic impact factor (PIF) were utilized to better understand intramuscular adipogenesis. A total of 16,101 genes were analyzed in both groups (high (H) and low (L) GEBV) and 77 DEG (FDR 10%) were identified between the two groups. Pathway Studio software identified 13 significantly over-represented pathways, functional classes and small molecule signaling pathways within the DEG list. PCIT analyses identified genes with a difference in the number of gene-gene correlations between H and L group and detected putative regulatory factors involved in IMF content. Candidate genes identified by PCIT include: ANKRD26, HOXC5 and PPAPDC2. RIF and PIF analyses identified several candidate genes: GLI2 and IGF2 (RIF1), MPC1 and UBL5 (RIF2) and a host of small RNAs, including miR-1281 (PIF). These findings contribute to a better understanding of the molecular mechanisms that underlie fat content and energy balance in muscle and provide important information for the production of healthier beef for human consumption.
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.
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.
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.
MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global ...microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult.
Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate.
These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.
Pregnancy success results from the interaction of multiple factors, among them are folliculogenesis and early embryonic development. Failure during these different processes can lead to difficulties ...in conception. Alternatives to overcome these problems are based on assisted reproductive techniques. Extracellular vesicles are cell-secreted vesicles present in different body fluids and contain bioactive materials, such as messenger RNA, microRNAs (miRNAs), and proteins. Thus, our hypothesis is that extracellular vesicles from follicular fluid from 3-6 mm ovarian follicles can modulate bovine embryo development in vitro. To test our hypothesis follicular fluid from bovine ovaries was aspirated and small-extracellular vesicles (<200 nm) were isolated for further analysis. Additionally, small-extracellular vesicles (EVs) were utilized for functional experiments investigating their role in modulating messenger RNA, microRNA as well as global DNA methylation and hydroxymethylation levels of bovine blastocysts. EVs from 3-6 mm follicles were used for RNA-seq and miRNA analysis. Functional annotation analysis of the EVs transcripts revealed messages related to chromatin remodeling and transcriptional regulation. EVs treatment during oocyte maturation and embryo development causes changes in blastocyst rates, as well as changes in the transcription levels of genes related to embryonic metabolism and development. Supplementation with EVs from 3-6 mm follicles during oocyte maturation and early embryo development (until the 4-cell stage) increased the levels of bta-miR-631 (enriched in EVs from 3-6 mm follicles) in embryos. Interestingly, the addition of EVs from 3-6 mm follicles induced changes in global DNA methylation and hydroxymethylation levels compared to embryos produced by the standard in vitro production system. Our results indicate that the supplementation of culture media with EVs isolated from the follicular fluid of 3-6 mm follicles during oocyte maturation and early embryo development can partially modify metabolic and developmental related genes as well as miRNA and global DNA methylation and hydroxymethylation, suggesting that EVs play an important role during oocyte maturation and early embryo development in vitro.
Positively correlated with carcass weight and animal growth, the ribeye area (REA) and the backfat thickness (BFT) are economic important carcass traits, which impact directly on producer's payment. ...The selection of these traits has not been satisfactory since they are expressed later in the animal's life and multigene regulated. So, next-generation technologies have been applied in this area to improve animal's selection and better understand the molecular mechanisms involved in the development of these traits. Correlation network analysis, performed by tools like WGCNA (Weighted Correlation Network Analysis), has been used to explore gene-gene interactions and gene-phenotype correlations. Thus, this study aimed to identify putative candidate genes and metabolic pathways that regulate REA and BFT by constructing a gene co-expression network using WGCNA and RNA sequencing data, to better understand genetic and molecular variations behind these complex traits in Nelore cattle.
The gene co-expression network analysis, using WGCNA, were built using RNA-sequencing data normalized by transcript per million (TPM) from 43 Nelore steers. Forty-six gene clusters were constructed, between them, three were positively correlated (p-value< 0.1) to the BFT (Green Yellow, Ivory, and Light Yellow modules) and, one cluster was negatively correlated (p-value< 0.1) with REA (Salmon module). The enrichment analysis performed by DAVID and WebGestalt (FDR 5%) identified eight Gene Ontology (GO) terms and three KEGG pathways in the Green Yellow module, mostly associated with immune response and inflammatory mechanisms. The enrichment of the Salmon module demonstrated 19 GO terms and 21 KEGG pathways, related to muscle energy metabolism, lipid metabolism, muscle degradation, and oxidative stress diseases. The Ivory and Light yellow modules have not shown significant results in the enrichment analysis.
With this study, we verified that inflammation and immune response pathways modulate the BFT trait. Energy and lipid metabolism pathways, highlighting fatty acid metabolism, were the central pathways associated with REA. Some genes, as RSAD2, EIF2AK2, ACAT1, and ACSL1 were considered as putative candidate related to these traits. Altogether these results allow us to a better comprehension of the molecular mechanisms that lead to muscle and fat deposition in bovine.
Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) ...in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or ;1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.
Commercial cuts yield is an important trait for beef production, which affects the final value of the products, but its direct determination is a challenging procedure to be implemented in practice. ...The measurement of ribeye area (REA) and backfat thickness (BFT) can be used as indirect measures of meat yield. REA and BFT are important traits studied in beef cattle due to their strong implication in technological (carcass yield) and nutritional characteristics of meat products, like the degree of muscularity and total body fat. Thus, the aim of this work was to study the Longissimus dorsi muscle transcriptome of Nellore cattle, associated with REA and BFT, to find differentially expressed (DE) genes, metabolic pathways, and biological processes that may regulate these traits.
By comparing the gene expression level between groups with extreme genomic estimated breeding values (GEBV), 101 DE genes for REA and 18 for BFT (false discovery rate, FDR 10%) were identified. Functional enrichment analysis for REA identified two KEGG pathways, MAPK (Mitogen-Activated Protein Kinase) signaling pathway and endocytosis pathway, and three biological processes, response to endoplasmic reticulum stress, cellular protein modification process, and macromolecule modification. The MAPK pathway is responsible for fundamental cellular processes, such as growth, differentiation, and hypertrophy. For BFT, 18 biological processes were found to be altered and grouped into 8 clusters of semantically similar terms. The DE genes identified in the biological processes for BFT were ACHE, SRD5A1, RSAD2 and RSPO3. RSAD2 has been previously shown to be associated with lipid droplet content and lipid biosynthesis.
In this study, we identified genes, metabolic pathways, and biological processes, involved in differentiation, proliferation, protein turnover, hypertrophy, as well as adipogenesis and lipid biosynthesis related to REA and BFT. These results enlighten some of the molecular processes involved in muscle and fat deposition, which are economically important carcass traits for beef production.
The amount of intramuscular fat can influence the sensory characteristics and nutritional value of beef, thus the selection of animals with adequate fat deposition is important to the consumer. There ...is growing knowledge about the genes and pathways that control the biological processes involved in fat deposition in muscle. MicroRNAs (miRNAs) belong to a well-conserved class of non-coding small RNAs that modulate gene expression across a range of biological functions in animal development and physiology. The aim of this study was to identify differentially expressed (DE) miRNAs, regulatory candidate genes and co-expression networks related to intramuscular fat (IMF) deposition. To achieve this, we used mRNA and miRNA expression data from the Longissimus dorsi muscle of 30 Nelore steers with high (H) and low (L) genomic estimated breeding values (GEBV) for IMF deposition.
Differential miRNA expression analysis between animals with extreme GEBV values for IMF identified six DE miRNAs (FDR 10%). Functional annotation of the target genes for these microRNAs indicated that the PPARs signaling pathway is involved with IMF deposition. Candidate regulatory genes such as SDHAF4, FBXO17, ALDOA and PKM were identified by partial correlation with information theory (PCIT), phenotypic impact factor (PIF) and regulatory impact factor (RIF) co-expression approaches from integrated miRNA-mRNA expression data. Two DE miRNAs (FDR 10%), bta-miR-143 and bta-miR-146b, which were upregulated in the Low IMF group, were correlated with regulatory candidate genes, which were functionally enriched for fatty acid oxidation GO terms. Co-expression patterns obtained by weighted correlation network analysis (WGCNA), which showed possible interaction and regulation between mRNAs and miRNAs, identified several modules related to immune system function, protein metabolism, energy metabolism and glucose catabolism according to in silico analysis performed herein.
In this study, several genes and miRNAs were identified as candidate regulators of IMF by analyzing DE miRNAs using two different miRNA-mRNA co-expression network methods. This study contributes to the understanding of potential regulatory mechanisms of gene signaling networks involved in fat deposition processes measured in muscle. Glucose metabolism and inflammation processes were the main pathways found in silico to influence intramuscular fat deposition in beef cattle in the integrative mRNA-miRNA co-expression analysis.